我学习ansys将近两个月了,看了很多移动热源的例子,自己小试了一下,居然成功了,哈哈
/prep7
et,1,solid70 !定义单元类型
!定义材料属性
mp,kxx,1,1.31
mp,dens,1,2205
mp,c,1,766
mp,kxx,2,148
mp,c,2,712
mp,dens,2,2330
!建模
block,0,0.1524,0,0.1524,0,0.00635
block,-0.05,0.2,-0.05,0.2,-0.01,0
!划分网格
TYPE,1
MAT,1
VSEL,S,,,1
ASLV,S
LSLA,S
FLST,5,8,4,ORDE,2
FITEM,5,1
FITEM,5,-8
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,18, , , , ,1
!*
FLST,5,4,4,ORDE,2
FITEM,5,9
FITEM,5,-12
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,3, , , , ,1
vmesh,1
ALLSEL,ALL
TYPE,1
MAT,2
VSEL,S,,,2
ASLV,S
LSLA,S
FLST,5,8,4,ORDE,2
FITEM,5,13
FITEM,5,-20
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,15, , , , ,1
!*
FLST,5,4,4,ORDE,2
FITEM,5,21
FITEM,5,-24
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,2, , , , ,1
vmesh,2
allsel,all
!建立接触单元
ET,2,170
ET,3,174
KEYOPT,3,9,0
R,3,
RMORE,
RMORE,,150
RMORE,0
KEYOPT,3,1,2
! Generate the target surface
ASEL,S,,,8
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF,ALL
allsel,all
! Generate the contact surface
ASEL,S,,,1
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF,ALL
RMORE,,150
allsel,all
!接触完毕
!生成移动热源表格
*dim,flux1,table,19,19,13,X,Y,TIME
*do,k,1,13,1
*do,i,1,19,1
flux1(i,0,k)=(i-1)*0.0254/3
*enddo
*do,j,1,19,1
flux1(0,j,k)=(j-1)*0.0254/3
*enddo
*do,i,1,19,1
*do,j,1,19,1
*if,i,eq,k+3,then
*if,j,gt,3,then
*if,j,lt,17,then
flux1(i,j,k)=6210
*else
flux1(i,j,k)=0
*endif
*else
flux1(i,j,k)=0
*endif
*else
flux1(i,j,k)=0
*endif
*enddo
*enddo
flux1(0,0,k)=(k-1)*0.01725
*enddo
finish
!进入求解
/solu
ANTYPE,TRANS
TUNIF,20 !初始温度150度
vsel,s,,,2
nslv,s,1
d,all,temp,20
allsel,all
KBC,1
*SET,DELT,1
TIMINT,ON
AUTOTS,OFF
TT=0.01725
!移动热源加载
*DO,TM,1E-6,0.207,TT
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
/GO
!*
TIME,TM
AUTOTS,0
DELTIM, , , ,1
KBC,0
!*
TSRES,ERASE
/GO
!*
!*
asel,s,,,2
SFA,all,,HFLUX,%FLUX1%
/STATUS,SOLU
SOLVE
*ENDDO
FINISH |