24小时热门版块排行榜    

查看: 2992  |  回复: 21

微言1021

木虫 (小有名气)

[求助] ansys瞬态热分析求助

本人在做塑料管焊接的ansys瞬态热分析,命令流如下:
*set,d1,0.0654
*set,d2,0.08
*set,l1,0.005
*set,l2,0.1
/prep7
et,1,plane55
keyopt,1,3,1
mptemp,1,-20,100,125,150,200,250
mpdata,dens,1,1,945,950,955,960,965
mpdata,kxx,1,1,0.341,0.3,0.265,0.264,0.263
mpdata,enth,1,0,2e5,3.18e5,5.35e5,7e5,9e5
mpdata,c,1,1,2.87e3,5.56e3,2.65e3,2.67e3,3.03e3
rect,d1,d2,0,l1
rect,d1,d2,l1,l2
aglue,all
lplot,all
lsel,s,loc,y,0
lsel,a,loc,y,l1
lsel,a,loc,y,l2
lesize,all,,,20
lsel,s,length,,l1
lesize,all,,,10
lsel,s,length,,l2-l1
lesize,all,,,50
lsel,all
aplot
aatt,1,,1
mshape,0,2d
mshkey,1
amesh,all
/solu            
tunif,-20
nsel,s,loc,x,d2
sf,all,conv,37.24,-20
nsel,s,loc,x,d1
sf,all,conv,0.72,-20
nsel,s,loc,y,l2
sf,all,conv,0,-20
lsel,s,loc,y,0
nsll,s,all
d,all,temp,210
alls,all
eplot,all

antype,trans
eqslv,iter
solcontrol,on
lnsrch,on
kbc,1
timint,on
time,20                  
nsubst,20
autots,on
tintpr,,,,1.0
lsel,s,loc,y,0
nsll,s,all
d,all,temp,210
alls,all
outres,all,all
solve
time,160                 
timint,on
nsubst,140
autots,on
lsel,s,loc,y,0
nsll,s,all
d,all,temp,210
alls,all
outres,all,all
solve
time,164                 
timint,on
nsubst,14
autots,on
nsel,all
nsel,s,loc,y,0
ddele,all,temp
sf,all,conv,0.5,210
eplot,all
alls,all
outres,all,all
solve
time,644                     
timint,on
nsubst,700
autots,on
nplot,all
nsel,s,loc,y,0
sf,all,conv,0,150
eplot,all
alls,all
outres,all,all
solve
time,1244                     
timint,on
nsubst,600
autots,on
nplot,all
nsel,s,loc,y,0
sf,all,conv,0,-20
alls,all
outres,all,all
solve
做出的结果不对,在第一、第二载荷步中施加的温度有210℃,在第三载荷步删除温度载荷,第三步只有4s的时间,到第四载荷步的时候温度就下降到了-18℃,我觉得不符合实际情况,请各位高手解答一下是哪里错了?导致温度下降这么迅速。我设置的平均温度为-20℃。谢谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

言之命至,心随己愿
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
微言1021: 金币+5, ★★★很有帮助 2013-04-14 22:20:24
xiegangmai: 金币+3, 谢谢参与! 2013-04-15 18:43:09
你的问题在于第三步时,删除了底部的固定温度 210 度,代以对流换热
        sf,all,conv,0.5,210
    但是,其中的换热系数 0.5,太小了,所以导致你说的结果。假如将换热系数加大,比如改为 50,温度下降就不那么快了。
    分析了一下你的命令流,做了一些注解,你可以看一下:


fini
/clear
/filn,weiyan-1021

*set,d1,0.0654
*set,d2,0.08
*set,l1,0.005
*set,l2,0.1
!*
/prep7
et,1,plane55
keyopt,1,3,1
mptemp,1,-20,100,125,150,200,250
mpdata,dens,1,1,945,950,955,960,965
mpdata,kxx,1,1,0.341,0.3,0.265,0.264,0.263
mpdata,enth,1,0,2e5,3.18e5,5.35e5,7e5,9e5
mpdata,c,1,1,2.87e3,5.56e3,2.65e3,2.67e3,3.03e3
rect,d1,d2,0,l1
rect,d1,d2,l1,l2
aglue,all
lplot,all
lsel,s,loc,y,0
lsel,a,loc,y,l1
lsel,a,loc,y,l2
lesize,all,,,20
lsel,s,length,,l1
lesize,all,,,10
lsel,s,length,,l2-l1
lesize,all,,,50
lsel,all
aplot
aatt,1,,1
mshape,0,2d
mshkey,1
amesh,all

