ANSYS层合板命令

!采用shell99单元(线性层合单元)，模拟[0/45/-45/90/90/-45/45/0]层合板受单向拉伸,拉伸应力100MPa
!层合板尺寸：长90mm，宽25mm，厚1mm

/prep7
ET,1,SHELL99,,,,,2,4
KEYOPT,1,8,1
R,1,8,1
RMODIF,1,13,1,0, 0.125,1,45,0.125
RMODIF,1,19,1,-45,0.125,1,90,0.125
MP,EX,1,135000
MP,EY,1,8000
MP,EZ,1,8000
MP,PRXY,1,0.34
MP,PRYZ,1,0.34
MP,PRXZ,1,0.34
MP,GXY,1,4500
MP,GYZ,1,3700
MP,GXZ,1,4500

!建立层合板模型
rectng,-45,45,-12.5,12.5

!定义网格划分
LSEL,S,LENGTH,,90 !选择所有长度为90mm的线段(长度方向上的线段)
LESIZE,all, , ,10, , , , ,1 !长度方向等长划分10个单元
lsel,s,length,,25 !选择所有长度为25mm的线段(宽度方向上的线段)
lesize,all, , ,20,-10, , , ,1 !宽度方向划分20个单元，中间疏，两边密

!映射划分层合板
MSHAPE,0,3D
MSHKEY,1
AMESH,all

!定义约束
/solu
nsel,s,loc,x,0 !选择所有x=0的节点
D,ALL, , , , , ,UX !给选中的节点施加x方向的约束(无x方向的位移)
nsel,s,loc,y,0 !选择所有y=0的节点
D,ALL, , , , , ,UY !给选中的节点施加y方向的约束(无y方向的位移)
allsel,all


!施加拉伸载荷
lsel,s,loc,x,45
lsel,a,loc,x,-45
SFL,ALL,PRES,-100, !施加100MPa的拉应力(pressure)

!求解设置
allsel,all !求解前选择所有元素
/solu
solve

!计算结果分析(因为是shell单元，故只能计算面内应力，不能考虑层间应力，在边界处没有边缘效应)
/post1
layer,1 !取第一层(从下往上),即0度层结果
PRNSOL,S,COMP !list 结果
layer,2 !取第二层(从下往上),即45度层结果
PRNSOL,S,COMP !list 结果
layer,3 !取第三层(从下往上),即-45度层结果
PRNSOL,S,COMP !list 结果
layer,4 !取第四层(从下往上),即90度层结果
PRNSOL,S,COMP !list 结果
layer,5 !取第五层(从下往上),即90度层结果
PRNSOL,S,COMP !list 结果
layer,6 !取第六层(从下往上),即-45度层结果
PRNSOL,S,COMP !list 结果
layer,7 !取第七层(从下往上),即45度层结果
PRNSOL,S,COMP !list 结果
layer,8 !取第八层(从下往上),即0度层结果
PRNSOL,S,COMP !list 结果
!注：以上所得结果与理论解完全符合


RMODIF那里表面了0/45/-45/90/90/-45/45/0