24小时热门版块排行榜    

查看: 2941  |  回复: 0

823978808

新虫 (初入文坛)

[求助] 求助ansys温度场,附命令流

本人需要研究制动器摩擦温度场,思想如下:
将制动时间等分成N等分,在每个时间段内,对模型一部分施加热流密度,对其他面施加对流换热,时间推移,移动施加热流密度位置,进行多载荷步求解,得到全时间段温度结果。
以下代码是论证这个思想的实验,有以下简化:
1. 对流换热没有对所有非接触其他表面,只是选取部分,建成COMPONENT
2. 没有细究时间段内热流密度施加的位置,只是简单移动到相邻位置
3. 为简化计算,总时间为10s,时间间隔选1s
APDL如下:


/prep7
et,1,solid70
mp,ex,1,175e3
mp,prxy,1,0.3
mp,dens,1,7.228e-9
mp,alpx,1,11e-6               
uimp,1,reft,,,20               
mp,kxx,1,48
mp,c,1,480e6                  
!------------------------
rb=1000/2                              
rs=148/2                              
b_disc=30                              
x_theo=435                              
length_pan=255/2                     
width_pan=120/2                        
b_pan=35                              
numcirc=18                        
rfmin=x_theo-width_pan                 
rfmax=x_theo+width_pan                  
numrad=12                              
incrad=(rfmax-rfmin)/numrad         
numrsmesh_b=6                     
numrfminmesh=5                          
numradmesh=20   
*do,i,1,numcirc,1
cylind,rfmin,rs,0,b_disc,360*i/numcirc,360*(i+1)/numcirc
*enddo
!*do,i,1,numcirc,1
!*do,j,0,numrad-1,1
!cylind,rfmin+j*incrad,rfmin+(j+1)*incrad,0,b_disc,360*i/numcirc,360*(i+1)/numcirc
!*enddo
!*enddo
*do,i,1,numcirc,1
cylind,rfmin,rfmax,0,b_disc,360*i/numcirc,360*(i+1)/numcirc
*enddo
*do,i,1,numcirc,1
cylind,rfmax,rb,0,b_disc,360*i/numcirc,360*(i+1)/numcirc
*enddo
allsel,all
vglue,all
numcmp,area
FLST,5,18,4,ORDE,18
FITEM,5,5   
FITEM,5,655
FITEM,5,659
FITEM,5,663
FITEM,5,667
FITEM,5,671
FITEM,5,675
FITEM,5,679
FITEM,5,683
FITEM,5,687
FITEM,5,691
FITEM,5,695
FITEM,5,699
FITEM,5,703
FITEM,5,707
FITEM,5,711
FITEM,5,715
FITEM,5,719
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,6, , , , ,1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,8   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,10, , , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,726
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,5, , , , ,1  
FLST,5,1,4,ORDE,1   
FITEM,5,792
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,1, , , , ,1  
!*  
MSHAPE,1,3D
MSHKEY,0
!*  
FLST,5,54,6,ORDE,3  
FITEM,5,1   
FITEM,5,55  
FITEM,5,-107
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
!------------网分------
asel,s,,,2,4,2
asel,a,,,30,96,2
cm,area_conv,area
asel,s,,,104,164,4
asel,a,,,26,167,141
cm,fzone_conv,area
!----------
/sol
tunif,20,
/sol
w0=1780
zt=10
dt=1
arad=w0/zt
u=0.3
p=1.8
r=rfmax
conve=10        !对流换热系数
convet=20       !参考温度
subst1=10
numloop=18
*dim,t,array,10,,,,,
*do,i,1,10
t(i)=dt*i
*enddo
*dim,w,array,10,,,,,,
*do,i,1,10
w(i)=w0-arad*i*dt
*enddo
*dim,q,array,10,,,,
*do,i,1,10
q(i)=u*p*w(i)*r
*enddo
!------------
/sol
timint,on,therm
tref,20

sfa,area_conv,1,conv,conve,convet   !这句存疑,sfa命令能够在多载荷步中一直存在吗?
!-------
sfa,fzone_conv,1,conv,conve,convet
sfa,26,1,hflux,q(1)
time,dt
kbc,1
nsubst,subst1
antype,4
trnopt,full
toffset,273,
outpr,all,all
solve
sfadele,area_conv,1,conv,conve,convet
!-------------------
*do,i,1,8
sfa,fzone_conv,1,conv,conve,convet
sfa,104,1,hflux,q(i+1)
j=104+4*i
time,dt*(i+1)
kbc,1
nsubst,subst1
antype,4
trnopt,full
toffset,273,
outpr,all,all
solve
sfadele,area_conv,1,conv,conve,convet
*enddo
!--------由于加载的表面由多个部分粘结而成,
!--------编号不连续或呈等差关系,
!--------所以不能只用一次循环
sfa,fzone_conv,1,conv,conve,convet
sfa,167,1,hflux,q(10)
time,dt*10
kbc,1
nsubst,subst1
antype,4
trnopt,full
toffset,273,
outpr,all,all
solve
sfadele,area_conv,1,conv,conve,convet
finish                       
有以下问题:
1.  能够运算,但是结果显示好像施加热流密度的位置没有移动
2.  温度出现负值
3.  看到论坛有好多温度场模拟,请高手指教,望不吝给出命令流进行参考!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 823978808 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学求调剂 +4 奶油草莓. 2026-03-22 5/250 2026-03-23 19:37 by pswait
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[考研] 一志愿070300浙大化学358分,求调剂! +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
信息提示
请填处理意见