24小时热门版块排行榜

返回列表

823978808

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 101
帖子: 30
在线: 147.1小时
虫号: 2506546
注册: 2013-06-13
专业: 机械摩擦学与表面技术

[求助] 求助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.  看到论坛有好多温度场模拟，请高手指教，望不吝给出命令流进行参考！

回复此楼