第11章 LS-DYNA自适应网格划分
LS-DYNA使用指南中文版本
第一章引言ANSYS/LS-DYNA将显式有限元程序LS-DYNA和ANSYS^序强大的前后处理结合起来。
用LS-DYNA的显式算法能快速求解瞬时大变形动力学、大变形和多重非线性准静态问题以及复杂的接触碰撞问题。
使用本程序，可以用ANSYS建立模型，用LS-DYNA做显式求解，然后用标准的ANSYSH处理来观看结果。
也可以在ANSYS和ANSYS-LS-DYNA之间传递几何信息和结果信息以执行连续的隐式-显式/显式-隐式分析，如坠落实验、回弹、及其它需要此类分析的应用。
显式动态分析求解步骤概述显式动态分析求解过程与ANSYS程序中其他分析过程类似，主要由三个步骤组成：1:建立模型（用PREP7前处理器）2:加载并求解（用SOLUTIONS理器）3:查看结果（用POST1和POST26后处理器）本手册主要讲述了ANSYS/LS-DYN湿式动态分析过程的独特过程和概念。
没有详细论述上面的三个步骤。
如果熟悉ANSYS程序，已经知道怎样执行这些步骤，那么本手册将提供执行显式动态分析所需的其他信息。
如果从未用过ANSYS就需通过以下两本手册了解基本的分析求解过程：•ANSYS Basic Ana lysis Guide•ANSYS Modeli ng and Meshi ng Guide使用ANSYS/LS-DYNA寸，我们建议用户使用程序提供的缺省设置。
多数情况下，这些设置适合于所要求解的问题。
显式动态分析采用的命令在显式动态分析中，可以使用与其它ANSYS分析相同的命令来建立模型、执行求解。
同样，也可以采用ANSYS图形用户界面（GUI）中类似的选项来建模和求解。
然而，在显式动态分析中有一些独特的命令，如下：EDADAPT:激活自适应网格EDASMP:创建部件集合EDBOUND定义一个滑移或循环对称界面EDBVIS :指定体积粘性系数EDBX :创建接触定义中使用的箱形体EDCADAPT：指定自适应网格控制EDCGEN指定接触参数EDCLIST :列出接触实体定义EDCMORE为给定的接触指定附加接触参数EDCNSTR定义各种约束EDCONTACT指定接触面控制EDCPU : 指定CPU时间限制EDCRB :合并两个刚体EDCSC :定义是否使用子循环EDCTS :定义质量缩放因子EDCURVE定义数据曲线EDDAMP定义系统阻尼EDDC :删除或杀死/重激活接触实体定义EDDRELAX进行有预载荷几何模型的初始化或显式分析的动力松弛EDDUMP指定重启动文件的输出频率（d3dump）EDENERGY定义能耗控制EDFPLOT:指定载荷标记绘图EDHGLS:定义沙漏系数EDHIST :定义时间历程输出EDHTIME :定义时间历程输出间隔EDINT :定义输出积分点的数目EDIS :定义完全重启动分析的应力初始化EDIPART :定义刚体惯性EDLCS :定义局部坐标系EDLOAD:定义载荷EDMP:定义材料特性EDNB :定义无反射边界EDNDTSD清除噪声数据提供数据的图形化表示EDNROT:应用旋转坐标节点约束EDOPT:定义输出类型，ANSYS或LS-DYNA EDOUT:定义LS-DYNA ASCII输出文件EDPART :创建，更新，列出部件EDPC :选择、显示接触实体EDPL :绘制时间载荷曲线EDPVEL:在部件或部件集合上施加初始速度EDRC :指定刚体/变形体转换开关控制EDRD:刚体和变形体之间的相互转换EDREAD:把LS-DYNA的ASCII输出文件读入到P0ST26勺变量中EDRI :为变形体转换成刚体时产生的刚体定义惯性特性EDRST:定义输出RST文件的时间间隔EDSHELL :定义壳单元的计算控制EDSOLV:把“显式动态分析”作为下一个状态主题EDSP :定义接触实体的小穿透检查EDSTART:定义分析状态（新分析或是重启动分析）EDTERM:定义中断标准EDTP :按照时间步长大小绘制单元EDVEL :给节点或节点组元施加初始速度EDWELD定义无质量焊点或一般焊点EDWRITE:将显式动态输入写成LS-DYNA输入文件PARTSEL:选择部件集合RIMPORT :把一个显式分析得到的初始应力输入到ANSYSREXPORT:把一个隐式分析得到的位移输出到ANSYS/LS-DYNAUPGEOM相加以前分析得到的位移，更新几何模型为变形构型关于ANSYS命令按字母顺序排列的详细资料（包括每条命令的特定路径），请参阅《ANSYS Comma nds Refere nee》。
第11章 LS-DYNA自适应网格划分汇编
网格自适应 h-Adaptive
精确: 使用two-pass 方法
One pass 可选
LS-DYNA 网格自适应 — 网格细化 四边形单元细化 – 一级细化
三角形单元细化 — 一级细化
LS-DYNA 网格自适应 — 网格细化
网格细化中自适应级别由三条规则限制:
• 级别数由允许的最大自适应级别限制，一般设置为 3 或 4 • 网格自适应级别必须与邻近单元的自适应级别最多差别一个级别 • 总的单元树受到内存的限制，一旦内存容量不足，自适应终止
*CONTROL_ADAPTI总的角度相对于周围的单元改变: adpopt = 2 . 错误容限 = 0.1 度
. 细化级别: maxlvl = 2
*PART . blank (part 3) - h-adaptivity (adpopt = 1)
自适应网格粗化 在成型过程中小单元是必要的 • 在90°范围内至少4个单元才能有最好的回弹结果 • 在最佳的位置自适应网格产生单元 • 问题：在深冲压中侧壁的单元过多
自适应网格粗化
一新功能在发展，即在回弹前对自适应网格进行粗化 好处
• 速度
• 内存 缺点
• 丢失信息（精度？）
过程 • dynain文件
粗化能节省 CPU 时间，并无明显的精度损失
自适应网格粗化 精度考虑 统一的网格与自适应 + 粗化
LS-DYNA 关键字 *CONTROL_ADAPTIVE • 自适应时间间隔
• 类型
• 错误容限 • 细化级别
*CONTROL_ADAPSTEP
• 在自适应过程中控制接触界面力 *PART
• 激活自适应 parts
网格自适应
一边两分规则
LS-DYNA 理论及功能(简介)
LS-DYNA 理论及功能LS-DYNA 的理论及功能LS-DYNA 发展概况 (LS-DYNA Introduction)LS-DYNA是以显式为主、隐式为辅的通用非线性动力分析有限元程序，特别适合求解 各种二维、三维非线性结构的高速碰撞、爆炸和金属成形等非线性动力冲击问题，同时可以 求解传热、流体及流固耦合问题。
DYNA 程序系列最初是 1976 年在美国 Lawrence Livermore National Lab. 由 J.O.Hallquist 博士主持开发完成的，主要目的是为武器设计提供分析工具，后经 1979、1981、1982、1986、 1987、1988 年版的功能扩充和改进，成为国际著名的非线性动力分析软件，在武器结构设 计、内弹道和终点弹道、军用材料研制等方面得到了广泛的应用。
1988 年 J.O.Hallquist 创建 LSTC 公司，推出 LS-DYNA 程序系列，主要包括显式 LS-DYNA2D、LS-DYNA3D、隐式 LS-NIKE2D、LS-NIKE3D、热分析 LS-TOPAZ2D、 LS-TOPAZ3D、前后处理 LS-MAZE、LS-ORION、LS-INGRID、LS-TAURUS 等商用程序， 进一步规范和完善 DYNA 的研究成果，陆续推出 930 版（1993 年）、936 版（1994 年）、940 版（1997 年），950 版（1998 年）增加了汽车安全性分析（汽车碰撞、气囊、安全带、假人）、 薄板冲压成形过程模拟以及流体与固体耦合（ALE 和 Euler 算法）等新功能，使得 LS-DYNA 程序系统在国防和民用领域的应用范围进一步扩大，并建立了完备的质量保证体系。
1997 年LSTC公司将LS-DYNA2D、LS-DYNA3D、LS-TOPAZ2D、LS-TOPAZ3D等程序 合成一个软件包，称为LS-DYNA，PC版的前后处理采用ETA公司的FEMB，新开发的后处 理器为LS-POST。
ANSYSLS-DYNA网格划分
ANSYSLS-DYNA网格划分ANSYS LS-DYNA结构冲击动力学分析专题培训学习心得——网格控制心得：本次培训最大的收获在于利用workbench进行模型的前处理方面，尤其是网格划分控制上，前期我们进行分析的主要网格划分方式多为系统自动划分，对于结构形式复杂的模型，很多时候都不太会对网格进行控制。
在三维网格划分方面主要有以下几个方式1. 四面体网格划分2. 扫略网格划分3. 多区网格划分4. 六面体为主网格划分5. 自动网格划分算法区别这些我们都有了解，而网格划分算法中的Patch Conforming 算法和Patch Independent 算法的区别却不太清楚，其主要差别在于Patch Independent 算法有较强的几何容忍度，小于某一给定尺寸的几何形状会被忽略，但同时也带来了计算精度有偏差的问题，如何均衡这两者的关系需要根据工程实例情况来进行区分。
对于单个模型的多次网格控制对于一个模型可以添加多个不冲突的网格控制，尺寸控制几乎可以跟任何一种划分方式合用，这样可以保证模型网格的规则性，也可以控制单元数量（因为在DYNA中，单元数量及最小单元尺寸与计算效率及计算精度关系非常密切）。
局部网格控制局部网格控制同样也是以尺寸控制为基础，用单元尺寸、线份数、影响球等手段来实现所关注的局部网格质量。
这个控制的合理应用可以提高计算的效率和精度。
在高级尺寸共功能上，打开调整曲率功能可以调整曲率法向角，细化转角处网格，还可以通过控制狭缝间的网格层数来对细微部分进行细化处理。
单元数量控制DANY的分析对单元数量很敏感，普通的双核CPU的计算机，计算400万单元的模型已经是极限，大型计算机的并行计算也需要进行合理的网格数量及尺寸的控制，另外不同的分析类型对网格质量的要求也不同。
运用DYNA进行碰撞模拟时多采用的是显示动力学分析的模式，这是因为一方面，计算收敛程度高，另一方面，计算结果更为精确，只是由于显式分析计算量大，对网格的质量要求就很高。
LS-DYNA使用指南
第七章材料模型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来定义材料模型。
LS-DYNA自适应网格划分
ANSYS/LS-DYNA自适应网格划分在金属成形和高速撞击分析中，物体要经历很大的塑性变形。
单积分点显式单元，常用于大变形，但是在这种情况下，由于单元纵横比不合适可能给出不精确的结果。
为了解决这一问题，ANSYS/LS-DYNA程序可以在分析过程中自动重新划分表面来改善求解精度。
这一功能，即自适应网格划分，由EDADAPT 和EDCADAPT 命令控制。
EDADAPT 命令在一个指定的PART内激活自适应网格划分。
（用EDPART 命令创建或显示有效PART IDs），例如，为了给PART1打开自适应网格划分，可以执行下列命令：EDADAPT，1，ON注意—自适应网格划分功能仅对包含SHELL163单元的部件有效。
当此项功能打开时，分析中该部件的网格将自动重新生成。
从而保证在整个变形过程中有合适的单元纵横比。
自适应网格划分一般应用在大变形分析例如金属变形中（调节网格最典型的应用是板料）。
在一个模型中要在多个部件上应用此功能，必须对每个不同的PART ID执行EDADAPT 命令。
缺省时，该功能是关闭的。
在指定哪些部件重新划分后，必须用EDCADAPT 命令定义网格划分参数。
采用EDADAPT 命令定义需要网格划分的所有PART ID号，用EDCADAPT 命令对其设置控制选项。
EDCADAPT 命令控制的参数如下所示：·Frequency(FREQ)-调节自适应网格划分的时间间隔。
例如，假设FREQ设置为0.01，如果单元变形超过指定的角度容差，则其将每隔0.01秒被重新划分一次（假设时间单位为秒）。
因为FREQ的缺省值为0.0，所以在分析中应用自适应网格划分时必须指定此项。
·Angle Tolerance(TOL)-对于自适应网格划分（缺省值为1e31）有一个自适应角度公差。
TOL域控制着单元间的纵横比，它对保证结果的精度是非常重要的，如果单元之间的相对角度超过了指定的TOL值，单元将会被重新划分。
LS-DYNA(ANSYS)中文培训教程_2
就当前所选的节点使用 CM 命令创建节点组元 Page 6
加载，定义边界条件及刚体
...一般加载过程
最好给创建的集合定义一个意义明确的名字，这样可以避免将载荷加 到错误的结构上。
定义数组参数 :
在 ANSYS/LS-DYNA中，所有载荷都是按特定时间间隔施加 ，并成 组地定义时间数组参数及相应的载荷数组参数 (*DIM and *SET):
最后，你要指定分析阶段，默认值是瞬态阶段 “Transient only” ，如 果你的载荷是应用于动力松驰分析（用来模拟隐式静态求解），选 “Dynamic relax”。如果包括两种分析，请选 “Trans and Dynam” 。
Page 11
加载，定义边界条件及刚体
...一般加载过程
两种类型均需考虑(继续): 实际温度载荷=用户施加的温度-TREF
TREF 默认值为0
PLANE162, SHELL163, and SOLID164 单元支持
PLANE162 不支持 “drelax” 文件(参见第4-3章), 只有瞬态热载荷或 冲击热载荷适于单元 PLANE162 设置ALPX=0 ，仅定义与温度相关的机械力学性能(EX, PRXY, …) ，使 瞬态热－显式分析稳定。
使用EDLOAD命令加载:
完成组元（component）或Part和参数数组定义后可以使用 EDLOAD 命令来加载:
Solution > Loading Options > Specify Loads
先选择载荷选项中的（Add Loads）. 这一选项同时也可以列出载荷 （listing Loads） 和删除载荷（ deleting load）.
Page 5
自适应网格生成算法与应用
自适应网格生成算法与应用自适应网格生成算法是一种重要的数值计算方法，广泛应用于科学计算、工程模拟和物理仿真等领域。
本文将详细介绍自适应网格生成算法的原理和应用，并探讨其在不同领域中的实际应用案例。
一、自适应网格生成算法的原理自适应网格生成算法是一种基于网格重构的数值计算方法。
其核心思想是根据问题的特性和求解需求，在计算过程中动态调整网格的大小和形状，以提高数值计算的效率和精度。
自适应网格生成算法通常包括以下几个步骤：1. 初始网格生成：根据问题的几何形状和边界条件，生成初始的网格。
2. 误差估计：通过计算网格单元内部的数值误差或局部残差，评估当前网格的精度。
3. 网格划分：根据误差估计结果，确定需要细分的网格单元，并对其进行细分操作。
4. 网格合并：根据误差估计结果，确定需要合并的网格单元，并对其进行合并操作。
5. 网格重构：根据细分和合并操作的结果，对整个网格进行重构，以适应新的问题求解需求。
6. 计算求解：在重构后的网格上进行数值计算，并更新问题的解。
7. 收敛判断：通过对比前后两次计算结果，判断数值计算是否收敛，若未收敛，则返回第二步。
二、自适应网格生成算法的应用自适应网格生成算法在科学计算、工程模拟和物理仿真等领域得到了广泛应用，能够提高计算效率和精度，降低计算成本。
下面将分别从这几个领域进行具体介绍。
1. 科学计算：自适应网格生成算法在科学计算中的应用非常广泛。
例如，对于流体力学中的空气动力学问题，通过自适应网格生成算法，可以在物体表面和流动区域自动调整网格的密度，以捕捉流动细节，提高计算精度和效率。
2. 工程模拟：工程模拟是一个复杂的过程，需要对多个物理场耦合求解。
自适应网格生成算法可以根据不同物理场的特性，分别调整网格在不同区域的分布和密度，以适应多个物理场的求解需求，从而提高工程模拟的准确性和可信度。
3. 物理仿真：自适应网格生成算法在物理仿真中也有重要应用。
例如，在计算结构强度和变形时，通过自适应网格生成算法，可以根据不同部位的应力和应变分布，自动调整网格的大小，以更好地捕捉局部细节，提高仿真结果的准确性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
过程 dynain文件
q (度)
自适应网格粗化 例子: NUMISHEET ‘96 S-rail, aluminum, 10kN BHF
自适应网格粗化 粗化参数研究 Numisheet ‘96 S-Rail – 三级自适应 初始网格: 550 单元
例子— 深冲压
结束语
• 自适应细化 • 自适应粗化 •制造过程的生产模式 •在许多另外的领域有越来越多的应用（包括碰撞）
Pipe-on-Pipe 自适应
例子— 深冲压
接触 0.01 20 • *CONTACT_FORMING_ONE_WAY_SURFACE_TO_SURFACE blank 与 punch间接触 在 blank 与binder间接触 在 blank 与die间接触
parts
例子— 深冲压 punch (1) binder (2) blank (3) die (4)
例子— 深冲压
die (4)
• *MAT_RIGAD_RIGID_BODY punch (1) *MAT_RIGID *BOUNDARY_PRESCRIBED_MOTION_RIGID
比较：角度改变对逼近接触
自适应网格粗化 在成型过程中小单元是必要的 90°范围内至少4个单元才能有最好的回弹结果
自适应网格粗化
一新功能在发展，即在回弹前对自适应网格进行粗化 好处
• 速度
网格自适应 h-Adaptive
精确: 使用two-pass 方法
One pass 可选
LS-DYNA 网格自适应 — 网格细化 四边形单元细化 – 一级细化
三角形单元细化 — 一级细化
LS-DYNA 网格自适应 — 网格细化
网格细化中自适应级别由三条规则限制:
3或4
第十一章
LS-DYNA中的自适应网格划分
在非线性结构分析中两种主要具有冲击影响的趋势
自适应重划分 • 自适应重划分提供一种方法，用来精确的处理那些不易预测 的结构失效模式：
— 屈曲
— — Localization
分布式并行处理
—
例子— 深冲压
*CONTROL_ADAPTIVE
: adpfreq = 5
– 总的角度相对于周围的单元改变: adpopt = 2 = 0.1 度
: maxlvl = 2
*PART blank (part 3) - h-adaptivity (adpopt = 1)
网格自适应
一边两分规则
网格自适应 细化判定 • 基于角度
—
平面内 平面外
无变形
变形 — —
无变形
变形
定义板料的网格以适应模具是非常困难的 ，网格自适应使板料的网格非常简单化
基于逼近接触的自适应
—
自适应 优点
• 应用领域的的扩展
基本原理 LS-DYNA 自动修改网格 • 错误标准 — L Y — — — Impending contact
多工序冲压成型仿真
自适应冲压 网格自动重划分以便捕捉到模具的细节
自适应的类型
r-method,重新部署节点
— — — LS-DYNA 二维单元 — h-method, 调整单元的尺寸 h — LS-DYNA 壳单元 — alpha 测试中 p-method, 调整单元插值阶数 p — 在显式动力学中一般不使用
粗化能节省 CPU 时间，并无明显的精度损失
自适应网格粗化 精度考虑 统一的网格与自适应 + 粗化
LS-DYNA 关键字 *CONTROL_ADAPTIVE
细化级别
*CONTROL_ADAPSTEP
L *PART