材料模型与状态方程
不同外罩材料双层药型罩射流成型研究
231中国设备工程 2020.03 （上）为了提高小口径破甲弹的威力，美国提出了分离式装药的概念[1]；鲁修国[2]等研究了双层罩外罩材料及罩间距在F 装药结构中对射流特性的影响。
沈慧铭[3]等发现变壁厚双层药型罩形成的侵彻体可以提高战斗部的侵彻能力。
滕桃居[4]等将数据仿真和试验仿真两种方法进行结合应用研究，对不同的材料双层药型罩EFP 的成型和侵彻性能进行系统的研究。
王哲[5]等以微元爆轰驱动及碰撞理论为基础，搭建了双层药型罩串联EFP 速度分析模型。
针对现有环形装药结构存在的两点不足，一是对舱适应性较差，舱中的设备会使射流造成大量消耗，影响终点威力；二是舱的存在相当于增加了炸高，在作用过程中射流被充分拉伸，极易断裂，造成侵彻能力不稳定。
考虑建立双层药型罩结构，研究不同药型罩外罩材料对射流成形性能的影响。
1 数值模拟方案及有限元模型1.1 模型建立为了开展不同外罩材料双层药性罩射流成型研究，文中搭建了相应的计算模型，首先，假设空气、炸药、药型罩和装药壳体均为连续介质及整个爆炸过程为绝热过程。
并根据需要，将装药高度设计为120mm，其最大外径为520mm。
1.2 料模型与状态方程炸药选用B 炸药，RDX 和TNT 混合比为3比2，文中采用MHEB 模型和E-J状态方程开展相应的计算。
式中：p 为爆轰压力，E 为炸药比内能，υ为相对比容。
内层药型罩选用紫铜，外层药型罩选用2024T-351硬铝、钨以及某型硬质钢，均采用M-S 材料模型和E-G 状态方程。
在选材上，壳体采用高强度钢30CrMnSiNi2A，材料呈应变硬化规律；钢板选用舰船用某型钢，环形射流切割过程中钢靶材料应变率高，在这样高应变率情况下，材料的力学性能与准静态情况有较大的差异。
对于多介质ALE 算法而言，还需建立覆盖整个射流形成、炸药爆轰范围的空气网格，并在边界节点上施加压力流出边界条件，避免压力在边界上的反射。
2 双层药型罩射流形成过程内层紫铜药型罩经历了类似单层紫铜罩的变化过程：药不同外罩材料双层药型罩射流成型研究陆志毅1，苏卿1，李大超1，刘宏杰2（1.海装驻上海地区第十军事代表室，上海 200241；2.海军潜艇学院，山东 青岛 266044）摘要：环形聚能装药结构中，单层罩射流容易偏斜和飞散，不利于射流拉伸成型，最终影响射流的侵彻性能，双层罩射流中外罩材料形成的杵状物对内层射流起包覆作用，使得内层射流不易发生偏斜。
微声聚能穿孔装置优化设计及IED安全销毁实验
微声聚能穿孔装置优化设计及IED安全销毁实验徐全军;白帆;李裕春;龚自正;张庆明【摘要】Based on the Held initiation theory, and the fact that shaped charge jet against the charge can produce low velocity explosion, this paper presents a new design of micro-acoustic shaped charge device. The orthogonal analysis methods were applied to optimize the charge structure. The experiments indicate that the noise elimination ability and the anti-detonation ability of the device are fine. The simulator charge boxes were destroyed with low velocity detonation, which helps reduce the explosive harm and allows for destroying IED safely.%根据Held引爆理论,利用聚能射流使炸药只燃烧或发生不完全爆轰的特点,设计了微声聚能穿孔装置,并对装药结构进行了正交优化.实验结果表明,装置的消声能力良好,用于销毁模拟爆炸药盒时,装药仅发生低速爆轰,达到了安全销毁简易爆炸装置的目的.【期刊名称】《爆炸与冲击》【年(卷),期】2012(032)001【总页数】5页(P79-83)【关键词】爆炸力学;安全销毁;微声聚能穿孔装置;简易爆炸装置;正交分析【作者】徐全军;白帆;李裕春;龚自正;张庆明【作者单位】解放军理工大学工程兵工程学院,江苏南京210007;解放军理工大学工程兵工程学院,江苏南京210007;解放军理工大学工程兵工程学院,江苏南京210007;北京卫星环境工程所,北京100094;北京理工大学,北京100081【正文语种】中文【中图分类】O383简易爆炸装置（improvised explosive device，IED）是指利用就便器材临时制作的爆炸装置。
空气中TNT爆炸冲击波超压峰值的预测及数值模拟_杨鑫
ForecastandSimulationofPeakOverpressureofTNT ExplosionShockWaveintheAir
YANGXin, SHIShao-qing, CHENGPeng-fei (DeptofCivilEngineering, LogisticalEngineeringUniversity, Chongqing400041, China)
97
-0.0Z353
72
+0.0
00 Z4
6
25,
0.05
≤
Z≤
0 .3
Pso
=
0.619 38 Z
-0.0Z322 62
+0.2Z133 24,
0.3 ≤ Z≤ 1
0 .066 Z
2
+0.Z4205
+0.3Z238
8,
1 ≤ 绍了 TNT爆炸冲击波峰值超压 (单位 :MPa)的一种表达式[ 1] :
Abstract: Differencesbetweenseveralscholars′forecastsofpeakoverpressurehavebeencomparedbythecon-
trastiveanalysisofresearchonTNTexplosionshockwavespreadingintheair.Anewformulahasbeensummedup, andLS-DYNA hasbeenusedtosimulatethespreadingofshockwaveintheairgeneratedbyTNTexplosionondifferentconditions.Thenumericalvaluesarelessthanempiricalformula′s.Theresultsoftheexplosionwithheavier loadaremoreapproximatethantheforecastingresultsofempiricalformula.Soit′snecessarytodomoreresearchon explosionshockwaveparameters.
材料模型
第七章材料模型ANSYS/LS-DYNA包括40多种材料模型，它们可以表示广泛的材料特性，可用材料如下所示。
本章后面将详细叙述材料模型和使用步骤。
对于每种材料模型的详细信息，请参看Appendix B,Material Model Examples或《LS/DYNA Theoretical Manual》的第十六章（括号内将列出与每种模型相对应的LS-DYNA材料号）。
线弹性模型·各向同性（#1）·正交各向异性（#2）·各向异性（#2）·弹性流体（#1）非线弹性模型·Blatz-ko Rubber(#7)·Mooney-Rivlin Rubber(#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来定义材料模型。
第八讲 材料模型
型一样
Strength
Failure
AUTODYN附加材料模型
Ideal Gas 状 态方程
正交各向异性材料
正交各向异性实体
Two Phase 状 态方程 SESAME 表格
累积损伤模型
复合壳
高能炸药 (HE)
Beam 强化模型
破片分析
爆炸 爆轰产物(气体) 燃烧效应 点火和增长
刚性材料
用来定义刚性体 在标准材料输入窗口选择 “EOS Rigid” 用刚性材料填充任意非结构化网 格
不适用于结构化网格
刚性材料单元就是一个带有质量/ 惯量特性的单个刚体 质量/惯量特性定义方式
材料密度和单元体积 材料中有明确定义
使用一个刚体材料定义多个刚体
刚性材料
• 例子: 3D 斜侵彻
正交各向异性材料
层状复合壳
材料模型
与标准壳单元协调的材料模型可以用于复合壳单元 正交各向异性材料 需定义材料方向 Tsai-Wu, Hoffman 和 Tsai-Hill 失效准则 包含压缩和拉伸失效 仅体积失效 11 和22 位于壳平面 33 为厚度方向 材料轴选项 I-J-K (推荐) 默认 11 : 指标k增加的方向 q 为11 方向绕单元中心的转角 22 垂直于 11 在单元平面 X-Y-Z
JWL 状态方程
输入参数包括
状态方程参数 爆速 单位体积化学能
AUTODYN标准材料库 中包含大量的炸药数据 压缩比和燃烧前体积模量 定义压缩起爆
jc本构与点火增长状态方程
jc本构与点火增长状态方程JC本构与点火增长状态方程JC本构是指应力与应变之间的关系，是材料力学性质的基本描述。
在材料科学与工程中，JC本构模型被广泛应用于研究材料的力学行为以及破坏机理。
点火增长状态方程则是描述材料点火和燃烧过程的数学模型。
我们来了解一下JC本构模型。
JC本构模型是由美国学者J.C. Jaeger在20世纪40年代提出的。
该模型假设材料在弹性区域内的应力与应变之间存在线性关系，而在塑性区域内则存在一定的非线性关系。
具体来说，在弹性区域内，材料的应力与应变成正比，比例系数为弹性模量；而在塑性区域内，材料的应力与应变之间则存在一定的非线性关系，这种非线性关系可以通过材料的屈服强度来描述。
JC本构模型的核心是屈服准则，即确定材料何时进入塑性区域。
常见的屈服准则有屈服强度准则、能量准则和变形准则等。
其中，屈服强度准则是最为常用的，它假设材料的屈服发生在达到某个临界值的应力状态下。
根据相应的屈服准则，可以得到JC本构模型的应力-应变关系。
接下来，我们来了解一下点火增长状态方程。
点火增长状态方程是用于描述材料点火和燃烧过程的数学模型。
点火是指燃料与氧化剂之间的反应开始，并产生明显的火焰和燃烧现象。
点火过程包括引燃、燃料氧化和燃烧产物的生成等多个阶段。
点火增长状态方程的核心是热力学平衡条件和质量守恒方程。
热力学平衡条件可以用来描述燃烧过程中各个组分的浓度和温度之间的关系，而质量守恒方程则可以用来描述燃料和氧化剂之间的反应速率。
通过求解这些方程，可以得到点火过程中各个组分的浓度和温度随时间的变化规律。
JC本构与点火增长状态方程在材料科学与工程领域具有重要的应用价值。
在材料设计和工程实践中，通过建立适合材料特性的JC本构模型，可以预测材料在不同应力状态下的力学行为，从而指导材料的设计和使用。
而点火增长状态方程则可以用于研究材料的燃烧行为，对于燃料的选择和燃烧过程的控制具有重要意义。
JC本构与点火增长状态方程是材料科学与工程中两个重要的数学模型。
材料模型与状态方程
1 John-Cook 材料本构模型)1)(ln 1)((**mnp y T C B A -++=εεσ 式中，pε—— 等效塑性应变；*ε——0.10=εs -1的无量纲塑性比，0*εεε p =；*T——相对温度，roommelt roomT T T T T--=*A —— 屈服应力，Pa ；B —— 应变硬化系数，Pa ； n —— 应变硬化指数；C —— 应变率相关系数；m——温度相关系数。
表达式的第一项)(n B A ε+表示对于0.1*=ε和0*=T （等温状态）时的应力与应变的函数关系；表达式的第二项)ln 1(*εC +和第三项)1(*mT -分别表示应变和率温度的影响。
表 Johnson 和Cook 给出的值韩永要《弹道学报》第16卷第2期（断裂破坏时的）应变]1][ln 1][[*5*421*3T D D e D D D f +++=εεσ其中，D 1、D 2、D 3、D 4、D 5输入参数，σ*是压力与有效应力之比，eff p σσ/*=。
当破坏参数∑∆=fpD εε达到1时，发生破坏。
* Hirofumi Iyama, Kousei Takahashi, Takeshi Hinata, Shigeru Itoh ．Numerical Simulation of Aluminum Alloy Forming Using Underwater Shock Wave ．8thInternational LS-DYNA Users Conference2 Steinberg-Guinan 材料本构模型定义材料熔化前的剪切模量i m iE E fE c i e R E E h bpV G G --⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-'--+=300313/10p ——压力，V ——相对体积，Ec ——冷压缩能，Em ——熔化能AR R ρ='，R ——气体常数，A ——原子量 屈服强度ii m iE E fE c i y eR E E h pV b --⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-'--'+'=300313/10σσ如果Em 超过Ei ，[]np i y )(10εγβσσ++'= i γ——初始塑性应变，当0σ'超过m σ，设置0σ'等于m σ。
jc本构与点火增长状态方程
jc本构与点火增长状态方程JC本构与点火增长状态方程引言：JC本构模型是一种常用的材料模型，用于描述材料的本构行为。
而点火增长状态方程则是用来研究点火过程中的火焰传播速度与温度、压力等因素之间的关系。
本文将介绍JC本构模型和点火增长状态方程的原理和应用。
一、JC本构模型的原理JC本构模型是由James C. Simo和Thomas J.R. Hughes于1986年提出的一种材料模型，用于描述材料的弹性和塑性行为。
该模型基于弹塑性理论，将材料的应力-应变关系分为弹性和塑性两个部分。
在弹性阶段，材料的应变是可恢复的，而在塑性阶段，材料的应变是不可恢复的。
JC本构模型通过引入一个等效塑性应变来描述材料的塑性变形，并根据材料的本构参数来确定其应力-应变关系。
二、JC本构模型的应用JC本构模型广泛应用于材料的力学性能分析和仿真中。
例如，在有限元分析中，可以使用JC本构模型来模拟材料的变形行为，从而预测材料的应力分布和变形情况。
JC本构模型还可以用于模拟材料的破裂行为。
通过引入断裂准则，可以预测材料在受力下发生断裂的位置和形态，从而提高材料设计的可靠性和安全性。
三、点火增长状态方程的原理点火增长状态方程是用来描述点火过程中火焰传播速度与温度、压力等因素之间的关系的数学模型。
该方程基于化学动力学和热力学原理，通过考虑燃烧反应速率、热传导和质量传输等因素，来描述火焰的形成和传播过程。
点火增长状态方程通常由两个主要部分组成：点火延迟模型和火焰传播速度模型。
点火延迟模型用于描述点火的延迟时间，即燃料与氧气混合后到点火发生的时间。
火焰传播速度模型用于描述火焰在点火后的传播速度，即火焰的扩展速度。
四、点火增长状态方程的应用点火增长状态方程广泛应用于燃烧模拟和燃烧器设计中。
例如，在内燃机燃烧过程的模拟中，可以使用点火增长状态方程来预测燃烧的起始时间和速度，从而优化燃烧过程，提高燃烧效率和减少排放。
点火增长状态方程还可以用于火灾模拟和安全评估中。
空气中TNT爆炸的数值模拟
第31 卷第4 期2014 年12 月爆破BLASTINGV o l．31 N o．4D ec．2014d o i:10．3963 /j．i ss n．1001 －487X．2014．04．009空气中TNT 爆炸的数值模拟*胡兆颖，唐德高( 解放军理工大学国防工程学院，南京210007)摘要: 为了研究T N T 炸药爆炸产生的冲击波在空气中的传播规律和预测不同比例距离的超压峰值，应用LS-DYNA有限元软件模拟了7．5 k g T N T 爆炸的冲击波传播过程，揭示其能量衰减规律。
并用2 种经验公式计算不同比例距离的冲击波超压峰值。
对数值模拟结果、经验公式结果和已有的实验数据进行对比。
结果表明: 数值模拟结果与实验数据吻合较好，误差在10% 以内，证明了计算模型和参数的合理性。
2 种经验公式，叶晓华推荐的经验公式与实验数据的误差相对较小，距爆心3．5 m处，误差仅为0．66% 。
说明叶晓华公式相比H e n r yc h 公式更为可靠。
但随着爆距的增大，误差也明显增大。
建议此公式在比例距离小于2．6 m/k g1 /3 时采用。
关键词: 爆炸冲击波; 超压峰值; 有限元; LS-DYNA; 比例距离; 经验公式中图分类号: O625 文献标识码: A 文章编号: 1001 －487X( 2014) 04 －0041 －05Numerical Simulation of TNT Explosion in AirHU Zhao-ying，TANG De-gao( E n g i n ee r i n g I n s tit u t e o f N a ti o n a l D e f e n s e，PL A U n i v o f S c i＆ T ec h N a n ji n g 210007，C h i n a)A b s t r a c t: T o s t ud y t h e o f s h oc k w ave p r o p aga ti o n b e h av i o r o f T N T ex p l o s i o n i n t h e a i r a nd t o p r e d i c t t h e ove r-p r e ss u r e p ea k i n d iff e r e n t s ca l e d i s t a n ce，t h e fi n it e e l e m e n t s o ft w a r e LS-DYNA w a s u s e d t o s i m u l a t e t h e p r o p aga ti o no f s h oc k w ave b y7．5 k g T N T ex p l o s i o n a nd t o r evea l t h e a tt e nu a ti o n l a w s o f e n e r gy s h oc k w ave．E m p i r i ca l f o r m u l aw a s a l s o u s e d t o ca l c u l a t e t h e p ea k p r e ss u r e i n d iff e r e n t s ca l e d i s t a n ce．B y co m p a r i n g t h e r e s u lt s o f nu m e r i ca l s i m u l a-ti o n a nd e m p i r i ca l f o r m u l a w it h ex p e r i m e n t a l d a t a，t h e r e s u lt s s h o w t h a t t h e s i m u l a ti o n r e s u lt s ag r ee d w e ll w it h ex p e r-i m e n t a l d a t a a nd t h e e rr o r w a s l e ss t h a n 10% ．C o m p a r e d w it h t h e ex p e r i m e n t a l d a t a，t h e f o r m u l a pu t f o r w a r d b y Y EX i ao-hu a m a d e a s m a ll e r e rr o r，w h i c h w a s o n l y0．66% 3．5 m e t e rs a w ay f r o m ex p l o s i o n ce n t e r．Ｒe s u lt s h o w s Y E X i-ao-hu a f o r m u l a m o r e r e li a b l e t h a n H e n r yc h f o r m u l a．B u t a s t h e d i s t a n ce i n c r ea s i n g，t h e e rr o r i n c r ea s e do b v i o u s l y．Th e r e f o r e，t h e f o r m u l a co u l d b e a pp li e d w h e n t h e s ca l e d i s t a n ce w a s l e ss t h a n 2．6 m/k g1 /3 ．K ey wo r d s: ex p l o s i o n s h oc k w ave;p ea k ove r p r e s s u r e;fi n it e e l e m e n t;LS-DYNA;s ca l e d i s t a n ce;e m p i r i ca lf o r m u l a近年来，由于全球范围内极端势力和分裂势力的盛行，恐怖袭击层出不穷，给人员安全和建筑物造成了重大的损害。
某炸药库内爆作用下坑道内隔爆防护墙厚度估算
某炸药库内爆作用下坑道内隔爆防护墙厚度估算徐全军;渠银录;姜楠;夏裕帅;万仁红【摘要】文章以某炸药库安全评估工程的实际问题为背景,考虑到周围墙体的限制,坑道内的冲击波得以加强,以及附近坑道施工不受到影响为原则,通过数值模拟和简单公式,分析了10 t炸药库的爆炸破坏作用,提出了内爆作用下隔爆防护墙厚度的估算方法,并研究在一定的冲击波作用下隔爆防护墙的厚度,为安全评估提供参考.%In this paper, an actual project in explosive magazine safety evaluation is taken as the background. When the explosion happens in the magazine, consider the limit of the wall, the blast shock wave in the tunnel will strengthen, and it should not affect the adjacent tunnel project. Based on those above principles, the numerical simulation and the formula calculation are used to analyze the blast destruction effect, and a thickness calculation method of the explosion protection wall under the blast shock wave load is proposed. The thickness of the explosion protection wall under a certain blast shock wave is also studied, to provide references for safety evaluation.【期刊名称】《爆破器材》【年(卷),期】2012(041)002【总页数】4页(P34-36,40)【关键词】防护墙;冲击波;炸药库;坑道【作者】徐全军;渠银录;姜楠;夏裕帅;万仁红【作者单位】解放军理工大学工程兵工程学院江苏南京 210007;解放军理工大学工程兵工程学院江苏南京 210007;解放军理工大学工程兵工程学院江苏南京210007;深圳市蛇口招商港湾建设工程有限公司广东深圳 518067;解放军理工大学工程兵工程学院江苏南京 210007【正文语种】中文【中图分类】TU746.5;TD235.4引言随着我国经济的不断发展，炸药仓库周边的居民和工程施工越来越多，对炸药仓库安全设施的要求也越来越高。
ANSYS材料模型
第七章材料模型ANSYS/LS-DYNA包括40多种材料模型，它们可以表示广泛的材料特性，可用材料如下所示。
本章后面将详细叙述材料模型和使用步骤。
对于每种材料模型的详细信息，请参看App endix B,Materi al ModelExampl es或《LS/DYNA Theore tical Manual》的第十六章（括号内将列出与每种模型相对应的L S-DYNA材料号）。
线弹性模型·各向同性（#1）·正交各向异性（#2）·各向异性（#2）·弹性流体（#1）非线弹性模型·Blatz-ko Rubber(#7)·Mooney-Rivlin Rubber(#27)·粘弹性（#6）非线性无弹性模型·双线性各向同性（#3）·与温度有关的双线性各向同性（#4）·横向各向异性弹塑性（#37）·横向各向异性FLD（#39）·随动双线性（#3）·随动塑性（#3）·3参数Bar lat（#36）·Barlat各向异性塑性（#33）·与应变率相关的幂函数塑性（#64）·应变率相关塑性（#19）·复合材料破坏（#22）·混凝土破坏（#72）·分段线性塑性（#24）·幂函数塑性（#18）压力相关塑性模型·弹-塑性流体动力学（#10）·地质帽盖材料模型(#25)泡沫模型·闭合多孔泡沫(#53)·粘性泡沫(#62)·低密度泡沫(#57)·可压缩泡沫(#63)·Honeyc omb(#26)需要状态方程的模型·Bamman塑性（#51）·Johnso n-Cook塑性（#15）·空材料（#9）·Zerill i-Armstr ong(#65)·Steinb erg(#11)离散单元模型·线弹性弹簧·普通非线性弹簧·非线性弹性弹簧·弹塑性弹簧·非弹性拉伸或仅压缩弹簧·麦克斯韦粘性弹簧·线粘性阻尼器·非线粘性阻尼器·索（缆）（#71）刚性体模型·刚体（#20）7.1定义显示动态材料模型用户可以采用ANSYS命令 MP，MPTEMP，MPDATA，TB， TBTEMP和 TBDATA以及ANSY S/LS-DYNA命令 EDMP来定义材料模型。
3.3 材料模型种类及参数
3.3 材料模型种类及参数
用命令定义材料模型
线性 非线性 结合状态方程定义 离散单位模型 刚性材料 弹性材料 非线无弹性材料 泡沫材料 状态方程 刚性材料
3.3 材料模型种类及参数
1)线弹性模型 1)线弹性模型 线弹性模型: σ=Eε (σ<σb；或σ<σys) ε σ σ σ
2)非线性弹性模型 2)非线性弹性模型: 非线性弹性模型: σ=kεn =kε (σ<σb；或σ<σys)
3.3 材料模型种类及参数
3)理想刚塑模型： 3)理想刚塑模型： 理想刚塑模型 忽略弹性变形，也不考虑应变硬化。 忽略弹性变形，也不考虑应变硬化。 当σ<σys时， ε=0 σ (σ≤ (σ≤σys) 当ε>0时， 时 σ=σys σ 4)理想弹塑性模型： 4)理想弹塑性模型： 理想弹塑性模型 线弹性+理想塑性。 线弹性 理想塑性。 理想塑性 ε≤εys 当 ε≤ε 时, σ=Eε ε 当 ε>εys 时, σ=σys=Eεys ε σ ε (s≤sys) ≤
3.3 材料模型种类及参数
聚碳酸酯屈服应力与温度、应变速度之间关系 聚碳酸酯屈服应力与温度、
6061 T4 铝合金成形极限图 电液成形 铝合金成形极限图(电液成形 电液成形)
3.3 材料模型种类及参数 应力
真应力、 真应力、真应变
工程应力S 工程应变e: 工程应力S、工程应变e:
l0 ∆l σb l σys
σ=F/A=Fl /A0 l 0=(F/A0)[(l 0+∆ l )/ l 0]=S(1+e) >S ∆
隧道围岩爆破冲击损伤防护的轻气炮试验及数值模拟
ｐ ａ ｎ ａ ｇ ｆｒ ｃ ａ ｌ ’ｌｓｉ ｇ ｓ ｒａ ｅ；ａ ｄ ｗｉ ｈ ａ ｏ ｉ ｏ ｆ ｉ ａ ｔｖ ｌｃ ｔ ｔ ｅ ｔ ｃ ｅ ｈ ＶＣ— ｅ ｋ ａ ｄ ｄ ｍａ ｅ ｏ ｏ ｋ ｓ ｍｐ ｅｂ ａ ｔｎ ｕ ｆｃ ｎ ｔ ｔ ｅ ｓ ｍｅ ｃ ｎｄｔ ｎ ｏ ｍｐ ｃ ｅｏ ｉ ｈ ｉ ｙ，ｈ ｈｉｋ ｒｔ ｅ Ｐ Ｕ
ｍ ｔ ａ ｙｒ ｔｅｌ ｅ ｅｂａｔ ｇｓｒ ｃ ｒｓ ｅｋｖｌｅ ｂ ｓ ｇ ｍ ｐｏ ｃ ｖ ｙｒｔｅｉ ｕｙｔｌ ａｃ ｆ ｏｋ ａｅ ｌａｅ ， ｈ ｗ ｒ ｈ ｌｓｎ ｆ ｅｓ ｅｓｐａ ａ ； ｙｕｉ ｍ ｒ ｅｔ ｅｌ ｅ， ｊｒ ｅ ｎ ｅｏ ｃ ｉ ｒ ｌ ｏ ｔ ｉ ｕａ ｔ ｕ ｎ４ ｔ ｉ ａ ｈ ｎ ｏｒ ｒ
１ ．炮孔 ；．药包 ；．保 护层 材料 （ Ｖ — ２ ３ Ｐ Ｃ Ｕ） 图 １ 装药结构示 意图
Ｆ ｇ １ Ｓ ｅｃ ｐ ｏ ｄｆ ｄ ｃ ａ ｇ ｔ ｃ ｕ ｅ ｉ ． ｋ ｔ ｈ ｍａ ｆ ｍｏ ｉ ｅ ｈ ｒ ｅｓｒ ｔ ｒ ｉ ｕ
作 用 。该方 法在 实 际 中得 到 了应 用 ， 其 机 理 需 进 一 但
０ ８ ４ ２ 西南科技大学 环境与资源学 院， 四川 （． １ 同济大学 道路与交通工程教 育部重点试 验室 ， 上海 ２ １０ ；． 绵阳 ６ １１） ２０ ０
Ansys_AUTODYN计算求解一般步骤中文解析
-
C γ
eγ
ε
效
概
率
失效应变
• γ 由用户定义
• C 通过计算得到 • 分布类型：
– Fixed–每一时刻相同 – Random-任意 • 可用于许多材料的失效模型
一般问题的分析步骤－第二步 定义材料(续)
破片分析
选材总结－状态方程总结
常高 (大于105/秒)
• 其中， 为有效应变率；T为温度(K)；
• h v0 v 为压缩比 。
一般问题的分析步骤－第二步 定义材料(续) MO-Granular 强度模型
常用于干土、沙子、岩石、混凝土和陶瓷等材料。
压强硬化 10 点分段屈服应力-压强曲线
密度硬化 10 点分段屈服应力-密度曲线
避免出现压力低于10-6的单位制； 避免出现单元质量低于10-6的单位。
缺省单位：
长度 质量 时间 速度 力 压强 密度 能量
mm mg ms m/s mN kPa g/cm3 mJ
一般问题的分析步骤－第二步 定义材料
第一种方式：材料库选择材料模型 状态方程(体积改变)
Single phase、Multi-phase、Compaction、Explosive
P T1 T22 B0refe ( 0 ) 1
SHOCK EOS:copper Us c s1up s2u2p
“状态方程” 基本假定:容变律与畸变律解耦；忽略体积粘性。
一般问题的分析步骤－第二步 定义材料(续)
一些带孔的材料，变形过程中孔破裂导致不可逆 转的体积变形，
D
f
影响因素： f D1 D2eD3* 1 D4 ln * .1 D5T *
AUTODYN9材料模型
Material Models
刚性材料
• 用来定义刚性体
– 显式动力学(ANSYS)中定义刚性体的不同 方式
• 在标准材料输入窗口选择 “EOS Rigid” • 用刚性材料填充任意非结构化网格
– 不适用于结构化网格 • 刚性材料单元就是一个带有质量/惯量特性
的单个刚体
• 质量/惯量特性定义方式 – 材料密度和单元体积 – 材料中有明确定义
– 正交各向异性塑性
• 广义二次屈服面
– 正交各向异性失效
• 损伤失效 • 脆性失效
ANSYS, Inc. Proprietary © 2009 ANSYS, Inc. All rights reserved.
9-15
February 27, 2009 Inventory #002665
Material Models
© 2009 ANSYS, Inc. All rights reserved.
9-9
February 27, 2009 Inventory #002665
Material Models
Beam强化模型
• 例子: 比例为1/3 的拉覆试验
– 试验
• 失效载荷: 86kN 4KN
Training Manual
• 例子: 3D 斜侵彻
Training Manual
可变性弹丸
ANSYS, Inc. Proprietary © 2009 ANSYS, Inc. All rights reserved.
刚性弹丸
9-13
February 27, 2009 Inventory #002665
Material Models
• 使用一个刚体材料定义多个刚体
状态方程在材料科学中的应用
状态方程在材料科学中的应用材料科学是一门研究材料基本性质、性能和结构的学科，研究的范围包括金属、非金属、复合材料、纳米材料等各种材料。
在材料科学中，状态方程是一个非常重要的概念，因为它可以用来描述材料的性质和特性。
在本文中，我们将讨论状态方程在材料科学中的应用。
一、状态方程的定义状态方程是一组数学方程，它描述了材料在给定条件下的基本物理特性和状态。
这些条件可能包括：温度、压力、体积、熵等等。
状态方程可以用数学公式的形式表示。
其中，最常使用的状态方程之一就是理想气体状态方程。
二、理想气体状态方程理想气体状态方程是状态方程的一种常用形式。
它描述了理想气体在给定的温度、压力和体积下的状态。
其数学形式为：PV = nRT在这个公式中，P是气体的压力，V是气体的体积，n是气体的摩尔数，R是气体常数，T是气体的温度。
这个公式表明，当气体的温度升高时，它的体积扩大，压力也随之增加。
当压力升高时，气体的体积减小，温度也随之升高。
理想气体状态方程是理想气体模型的基石，是描述气体行为的基本工具。
三、状态方程在材料科学中的应用除了理想气体状态方程之外，还有很多其他的状态方程。
在材料科学中，这些状态方程被用来描述材料的性质和特性。
例如，状态方程可以用来描述材料的机械性质、热力学性质、电学性质等等。
以下是几个状态方程在材料科学中的应用案例。
1. 埃伦费斯特方程埃伦费斯特方程是一种用来描述材料的弹性力学性质的方程。
它包含了弹性恢复力、刚度、杨氏模量等参数，可以用来计算材料的弹性应变和弹性模量。
埃伦费斯特方程在研究材料的力学性质时有着广泛的应用。
2. 范德华方程范德华方程是用来描述材料间作用力的方程。
它是基于分子间作用模型而提出的，可以用来计算分子间的引力和斥力，从而描述材料的热力学性质。
范德华方程在研究材料的化学性质和相变规律时有着广泛的应用。
3. 内能方程内能方程是描述材料内部自身能量的方程。
它包含了材料的温度、压力、热容等参数，可以用来计算材料的内部能量变化和热力学性质。
AUTODYN_Chapter9_材料模型
? 正交正各交向各异性向失异效性: 材损伤料模型 – 失效开始后, 程序会计算损伤张量用来降低屈服面强度
? KEVLAR?-epoxy材料的静态拉伸实验
正交各向异性材料
? 例子: 破片碰撞 GFRP 靶
正交各向异性材料
? 层状复合壳
正交各向异性材料 – 主要用于结构载荷（非冲击载荷） 作用下的薄壳结构 – 层状复合壳通过 “Fill功”能建模
刚性材料
? 用来定义刚性体
– 显式动力学(ANSYS)中定义刚性体的不同 方式
? 在标准材料输入窗口选择 “EOS Rigid”
? 用刚性材料填充任意非结构 化网格
– 不适用于结构化网格
? 刚性材料单元就是一个带有 质量/惯量特性的单个刚体
? 质量/惯量特性定义方式
– 材料密度和单元体积
? 例子刚: 3性D 斜材侵料彻
f (?
ij )
?
a11?
2 11
?
a22?
2 22
?
a33?
2 33
?
2a12?
11?
22
?
2a23?
22?
33
?
2a13? 11?
33
?
2a44?
2 23
?
2a55?
2 31
?
2a66?
2 12
?
k
– 屈服面形状通过系数矩阵定义, aij – 硬化参数, k
– 广义形式简化为：
? Hills 正交各向异性屈服函数
– 多个用户子程序的有效结合
? 模 性型曲Be的线a强定m度 义强数化据通模过型10点分段线
– 轴向力- 轴向应变 – 力矩-曲率 – 力矩-曲率 – 力矩-曲率
LS-DYNA-5材料模式
材料模型 – 线弹性（续） • 各向异性: – 各向异性材料由 21 个独立参数和 DENS定义
1. 选择Linear Elastic > Anisotropic 2. 输入参数与密度 3. 输入坐标系 ID
ANSYS/LS-DYNA Training
材料模型 – 非线性弹性
• 在非线性弹性材料组中有3种模型 : – Blatz-Ko : 用于像橡胶一样的可压缩材料 – Mooney Rivlin : 用于定义不可压缩橡胶材料 – Viscoelastic : 定义玻璃类材料 非线性弹性材料可以经受大的可恢复的弹性变形
– 除了密度 DENS, 还需输入以下参数: • Go = The short term (origin) elastic shear modulus • G = The long term (infinity) elastic shear modulus • K = elastic bulk modulus • 1/ = decay constant ANSYS/LS-DYNA Training
•
Blatz-Ko Rubber:
• • •
Blatz-Ko 橡胶材料只用于受压橡胶材料 泊松比 (NUXY) ANSYS自动设为 .463, 因此只需要 DENS 和GXY 材料响应通过应变能量密度函数 (W)确定:
W Gxy I 2 I 2 I 3 5 2 3 w here 2 and 3 are the straininvariants .
材料模型 - 塑性类别 2
2b. 率敏感：
– 具有双线性等向硬化的塑性行 为 – 用强度系数K和硬化系数N定义 的幂法则 – 屈服函数如下定义:
1 P y 1 k e ε eff p C
一定深度水下爆破冲击波衰减规律的研究
一定深度水下爆破冲击波衰减规律的研究陈坤鑫;韩雪靖;杨嘉颖;谢浩【摘要】通过理论研究结合数值模拟的方式,采用ANSYS/LS-DYNA软件建立起水下爆破模型,得到水击波的传播图形,以及数值模拟波形产生锯齿状的结论和15m 水深处炸药爆炸的水击波峰值压力公式,可供相似条件下的水下爆破参考.【期刊名称】《矿业工程》【年(卷),期】2018(016)004【总页数】4页(P48-51)【关键词】水下爆破;水击波;衰减规律;数值模拟【作者】陈坤鑫;韩雪靖;杨嘉颖;谢浩【作者单位】福州大学紫金矿业学院,福建福州350116;福州大学紫金矿业学院,福建福州350116;福州大学紫金矿业学院,福建福州350116;福州大学紫金矿业学院,福建福州350116【正文语种】中文【中图分类】TD2350 引言我国航道和港口的修整和建造渐渐由浅海向深海发展，而使用炸药进行水下爆破的规模和起爆药量也逐步加大，爆破水击波对环境的威胁已经成为急需解决的问题。
水下爆破水击波是衡量水下爆破危害的重大因素，由于水下爆破情况较为复杂，水下地形、水流、水深、地质条件、施工方法以及起爆方式均会影响其传播。
因此，研究水击波的产生机理、衰减规律和防范措施显得尤为重要且十分迫切。
P.Cole[1]在大量试验研究的基础上，于1948年出版了《水下爆炸》的著作，通过一定的理论分析建立了一定范围内爆炸流场中冲击波压力峰值、比冲量及能量密度的计算公式。
前苏联学者B.B加尔基通过大量研究提出了水击波的一套计算公式[2]。
John F.Goarrnor等在前人的理论研究基础上研究了水面下、小深度处爆炸流场的解析解。
B.Mehaute和S.Wang汇集了关于水下爆破形成水中冲击波和其传播的理论基础及试验结果，分析了各种理论对应的限制和适用条件。
Sliflke J P研究了深水中爆破的压力特性，认为冲击波压力与深度等因素存在某种函数关系。
本文对水击波的传播机理进行系统阐述和深入研究。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1 John-Cook 材料本构模型)1)(ln 1)((**mnp y T C B A -++=εεσ 式中，pε—— 等效塑性应变；*ε——0.10=εs -1的无量纲塑性比，0*εεε p =；*T——相对温度，roommelt roomT T T T T--=*A —— 屈服应力，Pa ；B —— 应变硬化系数，Pa ； n —— 应变硬化指数；C —— 应变率相关系数；m——温度相关系数。
表达式的第一项)(n B A ε+表示对于0.1*=ε和0*=T （等温状态）时的应力与应变的函数关系；表达式的第二项)ln 1(*εC +和第三项)1(*mT -分别表示应变和率温度的影响。
表 Johnson 和Cook 给出的值韩永要《弹道学报》第16卷第2期（断裂破坏时的）应变]1][ln 1][[*5*421*3T D D e D D D f +++=εεσ其中，D 1、D 2、D 3、D 4、D 5输入参数，σ*是压力与有效应力之比，eff p σσ/*=。
当破坏参数∑∆=fpD εε达到1时，发生破坏。
* Hirofumi Iyama, Kousei Takahashi, Takeshi Hinata, Shigeru Itoh ．Numerical Simulation of Aluminum Alloy Forming Using Underwater Shock Wave ．8thInternational LS-DYNA Users Conference2 Steinberg-Guinan 材料本构模型定义材料熔化前的剪切模量i m iE E fE c i e R E E h bpV G G --⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-'--+=300313/10p ——压力，V ——相对体积，Ec ——冷压缩能，Em ——熔化能AR R ρ='，R ——气体常数，A ——原子量 屈服强度ii m iE E fE c i y eR E E h pV b --⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-'--'+'=300313/10σσ如果Em 超过Ei ，[]np i y )(10εγβσσ++'= i γ——初始塑性应变，当0σ'超过m σ，设置0σ'等于m σ。
材料熔化之后，y σ和G 设置为初始值的一半。
$ OFHC 为高导无氧铜，聚能装药药型罩常用材料 *MAT_STEINBERG$ MID R0 G0 SIGO BETA N GAMA SIGM2 8.93 0.477 0.120E-02 36.0 0.450 0.00 0.640E-02 $ B BP H F A TMO GAMO SA2.83 2.83 0.377E-03 0.100E-02 63.5 0.179E+04 2.02 1.50 $ PC SPALL RP FLAG MMN MMX ECO EC1 -9.00 3.00 0.00 0.00 0.00 0.00 0.00 0.00 $ EC2 EC3 EC4 EC5 EC6 EC7 EC8 EC9 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 *EOS_GRUNEISEN$ EOSID C S1 S2 S3 GAMAO A E0 2 0.394 1.49 0.00 0.00 2.02 0.470 0.00 $ V0 1.00M.Katayama, S.Kibe, T.Yamamoto ．Numerical and Experimental Study on the Shaped Charge for Space Debris Assessment ．Acta Astronauttca Vol.48,No.5-12,pp.363-372,2001W.H.Lee, J.W.Painter ．Material void-opening computation using particle method ．International Journal of Impact Engineering 22(1999)1-22二阶状态方程210221021)()(D C C B B B A A P +++++++=εεμεμμμμμ剪切模量G 与流体应力Y 间的本构关系⎥⎦⎤⎢⎣⎡--⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡-'-++=E x E fE R x E E h P b G G m )(exp 3003)(103/10η⎥⎦⎤⎢⎣⎡--⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡-'-+++=E x E gE R x E E h P qb e Y Y m n )(exp 3003)(1)1(03/10ηβmax 0)1(Y e Y n ≤+β，)('3)()(0x T R x E x E m m +=，α)1()2ex p()(0x ax T x T m m -=31)(20--=a γα，TAD R dx x P x E x '3)()(0-=⎰，ax ax TAD --=0)1()ex p(300γ3 Mie-Gruneisen 状态方程 定义压缩材料的压力为22332212020)1(1)1(12)21(1⎥⎦⎤⎢⎣⎡+-+---⎥⎦⎤⎢⎣⎡--+=μμμμμμμγμρS S S a C p定义膨胀材料的压力为E a C p )(020μγμρ++=其中：C 为u s -u p 曲线的截距，体积声速S 1、S 2、S 3是u s -u p 曲线斜率的系数，0γ是Gruneisen 常数，1/0-=ρρμ，a 是0γ的一阶体积修正。
(1)Hirofumi Iyama, Kousei Takahashi, Takeshi Hinata, Shigeru Itoh ．Numerical Simulation of Aluminum Alloy Forming Using Underwater Shock Wave ．8thInternational LS-DYNA Users Conference(2)M. Katayama, S. Kibe, T. Yamamoto ．Numerical and Experimental Study on the Shaped Charge for Space Debris Assessment ．Acta Astronauttca Vol.48,No.5-12,pp.363-372,2001(3)JingPing Lu, Helen Dorsett, David L. Kennedy ．Simulation of Aquarium Tests for PBXW-115(AUST) (4)S. Itoh, H. Hamashima ．Determination of JWL Parameters from Underwater Explosion Test (5)Katsuhiko Takahashi, Kenji Murata, Akio Torii, Yukio Kato ．Enhancement of Underwater Shock Wave by Metal Confinement4 多线性多项式状态方程 压力由下式定义E C C C C C C C P )(2654332210μμμμμ++++++=其中，1/0-=ρρμ，如果0<μ，则设置022=μC ，026=μC 。
当设置063210=====C C C C C ，154-==γC C 时，就可以用于符合γ律状态方程的气体，其中γ为比热系数。
E P 0)1(ρργ-=5空白材料在仿真计算中，水介质的材料模型可以选用空白材料（NULL），通过此材料来避免计算应力、应变。
在LS-DYNA中为材料模型9。
空白材料模型必须使用状态方程。
6炸药的材料模型在LS-DYNA中，炸药的材料模型一般都选用材料类型8，即MA T_HIGH_EXPLOSIVE_BURN。
需要定义的参数有，密度、爆速与C-J爆轰压力等。
此种材料类型必须与状态方程一块使用。
7 JWL 状态方程炸药爆轰产物的状态方程常采用JWL 方程。
此状态方程通常用于描述高能炸药及爆轰产物，其形式为VEe VR B e VR A p V R V R ϖϖϖ+-+-=--))1()1(2121(1) Hirofumi Iyama, Kousei Takahashi, Takeshi Hinata, Shigeru Itoh ．Numerical Simulation of Aluminum Alloy Forming Using Underwater Shock Wave ．8thInternational LS-DYNA Users Conference(2)M.Katayama, S.Kibe, T.Yamamoto ．Numerical and Experimental Study on the Shaped Charge for Space Debris Assessment ．Acta Astronauttca Vol.48,No.5-12,pp.363-372,2001(3)M.Katayama, S.Kibe ．Numerical Study of the Conical Shaped Charge for Space Debris Impact ．International Journal of Impact Engineering 26(2001) 357-368(4)Mark Z. Vulitsky, Zvi H. Karni ．Ship Structures Subject to High Explosive Detonation ．7thInternational LS-DYNA Users Conference(5)金乾坤等．3D Numerical Simulations of Penetration of Oil-Well Perforator into ConcreteTargets ．7thInternational LS-DYNA Users Conference(6)W.H.Lee, J.W.Painter ．Material void-opening computation using particle method ．International Journal of Impact Engineering 22(1999)1-22(7)JingPing Lu, Helen Dorsett, David L.Kennedy ．Simulation of Aquarium Tests for PBXW-115(AUST) (8)S. Itoh, H. Hamashima ．Determination of JWL Parameters from Underwater Explosion Testb.见LX-10-1炸药的各值，它们是类似的混合炸药James L. O’Daniel, Theodor Krauthammer, Kevin L. Koudela．An UNDEX response validation methodology．International Journal of Impact Engineering 27(2002) 919-937$铝*MAT_JOHNSON_COOK2 2.77 0.2763.37E-03 3.43E-03 0.410 0.100E-01 1.00 0.877E+03 300. 0.100E-050.875E-05 -9.00 3.00 0.00 0.00 0.00 0.00 0.000.00*EOS_GRUNEISEN2 0.5328 1.338 0.00 0.00 2.00 0.00 0.000.00$金属-铜*MAT_JOHNSON_COOK1 8.96 0.478 1.28 0.340000 0.0000000.900E-03 0.292E-02 0.310 0.250E-01 1.09 0.136E+04 294. 0.100E-050.383E-05 -9.00 3.00 0.00 0.00 0.00 0.00 0.000.00*EOS_GRUNEISEN1 0.394 1.49 0.00 0.00 2.02 0.470 0.000.00$炸药*MAT_HIGH_EXPLOSIVE_BURN2 1.7870001 0.8390000 0.3400000 0.0000000 0.0000000 0.0000000 0.0000000*EOS_JWL2 5.8140002 6.8010E-2 4.1000000 1.0000000 0.3500000 0.0900000 1.0000000$空气*MAT_NULL3 1.2250E-3 0.0 17.456E-6 0.0 0.0 0.0 0.0*EOS_LINEAR_POL YNOMIAL3 0.0 0.0 0.0 0.0 0.40 0.4 0.02.5E-6 1.0$水*MAT_NULL4 1.0 -1.00E-4 1E-6 0.0000000 0.0000000 0.0000000 0.0000000*EOS_GRUNEISEN4 .14840 1.7900000 .0000000 .0000000 0.110000 3.0000000 .00000001.0$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ MA TERIAL DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $*MAT_HIGH_EXPLOSIVE_BURN2 1.7870001 0.8390000 0.3400000 0.0000000 0.0000000 0.0000000 0.0000000*EOS_JWL2 5.8140002 6.8010E-2 4.1000000 1.0000000 0.3500000 0.0900000 1.0000000*INITIAL_DETONA TION2 0.0000000 0.0000000 0.0000000 0.0000000*MAT_NULL1 1.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000$水*MAT_0107 1.000E+00 0.2e-2 0.0 0.000 0.000E-00 2.00*EOS_GRUNEISEN$水$ EOSID C S1 S2 S3 GAMA0 A E06 0.1480 1.79 0 0.0000 1.65 0.0 2.895e-6$ V01.0*EOS_GRUNEISEN$水$ EOSID C S1 S2 S3 GAMA0 A E015 0.1647 1.921 -0.096 0.0000 0.35 0.0 2.895e-6$ V01.0*EOS_GRUNEISEN$水$ EOSID C S1 S2 S3 GAMA0 A E01 0.1480 2.56 -1.986 0.2268 0.50 0.0 2.895e-6$ V01.0$*MAT_JOHNSON_COOK$钨合金9 17.6 1.36 3.50 0.286 0.0000001.806E-02 0.177E-02 0.120 1.600E-02 1.00 1.723E+03 294. 0.100E-050.134E-05 -9.0 2.00 0.00 2.00 0.00 0.00 0.000.00$*EOS_GRUNEISEN$钨合金9 0.399 1.24 0.00 0.00 1.54 0.00 0.000.00$*MAT_PLASTIC_KINEMATIC$树脂8 1.19 0.780E-01 0.000000 0.800E-03 0.00 1.000.00 0.00 2.00$*MAT_PLASTIC_KINEMATIC$钨合金$ mid ro e pr sigy etan beta6 1.862E+01 1.170E+00 0.22 1.790E-02 1.0$ src srp fs0.8$$空气*MAT_NULL4 1.280E-03 0.000E+00 0.000E+00 0.000E+00 0.000E+00*EOS_LINEAR_POL YNOMIAL4 0.000E-00 1.000E-05 0.000E+00 0.000E+00 0.400 0.400 0.000E+00 0.000E+00 0.000E+00$$*MAT_JOHNSON_COOK$铁5 7.83000 0.777.920E-03 5.10E-03 0.260 0.140E-01 1.03 0.1793E+04 294 0.100E-05 0.477E-05 -9.00E+0 3.00 0.0 0.80 0.00 0.00 0.000.00*EOS_GRUNEISEN$铁5 0.4569 1.490 0.00 0.00 2.17 0.46 0.001.00$金属-铜*MAT_JOHNSON_COOK3 8.96 0.478 1.28 0.340000 0.0000000.900E-03 0.292E-02 0.310 0.250E-01 1.09 0.136E+04 294. 0.100E-05 0.383E-05 -9.00 3.00 0.00 0.00 0.00 0.00 0.000.00$金属-铜*EOS_GRUNEISEN3 0.394 1.49 0.00 0.00 2.02 0.470 0.00 0.00$Define Steel 4340*MAT_JOHNSON_COOK1, 7.8400E+00, 0.759E-00, 2.000E+00, 0.320E+00, 0.000E+000.793E-02, 0.510E-02, 2.60E-01, 1.400E-02, 1.030E+00, 1.793E+03, 2.930E+02, 1.000E-05 4.400E-06, 2.000E-02, 2.000E+00, 0.000E+00, 0.80E+00, 2.100E+00, -0.05E+00, 2.000E-03 0.610E+00*EOS_GRUNEISEN1, 4.578E-01, 1.330E+00, 0.000E+00, 0.000E+00, 1.670E+00, 0.430E+00, 0.000E+001.000E+00$ SI unit :cm-g-microsecond$ unit conversion factor$ 1 psi=6895 pa, 1 dyn=1e-5 n,1 bar=1e5 pa,1 g/cm**3=0.0361 lb/in**3,1 lb=0.454kg本参数援引University of Nevada Las Vegas。