*SET,tshell,8 !
*SET,ltb,1250 !
*SET,dshell,447 !
*SET,gbtk,50 !
*SET,d2,19 !
*SET,sp,32 !
*SET,t2,2 !
*SET,ny,6 !
*SET,JD,90
*DIM,nxz,,ny !
*SET,nxz(1),7,7,6,6,6,5
*SET,ps,1.2 !
*SET,pt,0.96 !
*SET,ttube,5 !
*SET,pi,3.1415926
*SET,nr0,60
*SET,nr1,3
*SET,nb1,2 !
*SET,nb2,6 !
*SET,nb3,2 !
*SET,ng1,30 !
*SET,ng2,1 !
*SET,s1,sp
*SET,ROL2,504/2
! ***************物性参数***************
/prep7
et,1,shell57 !定义辅助建模的壳单元类型
et,2,solid185 !进行热分析的体单元类型
!et,1,shell57 !
!et,2,solid70 !
*SET,rol,dshell/2
MPTEMP,,,,,,,,
MPTEMP,1,-45
MPTEMP,2,0
MPTEMP,3,20
MPTEMP,4,50
MPTEMP,5,200
MPTEMP,6,220
MPDATA,EX,1,,1.95E5
MPDATA,EX,1,,1.95E5
MPDATA,EX,1,,1.95E5
MPDATA,EX,1,,1.935E5
MPDATA,EX,1,,1.84E5
MPDATA,EX,1,,1.828E+005
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,-45
MPTEMP,2,0
MPTEMP,3,20
MPTEMP,4,50
MPTEMP,5,200
MPTEMP,6,220
UIMP,1,REFT,,,20
MPDATA,ALPX,1,,16.001E-6
MPDATA,ALPX,1,,16.28E-6
MPDATA,ALPX,1,,16.384E-6
MPDATA,ALPX,1,,16.54E-6
MPDATA,ALPX,1,,17.25E-6
MPDATA,ALPX,1,,17.318E-6
MPTEMP,1,-45
MPTEMP,2,0
MPTEMP,3,20
MPTEMP,4,100
MPTEMP,5,150
MPTEMP,6,200
MPTEMP,7,250
MPDATA,KXX,1,,16.3
MPDATA,KXX,1,,16.3
MPDATA,KXX,1,,16.3
MPDATA,KXX,1,,16.3
MPDATA,KXX,1,,17
MPDATA,KXX,1,,17.3
MPDATA,KXX,1,,18.5
pcirc,,rol,,JD
pcirc,ROL,rol2,,JD
*SET,ntol,0 !
*do,j,1,ny !
*SET,ntemp,j
*do,i,1,nxz(ntemp)
*if,ntemp,ne,nint(ntemp/2)*2,then !
wpoffs,sp*(i-1),sp*(j-1) !行数j为奇数时该行第i
*else
wpoffs,sp*(i-1),sp*(j-1) !行数j为偶数时该行第i
*endif
pcirc,,d2/2.0d0 !
*SET,ntol,ntol+1 !计数器加1
wpcsys,1,0 !
*enddo
*enddo
asel,u,area,,1 !
asel,u,area,,2
cm,base1,area !将上一行选中的所有以换热管外直径为直径的圆平面命名为base1
allsel !
asba,1,base1 !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!换热管下表面,
*do,j,1,ny !
*SET,ntempj,j
*do,i,1,nxz(ntempj)
*SET,ntempi,i
*SET,x2,(i-1)*sp
*SET,y2,(j-1)*sp
wpoffs,x2,y2
*if,ntempj,eq,1,and,ntempi,eq,1,then
pcirc,d2/2.0d0-t2,d2/2.0d0,,90
*else
*if,ntempj,eq,1,then
pcirc,d2/2.0d0-t2,d2/2.0d0,,90
pcirc,d2/2.0d0-t2,d2/2.0d0,90,180
*else
*if,ntempi,eq,1,then
pcirc,d2/2.0d0-t2,d2/2.0d0,270,360
pcirc,d2/2.0d0-t2,d2/2.0d0,,90
*else
pcirc,d2/2.0d0-t2,d2/2.0d0
*endif
*endif
*endif
wpcsys,1,0
*enddo
*enddo
lsel,s,radius,,rol !
lsel,A,radius,,rol2
lesize,all,,,nr0 !
allsel !
lsel,s,radius,,d2/2,d2/2-t2
lesize,all,,,nr1 !
allsel !
lsel,s,length,,t2 !
lesize,all,,,1 !换热管壁厚方向剖分1
allsel !
mshape,0,2d !
type,1
!esize,2 !
aglue,all !
numcmp,area !
numcmp,line !
lsel,s,radius,,rol !
asll,s,0 !
cm,gb1,area !给布管区管板的下表面命名为basplt
allsel !
asel,u,area,,gb1 !
cm,hrg1,area !给换热管下表面的圆环面命名为bastub
allsel
esize,2
mat,1
/AUTO,1
/REP,FAST
FLST,5,2,4,ORDE,2
FITEM,5,4
FITEM,5,6
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,5, , , , ,1
!*
FLST,5,12,4,ORDE,3
FITEM,5,286
FITEM,5,288
FITEM,5,-298
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,3, , , , ,1
!*
FLST,5,1,4,ORDE,1
FITEM,5,287
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,8, , , , ,1
!*
lsel,s,radius,,d2/2,d2/2-t2
lesize,all,,,nr1
ALLSEL,ALL
amesh,gb1
k,10001,0,0,0 !
k,10002,0,0,ttube !
k,10003,0,0,gbtk-tshell !
k,10004,0,0,gbtk !
numcmp,line !
l,10001,10002,nb1 !
l,10002,10003,nb2 !
l,10003,10004,nb3 !
*get,line0,line,,num,max!
type,2 !
vdrag,gb1,,,,,,line0-2,line0-1,line0 !
aclear,gb1 !
mshape,0,2d !
type,1 !
!mat,1 !
amesh,hrg1 !
k,10005,0,0,gbtk+ltb !
l,10004,10005 !
*get,line1,line,,num,max!
lesize,line1,,,ng1,ng2 !
type,2 !
vdrag,hrg1,,,,,,line0-2,line0-1,line0 ,line1 !
aclear,hrg1 !
FINISH
/SOL
lsel,s,radius,,d2/2
asll
asel,r,loc,z,0
aplot
cm,cma1,area
lsel,s,radius,,d2/2-t2
asll
asel,u,loc,z,0
asel,u,loc,z,ttube
asel,u,loc,z,gbtk-tshell
asel,u,loc,z,gbtk
asel,u,loc,z,gbtk+ltb
cm,cma2,area
LSEL,S,RADIUS,,D2/2
ASLL
asel,u,loc,z,0
asel,u,loc,z,0,gbtk
asel,u,loc,z,gbtk+ltb
CM,CMA3,AREA
LSEL,S,RADIUS,,ROL
Lsel,R,loc,z,gbtk
ASLL
ASEL,R,LOC,Z,GBTK
CM,CMA4,AREA
LSEL,S,RADIUS,,ROL2
Lsel,R,loc,z,gbtk
ASLL
CM,CMA5,AREA
CMSEL,S,CMA4
CMSEL,U,CMA5
CM,CMA4,AREA
cmsel,s,cma5
asel,r,loc,z,gbtk
cm,cma5,area
asel,s,loc,z,ltb+gbtk
cm,cma6,area
cmsel,s,cma1
cmsel,a,cma2
cm,cma1,area
EPLOT
/USER, 1
/VIEW, 1, 0.627452583991E-01, 0.627528599713E-01, 0.996054773149
/ANG, 1, 0.108530502480
/REPLO
/VIEW, 1, 0.672211728368E-01, 0.104293090943 , -0.992272273675
/ANG, 1, 26.6050182769
/REPLO
/VIEW, 1, -0.146419155306 , -0.174414320817 , -0.973725359459
/ANG, 1, 26.6923257035
/REPLO
/sol !!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*SET,bc1,8541.71/1000000
*SET,bc2,902.15/1000000
!!
allsel
DADELE,all,TEMP
cmsel,s,cma1
SFA,all,1,CONV,bc1,220
cmsel,s,cma3
SFA,all,1,CONV,bc2,20
cmsel,s,cma4
SFA,all,1,CONV,bc2,20
alls
LSWRITE,1,
!!
DADELE,all,TEMP
cmsel,s,cma1
SFA,all,1,CONV,bc1,220
cmsel,s,cma3
SFA,all,1,CONV,bc2,-45
cmsel,s,cma4
SFA,all,1,CONV,bc2,-45
alls
LSWRITE,2,
!!
DADELE,all,TEMP
cmsel,s,cma1
SFA,all,1,CONV,bc1,80
cmsel,s,cma3
SFA,all,1,CONV,bc2,20
cmsel,s,cma4
SFA,all,1,CONV,bc2,20
alls
LSWRITE,3,
!!
DADELE,all,TEMP
cmsel,s,cma1
SFA,all,1,CONV,bc1,210
cmsel,s,cma3
SFA,all,1,CONV,bc2,20
cmsel,s,cma4
SFA,all,1,CONV,bc2,20
alls
LSWRITE,4,
!!
DADELE,all,TEMP
cmsel,s,cma1
SFA,all,1,CONV,bc1,210
cmsel,s,cma3
SFA,all,1,CONV,bc2,-36
cmsel,s,cma4
SFA,all,1,CONV,bc2,-36
alls
LSWRITE,5,
LSSOLVE,1,5,1, |