/solu            
tunif,-20             !  均匀温度 -20
nsel,s,loc,x,d2
sf,all,conv,37.24,-20      !  右边界对流边界条件,流体温度 -20
!*
nsel,s,loc,x,d1
sf,all,conv,0.72,-20       !  左边界对流边界条件,流体温度 -20
!*
nsel,s,loc,y,l2
sf,all,conv,0,-20          !  顶部对流边界条件,流体温度 -20
!*
lsel,s,loc,y,0
nsll,s,all
d,all,temp,210
!*
alls,all
eplot,all

antype,trans
eqslv,iter
solcontrol,on
lnsrch,on
kbc,1
timint,on
time,20                  
nsubst,20
autots,on
tintpr,,,,1.0
lsel,s,loc,y,0          !  y=0, 固定温度 210。
nsll,s,all
d,all,temp,210
alls,all
outres,all,all
solve                   !  以上第一步, 3 处对流边界,一处固定温度。
!*
time,160                 
timint,on
nsubst,140
autots,on
lsel,s,loc,y,0
nsll,s,all
d,all,temp,210
alls,all
outres,all,all
solve                    !  以上第二步,边界条件与第一步相同。
!*
time,164                 
timint,on
nsubst,14
autots,on
nsel,all
nsel,s,loc,y,0
ddele,all,temp           !  删除温度边界条件,上面三个对流边界条件仍在。
sf,all,conv,0.5,210      !  增加 y=0 处对流边界条件,对流换热系数为 0.5。
eplot,all
alls,all
outres,all,all
solve                    !  第三步, 去掉底部固定温度条件,增加底部对流边界条件,但是由于换热系数较小,作用不大,因而模型中温度很快下降,接近上面三个对流换热的流体温度。
!*
time,644                     
timint,on
nsubst,700
autots,on
nplot,all
nsel,s,loc,y,0
sf,all,conv,0,150        !  换热系数为零,这一对流换热条件不起作用,因而上面三个对流为主,使整个模型的温度与流体相同。
eplot,all
alls,all
outres,all,all
solve                    !  第四步,只有上面三个对流边界条件起作用,是整个模型温度为 -20。
!*
time,1244                     
timint,on
nsubst,600
autots,on
nplot,all
nsel,s,loc,y,0
sf,all,conv,0,-20
alls,all
outres,all,all
solve                    !  第五步,同第四步。
2楼2013-04-14 16:36:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微言1021

木虫 (小有名气)

引用回帖:
2楼: Originally posted by htbbzzg at 2013-04-14 16:36:56
你的问题在于第三步时,删除了底部的固定温度 210 度,代以对流换热
        sf,all,conv,0.5,210
    但是,其中的换热系数 0.5,太小了,所以导致你说的结果。假如将换热系数加大,比如改为 50,温度下降就不那 ...

非常感谢,对流换热系数小说明传热量小,在这里忽略不计了是吧,我以为传热系数小就应该是降温慢呢。
言之命至,心随己愿
3楼2013-04-14 22:22:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 谢谢参与! 2013-04-15 18:43:16
底面的对流 (介质温度 210 度) 应该是用来使模型下面能有一定的温度,它的换热系数小,应而上面三个对流起主要作用,使整个模型的温度向该三个对流边界的介质温度 -20 靠近了。
4楼2013-04-15 08:56:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微言1021

木虫 (小有名气)

引用回帖:
4楼: Originally posted by htbbzzg at 2013-04-15 08:56:37
底面的对流 (介质温度 210 度) 应该是用来使模型下面能有一定的温度,它的换热系数小,应而上面三个对流起主要作用,使整个模型的温度向该三个对流边界的介质温度 -20 靠近了。

那如果使底面和上面绝热的话,该怎么做呢?那个换热系数是我算出来的,应该不能改。
言之命至,心随己愿
5楼2013-04-15 09:24:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微言1021

木虫 (小有名气)

引用回帖:
4楼: Originally posted by htbbzzg at 2013-04-15 08:56:37
底面的对流 (介质温度 210 度) 应该是用来使模型下面能有一定的温度,它的换热系数小,应而上面三个对流起主要作用,使整个模型的温度向该三个对流边界的介质温度 -20 靠近了。

而且,我发现用同样的命令流,只要把尺寸改大一点,就能得到我想要的结果,这是为什么呢?
言之命至,心随己愿
6楼2013-04-15 09:35:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 谢谢参与! 2013-04-15 18:43:23
不施加任何边界条件,就是绝热边界了。
不过,如果底面和顶面都是绝热的话,结果很可能是整个模型在 -20 度左右,因为只有对流散热,而没有热量输入,模型温度会逐渐和介质温度一样。
需要检查一下换热系数的计算是否有问题。
7楼2013-04-15 09:37:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微言1021

木虫 (小有名气)

