ABAQUS子程序UMAT里弹塑本构的实现

合集下载

5.本构模型-UMAT-JC模型

5.本构模型-UMAT-JC模型

[Kg/m3]
[MPa]

A [MPa]
B [MPa]
n
C
M
数值
2.7×103 68.0×103 0.33
66.562
108.853 0.238 0.029 0.5
有限元模拟SHPB实验
为了模拟SHPB实验,在二维模型和三维模型中均施加了必要 的边界条件。在对称轴或对称面上施加了对称性边界条件,同时 保证压杆和试件可以沿轴线方向自由无约束的运动。压杆和试件 之间的接触为硬接触,光滑无摩擦。 为了确定输入应力脉冲的时间,进行了简单的计算。弹性材 料中纵波波速的计算公式为:
1-6
泊松比
塑性耗散比
7-12 塑性应变
A
B
n
13
C
M
弹性应变
等效塑性应变
UMAT 流程图
3 SHPB实验
分离式 Hopkinson 压杆( Split Hopkinson Pressure Bar , 简称 SHPB )实验是从经典 Hopkinson 实验基础之上发展而来的一 种实验技术,用来测量材料的动态应力 - 应变行为。该实验技术 的理论基础是一维应力波理论,通过测量两根压杆上的应变来推 导试件上的应力-应变关系。
Cd E

由此可以计算输入应力波在压杆中的传播速度为 Cd 5048m s
应力波前沿到达试件时的Mises应力 (t=1.98×10-4 s)

统一弹塑性本构模型在ABAQUS中的开发与应用

统一弹塑性本构模型在ABAQUS中的开发与应用

2010年4月 Rock and Soil Mechanics Apr. 2010

收稿日期:2009-05-27

第一作者简介:潘晓明,男,1979年生,博士研究生,主要从事隧道及地下结构方面。E-mail: pxm155138@

文章编号:1000-7598 (2010) 04-1092-07

统一弹塑性本构模型在ABAQUS 中的开发与应用

潘晓明1, 2 ,孔 娟1, 2, 杨 钊1, 2, 刘 成1, 2

(1.同济大学 地下建筑与工程系,上海 200092;2.同济大学 岩土及地下工程教育部重点实验室,上海 200092)

摘 要:基于统一弹塑性本构模型的有限元理论格式,根据ABAQUS 的UMAT 格式要求,编制相应的接口程序,将统一弹塑性本构模型引入ABAQUS 中。采用退化的统一强度模型(0b =时,为Mohr-Coulomb 模型)与ABAQUS 自带的Mohr-Coulomb 模型,对单轴试验和圆形硐室进行弹塑性分析,验证所开发材料子程序的正确性及高效性。考虑到统一强度模型的一般情况(0b ≠)和屈服面硬化条件,对圆形硐室进行弹塑性计算,得到应力场的变化规律。所提出的研究思路具有普遍性,为采用ABAQUS 平台进行本构模型的二次开发提供了借鉴和参考。 关 键 词:统一强度理论;屈服面;流动矢量;奇异点;应力拉回算法;ABAQUS 中图分类号:TD 353.6 文献标识码:A

Secondary development and application of unified elastoplastic

ABAQUS-材料本构模型及编程

ABAQUS-材料本构模型及编程

材料本构模型及编程-ABAQUS-UMAT

材料本构模型及编程实现:简介

1、什么时候用用户定义材料(User-defined material,UMAT)?

很简单,当ABAQUS没有提供我们需要的材料模型时。所以,在决定自己定义一种新的材料模型之前,最好对ABAQUS已经提供的模型心中有数,并且尽量使用现有的模型,因为这些模型已经经过详细的验证,并被广泛接受。

2、好学吗?需要哪些基础知识?

