24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2948  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工学 309分求调剂 +7 Yin DY 2026-04-08 7/350 2026-04-13 22:50 by pies112
[考研] 279学硕食品专业求调剂院校 20+6 孤独的狼爱吃羊 2026-04-12 23/1150 2026-04-13 20:42 by biomen
[考研] 290求调剂 +16 luoziheng 2026-04-10 16/800 2026-04-13 20:34 by biomen
[考研] 材料299专硕求调剂 +13 +21 2026-04-09 13/650 2026-04-13 14:16 by 张zhihao
[考研] 366求调剂 +9 不知名的小卅 2026-04-11 9/450 2026-04-13 01:19 by 幸免 ..
[考研] 085410 273求调剂 +10 X1999 2026-04-09 10/500 2026-04-12 09:24 by 逆水乘风
[考研] 085404 293求调剂 +9 勇远库爱314 2026-04-08 9/450 2026-04-12 02:24 by 秋豆菜芽
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 一志愿东北大学控制工程085406数二英二385,求调剂 +8 Ezra_Zhang 2026-04-09 8/400 2026-04-11 09:15 by 猪会飞
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 材料专业344求调剂 +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 367求调剂 +10 hffQAQ 2026-04-09 10/500 2026-04-09 18:06 by lijunpoly
[考研] 机械工程313分找工科调剂 +3 双一流本科机械 2026-04-08 3/150 2026-04-08 20:41 by 土木硕士招生
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +10 努力奋斗112 2026-04-07 10/500 2026-04-08 15:01 by screening
[考研] 计算机11408 287 求调剂 +3 LiLe5 2026-04-07 3/150 2026-04-07 23:15 by shanqishi
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
信息提示
请填处理意见