非轴对称荷载下轴对称问题的有限元分析 - 道客巴巴 (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应力