(整理)ANSYS材料模型.

1、精品文档第七章材料模型ANSYS/LS-DYN怠才40多种材料模型，它们可以表示广泛的材料特性，可用材料如下所示。本 章后面将详细叙述材料模型和使用步骤。对于每种材料模型的详细信息，请参看Appendix B,Material Model Examples或LS/DYNA Theoretical Manual »的第十六章(括号内将列出与每种模型相对应 的LS-DYN阳料号)。线弹性模型 各向同性(#1) 正交各向异性(#2) 各向异性(#2) 弹性流体(#1)非线弹性模型 Blatz -ko Rubber(#7) Mooney-Rivlin Rubber(#27) 粘弹性(#6)非

2、线性无弹性模型 双线性各向同性(#3) 与温度有关的双线性各向同性(#4) 横向各向异性弹塑性(#37) 横向各向异性FLD (#39) 随动双线性(#3) 随动塑性(#3) 3 参数 Barlat (#36) Barlat各向异性塑性(#33) 与应变率相关的募函数塑性(#64) 应变率相关塑性(#19) 复合材料破坏(#22) 混凝土破坏(#72) 分段线性塑性(#24) 募函数塑性(#18)压力相关塑性模型 弹-塑性流体动力学(#10) 地质帽盖材料模型(#25)泡沫模型 闭合多孔泡沫(#53) 粘性泡沫(#62) 低密度泡沫(#57) 可压缩泡沫(#63) Honeycomb(#26)

3、需要状态方程的模型 BammaiH性(#51) Johnson-Cook 塑性(#15) 空材料(#9) Zerilli -Armstrong(#65) Steinberg(#11)离散单元模型 线弹性弹簧 普通非线性弹簧 非线性弹性弹簧 弹塑性弹簧 非弹性拉伸或仅压缩弹簧 麦克斯韦粘性弹簧 线粘性阻尼器 非线粘性阻尼器 索（缆）（#71）刚性体模型 刚体（#20）7.1定义显示动态材料模型用户可以采用ANSYS命令 MP , MPTEMP, MPDATA , TB , TBTEMP 和 TBDATA 以及ANSYS/LS-DYNAr令EDMP来定义材料模型。下一节 显动态材料模型的描述,说明

4、了怎样使用命令定义每种材料模型的特性。通过GUI路径定义材料模型比使用命令直接得多：1 .选择菜单路径 Main Menu>Preprocessor>Material Props>Material Models.Define Material Model Behavior对话框出现。注-如果不事先定义 ANSYS/LS-DYNAI元类型，那么就不能定义 ANSYS/LS-DYNA料模型。2 .在Material Models Available 窗口的右侧，双击LS-DYNA然后选择一种材料模型种类：线 性、非线性、状态方程、离散单元特性或刚体材料。3 .双击一种材料的子目录

5、。例如，在非线性材料中，有弹性、非弹性和泡沫材料模型。4 .继续双击下面的材料分类直到数据输入对话框出现。框中的选项包括所有的材料模型，它对所 选的材料模型都有效。5 .输入所需的值，单击 OK。然后在Materials Models Defined窗口左边就列出了材料模型的类型和号码。然后用户可以双击 Materials Models Defined窗口左边的材料模型使相关数据对话框出现。这样就可以修改其值。然后单击OK o用户可以选择Edit>Copy并指定新模型号来复制现有材料模型的内容，复制的材料模型以新模型 号列在Materials Models Defined 窗口左侧，其内

6、容与原材料模型内容相同。单击模型号选定它，然后选择Edit>Delete,可以删除材料模型。使用 GUI路径定义材料的详细信息，参看§1.2.4.4« ANSYS Basic Analysis Guide中的Material Model Surface,也可参看ANSYSOperations Guide» 的§ 4.2.1.10 Using Tree Structure 来获得材料模型界面结构层的详细信息。如果用户通过GUI路径来定义、修改、复制或删除材料模型，ANSYSI自动发出正确命令并将其写入log文件中。7.2显式动态材料模型的描述本节将详

7、细讲述每一种材料模型。每当提及“加载曲线 ID”时，就需要输入一条材料数据曲线 ID,用EDCURVEr令定义材料数据曲线，见第四章,Loading。当采用交互工作方式时，所有材料模型的可用特性都出现在材料模型对话框中。当使用批处理或命令流方式时，相应的命令都提供在这里。要保证定义材料属性为模型列出的，不要定义与模型无关的数据。7.2.1 线弹性模型7.2.1.1 各向同性弹性模型各向同性弹性模型。使用M喉令输入所需参数：MP, DEN秘密度MP, EXH弹性模量MP, NUX一泊松比此部分例题参看 B.2.1 , Isotropic Elastic Example:High Carbon S

8、teel 。7.2.1.2 正交各向异性弹性模型正交各向异性弹性模型。用 MP命令输入所需参数：MP, DENS密度MP, EX-弹性卞K量（EY, EZ）;需一值MP, NUX从泊松比（NUXY NUXZ ;需一值或MP, PRXF主泊松比（PRYZ PRXZ ;需一值MP , GX剪切IK量（GYZ,GXZ ;需一值当仅给定一个值时（例如，EX）其它值将自动定义（EY=EZ=EX o用EDLCS和EDMP , ORTHO命令定义材料坐标系统。如果没有给定材料坐标系统，材料特性将单元的I , J, L节点定义的材料轴保持正交各向异性（参看下图）。对于多层复合壳，用 TB , COMPT令代替

9、，并作为SHELL16弹元 实常数给定层性质。详细信息参看§7.2.3.11 Composite Damage Model.例题参看 B.2.2 Orthopic Elastic Example:Aluminum Oxide.7.2.1.3 各向异性弹性模型此种材料的描述需要全弹性矩阵。由于其对称性，仅需21种常数。这种材料仅对SOLID164单元和PLANE162I元有效（轴对称和平面应变问题）。用MP命令输入密度。用 TB , ANEL命令以上三角形式输入常数。用 EDLCS和EDMP, ORTHO 命令定义材料方向轴。如果没有定义材料坐标系，材料性质将与单元的I、J、L节点所定

10、义的材料轴保持正交各向异性（参看上面的单元坐标系图）。MP, DENS密度TB , ANELTBDATA, 1, C11,C12,C22,C13,C23,C33TBDATA, 7, C14,C24,C34,C44,C15,C25TBDATA, 13, C35, C45, C55, C16, C26, C36TBDATA, 19, C46, C56, C66当用户使用TBLIST显示材料类型的数据信息时，这些常数以下三角形式D出现而不是上三角形式C o这一矛盾不是计算错误；材料数据已准确传递给LS-DYNA!序。例题参看 B.2.3,Anisotropic Elastic Example:Cad

11、mium 。7.2.1.4 弹性流体模型使用此选项来模拟动态冲击载荷作用下盛满流体的容器。可以用MP命令输入密度（DENS ,用EDMP命令定义材料模型为弹性流体：MP, DENSEDMP, FLUID, MAT VAL1流体模型要求指定体积模量，可以在上述命令的VAL1域输入。除了使用EDMP外，用户也可用MP命令输入弹性模量（ER和泊才比（NUXY o然后程序将计算体积模量如下所示：MP, EXMP, NUXY3（1-2叨如果VAL1 （ EDMP内）、EX和NUXYTB指定了，VAL1将用作体积模量。7.2.2 非线性弹性模型7.2.2.1Blatz-ko弹性橡胶模型Blatz和ko定义

12、的超弹连续橡胶模型。该模型使用第二类Piola-Kirchoff 应力：1% =G产厂47乜其中，G-剪切模量，V一相对体积，v 一泊松比，一右柯西-格林应变张量，而4 Kronecker delta o用MP命令输入密度（DENS和剪切模量（GXY。例题参看 B.2.4,Blatz-Ko Example:Rubber 。7.2.2. 2 Mooney-Rivlin橡胶弹性模型不可压缩橡胶模型。它与 ANSYS勺Mooney-Rivlin 2-参数模型很相似。输入 G。，Gil和1?来定义应变能量密度函数:取，出一 3) + %出-3) + 0+ D区 _ I)22(1 - 2v)、和4是右柯

