24小时热门版块排行榜    

查看: 1095  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 殷娜的宽恕 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +4 Tide man 2026-08-11 4/200 2026-08-12 00:39 by liuqidi
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
信息提示
请填处理意见