!采用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