LS-DYNA常见问题汇总2.0 yuminhust2005 2008-10-12 17. 怎样进行二次开发

1.二次开发的过程和需要的文件

在回答这个问题前，给大家一点背景知识（针对vwindows平台）：

1、DYNA只是一个求解器（一个固定的可执行的exe文件），它可以读取K文件，并进行计算，最终输

出结果。

2、DYNA开发过程：

编译连接

DYNA源代码－－－－>DYNA.OBJ文件或DYNA.lib————>DYNA.exe

如果进行二次开发，实际上就是要修改原有的exe文件，在其中加入自己的程序（子材料什么的）。而不管你是―一般用户‖还是正版用户，都不太可能直接修改exe文件（除非你特别特别牛，能直接编写机器语言（比汇编语言还低一级别），开个玩笑，基本不可能实现）。所以目前DYNA开发商，给正版用户提供了下面的方式进行二次开发：

编译连接

自定义程序代码（F文件）————>OBJ文件 + DYNA.lib————>新的DYNA.exe文件（其中包含自定义的部分）

从上面可以看出二次开发实际上是重新生成了一个全新的可执行程序，这一过程需要F文件、 DYNA.lib和其它资源文件，而不需要原来的DYNA.EXE文件，所以不管你是正版还是别的，如果没有DYNA.lib文件等都是不可能进行二次开发的。而―一般用户 ‖和正版用户的区别在于：

1、正版用户可以直接要求开发商提供F文件、 DYNA.lib和其它资源文件，而一般用户只有通过其它渠道。

2、将K文件提交DYNA.EXE计算时，正版用户由于有license ，可以直接通过检查并进行计算；而―一般用户‖目前基本上是通过修改 DYNA.exe文件。跳过检查，最终实现计算的。

一般用户‖要实现二次开发就必需首先有F文件、 DYNA.lib和其它资源文件；通过编译连接后生成新的EXE文件；再修改EXE文件，跳过license检查（这个网上有说明具体过程，所以关键是要得到那几个二次开发文件）。

―一般用户‖完全有可能进行二次开发，关键是要通过其它渠道得到正版用户才有的F文件、DYNA.lib和其它资源文件 2.子程序的编写

见白金泽和赵海鸥的书。

第 30 页

LS-DYNA常见问题汇总2.0 yuminhust2005 2008-10-12 18. 如何实现程序的连续计算

通过批处理设置连续计算可以省下很多时间。首先创建一个记事本，将扩展名改成.bat比如run.bat 。编辑批处理文件，写下面的语句：

f:

cd F:\\111\\

d:\\lsdyna\\program\\ls971_s_7600_win32.exei=*.kmemory=内存大小 cd..

第一步： f：这句很重要（假如你的k文件目录在f盘），因为批处理就是在命令提示符（DOS）下运行的。在某个盘（除了c盘）建立一个目录，直接用cd 命令是进入不到这个相应的目录的。比如在f盘建立一个文件夹，命名为111。

用cd f：\\111这句命令是进不到这个目录的，还是在c：。

必须先输入f：

在输入 cd 111才能进入目录。

第二步：找到求解器所在的目录和命名，再设置好memory大小（也可以在k文件第一行*keyword 后面设置），CPU个数即可。

D:\\lsdyna\\program\\ls971_s_7600_win32.exei=*.kmemory=内存大小

C:\\Program Files\\LSDYNA\\program\\ls970_s_6763.exei=*.kmemory=内存大小 上面是ls971求解器的，下面是ansys自带的970求解器的。

第三步：复制第二步的内容，有几个k文件就复制几次，改掉相应的路径和k文件名即可。若想求解完毕后自动关机，还可以在最后一行写上 “shutdown -s\。

第四步：保存批处理文件，双击求解。

最后在提醒一点，每个k文件最好都提前试一下，运行无误后写入让它自动求解。免得设置好高高兴兴走了，回来一看有问题，那就白高兴了，哈哈！

祝大家学业有成，工作顺利！ Timothy1985

