求救：LS-DYNA 悬臂梁自由端加载位移，应变值偏

我用shell163单元模拟一个长40mm,宽7mm,克厚0.2mm的悬臂梁，在其自由端加载垂直位移，其变形情况符合挠度曲线，但是其应变值均比理论计算值小很多，请教高手指点。我用了S/R Hughes-Liu和Belytschko-wong算法，效果是一样的，积分点数是3。材料选用的各向同性弹性材料模型，密度输入7.8e-6，泊松比0.3，弹性模量2.1e5.单位采用的是mm-kg-s.附件是我的模型图