求助：请高手帮我看看，为什么条件对称的两边温度不一样呢？

caicai1959

这是我的命令流，我做的模型是一个包含六层的三维物体，每层热系数都不同，通过热传递来交换热量。最上面那层有空气的对流和由于洒水带走的热流密度，需要进行瞬态温度场分析。
不知道为什么，左右两边条件都是一模一样，但是最后出来的温度却差了好几度，左边明显比右边低。是不是和网格划分有关？请高手指点一下，命令流写得有点啰嗦，希望看得懂~~
finish
/clear
!热分析
/filename,themal-pavement
/title,con=24.6  
/unit,si  
Ha=0.06                 !厚度
Ht=1.28+Ha
Tair=20            !空气温度
Ta=140
Tdown=30        !下层温度
kxx1=2                 ！不同温度下的导热系数
kxx2=1.9
kxx3=1.75
kxx4=1.6
kxx5=1.45
con=24.6               ！与空气的对流系数
L=19.92
p=-13700
q=-6700
/PREP7               
!定义单元类型
TOFFST,273! 设定绝对零度
ET,1,SOLID70
et,2,surf152        
DOF,TEMP
KEYOPT,2,3,0
KEYOPT,2,4,1   
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,8,4      !根据固体表面温度计算对流系数
KEYOPT,2,9,0
!设置材料属性
MPTEMP,,60,80,100,120,140 !建立温度表
MPDATA,KXX,1,,kxx1,kxx2,kxx3,kxx4,kxx5
MPDATA,C,1,,920
MPDATA,DENS,1,,2565
MPDATA,KXX,2,,1.5
MPDATA,C,2,,920
MPDATA,DENS,2,,2550  
MPDATA,KXX,3,,1.5
MPDATA,C,3,,920
MPDATA,DENS,3,,2540
MPDATA,KXX,4,,2.44  
MPDATA,C,4,,860  
MPDATA,DENS,4,,2750
MPDATA,KXX,5,,2.2   
MPDATA,C,5,,860
MPDATA,DENS,5,,2650
MPDATA,KXX,6,,1.7   
MPDATA,C,6,,860
MPDATA,DENS,6,,1870
MPDATA,HF,7,,con   !空气对流系数
Save
!创建几何模型
k,1,-6.25,0,0,,,,   
k,2,6.25,0,0,,,,
k,3,6.25,0.5,0,,,,
k,4,-6.25,0.5,0,,,,   
k,5,-6.25,0.8,0,,,,
k,6,6.25,0.8,0,,,,   
k,7,6.25,1.1,0,,,,  
k,8,-6.25,1.1,0,,,,
k,9,-6.25,1.19,0,,,,
k,10,6.25,1.19,0,,,,
k,11,6.25,1.28,0,,,,
k,12,-6.25,1.28,0,,,,   
k,13,-6.25,Ht,0,,,,   
k,14,6.25,Ht,0,,,,
FLST,2,4,3  
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
A,P51X  
FLST,2,4,3  
FITEM,2,4   
FITEM,2,3   
FITEM,2,6   
FITEM,2,5   
A,P51X  
FLST,2,4,3  
FITEM,2,5   
FITEM,2,6   
FITEM,2,7   
FITEM,2,8   
A,P51X  
FLST,2,4,3  
FITEM,2,8   
FITEM,2,7   
FITEM,2,10  
FITEM,2,9   
A,P51X  
FLST,2,4,3  
FITEM,2,9   
FITEM,2,10   
FITEM,2,11   
FITEM,2,12   
A,P51X
FLST,2,4,3  
FITEM,2,12   
FITEM,2,11   
FITEM,2,14   
FITEM,2,13   
A,P51X
FLST,2,6,5,ORDE,2   
FITEM,2,1   
FITEM,2,-6
VEXT,P51X, , ,0,0,L,,,,  
FLST,2,6,6,ORDE,2   
FITEM,2,1   
FITEM,2,-6  
VGLUE,P51X   
Save
!定义各层属性
vsel,,,,6
vatt,1,,1,0 !指定体的单元属性
vsel,,,,5
vatt,2,,1,0
vsel,,,,4
vatt,3,,1,0
vsel,,,,3
vatt,4,,1,0
vsel,,,,2
vatt,5,,1,0
vsel,,,,1
vatt,6,,1,0
allsel
!sweep划分：长80份、宽20份、厚度3、3、3、6、6、10份
FLST,5,1,4,ORDE,1   
FITEM,5,25  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,80, , , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,22  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,60, , , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,48  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,3, , , , ,1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,43  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,3, , , , ,1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,38  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,3, , , , ,1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,33  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,6, , , , ,1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,28  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,6, , , , ,1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,21  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  LESIZE,_Y1, , ,10, , , , ,1
!*  
FLST,5,6,6,ORDE,2   
FITEM,5,1   
FITEM,5,-6  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
asel,s,loc,y,Ht     
type,2         !生成表面效应单元
amesh,all
esel,s,type,,2      !选择表面效应单元
sf,all,conv,-7,Tair !确定对流边界条件
asel,s,area,,1,5,1
asel,a,area,,7,27,2
asel,a,area,,12,24,4
asel,a,area,,8
da,all,temp,Tdown   !两边和底面绝热
allsel,all
finish
！第一遍钢轮去1
/solu
Antype,trans  !设定分析类型为瞬态分析
esel,s,mat,,2,6,1
nsle,s
/replot
IC,all,TEMP,Tdown, ！施加初始下层温度
esel,s,mat,,1
nsle,s
/replot
IC,all,TEMP,Ta, ！施加初始温度
allsel
kbc,1   
timint,on !打开瞬态功能
autots,on!打开自动时间步长
tim=0
*do,i,0,L-1.66,1.66
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-6.25,-4.27  
cm,hl1,nodes
sf,hl1,hflux,p
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,4.27,6.25
cm,hl2,nodes
sf,hl2,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,36
deltim,3
solve
finish
！第一遍回1
/solu
antype,,rest
tim=36   
*do,i,0,L/1.66-1,1
tim=tim+2
time,tim
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-6.25,-4.27
cm,hl1,nodes
sf,hl1,hflux,p
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,4.27,6.25
cm,hl2,nodes
sf,hl2,hflux,p  
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,72  
deltim,3
solve
finish
！第一遍去2
/solu
antype,,rest
tim=72
*do,i,0,L-1.66,1.66
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-4.5,-2.52  
cm,hl1,nodes
sf,hl1,hflux,p
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,2.52,4.5
cm,hl2,nodes
sf,hl2,hflux,p
allsel  
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,108   
deltim,3
solve
finish
！第一遍回2
/solu
antype,,rest
tim=108   
*do,i,0,L/1.66-1,1
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-4.5,-2.52
cm,hl1,nodes
sf,hl1,hflux,p
allsel  
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,2.52,4.5
cm,hl2,nodes
sf,hl2,hflux,p
allsel  
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,144   
deltim,3
solve
finish
！第一遍去3
/solu
antype,,rest
tim=144
*do,i,0,L-1.66,1.66
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-2.75,-0.77  
cm,hl1,nodes
sf,hl1,hflux,p
allsel  
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,0.77,2.75
cm,hl2,nodes
sf,hl2,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,180   
deltim,3
solve
finish
!第一遍回3
/solu
antype,,rest
tim=180   
*do,i,0,L/1.66-1,1
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-2.75,-0.77  
cm,hl1,nodes
sf,hl1,hflux,p
allsel  
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,0.77,2.75
cm,hl2,nodes
sf,hl2,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,216  
deltim,3
solve
finish
!第一遍去4
/solu
antype,,rest
tim=216
*do,i,0,L-1.66,1.66
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-0.99,0.99
cm,hl1,nodes
sf,hl1,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,252   
deltim,3
solve
finish
!第一遍回4
/solu
antype,,rest
tim=252   
*do,i,0,L/1.66-1,1
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-0.99,0.99
cm,hl1,nodes
sf,hl1,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,288
deltim,3
solve
finish

