*EOS_LINEAR_POLYNOMIAL $ EOSID C0 C1 C2 C3 C4 C5 C6

1 0.0000000 1.50000+9 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 $ E0 V0

0.0000000 1.0000000

$======================================================================== *PART

void portion in the box 2 1 1 1 0 0 0 0 *INITIAL_VOID_PART 2

$======================================================================== *PART

rigid box containing water

$ PID SECID MID EOSID HGID GRAV ADPOPT TMID 3 3 3 0 0 0 0 0 *SECTION_SOLID $ SECID ELFORM AET 3 0 *MAT_RIGID

3 2000.0000 1.00000+8 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000

0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 $======================================================================== *PART

rigid super-heavy platform

$ PID SECID MID EOSID HGID GRAV ADPOPT TMID 4 4 4

*SECTION_SHELL

$ SID ELFORM SHRF NIP PROPT QR/IRID ICOMP 4 0

$ T1 T2 T3 T4 NLOC 0.011 0.011 0.011 0.011 *MAT_ELASTIC

$ MID RHO E PR DA DB K 4 1000000.0 1.0000+14

$======================================================================== $ [8] BC's + IC's + BODY LOADS + FORCE FIELDS

$======================================================================== *INITIAL_VELOCITY $ NSID NSIDEX BOXID 0

$ VX VY VZ VXR VYR VZR 0.0 -20.0 0.0

$------------------------------------------------------------------------------- *LOAD_BODY_Y

$ LCID SF LCIDDR XC YC ZC 1 1.00 *DEFINE_CURVE

$ LCID SIDR SFO OFFA OFFO DATTYP

1

$ X=abcissa Y=ordinate 0.0 981.0 1.0 981.0

$======================================================================== $ [9] LAGRANGIAN CONTACTS CONSTRAINTS, ...

$======================================================================== $ SFS = scale fact on dflt SLAVE penal stifns (see CONTROLL_CONTACT) $ SFM = scale fact on dflt MASTER penal stifns (see CONTROLL_CONTACT) *CONTACT_AUTOMATIC_NODES_TO_SURFACE

$ SSID MSID SSTYP MSTYP SBOXID MBOXID SPR MPR 3 4 3 3

$ FS FD DC VC VDC PENCHK BT DT

$ SFS SFM SST MST SFST SFMT FSF VSF 100. 100.

$======================================================================== $ [10] EULERIAN & ALE CONTACTS CONSTRAINTS, ...

$======================================================================== *CONTROL_ALE

$ DCT NADV METH AFAC BFAC CFAC DFAC EFAC 2 1 4-1.0000000 0.0000000 0.0000000 0.0000000 $ START END AAFAC VFACT VLIMIT EBC 0.0000000 0.0000000 0.0000000 0.0 *ALE_REFERENCE_SYSTEM_GROUP

$ SID STYPE PRTYP PRID BCTRAN BCEXP BCROT ICOORD 1 0 5 1

$ XC YC ZC EXPLIM

*SET_PART_LIST $ SID DA1 DA2 DA3 DA4 1

$ PID1 PID2 PID3 PID4 PID5 PID6 PID7 PID8 1 2

*ALE_REFERENCE_SYSTEM_NODE $ NSID 1

$ N1 N2 N3 N4 N5 N6 N7 N8 5 6 7

2、 SPH算法

SPH算法作为DYNA中第一种无网格（meshfree）算法，在连续体的破碎或分离分析中得到了广泛的关注和应用。在解决极度变形和破坏类型的问题上SPH有着其他方法无法比拟的优势，可以说无网格算法正在成为数值分析领域的研究热点，具有很好的发展前景。 我们知道传统的有限单元法中，单元的形状对结果的精度影响很大，如果单元因为变形过大可能造成矩阵奇异，使得精度降低甚至无法计算下去。而SPH算法则是把每个粒子作为一个物质的插值点，各个粒子间通过规则的内插函数计算全部质点即可得到整个问题的解。

主要的关键字如下： *section_sph

提供算法选择，以及sph粒子的滑顺长度的定义； *control_sph

提供sph算法的控制，如粒子排序后的循环次数、计算空间、中止时间以及维数； 处理sph粒子与其它结构的相互作用采用接触算法。 下面给出某一算例的部分命令流： *KEYWORD *TITLE sph test $

*DATABASE_FORMAT 0

$units:cm,gm,us

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ CONTROL OPTIONS $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*CONTROL_PARALLEL