ANSYS-LS-DYNA-显式动力学-培训手册-第三天

1、材料的定义材料的定义第第3-1章章Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-2材料的定义材料的定义本章目标本章目标本章内容覆盖了本章内容覆盖了ANSYS/LS-DYNA中材料模式的选择与定义中材料模式的选择与定义标题标题:现有的材料模型现有的材料模型ANSYSLS-DYNA材料的图形操作界面（材料的图形操作界面（GUI）输入材料数据输入材料数据线弹性材料线弹性材料非线性弹性材料非线性弹

2、性材料塑性材料塑性材料应变率无关各向同性塑性材料应变率无关各向同性塑性材料应变率相关各向同性塑性材料应变率相关各向同性塑性材料应变率相关各向异性塑性材料应变率相关各向异性塑性材料压力相关塑性材料压力相关塑性材料温度敏感塑性材料温度敏感塑性材料A. 状态方程模型状态方程模型Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-3材料的定义材料的定义 本章目标本章目标标题标题 (续续):空材料空材料

3、损伤模型损伤模型泡沫材料泡沫材料离散单元性质离散单元性质索单元性质索单元性质刚性材料刚性材料一般材料定义指导一般材料定义指导材料定义练习材料定义练习Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-4材料的定义材料的定义A. 可用的材料模型可用的材料模型ANSYS/LS-DYNA支持比支持比ANSYS隐式更大的材料库，因此，它几乎能模拟隐式更大的材料库，因此，它几乎能模拟任何实际问题。任何实际

4、问题。ANSYS/LS-DYNA材料库提供许多特性，其中包括材料库提供许多特性，其中包括 :考虑应变失效的应变率相关塑性材料模型考虑应变失效的应变率相关塑性材料模型温度相关和温度敏感塑性材料模型温度相关和温度敏感塑性材料模型状态方程和空材料模型（鸟撞分析等）状态方程和空材料模型（鸟撞分析等）为了使用方便，这些材料模型分为为了使用方便，这些材料模型分为5组组GUI菜单菜单:LinearNonlinearEquation of StateDiscrete Element PropertiesRigid MaterialsExplicit Dynamics with ANSYS/LS-DYNA 6.

5、0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-5材料的定义材料的定义B. ANSYS/LS-DYNA材料图形用户界面材料图形用户界面材料材料GUI目录树图结构与隐式目录树图结构与隐式ANSYS一致一致.IIIIIIIVVExplicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #001630

6、3-6材料的定义材料的定义. ANSYS/LS-DYNA材料图形用户界面材料图形用户界面现有的现有的ANSYS/LS-DYNA材料库包括材料库包括: Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-7材料的定义材料的定义 ANSYS/LS-DYNA材料图形用户界面材料图形用户界面材料输入方法简单不易出错材料输入方法简单不易出错图形用户界面禁止密度随温度变化的输入图形用户界面禁止密度随温度变化

7、的输入不再使用不再使用MPMOD和和MPUNDO命令命令图形用户界面能接受直接批处理输入方式图形用户界面能接受直接批处理输入方式 (w/o MPMOD)Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-8材料的定义材料的定义C. 输入材料数据输入材料数据ANSYS/LS-DYNA中大多数材料需要输入密度中大多数材料需要输入密度 (DENS), 杨氏模量杨氏模量(EX)和泊和泊松比松比 (NUX

8、Y or PRXY) , 这些定义都使用这些定义都使用MP命令。命令。在在ANSYS/LS-DYNA中，一些材料模型需要输入载荷曲线。这些曲线用来中，一些材料模型需要输入载荷曲线。这些曲线用来定义材料的两个变量的相关性，例如屈服应力随塑性应变的变化。定义材料的两个变量的相关性，例如屈服应力随塑性应变的变化。使用两组数组和使用两组数组和EDCURVE 定义载荷曲线定义载荷曲线:Preprocessor Material Props Curve Options .通常，应力应变数据是指通常，应力应变数据是指真应力真应力与与真应变格式真应变格式Explicit Dynamics with ANSYS

9、/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-9材料的定义材料的定义. 输入材料数据输入材料数据通过通过EDCURVE, LIST或或EDCURVE, PLOT命令检查应力应变数据命令检查应力应变数据EDCURVE, LISTCurve starts at Yield PointExplicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Trai

