请问哪位有典型算例－水箱跌落的命令流啊？

franck

我这里只有sysh320 大师 留下的一部分 命令流，请问哪位有全部的命令流？？谢谢了！
最后给出一个典型算例－水箱跌落的部分关键字：
*KEYWORD
*TITLE
boxwater2.k: dropping a water box onto a rigid platform
$========================================================================
$ [1] EXECUTION CONTROLS
$========================================================================
*CONTROL_TERMINATION
$ ENDTIM ENDCYC DTMIN ENDENG ENDMAS
0.0500000 0 0.0000000 0 0.0000000
*CONTROL_TIMESTEP
$ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MS1ST
0.0000000 0.2000000 0 0.0000000 0.0000000 0 0 0
*CONTROL_ENERGY
$ HGEN RWEN SLNTEN RYLEN
2 2 2 2
$========================================================================
$ [3] OUTPUT CONTROLS
$========================================================================
*DATABASE_BINARY_D3PLOT
$ DT CYCL LCDT BEAM
0.0005000 0
*DATABASE_GLSTAT
0.0001000
$========================================================================
$ [5] |SECTIONS|PARTS| DEFs
$========================================================================
*PART
water in the box
$ PID SECID MID EOSID HGID GRAV ADPOPT TMID
1 1 1 1 0 0 0 0
*SECTION_SOLID_ALE
$ SECID ELFORM AET
1 12
$ AFAC BFAC CFAC DFAC START END AAFAC
0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
$-------------------------------------------------------------------------------
*MAT_NULL
$ MID RHO PC MU TEROD CEROD YM PR
1 1000.0000 -1.000+10 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
*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