LS_DYNA常见问题解答与技巧分享

一、如何得知 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 的求解方式，可以为您解决静力、准静力、掉落、冲击、模态分析、振动、预应力/预应变、流固耦合、流场、热应力....等等的力学分析，使用者不必再额外添购模块程序，徒增预算的浪费与软件学习的负担。

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

解答：

由于 LS-DYNA 是全球使用率最高的结构分析软件，因此其它的前处理软件与 LS-DYNA 是完全兼容的。

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

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

解答：

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

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

解答：

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

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

五、 FEMB 能够自动产生有限元素网格吗？

解答：

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

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

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

解答：

这样的问题事实上是属于「Engineer Judgement」的层次，而非软件层次的问题。

使用者当然可以将完全按照 CAD 几何来建立有限元素的模型，然而任何一位有经验的工程师一定会先将不必要的几何造型予以忽略，以提升分析效益，并将分析课题明确化。

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

势流科技除了提供软件分析技术之外，同时更致力于协助业界将 CAE 导入于生产流程的正确位阶，培养正确的 CAE 分析观念，以提升产品研发品质，缩短生产流程。

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

解答：

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

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

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

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

2) 若是发生降伏产生大变形，在厚度方向上即使分割 8 层以上的实体元素，其掌握材料塑性应变的收敛速度仍然相当缓慢，而付出的时间代价为 shell element 的 25 倍以上。厚壳元素的运算速度及精度与三层实体元素相当，然而就塑性应变的精度而言，厚壳元素的表现比实体元素要好太多。

八、在 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 XXXcoordinate Frame」才是我们要的，再将相对应的数值填入 IXX~IZZ 即可。

九、在 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 档案中所记录的位置资料，经过微分计算而得到。如果点位取得够密的话，理论上会与板壳中性平面的应变值相当接近。

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

解答：

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

因此，如果不希望有离心膨胀的效果的话，可以利用*Load_body_rz(or rx、 ry)，并配合*Define_curve 给定角速度对应时间的关系曲线。

几个要注意的地方如下：

步骤 1：起始速度的设定选项，请设定 Phase=1，以进行 dynamic relaxation。

步骤 2：定义角速度的曲线中，必须开启 stress initialization 的功能(SIDR=1)。

步骤 3：在 control card 的部份须开启 accuracy 及 dynamic relaxation 的功能。

步骤 4：参考范例请至下载区 download。

十一、 LS-DYNA 该如何进行振动模态的分析设定？

解答：

模态分析即为特征值分析，须采取 Implicit 的方式求解，您只需要在前处理程序 FEMB 里依照下列方骤进行设定即可。

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

步骤 2：将一般 Implicit 功能选项开启，其中一项即为「IMP.EIGENVALUE」。

步骤 3：最基本的设定，您只需要设定模态个数(NEIG)的值即可。

步骤 4： LS-DYNA 针对特征值，还有许多非常实用的输出控制，详细用法请参考 LS-DYNAKeyword User’s Manual Vol. I，或是直接与势流科技的工程师联络。

十二、究竟采用 Implicit 及 Explicit 的分析方式会有多大的差异性？

解答：

右图为材料试片拉伸试验的仿真，分别采用 Implicit 及 Explicit 的分析方式。

依照静力的观念，其破坏点必发生于试片中点；然而若是以瞬间拉扯试片的话，因为应力传递速度与路径的关系，其破坏点未必会发生在中点，一般的静力软件必定不能满足真实的结构行为。

事实上， ASTM 对于材料性质试验，都有相关的规范以定义外力加载的速率，也就是因为这个原因。

坊间的结构动力学及数值分析相关参考书，对于 Explicit / Implicit 将有更详细的说明。IMPLICIT EXPLICIT

十三、要如何设定才能让 LS-DYNA 自动进行 Implicit / Explicit 的切换？

解答：

在 LS-DYNA 960 的版本即有 Implicit / Explicit 切换的功能，使用者可在FEMBDynaMiscIMP.General 的第一个字段输入「-100」，表示 Implicit / Explicit 的切换时间是依据#100 的 load curve 来定义的。而在 LS-DYNA 970 最新版本更提供了自动切换的功能，做法如下：

