一个长2 m,宽1 m的矩形板,中心开孔半径r =50mm。水平两端受拉力,q =50 MPa。弹性模量E =210 GPa,泊松比μ =0.3。计算孔周围的应力分布和应力集中系数。
计算结果:
命令流:
/PREP7
ET,1,PLANE183 !单元类型
KEYOPT,1,3,0 !平面应力
MP,EX,1,2.1e5 !材料属性
MP,PRXY,1,0.3
RECTNG,,1000,,500, !几何体
PCIRC,50, ,0,360,
ASBA,1,2 !布尔操作面减去面
AESIZE,all,50, !面的单元尺寸
MSHAPE,0,2D
MSHKEY,0
AMESH,all !面网格划分
LREFINE,5, , ,3,1,1,1 !在线周围细化网格
NSEL,S,LOC,X,0 !由坐标位置选节点
D,all,UX !位移约束
NSEL,S,LOC,Y,0
D,all,UY
ALLSEL,ALL
SFL,2,PRES,-50, !拉力载荷
FINISH
/SOL
ALLSEL,ALL
SOLVE !求解
FINISH
/POST1
PLNSOL, S,X, 0,1.0 !直角坐标系SX应力
/EXPAND,4,POLAR,HALF,,90 !1/4扩展
PATH,p1,2,30,20, !定义路径
PPATH,1,0,,50,,0,
PPATH,2,0,,500,,0,
/PBC,PATH,1
/REPLOT
PATH,P1
PDEF,p1-sx,S,X,AVG !路径应力
PLPATH,P1-SX
RSYS,1 !结果柱坐标系
PLNSOL, S,Y, 0,1.0 !柱坐标系SY应力(环向)