LS-DYNA第七章：材料模型

第七章材料模型ANSYS/LS-DYNA包括40多种材料模型，它们可以表示广泛的材料特性，可用材料如下所示。本章后面将详细叙述材料模型和使用步骤。对于每种材料模型的详细信息，请参看AppendixB,MaterialModelExamples或《LS/DYNATheoreticalManual》的第十六章（括号内将列出与每种模型相对应的LS-DYNA材料号）。线弹性模型·各向同性（#1）·正交各向异性（#2）·各向异性（#2）·弹性流体（#1）非线弹性模型·Blatz-koRubber(#7)·Mooney-RivlinRubber(#27)·粘弹性（#6）非线性无弹性模型·双线性各向同性（#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)需要状态方程的模型·Bamman塑性（#51）·Johnson-Cook塑性（#15）·空材料（#9）·Zerilli-Armstrong(#65)·Steinberg(#11)离散单元模型·线弹性弹簧·普通非线性弹簧·非线性弹性弹簧·弹塑性弹簧·非弹性拉伸或仅压缩弹簧·麦克斯韦粘性弹簧·线粘性阻尼器·非线粘性阻尼器·索（缆）（#71）刚性体模型·刚体（#20）7.1定义显示动态材料模型用户可以采用ANSYS命令MP，MPTEMP，MPDATA，TB，TBTEMP和TBDATA以及ANSYS/LS-DYNA命令EDMP来定义材料模型。下一节显动态材料模型的描述，说明了怎样使用命令定义每种材料模型的特性。通过GUI路径定义材料模型比使用命令直接得多：1．选择菜单路径MainMenuPreprocessorMaterialPropsMaterialModels.DefineMaterialModelBehavior对话框出现。注--如果不事先定义ANSYS/LS-DYNA单元类型，那么就不能定义ANSYS/LS-DYNA材料模型。2．在MaterialModelsAvailable窗口的右侧，双击LS-DYNA，然后选择一种材料模型种类：线性、非线性、状态方程、离散单元特性或刚体材料。3．双击一种材料的子目录。例如，在非线性材料中，有弹性、非弹性和泡沫材料模型。4．继续双击下面的材料分类直到数据输入对话框出现。框中的选项包括所有的材料模型，它对所选的材料模型都有效。5．输入所需的值，单击OK。然后在MaterialsModelsDefined窗口左边就列出了材料模型的类型和号码。然后用户可以双击MaterialsModelsDefined窗口左边的材料模型使相关数据对话框出现。这样就可以修改其值。然后单击OK。用户可以选择EditCopy并指定新模型号来复制现有材料模型的内容，复制的材料模型以新模型号列在MaterialsModelsDefined窗口左侧，其内容与原材料模型内容相同。单击模型号选定它，然后选择EditDelete,可以删除材料模型。使用GUI路径定义材料的详细信息，参看§1.2.4.4《ANSYSBasicAnalysisGuide》中的MaterialModelSurface,也可参看《ANSYSOperationsGuide》的§4.2.1.10UsingTreeStructure来获得材料模型界面结构层的详细信息。如果用户通过GUI路径来定义、修改、复制或删除材料模型，ANSYS将自动发出正确命令并将其写入log文件中。7.2显式动态材料模型的描述本节将详细讲述每一种材料模型。每当提及“加载曲线ID”时，就需要输入一条材料数据曲线ID，用EDCURVE命令定义材料数据曲线，见第四章，Loading。当采用交互工作方式时，所有材料模型的可用特性都出现在材料模型对话框中。当使用批处理或命令流方式时，相应的命令都提供在这里。要保证定义材料属性为模型列出的，不要定义与模型无关的数据。7.2.1线弹性模型7.2.1.1各向同性弹性模型各向同性弹性模型。使用MP命令输入所需参数：MP，DENS—密度MP，EX—弹性模量MP，NUXY—泊松比此部分例题参看B.2.1，IsotropicElasticExample:HighCarbonSteel。7.2.1.2正交各向异性弹性模型正交各向异性弹性模型。用MP命令输入所需参数：MP，DENS—密度MP，EX—弹性模量（EY，EZ）；需一值MP，NUXY—从泊松比（NUXY，NUXZ）；需一值或MP，PRXY—主泊松比（PRYZ，PRXZ）；需一值MP，GXY—剪切模量（GYZ,GXZ）;需一值当仅给定一个值时（例如，EX）其它值将自动定义（EY=EZ=EX）。用EDLCS和EDMP，ORTHO命令定义材料坐标系统。如果没有给定材料坐标系统，材料特性将单元的I，J，L节点定义的材料轴保持正交各向异性（参看下图）。对于多层复合壳，用TB，COMP命令代替，并作为SHELL163单元实常数给定层性质。详细信息参看§7.2.3.11CompositeDamageModel.例题参看B.2.2OrthopicElasticExample:AluminumOxide.7.2.1.3各向异性弹性模型此种材料的描述需要全弹性矩阵。由于其对称性，仅需21种常数。这种材料仅对SOLID164单元和PLANE162单元有效（轴对称和平面应变问题）。用MP命令输入密度。用TB，ANEL命令以上三角形式输入常数。用EDLCS和EDMP，ORTHO命令定义材料方向轴。如果没有定义材料坐标系，材料性质将与单元的I、J、L节点所定义的材料轴保持正交各向异性（参看上面的单元坐标系图）。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]。这一矛盾不是计算错误；材料数据已准确传递给LS-DYNA程序。例题参看B.2.3,AnisotropicElasticExample:Cadmium。7.2.1.4弹性流体模型使用此选项来模拟动态冲击载荷作用下盛满流体的容器。可以用MP命令输入密度（DENS），用EDMP命令定义材料模型为弹性流体：MP，DENSEDMP，FLUID，MAT，VAL1流体模型要求指定体积模量，可以在上述命令的VAL1域输入。除了使用EDMP外，用户也可用MP命令输入弹性模量（EX）和泊松比（NUXY）。然后程序将计算体积模量如下所示：MP，EXMP，NUXY如果VAL1（EDMP内）、EX和NUXY都指定了，VAL1将用作体积模量。7.2.2非线性弹性模型7.2.2.1Blatz-ko弹性橡胶模型Blatz和ko定义的超弹连续橡胶模型。该模型使用第二类Piola-Kirchoff应力：其中，G—剪切模量，V—相对体积，ν—泊松比，—右柯西-格林应变张量，而—Kroneckerdelta。用MP命令输入密度（DENS）和剪切模量（GXY）。例题参看B.2.4,Blatz-KoExample:Rubber。7.2.2.2Mooney-Rivlin橡胶弹性模型不可压缩橡胶模型。它与ANSYS的Mooney-Rivlin2-参数模型很相似。输入，和来定义应变能量密度函数：，和是右柯西-格林张量不变量。用MP命令输入泊松比（）和密度。（泊松比的值要比推荐的大一些，太小的值不能工作。）用TB和TBDATA命令输入Mooney-Rivlin常数，只允许一种温度下的数据，并且必须放在数据表中的1和2位置。TB，MOONEY，，，，0TBDATA，1，TBDATA，2，如果不直接输入和，可以设这些常数为0，然后用载荷曲线提供表格式单轴数据。程序将根据TBDATA命令的3-6项所输入的实验数据来计算这些常数。使用这种输入法，必须设TB命令的TBOPT=2：TB，MOONEY，，，，2TBDATA，1，（设为0，应用实验数据）TBDATA，2，（设为0，应用实验数据）TBDATA，3，（试样测量长度）TBDATA，4，（试样测量宽度）TBDATA，5，（试样厚度）TBDATA，6，（载荷曲线ID）提供单轴数据的载荷曲线应使测量长度随相应力的变化而变化。在压缩中，力和长度变化须为负值。在拉伸中，力和测量长度变化须为正值。单轴方向的主拉伸比由下式给出：—初始长度，L—实际长度。或者可以通过设定测量长度、设置厚度和宽度为1.0，并且在测量长度变化处定义工程应变以及在有作用力的地方定义名义（工程）应力，从而输入应力-应变曲线。在ANSYS/LS-DYNA求解的初始阶段，用最小二乘法来处理输入的实验数据。7.2.2.3粘弹性模型Herrmann和Peterson提出的线性粘弹性模型。模型采用偏量特性：这里剪切松弛模量由下式给出：在模型中，由体积V计算增量积分压力时，需事先进行弹性体积假设，即V：p=K。用参数、、K（体积模量）和β来定义线粘弹性模型。用TB，EVISC和TBDATA命令的46、47、48和61项输入以上数据：TB，EVISCTBDATA，46，TBDATA，47，TBDATA，48，KTBDATA，61，1/β注--对于这种材料选项，必须用MP命令定义密度（DENS）。例题参看B.2.6,ViscoelasticExample：Glass。7.2.3非线性无弹性模型7.2.3.1双线性各向同性模型使用两种斜率（弹性和塑性）来表示材料应力应变行为的经典双线性各向同性硬化模型（与应变率无关）。仅可在一个温度条件下定义应力应变特性。（也有温度相关的本构模型；参看TemperatureDependentBilinearIsotropicModel）。用MP命令输入弹性模量（Exx），泊松比（NUXY）和密度（DENS），程序用EX和NUXY值计算体积模量（K）。用TB和TBDATA命令的1和2项输入屈服强度和切线模量：TB，BISOTBDATA，1，（屈服应力）TBDATA，2，（切线模量）例题参看B.2.7,BilinearIsotropicPlasticityExample：NickelAlloy。7.2.3.2与温度相关的双线性各向同性模型应变率无关、用两种斜率（弹性和塑性）来表示材料应力应变特性的经典双线性各向同性硬化模型。可以在六个不同的温度下定义应力应变行为。如果仅在一个温度下定义应力应变行为，就需假定双线性各向同性材料模型（与应变率和温度无关）。可以通过输入较大的屈服强度值来以该模型、表热-弹性材料。用MP命令输入密度（DENS）（温度无关）。用MPTEMP和MPDATA输入弹性模量（Exx）、泊松比(NUXY)和热胀系数(ALPX)（这些特性和温度有关）。用TB，BISO，，NTEMP、TBTEMP以及TBDATA命令的1和2项输入屈服强度和切线模量。屈服强度和切线模量必须相对于同一温度定义，在MPTEMP命令中输入。MP，DENSMPTEMP，1，，，…，MPDATA，EX，，1，，，…，MPDATA，NUXY，，1，，，…MPDATA，ALPX，，1，，，…，TB，BISO，，NTEMP（NTEMP可为2到6）TBTEMP，（第一个温度点）TBDATA，1，（屈服应力）TBDATA，2，（切线模量）TBTEMP，（第二个