LS-DYNA的状态方程模型

17. 状态方程模型

17.1 状态方程形式1Linear Polynomial

这个多项式状态方程,单位初始体积的内能呈线性,E

LS-DYNA的状态方程模型的图1 17.1.1

其中C0C1C2C3C4C5C6是用户定义常数。

LS-DYNA的状态方程模型的图2 17.1.2

V是相对体积,在膨胀单元中,的系数设为零,即:

LS-DYNA的状态方程模型的图3

线性多项式状态方程可用伽马定律状态方程来模拟气体。这可以通过设置来实现:

LS-DYNA的状态方程模型的图4

LS-DYNA的状态方程模型的图5

其中是比热的比率。压力则由下式给出:

LS-DYNA的状态方程模型的图6

请注意,E的单位是压力的单位

17.2 状态方程形式2JWL High Explosive

JWL状态方程将压力定义为相对体积,V,以及单位初始体积的内能E,的函数:

LS-DYNA的状态方程模型的图7 17.2.1

其中,ω、ABR1R2为用户定义的输入参数。这个状态方程通常用于在涉及金属加速度的应用中确定烈性炸药的爆轰产物的压力。该方程的输入参数由Dobratz [1981]给出了各种高爆炸材料的输入参数该状态方程与爆炸燃烧(材料模型8)材料模型一起使用,该模型决定了爆炸单元点火时间。

17.3 状态公式形式3Sack “Tuesday” High Explosives

爆轰产物的压力根据相对体积V单位初始体积的内能E给出,如[Woodruff 1973]

LS-DYNA的状态方程模型的图8 17.3.1

其中A1A2A3B1B2为用户定义的输入参数

该状态方程与爆炸燃烧(材料模型8)材料模型一起使用,该模型决定了爆炸单元点火时间。

17.4 状态方程形式4Gruneisen

具有立方激波速度-粒子速度的Gruneisen态方程将压缩材料的压力定义为

LS-DYNA的状态方程模型的图9 17.4.1

其中E单位初始体积的内能,Cvsvp剪切-压缩波速曲线的截距,S1S2S3vsvp曲线的斜率系数,γ0Gruneisen gammaa为对γ0的一阶体积修正。常数CS1S2S3γ0a均为输入参数。压缩根据相对体积V定义为:

LS-DYNA的状态方程模型的图10

对于膨胀材料,其压力的定义为:

LS-DYNA的状态方程模型的图11 17.4.2

17.5 状态方程形式5Ratio of Polynomials

多项式状态方程比率将压力定义为:

LS-DYNA的状态方程模型的图12 17.5.1

其中

LS-DYNA的状态方程模型的图13 17.5.2

LS-DYNA的状态方程模型的图14

LS-DYNA的状态方程模型的图15

在扩展区域中,F1F’1=F1+βµ2替代,常数Aijαβ是用户输入。

17.6 状态方程形式6Linear With Energy Deposition

这个多项式状态方程,单位初始体积,E,的内能呈线性,由

LS-DYNA的状态方程模型的图16 17.6.1

其中C0C1C2C3C4C5C6是用户定义常数。

LS-DYNA的状态方程模型的图17 17.6.2

V是相对体积,在膨胀单元中,的系数设为零,即:

LS-DYNA的状态方程模型的图18

内能E根据能量沉积速率随时间曲线的增加,其ID在输入中定义。

17.7状态方程形式7Ignition and Growth Model

JWL状态方程将未反应高爆物的压力定义为

LS-DYNA的状态方程模型的图19 17.7.1

其中Ve为相对体积,Ee为内能,常数AeBeωeR1eR2e为输入常数。同样地,反应产物中的压力也由另一种JWL形式来定义

LS-DYNA的状态方程模型的图20 17.7.2

未反应的炸药和反应产物的混合物由反应的分数F定义(F = 0表示没有反应,F = 1表示从炸药完全转化为产物)。假设压力和温度处于平衡状态,并假设相对体积为添加剂:

LS-DYNA的状态方程模型的图21 17.7.3

反应的速率被定义为

LS-DYNA的状态方程模型的图22 17.7.4

其中,IGHxyzm(通常是m = 0)是输入常数。

PBX-9404RX-03-BBPETN和铸造TNT四种炸药的JWL状态方程和反应速率拟合了一维和二维冲击起始和爆炸数据。计算方法的细节由CochranChan [1979]描述。LeeTarver [1980]给出了这四种炸药的详细的一维计算和参数。TarverHallquist [1981]讨论了PBX 9404LX- 17对该模型的二维计算。

