轴对称热应力中如何获得周向应力
我做两个不同材料的圆柱形结合在一起,采用轴对称建模(XY平面内),但是结果中无法获得沿圆周方向的应力?程序如下:
/units,si !取用国际单位制
Rad=20e-3
H=1e-3
L=1e-3+0.2e-3
/prep7
et,1,plane13 !定义单元类型
KEYOPT,1,1,4 !设置单元基本选项,即设置options对话框
KEYOPT,1,3,1 !输出变量控制
MP,EX,1,1050e9 !设置材料属性
MP,PRXY,1,0.2
MP,DENS,1,3e3
mp,c,1,1.47e-3
mp,kxx,1,2000
mptemp,1,300,500,1173
mpdata,alpx,1,,1e-6,2.7e-6,4.4e-6
MP,PRXY,2,0.324
mp,c,2,272.35
MP,DENS,2,10.2e3
MP,ex,2,315.882e9
MP,KXX,2,1470
MPTEMP,1,300,500,1173
MPDATA,alpx,2,,5.8e-6,6.0e-6,6.2e-6
rectng,0,Rad,0,H!建立几何模型
rectng,0,Rad,H,L
AGLUE,all
NUMCMP,area
/PNUM,area,1
/PNUM,line,1
aplot
asel,s,area,,1!分配单元属性
aatt,1,1,1
asel,s,area,,2
aatt,2,1,1
allsel
/PNUM,MAT,1
aplot
AESIZE,all,100e-6, !网格化(CVD)
MSHAPE,0,2D
MSHKEY,0
amesh,all
ALLSEL,ALL
tref,0 !加载
ic,all,temp,1173
***CTRAN
nsel,s,loc,y,0
d,all,uy,0
sf,all,conv,120,300
ALLSEL,ALL
lsel,s,loc,x,0
dl,all,,symm !对称加载
allsel
***CTRAN !将实体载荷转化为有限元模型载荷
allsel,all
save
finish
/solu
antype,4
TIMINT,ON
time,120
deltim,2,1,5,1
kbc,1
autots,1
outres,ESOL
save
solve
finish
/post1
set,last
PLNSOL,TEMP, ,1,
plnsol,s,x,1,
plnsol,s,y,1,
plnsol,s,z,1,
finish
请各位大哥和斑竹多多帮忙!急!!!
Re:轴对称热应力中如何获得周向应力
有一事不明:***CTRAN
请楼主解释一下
Re:轴对称热应力中如何获得周向应力
***CTRAN是我的错误,它因改为:***CTRAN还有在AGLUE,al后l对面积编号有重新整理:NUMCMP,area
编号为3的面积转为编号为2的面积!
我现在的疑惑的是结构中沿圆周方向的热应力基本与径向一致,这和理论分析
不符合!
能不能帮我分析一下如何产生正确的圆周方向的热应力
Re:轴对称热应力中如何获得周向应力
我怎么觉得模型简化不是很清晰其中的约束怎么看也觉得不明白,对称约束是否加错地方了?:
***CTRAN
nsel,s,loc,y,0
d,all,uy,0
sf,all,conv,120,300
ALLSEL,ALL
lsel,s,loc,x,0
dl,all,,symm !对称加载
allsel
***CTRAN !将实体载荷转化为有限元模型载荷
我觉得对称约束应该加在y=0的位置,不知楼主觉得如何?
Re:轴对称热应力中如何获得周向应力
我对此也不是很肯定,请版主帮忙确定一下,对称约束加载在那边?我也曾经找过相关的对称加载,但是总是不明确!
希望版主能帮忙弄清楚!
页:
[1]