步骤 1：进入 FEMBDynaMiscCONTROL

步骤 2：进入 IMP.AUTO 选项，给定 DTEXP 的值，此值即为 explicit 分析所需的计算步长。

步骤 3：进入 IMP.General 选项，设定 IMFLAG=4 or 5 即可。

十四、您好，我们评估了很多结构软件，发现 LS-DYNA 是唯一能够将 Implicit 与 Explicitcode 整合在一起的软件，这对实际的分析工作有何帮助？

解答：

LS-DYNA 将 Implicit 及 Explicit 两种不同的求解方式整合在一起，对于初学者而言可以降低

学习的困扰、缩短软件学习的时程、省去面临不同软件的适应问题；对进阶的使用者而言，

可以从事更详实的力学分析。最实际的效益在于：

1) 预变形及预应力的仿真：真实的结构行为通常都是以连续作动的方式在进行，因为第一阶段的力学行为而衍生出后续阶段的变形。传统的作法，使用者必须凭借个人经验或是以简化的方式在结构体上施予一假设的预应力(预变形)量，然而所施加的量值却未必能符合真实的数值。 LS-DYNA 整合 Implicit 及 Explicit 解法，可以忠实的仿真各种结构反应。以卡勾卡合的动作机制而言，在第一阶段可利用 Implicit 解法指定卡勾运动方向，接着可马上切换至 Explicit 解法，仿真卡合瞬间的反应。

2) 相信以往从事 Implicit 分析的工程师都有调整收敛参数的不愉快经验，而 LS-DYNA970 可以视运算过程收敛性的难易程度，自动将 solver 由 Implicit 切换到 Explicit，以克服收敛性的问题，这个功能绝对是非常吸引人的！！

十五、请问以 Explicit 方式进行分析，影响运算时间的因素有哪些？

解答：

采用 Explicit 方式最直接影响运算时间的因素是「计算步长」 (time step)的大小， LS-DYNA决定计算步长的方式为：另外还有许多因素会影响运算的时间：

1) 分析历时(Termination Time)，分析历时较短者，可缩减运算时间。

2) 接触算法： LS-DYNA 采用最有效率的接触演算理论。

3) 机器等级：计算机等级高者，运算时间较短。

4) 元素理论：采用较多高斯积分点的元素，其运算时间会增加。

5) 材料模式之复简：材料性质直接影响计算步长。

6) 模型的复简：适当的简化 CAD model 可大幅提升运算品质。

所以当您在评估软件时，若有任何关于 LS-DYNA 计算速度的疑问，都欢迎与势流科技联络，我们将提供最正确的 CAE 信息与您分享。

十六、势流科技您好，请问弹簧/阻尼系统该采用什幺元素来仿真比较恰当？

解答：

若您是古早版本(LS-DYNA 940)的使用者，可以直接选用 discrete element 来仿真弹簧或是阻尼元素。

不过 LSTC 原厂建议采用最新版本 LS-DYNA 970 所提供的 Beam Element(#6)来仿真 Spring /Damper，配合*MAT(#66、 #67)可一次定义多轴的劲度与阻尼值，并且不用额外定义 LumpMass。

例如右图即是采用单一 Beam Element(#6)来仿真单自由度(SDOF)系统在 No Damped、 UnderDamped、 Criticle Damped 及 Over Damped 状况下的振动行为。

十七、请问安装 LS-DYNA(MPP)版本时，操作系统需有哪些额外的设定？

解答：

若你所采用的操作系统是 Window2000，可以请贵单位的网管人员协助下列设定：

1) 确定参与平行运算的各机器皆处于同一网域内。

2) 各机器须提供公用帐号及密码给予 MPP 的使用者登入。

3) 于每一台机器上安装最新版本的 MPICH 程序(Mpich-1.2.5)，利用所附的 mpiregister 程序设定 MPP 的使用者帐号及密码。

若您所用的系统是 UNIX / Linux，亦请贵单位的网管人员协助以下的设定：