13、西-格林张量不变量精品文档用MP命令输入泊松比（I）和密度。（泊松比的值要比推荐的大一些，太小的值不能工作。）用TB和TBDATA命令输入Mooney-Rivlin常数，只允许一种温度下的数据， 并且必须放在数据表中 的1和2位置。TB , MOONE丫，0TBDATA, 2,TBDATA, 1,01如果不直接输入 to和C。,可以设这些常数为 0,然后用载荷曲线提供表格式单轴数据。程序将根据TBDATA命令的3-6项所输入的实验数据来计算这些常数。使用这种输入法，必须设TB命令的 TBOPT=2TB , MOONEY, , , 2TBDATA, 1, C1。（设为0,应用实验数据）TBDAT

14、A, 2,设为0,应用实验数据）TBDATA, 3,（试样测量长度TBDATA, 4,（试样测量宽度）TBDATA, 5,（试样厚度）TBDATA, 6,（载荷曲线ID）提供单轴数据的载荷曲线应使测量长度 AL随相应力的变化而变化。在压缩中，力和长度变化 须为负值。在拉伸中，力和测量长度变化须为正值。单轴方向的主拉伸比|儿由下式给出：'。一初始长度，L实际长度。或者可以通过设定测量长度、设置厚度和宽度为1.0 ,并且在测量长度变化处定义工程应变以及在有作用力的地方定义名义（工程）应力，从而输入应力 -应变曲线。在ANSYS/LS-DYNA：解的初始阶段，用最小二乘法来处理输入的实验数据

