圆锥罩聚能射流LSDYNA仿真计算(图解)
说明
LS-DYNA 是世界上著名的动力学分析有限元程序，可以精确可靠地处 理各种高度非线性问题，如碰撞分析、爆炸分析、冲压成型分析和常规武 器设计等。LS-DYNA 以 Lagrange 算法为主，兼有 ALE 和 Euler 算法；以 显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体 —结构耦合功能；以非线性动力分析为主，兼有静力分析功能；它计算的 可靠性已经被无数次试验所证明。
本文主要对涉及到聚能射流的 LS-DYNA 显式动力学分析准备内容学 习，并附有圆锥罩聚能射流的二维模拟实例，期间介绍了相关的内容与知 识点，如建模、分网和关键字的含义等。
由于作者水平有限，文中错误在所难免，敬请各位朋友批评指正！ 文中特殊提示部分用蓝色进行标识；表格部分用深蓝色标识；注释部 分用红色进行标识。
CP, PC, SPALL, IT, D1, D2, D3, D4
D5
其中：
MID 为材料模型的 ID 号；
RO 为材料的密度；
G 为材料的剪切模量；
E 为材料的杨氏模量（仅用于壳单元）；
PR 为材料的泊松比（仅用于壳单元）；
TM 为熔化温度；
TR 为室内温度；
CP 为材料比热；
SPALL 为材料的破碎形式；
QQ：759195361
第 8 页 共 44 页 Email：qingzhouwang@
ANSYS/LS-DYNA
材料模型的级联菜单，每个菜单选项对应一种材料模型。图示如下：
对于每种单元类型，未必能够使用所有的材料模型，每种材料模型并不是适合所 有的单元类型，因此使用时要查阅单元手册来确认可以用哪种模型。
QQ：759195361
第 1 页 共 44 页 Email：qingzhouwang@
ANSYS/LS-DYNA
第1部分 基础知识
本部分涉及定义单元类型、材料属性、实体模型的建立、有限元模型 的建立、PART 的定义、求解设置、输出 K 文件、求解和 LS-PREPOST 等聚能射流分析用学习资料。
SOLID164 是用于 3 维的结构实体单元，由 8 节点构成。下图描述了 SOLID164 几何特性、节点位置和坐标系。这个单元只用在动力分析，它支持所有许可的非线性 特性。
使用这种单元时应该注意以下几点： (1) 这个单元必须有 8 个节点。它可以退化为楔形、四面体和棱 锥形。但是一个 0 体积单元是不允许的。单元不可以被扭曲以至 有两个分离的体积，当单元没有正确编号时，这是经常会发生的。 (2) 该单元可供设置的选项包括：单元算法、单元连续性特征。缺省情况下，SOLID164 实体单元是用一点的积分加上粘性沙漏控制来加快单元的方程式。可以利用 Lagrangian(缺省)和 Arbitrary Lagrangian-Eulerian(ALE)网格。为了使 ALE 方程启用， (3) 除了设置 KEYOPT(5)=1 外，还必须用 EDALE 和 EDGCALE 设置正确的参数。 (4) 该单元没有实常数。 (5) 压力能通过对单元表面施加面荷载，正的压力施加在单元上。 定义单元类型及其算法的菜单操作步骤如下： 选择 Main Menu > Preprocessor > Element Type > Add/Edit/Delete 命令，弹 出 Element Type 对话框。单击 Add… 按钮，弹出 Library of Element Types 对话框。
ANSYS/LS-DYNA
选择 Main Menu > Preprocessor > Material Props > Material Models 命令，弹 出 Define Material Model Behavior 对话框。左侧框图中，显示为材料模型代号，有多 种材料时，点击菜单 Material > New Model… 命令，弹出 Define Material ID 对话框， 输入材料 ID，即可新建材料。右侧框图的 Material Models Available 选择栏对应的是
I2 部分材料的 MAT_JOHNSON_COOK 参数
材料
洛氏 密度 比热 熔温 A
B
硬度 g/cm3 J/kg*K K MPa MPa N
C
M
高导无氧铜 F-30 8.9
383 1356 90 292 0.31 0.025 1.09
药筒黄铜
使用 SHELL163 单元时应该注意以下问题： （１）SHELL163 提供了 12 种算法； （２）积分规则（常用于标准积分规则，即积分规则标识为标准积分），它包括：高 斯规则（最多允许五个积分点）、梯形规则（最多允许 100 个积分点）两种规则。 （３）标记层复合材料模式，它有两种模式，即非复合材料模式和复合材料模式。
QQ：759195361
第 7 页 共 44 页 Email：qingzhouwang@
ANSYS/LS-DYNA
4340 钢
C-30 7.83 477 1793 792 510 0.26 0.014 1.03
*MAT_NULL
//注：该关键字用于空气参数设定
MID, RO, PC, MU, TEROD, CEROD, YM, PR MID 为材料模型 ID 号； RO 为材料的密度； PC 为材料的剪切应力； MU 为材料的泊松比； TEROD 为材料扩张侵蚀过程的相对体积，取 0 表示被忽略； CEROD 为材料压缩侵蚀过程的相对体积，取 0 表示被忽略； YM 为材料的杨氏模量； PR 为材料的动力粘性系数。 表 3 部分材料的 MAT_NULL 参数
*MAT_HIGH_EXPLOSIVE_BURN //注：该关键字用于炸药参数设定，须通过编辑、修改 DYNA 的数据输入文件 Jobname.K 来进行加入。
MID, RO, D, PCJ, BETA, K, G, SIGY
其中：
MID 为材料模型的 ID 号；
RO 为材料的密度；
D 为材料的爆速；
QQ：759195361
第 2 页 共 44 页 Email：qingzhouwang@
ANSYS/LS-DYNA
1.1 定义单元类型
单元类型及其算法的选择、实常数的定义所对应的菜单操作和 ANSYS 命令流代 码，对各种类型的分析几乎都是同一种模式。
LSDYNA 中可供选择的单元有： 1. LINK160 杆单元； 2. BEAM161 梁单元； 3. PLANE162 平面实体单元； 4. SHELL163 薄壳单元； 5. SOLID164 实体单元； 6. COMBIN165 弹簧阻尼单元； 7. MASS166 质量单元； 8. LINK167 仅拉伸单元（索单元）； 9. SOLID168 10 节点四面体单元。 本次中仅用到 SHELL163 薄壳单元和 SOLID164 实体单元。 SHELL163 是一个 4 节点显式结构薄壳单元，有弯曲和膜特征，可加平面和法向 载荷。单元在每个节点上有 12 个自由度：在节点 X、Y 和 Z 方向的平动，加速度， 速度和绕 X、Y 和 Z 轴的转动。该单元支持显式动力学分析所有非线性特性。
ANSYS/LS-DYNA
1.2 定义材料属性
材料属性指材料的杨氏模量、泊松比、密度等参数。 本文中用到的材料有： *MAT_HIGH_EXPLOSIVE_BURN； *MAT_JOHNSON_COOK； *MAT_NULL； *MAT_PLASTIC_KINEMATIC； 材料的关键字介绍如下：
QQ：759195361
第 3 页 共 44 页 Email：qingzhouwang@
ANSYS/LS-DYNA
（４）积分规则标识：它包括标准积分选项和用户定义积分规则选项。 （５）SHELL163 提供了如下实常数：SHRF 是剪切因数，推荐 5/6，如果不指定则默 认为 1；NIP 是通过单元厚度的积分点数值，最大值为 100。如果 NIP 输入值为 0 或 空，ANSYS 会默认积分值为 2；T1~T4 是四个节点中每个节点处的壳厚度。
特别提示 在定义材料属性时，确保使用一致的单位制，不正确的 单位制不仅会影响材料响应，而且影响最终的计算结果。 不要低估准确材料数据对结果的重要性，尽量花费额外 的时间和金钱去获得更为准确的材料数据。
QQ：759195361
第 9 页 共 44 页 Email：qingzhouwang@
材料
RO
PC
MU TEROD CEROD YM
PR
水
1
-1e-4
1e-6
0
0
0
0
1.225e-3
0
17.456e-6
0
0
0
0
空气
1.28e-3
0
0
0
0
*MAT_PLASTIC_KINEMATIC
//注：该关键字用于靶板参数设定
MID, RO, E, PR, SIGY, ETAN, BETA
SRC, SRP, FS, VP MID 为材料模型的 ID 号； RO 为材料的密度； E 为材料的杨氏模量； PR 为材料的泊松比； SIGY 为材料的屈服强度。 定义材料的菜单操作步骤如下：
-5.6
RDX
1.00
6.048 5.98
+1.1
-
-
-
QQ：759195361
第 6 页 共 44 页 Email：qingzhouwang@
ANSYS/LS-DYNA
φ=6.784
1.20
6.731 6.750 -0.3 152.2 152
+0.1
1.59
8.063 8.100 -0.5 267.1 287
-6.9
HMX
1.9
φ=6.772
基于LS-DYNA的智能植保车喷雾仿真分析
|四川农业与农机/2023年5期|>>>我国是世界上果园种植面积最大的国家[1]，随着人民生活水平的日益提高，鲜果需求量极大。
据国家统计局数据，我国果园面积呈逐年增长趋势，2021年，我国果园种植面积达1280.80万hm 2[2]，预计2023年将达1351.06万hm 2[3]。
我国丘陵山区果园约占果园总面积的65%[4]，其主要种植方式为多植株密集型果园种植模型[5]。
高效、精细的果园管理是水果增产、增收的关键，病虫害防治是当前果园植保的关键环节，其中，化学防治是果园植保病虫害防治的主要方法之一[6]。
但过度使用化学农药，会造成环境污染和农药残留，危害人体健康。
为有效减少化学农药的使用，提高农药喷洒效率，降低人工成本，实现农药精准喷洒，本文研制了新能源农业智能植保车，并对智能植保车喷雾进行了仿真分析，以期实现化学农药的精确喷洒，降低农药使用量，减少环境污染，提高生产效率。
在无人植保车喷雾性能方面，国内外学者进行了大量科学研究，主要研究方法为试验法和计算流体动力学方法[7]。
石河子大学李瑞敏等通过试验方法，分析了扇形雾喷头结构参数对喷雾性能的影响[8]；JING D 、ZREMBA M 等通过流体动力学计算分析了喷嘴雾化特性[9-10]；聂涛等分析了压力旋流喷嘴雾化特性[11]，蒋小平等通过流体动力学计算分析了喷射参数对扇形喷嘴雾化特性的影响[7]。
试验法周期长、成本大，计算流体动力学方法求解时间长，资源消耗大。
基于以上分析，本文运用时空守恒元/解元方法，基于LS-DYNA 对无人植保车喷雾仿真流固耦合模型进行分析，从而获得智能无人植保车喷雾特性，以期为研究无人植保车喷雾结构及优化参数提供理论支持。
1智能植保车结构根据果园植保管理要求，智能植保车的结构设计综合考虑了果园种植行间距、丘陵山区地形坡度等因素，具有药液载量大、续航久、动力足、通过性能好以及爬坡能力强等特点。
智能植保车主要由动力传动系统、液压驱动系统、喷药机构以及履带式行走机构组成，如图1所示。
聚能射流形成及破甲过程的数值模拟分析
23 期
韩秀清, 等: 聚能射流形成及破甲过程的数值模拟分析
69 63
69 64
科学技术与工程
9卷
以上各图的标度条为图中各单元的应力大小, 以不同颜色清晰的表示。最上边的红色为最大值, 最下边的蓝色为最小值, 中间从下到上表示的应力 值逐渐增大。
3 炸高对聚能射流侵彻作用的影响
石油射孔弹药型 罩底径端面到靶 板之间的垂 直距离, 即所谓的炸高, 影响着聚能射流对靶板 ( 这 里是射孔枪壁 ) 的侵彻深度, 也就是 破甲深度。作 为算例, 在前述模型数据的基础上, 使用实验参数, 以期找到炸高这一参数 对侵彻深度的 影响。具体 模拟结果如表 1所示。
由于爆炸实 验费 用昂 贵, 试验 过程也 比较 漫 长和复杂, 测 量手段 和观 测条件 在很大 程度 上限 制了爆破的试验研究。利用流体力学和 动力学模 型来描述射流 的过程 需要综 合化学 反应 方程、反 应率方程、热 传导 方程、材料 本构关 系等 方程, 这 些方程只有在 很多假 设、极其简 化的情 况下 才能
1. 3 材料模型 在聚能装药射流形成及破甲的数值模拟中, 涉
及到炸药、药型罩和射孔枪 壁三种材料模 型, 其具 体参数在 K 文件里可显示。 1. 3. 1 炸药
程序提供了 MAT _H IGH _EXPLOS IVE _BURN 高 能炸药材 料模型, 爆 轰产物 压力 与体积 关系 采用 JW L 状态方程, 其一般形式为 [ 3] :
侵彻深度。因此, 最佳炸高为 22. 92 mm。
4 结语
运用大型有限元软件做前处理, 利用通用显式动 力分析程序对圆锥型聚能装药射流进行了数值模拟。 这样可以在不引入其他干扰因素的情况下, 对侵彻过 程进行计算, 遇到特殊点还可以任意跟踪和放大观 察[ 5] , 从而得到试验无法得到的数据。数值计算结果 符合石油射孔弹聚能装药射流的物理现象和规律, 存 在一定程度的误差, 但在工程允许的范围之内。
LS-DYNA精编教程
LS-DYNA软件1.1 LS-DYNA 简介LS-DYNA 是世界上最著名的通用显式动力分析程序，能够模拟真实世界的各种复杂问题，特别适合求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成型等非线性动力冲击问题，同时可以求解传热、流体及流固耦合问题。
在工程应用领域被广泛认可为最佳的分析软件包。
与实验的无数次对比证实了其计算的可靠性。
由J.O.Hallquist主持开发完成的DYNA程序系列被公认为是显式有限元程序的鼻祖和理论先导，是目前所有显式求解程序（包括显式板成型程序）的基础代码。
1988年J.O.Hallquist创建LSTC公司，推出LS-DYNA程序系列，并于1997年将LS-DYNA2D、LS-DYNA3D、LS-TOPAZ2D、LS-TOPAZ3D等程序合成一个软件包，称为LS-DYNA。
LS-DYNA的最新版本是2004年8月推出的970版。
1.1.1 LS-DYNA功能特点LS-DYNA程序是功能齐全的几何非线性（大位移、大转动和大应变）、材料非线性（140多种材料动态模型）和接触非线性（50多种）程序。
它以Lagrange 算法为主，兼有ALE和Euler算法；以显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体-结构耦合功能；以非线性动力分析为主，兼有静力分析功能（如动力分析前的预应力计算和薄板冲压成型后的回弹计算）；军用和民用相结合的通用结构分析非线性有限元程序。
LS-DYNA功能特点如下：1.分析能力：z非线性动力学分析z多刚体动力学分析z准静态分析（钣金成型等）z热分析z结构-热耦合分析z流体分析：欧拉方式任意拉格郎日-欧拉（ALE）流体-结构相互作用不可压缩流体CFD分析z有限元-多刚体动力学耦合分析（MADYMO，CAL3D）z水下冲击z失效分析z裂纹扩展分析z实时声场分析z设计优化z隐式回弹z多物理场耦合分析z自适应网格重划z并行处理（SMP和MPP）2.材料模式库(140多种)z金属z塑料z玻璃z泡沫z编制品z橡胶（人造橡胶）z蜂窝材料z复合材料z混凝土和土壤z炸药z推进剂z粘性流体z用户自定义材料3.单元库z体单元z薄/厚壳单元z梁单元z焊接单元z离散单元z束和索单元z安全带单元z节点质量单元z SPH单元4.接触方式(50多种) z柔体对柔体接触z柔体对刚体接触z刚体对刚体接触z边-边接触z侵蚀接触z充气模型z约束面z刚墙面z拉延筋5.汽车行业的专门功能 z安全带z滑环z预紧器z牵引器z传感器z加速计z气囊z混合III型假人模型6.初始条件、载荷和约束功能z初始速度、初应力、初应变、初始动量（模拟脉冲载荷）；z高能炸药起爆；z节点载荷、压力载荷、体力载荷、热载荷、重力载荷；z循环约束、对称约束（带失效）、无反射边界；z给定节点运动（速度、加速度或位移）、节点约束；z铆接、焊接（点焊、对焊、角焊）；z二个刚性体之间的连接－球形连接、旋转连接、柱形连接、平面连接、万向连接、平移连接；z位移/转动之间的线性约束、壳单元边与固体单元之间的固连；z带失效的节点固连。
聚能射流对靶板侵彻的数值仿真
兵 工 自 动 化
Ｏｒ ｄｎａｎｃ ｎｄｕｓ ｒ ｕｔ ｅＩ ｔｙ Ａ ｏｍ ａｔｏｎ ｉ
２０１ ０ｌ ０．
２ （） ９ １
ｄ ｉ １．９ ９ｊ ｓｎ１０ ・５ ６２ １．１ ０ ｏ： ０３ ６ ／ｉｓ．０ ６１ ７ ．０ Ｏ０ ． ７ ． ０
聚能射 流对靶 板 侵彻 的数 值仿 真
崔 军 ， 徐 峰 ， 李 向 荣 （ 甲 兵 工 程 学 院 兵 器 工 程 系 ， 北 京 １０ ７ ） 装 ００ ２
摘要 ：利 用 ＡＵＴ ＤＹＮ非 线性 动 力学分析 软 件 ，采 用 Ｊ ｈ ｓｎＣ ｏ 动 态本 构模 型和 多物 质 Ｅ ｌｒ算法 ，对 某子 Ｏ ｏ ｎｏ ｏ ｋ ｕｅ 母 弹的 子弹 聚 能射 流 形成 、侵 彻 钢板 过程 进行 了数 值 仿真 ，得到 与 试验 结 果相 近的 聚 能射 流形 成和 侵彻 的 物理 现 象 和规律 ，验 证 了该 模 型和数 值 模拟 的合 理 性 ，为 该弹 的测 试 和鉴 定提 供 一些 理论依 据 ，对 聚能破 甲战斗部 的工程设
ANSYSLS-DYNA数值模拟霍普金森压杆试验资料
A N S Y S L S-D Y N A数值模拟霍普金森压杆试验ANSYS/LS-DYNA数值模拟霍普金森压杆试验1 功能概述大多数材料在强度等力学性质方面都表现出某种程度的加载率或应变率敏感性，高幅值短持续时间脉冲和荷载所引起材料力学性质的应变率效应，对于抗动载的结构设计和分析是非常重要的。
这些动载来至常规武器侵彻与爆炸、偶然爆炸和高速撞击等许多军事和民用事件，对于这些事件的理论分析和数值模拟必须知道材料的高应变率强度、断裂特性和应力-应变关系等本构性质。
要研究材料在脉冲动载作用下的力学性质的实验设备和实验必须模拟类似现场的应变率条件，分离式霍普金森杆被公认为是最常用最有效的研究脉冲动载作用下材料力学性质的实验设备。
数值模拟是一种依靠电子计算机对工程问题和物理问题乃至自然界各类问题进行研究的技术。
它利用材料的本构函数，结合有限元或有限容积的概念，采用数值计算和图像显示的方法，因此具有如下优势：（1）检验理论结果是否正确；（2）弥补实验与观测得不足；（3）利用模拟结果，了解非线性过程中的因果关系与主要物理机制；（4）预测在不同初始条件与边界条件下非线性过程的发展情形；（5）数值模拟成本低，可以带来巨大社会经济效益。
由于很多材料的本构性质已经知道，因此在设计产品时，可以利用材料的本构性质通过仿真来模拟复杂的系统。
ANSYS/LS-DYNA数值模拟霍普金森压杆试验，就是通过ANSYS/LS-DYNA软件来模拟霍普金森压杆实验，通过设置弹丸不同速度，对试件进行研究。
霍普金森压杆实验分为自由式和分离式两种，本仿真采用分离式的办法。
2 原理简介2.1 霍普金森压杆实验简介霍普金森杆实验装置的基本原型最早是由Hopkinson提出的，它可用于测量冲击载荷的脉冲波形。
1949年Kolsky将压杆分成两段，试件置于输入杆和输出杆中间，通过加速的质量块、短杆撞击或炸药爆轰产生加速脉冲，利用这一装置可测量材料在冲击载荷作用下的应力-应变关系。
多模成型装药的药型罩结构设计
多模成型装药的药型罩结构设计
多模成型装药是指根据不同的毁伤目标而自适应选择不同模式毁伤元,如大炸高小穿深的爆炸成型弹丸(Explosive Formation Penetration—简称EFP)、具有大穿深的聚能高速金属射流和聚能杆式侵彻体(国外称之为JPC, Jetting Projectile Charge)。
为了使得同一弹药能对付战场中出现的多种目标,多模成型装药成为当前成型装药技术研究的热点问题之一。
本文主要针对多模成型装药研究的关键技术之一——药型罩结构参数与多模毁伤元的匹配关系研究。
基于成型装药基本理论及多模毁伤元形成理论,采用LS-DYNA仿真软件,数值仿真研究了成型装药毁伤元形成过程,得到了次口径药型罩结构参数与毁伤元成型参数之间的变化规律曲线。
建立多模成型装药结构模型,选取药型罩结构参数曲率半径和壁厚仿真计算方案,分析计算结果,并进行了Matlab三维作图,找出了形成多模毁伤元的药型罩结构参数范围,当药型罩曲率半径为0.7～0.9倍装药口径,药型罩壁厚为
0.40～0.56倍装药口径时,此时所形成的多模毁伤元EFP和JPC均满足课题试验项目的指标要求。
最后,进行了多模毁伤元成型及侵彻的试验研究,试验结果与仿真结果较一致,为今后进一步研究多模战斗部提供了参考依据。
基于LS-DYNA的弹体入水动力学仿真研究
摘
要： 运 用动 力学有限元软件 Ｌ Ｓ ． Ｄ ＹＮ Ａ，采用 Ａ Ｌ Ｅ算法分别对平头弹体和 圆头弹体在速度 ７ ５ ． ４～ １ ４ ２ ． ７ ｍ／ ｓ 范围
入水 角 ３ Ｏ 。 、６ ０ 。 和９ ０ 。 环境下进行 了数值模拟 ，得到 了两种头型弹体在入水初期 的弹道轨迹和 空泡 形状 、 自由液面变化 及 飞溅效应 。研 究 了头型及 入水角度 对弹体轨迹、速度 的影响 ，分 析了入水压力峰值随速度的变化 ，给 出了速度 与阻力 系数 的关 系式。通 过对 比实验数据，模型 的计算结果和实验结果 吻合得较好 ，为入水弹体头型设计 、抗冲击 设计 和材 料 选择提供 了有价值 的参考 。 关键 词：流固耦合 ：入 水；数值 模拟 ；ＡＬ Ｅ算法 中图分类号：Ｔ Ｊ ６ 文献标志码 ：Ａ 中图分类号 ：１ ６ ７ ２ － ４ ８ ０ １ （ ２ ０ １ ３ ） ０ ５ － ０ １ ６ － ０ ３ ２
２ ． ０
０ 一 。 ７ ～ 一 １ ｊ ￣ 一 ＋ 一 ９ 晕 ５ ｓ ． ４ ｐ ｍ ／ ｓ Ｑ ｐ Ｉ 一 ≥ ～ ｊ … ≮ 一 。 ‘ ｙ … … … ｔ 一 一 － ’ － ． ／ ＇ ／ 。 ～ ｊ ｒ 一 ’ ‘
本 文 水 和 空气 的状 态 方 程 都 选 用 Ｇ ｒ ｕ ｎ ｅ ｉ ｓ ｅ ｎ 状 态方 程 。该 方程 定 义 空气域 和 水域 压 缩材 料 的 压 力为 ：
作 者 简介 ：宋 盼盼 （ １ ９ ８ ６ 一） ，女 ，在读 硕 士研 究生 ，主 要研 究方 向为 入水 动力 学 。
１ ６
机 电技术
２ ０ １ ３ 年ｌ ０ 月
基于 Ｌ Ｓ － Ｄ Ｙ ＮＡ的弹体入水 动力学仿 真研 究
宋盼盼 赵捍 东 吴建萍２
射流侵彻混凝土介质数值模拟及影响因素研究_郑平泰
纲应变率 , E为实际应变率, E 0 = 1 1 0 s- 1 为参照应
116
弹 箭与 制 导学 报
- 3 3
2007ቤተ መጻሕፍቲ ባይዱ年
变率 ; A 、 B、 C、 N 为材料常数, A 为规一化内聚强 度, B 为规一化压力硬化系数 , N 为压力硬化指 数, C 为应变率系数。 H JC 模型关于混凝土材料的损伤因子 D( 0 [ D [ 1 ) 由等效塑性应变和塑性体积应变累加 得到 , D 的表达式为 : D=
结果与实验数据基本吻合。通过数值模拟计算重点分析了药型罩锥角、 药型罩厚度、 壳体厚度、 炸高和装药长 度对聚能射流侵彻效能的影响规律。研究结果表明 : 药 型罩锥角、 药型罩 厚度和炸高 对聚能射 流的侵彻 深度 和侵彻孔径影响较大 , 而壳体厚度和装药长度对 侵彻效能影响较小。 关键词 : 聚能射流 ; 侵彻 ; 混凝土介质 ; 数值模拟 ; 影响因素 中图分类号 : TJ4101 6 文献标志码 : A
为: 密度 : 10 g/ mm , 强度或屈服应力 : M Pa , 杨 氏模量 E: GPa , 剪切模量 G: GPa , 体积模量 K : GPa , K 1 ~ K 3 : GP a, k: 10 # m m , 其它参数的
12 -3
量纲均由基本物理量导出。
表3 Q 0 G 14. 86 P crush 16 C 0 . 007 混凝土材料模型参数 f cs 44 E v, crush 0 . 001 EF min 0 . 01 S max 7. 0 E v , l ock 0. 1 N 0 . 61 T 4. 1 P l ock 800 K3 208 A 0. 75 K1 85 B 1 . 65 K2 - 171
紧凑型聚能装药射流成型的数值模拟
紧凑型聚能装药射流成型的数值模拟姚志华;王志军;王向东;李德战;刘洋【期刊名称】《弹箭与制导学报》【年(卷),期】2012(032)005【摘要】为研究药型罩结构参数对K装药射流成型的影响,利用ANSYS/LS-DYNA软件对其成型过程进行数值模拟,得到了结构参数对射流的头部速度、长度、动能和能量利用率的影响规律,并对有无波形控制器的情况进行了分析.研究表明:药型罩高度40 mm ＜H＜60 mm、曲率半径300mm＜R＜500mm、厚度3mm＜δ＜5 mm、装药长径比0.7＜λ＜0.9时,形成射流性能较好.与无波形调整器的装药结构相比,射流的能量利用率提高17.6％,速度提高42.5％,长度增加1倍,能量增加2.5倍.【总页数】4页(P79-82)【作者】姚志华;王志军;王向东;李德战;刘洋【作者单位】中北大学机电工程学院,太原030051;空军场务技术试验中心,山东济宁272000;中北大学机电工程学院,太原030051;晋西工业集团有限责任公司,太原030027;空军场务技术试验中心,山东济宁272000;空军场务技术试验中心,山东济宁272000【正文语种】中文【中图分类】TJ410.34【相关文献】1.椭圆形罩线型聚能装药射流成型过程数值模拟研究 [J], 武双章;顾文彬2.装药长径比对半球形聚能装药射流成型的数值模拟 [J], 张万君;李国辉;王凯琳;吴晓颖3.环形聚能装药尼龙射流成型及侵彻过程数值模拟研究 [J], 吉庆;王志军;伊建亚;董理嬴;汤雪志4.环形聚能装药尼龙射流成型及侵彻过程数值模拟研究 [J], 吉庆;王志军;伊建亚;董理嬴;汤雪志5.新型聚能装药爆炸成型杆式射流数值模拟及试验研究 [J], 徐文龙;王成;徐斌;胡平因版权原因，仅展示原文概要，查看原文内容请购买。
第十七章 LS-DYNA 的隐式求解讲解学习
第十七章L S-D Y N A 的隐式求解第十七章 LS-DYNA 的隐式求解LS-DYNA 作为著名的显式求解程序只能求解瞬态动力问题，对于时间历程较长的静、动力问题， LS-DYNA 的显式中心差分法有它的局限性，而一些与瞬态动力分析紧密相关的问题要求LS-DYNA 也能够求解，如：➢ 冲压成型后的回弹计算➢ 应力初始化➢ 冲击后常时间低频动力响应➢ 静力分析➢ 特征值分析实际上从950版本开始，LS-DYNA 已增加了隐式求解功能。
刚开始的应用主要在冲压成型后的回弹计算，经过960版到970版的发展，LS-DYNA 的隐式求解功能已大增强，已经能满足以上的求解需要。
17.1显式与隐式的区别：17.1.1 LS-DYNA 显式求解：采用中心差分方法进行显式时间积分int n ext n n f fMa -=● 方程非耦合，可以直接求解(显式)● 但需要常小的时间步保持稳定状态● 不需要求解刚度矩阵 ● 适合冲击、穿透等高频非线性动力响应问题17.1.2 LS-DYNA 隐式求解：采用Newmark 隐式时间积分n n ext n n n Ma f fu K a M --=∆+∆++int 11● 对于线性问题，无条件稳定● 可以采用大的时间步 ● 对于非线性问题，需要一系列线性逼近（Newton-Raphson ）叠代求解 ● 需要求解刚度矩阵● 适合静力问题、低频动力问题及特征值分析。
17.2 LS-DYNA 中隐式分析的激活及相关关键字在LS-DYNA 中，缺省的求解是显式求解，可以通过下面的关键字来激活隐式求解：*CONTROL_IMPLICIT_GENERAL*CONTROL_IMPLICIT_GENERAL$ imflag dt0 iefs nstepsb igso1 0.01 0 0 0其中参数imflag=1激活全隐式求解imflag=0为缺省的显式求解。
imflag=2为显式求解后无缝进行隐式回弹求解。