10、ning ManualMarch 7, 2002Inventory #0016303-10材料的定义材料的定义D. 线弹性材料线弹性材料线弹性线弹性材料中包括四种不同的材料模型材料中包括四种不同的材料模型:流体流体 : 充满流体的容器在冲击载荷下的弹性性质充满流体的容器在冲击载荷下的弹性性质各向同性各向同性 : 材料性质各个方向都相同材料性质各个方向都相同正交各向异性正交各向异性 : 3个相互正交对称平面上的性质不同个相互正交对称平面上的性质不同各向异性各向异性 : 材料中一点的性质与该点在材料中的位置无关材料中一点的性质与该点在材料中的位置无关线弹性材料没有塑性变形而完全由线弹性材料没有塑性

11、变形而完全由虎克定律虎克定律来定义来定义:流体流体 : EDMP, FLUID, MAT, K弹性流体模型要求输入弹性流体模型要求输入DENS（密度）、密度）、K（体积模量）体积模量）体积模量可以通过体积模量可以通过EDMP命令直接输入或者由命令直接输入或者由EX和和NUXY的值用下式自动计算的值用下式自动计算:)21 (3EKjijIc Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-11

12、材料的定义材料的定义. 线弹性材料线弹性材料各向同性各向同性 :大多数工程材料（例如钢铁）都是各向同性的大多数工程材料（例如钢铁）都是各向同性的通过通过DENS、EX和和NUXY定义即可定义即可正交各向异性正交各向异性 : EDMP, ORTHO正交各向异性材料通过正交各向异性材料通过9个独立的常数和个独立的常数和DENS定义定义横向各向同性材料（正交各向异性的一种特例）通过横向各向同性材料（正交各向异性的一种特例）通过5个独立的常数个独立的常数(EXX, EZZ, NUXY, NUXZ, GXY)和和 DENS定义定义正交各向异性材料的定义与正交各向异性材料的定义与EDLCS命令所定义的坐标

13、系编号有关命令所定义的坐标系编号有关: Preprocessor: Material Props Local CS Create Local CS各向异性各向异性 : EDMP, ORTHO 和和 TB, ANEL各向异性材料通过各向异性材料通过21个相互独立的参数和个相互独立的参数和DENS定义定义使用局部坐标系使用局部坐标系 (EDLCS) 和数据表和数据表 (TB,ANEL)Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002

14、Inventory #0016303-12材料的定义材料的定义E. 非线性弹性材料非线性弹性材料非线性弹性非线性弹性材料包括三种材料模型材料包括三种材料模型:Blatz-Ko模型模型: 可压缩泡沫材料（例如聚氨酯橡胶）可压缩泡沫材料（例如聚氨酯橡胶）Mooney-Rivlin模型模型 : 不可压橡胶材料不可压橡胶材料粘弹性材料粘弹性材料 : 玻璃和类玻璃材料玻璃和类玻璃材料非线性弹性材料能够经受大的可恢复弹性变形。所有的非线性弹性材料能够经受大的可恢复弹性变形。所有的超弹超弹材料材料(Blatz-Ko 和和 Mooney-Rivlin)应变是可逆的，但是应变是可逆的，但是粘弹性粘弹性材料粘性部

15、分的应变是不可逆材料粘性部分的应变是不可逆的，弹性应变部分是可逆的。的，弹性应变部分是可逆的。Blatz-Ko 超弹模型超弹模型 :Blatz-Ko 材料模型只用于压缩下的橡胶材料材料模型只用于压缩下的橡胶材料ANSYS/LS-DYNA自动设定泊松比自动设定泊松比 (NUXY)为为 0.463，因此因此仅仅需要输入需要输入DENS和和GXY材料响应通过应变能密度函数材料响应通过应变能密度函数W定义定义,是应变不变量和32 522332IIIGxyWExplicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DY

16、NA 6.0Training ManualMarch 7, 2002Inventory #0016303-13材料的定义材料的定义. 非线性弹性材料非线性弹性材料Mooney-Rivlin 超弹材料超弹材料 : TB, MOONEY, , , , TBOPT用于定义不可压缩橡胶的材料响应用于定义不可压缩橡胶的材料响应与与ANSYS隐式双参数模型基本相同隐式双参数模型基本相同需要输入需要输入DENS, NUXY, 和常数和常数C10 和和 C01为保证不可压缩性质为保证不可压缩性质, NUXY 必须在必须在 0.49和和0.50间间取值取值Mooney-Rivlin系数可以直接输入系数可以直接输

17、入(TBOPT=0)或者通过或者通过LS-DYNA对测试数据的计对测试数据的计算得到算得到 (TBOPT=2)。对后一种情况，举例如下对后一种情况，举例如下: TB, MOONEY, 1, , , 2 ! 材料材料1的计算数据的计算数据 TBDATA, 1, 0.0, 0.0 ! C10 和和 C01 设为零设为零 TBDATA, 3, L0 , w, t ! 样品的原始长度、宽度和厚度样品的原始长度、宽度和厚度 TBDATA, 6, LCID ! 长度变化与力载荷曲线数据长度变化与力载荷曲线数据 注注: 如果如果L0 , w和和t都等于都等于1.0, 那么那么LCID =工程应力工程应力vs

18、. 应变应变 材料响应通过应变能密度函数材料响应通过应变能密度函数W定义：定义： 1 , 2 和和 3 是应变不变量，是应变不变量，K 是体积模量。是体积模量。2320111015 .3C3CWExplicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-14材料的定义材料的定义. 非线性弹性材料非线性弹性材料粘弹性材料粘弹性材料 : TB, EVISC这类材料具有这类材料具有弹性弹性（变形可恢复）部分和

19、（变形可恢复）部分和粘性粘性（变形不可恢复）部分（变形不可恢复）部分用于模拟类玻璃材料和火箭固体燃料等用于模拟类玻璃材料和火箭固体燃料等剪切关系由下式表达剪切关系由下式表达:除了密度除了密度 DENS, 还需输入以下参数还需输入以下参数: Go = 短期短期 (初始初始) 弹性剪切模量弹性剪切模量 G = 长期长期 (无限无限)弹性剪切模量弹性剪切模量 K = 弹性体积模量弹性体积模量 1/ = 衰减常数衰减常数 toe)GG(GtGExplicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0T

20、raining ManualMarch 7, 2002Inventory #0016303-15材料的定义材料的定义F. 塑性材料塑性材料塑性材料塑性材料模型包含模型包含ANSYS/LS-DYNA中大多数非线性中大多数非线性非弹性非弹性材料。要根据所材料。要根据所分析材料的类型、应用领域和材料常数的可获取性来选择某个特定塑性模型分析材料的类型、应用领域和材料常数的可获取性来选择某个特定塑性模型。塑性模型可以分为五类塑性模型可以分为五类:类别类别 1：应变率和塑性无关的各向同性材料应变率和塑性无关的各向同性材料类别类别 2: 应变率和塑性相关的各向同性材料应变率和塑性相关的各向同性材料类别类别

21、3: 应变率和塑性无关的各向异性材料应变率和塑性无关的各向异性材料类别类别 4: 压力相关塑性压力相关塑性类别类别 5: 温度敏感塑性温度敏感塑性为分析材料选择正确的类别为分析材料选择正确的类别非常重要非常重要， 在某类别内选择特定的模型相比就不在某类别内选择特定的模型相比就不那么重要了，这通常取决于材料数据的可获取性。那么重要了，这通常取决于材料数据的可获取性。Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory

22、 #0016303-16材料的定义材料的定义. 塑性材料塑性材料大多数非线性有限元分析精确性的关键在于材料常数的质量。为了得到最好大多数非线性有限元分析精确性的关键在于材料常数的质量。为了得到最好的结果，应该从材料供应者那里得到材料常数或者去做材料特性分析。的结果，应该从材料供应者那里得到材料常数或者去做材料特性分析。一些塑性模型需要输入附加的状态方程一些塑性模型需要输入附加的状态方程 (EOS)。这些方程将会在介绍完所有这些方程将会在介绍完所有的塑性模型后进行详细的讨论。的塑性模型后进行详细的讨论。Explicit Dynamics with ANSYS/LS-DYNA 6.0Explici

23、t Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-17材料的定义材料的定义G. 应变率无关各向同性塑性材料应变率无关各向同性塑性材料类别类别1: 应变率和塑性无关的应变率和塑性无关的 各向同性材料各向同性材料1a: 经典的经典的双线性随动硬化双线性随动硬化 (BKIN)1b: 经典的经典的双线性等向硬化双线性等向硬化 (BISO)BKIN 和和 BISO : TB, BKIN 和和 TB, BISO这些模型通常用于整个成形持续时间相对较长的过程（例如金属板冲压），这些模型通常用于整个成

24、形持续时间相对较长的过程（例如金属板冲压），以及大多数工程材料中（钢、铝、铸铁等）。以及大多数工程材料中（钢、铝、铸铁等）。两种模型都有两个斜率，弹性模量两种模型都有两个斜率，弹性模量(EX)和切向模量和切向模量(ETAN) 来表达材料的应力来表达材料的应力应变行为。应变行为。 BKIN和和BISO模型所需输入的参数相同模型所需输入的参数相同 : MP命令定义命令定义DENS, EX和和NUXYTB和和TBDATA命令定义屈服应力和切向模量命令定义屈服应力和切向模量Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with AN

25、SYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-18材料的定义材料的定义.应变率无关各向同性塑性材料应变率无关各向同性塑性材料BKIN和和BISO模型的差别仅在于硬化假设不同。模型的差别仅在于硬化假设不同。 随动硬化随动硬化假设二次屈服发生在假设二次屈服发生在2 y 等向硬化等向硬化假设二次屈服发生在假设二次屈服发生在2 maxBISO模型允许有模型允许有温度相关温度相关:预加热载荷分析和热瞬态分析预加热载荷分析和热瞬态分析 (见见第第2-1章章)指定六个温度条件下的材料数据指定六个温度条件下的材料数据使用使用MPTE

26、MP, MPDATA, TBTEMP和和TBDATA设定高屈服应力模拟热弹性材料设定高屈服应力模拟热弹性材料Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-19类别类别2:应变率和塑性相关的应变率和塑性相关的 各向同性材料各向同性材料2a: 塑性随动塑性随动: 带有失效应变的带有失效应变的Cowper-Symonds模型模型2b:幂率硬化幂率硬化: 带有强度和硬化系数的带有强度和硬化系数的C

27、owper-Symonds模型模型2c: 分段线性分段线性: 带有多线性曲线和失效应变的带有多线性曲线和失效应变的Cowper-Symonds模型模型2d: 率相关率相关: 通过载荷曲线和失效应力定义应变率通过载荷曲线和失效应力定义应变率2e: 应变率敏感应变率敏感: 超塑性成形的超塑性成形的Ramburgh-Osgood模型模型模型模型2a - 2d可用于一般材料和各向同性材料的塑性成形分析可用于一般材料和各向同性材料的塑性成形分析模型模型2a - 2c利用利用Cowper-Symonds 模型，模型的屈服应力与应变率因子有关模型，模型的屈服应力与应变率因子有关:P1C1 材料的定义材料的定

28、义H.应变率相关各向同性塑性材料应变率相关各向同性塑性材料C和和P是是Cowper-Symonds应变率参数。应变率参数。Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-20塑性随动塑性随动 : TB, PLAW, , , , 1双线性硬化塑性双线性硬化塑性 ( y 和和 ETAN)硬化参数硬化参数 在在 0 (kinematic) 和和 1 (isotropic)之间之间失效应变决定删除失

29、效单元失效应变决定删除失效单元屈服应力屈服应力:材料的定义材料的定义.应变率相关各向同性塑性材料应变率相关各向同性塑性材料其中其中 0 是初始屈服应力是初始屈服应力, peff 是有效塑性应变是有效塑性应变,Ep 是塑性硬化模量，由此公式确定是塑性硬化模量，由此公式确定 :effo1EC1pPPy tantanPEEEEEExplicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-21材料的定义材料的定

30、义.应变率相关各向同性塑性材料应变率相关各向同性塑性材料幂率硬化幂率硬化 : TB, PLAW, , , , 4带有双线性等向强化的塑性行为带有双线性等向强化的塑性行为强度系数强度系数k和硬化系数和硬化系数n定义的幂率硬化定义的幂率硬化屈服应力屈服应力:npePyeff1kC1 其中其中 e 是弹性应变。是弹性应变。Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-22材料的定义材料的定义.应

31、变率相关各向同性塑性材料应变率相关各向同性塑性材料分段线性分段线性 : TB, PLAW, , , , 8模型在求解时非常有效，通常用于碰撞分析模型在求解时非常有效，通常用于碰撞分析与与ANSYS 隐式中的隐式中的TB, MISO模型类似模型类似用有效真应力与有效真应变载荷曲线定义应力应变行为用有效真应力与有效真应变载荷曲线定义应力应变行为输入失效应变，以确定需要删除的单元输入失效应变，以确定需要删除的单元屈服面由屈服面由 Cowper-Symonds 模型针对率相关进行缩放模型针对率相关进行缩放Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit

32、Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-23材料的定义材料的定义.应变率相关各向同性塑性材料应变率相关各向同性塑性材料率相关率相关 : TB, PLAW, , , , 5这是最普通的应变率相关的塑性模型，因为弹性模量这是最普通的应变率相关的塑性模型，因为弹性模量(E) 、屈服应力、屈服应力( y) 、切向模量切向模量(ETAN)和失效应力和失效应力( FAIL)均可以应变的函数形式输入。均可以应变的函数形式输入。给定塑性应变率后的屈服应力如下定义给定塑性应变率后的屈服应力如下定义

33、:effpipioiyE其中其中LCID 1 = defines y as a function of LCID 2 = defines E as a function ofLCID 3 = defines ETAN as function ofLCID 4 = defines effective von Mises stress at failure as a function ofeff eff eff eff Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training M

34、anualMarch 7, 2002Inventory #0016303-24材料的定义材料的定义.应变率相关各向同性塑性材料应变率相关各向同性塑性材料nmyk 应变率敏感应变率敏感 : TB, PLAW, , , , 2是专门用于超塑性成形的特定模型。是专门用于超塑性成形的特定模型。Ramburgh-Osgood本构关系的屈服应力如下本构关系的屈服应力如下 :其中其中k 是材料系数是材料系数, m 是硬化系数是硬化系数, n 是应变率参数是应变率参数, 是应变率。是应变率。 Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics w

35、ith ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-25材料的定义材料的定义I.应变率相关各向异性塑性材料应变率相关各向异性塑性材料类别类别 3:应变率相关各向异性塑性材料应变率相关各向异性塑性材料3a: 横向各向异性横向各向异性: 应变率相关应变率相关Hills屈服准则屈服准则3b: 3参数参数Barlat: 用于铝板成形等正交各向异性模型用于铝板成形等正交各向异性模型3c: Barlat各向异性各向异性: 三维连续体成形各向异性模型三维连续体成形各向异性模型3d: 横向各向异性成形极限图硬化横向各向异性成形

36、极限图硬化模型模型3a用于模拟常规各向异性材料的高应变率成形过程用于模拟常规各向异性材料的高应变率成形过程模型模型3b和和3c由由ALCOA开发专用于金属铝的加工成形开发专用于金属铝的加工成形模型模型3d专用于钣金成形专用于钣金成形Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-26材料的定义材料的定义.应变率相关各向异性塑性材料应变率相关各向异性塑性材料横向各向异性横向各向异性 : TB,

37、 PLAW, , , , 7通常用于各向异性的钣金成形通常用于各向异性的钣金成形可以通过定义载荷曲线参数来表达有效屈服应力和有效塑性应变的关系可以通过定义载荷曲线参数来表达有效屈服应力和有效塑性应变的关系屈服应力屈服应力:2122211222211112212RRRRyp33p22R 各向异性硬化参数各向异性硬化参数R用面内塑性应变率与面外塑性应变率的比值来用面内塑性应变率与面外塑性应变率的比值来定义定义:Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training Manual

38、March 7, 2002Inventory #0016303-27材料的定义材料的定义.应变率相关各向异性塑性材料应变率相关各向异性塑性材料3参数参数Barlat : TB, PLAW, , , , 3开发用于开发用于铝板铝板在平面应力状态下的钣金成形在平面应力状态下的钣金成形对于线性硬化准则对于线性硬化准则, 输入输入 y 和和 ETAN 对于指数硬化准则对于指数硬化准则, 输入输入 n 和和 m 推荐的推荐的Barlat指数指数:体心立方晶格（体心立方晶格（BCC）金属取金属取m=6，面心立方晶格（面心立方晶格（FCC）金属取金属取 m=8正交各向异性正交各向异性Lankford系数用于

39、长度与厚度比系数用于长度与厚度比正交各向异性材料的坐标系通过正交各向异性材料的坐标系通过EDMP命令输入命令输入Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-28材料的定义材料的定义.应变率相关各向异性塑性材料应变率相关各向异性塑性材料Barlat各向异性各向异性 : TB, PLAW, , , , 6用于三维连续材料的钣金成形过程中的材料模型，特别是铝用于三维连续材料的钣金成形过程中的材

40、料模型，特别是铝大多用于实体材料（例如非板材）大多用于实体材料（例如非板材）实验确定实验确定6个各向异性参数个各向异性参数: a,b,c,f,g,h推荐的推荐的Barlat指数指数:体心立方晶格体心立方晶格金属取金属取m=6，面心立方晶格面心立方晶格金属取金属取 m=8屈服应力屈服应力: y=k( o+ p)n其中其中 o 和和 p 是初始屈服应变和塑性应变是初始屈服应变和塑性应变Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002

41、Inventory #0016303-29材料的定义材料的定义.应变率相关各向异性塑性材料应变率相关各向异性塑性材料横向各向异性横向各向异性 FLD : TB, PLAW, , , , 10用于模拟横向各向异性金属的钣金成形用于模拟横向各向异性金属的钣金成形仅用于壳单元仅用于壳单元屈服行为可以通过屈服行为可以通过 y和和ETAN或有效应力和塑性应变曲线来定义或有效应力和塑性应变曲线来定义成形极限图也可以用载荷曲线输入，来计算最大应变率成形极限图也可以用载荷曲线输入，来计算最大应变率Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics

42、 with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-30tthEEEEE应力应变行为也可以通过多达应力应变行为也可以通过多达16个数据点来定义。个数据点来定义。线性多项式状态方程也必须给定线性多项式状态方程也必须给定 (EOSOPT)。phoyE材料的定义材料的定义J.压力相关塑性材料压力相关塑性材料类别类别4: 压力相关塑性压力相关塑性4a: Elastic-Plastic Hydrodynamic: 在大应变下会失效的材料在大应变下会失效的材料4b: Geological Cap Model: 地质盖帽

43、模型地质盖帽模型Elastic-Plastic Hydrodynamic : TB, PLAW, , , , 9, EOSOPT如果没有定义有效真应力和真应变，就假定为等向强化，而且必须指定如果没有定义有效真应力和真应变，就假定为等向强化，而且必须指定Sy和和ETAN来定义屈服强度来定义屈服强度 :其中塑性硬化模量其中塑性硬化模量 Eh通过通过E 和和 ETAN定义定义Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Invent

44、ory #0016303-31材料的定义材料的定义.压力相关塑性材料压力相关塑性材料Geological Cap Model : TB, GCAP无粘滞性双不变量地质盖帽材料模型无粘滞性双不变量地质盖帽材料模型地质力学问题或混凝土类材料地质力学问题或混凝土类材料双不变量盖帽理论扩展至包含随动强化双不变量盖帽理论扩展至包含随动强化采用剪切模量代替弹性模量采用剪切模量代替弹性模量模型的详细信息请参考模型的详细信息请参考LS-DYNA理论手册理论手册（在与（在与ANSYS6.0理论手册的理论手册的同一张光碟上）同一张光碟上）。Explicit Dynamics with ANSYS/LS-DYNA

45、6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-32材料的定义材料的定义K.温度敏感塑性材料温度敏感塑性材料类别类别5: ANSYS/LS-DYNA中有中有四种塑性模型四种塑性模型来来描述温度效应描述温度效应. 大部分模型需大部分模型需要定义附加状态方程要定义附加状态方程(EOS) :5a: Bamman: 使用内部状态方程变量的复杂材料模型使用内部状态方程变量的复杂材料模型5b: Johnson-Cook: 高应变率和温度相关问题高应变率和温度相关问题5c: Z

46、erilli-Armstrong: 高速冲击和某些金属成形过程高速冲击和某些金属成形过程5d: Steinberg: 极高应变率条件下极高应变率条件下(105)材料成形的模拟材料成形的模拟Bamman : TB, EOS, , , , 4用于金属成形的复杂应变率相关材料模型用于金属成形的复杂应变率相关材料模型内部状态方程变量通过常数内部状态方程变量通过常数Ai直接输入模型，无需定义附加直接输入模型，无需定义附加状态方程状态方程。 Bamman模型需要输入流动应力参数模型需要输入流动应力参数CiExplicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dy

47、namics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-33A, B, C, m, 和和n 是用实验方法确定的常数，是用实验方法确定的常数， p是有效塑性应变是有效塑性应变也需要输入有效塑性应变率，由下式给出也需要输入有效塑性应变率，由下式给出:m_npyT1lnC1BAroommeltroomTTT-TT :is etemperatur homologous the and op材料的定义材料的定义.温度敏感塑性材料温度敏感塑性材料Johnson-Cook : TB, EOS, , , , 1,

48、EOSOPTJohnson-Cook模型主要用于高应变率过程，如伴随较大温升的机加工模型主要用于高应变率过程，如伴随较大温升的机加工模型的开发源自弹道学模型的开发源自弹道学需要输入需要输入 DENS, EX, 和和 NUXY屈服应力屈服应力:Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-34温度计算需要输入比热、熔点和室温温度计算需要输入比热、熔点和室温通过下式中的失效常数通过下式中的失效

49、常数D1-D5 ，失效应变可以并入到模型中失效应变可以并入到模型中:其中其中 Johnson-Cook参数输入后，要输入状态方程常数，可以是线性多项式或参数输入后，要输入状态方程常数，可以是线性多项式或Gruneisen模型模型(EOSOPT 后面讨论后面讨论)。TD1lnD1expDD54D21f3材料的定义材料的定义.温度敏感塑性材料温度敏感塑性材料.压力与有效应力之比 effExplicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002

50、Inventory #0016303-35材料的定义材料的定义.温度敏感塑性材料温度敏感塑性材料Zerilli-Armstrong : TB, EOS, , , , 3, EOSOPT Zerilli-Armstrong模型用于金属成形过程和高速冲击等这些应力与应变、模型用于金属成形过程和高速冲击等这些应力与应变、应变率及温度相关的应用领域应变率及温度相关的应用领域 。这种模型必须使用状态方程。这种模型必须使用状态方程。 Zerilli-Armstrong模型需输入流动应力模型需输入流动应力 (Ci), 温度温度 (Bi) 和热容和热容(Gi)系数。系数。)293()()(:)293()()(

51、:65)ln(215)ln(5 .214343TCCeCCTCeCCnpTeCCyTeCCpy体心立方晶格金属面心立方晶格金属Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-36材料的定义材料的定义.温度敏感塑性材料温度敏感塑性材料Steinberg : TB, EOS, , , , 5, EOSOPT带失效的高应变率带失效的高应变率(105)的固体变形的固体变形非常适合于加工过程和高速冲击

52、分析非常适合于加工过程和高速冲击分析屈服强度是温度和压力的函数屈服强度是温度和压力的函数状态方程状态方程(EOS)决定压力决定压力Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-37材料的定义材料的定义L. 状态方程状态方程 (EOS)ANSYS/LS-DYNA中有三种不同类型的中有三种不同类型的 EOS :线性多项式线性多项式 : EOS 是内能的线性方程是内能的线性方程Gruneisen

53、 : EOS 通过两种方式定义压力体积关系通过两种方式定义压力体积关系列表格式（列表格式（ Tabulated ）: 简化的简化的EOS ，也是内能的线性方程也是内能的线性方程线性多项式线性多项式 : EOSOPT = 1EOS 是内能的线性方程是内能的线性方程压力由压力由 和线性系数和线性系数Ci决定决定:P = C0 + C1 + C2 2 + C3 3 + (C4 + C5 +C6 2) E其中其中 = / 0 1 ， 和和 0 是当前和初始密度是当前和初始密度Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with A

54、NSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-38EaSSSaCpooo2332212211112211材料的定义材料的定义.状态方程状态方程 (EOS)Gruneisen : EOSOPT = 2压力体积关系与材料是压缩还是膨胀有关。带有立方冲击速度压力体积关系与材料是压缩还是膨胀有关。带有立方冲击速度-质点速度的质点速度的Gruneisen状态方程通过状态方程通过 和和Gruneisen系数系数C, a, S1, S2, S3及及 0定义压力，定义压力，对压缩材料：对压缩材料：列表格式（列表格式（Tabulat

55、ed） : EOSOPT = 3EOS是内能的线性方程是内能的线性方程. 压力为压力为: P = Ci ( vi) + Ti ( vi) E 其中其中 Ci 和和 Ti 分别是体积压力和温度分别是体积压力和温度, vi 是体积应变是体积应变Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-39材料的定义材料的定义M. 空材料空材料空材料空材料 : TB, EOS, , , , 2, EOSOP

56、T空材料允许在不计算偏应力的条件下考虑状态方程空材料允许在不计算偏应力的条件下考虑状态方程空材料大多用于鸟撞分析空材料大多用于鸟撞分析需要输入需要输入密度密度DENS和和截止压力截止压力。对空的梁和壳需要输入对空的梁和壳需要输入EX 和和 NUXY也可定义拉压下的粘性和侵蚀（可选）也可定义拉压下的粘性和侵蚀（可选）EOSOPT 需设定为需设定为1、2或或3来指定参考的状态方程来指定参考的状态方程Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch

57、7, 2002Inventory #0016303-40材料的定义材料的定义N. 损伤模型损伤模型ANSYS/LS-DYNA有两种损伤模型有两种损伤模型:混凝土损伤混凝土损伤: 冲击载荷下的钢筋混凝土冲击载荷下的钢筋混凝土复合材料损伤复合材料损伤: 用于吸能的复合材料失效分析用于吸能的复合材料失效分析混凝土损伤混凝土损伤: TB, CONCR, , , , 2仅用于仅用于SOLID164单元单元, 需要输入泊松比、密度、混凝土和加强钢筋常数需要输入泊松比、密度、混凝土和加强钢筋常数必须使用列表格式状态方程必须使用列表格式状态方程复合材料损伤复合材料损伤: TB, COMP必须输入每个方向上的弹

58、性模量、剪切模量和泊松比（无默认值）必须输入每个方向上的弹性模量、剪切模量和泊松比（无默认值）可以输入剪切、纵向拉伸、横向拉伸、横向压缩强度来定义失效可以输入剪切、纵向拉伸、横向拉伸、横向压缩强度来定义失效压缩失效需要定义材料的体积模量压缩失效需要定义材料的体积模量Explicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-41材料的定义材料的定义O. 泡沫材料泡沫材料ANSYS/LS-DYNA有有5种

59、不同的泡沫材料模型种不同的泡沫材料模型:闭室泡沫闭室泡沫: 低密度泡沫材料低密度泡沫材料 (例如聚氨酯泡沫例如聚氨酯泡沫)低密度泡沫低密度泡沫: 高度可压缩性材料高度可压缩性材料 (例如座垫例如座垫)粘性泡沫粘性泡沫: 用在碰撞跌落模拟中的能量吸收泡沫用在碰撞跌落模拟中的能量吸收泡沫可压碎泡沫可压碎泡沫: 永久性压碎材料永久性压碎材料 (例如聚苯乙烯例如聚苯乙烯)蜂窝材料蜂窝材料: 正交各向异性压碎泡沫正交各向异性压碎泡沫特定模型的选择取决于所分析材料的类型特定模型的选择取决于所分析材料的类型ANSYS/LS-DYNA所有泡沫材料模型主要用于汽车碰撞分析所有泡沫材料模型主要用于汽车碰撞分析Ex

60、plicit Dynamics with ANSYS/LS-DYNA 6.0Explicit Dynamics with ANSYS/LS-DYNA 6.0Training ManualMarch 7, 2002Inventory #0016303-42材料的定义材料的定义. 泡沫材料泡沫材料闭室泡沫闭室泡沫: TB, FOAM, , , , 1开发用于开发用于低密度聚氨酯低密度聚氨酯 (常用于船用集装箱及汽车设计中的冲击限制物建模常用于船用集装箱及汽车设计中的冲击限制物建模)要输入要输入DENS, EX, 初始泡沫压力初始泡沫压力Po, 泡沫聚合密度泡沫聚合密度比率比率 包括封闭气压效应包括