15、。例题参看 B.2.5 Mooney-Rivlin Example : Rubber。7.2.2 . 3粘弹性模型Herrmann和Peterson提出的线性粘弹性模型。模型采用偏量特性：f del （）|这里剪切松弛模量由下式给出：/）=，+（鸟在模型中，由体积V计算增量积分压力时，需事先进行弹性体积假设，即V: p=K 。用参数G©、 G。、K （体积模量）和 0来定义线粘弹性模型。用 TB , EVISC和TBDATA命令的46、47、48和61 项输入以上数据：TB , EVISCTBDATA, 46,-i.TBDATA, 47TBDATA, 48, KTBDATA, 61,

16、 1/0注-对于这种材料选项，必须用 MP命令定义密度（DENS。例题参看 B.2.6, Viscoelastic Example: Glass。7.2.3 非线性无弹性模型7.2.3.1 双线性各向同性模型使用两种斜率（弹性和塑性）来表示材料应力应变行为的经典双线性各向同性硬化模型 （与应变 率无关）。仅可在一个温度条件下定义应力应变特性。 （也有温度相关的本构模型； 参看Temperature Dependent Bilinear Isotropic Model ）。用 MP 命令输入弹性模量（Exx）,泊松比（NUXY 和密 度（DENS ,程序用EX和NUXY1计算体积模量（K）。用T

17、B和TBDATA命令的1和2项输入屈服 强度和切线模量：TB , BISOTBDATA, 1,（屈服应力）TBDATA, 2, 几（切线模量）例题参看 B.2.7,Bilinear Isotropic Plasticity Example: Nickel Alloy 。7.2.3.2 与温度相关的双线性各向同性模型应变率无关、用两种斜率（弹性和塑性）来表示材料应力应变特性的经典双线性各向同性硬化模 型。可以在六个不同的温度下定义应力应变行为。如果仅在一个温度下定义应力应变行为，就需假定 双线性各向同性材料模型（与应变率和温度无关）。可以通过输入较大的屈服强度值来以该模型、表 热-弹性材料。用M

18、P命令输入密度（DENS （温度无关）。用 MPTEM评口 MPDATA输入弹性模量（Exx）、泊 松比（NUXY却热胀系数（ALPX）（这些特性和温度有关）。用TB , BISO, , NTEMPTBTEMP以及TBDATA 命令的1和2项输入屈服强度和切线模量。屈服强度和切线模量必须相对于同一温度定义，在MPTEMP命令中输入。MP, DENSMPTEMP 1, T的%, TE峭，阳明网？MPDATA EX, ,1,第2 ,，MPDATA NUXY , 1,而弭映啊MPDATA ALPX , 1, ALPX,，"PX滥股TB , BISO）, , NTEMPNTEM阿为 2 至U

19、 6）TBTEMP,（第一个温度点）TBDATA, 1, 4 （屈服应力）ffTBDATA, 2, 皿（切线模量）TBTEMP, TEMR （第二个温度点）TBDATA, 1,内（屈服应力）TBDATA, 2, 口恤！（切线模量）（重复此形式NTEM欣）TBTEMP,匕叫棋翻p （最后一个温度点）TBDATA, 2,TBDATA, 1,（屈服应力）（切线模量）注-对于这些材料模型，必须提供足够范围的温度数据，使之能够覆盖分析中的实际温度。否 则，分析将会中止。7.2.3.3 横向各向异性硬化模型仅供壳单元和2-D单元使用的全迭代各向异性塑T模型。在此模型中，由 HILL给出的屈服函数 在平面应

20、力情况下简化如下：12 2R_ 2 衣 +1F= + 0厂 +2%VA t1式十回1恒1这里R为各向异性硬化参数，它为平面内的塑性应变率22 和平面外应变率 出 之比：rJIl 稚应力应变仅定义在一个温度下。用 MP命令输入弹性模量（Exx）,密度（DENS加泊才H匕（NUXY卜 用TB , PLAW , , , 7和TBDATA命令的1-4项输入屈服应力，切线模量，各向异性硬化参数以及 有效屈服应力相对于有效塑性应变的载荷曲线 ID号：TB , PLAW , , , 7TBDATA,1,（屈服应力）TBDATA, 3, R （各向异性硬化参数）TBDATA, 2,（切线模量）TBDATA,

21、4, LCID （屈服应力和塑性应变的载荷曲线ID）例题参看 B.2.8,Transversely Anisotropic Elastic Plastic Example: 1010 Steel。7.2.3.4 横向各向异性FLD硬化模型这种材料模型用于模拟各向异性材料的板料成形。仅考虑横向各向异性材料。对于此模型，可以用定义的载荷曲线来模拟流动应力和有效塑性应变的关系（EDCURVE） o另外，也可以定义成形极限图（也可用EDCURVE,如下图所示）。ANSYS/LS-DYNA!序用此图来计算材料所承受的最大应变 比。这一塑性模型仅在壳单元和 2-D单元中使用。这一模型遵循前边所述的横向各向

