输入的本构曲线和输出的结果有比较大的差异？ 75

       采用混凝土损伤塑性模型做简单的单轴压缩试验。当采用平面应力单元验证本构的时候，输入曲线和输出曲线吻合（图1）；但是换成3维模型后（100*100*100mm立方体、150*150*300mm棱柱体、φ50*100mm圆柱体），输出的应力-应变曲线和输入曲线不一致（图2）。

       模拟中约束了试件底部的竖向自由度，然后在试件顶部施加竖向位移（图3）。最后提取试件底部所有节点反力来计算应力，提取试件顶部竖向位移来计算应变，以此绘制应力-应变曲线。

       1.通过对比输入曲线可以看到（图2，模拟得到的应力-应变曲线峰值强度略低，且承载力下降得较快。

       2.立方体试件的曲线下降得最快，而圆柱体试件的曲线与棱柱体接近，但有一个应力下降又上升的现象。

       请问这是怎么回事呢？