用dyna计算金属切削时能量不守恒的问题？

tangrocpeng

    模型由3部分组成：包括工件(mat1)、刀具(mat2)、夹具(mat3)。对应part为1、3、2。切削原理就就是给夹具施加10m/s的初速度从而带动刀具切削工件。切削深度为1mm。单位制为mm，kg，pa。
http://home.simwe.com/attachment/201007/19/473778_12795220698R2g.jpg
几个k文件关键参数如下：


1、材料定义【工件为45#钢(500*50*15)，刀具为硬质合金YT15(30*10*20)，夹具为自定义材料(100*50*20),密度自定义，质量50kg，保证初始动能是2500J】
*MAT_PLASTIC_KINEMATIC

1 0.786E+04 0.210E+12
0.270000 0.355E+09 0.768E+09 0.300


40.0
5.00
0.750

*MAT_ELASTIC

2 0.115E+05 0.510E+12
0.300000
0.0
0.0
0.0
*MAT_ELASTIC

3 0.500E+06 0.210E+12
0.300000
0.0
0.0
0.0
*HOURGLASS

1
6 0.100
0
1.50
0.600E-01
0.00
0.00

*HOURGLASS

2
6 0.100
0
1.50
0.600E-01
0.00
0.00

*HOURGLASS

3
6 0.100
0
1.50
0.600E-01
0.00
0.00
*PART
Part
1 for Mat
1 and Elem Type
1

1
1
1
0
1
0
0
$
*PART
Part
2 for Mat
3 and Elem Type
1

2

1
3
0
3
0
0
$
*PART
Part
3 for Mat
2 and Elem Type
1

3
1
2
0
2
0
0



2、刀具与工件为侵蚀接触：静、动摩擦系数分别为0.15和0.1
*CONTACT_ERODING_SURFACE_TO_SURFACE

1
3
3
3
0
0
0
0
0.15000
0.1000
0.000
0.000
0.000
0 0.000
0.1000E+08

1.000
1.000
0.000
0.000
1.000
1.000
1.000
1.000


0
1
0


3、控制选项如下

*CONTROL_ENERGY

2
2
2
2
*CONTROL_HOURGLASS

1 0.100

*CONTROL_SHELL

20.0
1
-1
1
2

2
1
*CONTROL_BULK_VISCOSITY

1.50
0.600E-01
*CONTROL_CONTACT

0.10000
0.00000
2
0
1
1
1

0
0
10
0
4.00

*CONTROL_TIMESTEP

0.0000
0.9000
0
0.00
0.00

*CONTROL_TERMINATION

0.500E-01
0
0.00000
0.00000
0.00000

      按理来说，模型是相当的简单，因为不需要考虑切削质量等问题，只是考虑初始2500J的动能转化为内能问题，但是就是不知道为什么系统计算后能量会增大。请大家看看，到底是哪里的为题？实在感激不尽！