22、异性弹塑性模型 所介绍的塑性理论。理论基础可参考该模型。使用横向各向异性FLD模型，需用MP命令输入密度（DENS,弹性模量（Exx）和泊松比（NUXY 如下所示，可以用 TB , PLAW, , , 10和TBDATA命令中的1-5项定义其它参数。TB , PLAW , , , 10TBDATA, 1, 口？（屈服应力）ffTBDATA, 2, 口（切线模量）TBDATA, 3, R （各向异性硬化参数）TBDATA, 4, LCID1 （有效应力和塑性应变的载荷曲线）TBDATA, 5, LCID2 （定义FLD的载荷曲线）例题参看 B.2.9 , Transversely Anisotr

23、opic FLD Example : Steinless SteelPLANE STRAIN %7.23.5 双线性随动模型（与应变率无关）经典的双线性随动硬化模型，用两个斜率（弹性和塑性）来表示材料的应力应 变特性。用MP命令输入弹性模量（Exx）,密度（DENS和泊才H匕（NUXY。可以用TB , BKIN和 TBDATA命令中的1-2项输入屈服强度和切线模量：TB , BKINTBDATA,1,（屈服应力）TBDATA,2,抽（切线模量）例题参看 B.2.10 , Bilinear Kinematic Plasticity Example: Titanium Alloy7.2.3.6塑性

24、随动模型各向同性、随动硬化或各向同性和随动硬化的混合卞K型，与应变率相关，可考虑失效。通过在 0（仅随动硬化）和1 （仅各向同性硬化）间调整硬化参数0来选择各向同性或随动硬化。应变率用Cowper-Symonds模型来考虑，用与应变率有关的因数表示屈服应力，如下所示：这里1+ C（为十您*的一初始屈服应力，E一应变率，C和P-Cowper Symonds为应变率参数。一有效塑%性应变,塑性硬化模量，由下式给出:应力应变特性只能在一个温度条件下给定。用MP命令输入弹性模量（Exx）,密度（DENS和泊松比（NUXY。用TB , PLAW , , , 1和TBDATA命令中的1-6项输入屈服应力，

25、切线斜率，硬化 参数，应变率参数 C和P以及失效应变：如下所示，可以用TB , PLAW, , , 10和TBDATA命令中的1-5项定义其它参数。TBDATA,1,（屈服应力）TBDATA,2,侬（切线模量）TBDATA,3,B （硬化参数）TBDATA,4,C （应变率参数）TBDATA,5,P （应变率参数）TBDATA,6,/ （失效应变）例题参看B.2.11 , Plastic Kinematic Example : 1018 Steel。7.2.3.7 3-参数Barlat模型由 BarlatLian提出的各向异性塑性模型，用于平面应力条件下的铝质薄板模型。使用了指数和线性硬化法则

26、。平面应力下各向异性屈服准则定义为:2（5）常二*1+舄+M& -+ ,2占“丁一屈服应力，a和c各向异性材料模型，m-Barlat常数，和舄 定义为勺二阳+户用其中h和p为附加各向异性材料常数。对于指数硬化选项，材料屈服强度给定如下:k屈服系数，卜初始屈服应变，耳塑性应变，n硬化系数。所有各向异性材料常数, 除p隐含定义外，都由Barlat和Lian定义的宽厚应变比（R）决定：c=2-a对于任意角4的宽厚应变比可由下式计算;MP命令输入沿4万向的单轴向拉伸应力。仅在同一个温度下给定应力应变特性。用 弹性模量（Exx）,密度（DENS和泊才H匕（NUXY。硬化7B则类型HR （线T为1

27、或指数为2）,切 线模量（HR=D或屈服系数（HR=2 ,屈服应力（HR=1或硬化系数（HR=2 , Barlat指数，m,厚度和宽度方向的应变比,和小0以及正交各向异性材料轴，用TB , PLAW, , ,3和TBDATA命令的1-8项输入:最后一项CSID有两个有效值：0 （缺省）和2,如果CSID=0,局部坐标系由单元节点I,J和L 定义（如上图所示）；如果 CSID=2,材料轴由EDLCS命令给定的局部坐标系决定（对于确定轴向的 详细信息，请参看命令的描述）。在定义材料特性之前，必须用EDLCS定义局部坐标系，然后执行EDMPRTHO,VAL1,其中VAL1值为EDLCS命令定义的坐标

