24小时热门版块排行榜    

查看: 3020  |  回复: 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.  看到论坛有好多温度场模拟,请高手指教,望不吝给出命令流进行参考!
回复此楼

» 收录本帖的淘帖专辑推荐

ANSYS命令流

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

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

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 823978808 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 范进中举一文的中心思想 +4 炎黄贵胄 2026-08-22 5/250 2026-08-24 02:09 by 流流伤
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见