lsdyna笔记摘录--sph一些知识大杂烩
1，sph泥石流材料
*MAT_ELASTIC_PLASTIC_HYDRO
$# mid ro g sigy eh pc fs charl
3 1.13 0.027 1.2E-3 4.5E-5 0 0.0 0.000
$# eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
$# eps9 eps10 eps11 eps12 eps13 eps14 eps15 eps16
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
$# es1 es2 es3 es4 es5 es6 es7 es8
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
$# es9 es10 es11 es12 es13 es14 es15 es16
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
*EOS_GRUNEISEN
3 0.257 1.85 1.07 0.0 1.07 0.0 0.0
1.0
2，sph-lagrange绑定
*CONTACT_TIED_SHELL_EDGE_TO_SOLID_OFFSET_ID
$# cid title
2lingbujian-gudingjiechu
$# surfa surfb surfatyp surfbtyp saboxid s bboxid sapr s bpr
1 2 4 3 1 1
$# fs fd dc vc vdc penchk bt dt
0.5 0.5 0.0 0.0 20.0 0.01.00000E20
$# sfsa sfs b sast s bst sfsat sfs bt fsf vsf
1.0 1.0 1.0 1.0 1.0 1.0
3，sph-lag接触---
*CONTACT_AUTOMATIC_NODES_TO_SURFACE_THERMAL_ID
$# cid title
1contact
$# surfa surfb surfatyp surfbtyp saboxid s bboxid sapr s bpr
1 1 4 3 1 1
$# fs fd dc vc vdc penchk bt dt
0.5 0.5 0.0 0.0 20.0 0.01.00000E20
$# sfsa sfs b sast s bst sfsat sfs bt fsf vsf
1.0 1.0 1.0 1.0 1.0 1.0
$# k frad h0 lmin lmax ftosa bc_ fl g algo
6.70000E-6 0.2 1000.0 0.9
4，2d-sph-lag-sph2d-lagrange接触
使用part集 合定义；
*CONTACT_2D_NODE_TO_SOLID_TITLE
$# cid title
1
$# sph solid tbirth tdeath
-2 -1
$# soft unused vc offd pen fs fd dc
0 1.0 0.24 0.24 0.0
*SET_PART_LIST_TITLE
tool-set
$# sid da1 da2 da3 da4 solver
1 0.0 0.0 0.0 0.0MECH
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1 0 0 0 0 0 0 0
*SET_PART_LIST_TITLE
sph-part-set
$# sid da1 da2 da3 da4 solver
2 0.0 0.0 0.0 0.0MECH
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
2 0 0 0 0 0 0 0
5，sph-euler耦合
不支持
6，fem-sph绑定
【lspp，sph边界的节点集 合，fem边界的sgment集 合，tie_nodes_surface】
fem-sph绑定结果如下：
7，Sph-sph接触
*DEFINE_SPH_TO_SPH_COUPLING and sets the parameter ISOFT to 1 (Soft constraint formulation).
8，sph不飞溅？
如果定义了材料失效，就不会飞溅了。
反之，不定义失效，则具有飞溅效果
9，其他
17，sph，控制sph，sph-fem渗透，sph总结
1）尽量定义box【但是又不要太小，有模型的地方都包裹住】
2）Memory500合适
3）Form=15用于流体物质的sph
4）Maxv=超过该速度的sph粒子被停用/抑制。
5）如果sph是基于实体中心创建的，那么使用rigidwall的坐标=0即可，不用添加偏置。
如下图，刚性墙中心坐标为0，则刚性墙于sph粒子之间已经有间隙了，要求间隙大于粒子半径即可。如果sph粒子是紧贴着对称面，也就是sph粒子在x=0面上存在，那么刚性墙与sph粒子就需要间隙了。
6）当然也可以使用sph对称关键字。
7）Sph-fem接触如下：直接使用part或者part集 合即可。
8）sph-sph接触无需定义即可；
9）Sph对称边界的sph粒子，可以定义一下spc约束；
1）sph-fem尽量使用绑定接触，sph粒子和fem的segment定义tied绑定；
2）sph-fem中fem尽量使用刚体；
3）Sph单元