28、系标号。7.2.3.8 Barlat各向异性塑性模型由Barlat , lege和Berm发展的各向异性塑性模型，用于模拟成形过程的材料特性， 各向异性屈 服函数定义如下：这里m为流动指数；研为对称矩阵 U的主值，这里a、b、c、f、g和h代表各向异性材料常数，当 a=b=c=f=g=h=1,就会模拟各向同性材料行 为，而屈服表面就会简化为 Tresca表面（m=1和Von Mises表面（m=2或4）,对于此材料选项， 屈服强度由下式给出：取 T/3H%S犷= 1，3b （% - %- bjy）工=1咻仁一%）75-%）L二凡这里k是强度系数，是塑性应变，*0是初始屈服应变，n是硬化系数，仅

29、在同一温度下定义应力、应变特性。用 MP命令输入弹性模量（Exx）,密度（DENS和泊才比（NUXY ,强度系数， 初始屈服应变，硬化系数，流动指数和Barlat各向异卜常数a-h ,用TB ,PLAW,6 和TBDATA命令的第1 10项输入。TB , PLAW,6TBDATA 1, k （强度系数）TBDATA 2,% （初始应变）TBDATA 3, n （硬化系数）TBDATA 4, m （流动指数（Barlat）TBDATA 5, aTBDATA 6, bTBDATA 7, cTBDATA 8, fTBDATA 9, gTBDATA 10, h例题参看 B.2.13 , Barlat

30、Anisotropic Plasticity Example:2008-T4 Aluminum。7.2.3.9 应变率敏感的募函数式塑性模型与应变率相关的塑性模型，主要用于超塑性成形分析，该模型遵循Ramburgh-Osgood本构关系:.这里£ 应变；应变率；m-硬化系数；k 材料常数；n应变率敏感系数。应力应变关系只能定义于一个温度下。用 MP命令输入弹性模量（EXX ,密度（DENS和泊才比（NUXY 用TB ,PLAW,4和TBDATA命令的第14项定义材料常数、硬化系数、应变率敏感系数及初始应 变率。TB , PLAW,4TBDATA 1, k （材料常数）TBDATA 2

31、, m （硬化系数）TBDATA 3, n （应变率灵敏系数）TBDATA 4,瓦（初始应变率）例题参看 B.2.14 , Rate Sensitive Powerlaw Plasticity Example:A356 Aluminum 。7.2.3.10 应变率相关的塑性模型应变率相关各向同性塑性模型主要用于金属和塑性成形分析，在此模型中，载荷曲线用来描述初始屈服强度 °。与有效应变率之间的函数关系。屈服应力定义如下：%二CT温+E成f式中 讶初始屈服强度，声有效应变率，严有效塑性应变，应力应变特性仅定义于同一温度下。用MP命令输入弹性模量（EXX 、密度（DENS和泊松比（NUX

32、Y o定义初始屈服应力和有效应变率的载荷曲线号，切线模量，定义弹性模量和有效应变率 的载荷曲线号，定义切线模量和有效应变率的载荷曲线ID,定义Von Misess失效应力和有效应变率的载荷曲线号，用TB ,PLAW,5和TBDATA命令的第1 5项输入。对于壳单元，可在第 6项中给 定Mn Time,取代第5项中的LCID4定义材料失效。Mn Time为自动删除单元的最小步长。TB , PLAW,5TBDATA 1, LCID1 （定义初始屈服应力和有效应变率的载荷曲线ID）TBDATA 2, E tan （切向（塑性硬化）模量 ）TBDATA 3, LCID2 （定义弹性模量和有效应变率的载

33、荷曲线ID）TBDATA 4, LCID3 （定义切线模量和有效应变率的载荷曲线ID）TBDATA 5, LCID4 （定义Von Misess 失效应力和有效应变率的载荷曲线ID）TBDATA 6, Mn Time （自动删除单元的最小步长，仅用于壳单元）例题参看 B.2.15 , Strain Rate Dependent Plasticity Example;4140 Steel。7.2.3.11 复合材料破坏模型此材料模型是由Chang & Chang发展的复合材料失效模型，模型采用如下5个参数：S11由向拉伸强度S2= 黄向拉伸强度512=剪切强度C2= 黄向压缩强度二=非线

34、性剪切应力参数所有参数均由实验确定，用MP命令输入弹性模量（Exx,Eyy,Ezz ）、剪切模量（Gxy,Gyz,Gxz）, 密度（DENS和泊松比（NUXY,NUYZ,NUXZ）压缩失效时的体积模量、剪切强度、轴向拉深强度、横向 拉深强度、横向压缩强度以及非线性剪切应力参数用TB ,COMP和TBDATA命令的第16项输入：TB , COMPTBDATA 1, KFAIL （压缩失效时的体积模量）TBDATA 2, S12 （剪切强度）TBDATA 3, S1 （轴向拉伸强度）TBDATA 4, S2 （横向拉伸强度）TBDATA 5, C2 （横向压缩强度）TBDATA 6,Q'（