!第二遍去1
/solu
antype,,rest
tim=288
*do,i,0,L-1.66,1.66
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-6.25,-4.27  
cm,hl1,nodes
sf,hl1,hflux,p
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,4.27,6.25
cm,hl2,nodes
sf,hl2,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,324     
deltim,3
solve
finish
!第二遍回1
/solu
antype,,rest
tim=324   
*do,i,0,L/1.66-1,1
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-6.25,-4.27
cm,hl1,nodes
sf,hl1,hflux,p
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,4.27,6.25
cm,hl2,nodes
sf,hl2,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,360  
deltim,3
solve
finish
!第二遍去2
/solu
antype,,rest
tim=360
*do,i,0,L-1.66,1.66
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-4.5,-2.52  
cm,hl1,nodes
sf,hl1,hflux,p
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,2.52,4.5
cm,hl2,nodes
sf,hl2,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,396   
deltim,3
solve
finish
!第二遍回2
/solu
antype,,rest
tim=396  
*do,i,0,L/1.66-1,1
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-4.5,-2.52
cm,hl1,nodes
sf,hl1,hflux,p
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,2.52,4.5
cm,hl2,nodes
sf,hl2,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,432
deltim,3
solve
finish
!第二遍钢轮去3
/solu
antype,,rest
tim=432
*do,i,0,L-1.66,1.66
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-2.75,-0.77  
cm,hl1,nodes
sf,hl1,hflux,p
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,0.77,2.75
cm,hl2,nodes
sf,hl2,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,468   
deltim,3
solve
finish
!第二遍钢轮回3
/solu
antype,,rest
tim=468     
*do,i,0,L/1.66-1,1
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-2.75,-0.77  
cm,hl1,nodes
sf,hl1,hflux,p
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,0.77,2.75
cm,hl2,nodes
sf,hl2,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,504  
deltim,3
solve
finish
!第二遍钢轮去4
/solu
antype,,rest
tim=504
*do,i,0,L-1.66,1.66
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-0.99,0.99
cm,hl1,nodes
sf,hl1,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,540   
deltim,3
solve
finish
!第二遍钢轮回4
/solu
antype,,rest
tim=540     
*do,i,0,L/1.66-1,1
tim=tim+2
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-0.99,0.99
cm,hl1,nodes
sf,hl1,hflux,p
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,576
deltim,3
solve
finish

