Skip to content

ElemTest-未改进wilson非协调元PS4I&CPS4I对比

测试对象:FE3d/src/ps4.py的PS4I单元

测试案例:悬臂梁弯曲

abaqus不支持输出CPS4I单元的单元刚度矩阵无法直接对比,所以对比悬臂梁的位移。

材料参数:E=1000mPA,NU=0.15,t=1.0

悬臂梁:1.0mm*2.0mm*40.0mm

对比结果:位移计算结果小于1%

abaqus模型:

Img

abaqus-CPS4I单元计算结果:

Img

Img

Img

自编PS4I单元计算结果:

Img

Img

Img