圆管受压里的问题

wugang2011

用ls-dyna结构分析圆管受压模型，采用了实体单元。在后处理器中提取了圆管的内能、动能、沙漏能和位移能，这四项加在一起应该就是圆管在压缩过程中吸收的总能量吧，然后提取了冲击板收到的反力，通过反力对位移的积分算出冲击板对圆管做的功。现在遇到的问题是，这冲击板对圆管做的功不等于前面说的四项能量相加所得到的的总能量！？？？
为什么这两者不相等？和我采用实体单元有关吗？实体单元有什么地方要设置吗？整个压缩过程沙漏能和位移能都较小。（其实我也用过壳单元建立这个模型，当时的总能量和功是相同的，我不知道为什么换成实体单元就不同了~~~~~）

求助！

zhangyanbin1236

原文件在哪里啊