35、非线性剪切应力参数）注-关于LS-DYNA材料模型# 22 （复合破坏）的详细信息，请参考 LS-DYNA TheoreticalManual。即使不使用失效特性，多层复合薄片也要求此模型。薄片特性定义为SHELL 163的实常数。7.2.3.12 混凝土破坏模型此模型用于分析承受混合冲击载荷的刚劲加强混凝土材料。这一模型要求混凝土和加强材料常数以及状态方程（有关状态方程的详细信息参考7.2.6,Equation of State Models ）。用MP命令输入密度（DENS和泊松比（NUXY用TB , CONC R, , , 2命令和TBDATA命令的1-78项输入下列 值：TB ,CON

36、CR,2TBDATA1, 0f （失效的最大主应力）TBDATA,2,4 （内聚常数）TBDATA3,（压力硬化系数）TBDATA4,（压力硬化系数）TBDATA5,TBDATA6,TBDATA7,TBDATA8,TBDATA,TBDATA,TBDATA,TBDATA,TBDATA,TBDATA,TBDATA,TBDATA,TBDATA,TBDATA,TBDATA,TBDATA,（屈服内聚力）10,11,12,13,14,15,16,18,19,（屈服极限的压力硬化系数）（屈服极限的压力硬化系数）（失效材料的压力硬化系数）9,（失效材料的压力硬化系数）（破坏比例因子）（单轴拉伸的破坏比例因子）

37、（三轴拉伸的破坏比例因子）PRE （加强筋的百分比） （加强筋的弹性模量）肌叫（加强筋的泊松比）（初始屈服应力）回（切线模量）17,LCP （主材料速率敏感度的载荷曲线ID）LCR （加强筋速率敏感度的载荷曲线ID）20-32,'13 （破坏函数 1-13）精品文档TBDATA,3,C （应变率参数）精品文档TBDATA,33-45 ,（比例因子1-13）TBDATA,46,GAMA温度常数）它是一个很常用的塑性准则，特Cowper-Symbols模型考虑应变率7.2.3. 13分段线性塑性模型多线性弹塑性材料模型，可输入与应变率相关的应力应变曲线 别用于钢。采用这个材料模型，也可根据

38、塑性应变定义失效。采用 的影响，它与屈服应力的关系为：i+（“ +这里£_有效应变率，c和p应变率参数,常应变率处的屈服应力，而是基于有效塑性应变的硬化函数。用MP命令输入弹性模量（Exx）,密度（DENSW口泊才H匕（NUXY）。用TB , PLAW , , , 8和TBDATA命令的1-7项输入屈服应力、切线模量、失效的有效真实塑性应 变、应变率参数 G应变率参数P、定义有效全应力相对于有效塑性真应变的载荷曲线ID以及定义应变率缩放的载荷曲线ID。TB ,PLAW, 8TBDATA2, 月血（切线模量）TBDATA3,（失效时的有效塑性真应变）TBDATA4,C （应变率参数）T

39、BDATA5,P （应变率参数）TBDATA6,LCID1 （定义全真应力相对于塑性真实应变的载荷曲线）TBDATA7,LCID2 （关于应变率缩放的载荷曲线）注-如果采用载荷曲线LCID1,则用TBDATA命令输入的屈服应力和切线模量将被忽略。 另外， 如果C和P设为0,则略去应变率影响。如果使用 LCID2,用TBDATA命令输入的应变率参数 C和P 将被覆盖。只考虑真实应力和真实应变数据。在数据曲线一节中讲述了此种类型的例题。注-例题参看 B.2.16 , Piecewise Linear Plasticity Example : High Carbon Steel 。7.2.4. 14

40、募函数塑性模型用于金属和塑性成形分析的与应变率有关的塑性模型。该模型提供各向同性硬化的弹塑性行为。 并且它用一个包括Cowper-Symbols乘子的募函数本构关系来描述应变率的影响：11 rm i+ - c /+s"其中 应变率，C, PCowper-Symbols应变率参数， 弹性应变， D 有效塑性应变，k强度系数，n硬化系数，仅能在一个温度下指定应力应变特性。用 MP命令输 入弹性卞K量（Exx）、密度（DENSW口泊松比（NUXY）。用TB , PLAW , , , 2和TBDATA命令的1-4 项输入强度系数、硬化系数和应变率参数 C和P:TB, PLAW,2TBDATA

41、1,k（强度系数）TBDATA2,n （硬化系数）精品文档TBDATA4,P （应变率参数）例题参看 B.2.17 , Powerlaw Plasticity Example:Aluminum 1100。7.2.5. 力相关的塑性模型7.2.5.1 弹塑性流体动力学模型该模型用于模拟承受大应变的材料，这里塑性特性可以由一系列数据点定义或屈服应力和切线模 量定义。如果不指定有效真实塑性应变和有效真实应力数据，屈服强度将按下式计算（依据各向同性 硬化）：|。/生+的根据杨氏模量和切线模量可计算塑性硬化模量£%:纥=E-Et如果指定了有效真实塑性应变和应力值，应力应变特性可以由有效真实应力

