24小时热门版块排行榜    

查看: 1065  |  回复: 0

殷娜的宽恕

新虫 (初入文坛)

[求助] ansys瞬态热分析热膨胀变化速度异常

ansys瞬态热分析进行一个铝板传热的仿真。
一开始热膨胀变化很快,到后面热膨胀变化变得很慢,一直达不到热平衡状态。仿真时间达到90000s以上,热变形还在继续。导致无法得到预期结果。
下面是命令流
/prep7
lengh1=0.050
lengh2=0.050
lengh3=0.050
lengh4=0.050
lengh5=0.050
lengh6=0.050
lengh7=0.050
lengh8=0.050
lengh9=0.050
lengh10=0.050
llengh=0.5
!带宽分段水冷

glengh1=0.040
glengh2=0.260
glengh=0.8
!辊身其余部分
blengh1=0.256
blengh2=0.144!轴承作用范围
blengh=0.4
r0=0.09750
r1=0.090
r2=0.040
r3=0.325

ed1=0.005        !网格分段1
ed2=0.005        !网格分段2
ed3=0.005        !网格分段3
ed4=0.005        !网格划分4
ed5=0.005        !网格划分5
ed6=0.005        !网格划分6
ed7=0.005
edr1=0.005
edr2=0.002         

heatdens=5.1168e4*1.04       !热流密度热源
heatcon_bear=1.0111e4*0.03!轴承热流密度热源
con_water=1725          !水冷散热系数
con_water1=1725         !强水冷散热系数        
con_water2=250          !弱水冷散热系数
con_air=40             !空冷散热系数
con_bear=250
tem_water=38            !冷却液环境温度
tem_air=26              !空冷环境温度
tem_workroll=26         !工件温度

et,1,13,4,,1,   
et,2,151,,,1,
KEYOPT,2,8,2
et,3,151,,,1,
KEYOPT,3,8,1
MP,EX,1,2.36e11   
MP,PRXY,1,0.3  
MP,kxx,1,44   
MP,dens,1,7810
MP,c,1,460
UIMP,1,REFT,,,20
MP,ALPX,1,1.15e-5  



!参数化建模
RECTNG,,r0,,lengh1,
RECTNG,,r0,lengh1,lengh1+lengh2,
RECTNG,,r0,lengh1+lengh2,lengh1+lengh2+lengh3,
RECTNG,,r0,lengh1+lengh2+lengh3,lengh1+lengh2+lengh3+lengh4,
RECTNG,,r0,lengh1+lengh2+lengh3+lengh4,lengh1+lengh2+lengh3+lengh4+lengh5,
RECTNG,,r0,lengh1+lengh2+lengh3+lengh4+lengh5,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6,
RECTNG,,r0,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6, lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7
RECTNG,,r0,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7, lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8
RECTNG,,r0,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8, lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9
RECTNG,,r0,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9, lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9+lengh10

RECTNG,,r0,llengh,llengh+glengh1,
RECTNG,,r0,llengh+glengh1,llengh+glengh1+glengh2,

RECTNG,,r0,glengh,glengh+blengh1,
RECTNG,,r0,glengh+blengh1,glengh+blengh1+blengh2,


RECTNG,r0,r1+r2,,lengh1,
RECTNG,r0,r1+r2,lengh1,lengh1+lengh2,
RECTNG,r0,r1+r2,lengh1+lengh2,lengh1+lengh2+lengh3,
RECTNG,r0,r1+r2,lengh1+lengh2+lengh3,lengh1+lengh2+lengh3+lengh4,
RECTNG,r0,r1+r2,lengh1+lengh2+lengh3+lengh4,lengh1+lengh2+lengh3+lengh4+lengh5,
RECTNG,r0,r1+r2,lengh1+lengh2+lengh3+lengh4+lengh5,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6,
RECTNG,r0,r1+r2,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6, lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7
RECTNG,r0,r1+r2,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7, lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8
RECTNG,r0,r1+r2,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8, lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9
RECTNG,r0,r1+r2,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9, lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9+lengh10

RECTNG,r0,r1+r2,llengh,llengh+glengh1,
RECTNG,r0,r1+r2,llengh+glengh1,llengh+glengh1+glengh2,
AGLUE,all

!对应的网格分段设置
lsel,s,loc,y,lengh1*0.5
LESIZE,all,ed1
lsel,s,loc,y,lengh1+lengh2*0.5
LESIZE,all,ed2
lsel,s,loc,y,lengh1+lengh2+lengh3*0.5
LESIZE,all,ed3
lsel,s,loc,y,lengh1+lengh2+lengh3+lengh4*0.5
LESIZE,all,ed4
lsel,s,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5*0.5
LESIZE,all,ed5
lsel,s,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6*0.5
LESIZE,all,ed6
lsel,s,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7*0.5
LESIZE,all,ed7
lsel,s,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8*0.5
LESIZE,all,ed7
lsel,s,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9*0.5
LESIZE,all,ed7
lsel,s,loc,y,llengh-lengh10*0.5
LESIZE,all,ed7
lsel,s,loc,y,llengh+glengh1*0.5
LESIZE,all,ed7
lsel,s,loc,y,llengh+glengh1+glengh2*0.5
LESIZE,all,ed7
lsel,s,loc,y,glengh+blengh1*0.5
LESIZE,all,ed7
lsel,s,loc,y,glengh+blengh1+blengh2*0.5
LESIZE,all,ed7

