非軸對稱荷載下軸對稱問題的有限元分析 - 道客巴巴 (doc88.com)

最近做了個例子，關于軸對稱結構受非軸對稱載荷情況下采用二維模型進行計算的，文中說對每個截面的應力分布均可以求得，但我實際計算只求得0度截面的應力分布同三維模型結果接近，45度的計算結果差距較大，想請教大家，如何求得任意角度截面的應力分布？其原理又是啥？

2Dmesh.rar

給出模型及求解及后處理的命令流

!靜力計算，直線工況，只加垂向力

/FILNAME,file,0

/solu

allsel,all

lsclear,all

DL,2,,ALL, 

mode,0,1

SFL,4,PRES,1e7/30    

solv

pi=3.1415926

FU=195     !傅里葉展開的項數

*do,i,1,FU

mode,i,1

zaihe=2e7/i/pi*sin(i*pi/30)*cos(i*0)

SFL,4,PRES,zaihe  

solv

*enddo

FINISH  

!載荷工況最多可設置99個，將總載荷步數按99一組進行分組

!第一組99個

/post1

SET,1,,1,,,0  

LCDEF,1,1,,    !設定工況1=載荷步1，第一個1為工況編號，第二個1為載荷步編號

LCASE,1        !讀入工況1，database＝1

*do,i,2,99 

SET,i,,1,,,0    

LCDEF,i,i,, 

LCASE,i 

*enddo          !依次將各個工況讀入，database＝i

!荷載組合，組合方式為工況i+工況1+工況2+...+工況i-1

*do,j,1,98      !荷載組合，database＝database＋1

LCOPER,ADD,j,,,   

*enddo  

plnsol,s,EQV    !查看該組合工況下的結構mises應力

!第二組96個，在第二組的基礎上直接疊加

*do,n,99+1,FU  

*do,m,1,FU-99

SET,n,,1,,,0

LCDEF,m,n,, 

LCASE,m 

*enddo

*enddo          

*do,j,1,FU-99     !荷載組合，database＝database＋1

LCOPER,ADD,j,,,   

*enddo  

plnsol,s,EQV    !查看該組合工況下的結構mises應力