17.8状态方程形式8Tabulated Compaction

压缩压力为正,容积应变εV为正。表格压缩模型与单位体积内能呈线性关系。压力定义为

LS-DYNA的状态方程模型的图23 17.8.1

增加加载(压缩)。卸载发生在峰值(最压缩)体积应变处的体积模量对应的斜率处,如图17.1所示。重新加载沿着卸载路径一直进行到开始卸载的点,然后继续按照式(17.8.1)描述的加载路径进行

17.9状态方程形式9Tabulated

表中的状态方程模型与内能呈线性关系。压力为

LS-DYNA的状态方程模型的图24 17.9.1

体积应变εV由相对体积的自然算法给出。在定义表格函数时,最多可以使用10个点,最少可以用2个点。如有必要,可以外插出该压力。与状态方程8不同,加载和卸载是沿同一曲线。

LS-DYNA的状态方程模型的图25

17.1.8状态方程的压力与体积应变曲线。在压实状态下,体积卸载模量取决于峰值体积应变。

17.9状态方程形式10Propellant-Deflagration

爆燃(燃烧速率)反应性流动模型需要一个未反应的固体状态方程、一个反应产物状态方程、一个反应速率定律和两种(或更多)物质的混合规则。标准点火和生长模型的混合规则[LeeTarver 1980]假设,随着反应的进行,压力和温度都是完全平衡的。然而,可以修改混合规则,不允许反应产物气体热传导或部分加热固体。对于这种相对缓慢的安全气囊推进剂燃烧过程,热平衡和压力平衡假设是有效的。目前在燃烧模型中使用的状态方程是JWLGruneisen、范德瓦尔斯共体积和完美气体定律,但其他的状态方程可以很容易地实现。在这种推进剂燃烧中,燃烧的叠氮化钠产生的气体氮在充满安全气囊的气中时,遵循完美的气体定律,但可能必须被建模为在推进室内产生的高压和温度下的范德瓦尔气体化学反应速率定律与压力、颗粒几何形状和表面积有关,以及大多数高压燃烧过程。当反应体系的温度分布已知时,可以使用与温度相关的阿伦尼乌斯化学动力学。

由于安全气囊推进剂的组成和性能数据是公司的私人信息,因此很难获得燃烧速率建模所需的信息。然而,帝国化学工业(ICI)公司提供了压力指数、粒子几何形状、堆积密度、反应热和大气压力燃烧率数据,这使我们能够开发他们的叠氮化钠+三氧化二铁驱动安全气囊推进剂的数值模型。[Hallquist等人,1990]描述了爆燃模型及其实现和ICI推进剂的结果。

未反应的推进剂和反应产物状态方程的形式均为:

LS-DYNA的状态方程模型的图26

其中p为压力(以Mbars表示),V为相对比体积(相对密度倒数),ωGruneisen系数,Cv为热容(以Mbars-cc/cc°K表示),T为温度°Kd为共体积,ABR1R2为常数。设置A=B=0。得到范德瓦尔的上体积状态方程。JWL状态方程通常在超过几千的压力下有用,而范德瓦尔方在压力低于该范围和高于完美气体定律成立的范围时有用。当然,设置A=B=d=0可以得到完美的气体定律。如果使用准确的ωCv值加上压缩和内能之间的正确分布,计算出的温度非常合理,因此可以用来检查推进剂的性能。

推进剂爆燃过程所用反应的反应速率如下:

LS-DYNA的状态方程模型的图27

