24小时热门版块排行榜    

查看: 2942  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南航材料专317分求调剂 +3 炸呀炸呀炸薯条 2026-03-23 3/150 2026-03-23 20:47 by pswait
[考研] 环境学硕288求调剂 +7 皮皮皮123456 2026-03-22 7/350 2026-03-23 20:47 by baobaoye
[考研] 306求0703调剂一志愿华中师范 +6 纸鱼ly 2026-03-21 7/350 2026-03-23 19:06 by macy2011
[考研] 生物学一志愿985,分数349求调剂 +6 zxts12 2026-03-21 9/450 2026-03-23 18:37 by macy2011
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 310求调剂 +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
信息提示
请填处理意见