42、与真实塑性应变曲线 的数据点定义。可以最多定义16个数据点。如果应变值超过了最大输入值，将使用线性插值；因此， 需输入其它值来覆盖分析中所有的应变值。用MP命令输入密度（DENS）、弹性卞K量（EX）和剪切模量（GXY。用TB ,PLAW,9 和TBDATA命令的第1-45项输入下列参数：TB ,PLAW,9TBDATA1, % （初始屈服应力）TBDATA2, 口 k （硬化模量）TBDATA3,PC （截断压力值）TBDATA4, 1（失效应变）TBDATA5-20, W '16 （有效应变数据曲线值）TBDATA21-36, "1 "16 （有效应力数据曲线值

43、）TBDATA37,（状态常数的线性多项式方程）TBDATA38,（状态常数的线性多项式方程）TBDATA39, G （状态常数的线性多项式方程）TBDATA40, C?（状态常数的线性多项式方程）TBDATA41,（状态常数的线性多项式方程）TBDATA42, 5 （状态常数的线性多项式方程）TBDATA43, 6 （状态常数的线性多项式方程）TBDATA44, % （初始内能）TBDATA45,匕（初始相又t体积）注意TBDATA命令指定的37-45的常数与状态模型的线性多项式方程相同。详细信息请参看 § 7.2.6 , Equation of State Models 。12.

44、42 地质帽盖模型该模型是一种用于地质力学问题或诸如混凝土材料分析的非粘性、双常量材料模型。该模型中， 双常量帽盖理论又被扩展到包括非线性随动硬化。下面将讨论扩展的帽盖模型及其参数。精品文档Sur&jce of tine Two-invariant C理 Model.41XT0图7-1用应力张量不变量来描述帽盖模型。从偏量应力得出偏应力张量第二不变量的平方根如下所示:1此为变形或剪切力的客观标量尺寸。应力/1-第一不变量是应力张量的轨迹。帽盖模型包括压力人空间的三个表面，如图 7-1Surface of the Two-invariant CapModel所示。表面失效包络,?是极限表

45、面，而是拉伸中止值。的函数形式如下:Z 二 M(兄(J), 丁冠卿)这里& 士由下式给出:月伉)三好"型(-4)+跖。这一失效包络面固定在1空间，因此，如果不存在随动硬化就不会硬化，接着，在图中，有一个帽盖表面下式给出:这里%由下式给出"（人出三"也取厂与耳.£君是帽盖表面和A轴的交叉:x（旷上+叫尬而L （k）定义为由硬化准则，硬化参数 k和塑性体积的变化 有关，暧二印1一E即一。（工（同一工4|在几何上，认为R为帽盖表面和失效表面交叉处的坐标，最后，有一个截止拉伸表面，在图中表示为 h,函数 h由下式给出式中T为输入的材料参数，它来定义材料所

46、支持的最大静水张力，区;一1处的弹性区域由上面的失效包络面，左边的拉伸截止表面和右边的帽盖表面来定义边界。用MP命令输入密度（DENS和剪切模量（GXY。用TB , GCA脸令和TBDATA命令的1-13项输入下列参数。TB , GCAPTBDATA, 1, K （体积模量）TBDATA, 2, a （失效包络参数）TBDATA, 3, 0 （失效包络线性系数）TBDATA, 4, 丫（失效包络指数系数）TBDATA, 5, 0 （失效包络指数）TBDATA, 6, R （帽盖表面中心线比率）TBDATA, 7, D （硬化率指数）TBDATA, 8, W（硬化率系数）TBDATA, 9, X

47、。（硬化率指数）TBDATA, 10, C （动态硬化系数）TBDATA, 11, N （动态硬化参数）TBDATA, 12, Ftype （公式标志：1表示土和混凝土， 2表示石头）TBDATA, 13,Toff （拉伸截止值；Toff0,在压缩中为正）对于该种材料的详细信息请参看 LS-DYNA Theoretical Manual »。例题参看 B.2.18 , Geological Cap Example : SRI Dynamic Concrete 。7.2.5泡沫模型1.1.1.1 1闭合多孔泡沫模型刚性、闭合多孔、低密度聚氨酯泡沫材料模型通常用于汽车设计的撞击限制器模型

48、。该模型与 honeycomb很相似，在体积压缩达到之前，所有应力张量的分量都不耦合。但与honeycomb不同的是, 闭合多孔泡沫是各向同性的，还受约束的空气压力的影响，材料模型定义应力为：这里一轮廓应力，FQ初始泡沫应力，状一泡沫与聚合物密度之比，Kroneckerdelta,体积应变定义如下：了 一 'T + 10这里尸相对体积,一初始体积应变。屈服条件使用试探主应力，定义如下:这里a,b,c为用户自定义常数。只能在同一温度下定义应力应变特性。用 MP命令输入弹性模 量（Exx）和密度（DENS）,假设该模型的泊松比为 0,用TB ,FOAM,1和TBDATA命令的1-6项输入