先看一下ABAQUS手册(ABAQUS Analysis User's Manual)里的一段话:

Warning:The use of this option generally requires considerable expertise.The user is cautioned that the implementation of any realist

ic constitutive model requires extensive development and testing.Initial testing on a single element model with prescribed traction lo ading is strongly recommended.

但这并不意味着非力学专业,或者力学基础知识不很丰富者就只能望洋兴叹,因为我们的任务不是开发一套完整的有限元软件,而只是提供一个描述材料力学性能的本构方程(Constitutive equation)而已。当然,最基本的一些概念和知识还是要具备的,比如

abaqus1用户材料子程序

abaqus1用户材料子程序

19 ABAQUS用户材料子程序(UMAT)

虽然ABAQUS为用户提供了大量的单元库和求解模型,使用户能够利用这些模型处理绝大多数的问题;但是现实世界毕竟十分复杂,ABAQUS不可能把所有可能出现的问题都包含进去。所以ABAQUS提供了大量的用户子程序(User Subroutine)。用户子程序允许用户在找不到合适模型的情况下自行定义符合自己问题的模型。这些用户子程序涵盖了建模从载荷到单元的几乎各个部分。

ABAQUS为用户提供的这个接口,允许用户通过自定义的子程序定制ABAQUS,以实现特定的功能。用户子程序具有以下的功能和特点:(1)如果ABAQUS的一些固有选项模型功能有限;用户子程序可以提高ABAQUS中这些选项的功能;(2)通常用户子程序是用FORTRAN语言的代码写成;(3)它可以以几种不同的方式包含在模型中;(4)由于它们没有存储在restart文件中,如果需要的话,可以在重新开始运行时修改它;(5)在某些情况下它可以利用ABAQUS允许的已有程序。

要在模型中包含用户子程序,可以利用ABAQUS执行程序,在abaqus执行程序中应用user选项指明包含这些子程序的FORTRAN源程序或者目标程序的名字。

提示:ABAQUS的输入文件除了可以通过ABAQUS/CAE的作业模块中提交运行外,还可以在ABAQUS Command窗口中输入ABAQUS执行程序直接运行:

ABAQUS job=输入文件名 user=用户子程序的Fortran文件名

ABAQUS/Standard和ABAQUS/Explicit都支持用户子程序功能,但是他们所支持的用户子程序种类不尽相同,读者在需要使用时请注意查询手册。

UMAT全过程技术篇

UMAT全过程技术篇

UMAT全过程-"技术篇"

[写在前面:

这篇文章是UMAT全过程-感想篇的姊妹篇,答应要给大家写的一篇帖子,同时也是为了记录自己的学习过程,与大家分享!首先指出,俺的"技术篇"--是加了引号的,因为确实称不上有多么大的技术含量,还

望大家莫笑偶!只不过一是跟那个感想篇形成一个对照,同时主要内容为自己编子程序过程中涉及的技术边边上的小问题的一些解决方法,供仿友们参考!偶不是谦虚,也不是一个低调的人,大家谢谢和支持的话,我先行谢过啦!更希望大家能提出质疑或者别的更好的办法,大家相互

交流,共同进步!]

----------------------------------------------------------------------------

*转*入*正*题*

第一部分:相关知识

[特别声明,这部分来自于华中科技大学杨曼娟同学的硕士学位论文,在此对作者表示感谢!--大家可以去知网下载]

----------------------------------------------------------------------------

1.ABAQUS中材料非线性问题的处理

ABAQUS中材料非线性问题用Newton-Raphson法来求解。首先将载荷分为若干个微小增量,结构受到一个微小增量△P。ABAQUS用与初始结构位移相对应的初始刚度矩阵K0和荷载增量△P计算出结构的在这一步增量后的位移修正Ca、修正后的位移值Ua和相应的新的刚度矩阵Ka。ABAQUS用新的刚度矩阵计算结构的内力Ia,荷载P和Ia的差值为迭代的残余力Ra,即Ra=P-Ia。如果Ra在模型内的每个自由度上的值都为零,如图2-2中的a点,则结构处于平衡状态。但在非线性问题中,通常Ra是不可能为零,ABAQUS为此设置了一个残余力容差。如果Ra小于这个数字,ABAQUS就认为结构的内外力是平衡的。一般这个缺省值取为平均内力的0.5%(如图2-2)。

基于Abaqus子程序的高分子材料本构关系实现

基于Abaqus子程序的高分子材料本构关系实现
mo d e l t o s i mul a t e t h e p o l y me r s .Th e UMAT/ VUMAT s ub r o u t i n e i s pr o v i d e d b y Aba q us t o he l p u s e r t o d e in f e ma t e r i a l c o n s t i t u t i v e mo d e l a c c o r di ng t o t he a c t u a l n e e d s . Th e s t r e s s - s t r a i n r e l a t i o n s h i p o f p o l y me r s l c a n b e i mp l e me n t e d p r e c i s e l y b y de in f e d VUMAT,a n d S O t h e d i f f e r e n c e be t we e n t he s i mu l a t i o n a n d t h e t e s t c a n b e r e d u c e d. Me a n whi l e,t h e d i f f e r e n t u n l o a d pa t h s c a n b e d e in f e d b y di f f e r e n t u n l o a d la f g s . Ke y wo r d s:p o l y me r s;c o ns t i t u t i v e r e l a t i o n;Ab a q us ;UMAT;VUMAT

非线性本构关系在abaqus中的实现

非线性本构关系在abaqus中的实现

非线性本构关系在abaqus中的实现

ABAQUS中非线性本构关系的实现可以通过定义一个UMAT子程序,通过计算输入的应力和应变,来实现非线性本构关系。具体步骤如下:

1. 在ABAQUS中定义一个UMAT子程序,用于计算应力和应变之间的非线性本构关系;

2. 在ABAQUS中定义材料模型,指定使用的UMAT子程序;

3. 在ABAQUS中定义材料参数,并将其传递给UMAT子程序

4. 在UMAT子程序中,根据输入的应力和材料参数,计算应变;

5. 将计算的应变传递给ABAQUS,ABAQUS根据计算的应变计算应力;

6. 将应力传递给UMAT子程序,重复步骤4和5,直到应力和应变收敛为止。

ABAQUS子程序UMAT里弹塑本构的实现

ABAQUS子程序UMAT里弹塑本构的实现

前言

有限元法是工程中广泛使用的一种数值计算方法。它是力学、计算方法和计算机技术相结合的产物。在工程应用中,有限元法比其它数值分析方法更流行的一个重要原因在于:相对与其它数值分析方法,有限元法对边界的模拟更灵活,近似程度更高。所以,伴随着有限元理论以及计算机技术的发展,大有限元软件的应用证变得越来越普及。

ABAQUS软件一直以非线性有限元分析软件而闻名,这也是它和ANSYS,Nastran等软件的区别所在。非线性有限元分析的用处越来越大,因为在所用材料非常复杂很多情况下,用线性分析来近似已不再有效。比方说,一个复合材料就不能用传统的线性分析软件包进行分析。任何与时间有关联,有较大位移量的情况都不能用线性分析法来处理。多年前,虽然非线性分析能更适合、更准确的处理问题,但是由于当时计算设备的能力不够强大、非线性分析软件包线性分析功能不够健全,所以通常采用线性处理的方法。

这种情况已经得到了极大的改善,计算设备的能力变得更加强大、类似ABAQUS这样的产品功能日臻完善,应用日益广泛。

非线性有限元分析在各个制造行业得到了广泛应用,有不少大型用户。航空航天业一直是非线性有限元分析的大客户,一个重要原因是大量使用复合材料。新一代波音 787客机将全部采用复合材料。只有像 ABAQUS这样的软件,才能分析包括多个子系统的产品耐久性能。在汽车业,用线性有限元分析来做四轮耐久性分析不可能得到足够准确的结果。分析汽车的整体和各个子系统的性能要求(如悬挂系统等)需要进行非线性分析。在土木工程业, ABAQUS能处理包括混凝土静动力开裂分析以及沥青混凝土方面的静动力分析,还能处理高度复杂非线性材料的损伤和断裂问题,这对于大型桥梁结构,高层建筑的结构分析非常有效。

5.本构模型-UMAT-JC模型解析

5.本构模型-UMAT-JC模型解析
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,DRPLDE, 1 DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME, 2 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT, 3 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) C INCLUDE 'ABA_PARAM.INC' C CHARACTER*80 CMNAME DIMENSION STRESS(NTENS),STATEV(NSTATV), 1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS), 2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1), 3 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3) user coding to define DDSDDE, STRESS, STwenku.baidu.comTEV, SSE, SPD, SCD and, if necessary, RPL, DDSDDT, DRPLDE, DRPLDT, PNEWDT

