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。
1.1.1 LS-DYNA 功能特点
LS-DYNA 程序是功能齐全的几何非线性(大位移、大转动和大应变)、材料非线性(140 多种材料动态模型)和接触非线性(50 多种)程序。它以 Lagrange 算法为主,兼有 ALE 和 Euler 算法;以显式求解为主,兼有隐式求解功能;以结构分析为主,兼有热分析、流体-结构耦合功能;以非线性动力分析为主,兼有静力分析功能(如动力分析前的预应力计算和薄板冲压成型后的回弹计算);军用和民用相结合的通用结构分析非线性有限元程序。LS-DYNA 功能特点如下:
1.分析能力:
(1) 非线性动力学分析
(2) 多刚体动力学分析
(3) 准静态分析(钣金成型等)
(4) 热分析
(5) 结构-热耦合分析
(6) 流体分析:
a. 欧拉方式
b. 任意拉格郎日-欧拉(ALE)
c. 流体-结构相互作用
d. 不可压缩流体CFD分析
(7) 有限元-多刚体动力学耦合分析 (MADYMO,CAL3D)
(8) 水下冲击
(9) 失效分析
(10) 裂纹扩展分析
(11) 实时声场分析
(12) 设计优化
(13) 隐式回弹
(14) 多物理场耦合分析
(15) 自适应网格重划
(16) 并行处理(SMP 和 MPP)
2.材料模式库(140 多种)
(1) 金属
(2) 塑料
(3) 玻璃
(4) 泡沫
(5) 编制品
(6) 橡胶(人造橡胶)
(7) 蜂窝材料
(8) 复合材料
(9) 混凝土和土壤
(10) 炸药
(11) 推进剂
(12) 粘性流体
(13) 用户自定义材料
3.单元库
(1) 体单元
(2) 薄/厚壳单元
(3) 梁单元
(4) 焊接单元
(5) 离散单元
(6) 束和索单元
(7) 安全带单元
(8) 节点质量单元
(9) SPH 单元
4.接触方式(50 多种)
(1) 柔体对柔体接触
(2) 柔体对刚体接触
(3) 刚体对刚体接触
(4) 边-边接触
(5) 侵蚀接触
(6) 充气模型
(7) 约束面
(8) 刚墙面
(9) 拉延筋
5.汽车行业的专门功能
(1) 安全带
(2) 滑环
(3) 预紧器
(4) 牵引器
(5) 传感器
(6) 加速计
(7) 气囊
(8) 混合III型假人模型
6.初始条件、载荷和约束功能
(1) 初始速度、初应力、初应变、初始动量(模拟脉冲载荷);
(2) 高能炸药起爆;
(3) 节点载荷、压力载荷、体力载荷、热载荷、重力载荷;
(4) 循环约束、对称约束(带失效)、无反射边界;
(5) 给定节点运动(速度、加速度或位移)、节点约束;
(6) 铆接、焊接(点焊、对焊、角焊);
(7) 二个刚性体之间的连接-球形连接、旋转连接、柱形连接、平面连接、万向连接、平移连接;
(8) 位移/转动之间的线性约束、壳单元边与固体单元之间的固连;
(9) 带失效的节点固连。
7.自适应网格剖分功能
自适应网格剖分技术通常用于薄板冲压变形模拟、薄壁结构受压屈曲、三维锻压问题等大变形情况,使弯曲变形严重的区域皱纹更加清晰准确。
对于三维锻压问题,LS-DYNA 主要有两种方法:自适应网格剖分和任意拉格朗日-欧拉网格(ALE)网格进行 Rezoning),三维自适应网格剖分采用的是四面体单元。
8 ALE 和 Euler 列式
ALE 列式和 Euler 列式可以克服单元严重畸变引起的数值计算困难,并实现流体-固体耦合的动态分析。在 LS-DYNA 程序中 ALE 和 Euler 列式有以下功能:
(1) 多物质的 Euler 单元,可达 20 种材料;
(2) 若干种 Smoothing 算法选项;
(3) 一阶和二阶精度的输运算法;
(4) 空白材料;
(5) Euler 边界条件:滑动或附着条件;
(6) 声学压力算法;
(7) 与 Lagrange 列式的薄壳单元、实体单元和梁单元的自动耦合。
9.SPH 算法
SPH(Smoothed Particle Hydrodynamics)光顺质点流体动力算法是一种无网格 Lagrange 算法,最早用于模拟天体物理问题,后来发现解决其它物理问题也是非常有用的工具,如连续体结构的解体、碎裂、固体的层裂、脆性断裂等。SPH 算法可以解决许多常用算法解决不了的问题,是一种非常简单方便的解决动力学问题的研究方法。由于它是无网格的,它可以用于研究很大的不规则结构。
SPH 算法适用于超高速碰撞、靶板贯穿等过程的计算模拟。
10.边界元法
LS-DYNA 程序采用边界元法 BEM(Boundary Element Method)求解流体绕刚体或变形体的稳态或瞬态流动,该算法限于非粘性和不可压缩的附着流动。
11.隐式求解
用于非线性结构静动力分析,包括结构固有频率和振型计算。LS-DYNA 中可以交替使用隐式求解和显式求解,进行薄板冲压成型的回弹计算、结构动力分析之前施加预应力等。
12.热分析
LS-DYNA 程序有二维和三维热分析模块,可以独立运算,也可以与结构分析耦合,可进行稳态热分析,也可进行瞬态热分析,用于非线性热传导、静电场分析和渗流计算。
热传导单元:8 节点六面体单元(3D),4 节点四边形单元(2D);
材料类型:各向同性、正交异性热传导材料,可以与温度相关,以及各向同性热传导材料的相变;
边界条件:给定热流 flux 边界,对流 convection 边界,辐射 radiation 边界,以及给定温度边界,它们可随时间变化;给定初始温度,可计算二个物体接触界面的热传导和热辐射,给定材料内部热生成(给定热源);
热分析采用隐式求解方法,过程控制有:
(1) 稳态分析还是瞬态分析;
(2) 线性问题还是非线性问题;
(3) 时间积分法:Crank-Nicholson 法(a=0.5)和向后差分法( a=1);
(4) 求解器:直接法或迭代法;
自动时步长控制。
LS-DYNA 不可压缩流求解器是 970 版新增加的功能,用于模拟分析瞬态、不可压、粘性流体动力学现象。求解器中采用了超级计算机的算法结构,在确保有限元算法优点的同时计算性能得到大幅度提高,从而在广泛的流体力学领域具有很强的适用性。
(1) 多种控制选项和用户子程序使得用户在定义和分析问题时有很大的灵活性。
(2) 输入文件可分成多个子文件;
(3) 用户自定义子程序;
(4) 二维问题可以人工控制交互式或自动重分网格(REZONE);
(5) 重启动;
(6) 数据库输出控制;
(7) 交互式实时图形显示;
(8) 开关控制-可监视计算过程的状态;
(9) 对 32 位计算机可进行双精度分析。
LS-DYNA 利用 ANSYS、LS-INGRID、ETA/FEMB、TrueGrid、LS-POST 和 LS-PREPOST 强大的前后处理模块,具有多种自动网格划分选择,并可与大多数的 CAD/CAE 软件集成并有接口。
后处理:结果的彩色等值线显示、梯度显示、矢量显示、等值面、粒子流迹显示、立体切片、透明及半透明显示;变形显示及各种动画显示;图形的 PS、TIFF 及 HPGL 格式输出与转换等。
LS-DYNA 970 版的 SMP 版本和 MPP 版本是同时发行的。MPP 版本使一项任务可同时在多台分布计算机上进行计算,从而最大限度地利用已有计算设备,大幅度减少计算时间。计算效率随计算机数目增多而显著提高。
LS-DYNA 970 版的 SMP 版本和 MPP 版本可以在 PC 机(NT、LINUX 环境)、UNIX 工作站、超级计算机上运行。
1.1.2 LS-DYNA应用领域
1 汽车工业
(1) 碰撞分析
(2) 气囊设计
(3) 乘客被动安全
(4) 部件加工
2 航空航天
(1) 鸟撞
(2) 叶片包容
(3) 飞机结构冲击动力分析
(4) 碰撞,坠毁
(5) 冲击爆炸及动态载荷
(6) 火箭级间分离模拟分析
(7) 宇宙垃圾碰撞
(8) 特种复合材料设计
3 制造业
(1) 冲压
(2) 锻造
(3) 铸造
(4) 切割,等
4 建筑业
(1) 地震安全
(2) 混凝土结构
(3) 爆破拆除
(4) 公路桥梁设计
5 国防
(1) 内弹道和终点弹道;
(2) 装甲和反装甲系统;
(3) 穿甲弹与破甲弹设计;
(4) 战斗部结构设计;
(5) 冲击波传播;
(6) 侵彻与开坑;
(7) 空气,水与土壤中爆炸;
(8) 核废料容器设计,等
6 电子领域
(1) 跌落分析
(2) 包装设计
(3) 热分析
(4) 电子封装
7 石油工业
(1) 液体晃动
(2) 完井射孔
(3) 管道设计
(4) 爆炸切割
(5) 事故模拟
(6) 海上平台设计
8 其它应用
(1) 玻璃成型
(2) 生物医学
(3) 体育器材(高尔夫杆,高尔夫球,棒球杆,网球拍、头盔)。