1) 确定所有参与 MPP 计算的机器皆位于同一网域内。

2) 建立 NFS server 让 cluster 内的所有机器将 LS-DYNA 的所在目录挂载进来。

3) 进行相关安全设定(/etc/hosts、 /etc/hosts.equiv、 .rhosts....etc)，让 MPP 的使用帐号可以利用 rsh 或是 ssh 的方式在 Cluster 的机器间互相登入而无须输入密码。

4) 安装 MPICH 或是 LAM 最新版本(LAM-6.5.6)

十八、在 Linux 系统下如何指定 MPP 的 computing node？

解答：

若你是采用 MPICH 的话，请修改/usr/local/mpich/share/machines.LINUX 的档案即可，参考格式如下：

node1 :2

node2 :2

若你是采用 LAM 的话，请先编辑 hostfile，参考格式如下：

node1 cpu=2

node2 cpu=2

执行运算前请别忘了要先激活 LAM，请执行 lamboot –v hostfile；若是激活成功，屏幕上会显示参与 MPP 计算的 Node number。

十九、请问在执行 MPP 运算时，有没有类似 PC 版本 Ctrl+C 的功能以进行不同的 switch？

解答：

当然是有的。你可以利用另一个 terminal 窗口在工作目录产生名为「d3kil」的档案，将正确的 switch(如 sw2)填写在第一行， LS-DYNA 程序会周期性的读取 d3kil 档案，并进行相关的switch report。

若是 Linux / UNIX 系统的话，还可以有另一种方式：请另外开启新的 terminal，执行 sourcebg_switch，此时正在运算的程序会先暂停，并且等待你的 switch 讯号。

二十、如果我想要了解特定截面的力量分布情形，该设定哪些项目呢？

解答：

您所想获得的信息应该是属于 Section force 的输出结果。其做法如下：

1) 进入 Set 选单，将特定截面上的所有结点设成同一个群组(Node Set)。

2) 将这些结点所在的元素(在截面的同一侧)也设成同一个群组(Element Set)。

3) 进入 DynaMisc 选单点选 DatabaseCross section建立 Cross section set，将步骤(1)(2)所设定的结点及元素群组点选进来即可。

4) 最后再到 ACSII 选单定义 secforc 的记录输出时距即可

至于资料的读取方式与一般利用 LS-Pre/Post 读取 ACSII data 的作法一样。

一般常见的应用在于获得支承的扭转性质(如 LCD、铰接托架等结构)。

二十一、请问 LS-DYNA 在进行 Implicit 分析时，有没有什幺该注意的地方？

解答：

您这个问题好象有点「蛮广泛」

在这边先跟您说明一下大原则： Implicit method 分析主要就是在求解矩阵式，给定的力学条件不外乎是「Force Prescribed」与「Displacement Prescribed」，一般最常犯的毛病就是边界条件给定不足，造成 model 发生「刚体运动」的情形，您可以将这种情形想象成物体浮在空

中一样，宜尽量避免这种情形产生。在此提供您一个判断的方式：您可以先利用 LS-DYNA eigenvalue analysis 的功能先进行模态分析，藉后处理程序检查有无刚体模态，以判断 Implicit model 的束制条件是否充足。

二十二、请教一下， LS-DYNA 有没有什幺功能可以让上/下表面的网格移到中性平面来？

解答：

有的，请依照下列步骤进行：

1) 利用 Check 的功能，将元素法向量调整一致。

2) 进入 Element 选单，点选 transform。

3) 接着会出现向右图一样的窗口，请选择 NORMAL OFFSET。

4) 输入半个壳厚即可。

Normal offset 与直接平移的不同是在于 Normal offset 会依照壳厚修正曲率，能够让网格正确的移至中性平面的位置。

所以你大可放心在 CAD 的外表面产生 iges surface，再利用 Normal offset 的功能让网格产生在中性平面上，以加快建构网格的时间。

二十三、请问我该如何在后处理程序 LS-Pre/Post 里头看到壳元素的厚度？

解答：