(完整word版)Abaqus弹塑性分析简单实例

(完整word版)Abaqus弹塑性分析简单实例

(完整word版)Abaqus弹塑性分析简单实例

Abaqus弹塑性分析简单实例

ABAQUS默认的塑性材料特性应用金属材料的经典塑性理论,采用MISES屈服面来定义各向屈服。金属材料的弹塑性行为可以简述如下:在小应变时,材料性质基本为线弹性,弹性模量E为常数;应力超过屈服应力后,刚度会显著下降,此时材料的应变包括塑性应变和弹性应变两部分;在卸载后,弹性应变消失,而塑性应变是不可恢复的;如果再次加载,材料的屈服应力会提高,即所谓的加工硬化。在abaqus中,等效塑性应变PEEQ大于0表明材料发生了屈服。在工程结构中,等效塑性应变一般不应超过材料的破坏应变。对于金属成形等大变形问题,应根据生产工艺要求来确定许可的等效塑性应变量。需要注意的是在比例加载时,大多数材料的PEMAG和PEEQ相等。这两个量的区别在于,PEMAG描述的是变形过程中某一时刻的塑性应变,与加载历史无关,而PEEQ是整个变形过程中塑性应变的累积结果。下面我们以单向压缩过程的模拟来演示ABAQUS弹塑性仿真设置。模型如图所示,压头用解析刚体来模拟,试样用SHELL来模拟。采用轴对称模型。试样的截面属性设置如下图所示,注意塑性应变必须从0开始。