/solu
antype,,rest
deltim,40
time,1656
outres,all,all
solve
Finish
!第三遍去1
/solu
antype,,rest
tim=1656
*do,i,0,L-1.66,1.66
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-6.25,-4.27  
cm,hl1,nodes
sf,hl1,hflux,q
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,4.27,6.25
cm,hl2,nodes
sf,hl2,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1680
deltim,4
solve
finish
!第三遍钢轮回1
/solu
antype,,rest
tim=1680     
*do,i,0,L/1.66-1,1
tim=tim+1
time,tim
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-6.25,-4.27
cm,hl1,nodes
sf,hl1,hflux,q
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,4.27,6.25
cm,hl2,nodes
sf,hl2,hflux,q  
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1704
deltim,4
solve
finish
!第三遍钢轮去2
/solu
antype,,rest
tim=1704
*do,i,0,L-1.66,1.66
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-4.5,-2.52  
cm,hl1,nodes
sf,hl1,hflux,q
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,2.52,4.5
cm,hl2,nodes
sf,hl2,hflux,q
allsel  
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1728   
deltim,4
solve
finish
!第三遍钢轮回2
/solu
antype,,rest
tim=1728   
*do,i,0,L/1.66-1,1
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-4.5,-2.52
cm,hl1,nodes
sf,hl1,hflux,q
allsel  
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,2.52,4.5
cm,hl2,nodes
sf,hl2,hflux,q
allsel  
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1752  
deltim,4
solve
finish
!第三遍钢轮去3
/solu
antype,,rest
tim=1752
*do,i,0,L-1.66,1.66
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-2.75,-0.77  
cm,hl1,nodes
sf,hl1,hflux,q
allsel  
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,0.77,2.75
cm,hl2,nodes
sf,hl2,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1776   
deltim,4
solve
finish
!第三遍回3
/solu
antype,,rest
tim=1776  
*do,i,0,L/1.66-1,1
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-2.75,-0.77  
cm,hl1,nodes
sf,hl1,hflux,q
allsel  
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,0.77,2.75
cm,hl2,nodes
sf,hl2,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1800
deltim,4
solve
finish
!第三遍去4
/solu
antype,,rest
tim=1800
*do,i,0,L-1.66,1.66
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-0.99,0.99
cm,hl1,nodes
sf,hl1,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1824   
deltim,4
solve
finish
!第三遍回4
/solu
antype,,rest
tim=1824   
*do,i,0,L/1.66-1,1
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-0.99,0.99
cm,hl1,nodes
sf,hl1,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1848
deltim,4
solve
finish

