Ls-dyna中无网格法的一点心得

1、Ls-dyna 中无网格法的一点心得前处理：. 采用 ls-prepost在 ls-prepost 中第 7 页有 sphgen 的功能， 利用这个功能可以建立一些简单的无网格模型。操作具体步骤是：1 选择 sphgen，出现如图1 的菜单，选择其中的new，表示新建一个无网格的命令。可以看到下面有几个选项：box， sphere， cylinder ， sketch图 12 选择 box，出现如图2 所示的菜单，即为创建一个具有长方体形的无网格模型。图 2可以看到有坐标x，y，z，而且标号为x1， x2， y1，y2，z1，z2，表示以p1 点（ x1，y1，z1）起点，以向量（由（x1，

2、y1，z1）指向（ x2， y2， z2）为对角线创建一个长方形的box。NX,NY,NZ 表示 X,Y,Z 方向的 sph 粒子的数量。 Box name 和 part id 都需要用户制定，当给定part的 id，比如“1”后，即表明创建的sph 单元属于part1 ，然后点“enter ”，可以看到material id也被激活，可以看到如图3 的标识，点击这个按钮，就可以为这个part选择相应的材料了。图 3点这个按钮出现如图4 的对话框，选择一种材料，如elastic，点击 edit ，出现如图5 的窗口图 4图 5可以选择add，或者 newid ，新建一种材料， 比如点击add，

3、可以看到mid 的选项下面出现了 1 的标号，代表材料id 为“ 1”，然后填写相应的属性，填写完成选择“accept”，表示接受这种材料类型。如果还要创建，那么选择“add”，如果想保持当前的数据，只是想把材料id 改为别的数值，那么点击“newid ”。创建了材料之后，选择done 推出当前窗口，再选择done 结束材料的创建。在 sph 的主窗口下，选择 apply ，可以看到我们创建的中，如图 6 所示。如果不选择“ done ”，可以继续创建box 形的 sph sph 模型。模型，出现在窗口图 63 我们现在创建的box 使用的缺省尺寸1，1，1。用户可以根据需要采用不同的尺寸和s

4、ph的数量。4 现在所创建的仅仅是sph 的几何模型，对于属性的定义还要采用其他的关键字，比如section_sph 。下面就section_sph 的一些参数谈一下应用。在 ls prepost 的第三页有各种关键字的定义。选择“*section ”，进行单元属性的定义。在弹出的菜单中选择sph，点击 edit ，我们进行sph 单元属性的定义。与材料定义一样，出现了如图7 的窗口， 点击add，表示新建一个*section_sph的属性。图 7需要填写的参数有：CLSH：这是用于确定 smooth length 的一个常数。 通常用缺省的就行。这个常数的用途可能是这样的： 在开始的时候，

5、程序根据当前的粒子的间隔确定一个例子间隔的最大值，把它作为sph 法的最初的光滑长度， CLSH 是最初的光滑长度的系数，有点类似于显示积分中步长的确定，如果实际步长大于初始光滑长度，那么采用CLSH给初始光滑长度增加倍数，以满足程序的需要。ii.Hmin，Hmax 和CLSH的原理是一样的，都是调整smooth length的参数。SPHINI是强制的初始步长， 如果采用这个步长， 那么程序中就不会自动计算步长，而是一直使用这个 SPHINI的大小作为 SMOOTH LENGTH。Death： sph 粒子消亡的时间Start： sph 粒子开始的时间这两个参数， 是设定 sph 算法的作用

6、时间的，可以设定什么时候sph 法有作用，什么时候sph 法结束。设定完这些参数后，点击“accept ”，表示接受这些设定。然后点击done推出section_sph 的定义。5 把所创建的 *section_sph 付给相应的part 。还是在当前页，选择“*part ”，如图 8。选择 part 【 1】，点 edit ，弹出如下窗口，如图9。图 8图 9把 secid 给值为 1，表示将刚才创建的section 付给了 part1 。出现提示，点击replace。这样我们就完成了一个sph 有限元模型的创建过程，包含几何， 单元属性， 材料属性的有限元模型。6 Control_sph

7、的用法在当前页点击“*control ”，在当前窗口中选择sph，点击 cre/edit ，进入“ control_sph ”的创建菜单，如图10 所示。图 10关于其中参数的说明：NCBS，表示搜索 sph 时的循环次数Boxid，表示在这个 box 里面 sph 参与计算，而出去了就不对模型有贡献。这两个的意思有些类似于 contact 的搜索过程iii.Dt 和 section_sph 中的 death 一样iv.Dim，设定有限元模型的维数，3 表示三维有限元模型，2 表示平面问题，2 表示2 维轴对称。对于2d 的应用而言，模型必须建立在xy 平面中。Start 和 section_

8、sph 中的功能一样MAXV，粒子允许的最大速度，如果超过了这个速度，粒子就会失效。7 在 ls-dyna970 中还增加了 symmetry （对称条件的定义） ，sph_flow 的设定， 我没有试过，以后有时间我会试一下的。二后处理Sph 的后处理可以以下几步实现。1 定义 ASCII关键字中的sphout2 通过 *database_history_sph(_set) ，定义输出那些sph 单元。3 利用 ls-prepost 中第一页的ASCII菜单，把 sph 调入到后处理中来。如图11， 12 所示，在 file 右侧的列表中选择sphout ，然后点load ，载入 sph 的结果文件。其