在压头与试样之间定义无摩擦的接触。固定对称轴

上的径向位移U1和底边的轴向位移U2。压头是轴对称刚体,U2边界条件需要施加在压头的参考点上。设定两个分析步,第一个分析步让压头与试样建立平稳的接触,设置压头下移-5.001mm。第二个分析步,设定压头下移20mm。具体如下图所示:

提交分析,结果如下图所示:有限元在线因为专注所以卓越

ABAQUS使用解答_

ABAQUS使用解答_

Q: abaqus的图形如何copy?

A:file>print>file格式为png,可以用Acdsee打开。

Q: 用Abaqus能否计算[Dep]不对称的问题?

A: 可以,并且在step里面的edit step对话框other里面的matrix solver有个选项。

Q: 弹塑性矩阵【D】与ddsdde有何联系?

A: stress=D*stran;d(stress)=ddsdde*d(stran)。

Q: 在abaqus中,如果采用umat,利用自己的本构,如何让abaqus明白这种材料的弹塑性应变,也就是说,如何让程序返回弹性应变与塑性应变,好在output中输出,我曾想用最笨地方法,在uvarm中定义输出,利用getvrm获取材料点的值,但无法获取增量应力,材料常数等,研究了帮助中的例子,umatmst3.inp,umatmst3.for,他采用mises J2 流动理论,我在output history 显示他已进入塑性状态,但他的PE仍然为0!!?

A: 用uvar()勉强成功。

Q: 本人在用umat作本构模型时,

*static,

1,500,0.000001,0.1 此时要求的增量步很多,即每次增量要很小,

*static

1,500 时,在弹性向塑性过度时,出现错误,增量过大,出现尖点。?

A: YOU CAN TRY AS FOLLOWS:

*STEP,EXTRAPOLATION=NO,INC=2000000

*STATIC

0.001,500.0,0.00001,0.1。

混凝土塑性损伤模型及其ABAQUS子程序开发

混凝土塑性损伤模型及其ABAQUS子程序开发

混凝土塑性损伤模型及其ABAQUS子程序开发

一、本文概述

混凝土作为一种广泛使用的建筑材料,其力学行为一直是工程领域的研究热点。混凝土塑性损伤模型(Concrete Plasticity Damage Model)作为一种能够模拟混凝土在复杂应力状态下的非线性、弹塑

性及损伤行为的本构模型,对于准确预测混凝土结构的力学响应和破坏过程具有重要意义。本文旨在介绍混凝土塑性损伤模型的基本理论,以及如何利用ABAQUS软件的子程序开发功能,实现该模型在数值模

拟中的应用。

