【求助】金属切削过程中能量不守恒问题？有图有K文件。

tangrocpeng

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

1、材料定义【工件为45#钢(尺寸：500*50*15mm)，刀具为硬质合金YT15(尺寸：30*10*20mm)，夹具为自定义材料(尺寸：100*50*20mm),密度自定义，质量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

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的动能作用下，直到刀具停止有多少动能转化为内能的问题，但是就是不知道为什么系统计算后能量会增大。请大家看看，到底是哪里的问题？实在感激不尽！


切削模型及K文件如下：

tangrocpeng

我仔细查了一下，应该就是滑动能出现负值，造成内能增大的问题。搜索论坛以前的帖子，也有人碰到过这个情况，但是没人说明一下怎么解决啊？请版主们指点一二。

tangrocpeng

这几天论坛里关注“能量”方面的有很多。一般介绍的方法就是更改接触参数，控制穿透。但是对于切削方面来说，本来就是侵蚀接触，要删除失效单元的，那这个接触参数该如何调整呢？恳请高人解惑……