求助：关于轴对称模型施加非轴对称载荷分析的问题？ 50
非轴对称荷载下轴对称问题的有限元分析 - 道客巴巴 (doc88.com)
最近做了个例子，关于轴对称结构受非轴对称载荷情况下采用二维模型进行计算的，文中说对每个截面的应力分布均可以求得，但我实际计算只求得0度截面的应力分布同三维模型结果接近，45度的计算结果差距较大，想请教大家，如何求得任意角度截面的应力分布？其原理又是啥？
给出模型及求解及后处理的命令流
!静力计算，直线工况，只加垂向力
/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应力
