《AUTODYNChapter9材料模型》由会员分享,可在线阅读,更多相关《AUTODYNChapter9材料模型(48页珍藏版)》请在金锄头文库上搜索。
1、AUTODYN-Chapter-9-AUTODYN-Chapter-9-材料材料模型模型Material Models9-2ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual显式动力学显式动力学 (ANSYS)中的材料模型中的材料模型Failure ModelStrength ModelEquation of StateAUTODYNMaterial Models9-3ANSYS, Inc. Proprietary 2009
2、ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualAUTODYN 材料模型材料模型 AUTODYN 材料模型包含三个部分状态方程 (EOS)强度模型失效模型 EOS Strength Failure和显式动力学和显式动力学(ANSYS)中的模中的模型一样型一样Material Models9-4ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665T
3、raining ManualAUTODYN附加材料模型附加材料模型Ideal Gas 状态方程Two Phase 状态方程SESAME 表格累积损伤模型Beam 强化模型破片分析刚体材料 正交各向异性材料正交各向异性实体复合壳高能炸药 (HE)爆炸 爆轰产物(气体)燃烧效应点火和增长慢燃药用户自定义材料Material Models9-5ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual带能量项的 状态方程EOSg g = 理
4、想气体参数r r = 密度, e = 比内能绝热常数, C输入一个不等于0的值计算绝热反应P/g g = C压力换算减去大气压力Ideal Gas 状态方程状态方程Material Models9-6ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual用来表达过热流体材料的膨胀和气化行为例如. 反应堆冷却和 可压缩状态方程一起使用 单相区域使用Gruneisen状态方程饱和曲线作为参考曲线材料饱和曲线可以再用户子程序 EXTAB
5、中定义 AUTODYN提供水的饱和曲线Two Phase状态方程状态方程比容流体区域两相区域流体和气体气体区域压力Material Models9-7ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualSesame 库库 Sesame库作为存储状态数据的表格,并不是状态方程包含200多种材料,有金属、 矿物材料, 聚合物和混合物应用于一些特殊材料,密度和内能均有很大的变化范围,以满足特殊的变化需要谨慎使用Sesame 库属于美国
6、政府控制出口的产品不包含在标准的发行包里如果允许的话可以从ANSYS 获得可以从 LANL直接获得Material Models9-8ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual考虑积累效应,强度逐步降低 早期模型主要用于脆性材料的先前的Johnson-Holmquist 模型使用用户子程序建立模型多个用户子程序的有效结合累积损伤失效模型累积损伤失效模型Material Models9-9ANSYS, Inc. Prop
7、rietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual模型的强度数据通过10点分段线性曲线定义轴向力- 轴向应变 轴11力矩-曲率 轴 11力矩-曲率 轴 22力矩-曲率 轴 33试验得到的载荷-位移曲线直接添加到模型里面,以便取得理想的结构响应定义单元轴向、扭矩和弯矩响应的四条曲线数据不存在内在联系Beam强化模型强化模型Material Models9-10ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights
8、 reserved.February 27, 2009Inventory #002665Training ManualCourtesy of AWE (A), UKBeam强化模型强化模型例子: 比例为1/3 的拉覆试验试验失效载荷: 86kN 4KN模拟失效载荷: 83kN 5KNMaterial Models9-11ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualCourtesy Sandia National Lab.
9、查看和列表显示破片的形成特性例子: 子弹出膛后爆炸破片分析破片分析Material Models9-12ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual用来定义刚性体显式动力学(ANSYS)中定义刚性体的不同方式在标准材料输入窗口选择 “EOS Rigid” 用刚性材料填充任意非结构化网格不适用于结构化网格刚性材料单元就是一个带有质量/惯量特性的单个刚体质量/惯量特性定义方式材料密度和单元体积材料中有明确定义使用一个刚体材料
10、定义多个刚体刚性材料刚性材料Material Models9-13ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual可变性弹丸 刚性弹丸 例子: 3D 斜侵彻刚性材料刚性材料Material Models9-14ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manua
11、l例子: 金属板成型刚性冲头和模具非结构壳 (主四边形) 工件冲头工件模具刚性材料刚性材料Material Models9-15ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualAUTODYN 可以为模型定义大载荷谱条件下的各向异性材料模型正交各向异性线弹性响应 (结构加载)正交各向异性弹性刚度矩阵线性体积响应正交各向异性线弹性响应 与非线性状态方程的耦合 (瞬态冲击加载)改进的正交各向异性弹性刚度矩阵非线性体积响应正交各向异
12、性塑性广义二次屈服面正交各向异性失效损伤失效脆性失效正交各向异性材料正交各向异性材料Material Models9-16ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual正交各向异性 状态方程正交各向异性软化正交各向异性屈服联合使用正交各向异性 状态方程, 屈服和软化模型得到完整的响应正交各向异性材料正交各向异性材料Material Models9-17ANSYS, Inc. Proprietary 2009 ANSYS,
13、Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual123等价的正交各向异性材料 - 只显示单个体压层复合材料使用固体单元描述正交各向异性材料或或正交各向异性材料正交各向异性材料Material Models9-18ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualS = C-1 =C =正交各向异性线弹性响应线性状态方程近似求出体
14、积响应正交各向异性材料正交各向异性材料Material Models9-19ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual正交各向异性线弹性响应 与非线性状态方程的耦合PolynomialShockPorous正交各向异性材料正交各向异性材料Material Models9-20ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27,
15、 2009Inventory #002665Training Manual屈服面形状通过系数矩阵定义, aij硬化参数, k广义形式简化为:Hills 正交各向异性屈服函数Von-mises 屈服函数正交各向异性塑性使用广义二次屈服面正交各向异性材料正交各向异性材料Material Models9-21ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual正交各向异性失效 : 脆性 Failure三个正交各向异性脆性失效准则材料应
16、力材料应变材料应力/应变可以在每个主材料方向指定拉应力/应变和剪切应力/应变正交各向异性材料正交各向异性材料Material Models9-22ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual正交各向异性失效 : 损伤模型模型的失效初始准则 (面) 正交各向异性材料正交各向异性材料Material Models9-23ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights r
17、eserved.February 27, 2009Inventory #002665Training Manual正交各向异性失效 : 损伤模型失效开始后, 程序会计算损伤张量用来降低屈服面强度正交各向异性材料正交各向异性材料Material Models9-24ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual KEVLAR-epoxy材料的静态拉伸实验正交各向异性材料正交各向异性材料Material Models9-25A
18、NSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual例子: 破片碰撞 GFRP 靶正交各向异性材料正交各向异性材料Material Models9-26ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual层状复合壳主要用于结构载荷(非冲击载荷)作用下的薄壳结构层状复
19、合壳通过 “Fill” 功能建模选择 Composite 按钮每一层连同法线初始条件被应用到网格层数可以使任意的,可以选择和存储每一层可以使各项同性和正交各向异性对正交各向异性材料必须指定 11 方向每一层有一个厚度每一层能独立查看正交各向异性材料正交各向异性材料Material Models9-27ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual层状复合壳材料模型与标准壳单元协调的材料模型可以用于复合壳单元 正交各向异性材
20、料需定义材料方向Tsai-Wu, Hoffman 和 Tsai-Hill 失效准则包含压缩和拉伸失效仅体积失效材料方向11 和22 位于壳平面33 为厚度方向材料轴选项I-J-K (推荐)默认 11 : 指标k增加的方向q 为11 方向绕单元中心的转角22 垂直于 11 在单元平面X-Y-Z正交各向异性材料正交各向异性材料Material Models9-28ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual例子: 鸟撞机翼
21、(机翼使用复合壳)正交各向异性材料正交各向异性材料Material Models9-29ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual爆炸过程时间起爆初始起爆点/平面压缩起爆不推荐不完善的物理过程用点火和增长(ignition and growth) 模型代替爆炸产物的膨胀 (气体)JWL 状态方程 (Jones, Wilkins, Lee)高能炸药高能炸药Material Models9-30ANSYS, Inc. Pr
22、oprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualDetonation FrontsInitiation NodeCellT1T2S1S2T1 = S1 / DT2 = S2 / D时间起爆爆炸发生在点或平面 (使用者定义)爆轰波波前以爆速 D传播 T1 开始燃烧 T2结束燃烧化学能从 T1 到 T2线性释放燃烧分数随时间从 0.0 增加到1.0 单元变量 alpha = -T1, TT1高能炸药高能炸药 爆炸过程爆炸过程Material Models9
23、-31ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual时间起爆直接起爆起爆时间通过起爆点到单元中心的直线距离来计算起爆时间通过起爆点到单元中心的直线距离来计算 间接起爆起爆路径绕过惰性材料的最短有效路径决定起爆路径绕过惰性材料的最短有效路径决定直接起爆不能用直接起爆高能炸药高能炸药 爆炸过程爆炸过程Material Models9-32ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All
24、rights reserved.February 27, 2009Inventory #002665Training Manual时间起爆间接多点起爆 起爆点2是准确计算shadow区域起爆时间的关键高能炸药高能炸药 爆炸过程爆炸过程Material Models9-33ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual直接起爆间接起爆,1个起爆点间接起爆,2个起爆点间接起爆,3个起爆点高能炸药高能炸药 起爆路径起爆路径Mat
25、erial Models9-34ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualJWL状态方程用来描述炸药爆炸产物 (气体)的高速膨胀JWL 状态方程是一个经验公式,所需的数据均来源于数值和物理实验提供了很多常用炸药的数据膨胀气体的压力通过下式计算: A, B, R1, R2, w 为经验常数 , r = 密度, r0 = 参考密度, h = r / r0, e = 比内能log p高能炸药高能炸药 爆炸产物的膨胀爆炸产物的
26、膨胀log vMaterial Models9-35ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualJWL 状态方程输入参数包括状态方程参数爆速单位体积化学能 AUTODYN标准材料库中包含大量的炸药数据压缩比和燃烧前体积模量定义压缩起爆不推荐,默认值为0自动转变为理想气体考虑精度,推荐使用高能炸药高能炸药 爆炸产物的膨胀爆炸产物的膨胀Material Models9-36ANSYS, Inc. Proprietary 20
27、09 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualwhereQ = 附加的比能,a = 能量释放常数,m = 能量释放指数,n = 压力指数JWL 状态方程 Miller Extension非理想炸药,比如包含有铝 (Al) 或高锰酸铵 (AP) ,铝 (Al) 或高锰酸铵 (AP) 颗粒燃烧后会释放出更多的能量 Miller extension 模型表达了这种能量释放高能炸药高能炸药 爆炸产物的膨胀爆炸产物的膨胀Material Models9-37ANSYS, Inc. P
28、roprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualJWL 状态方程 - 能量释放范围的扩大 温压弹可以比传统的炸药释放更多的能量典型的为铝基炸药起爆后会消耗空气中的氧气 JWL 状态方程中附加的能量选项这个随时间变化的能量释放能量在指定的时间间隔内积聚高能炸药高能炸药 爆炸产物的膨胀爆炸产物的膨胀Material Models9-38ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved
29、.February 27, 2009Inventory #002665Training ManualJWL 状态方程 - 能量释放范围的扩大 10kg TNT球形装药,在0.12 和 0.55 msec之间附加能量2.15MJ/kg得到的结果更长的压力持续时间和更大的压力脉冲高能炸药高能炸药 爆炸产物的膨胀爆炸产物的膨胀Material Models9-39ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualLee-Tarver
30、 Ignition & Growth Model用于高能炸药 (HE)初始状态研究假设点火在局部过热点, 并从这些点向外扩张由三个部分组成:惰性炸药状态方程 (Shock或 JWL )JWL -反应产物的状态方程反应率方程描述燃烧的点火、生长和完成高能炸药高能炸药Material Models9-40ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual 0.5 km/s 0.7 km/s 1.0 km/sLee-Tarver I
31、gnition & Growth Model 例子: 冲击起爆高能炸药高能炸药Material Models9-41ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualPowder Burning 模型模拟主要物理特征为爆燃 的材料的燃烧(纵火器、弹药)两相模型单元内气体和固体同时存在固体相: Linear/Compaction 状态方程气体相: JWL/Exponential燃速, c, 和气体压力, Pg相关燃烧率和气体压力
32、 , Pg 以及燃烧分数, F有关有关公式: A Atwood, EK Friis and JF Moxnes, 一种带能量粉末材料燃烧的数学模型, 34th International Annual Conference of ICT, June 24-27, 2003, Karlsruhe Federal Republic of Germany慢燃炸药慢燃炸药GasSolid ParticlesNumerical Cell of Volume VMaterial Models9-42ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights re
33、served.February 27, 2009Inventory #002665Training ManualPowder Burn模型例子: 枪膛内发射带弹带的弹丸慢燃炸药慢燃炸药Material Models9-43ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual AUTODYN自带的材料数据的集合通过 Load,进入材料库材料可以通过 名称, 状态方程,强度模型和失效模型分类 所有的材料均有 EOS , 大多数有强度
34、模型,部分有失效模型 可以增加和修改材料库中的数据,并且可以自己建立材料库模型选择以后,数据会自动使用当前的单位制材料库材料库Material Models9-44ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual怎么使用材料模型怎么使用材料模型?我们如何对一个具体的材料进行模型选择呢?我们如何对一个具体的材料进行模型选择呢?根据材料本身,可以相对容易地大致确定材料种类根据材料本身,可以相对容易地大致确定材料种类液体或固体?液
35、体或固体?各项同向或各向异性?各项同向或各向异性?惰性惰性/活性?活性?多孔或没有多孔或没有?易延展或易碎易延展或易碎?压强是否依赖强度或粘性压强是否依赖强度或粘性?实际使用的材料很大程度依赖于使用的领域和可以得到的材料数据实际使用的材料很大程度依赖于使用的领域和可以得到的材料数据原则:材料模型近可能简单原则:材料模型近可能简单应该知道参数怎样影响结果,要得到好的结果那些是关键的参数Material Models9-45ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory
36、#002665Training Manual模块化的用户材料子程序:建立输入 GUI检查输入参数的连续性映射输入参数到求解参数写求解方程用 Fortran 90编写MDEOS_USER_1状态方程MDSTR_USER_1强度模型MDFAI_USER_1失效准则MDERO_USER_1侵蚀准则材料用户子程序材料用户子程序Material Models9-46ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training Manual子程序例子 : 强度模
37、型Module STR_USER_1声明标量和数组INIT_STR_USER_1定义输入参数并创建输入菜单SET_STR_USER_1复制输入参数到求解器的标量/数组变量CHECK_STR_USER_1检查输入参数是否有效SOLVE_STR_USER_1强度模型求解Material Models9-47ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.February 27, 2009Inventory #002665Training ManualGlobal and Material Erosion侵蚀是一个数值过程,在模拟过程中自动删除单元. 在退化前删除极度扭曲的单元 (degenerate). 保证时间步合理和足够大.确保求解能够完成. 能够用于材料破碎、切削和穿透在Explicit Dynamics (ANSYS), 能够为所有模型指定侵蚀准则在AUTODYN中, 可以为每一种材料定义侵蚀准则侵蚀不是一个物理结果 (或材料特性). 是为了防止网格畸变AUTODYN中有五种定义侵蚀的选项Geometric StrainPlastic StrainTimestepFailureUser Erosion用户子程序 EXEROD结束结束