《LSDYNA培训手册》PPT课件

1、材料的定义第3-1章材料的定义本章目标本章内容覆盖了ANSYS/LS-DYNA中材料模式的选择与定义标题:现有的材料模型ANSYSLS-DYNA材料的图形操作界面（GUI）输入材料数据线弹性材料非线性弹性材料塑性材料应变率无关各向同性塑性材料应变率相关各向同性塑性材料应变率相关各向异性塑性材料压力相关塑性材料温度敏感塑性材料状态方程模型培训ppt材料的定义 本章目标标题 (续):空材料 损伤模型泡沫材料离散单元性质索单元性质刚性材料一般材料定义指导材料定义练习培训ppt材料的定义A. 可用的材料模型ANSYS/LS-DYNA支持比ANSYS隐式更大的材料库，因此，它几乎能模拟任何实际问题。AN

2、SYS/LS-DYNA材料库提供许多特性，其中包括 :考虑应变失效的应变率相关塑性材料模型温度相关和温度敏感塑性材料模型状态方程和空材料模型（鸟撞分析等）为了使用方便，这些材料模型分为5组GUI菜单:LinearNonlinearEquation of StateDiscrete Element PropertiesRigid Materials培训ppt材料的定义B. ANSYS/LS-DYNA材料图形用户界面材料GUI目录树图结构与隐式ANSYS一致.IIIIIIIVV培训ppt材料的定义. ANSYS/LS-DYNA材料图形用户界面现有的ANSYS/LS-DYNA材料库包括: 培训ppt

3、材料的定义 ANSYS/LS-DYNA材料图形用户界面材料输入方法简单不易出错图形用户界面禁止密度随温度变化的输入不再使用MPMOD和MPUNDO命令图形用户界面能接受直接批处理输入方式 (w/o MPMOD)培训ppt材料的定义C. 输入材料数据ANSYS/LS-DYNA中大多数材料需要输入密度 (DENS), 杨氏模量(EX)和泊松比 (NUXY or PRXY) , 这些定义都使用MP命令。在ANSYS/LS-DYNA中，一些材料模型需要输入载荷曲线。这些曲线用来定义材料的两个变量的相关性，例如屈服应力随塑性应变的变化。使用两组数组和EDCURVE 定义载荷曲线:Preprocessor

4、 Material Props Curve Options .通常，应力应变数据是指真应力与真应变格式培训ppt材料的定义. 输入材料数据通过EDCURVE, LIST或EDCURVE, PLOT命令检查应力应变数据EDCURVE, LISTCurve starts at Yield Point培训ppt材料的定义D. 线弹性材料线弹性材料中包括四种不同的材料模型:流体 : 充满流体的容器在冲击载荷下的弹性性质各向同性 : 材料性质各个方向都相同正交各向异性 : 3个相互正交对称平面上的性质不同各向异性 : 材料中一点的性质与该点在材料中的位置无关线弹性材料没有塑性变形而完全由虎克定律来定义:

5、流体 : EDMP, FLUID, MAT, K弹性流体模型要求输入DENS（密度）、K（体积模量）体积模量可以通过EDMP命令直接输入或者由EX和NUXY的值用下式自动计算:培训ppt材料的定义. 线弹性材料各向同性 :大多数工程材料（例如钢铁）都是各向同性的通过DENS、EX和NUXY定义即可正交各向异性 : EDMP, ORTHO正交各向异性材料通过9个独立的常数和DENS定义横向各向同性材料（正交各向异性的一种特例）通过5个独立的常数(EXX, EZZ, NUXY, NUXZ, GXY)和 DENS定义正交各向异性材料的定义与EDLCS命令所定义的坐标系编号有关:Preprocesso

6、r: Material Props Local CS Create Local CS各向异性 : EDMP, ORTHO 和 TB, ANEL各向异性材料通过21个相互独立的参数和DENS定义使用局部坐标系 (EDLCS) 和数据表 (TB,ANEL)培训ppt材料的定义E. 非线性弹性材料非线性弹性材料包括三种材料模型:Blatz-Ko模型: 可压缩泡沫材料（例如聚氨酯橡胶）Mooney-Rivlin模型 : 不可压橡胶材料粘弹性材料 : 玻璃和类玻璃材料非线性弹性材料能够经受大的可恢复弹性变形。所有的超弹材料(Blatz-Ko 和 Mooney-Rivlin)应变是可逆的，但是粘弹性材料粘

7、性部分的应变是不可逆的，弹性应变部分是可逆的。Blatz-Ko 超弹模型 :Blatz-Ko 材料模型只用于压缩下的橡胶材料ANSYS/LS-DYNA自动设定泊松比 (NUXY)为 0.463，因此仅需要输入DENS和GXY材料响应通过应变能密度函数W定义,培训ppt材料的定义. 非线性弹性材料Mooney-Rivlin 超弹材料 : TB, MOONEY, , , , TBOPT用于定义不可压缩橡胶的材料响应与ANSYS隐式双参数模型基本相同需要输入DENS, NUXY, 和常数C10 和 C01为保证不可压缩性质, NUXY 必须在 0.49和0.50间取值Mooney-Rivlin系数可