引用回帖:
7楼: Originally posted by htbbzzg at 2013-04-15 09:37:43
不施加任何边界条件,就是绝热边界了。
不过,如果底面和顶面都是绝热的话,结果很可能是整个模型在 -20 度左右,因为只有对流散热,而没有热量输入,模型温度会逐渐和介质温度一样。
需要检查一下换热系数的计算 ...

我上面的命令流是做低温大风环境下的温度分布,现在我又改成在常温下的,换热系数小的边界我觉得可以类似于绝热。
我发现只要把尺寸改大了,结果就正确了。这个是为什么呢?把命令流发给你,你看看:
/filname,true
/title,pipeanalysis
*set,d1,0.654
*set,d2,0.8
*set,l1,0.05
*set,l2,1
/prep7
et,1,plane55
keyopt,1,3,1
mptemp,1,20,100,125,150,200,250
mpdata,dens,1,1,945,950,955,960,965
mpdata,kxx,1,0.445,0.341,0.3,0.265,0.264,0.263
mpdata,enth,1,0,2e5,3.18e5,5.35e5,7e5,9e5
mpdata,c,1,1,2.87e3,5.56e3,2.65e3,2.67e3,3.03e3
rect,d1,d2,0,l1
rect,d1,d2,l1,l2
aglue,all
lplot,all
lsel,s,loc,y,0
lsel,a,loc,y,l1
lsel,a,loc,y,l2
lesize,all,,,20
lsel,s,length,,l1
lesize,all,,,10
lsel,s,length,,l2-l1
lesize,all,,,50
lsel,all
aplot
aatt,1,,1
mshape,0,2d
mshkey,1
amesh,all
/solu            
antype,trans
time,0.01
nsubst,5
nropt,full,
eqslv,front
solcontrol,on
lnsrch,on
cnvtol,heat,,,,1.0e-6

tunif,20
nsel,s,loc,x,d2
sf,all,conv,23.43,20
nsel,s,loc,x,d1
sf,all,conv,1.149,20
nsel,s,loc,y,0
d,all,temp,210
alls,all
eplot,all
solve

time,20                  
timint,on
nsubst,20
autots,on
tintpr,,,,1.0
nsel,s,loc,x,d1
sf,all,conv,1.149,24
alls,all
eplot,all
kbc,1
outres,all,all
solve


time,160                 
timint,on
nsubst,140
autots,on
nsel,s,loc,x,d1
sf,all,conv,1.149,23
eplot,all
alls,all
kbc,1
outres,all,all
solve


time,164                 
timint,on
nsubst,40
autots,on
nsel,all
nsel,s,loc,y,0
ddele,all,temp
sf,all,conv,1.493,25
nsel,s,loc,x,d1
sf,all,conv,1.149,23
eplot,all
alls,all
kbc,1
outres,all,all
solve


time,644                     
timint,on
nsubst,480
autots,on
nplot,all
nsel,s,loc,y,0
sfdele,all,conv
nsel,s,loc,x,d1
sf,all,conv,1.149,21
eplot,all
alls,all
kbc,1
outres,all,all
solve


time,1244                     
timint,on
nsubst,600
autots,on
nsel,s,loc,x,d1
sf,all,conv,1.149,20
alls,all
kbc,1
outres,all,all
solve
言之命至,心随己愿
8楼2013-04-15 09:50:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

★ ★
xiegangmai: 金币+2, 谢谢参与! 2013-04-15 18:43:32
没有很多时间仔细研究你的命令流。简单看了一下,你这次的命令流和开始的命令流有很大差别,不仅是尺寸问题。
主要的,这次定义的几个对流换热边界条件,其介质温度基本是正的,而开始的是负值;此外,对流换热系数也变化较大,所以这次结果和开始的结果不同是必然的。如果这次的结果是正确的,那就说明首次的模型确实有问题。
9楼2013-04-15 18:15:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微言1021

木虫 (小有名气)

引用回帖:
9楼: Originally posted by htbbzzg at 2013-04-15 18:15:31
没有很多时间仔细研究你的命令流。简单看了一下,你这次的命令流和开始的命令流有很大差别,不仅是尺寸问题。
主要的,这次定义的几个对流换热边界条件,其介质温度基本是正的,而开始的是负值;此外,对流换热系数 ...

但是,我的实际尺寸是
*set,d1,0.0654
*set,d2,0.08
*set,l1,0.005
*set,l2,0.1
用这个尺寸算温度下降又很快了,我的问题还是没有解决。在164s时加热端温度为75℃,165s时加热端温度就变成了52℃。下降还是很快。麻烦您在帮忙分析一下,谢谢了。
言之命至,心随己愿
10楼2013-04-16 10:49:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员8JxXkC 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见