文章首先将对混凝土塑性损伤模型的基本原理进行阐述,包括模型的损伤演化方程、塑性流动法则以及相关的材料参数。随后,将详细介绍在ABAQUS软件中开发混凝土塑性损伤模型子程序的步骤和关

键技术,包括用户子程序的编写、模型参数的输入和输出处理等。通过具体的算例分析,文章将展示所开发子程序在模拟混凝土结构力学行为方面的应用效果,并与其他常用模型进行对比分析,以验证所开发子程序的准确性和可靠性。

本文旨在为从事混凝土结构数值模拟的研究人员和工程师提供

一套有效的混凝土塑性损伤模型子程序开发方法,以推动混凝土结构数值模拟技术的发展和应用。

二、混凝土塑性损伤模型的基本理论

混凝土塑性损伤模型是一种基于塑性力学和损伤力学的本构模型,用于描述混凝土在复杂应力状态下的力学行为。该模型能够考虑混凝土的塑性变形、刚度退化以及损伤演化,因此在结构分析和数值模拟中得到了广泛应用。

塑性流动理论:混凝土在受力过程中会发生塑性变形,这种变形是不可逆的。塑性流动理论通过引入塑性势函数和流动法则,描述了混凝土在塑性状态下的应力-应变关系。塑性势函数用于确定塑性应

abaqus umat于计算力学之应用

abaqus umat于计算力学之应用

abaqus umat是一种在计算力学中广泛应用的有限元分析软件。它可以通过用户自定义的子程序(也称为umat)进行材料本构关系的定义,使得在模拟复杂材料行为时能够更加精确地描述材料的非线性和非均

匀性等特性。abaqus umat能够有效地模拟材料的机械性能,并在工程领域具有广泛的应用。

1. 什么是abaqus umat?

abaqus umat是abaqus软件中用于用户自定义材料本构关系的子程序。它可以实现对材料行为的精确描述,包括材料的非线性、非均匀

性等特性。通过abaqus umat,用户可以自定义材料的本构关系和材料参数,以满足对于各种材料行为的精确模拟需求。

2. abaqus umat的实现原理

abaqus umat的实现依赖于有限元分析方法。用户可以通过编写程序,在abaqus中调用该程序来定义材料的本构关系。在有限元分析中,

材料的本构关系是描述材料应力和应变之间关系的重要参数,通过用

户自定义的umat程序,可以实现对材料行为的更为精确的描述。

3. abaqus umat的应用领域

abaqus umat在工程领域有着广泛的应用。例如在航空航天领域,abaqus umat可以用于模拟飞机结构的材料行为,预测飞机在不同载荷下的应力应变分布,进行疲劳分析等。在汽车工业中,abaqus umat可以用于模拟汽车结构在碰撞时的材料行为,以及进行车身强度

分析等。abaqus umat还被广泛应用于建筑、船舶、能源等领域,在模拟复杂材料行为时发挥着重要作用。

4. abaqus umat的优势

岩土弹塑性模型的回映算法及ABAQUS子程序开发

岩土弹塑性模型的回映算法及ABAQUS子程序开发

岩土弹塑性模型的回映算法及ABAQUS子程序开发

郭德伟

【摘要】有限元法被广泛应用于各类岩土工程的数值分析中。为工程的过程分析

和方案优化提供必要的计算支持。但岩土材料大多是弹塑性体,在进行有限元分析时,需要采用迭代方法对本构模型进行积分计算,以求出与实际边值问题相适应的应力场和位移场。一个优良的迭代方法可以加快数值收敛,得到更精确的计算成果。文章在前人的研究基础上,采用回映算法的数学原理,推导了岩土材料理想弹塑性模型有限元法的一般数值格式,并且将其应用于常用的D—P模型中,在ABAQUS有限元程序上进行了UMAT的二次开发。最后用一个数值算例验证了迭代算法的正确性。

【期刊名称】《四川建筑》

【年(卷),期】2012(032)006

【总页数】3页(P64-66)

【关键词】弹塑性模型;回映算法;D—P模型;UMAT子程序

【作者】郭德伟

【作者单位】中国土木工程集团有限公司,北京100038

【正文语种】中文

【中图分类】TU432

数值分析是岩土工程中最重要的分析手段之一,为岩土工程的设计和建设提供极大