第 31 页

LS-DYNA常见问题汇总2.0 yuminhust2005 2008-10-12 19. 如何得知LS-DYNA为一套完整的软件？

解答：LS-DYNA一直以来都是独立软件，原因很简单：

(1)美国原厂Livermore Software Technology Corp.(简称LSTC)授权「势流科技」为台湾地区LS-DYNA软件的独家代理商。

(2)LS-DYNA拥有LSTC官方的前、后处理程序：FEMB及LS-Pre/Post，无需另外增购。 LS-DYNA强调「One Model、One Code、Multi-results」，单一核心程序即可进行Implicit、Explicit及混合Implicit/Explicit的求解方式，可以为您解决静力、准静力、掉落、冲击、模态分析、振动、预应力/预应变、流固耦合、流场、热应力....等等的力学分析，使用者不必再额外添购模块程序，徒增预算的浪费与软件学习的负担。

20. LS-DYNA与市面上其它的前处理软件兼容性如何？

解答：由于LS-DYNA是全球使用率最高的结构分析软件，因此其它的前处理软件与LS-DYNA是完全兼容的。在此要强调的是：LS-DYNA的官方前处理程序为FEMB，因为FEMB是专门为LS-DYNA量身订作的前处理程序，有许多设定条件及架构逻辑是其它前处理软件所难望其项背的，为了避免在学习LS-DYNA的过程及操作上产生困扰，强烈建议使用者采用原厂出品的FEMB来做为LS-DYNA的前处理工具，使用者必定更能体会LS-DYNA直觉式的设定与强大的分析能力。

21. LS-DYNA似乎很重视「Contact Algorithm」，这是为什幺？

解答：是的，LS-DYNA很早以前就已经发展「接触算法」，这是因为基础力学所分析的对像均只考虑「力的受体」，故输入条件皆为外力量值。然而在真实情况下，物体受力通常是因为与其它的物体发生「接触」（Contact）才受力，此时外力量值是无法预期的，应该输入的条件往往都是几何上的接触条件。因为有完备的接触力学演算方式，LS-DYNA才得以忠实的仿真现实环境的复杂结构行为。

22. 如果要利用LS-DYNA进行MPP(平行运算)的计算，硬件配备及操作系统有无特殊需求？

解答：不论是PC cluster、工作站及一般的PC环境，都适合执行LS-DYNA的MPP平行运算功能，一般我们还是会建议要用来执行平行运算的计算机群组，彼此的等级宜尽量一致；操作系统方面并无特别需求，以一般的windows 2000、LINUX或是UNIX皆可执行。

国外已有很多厂商利用非办公时间，将办公室内的计算机串连在一起，配合LS-DYNA来分析问题，宛如一部超级计算机，不仅可以有效提升研发的竞争力，同时亦可弹性地运用计算机资源，一举数得。

23. FEMB能够自动产生有限元素网格吗？

解答：FEMB当然可以自动产生有限元素网格，使用者再也不必费心在每个几何边界上指定结点数量，仅需要输入元素尺寸的参考值，FEMB便会依此产生网格。当然，如果能够先在CAD进行合理的简化，auto-mesh得到的元素品质会更好。

在国外CAE领域有两句名言：「Well Done is Quickly Done」、「Quality Mesh leads to Quality Analysis」，因此势流科技建议使用者应该多花点心思在前处理上，这对后续的分析工作有着莫大的帮助。

24. 实际的产品CAD图文件，有许多复杂的几何造型，这些feature是否都该纳入分析考量？

解答：这样的问题事实上是属于「Engineer Judgement」的层次，而非软件层次的问题。 使用者当然可以将完全按照CAD几何来建立有限元素的模型，然而任何一位有经验的工程师一定会先将不必要的几何造型予以忽略，以提升分析效益，并将分析课题明确化。

以Block为例，有些工程师会以Solid element来建构完整的模型，有些则采用Shell element，甚至有些人会采用1D element...，这些方式都无伤大雅，重点是工程师如何从分析结果获得充份的信息。

