ls dyna

LS-DYNA爆炸模拟资料及技术问答汇总专帖

[i=s] 本帖最后由 wyc412721 于 2010-10-30 22:38 编辑 [/i]

[b][size=4][/size][/b]

{.@

L$s

O

:r t!o s

q-c9Z:z

f,D

[b][size=4][/size][/b]

[b][size=4][color=red]注：本帖属于资料及技术问题讨论帖，感谢辛苦之类的话就不要留言了。[/color][/size][/b]

-?%Q6X

t'N

[b][size=4][/size][/b] [b][size=4][/size][/b][b][size=4]

%{8X

X,^+K ~

6J

^

q.o

E+V

[b][size=4]1、爆炸模拟的一般过程[/size][/b]

[/size][/b]（1）首先要有一个总体思路，怎样[font=宋体]把爆炸程转化为程序模型（明白该怎样建立自己的模型，选取材料模型，采用哪儿种算法计算求解），这一步很重要，直接关系计算结果的合理性。[/font] （2）前处理，包括：画图、建模、划分网格、定义边界及约束条件、设定分析控制选项、输出K文件等。十分重要，需要注意）。

[b][color=red][/co

1 如何处理LS-DYNA中的退化单元？

在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1，N2，N3，N4，N4，N4，N4，N4和N1，N2，N3，N4，N5，N5，N6，N6。这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。所以为网格划分的方便起见，我们还是在同一Part ID下划分网格，通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。

2 LS-DYNA中对于单元过度翘曲的情况有何处理方法 有两种方法：

1. 采用默认B-T算法，同时利用*control_shell控制字设置参数BWC＝1，激活翘曲刚度选项；

2. 采用含有翘曲刚度控制的单元算法，第10号算法。该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。

除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。

3 在ANSYS计

！！！Dyna批处理！！！

我都是自己编辑批处理来进行计算的，比ansys的好，我自己觉得，虽然不算高科技，总有人不知道，贴出来大家看看吧，具体办法如下：

1、设置环境变量。右键单击‘我的电脑’，选属性----高级-----环境变量，找到path变量，在变量值（很长很长）的最后添加一个分号‘;’，然后查找你的ansys安装目录下一个叫ls970的文件，找到之后把路径复制下来，拷贝到刚才分号的后面，确定再确定，设置环境变量的目的是使系统认识ls970命令在什么地方；环境变量设置成功与否，可这样检验，在运行里面键入ls970，回车，有窗口调出就行了。1 y Y$ Z# }/ c2 \\7 P! b 2、编辑批处理文件。在你ansys工作目录（就是你要进行dyna计算的目录）新建文本文件，将后缀名从txt改变为bat，中间会询问你是否改变后缀名，是！！！右键单击之，选择编辑，出现空白的编辑界面：9 H! {; \\, [0 ?$ V& R

情况A：你只想进行某一个单独的计算，只需要在批处理文件中键入; g$ {( v- ~\ ls970 i=jobname.k

保存，然后令批处理文件与k文件同目录，双击批处理，计算开始！！！

情况B：想进行系列计

ANSYS/LS-DYNA前处理操作步骤

在ANSYS Launcher界面中，选择ANSYS Mechanical/LS-DYNA 1、菜单过滤 Main Menu→Preprocessor→LD-DYNA Explicit→OK 2、设置文件名及分析标题 Utility Menu→File→change Jobname→cutting→New log and error file ：YES→OK Utility Menu→File→change Title→cutting analysis →OK

3．定义单元类型

Main Menu| Preprocessor| Element Type| Add/Edit/Delete

选择3D SOLID164 →OK→options→选择const.stress；Lagrangian→OK 4．输入材料性能

Main Menu| Preprocessor| | Material Props| Material Models 在窗口中选择Rigid Material定义刀片材料 在DENS中 输入5.2e3 EX中 输入4.1e11

第二部分 ANSYS/LS-

DYNA程序的使用方法

1 概述

ANSYS/LS-DYNA 程序系统是将非线性动力分析程序LS-DYNA显式积分部分与

ANSYS程序的前处理PREP7和后处理POST1、POST26连接成一体。这样既能充分运用LS-DYNA程序强大的非线性动力分析功能，又能很好地利用ANSYS程序完善的前后处理功能来建立有限元模型与观察计算结果，它们之间的关系如下。

ANSYS前处理PREP7 Jobname.DB Jobname.k ANSYS/LS-DYNA Jobname.RST d3plot Jobname.HIS d3thdt

初始穿透问题

Warning Initialization - node will not be moved Penetrating node is too far below surface

single surface interface # = 1 penetrating node # = 990

segment # = 140

penetration = 0.1229370E+00

s,t = -0.9817361E+00 0.9780035E+00

方法一：在建立模型时就要注意该问题，尽量保持接触对中的接触空隙（如在空间位置中考虑壳单元的厚度），但对于复杂的模型，不可避免会出现初始穿透，此时可根据第

一次递交后程序给出的穿透信息（在MESSAGE 和D3HSP 文件中有详细记录），按照提示移动相关节点，消除穿透。但该方法工作量大，且会在局部改变模型的形状。9 v. T1

z6 ]! D1 i

方法二：对于LS-DYNA接触问题，有一接

LS-DYNA常见问题汇总2.0

资料整理：yuminhust2005

职 务：A2:LS-DYNA 版主

资料来源：网络和自己的总结

更新时间：2008-10-12

LS-DYNA常见问题汇总2.0 yuminhust2005 2008-10-12 目录

1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. 22. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. 33. 34. 35. 36. 37. 38. 39. 40. 41. 42. 43. 44. 45.

LS-DYNA 简介 ...................................................................................................................................................... 1 单位制度 ................

1 如何处理LS-DYNA中的退化单元？

在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1，N2，N3，N4，N4，N4，N4，N4和N1，N2，N3，N4，N5，N5，N6，N6。这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。所以为网格划分的方便起见，我们还是在同一Part ID下划分网格，通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。

2 LS-DYNA中对于单元过度翘曲的情况有何处理方法 有两种方法：

1. 采用默认B-T算法，同时利用*control_shell控制字设置参数BWC＝1，激活翘曲刚度选项；

2. 采用含有翘曲刚度控制的单元算法，第10号算法。该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。

除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。

3 在ANSYS计

LS-DYNA_关键字文件

LS-DYNA关键字文件

LS-DYNA_关键字文件

内容 板料 (Blank) 凸模 (Punch) 凹模 (Die) 压边圈 (Binder) 控制参数 FEM模型

LS-DYNA_关键字文件

关键字文件 Dynaform生成的输入文件，由LS-DYNA读取并运行 文本格式，可用文本编辑器打开和修改 文件后缀：*.dyn, *.mod– *.dyn: 定义板料成形的所有工艺参数 – *.mod: 定义模型的节点、单元等有限元模型信息

LS-DYNA_关键字文件

关键字文件 以$开始的是注释行$ ETA/DYNAFORM : DYNA3D(960) INPUT DECK$ DATE $ : Nov 3, 2003 at 17:19:47

输入文件的第一个关键词都是*KEYWORD*KEYWORD

该分析项目的标题*TITLENumisheet 96: S-rail Stamping

LS-DYNA_关键字文件

关键字文件基本结构

LS-DYNA_关键字文件

板料 (Blank)$---+----1----+----2----+----3----+----4----+----5----+----6----

1 如何处理LS-DYNA中的退化单元？

在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1，N2，N3，N4，N4，N4，N4，N4和N1，N2，N3，N4，N5，N5，N6，N6。这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。所以为网格划分的方便起见，我们还是在同一Part ID下划分网格，通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。

2 LS-DYNA中对于单元过度翘曲的情况有何处理方法 有两种方法：

1. 采用默认B-T算法，同时利用*control_shell控制字设置参数BWC＝1，激活翘曲刚度选项；

2. 采用含有翘曲刚度控制的单元算法，第10号算法。该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。

除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。

3 在ANSYS计