的计算支持,因此,愈来愈受到工程界的重视和应用。岩土工程弹塑性分析的难点之一是岩土材料非线性弹塑性本构模型在计算机上的程序化,这需要对本构方程进行积分,以便得到新的应力增量。在复杂的应力路径加载条件下,很难给出材料(包括岩土材料)弹塑性应力增量的显示积分解析式,而只能进行数值积分,在过去的数十年,在这方面有大量的文献报道[1-8]。这些文献报道的方法大致分为两类:显示积分方法[1-2]和隐式向后欧拉算法[3-8]。隐式向后欧拉算法最早由 Krieg 等人提出[3],经过发展,目前已成为应用最广的算法。

[转载]abaqus用户子程序

[转载]abaqus用户子程序

[转载]abaqus⽤户⼦程序

原⽂地址:abaqus⽤户⼦程序作者:热爱结构设计

当⽤到某个⽤户⼦程序时,⽤户所关⼼的主要有两⽅⾯:⼀是ABAQUS提供的⽤户⼦程序的接⼝参数。有些参数是ABAQUS 传到⽤户⼦程序中的,例如SUBROUTINE DLOAD中的KSTEP,KINC,COORDS;有些是需要⽤户⾃⼰定义的,例如F。⼆是ABAQUS何时调⽤该⽤户⼦程序,对于不同的⽤户⼦程序ABAQUS调⽤的时间是不同的。有些是在每个STEP的开始,有的是STEP结尾,有的是在每个INCREMENT的开始等等。当ABAQUS调⽤⽤户⼦程序是,都会把当前的STEP和INCREMENT利⽤⽤户⼦程序的两个实参KSTEP和KINC传给⽤户⼦程序,⽤户可编个⼩程序把它们输出到外部⽂件中,这样对ABAQUS何时调⽤该⽤户⼦程序就会有更深的了解。

(⼦程序中很重要的就是要知道由abaqus提供的那些参量的意义,如下)

⾸先介绍⼏个⼦程序:

⼀.SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS, JLTYP,SNAME)

参数:

1. F为⽤户定义的是每个积分点所作⽤的荷载的⼤⼩;

2. KSTEP,KINC为ABAQUS传到⽤户⼦程序当前的STEP和INCREMENT值;

3. TIME(1),TIME(2)为当前STEP TIME和INCREMENT TIME的值;

4. NOEL,NPT为积分点所在单元的编号和积分点的编号;

5. COORDS为当前积分点的坐标;

ABAQUS_Fortran二次开发

ABAQUS_Fortran二次开发

目录

摘要............................................................... ABSTRACT ........................................................... I

1.绪论 0

1.1.课题的研究背景 0

1。2.............................................. 本文的研究内容和方法1 2。基于ABAQUS软件的二次开发 (2)

2。1。...................................................... ABAQUS介绍2

2.2。................................................ ABAQUS各模块简介2

2。3............................................. ABAQUS的二次开发平台4 2。4。........................................... ABAQUS的二次开发语言5 3。用户材料子程序UMAT (7)

3.1。................................................ UMAT开发环境设置7

3.2.UMAT注意事项 (8)

3。3.................................................... UMAT接口的原理9