8、以直接输入(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. 应变材料响应通过应变能密度函数W定义：1 , 2 和 3 是应变不变量，K 是体积模量。培训ppt材料的定义. 非线性弹性材料粘弹性材料 :

9、TB, EVISC这类材料具有弹性（变形可恢复）部分和粘性（变形不可恢复）部分用于模拟类玻璃材料和火箭固体燃料等剪切关系由下式表达:除了密度 DENS, 还需输入以下参数:Go = 短期 (初始) 弹性剪切模量G = 长期 (无限)弹性剪切模量K = 弹性体积模量1/= 衰减常数培训ppt材料的定义F. 塑性材料塑性材料模型包含ANSYS/LS-DYNA中大多数非线性非弹性材料。要根据所分析材料的类型、应用领域和材料常数的可获取性来选择某个特定塑性模型。塑性模型可以分为五类:类别 1：应变率和塑性无关的各向同性材料类别 2: 应变率和塑性相关的各向同性材料类别 3: 应变率和塑性无关的各向异性

10、材料类别 4: 压力相关塑性类别 5: 温度敏感塑性为分析材料选择正确的类别非常重要， 在某类别内选择特定的模型相比就不那么重要了，这通常取决于材料数据的可获取性。培训ppt材料的定义. 塑性材料大多数非线性有限元分析精确性的关键在于材料常数的质量。为了得到最好的结果，应该从材料供应者那里得到材料常数或者去做材料特性分析。一些塑性模型需要输入附加的状态方程 (EOS)。这些方程将会在介绍完所有的塑性模型后进行详细的讨论。培训ppt材料的定义G. 应变率无关各向同性塑性材料类别1: 应变率和塑性无关的 各向同性材料1a: 经典的双线性随动硬化 (BKIN)1b: 经典的双线性等向硬化 (BISO

11、)BKIN 和 BISO : TB, BKIN 和 TB, BISO这些模型通常用于整个成形持续时间相对较长的过程（例如金属板冲压），以及大多数工程材料中（钢、铝、铸铁等）。两种模型都有两个斜率，弹性模量(EX)和切向模量(ETAN) 来表达材料的应力应变行为。 BKIN和BISO模型所需输入的参数相同 : MP命令定义DENS, EX和NUXYTB和TBDATA命令定义屈服应力和切向模量培训ppt材料的定义.应变率无关各向同性塑性材料BKIN和BISO模型的差别仅在于硬化假设不同。 随动硬化假设二次屈服发生在2y 等向硬化假设二次屈服发生在2maxBISO模型允许有温度相关:预加热载荷分析和

12、热瞬态分析 (见第2-1章)指定六个温度条件下的材料数据使用MPTEMP, MPDATA, TBTEMP和TBDATA设定高屈服应力模拟热弹性材料培训ppt类别2:应变率和塑性相关的 各向同性材料2a: 塑性随动: 带有失效应变的Cowper-Symonds模型2b:幂率硬化: 带有强度和硬化系数的Cowper-Symonds模型2c: 分段线性: 带有多线性曲线和失效应变的Cowper-Symonds模型2d: 率相关: 通过载荷曲线和失效应力定义应变率2e: 应变率敏感: 超塑性成形的Ramburgh-Osgood模型模型2a - 2d可用于一般材料和各向同性材料的塑性成形分析模型2a -

13、 2c利用Cowper-Symonds 模型，模型的屈服应力与应变率因子有关:材料的定义H.应变率相关各向同性塑性材料C和P是Cowper-Symonds应变率参数。培训ppt塑性随动 : TB, PLAW, , , , 1双线性硬化塑性 (y 和 ETAN)硬化参数 在 0 (kinematic) 和 1 (isotropic)之间失效应变决定删除失效单元屈服应力:材料的定义.应变率相关各向同性塑性材料其中0 是初始屈服应力,peff 是有效塑性应变,Ep 是塑性硬化模量，由此公式确定 :培训ppt材料的定义.应变率相关各向同性塑性材料幂率硬化 : TB, PLAW, , , , 4带有双线

14、性等向强化的塑性行为强度系数k和硬化系数n定义的幂率硬化屈服应力:其中e 是弹性应变。培训ppt材料的定义.应变率相关各向同性塑性材料分段线性 : TB, PLAW, , , , 8模型在求解时非常有效，通常用于碰撞分析与ANSYS 隐式中的TB, MISO模型类似用有效真应力与有效真应变载荷曲线定义应力应变行为输入失效应变，以确定需要删除的单元屈服面由 Cowper-Symonds 模型针对率相关进行缩放培训ppt材料的定义.应变率相关各向同性塑性材料率相关 : TB, PLAW, , , , 5这是最普通的应变率相关的塑性模型，因为弹性模量(E) 、屈服应力(y) 、切向模量(ETAN)和

15、失效应力(FAIL)均可以应变的函数形式输入。给定塑性应变率后的屈服应力如下定义:其中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 of培训ppt材料的定义.应变率相关各向同性塑性材料应变率敏感 : TB, PLAW, , , , 2是专门用于超塑性成形的特定模型。Ramburgh-Osgoo

16、d本构关系的屈服应力如下 :其中k 是材料系数, m 是硬化系数, n 是应变率参数, 是应变率。培训ppt材料的定义I.应变率相关各向异性塑性材料类别 3:应变率相关各向异性塑性材料3a: 横向各向异性: 应变率相关Hills屈服准则3b: 3参数Barlat: 用于铝板成形等正交各向异性模型3c: Barlat各向异性: 三维连续体成形各向异性模型3d: 横向各向异性成形极限图硬化模型3a用于模拟常规各向异性材料的高应变率成形过程模型3b和3c由ALCOA开发专用于金属铝的加工成形模型3d专用于钣金成形培训ppt材料的定义.应变率相关各向异性塑性材料横向各向异性 : TB, PLAW, ,

17、 , , 7通常用于各向异性的钣金成形可以通过定义载荷曲线参数来表达有效屈服应力和有效塑性应变的关系屈服应力:各向异性硬化参数R用面内塑性应变率与面外塑性应变率的比值来定义:培训ppt材料的定义.应变率相关各向异性塑性材料3参数Barlat : TB, PLAW, , , , 3开发用于铝板在平面应力状态下的钣金成形对于线性硬化准则, 输入 y 和 ETAN 对于指数硬化准则, 输入 n 和 m 推荐的Barlat指数:体心立方晶格（BCC）金属取m=6，面心立方晶格（FCC）金属取 m=8正交各向异性Lankford系数用于长度与厚度比正交各向异性材料的坐标系通过EDMP命令输入培训ppt材

18、料的定义.应变率相关各向异性塑性材料Barlat各向异性 : TB, PLAW, , , , 6用于三维连续材料的钣金成形过程中的材料模型，特别是铝大多用于实体材料（例如非板材）实验确定6个各向异性参数: a,b,c,f,g,h推荐的Barlat指数:体心立方晶格金属取m=6，面心立方晶格金属取 m=8屈服应力: y=k(o+ p)n其中 o 和 p 是初始屈服应变和塑性应变培训ppt材料的定义.应变率相关各向异性塑性材料横向各向异性 FLD : TB, PLAW, , , , 10用于模拟横向各向异性金属的钣金成形仅用于壳单元屈服行为可以通过sy和ETAN或有效应力和塑性应变曲线来定义成形极

19、限图也可以用载荷曲线输入，来计算最大应变率培训ppt应力应变行为也可以通过多达16个数据点来定义。线性多项式状态方程也必须给定 (EOSOPT)。材料的定义J.压力相关塑性材料类别4: 压力相关塑性4a: Elastic-Plastic Hydrodynamic: 在大应变下会失效的材料4b: Geological Cap Model: 地质盖帽模型Elastic-Plastic Hydrodynamic : TB, PLAW, , , , 9, EOSOPT如果没有定义有效真应力和真应变，就假定为等向强化，而且必须指定Sy和ETAN来定义屈服强度 :其中塑性硬化模量 Eh通过E 和 ETAN

20、定义培训ppt材料的定义.压力相关塑性材料Geological Cap Model : TB, GCAP无粘滞性双不变量地质盖帽材料模型地质力学问题或混凝土类材料双不变量盖帽理论扩展至包含随动强化采用剪切模量代替弹性模量模型的详细信息请参考LS-DYNA理论手册（在与ANSYS6.0理论手册的同一张光碟上）。培训ppt材料的定义K.温度敏感塑性材料类别5: ANSYS/LS-DYNA中有四种塑性模型来描述温度效应. 大部分模型需要定义附加状态方程(EOS) :5a: Bamman: 使用内部状态方程变量的复杂材料模型5b: Johnson-Cook: 高应变率和温度相关问题5c: Zerill

21、i-Armstrong: 高速冲击和某些金属成形过程5d: Steinberg: 极高应变率条件下(105)材料成形的模拟Bamman : TB, EOS, , , , 4用于金属成形的复杂应变率相关材料模型内部状态方程变量通过常数Ai直接输入模型，无需定义附加状态方程。 Bamman模型需要输入流动应力参数Ci培训pptA, B, C, m, 和n 是用实验方法确定的常数， p是有效塑性应变也需要输入有效塑性应变率，由下式给出:材料的定义.温度敏感塑性材料Johnson-Cook : TB, EOS, , , , 1, EOSOPTJohnson-Cook模型主要用于高应变率过程，如伴随较大

22、温升的机加工模型的开发源自弹道学需要输入 DENS, EX, 和 NUXY屈服应力:培训ppt温度计算需要输入比热、熔点和室温通过下式中的失效常数D1-D5 ，失效应变可以并入到模型中:其中 Johnson-Cook参数输入后，要输入状态方程常数，可以是线性多项式或Gruneisen模型(EOSOPT 后面讨论)。材料的定义.温度敏感塑性材料培训ppt材料的定义.温度敏感塑性材料Zerilli-Armstrong : TB, EOS, , , , 3, EOSOPTZerilli-Armstrong模型用于金属成形过程和高速冲击等这些应力与应变、应变率及温度相关的应用领域 。这种模型必须使用状

23、态方程。Zerilli-Armstrong模型需输入流动应力 (Ci), 温度 (Bi) 和热容(Gi)系数。培训ppt材料的定义.温度敏感塑性材料Steinberg : TB, EOS, , , , 5, EOSOPT带失效的高应变率(105)的固体变形非常适合于加工过程和高速冲击分析屈服强度是温度和压力的函数状态方程(EOS)决定压力培训ppt材料的定义L. 状态方程 (EOS)ANSYS/LS-DYNA中有三种不同类型的 EOS :线性多项式 : EOS 是内能的线性方程Gruneisen : EOS 通过两种方式定义压力体积关系列表格式（ Tabulated ）: 简化的EOS ，也是

24、内能的线性方程线性多项式 : EOSOPT = 1EOS 是内能的线性方程压力由和线性系数Ci决定:P = C0 + C1 + C22 + C33 + (C4 + C5 +C62) E其中 = /0 1 ， 和 0 是当前和初始密度培训ppt材料的定义.状态方程 (EOS)Gruneisen : EOSOPT = 2压力体积关系与材料是压缩还是膨胀有关。带有立方冲击速度-质点速度的Gruneisen状态方程通过和Gruneisen系数C, a, S1, S2, S3及0定义压力，对压缩材料：列表格式（Tabulated） : EOSOPT = 3EOS是内能的线性方程. 压力为:P = Ci

25、(vi) + Ti (vi) E其中 Ci 和 Ti 分别是体积压力和温度, vi 是体积应变培训ppt材料的定义M. 空材料空材料 : TB, EOS, , , , 2, EOSOPT空材料允许在不计算偏应力的条件下考虑状态方程空材料大多用于鸟撞分析需要输入密度DENS和截止压力。对空的梁和壳需要输入EX 和 NUXY也可定义拉压下的粘性和侵蚀（可选）EOSOPT 需设定为1、2或3来指定参考的状态方程培训ppt材料的定义N. 损伤模型ANSYS/LS-DYNA有两种损伤模型:混凝土损伤: 冲击载荷下的钢筋混凝土复合材料损伤: 用于吸能的复合材料失效分析混凝土损伤: TB, CONCR, ,

26、 , , 2仅用于SOLID164单元, 需要输入泊松比、密度、混凝土和加强钢筋常数必须使用列表格式状态方程复合材料损伤: TB, COMP必须输入每个方向上的弹性模量、剪切模量和泊松比（无默认值）可以输入剪切、纵向拉伸、横向拉伸、横向压缩强度来定义失效压缩失效需要定义材料的体积模量培训ppt材料的定义O. 泡沫材料ANSYS/LS-DYNA有5种不同的泡沫材料模型:闭室泡沫: 低密度泡沫材料 (例如聚氨酯泡沫)低密度泡沫: 高度可压缩性材料 (例如座垫)粘性泡沫: 用在碰撞跌落模拟中的能量吸收泡沫可压碎泡沫: 永久性压碎材料 (例如聚苯乙烯)蜂窝材料: 正交各向异性压碎泡沫特定模型的选择取决

27、于所分析材料的类型ANSYS/LS-DYNA所有泡沫材料模型主要用于汽车碰撞分析培训ppt材料的定义. 泡沫材料闭室泡沫: TB, FOAM, , , , 1开发用于低密度聚氨酯 (常用于船用集装箱及汽车设计中的冲击限制物建模)要输入DENS, EX, 初始泡沫压力Po, 泡沫聚合密度比率包括封闭气压效应NUXY 被近似地设为零屈服条件: y = a + b(1+c )其中 a, b和c 为实验确定的参数, = V/ Vo + o - 1其中 V 是当前体积, Vo 是初始体积, o 是初始体积应变培训ppt材料的定义. 泡沫材料低密度泡沫 : TB, FOAM, , , , 2主要用于模拟汽

28、车座垫需要输入密度DENS和弹性模量EX 采用载荷曲线标识号(LCID)输入应力应变行为对于压缩，模型用可能的能量耗散来假设材料的滞后行为对于拉伸，直到拉伸截止应力以前模型均表现为线性行为NUXY 近似的设为零使用滞后卸载时，如果衰减常数 = 0 ，重新加载将会沿卸载曲线进行如果非零, 初始载荷由1- e-t决定，可用粘性系数(0.05-0.5)模拟阻尼效应体积粘度如果标记为1，则被激活滞后卸载使用形状卸载因数，数值上小于1能量耗散输入的滞后卸载(HU)因数在0到1之间，如果HU=1 ，就没有能量耗散培训ppt材料的定义. 泡沫材料粘性泡沫 : TB, FOAM, , , , 3在碰撞模拟中吸

29、能泡沫常用于模拟吸能材料（例如假人）仅用于压缩载荷下的实体模型由平行的非线性弹簧和粘性阻尼构成需要有DENS, EX (初始杨氏模量E1) 和 NUXY弹性刚度E 定义为: E = E1 V-n1其中 n1 是弹性刚度幂率粘性系数V 定义为: V = V2 |1-V|n2其中 V2 = 初始粘度系数 n2 = 粘度系数幂率培训ppt材料的定义. 泡沫材料可压碎泡沫 : TB, FOAM, , , , 4用于可永久压碎的材料（如聚苯乙烯固体）可选择粘性阻尼和截止张力（撕裂）卸载认为是完全弹性过程拉伸按照完全弹塑性处理需DENS, EX, NUXY需输入定义屈服应力与体积应变的载荷曲线 体积应变

30、为: = 1 V其中 V 为当前体积和初始体积之比培训ppt材料的定义. 泡沫材料蜂窝材料 : TB, FOAM, , , , 5正交各向异性可压碎泡沫为侧面冲击缓冲器前端材料和航空结构而开发可以分别定义法向应力和剪切应力的非线性行为需要DENS, EX, NUXY和粘性系数需输入完全压实的蜂窝的屈服应力和体积需输入弹性模量和每个正交方向上的应力与相对体积或体积应变载荷曲线培训ppt材料的定义P. 离散单元性质COMBI165 离散单元要求有弹簧或阻尼性质弹簧需要有刚度系数或者力与变形曲线阻尼需要有阻尼系数或者力与变形速率曲线转动性质可以用来代替平动性质离散弹簧性质的六个选项:线弹性弹簧 (T

31、B, DISCRETE, , , , 0)通用非线性 (TB, DISCRETE, , , , 5)非线性弹性弹簧 (TB, DISCRETE, , , , 3)弹性塑料弹簧 (TB, DISCRETE, , , , 2)非弹性拉压 (TB, DISCRETE, , , , 7)Maxwell粘性弹簧 (TB, DISCRETE, , , , 6)离散阻尼性质的两个选项:线性粘性阻尼 (TB, DISCRETE, , , , 1)非线性粘性阻尼 (TB, DISCRETE, , , , 4)培训ppt材料的定义Q. 索单元性质LINK167 单元需要有索材料属性单元压缩时（松弛状态）不承担载荷

32、需要定义密度 (DENS)和弹性模量(EX)索链刚度:K = (E) (Area0) / (L0 Offset)或者通过指定工程应力-应变（载荷）曲线，而不考虑用杨氏模量来计算索链刚度培训ppt材料的定义R. 刚性材料声明一种材料是刚性的，从而这种材料构成的梁、壳和实体单元都成为刚体. 下述几点很重要:刚体材料的杨氏模量不能任意大。LS-DYNA 用杨氏模量计算接触罚刚度，而接触罚刚度决定了接触穿透。如果材料声明为刚性，那么任何属于这种材料的单元必须属于同一刚体。因此，分配材料属性时必须非常小心。正如前一章讨论的，最好考虑使用Parts，采用完全相同的单元类型、实常数集和材料来定义各种不同的P

33、arts.EDMP, RIGID命令不仅声明一种材料是刚性的，同时也约束了刚体的运 动特性培训ppt材料的定义S. 一般材料定义指导并非所有材料都支持每一种单元类型。查阅单元手册来决定采用哪种模型。一些模型需要说明状态方程来完成材料的定义。对每种材料模型，并非所有的常数和选项都需要输入。例如带有Cowper-Symonds常数的应变率相关塑性材料的常数设为零，就可以用作率无关模型。这可以用作处理允许的失效应变定义材料性质时单位要一致。不正确的单位不仅影响材料响应，也会影响接触刚度。单位的错误甚至可能导致模型无法求解。不要低估精确材料数据的重要性。多花些精力以获得精确的数据。培训ppt材料的定义

34、T. 材料定义练习本练习讨论的问题是:练习3-1. 弹簧质量阻尼系统响应模拟培训ppt接触面第 3-2章接触面本章目标本章是关于接触面的内容. 介绍了不同的算法及推荐使用的算法。主题:接触面概述单面接触算法点对面接触算法面对面接触算法一般接触类自动接触类刚性接触类固连接触类固连失效接触类培训ppt接触面 本章目标主题 (续):侵蚀接触类边接触类拉延筋接触类成形接触类2-D 接触类定义基本接触接触的列表和显示接触的删除/撤消/激活高级接触控制一般接触指导接触面的练习培训ppt 接触表面A. 接触表面概述在ANSYS/LS-DYNA中，定义接触有很多方法. 不同之处包括接触面如何描述，接触罚函数如

35、何表达, 以及不同算法间存在的优缺点。 对一些接触模型， 接触面被用来定义接触体的表面. 这与隐式ANSYS 中的新的面对面算法是类似的, 只是用户不必定义这些接触面本身， ANSYS/LS-DYNA 会通过指定的节点组元自动生成接触面。另外一些接触模型允许模型的任何表面与其他表面接触， 包括它本身。 这种完全任意接触实际上是最容易定义的, 而且在预先不知接触表面（如整车碰撞模拟）时是非常有用的。培训ppt 接触表面. 接触表面概述一般来说，当一个接触（从）节点或接触面穿透目标（主）面时, 恢复力（罚力）会迫使之返回边界。接触罚刚度由LS-DYNA程序基于接触体的杨氏模量自动计算而得。 这就是

36、在定义刚体（材料）时，必须定义真实的杨氏模量的重要原因。ANSYS/LS-DYNA 可以模拟很大范围的接触状态: 表面抛光 通过定义带有剪切失效应力的速度相关的摩擦来实现。侵蚀接触 当表面单元失效时，允许接触表面延伸到内部单元。边缘接触 允许一个壳单元的边检测另一个壳的边不同于面面接触的一种特性。培训ppt 接触表面. 接触表面概述在ANSYS/LS-DYNA中定义接触时, 只需要简单地指出接触表面（不总是必需的）、接触类型以及与给定接触类型相关的任何特定参数。由于许多不同的接触类型可以使用，因此确定哪一种接触类型能最准确的描述所建模型常常是非常困难的。理解ANSYS/LS-DYNA所提供的不

37、同接触算法和接触类型对正确选择给定模型的接触面是非常重要的。接触算法是程序处理接触面的方法。有三种不同的算法：单面接触点对面的接触面对面的接触培训ppt接触表面 接触表面概述接触类型是指具有某些特定的相似属性的一系列接触类型 。这里有十种不同的接触类:一般接触自动接触刚性接触固连接触固连失效接触 侵蚀接触边接触拉延筋接触成形接触二维接触ANSYS/LS-DYNA 支持 24 种此类接触算法和接触类型的组合。培训ppt每种接触算法和接触类型都将详细介绍，以帮助我们选择能够精确体现所模拟物理现象的最合理的接触模型。接触表面 接触表面概述培训ppt接触表面B. 单面接触算法当一个体的外表面接触其自身

38、或其他体的外表面时，可用单面接触算法建立接触。单面接触是最常见的接触类型 这是由于ANSYS/LS-DYNA 程序自动搜索模型的外部表面以确定是否发生穿透。由于包括了所有的外部表面, 因此不需要定义接触和目标表面。单面接触对于处理接触区域不能提前预知的自接触或大变形问题是非常有效的。培训ppt接触表面. 单面接触算法ANSYS/LS-DYNA中的单面接触算法处理有限的接触表面时只会引起CPU时间很少量的增加。大多数冲击和碰撞的问题需要定义单面接触，因为有的接触表面是不能预知的。当接触单元的接触穿透超过单元厚度的40%，单面接触算法将被自动解除。 这对以下状态存在潜在的问题:过薄的部件具有低刚度

39、值的柔性材料非常高速运动的物体之间的接触培训ppt接触表面. 单面接触算法这些状况会导致接触节点超出40%的深度条件 。接触算法将假设表面不再接触，而且将会发生材料渗透，最终接触节点残留在目标表面后面。单面接触算法不在ASCII rcforc的文档中记录总的接触力 。如果想得到接触力，则应该采用点对面或面对面算法。单面接触算法包括以下接触模型: 单面 (SS)， 自动单面 (ASSC)，自动一般 (AG)， 侵蚀单面 (ESS)，单边 (SE)，和二维自动单面 (ASS2D) 。培训ppt当接触节点穿透目标表面时，采用点对面接触算法建立接触。因其不对称性，此法是最快捷的算法， 此算法只处理冲击

40、目标表面的接触节点。对于点对面接触算法, 必须定义接触和目标表面的节点的组元或part （part集）号，这类似于ANSYS 隐式方法。点对面接触算法对于接触区域相对较小而且接触区域已提前预知的情况非常有效。并对节点与刚体的接触也非常有效。接触表面C. 点对面接触算法培训ppt以下是用于点对面接触的指南 : TARGET SURFACE CONTACT NODES平面或凹面应作为目标面而凸面应作为接触面较粗网格应为目标面 而较细网格应为接触面3. 对拉延筋接触，拉延筋总是节点接触表面而板料则是目标面。接触表面 .点对面接触算法点对面接触算法在ASCII rcforc的文档中记录接触合力 。点对

41、面接触算法包括以下接触模型: NTS, ANTS, RNTR, TDNS, TNTS, ENTS, DRAWBEAD, FNTS培训ppt当一个体的表面穿透另一个体的表面时，采用面对面接触算法建立接触。面对面接触是完全对称的，因此接触面与目标面的选择是任意的。对面对面接触，需要定义接触面&目标面节点组元或part （或part集）号。节点可以属于多个接触表面。面对面接触是一种普遍的算法，常应用到具有大的接触区域且接触表面已知的情况。接触表面D. 面对面接触算法培训ppt面对面接触算法在ASCII rcforc的文档中记录接触合力 。面对面接触算法包括以下接触模型: STS, OSTS, AST

42、S, ROTR, TDSS, TSES, TSTS, ESTS, FSTS, FOSSCONTACT AND TARGET DEFINITIONS ARBITRARYV面对面接触对产生大量相对滑移的接触（如一个木块在平面上的滑移）非常有效。接触表面. 面对面接触算法培训pptShell TOP surfaceShell BOTTOM surfacepenetrating nodeGeneral ContactContactRestoring Force接触表面E. 一般接触类一般接触 只考虑壳体单元一侧接触。对实体单元（brick） ， ANSYS/LS-DYNA 总是用外法线创建目标面 ，因

43、此所期望的接触总能被识别。恢复力（即抵抗节点穿透的惩罚力） 将随着接触节点穿透目标表面而持续增加。节点经过一定壳单元厚度后，此力不会被消除。一般接触类型包括:SSNTSSTSOSTS培训ppt接触表面F. 自动接触类自动接触考虑壳体单元两侧的接触。壳体接触表面的方向是自动确定的。恢复力将随着接触节点穿透目标表面而持续增加，但只能增加到一点， 这是由于壳单元的两个面都需要检测。自动接触模型包括:ASSCAG (includes SE)ASS2DANTSASTSpenetrating nodeShell TOP surfaceShell BOTTOM surfaceAutomatic Contac

44、tContactRestoring Force培训ppt膝盖与挡板的接触接触表面G. 刚性接触类刚性接触 RNTR 和ROTR 类似于 NTS 和 OSTS 接触，但它不是用线性刚度来阻止穿透，而是定义力-变形曲线来实现。这种接触通常用于多刚体动力学， 它们允许在不必模拟变形单元 的情况下，吸收能量。刚体对可变形体的接触必须用一般、自动或侵蚀接触来定义。培训ppt接触表面H. 固连接触类固连接触用于连接两个不相似的网格，在很大程度上。这与ANSYS隐式中所用的绑定接触几乎是一样的。 接触节点粘接在目标表面上。两个面必须是初始共面的。目标面可以变形，迫使接触节点随之变形。固连接触只影响平移 DO

45、Fs (UX, UY, UZ) 固连接触模型包括： TDNS, TDSS, and TSES 接触。培训ppt接触表面I. 固连断开接触固连断开接触 实质是失效的 固连接触 。常被用来模拟点焊或螺栓连接 。 一旦符合失效方程, 接触节点就能够从目标面上滑移或分离出来。TNTS 失效方程基于法向力或剪切力, 而TSTS失效方程基于法向或切向应力。fstied失效 (联接解除)fnffffnnfailmssfailm,+121培训ppt接触表面J. 侵蚀接触类型这些接触模型 (ESS, ENTS, ESTS) 应用于外表面上的实体单元发生失效的情况 (例如, 由于超过允许应变值)，从而需要内部实体

46、单元承担抵抗穿透的作用。培训ppt接触表面K. 边接触类型边接触用于壳单元的面法线与碰撞方向正交的情况。壳边接触（SE）自动选定所有的边线。 SE 接触也包含在自动一般 (AG) 接触中。培训pptF = Ffriction+ FbendingDepth of drawbead接触平面L. 拉延筋接触类型拉延筋接触一般用于对坯料有特殊要求的金属成形过程 。例如，冲压过程中经常由于坯料不能贴紧模具而引起了坯料的褶皱。拉延筋接触通过包含确保坯料在整个拉伸过程中与模具始终接触的弯曲和摩擦限制力来模拟实际的拉伸过程。培训ppt接触表面M. 成形接触类型成形接触选项用于金属成形过程中点对面 (FNTS)

47、成形, 面对面 (FSTS)成形以及单向面对面 (FOSS)成形。对于这些接触类型， 冲头和模具一般定义为目标面 而工件则定义为接触面。这些接触类型网格无需连通, 因此减小了接触定义的复杂性。冲件网格的方向必须一致。成形接触选项基于自动接触类型，因此功能十分强大。培训ppt接触表面N. 2-D 接触类由 PLANE162 单元组成的模型, 只能定义为二维接触。由 ANSYS/LS-DYNA支持的2D 接触选项是ASS2D.默认的接触表面是整个模型 (即不采用目标面)。接触表面 可以限定为一个part集 (EDASMP).2D 接触只支持静摩擦系数 (FS)。培训ppt接触表面O. 定义基本接触

48、ANSYS/LS-DYNA 中定义接触的三个基本步骤是:选择接触模型 能最好地模拟你的物理系统识别接触实体 (单面接触不需要)定义接触 及任何必要的其它输入参数选择接触模型 :选则最适合于给定条件的接触模型 ， 参考 ANSYS/LS-DYNA Users Guide并回顾本章的内容。对大多数分析，推荐使用自动一般 (AG), 点对面 (NTS), 和面对面 (STS) 接触，它们都非常强大，并适合于大多数应用。培训ppt接触平面. 定义基本接触识别接触实体 :除了单面接触模型外 （SS, ASSC, AG, ESS, SE, and ASS2D），所有ANSYS/LS-DYNA 接触模型都需

49、要定义接触面和目标面 。虽然单面接触模型不必定义, 但通过定义接触面它能够限定部分模型的接触 。对于缺省值， 此算法考虑了整个模型的可能接触 。接触面和目标面能通过节点组元（CM command）或 part号/part 集 （EDPART and EDASMP commands）来定义。 如第1-3章讨论到的，如果在改变模型后重建part ，则part列表会改变，因此有必要在模型修改后经常更新part列表。否则，用以前的part号来定义接触将得到错误的分析结果。 最好的做法是在执行任何ANSYS/LS-DYNA的特定命令前完成所有的ANSYS几何模型的建立 。培训ppt接触表面. 定义基本接

50、触如果利用part号 (或part集) 识别接触面与目标面 ，可以用box ID 对接触面进一步限制。 Box ID 由 EDBOX command定义:Preprocessor LS-DYNA Options Contact Define Box这些方形体积不能用于 节点组元。限制接触处理的量在大型模型中更重要一些 ，但在包括整个模型可能导致接触行为不合适的小模型中也是很有用的。培训ppt接触表面. 定义基本接触定义接触 :EDCGEN 命令可在接触面和目标面之间自动生成接触 :Preprocessor LS-DYNA Options Contact Define Contact 先选择需要

51、的接触算法和相应的模型: 接触模型 接触算法培训ppt接触表面. 定义基本接触接下来规定静态和动态摩擦系数 以及指数衰减系数，粘性摩擦和粘性阻尼。ANSYS/LS-DYNA中，摩擦系数可以与速度相关并且可以使用一个限制值来定义总的摩擦力。缺省值定义模型为无摩擦接触 。培训ppt接触表面. 定义基本接触摩擦系数， mc, 的定义是:DC 或v = 0类似于 mc = ms 最大摩擦力 Flim可以由粘性摩擦系数VC和接触段的面积Asegment来定义: Flim = VC Asegment极限摩擦系数常用于接触引起塑性流动的情况。 推荐其值为 VC = 0.577 so, 这里 so 表示接触材

52、料的屈服应力. mc = md + (ms - md )e -(DC) (v) wherems = 静摩擦系数md = 动摩擦系数 DC = 指数衰减系数v = 接触面间的相对速度培训ppt接触表面. 定义基本接触为避免接触中产生不期望的振动，可利用粘性阻尼系数VDC施加垂直于接触面的接触阻尼 。输入值 VDC，需要一个阻尼值作为临界百分比:x = (VDC/100.0) xcrit = (VDC/100.0) (2m w)这里m = 质量w = 接触段的自然频率在金属成形分析中，粘性阻尼常用于抑制法向振动 。培训ppt接触表面. 定义基本接触定义摩擦数据后要定义接触起始和终止时间:利于多级成

53、形过程中设定后一阶段使用的工具到相应位置而不会与坯料互相影响。默认值在整个运算时间激活接触。如果运用已定义 的Part 号或 Part 集，给定 Box ID 用来进一步限制接触面或目标面的定义。由于当前 UIDL的局限性，只有在整个EDCGEN命令给出后， GUI 才能识别 Box IDs 是否是有效。 因此只有当Box IDs (用EDBOX命令) 、Part IDs 或 集合已经被定义才可以设定Box IDs。培训ppt接触表面. 定义基本接触在第一个 EDCGEN 对话框点击 OK后, 又弹出一个对话框以区分接触面和目标面。对于单面接触模型, 第二个框一般不出现，这是因为模型不需要区分

54、接触面和目标面 。对单面接触模型，EDCGEN 命令需直接键入 (即批处理)来获得进一步限制接触面的特性。有些接触模型需要输入附加输入参数。对这些模型，相应的将出现特定的模型对话框 培训ppt接触表面. 定义基本接触需要输入附加数据的接触系列是 (EDCGEN fields V1-V4):接触的模型指定输入进一步定义了接触。 在左边介绍的 Tiebreak点对面接触中 (TNTS) ，会出现对话框，要求规定两个面何时脱开。培训ppt接触表面P. 接触的列表和显示由于接触单元不能创建，因此在求解前列出定义了的接触表面对确保接触的正确定义是非常重要的。 这用 EDCLIST command来完成:

55、Preprocessor LS-DYNA Options Contact List Entities所有定义过的实体与指定的摩擦参数一起自动列出。 注意每种接触定义有专门的参考号培训ppt接触表面. 接触的列表和显示利用由EDCLIST输出的已定义的接触 (reference) 号，可以由EDPC 命令显示接触对 :Preprocessor LS-DYNA Options Contact Select and Plot.Contact #1 Contact #2 Contact #3 EDPC既选择又显示出接触对。在求解前不要忘记重选需要的节点和单元 培训ppt接触表面Q. 删除/撤消/恢复接

56、触用 EDDC 命令删除定义错误的接触实体:Preprocessor LS-DYNA Options Contact Delete Entity 接触实体只能在新分析中删除。对于EDDC 命令， GUI 自动提供 “DELE” 选项, 但仍需要指定接触和目标组元和/或Part ID/集，单面接触一般都不需要这些表面。 删除模块中所有接触, 键入命令 EDDC, DELE, ALLEDDC, DELE 需要:接触模型接触part或组元目标Part 或组元培训ppt接触表面.删除/撤消/恢复接触小型重启动分析中，可以撤消(EDDC, DACT)或激活(EDDC, RACT) 一个定义了的接触实体。

57、Preprocessor LS-DYNA Options Contact Deactivate Entity Preprocessor LS-DYNA Options Contact Activate Entity 接触的撤消和激活在多阶段成形模拟很有用 。EDDC 命令不支持完全重启动。培训ppt接触表面R. 高级接触控制ANSYS/LS-DYNA 有许多高级接触控制选项大多数控制可以由 EDCONTACT命令定义，而且这些控制是全局的 (即适合所有定义了的接触实体):Preprocessor LS-DYNA Options Contact Advanced Control可以用EDSP 或

58、者EDCMORE 命令来定义某些局部(如对特定的实体)的高级控制 :Preprocessor LS-DYNA Options Contact Advanced Control Preprocessor LS-DYNA Options Contact Additional Parms 以下 全局选项 一般用于控制接触:控制接触搜索方法控制接触深度控制接触刚度培训ppt接触表面. 高级接触控制控制接触的搜索方法 :ANSYS/LS-DYNA中,两个不同 接触搜索方法 :网格连贯性搜索块排序方法网格连贯性搜索中, 接触搜索算法利用邻近单元的共享节点来搜索接触。当目标段与接触节点失去接触时， 检查邻近

59、段。网格连贯性方法是极快的，但是由于它要求接触表面具有连续的网格，因此又有局限性。过去，网格连贯性方法用于接触模型STS, NTS, OSTS, TSTS, TNTS, and TDNS 的缺省，但是新的块排序算法逐渐成为某些模块的缺省算法。培训ppt接触表面. 高级接触控制在块排序方法中,由目标面所占据的三维空间被分为许多立方体（buckets）。节点可以与同一块或相连的块中的任一段相接触。 bucket sort 算法功能十分强大, 但是在某种程度上比网格连贯性搜索要慢, 尤其针对大模型。 网格连通方法是ASTS, ANTS, ESTS, ENTS, 及所有的单面接触模型的缺省方法。由于许

60、多模型含有不连续网格, 最好通过EDCONTACT 命令的壳体厚度选项(SHTK) 将搜索算法转变为块排序方法 将 EDCONTACT 的SHTK 域设置为 Thk Incl-Exc Rgd (利用除刚体外的壳厚度)或 Thk Incl (包括刚体在内的壳厚度)将导致如 NTS的接触模型使用bucket sort算法。培训ppt接触表面. 高级接触控制控制接触深度 :对于接触选项STS, NTS, TNTS, TSTS, and OSTS， ANSYS/LS-DYNA假定搜索的接触深度 为1e10 (如果不考虑厚度)。当接触节点穿过目标面时，就会产生一个与接触深度成比例的接触（恢复）力。在动力