其中F是反应的分数(F = 0意味着没有反应,F = 1是完整的反应),t是时间,p是压力(以Mbars),rsuwxyy、极限1和极限2是用于描述反应速率的压力依赖性和表面积依赖性的常数。如果推进剂是混合物或在某种压力或温度下表现出反应速率的急剧变化,则包括两个(或更多)与压力有关的反应速率。燃烧表面积的依赖关系可以用(1-Fy Fx项来近似。其他形式的反应速率定律,如与阿伦尼乌斯温度相关的e-E/RT型速率,也可以使用,但这些需要非常精确的温度计算。虽然LS-DYNA理论手册17.8压力的理论证明还不完整,但大量的实验燃烧速率与压力数据确实证明了这种效应,使用压力依赖燃烧的水动力计算准确地模拟了这种实验。

爆燃反应流模型被反应材料中一个或多个区域边界上的任何压力或粒子速度的增加所激活。这种增加在这些区域产生了压力,并开始分解。如果压力减轻,反应速率就会降低,并可以降为零。这一特征对于短时间的部分分解反应很重要。如果压力保持不变,反应的分数最终达到1,材料就完全转化为产物分子。用该模型计算出的几种推进剂通过推进剂的爆燃前推进速率与实验观察到的燃烧速率与压力曲线非常接近。

为了获得与实验爆燃数据的良好一致性,该模型需要准确地描述未反应的推进剂状态方程,要么对实验压缩数据进行分析拟合,要么基于以前使用类似材料的经验进行估计拟合。对于反应产物的状态方程也是如此。现有的实验燃烧速率、压力产生和能量传递数据越多,反应速率方程的形式和常数就越好。因此,在燃烧子例程中使用的压力方程

LS-DYNA的状态方程模型的图28

其中,VuTu分别为未反应的推进剂的相对体积和温度。相对密度显然是相对体积的倒数。反应产物中的压力Pp为:

LS-DYNA的状态方程模型的图29

随着反应的进行,假设未反应和产物压力和温度是平衡的(Tu=Tp=T p=Pu=Pp),相对体积是相加的:

LS-DYNA的状态方程模型的图30

其中,V为总相对体积。其他的混合假设可以并且已经在不同版本的DYNA2D/3D中使用过。反应速率定律的形式:

LS-DYNA的状态方程模型的图31

如果F超过FMXGR,则将GROW1项设为零,如果F小于FMNGR,则GROW2项为零。因此,两个独立(或重叠)的燃烧速率可以用来描述推进剂分解的速率。

这个状态方程子程序与一个材料模型一起用于描述推进剂。在安全气囊推进剂的情况下,可以使用零材料模型(type#10)。当确定剪切模量和屈服强度时,材料类型#10通常用于固体推进剂或炸药。推进剂材料由材料模型和反应开始前的无反应状态方程来定义。使用计算出的混合态,直到反应完成,然后使用反应产物的状态方程。假设反应热ENQ是一个常数,在F的所有值下都是相同的,但可以实现更复杂的能量释放定律。


参考文献:

[1] Hallquist J O. LS-DYNA Theory manual. 2006.





ANSYS LS-DYNA状态方程模型

LS-DYNA的状态方程模型的评论1条

LS-DYNA的状态方程模型的相关案例教程

1、数值计算模型 运用LS-DYNA有限元软件建立二维计算模型。模型中战斗部布置在空气介质中,水介质前端面为铝合金板;在铝合金板和间隔靶之间填充60cm厚的水介质,间隔靶由3层厚度均为1cm的45钢板组成,靶板间距为5cm,边缘采用全约束方式进行固定。采用Euler算法描述水、空气、炸药和药型罩,间隔靶、壳体和铝合金板采用Lagrange算法,Euler单元和Lagrange单元采用流固耦合算法进
*MAT_RHT 这是材料类型272。该模型用于分析受脉冲载荷影响的混凝土结构;参见Riedel et al.(1999)和Riedel(2004)。 卡片格式 卡片定义: MID——材料标识。必须指定一个不超过8个字符的唯一数字或标签。 RO ——质量密度 SHEAR——剪切模量 ONEMPA——单位
LS-DYNA是世界最著名的显式非线性动力学计算程序,广泛应用于冲击侵彻、高速碰撞、金属成形等固体结构非线性动力学分析。近年来,LS-DYNA 中加入了爆炸流场计算功能,可以用于计算爆炸容器的内部爆炸流场和壁面爆炸载荷。 0 1 爆轰模型 炸药的起爆和爆炸过程是一种快速的化学反应过程,对于该过程的描述,主要存在 CJ 模型和 ZND 模型两种。LS-DYNA 中包含两种炸药爆轰模型:高能燃烧模型和
1.LS-DYNA简介 LS-DYNA 是世界上最著名的通用显式动力分析程序,能够模拟真实世界的各种复杂问题,特别适合求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成型等非线性动力冲击问题,同时可以求解传热、流体及流固耦合问题。在工程应用领域被广泛认可为最佳的分析软件包。与实验的无数次对比证实了其计算的可靠性。 由J.O.Hallquist主持开发完成的DYNA程序系列被公认为是显式有限元程序
1.LS-DYNA简介 LS-DYNA 是世界上最著名的通用显式动力分析程序,能够模拟真实世界的各种复杂问题,特别适合求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成型等非线性动力冲击问题,同时可以求解传热、流体及流固耦合问题。在工程应用领域被广泛认可为最佳的分析软件包。与实验的无数次对比证实了其计算的可靠性。 由J.O.Hallquist主持开发完成的DYNA程序系列被公认为是显式有限元程序
博士/博士后
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    1 13