内容发布更新时间 : 2024/4/24 14:18:27星期一 下面是文章的全部内容请认真阅读。

materials or where the mesh densities of the two contacting surfaces are dissimilar.

When SOFT=1, we use the max of the stiffness as calculated by soft=0 and soft=1. Therefore reducing SOFSCL has no effect if the soft=0 stiffness is larger.

k = max(SLSFAC*SFS*k0, SOFSCL*k1)

where k is the penalty stiffness SLSFAC is user input on *CONTROL_CONTACT SFS is user input on *CONTACT card 3 SOFSCL is user input on *CONTACT optional card A k0 is the stiffness calculated from material bulk modulus and element dimensions k1 is the stiffness calculated from nodal masses and the solution time step.

Note: For two way contact like

*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE, replace SFS with SFM (user input on *CONTACT card 3) in the above equation.

13.LS-DYNA中夹层板(sandwich)的模拟

在LS-DYNA中夹层板(sandwich plate)通常有三种方式来模拟

1. 采用几层体单元来模拟，当板比较薄，沿厚度方向单元尺寸比较小，时间步很难提高，

2. 中间的夹层材料用体单元模拟，外面的包层材料用shell单元模拟，在体单元和壳单元之间用Contact_tied (_offset)卡建立粘接关系，

3. 用一层壳单元来模拟。对材料*Mat_composite_damage和*mat_enhanced_composite_damage打开层板理论选项(*control_shell卡LAMSHT项设置为1)

对最后一种方式，需要定义*intergration_shell卡片。例如下面这个层合板：

*INTEGRATION_SHELL 1,8,0

-.9722, .02778, 1 -.9167, .02778, 1 -.6667, .22222, 2 -.2222, .22222, 2 .2222, .22222, 2

.6667, .22222, 2 .9167, .02778, 3 .9722, .02778, 3 *PART material 1 1 11 *PART material 2 1 12 *PART material 3 1 11

*SECTION_SHELL

1 2 0.0000000 8.0000000 0.00000 -1.000000 0 18.000000 18.000000 18.000000 18.000000 0.0000000 *mat_composite_damage

11, 2.7e-6, 73.4, 73.4, 73.4, 0.32, 0.32, 0.32 27.8, 27.8, 27.8 1.e9, 1.e9, 1.e9, 1.e9 *mat_composite_damage

12, 6.3e-7, 0.286, 0.286, 0.286, .3, .3, .3 0.11, 0.11, 0.11 1.e9, 1.e9, 1.e9, 1.e9

－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－ 其中*intergration_shell卡片中第一项S表示板厚度方向积分点的规一化坐标值，从-1到

1。

积分点规一化坐标＝(积分点全局Z坐标－中面Z坐标)/(thickness/2)

第二项WF(Weighting factors)为与积分点关联的厚度与整个板厚度的比值，WF

和应该为1

。第三项指向各层对应的part ID号，这个part ID号仅用来指定这一层的材料。

完整的一个实例如下：shell8lam22.rar 14. LS-DYNA隐式－显式自动转换例子

在LS-DYNA中可以通过卡片*Control_implicit_general的选项imflag指定为一个负值(这个负值的绝对值为一条*define_curve ID)，来自动实现在指定的时间内作隐式或显式分析。

这里的例子为一个一端固定平板，另一端施加从0-0.01秒从0增加到10.0的力，之后释放。曲线ID3 指定从0-0.01秒内用隐式计算，增量步为20步。之后自动转为显式继续求解。

下载input deck－》implicit_explicit_switchk.txt

注意请用5434,5434a or 971版本运行上面的k文件，不要用6763单精度或者双精度版本(这两个版本在implicit求解方面似乎存在bug，隐式分析部分无法完成 )

14. 怎样进行二次开发

14.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和其它资源文件；通过编译连接后生成新的