相关主题
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1.1. 课题的研究背景..................................................................................................... 1 1.2. 本文的研究内容和方法......................................................................................... 2 2. 基于 ABAQUS 软件的二次开发..................................................................................... 3 2.1. ABAQUS 介绍 ....................................................................................................... 3 2.2. ABAQUS 各模块简介 ........................................................................................... 3 2.3. ABAQUS 的二次开发平台.................................................................................... 5 2.4. ABAQUS 的二次开发语言.................................................................................... 6 3. 用户材料子程序 UMAT................................................................................................... 8 3.1. UMAT 开发环境设置 ............................................................................................ 8 3.2. UMAT 注意事项 .................................................................................................... 9 3.3. UMAT 接口的原理 ............................................................................................... 10 3.4. UMAT 的使用方法 ............................................................................................... 12 4. 材料非线性问题 ............................................................................................................. 14 4.1. 材料的弹塑性本构关系....................................................................................... 14 4.2. 非线性有限元算法理论....................................................................................... 17 4.3. 增量理论常刚度法公式推导............................................................................... 20 4.4. 增量理论切线刚度法公式推导........................................................................... 21 5. UMAT 程序设计和编码 ................................................................................................. 25 5.1. 本构关系描述....................................................................................................... 25 5.2. 常刚度法程序设计............................................................................................... 27 5.3. 常刚度法程序编码............................................................................................... 29 5.4. 切线刚度法程序设计........................................................................................... 32 5.5. 切线刚度法程序编码........................................................................................... 36 5.6. 程序的调试........................................................................................................... 39 6. 程序验证 ......................................................................................................................... 40
这种情况已经得到了极大的改善,计算设备的能力变得更加强大、类似 ABAQUS 这 样的产品功能日臻完善,应用日益广泛。
非线性有限元分析在各个制造行业得到了广泛应用,有不少大型用户。航空航天 业一直是非线性有限元分析的大客户,一个重要原因是大量使用复合材料。新一代波 音 787 客机将全部采用复合材料。只有像 ABAQUS 这样的软件,才能分析包括多个子 系统的产品耐久性能。在汽车业,用线性有限元分析来做四轮耐久性分析不可能得到 足够准确的结果。分析汽车的整体和各个子系统的性能要求(如悬挂系统等)需要进 行非线性分析。在土木工程业, ABAQUS 能处理包括混凝土静动力开裂分析以及沥青混 凝土方面的静动力分析,还能处理高度复杂非线性材料的损伤和断裂问题,这对于大 型桥梁结构,高层建筑的结构分析非常有效。
瞬态、大变形、高级材料的碰撞问题必须用非线性有限元分析来计算。线性分析 在这种情况下是不适用的。以往有一些专门的软件来分析碰撞问题,但现在 ABAQUS 在 通用有限元软件包就能解决这些问题。所以,ABAQUS 可以在一个软件完成线性和非线 性分析。
ABAQUS 给用户提供了强大二次开发接口,尤其是在材料本构方面,给用户开发符 合实际工程的材料本构模型提供了强大帮助,本文将针对其用户材料子程序展开研究, 总结常用材料模型的开发方法。
1
6.1. 问题描述............................................................................................................... 41 6.2. 本构关系............................................................................................................... 42 6.3. ABAQUS 自带材料模型计算 ............................................................................. 42 6.4. 常刚度法的 UMAT 验证 ..................................................................................... 44 6.5. 切线刚度法的 UMAT 验证 ................................................................................. 46 6.6. 两种算法的比较分析........................................................................................... 48 7. 结论与展望 ..................................................................................................................... 52 7.1. 结论....................................................................................................................... 52 7.2. 展望....................................................................................................................... 52 致 谢 ...................................................................................................................................... 54 参考文献 ................................................................................................................................ 55 附 1:ABAQUS 自带弹塑性材料验证的 INP 文件 ........................................................... 56 附 2:用于算法验证的 INP 文件......................................................................................... 62
前言
有限元法是工程中广泛使用的一种数值计算方法。它是力学、计算方法和计算机 技术相结合的产物。在工程应用中,有限元法比其它数值分析方法更流行的一个重要 原因在于:相对与其它数值分析方法,有限元法对边界的模拟更灵活,近似程度更高。 所以,伴随着有限元理论以及计算机技术的发展,大有限元软件的应用证变得越来越 普及。
目录
摘 要 .........................................................................................................................................I ABSTRACT .......................................................................................................................ຫໍສະໝຸດ Baidu.... II 1. 绪论 ................................................................................................................................... 1
ABAQUS 软件一直以非线性有限元分析软件而闻名,这也是它和 ANSYS,Nastran 等 软件的区别所在。非线性有限元分析的用处越来越大,因为在所用材料非常复杂很多 情况下,用线性分析来近似已不再有效。比方说,一个复合材料就不能用传统的线性 分析软件包进行分析。任何与时间有关联,有较大位移量的情况都不能用线性分析法 来处理。多年前,虽然非线性分析能更适合、更准确的处理问题,但是由于当时计算 设备的能力不够强大、非线性分析软件包线性分析功能不够健全,所以通常采用线性 处理的方法。
相关文档
最新文档