lsel,s,loc,x,r0*0.5
LESIZE,all,edr1
lsel,s,loc,x,r0+r3*0.5
LESIZE,all,edr2

MSHAPE,0,2D
MSHKEY,0
AMESH,all

!以下为施加热部分
!定义载荷的数组
*DIM,heatload,ARRAY,5,1,1, , ,  
*SET,HEATLOAD(1,1,1) , heatdens       ! 热流密度
*SET,HEATLOAD(2,1,1) , con_water      !水冷
*SET,HEATLOAD(3,1,1) , con_water1     !强水冷
*SET,HEATLOAD(4,1,1) , con_water2     !弱水冷
*SET,HEATLOAD(5,1,1) , con_air        !空冷



nsel,s,loc,x,r1+r2
TYPE, 2   
MAT,  1
ESURF,0

TYPE, 3   
MAT,  1
ESURF,0
nsel,s,loc,x,r0
nsel,r,loc,y,glengh,glengh+blengh
TYPE, 3  
MAT,  1
ESURF,0
!surf单元预处理


nsel,s,loc,x,r1+r2
nsel,r,loc,y,0,llengh
ESEL,S,TYPE,,3
ESLN,R,1
SFE,all,,HFLUX, ,HEATLOAD(1,1,1), , ,      ! 热流密度   
allsel
!轧辊与板带施加热部分
nsel,s,loc,x,r0
nsel,r,loc,y,glengh+blengh1,glengh+blengh
ESEL,S,TYPE,,3
ESLN,R,1
SFE,all,,HFLUX, ,heatcon_bear, , ,      ! 热流密度1
allsel
!轧辊与轴承施加热部分


nsel,s,loc,x,r1+r2
nsel,r,loc,y,0,lengh1
ESEL,S,TYPE,,2
ESLN,R,1  
SFE,all,1,CONV, ,HEATLOAD(3,1,1)       !1段水冷     改变位置 1   
SFE,all,1,CONV,2,tem_water
allsel
nsel,s,loc,x,r1+r2
nsel,r,loc,y,lengh1,lengh1+lengh2  
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(2,1,1)      !2 段水冷   
SFE,all,1,CONV,2,tem_water      
allsel

nsel,s,loc,x,r1+r2
nsel,r,loc,y,lengh1+lengh2,lengh1+lengh2+lengh3
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(2,1,1)         ! 3段水冷     
SFE,all,1,CONV,2,tem_water                 
allsel

nsel,s,loc,x,r1+r2
nsel,r,loc,y,lengh1+lengh2+lengh3,lengh1+lengh2+lengh3+lengh4
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(2,1,1)            ! 4段水冷   
SFE,all,1,CONV,2,tem_water
allsel

nsel,s,loc,x,r1+r2
nsel,r,loc,y,lengh1+lengh2+lengh3+lengh4,lengh1+lengh2+lengh3+lengh4+lengh5
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(2,1,1)          ! 5段水冷   
SFE,all,1,CONV,2,tem_water
allsel


nsel,s,loc,x,r1+r2
nsel,r,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(2,1,1)          ! 6段水冷   
SFE,all,1,CONV,2,tem_water
allsel

nsel,s,loc,x,r1+r2
nsel,r,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(2,1,1)            !7段水冷   
SFE,all,1,CONV,2,tem_water
allsel
nsel,s,loc,x,r1+r2
nsel,r,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(2,1,1)          !8段水冷   
SFE,all,1,CONV,2,tem_water
allsel
nsel,s,loc,x,r1+r2
nsel,r,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(2,1,1)            ! 9段水冷   
SFE,all,1,CONV,2,tem_water
allsel

nsel,s,loc,x,r1+r2
nsel,r,loc,y,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9,lengh1+lengh2+lengh3+lengh4+lengh5+lengh6+lengh7+lengh8+lengh9+lengh10
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(2,1,1)            ! 10段水冷   
SFE,all,1,CONV,2,tem_water
allsel


nsel,s,loc,x,r1+r2
nsel,r,loc,y,llengh,llengh+glengh1
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(5,1,1)            ! 5 段空冷  
SFE,all,1,CONV,2,tem_air
allsel

nsel,s,loc,x,r1+r2
nsel,r,loc,y,llengh+glengh1,glengh
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(5,1,1)            ! 6段空冷
SFE,all,1,CONV,2,tem_air
allsel
nsel,s,loc,y,glengh
nsel,r,loc,x,r0,r1+r2
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(5,1,1)            ! 辊肩空冷
SFE,all,1,CONV,2,tem_air
allsel

nsel,s,loc,x,r0
nsel,r,loc,y,glengh,glengh+blengh1
ESEL,S,TYPE,,2
ESLN,R,1
SFE,all,1,CONV, ,HEATLOAD(5,1,1)            ! 辊颈空冷
SFE,all,1,CONV,2,tem_air
allsel

!约束
nsel,s,loc,y,glengh+blengh
D,all,UY,0
allsel

nsel,s,loc,y,0
D,all,UY,0
allsel


nsel,s,loc,x,r0,r1+r2
nsel,r,loc,y,glengh
D,all,UY,0
allsel

FINISH  
/SOL
  
ANTYPE,TRANS               
KBC,1                     
OUTRES,,ALL               
AUTOTS,ON                  
DELTIM,50,50,100         
TIME,36000   

TUNIF,26  
            

solve

FINISH
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 殷娜的宽恕 的主题更新
信息提示
请填处理意见