非线性有限元经典例题

非线性有限元经典例题

厚壁圆筒受内压力P＝180Mpa作用，假设为平面应变问题，圆弹模E 2.1×10E5Mpa，泊松比 0.3 单轴屈服应力 240MPa，应变硬化参数 H＝0，采用MISES 屈服准则进行计算。 /prep7

ET,1,plane42

mp,ex,1,2.1e4

mp,nuxy,1,0.300000

!！ 你用tb命令输入，很好用的，我可以给你一个：

!tb,bkin, !经典的双线性随动强化

TB,BISO,1, , , , !双线性等向强化

tbdate,,24,0 !tbdate,,屈服应力,斜率

alls

CYL4,0,0,200,0,100,180

MSHAPE,0,2D

MSHKEY,1

LSEL,S,LENGTH,,100

LESIZE,ALL, , ,16, , , , ,1

ALLS

LSEL,U,LENGTH,,100

LESIZE,ALL, , ,30, , , , ,1

amesh,all

CSYS,1

DDELE,ALL,ALL

NSEL,S,LOC,Y,0

NSEL,A,LOC,Y,180

DSYM,SYMM,Y,1, !定义对称边界，施加约束

NSEL,S,LOC,Y,90

!NSEL,A,LOC,Y,180

!NSEL,R,LOC,X,200

D,ALL, , , , , ,UX, !定义不要侧向移动

DDELE,ALL,ALL

NSEL,S,LOC,Y,90

D,ALL, , , , , ,UX,

NSEL,S,LOC,Y,0

NSEL,A,LOC,Y,180

D,ALL, , , , , ,UY,

/solu

ALLS

PRED,ON !预测器选项

OUTRESS,ALL,ALL !输出控制

!SFCUM,PRES,ADD,1,1, !荷载迭加开关

SAVE

CSYS,1

NSEL,S,LOC,X,100

SF,ALL,PRES,18 !施加内压

ALLS

ANTYPE,0!定义结构

!solcontrol，OFF !指定是否使用一些非线性求解缺省值

NLGEOM,ON !考虑大变形影响0,NO

AUTOTS,0!使用自动时间步长

NSUBST,3,10,3 !指定此荷载步的子步数

TIME,1

SOLVE

结果：径向应力和环向应力分别为

总结：与理论解相当的接近。