!第四遍去1
/solu
antype,,rest
tim=1848
*do,i,0,L-1.66,1.66
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-6.25,-4.27  
cm,hl1,nodes
sf,hl1,hflux,q
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,4.27,6.25
cm,hl2,nodes
sf,hl2,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1872     
deltim,4
solve
finish
!四遍回1
/solu
antype,,rest
tim=1872   
*do,i,0,L/1.66-1,1
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-6.25,-4.27
cm,hl1,nodes
sf,hl1,hflux,q
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,4.27,6.25
cm,hl2,nodes
sf,hl2,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1896
deltim,4
solve
finish
!第四遍去2
/solu
antype,,rest
tim=1896
*do,i,0,L-1.66,1.66
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-4.5,-2.52  
cm,hl1,nodes
sf,hl1,hflux,q
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,2.52,4.5
cm,hl2,nodes
sf,hl2,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1920   
deltim,4
solve
finish
!第四遍回2
/solu
antype,,rest
tim=1920  
*do,i,0,L/1.66-1,1
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-4.5,-2.52
cm,hl1,nodes
sf,hl1,hflux,q
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,2.52,4.5
cm,hl2,nodes
sf,hl2,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1944
deltim,4
solve
finish
!第四遍去3
/solu
antype,,rest
tim=1944
*do,i,0,L-1.66,1.66
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-2.75,-0.77  
cm,hl1,nodes
sf,hl1,hflux,q
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,0.77,2.75
cm,hl2,nodes
sf,hl2,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1968   
deltim,4
solve
finish
！第四遍回3
/solu
antype,,rest
tim=1968     
*do,i,0,L/1.66-1,1
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-2.75,-0.77  
cm,hl1,nodes
sf,hl1,hflux,q
allsel
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,0.77,2.75
cm,hl2,nodes
sf,hl2,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,1992
deltim,4
solve
finish
！第四遍去4
/solu
antype,,rest
tim=1992
*do,i,0,L-1.66,1.66
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,i,i+1.66
nsel,r,loc,x,-0.99,0.99
cm,hl1,nodes
sf,hl1,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,2016
deltim,4
solve
finish
！第四遍回4
/solu
antype,,rest
tim=2016   
*do,i,0,L/1.66-1,1
tim=tim+1
time,tim
nsubst,1
esel,s,type,,1
nsel,r,loc,y,Ht
nsel,r,loc,z,L-1.66*(i+1),L-1.66*i
nsel,r,loc,x,-0.99,0.99
cm,hl1,nodes
sf,hl1,hflux,q
allsel
outres,all,all
solve
SFDELE,all,hflux
*enddo
time,2040
deltim,4
solve
finish

蓝莲花

你这个太长了吧，关键问题交待清楚就行了，网格化分和约束都有可能导致这样

caicai1959

谢谢回复~

kt266a

1# caicai1959

您不能给幅图说明一下啊？