第 32 页

LS-DYNA常见问题汇总2.0 yuminhust2005 2008-10-12 势流科技除了提供软件分析技术之外，同时更致力于协助业界将CAE导入于生产流程的正确位阶，培养正确的CAE分析观念，以提升产品研发品质，缩短生产流程。

25. LS-DYNA是否只能进行壳元素的分析？

解答：LS-DYNA不只能分析实体元素(Tetra、Hexa)、薄/厚壳元素(Quad.、Tri)，同时还有梁元素、质点、Spring/Damper、spotweld等元素型态，另外还有SPH element。

之所以建议台湾电子产业尽量以壳元素来仿真壳件产品，是由于电子产品的组件当中，板壳类占了90%以上，基于结构行为、准确性及分析效能的考量。除非使用者有特殊的分析需求，否则并不建议采用实体元素来仿真板壳类型的产品组件。

我们藉由右边的简单案例，以了解不同元素型态在仿真壳件产品的结果表现：

(1)实体元素在板壳厚度方向上至少需要分割三层以上，其位移变形量才能收敛至一定值；而壳元素及厚壳元素仅需一层，在位移量上即可达到理想精度。

(2)若是发生降伏产生大变形，在厚度方向上即使分割8层以上的实体元素，其掌握材料塑性应变的收敛速度仍然相当缓慢，而付出的时间代价为shell element的25倍以上。

厚壳元素的运算速度及精度与三层实体元素相当，然而就塑性应变的精度而言，厚壳元素的表现比实体元素要好太多。

26. 在LS-DYNA里该如何强制刚体运动？

解答：CAE分析常需要藉由刚体运动，以简化复杂物理行为。应用技巧有二：给定初始速度或是强制运动轨迹(prescribed motion)。

在LS-DYNA中，预设的刚体旋转轴位置为刚体的质心，若要强制刚体绕特定轴旋转，只要利用FEMB在PART选单里设定INERTIA相关的数值即可。

其中的XC、YC、ZC为旋转中心所在位置，即旋转轴通过的地方(内定值为刚体质心)。TM为Translational Mass，一般状况输入质量即可。

IXX~IZZ为此刚体质心位置相对于旋转中心的惯性矩，一般CAD软件应已具备此项计算功能，因此使用者无须再额外进行人工运算。

以Pro/E的操作为例：

1.确认图文件的单位(公制、英制)。 2.于旋转中心建立坐标系。

3.选择「Analysis」?「Model Analysis」?「Model Mass Properties」。 4.选择步骤2.产生之坐标系。

Pro/E计算结果会有二组Inertia Tensor，其中「Inertia at Center of Gravity with respect to XXX coordinate Frame」才是我们要的，再将相对应的数值填入IXX~IZZ即可。

27. 在LS-POST如何显示应变分布图？

解答：LS-DYNA内定的输出控制并不包含应变值。如果您希望能有应变值的输出的话，请在FEMB里依照下列步骤完成设定即可。

步骤1：点选「DynaMisc」，进入「DATABASE」选项。

步骤2：请在「Extent」选项中新增一个「Binary」的输出控制。 步骤3：设定STRFLG=1即可。

步骤4：分析完成后，使用者可直接在LS-Pre/Post中，点选 「Fcomp」?「Strain」观察计算后的应变分布图。LS-Pre/Post还有其它关于应变的显示按键，如：「Infin」、「Green」及「Almans」等，基本上这些应变量是根据d3plot档案中所记录的位置资料，经过微分计算而得到。如果点位取得够密的话，理论上会与板壳中性平面的应变值相当接近。

28. 强制结构体绕着特定轴旋转的时候，为何会出现奇怪的变形？

解答：在定义物体旋转自由度的时候，常会以*Initial_velocity_generation来设定角速度，然而这样的方式会将离心的效果纳入考量，特别是当角速度值很大的时候，，旋转物体会由旋转中心沿半径方向膨胀变形。

第 33 页