橡胶圆筒轴向受力的问题，求助~

作橡胶材料分析，关键是获取正确的材料常数，我采用的双参数M-R模型，对于边界条件就是把外和内边界节点旋转到柱坐标系

/PREP7
!*
ET,1,PLANE182
!*
!*
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
!*
!*
TB,HYPE,1,1,2,MOON
TBTEMP,0
TBDATA,,0.4,1.3,1,,,
CYL4,0,0,0.1,0,0.25,360
ESIZE,0.01,0,
amesh,1
/SOL
LSEL,ALL
FLST,5,8,4,ORDE,2
FITEM,5,1
FITEM,5,-8
LSEL,S, , ,P51X
NSEL,ALL
NSLL,S,1
NPLOT
CSYS,1
FLST,2,160,1,ORDE,2
FITEM,2,1
FITEM,2,-160
!*
/GO
D,P51X, , , , , ,UX, , , , ,
LSCLEAR,ALL
FLST,2,160,1,ORDE,2
FITEM,2,1
FITEM,2,-160
FINISH
/PREP7
FLST,2,224,1,ORDE,2
FITEM,2,1
FITEM,2,-224
NROTAT,P51X
FINISH
/SOL
FLST,2,160,1,ORDE,2
FITEM,2,1
FITEM,2,-160
!*
/GO
D,P51X, , , , , ,UX, , , , ,
FLST,2,4,4,ORDE,2
FITEM,2,5
FITEM,2,-8
/GO
!*
SFL,P51X,PRES,3e2,
ALLSEL,ALL
NLGEOM,1
/STATUS,SOLU
SOLVE

应该在圆柱坐标系中求解，以便设置外径固定的约束；此外地面应该设置 Z 向的约束，以消除刚体运动。
这是一个轴对称问题，如果采用轴对称模型，问题将大大简化，如果应力－应变关系不太复杂的话，估计也可以采用弹性理论导出相应的理论解。

这个就在软件中按你的实际情况加载就可以了，建议参考一下类似的例子