【急求助】一个Dyna的算例，谢谢！

模型是一根简化的直梁结构，分为五段，前四段长度为56.2mm，最后一段为13.2mm。每段梁之间用零长度的弹簧连接，弹簧属性赋予了一个力与位移的曲线。想要不计梁单元的变形，因此设为刚性。整个模型与运动的刚性墙碰撞，我的想法是这样（参考他人论文）：应该是仅有弹簧属性中的力与位移的曲线发生作用，即第一根弹簧的力达到73500N时开始变形，当压缩长度到达56.2mm时失效，第二根弹簧开始起作用，以此吸收碰撞能。
一直计算不出结果。困扰本人多时了，在此谢过！
K文件如下：
$$ HM_OUTPUT_DECK created 10:24:12 12-22-2009 by HyperMesh Version 10.0build6
$$ Ls-dyna Input Deck Generated by HyperMesh Version : 10.0build60
$$ Generated using HyperMesh-Ls-dyna 971 Template Version : 10.0build60
*KEYWORD
*CONTROL_TERMINATION
$$ ENDTIM ENDCYC DTMIN ENDENG ENDMAS
0.03
*CONTROL_TIMESTEP
$$ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MSIST
0.0 0.67 4.4450E-07 0 0 0
1
*CONTROL_CONTACT
$$ SLSFAC RWPNAL ISLCHK SHLTHK PENOPT THKCHG ORIEN ENMASS
0.1 0.0 2 1 1 2
$$ USRSTR USRFRC NSBCS INTERM XPENE SSTHK ECDT TIEDPRJ
4.0
$$ SFRIC DFRIC EDC INTVFC TH TH_SF PEN_SF

$$ IGNORE FRCENG SKIPRWG OUTSEG SPOTSTP SPOTDEL SPOTHIN
1
*CONTROL_CPU
$$ CPUTIM
1.0000E+09
$$DATABASE_OPTION -- Control Cards for ASCII output
*DATABASE_GLSTAT
1.0000E-03 1
*DATABASE_MATSUM
1.0000E-03 1
*DATABASE_NCFORC
1.0000E-03 1
*DATABASE_NODFOR
1.0000E-03 1
*DATABASE_NODOUT
1.0000E-03 1
*DATABASE_RCFORC
1.0000E-03 1
*DATABASE_RWFORC
1.0000E-03 1
*DATABASE_BINARY_D3PLOT
$$ DT/CYCL LCDT BEAM NPLTC
1.0000E-03
0
*NODE
1 0.0 0.0 0.0
2 56.2 0.0 0.0
3 56.2 0.0 0.0
4 112.4 0.0 0.0
5 112.4 0.0 0.0
6 168.6 0.0 0.0
7 168.6 0.0 0.0
8 224.8 0.0 0.0
9 224.8 0.0 0.0
10 238.0 0.0 0.0
11 -2.0 0.0 0.0
12 0.0 5.62 0.0
13 56.2 5.62 0.0
14 112.4 5.62 0.0
15 168.6 5.62 0.0
16 224.8 1.32 0.0
*MAT_RIGID
$HMNAME MATS 1beam
17.8000E-06 210.0 0.3
0.0
0
*MAT_SPRING_INELASTIC
$HMNAME MATS 2spring
2 1 1.0
*PART
$HMNAME COMPS 1beam
$HWCOLOR COMPS 1 8

1 1 1 0
$HMNAME COMPS 2spring
$HWCOLOR COMPS 2 7

2 2 2
*SECTION_BEAM
$HMNAME PROPS 1beam
1 10.83333333 -3 0.0 0.0 0.0
2.3 2.3 2.3 2.3 0.0 0.0
*SECTION_DISCRETE
$HMNAME PROPS 2spring
2 0 0.0 0.0 0.0 -56.2
56.2 0.0
*INTEGRATION_BEAM
$HMNAME PROPS 3int
3 12 0.1046 0 1
-52.5 29.85 0.083333 0
-17.5 29.850.08333333 0
17.5 29.85 0.083333 0
52.5 29.85 0.083333 0
-52.5 -29.85 0.083333 0
-17.5 -29.85 0.083333 0
52.5 -29.85 0.083333 0
52.5 -29.85 0.083333 0
-68.85 14.35 0.083333 0
-68.5 -14.35 0.083333 0
68.85 14.35 0.083333 0
68.85 -14.35 0.083333 0
*DEFINE_BOX
$HMNAME BLOCKS 1box
$HWCOLOR BLOCKS 1 11
1 -10000.0 10000.0 -10000.0 10000.0 -10000.0 10000.0
*CONTACT_AUTOMATIC_GENERAL_ID
$HMNAME GROUPS 3inter
$HWCOLOR GROUPS 3 3
3
0 5 1


*RIGIDWALL_PLANAR_MOVING_ID
$HMNAME GROUPS 1rigidwall
$HWCOLOR GROUPS 1 4
1
0 1
-2.0 0.0 0.0 -1.0 0.0 0.01.0000E-03
500.0 10000.0
*ELEMENT_DISCRETE
19 2 2 3
20 2 4 5
22 2 6 7
23 2 8 9
*ELEMENT_BEAM
14 1 1 2 12
15 1 3 4 13
16 1 5 6 14
17 1 7 8 15
18 1 9 10 16
$$
$$ Sets Defined In HyperMesh
$$
*SET_PART_LIST
$HMSET
$HMNAME SETS 1Set_1
1
1 2
*BOUNDARY_SPC_NODE
$HMNAME LOADCOLS 1spc
$HWCOLOR LOADCOLS 1 3
10 0 1 1 1 1 1 1
*DATABASE_HISTORY_NODE
$HMNAME OUTPUTBLOCKS 1rigid wall
11
*DEFINE_CURVE
$HMNAME CURVES 1curve1
$HWCOLOR CURVES 1 3
$HMCURVE 1 2 curve1
1 0 1.0 1.0 0.0 0.0 0
-56.2 -73500.0
-0.01 -73500.0
0.0 0.0
0.01 73500.0
56.2 73500.0
*END
$HMNAME BEAMSECTCOLS 1boxsection
$HMNAME BEAMSECTS
$HMCONT 1 1boxsection
$HMCONT 2 7 1 01.0 1.0 0.0 0.0 0.0
$HMCONT 0.0 0
$HMNAME BEAMSECTS BEAMSECTIONSTANDARD 0 6 0HMBox
$HMNAME BEAMSECTS BEAMSECTIONSTANDARD PARAMETERS 2.3 0.5 2.3
$HMNAME BEAMSECTS BEAMSECTIONSTANDARD PARAMETERS 140.0 10.0 140.0
$HMNAME BEAMSECTS BEAMSECTIONSTANDARD PARAMETERS 62.0 10.0 62.0
$HMNAME BEAMSECTS BEAMSECTIONSTANDARD PARAMETERS 2.3 0.5 0.5
$HMNAME BEAMSECTS BEAMSECTIONSTANDARD PARAMETERS 2.3 0.5 0.5
$HMNAME BEAMSECTS BEAMSECTIONSTANDARD PARAMETERS 2.3 0.5 0.5
$HMNAME BEAMSECTS END