49、屈服应力常数a,b和c,初始泡沫压力，泡沫与聚合物密度之比以及初始体积应变，如下示：TB ,FOAM,1TBDATA1,aTBDATA2,bTBDATA3,cTBDATA4, P0（初始泡沫压力）TBDATA5,夕（泡沫与聚合物密度之比）TBDATA6, （初始体积应变）1.1.1.2 粘性泡沫模型用于撞击模型的能量吸收泡沫材料。该模型包括并行的非线性弹性刚度和粘性阻尼。 在粘性吸收 能量同时使用用弹性刚度限定整体撞击。弹也刚度 V ,初始粘性系数片都是相对体积的非线性函 数：耳二旷£1初始弹性刚度， 匕一初始粘性系数，/1,分别为弹性刚度和粘性系数的募指数。仅能在同一种温度下定义应

50、力应变曲线。用MP命令输入弹性刚度（Exx）,泊松比（NUXY和密度（DENS。用TB ,FOAM,3 和TBDATA命令的1-4项输入弹性刚度募指数，初始粘性系数，粘性弹性刚度（防止产生时间步问题）和粘性募函数：TB ,FOAM,3TBDATA, 1,1 （弹性刚度的募指数）TBDATA,2,（初始粘性系数）TBDATA, 3, 耳（粘性弹性刚度）TBDATA, 4,（粘性系数的募指数）1.1.1.3 低密度泡沫模型高度可压缩泡沫材料模型，常常用于衬垫材料如椅子坐垫。在压缩中，该模型假设存在伴随能量耗散的滞后卸载特性。拉伸过程中撕裂发生之前，该材料模型呈线性。对于单轴载荷，该模型假设在 横向

51、方向上无耦合。采用输入形状因子控制（滞后卸载因子（HU ,延迟常数（0 ）和卸载形状因子）： 就可以近似估计泡沫的卸载特性。仅可在同一温度下定义应力应变特性。用 MP命令输入弹性模量（Exx）和密度（DENS。用TB ,FOAM,2和TBDATA命令的1-8项输入滞后卸载因子，延迟常数， 粘性系数，卸载形状因子，达到中止应力时的失效选项以及体积粘度标志：TB ,FOAM,2TBDATA1,LCID （应力应变行为的载荷曲线 ID）TBDATA2,TC （拉伸截断应力，缺省 =1E20）TBDATA3,HU （滞后卸载因子：1.0-无能量耗散；0.0 -全部能量耗散）TBDATA4, 0 （延迟

52、常数）TBDATA5,DAMP （粘性系数，推荐值为 0.05至U 0.5）TBDATA6,SHAPE （形状卸载因子，缺省值=1）TBDATA7,FAIL （达到截断应力时的失效选项：0.0 -截断值处的拉伸应力；1.0-拉伸应力设为0）TBDATA8,BVFLAG（体积粘度特性标志：0.0 -没有体积粘度（推荐值），1.0-激活体积粘度）1.1.1.4 可压缩泡沫模型该模型用于边侧撞击的可压缩泡沫或其它周期效应不太重要的应用。该模型与应变率有关并且在单向压缩时泊松比为0。在公式中，弹性模量认为是常数且采用弹性特性修正应力：%应变率，E弹性木量，t时间，该模型包括在拉伸载荷作用下定义失效的拉

53、伸应力截断值。对于拉伸截断值以下的应力，该模型在拉伸和压缩载荷作用下有相同的反应。重要的是该截断值需有非零值防止材料在小拉伸载荷下产生破坏。用 MP命令输入弹性模量（Exx）,密度（DENS）和泊松比 （NUXY）。用TB , FOAM，- 4和TBDATA命令的1-3项输入应力体积应变曲线，拉伸截断值和粘 性阻尼系数。TB ,FOAM,4TBDATA1,LCID （应力体积应变曲线ID）TBDATA2,TC（拉伸截断值）TBDATA3,DAMP（粘性阻尼系数,）1.1.1.5 Honeycomb 泡沫模型用于Honeycomb#料的正交各向异性材料模型。在压缩之前，材料为正交异性的，应力张量分量 不耦合，弹性模量呈线性分布与相对体积关系如下：%=% +舶-%)G以+以*%）Honeycomb材料的弹性剪切模量。并且，V一相对体积（定义为当前体积与原始体积之比）4 全压缩Honeycomb的相对体积载荷曲线用于表示平均应力幅值随相对体积的变化。 每个曲线必须有相同的横坐标值。曲线可以 定义为相又t体积（V或体积应变（1-V）的函数。用MP命令输入弹性模量（Exx）,密度（DENS 和泊松比（NUXY。用TB ,HONEY和TBADATA命令的1-17项输入下列