耦合场分析求助

mayanyun

小弟做耦合场分析，计算一段波导由于高频损耗引起的热分布，命令流如下：
现在的问题是热载荷无法转到热分析中，我查了一下，hf119单元的热流输出是hfluxavg1、hfluxavg2、hfluxavg3，但是solid87单元没有相应的自由度，而且SFA命令也没有hfluxavg1这个lab。而且hfluxavg1等必须先计算powerh,后才能有输出，不然是0，这些问题都是互相矛盾的，真的不知道怎么解决。
我的信箱：mayanyun@ihep.ac.cn
望您在百忙之中给予关注，谢谢。
/PREP7  
/UNITS,mks
smrt,off
*SET,CH,0.01                         ! WIDTH
*SET,CW,0.03                         ! HEIGHT   
*SET,CL,0.048                        ! LENGTH EQUAL TO ONE WAVELENGTH   
*SET,FREQ,8E9                        ! FREQUENCY (HZ.)  
*SET,LAMDA,3E8/FREQ                  ! CALCULATE WAVELENGTH
BLOCK,-CW/2,CW/2,0,CH,0,CL      ! CREATE WAVEGUIDE  
BLOCK,-CW/2-0.003,CW/2+0.003,0-0.003,CH+0.003,0,CL      ! CREATE copper
VOVLAP,all  
!!!!!!!!
ET,1,HF119                      ! HIGH FREQUENCY SOLID TETRAHEDRAL  
MP,MURX,1,1.                    ! REALTIVE PERMEABILITY
MP,PERX,1,1.                    ! RELATIVE PERMITTIVITY
!!!!!!!!!   
VSEL, , , ,1
VATT,       1, ,   1,0  
allsel,all  
ESIZE,.0015                      ! SET ELEMENT DIVISIONS
VMESH,1                        ! MESH VOLUME  
!!!!!!!!!   
ASEL,S,LOC,X,-CW/2  
ASEL,A,LOC,X,CW/2   
ASEL,A,LOC,Y,0  
ASEL,A,LOC,Y,CH
NSLA,S,1
SF,all,SHLD,5.8e+007,1,  
allsel,all  
!!!!!!!!!!  
!the high frequency simulation  
!!!!!!!!!!!
ASEL,S, , ,1
NSLA,S,1
SF,all,PORT,1  
HFPORT,1,RECT,0,TE10,IMPD,34.04/1000,36.07/1000+36.07/1000,1000000,0,90000000,0,
ALLSEL,ALL
!!!
ASEL,S, , ,2
NSLA,S,2
SF,all,PORT,2
HFPORT,2,RECT,0,TE10,IMPD,34.04/1000,36.07/1000+36.07/1000,,,,,
ALLSEL,ALL  
ASUM, DEFAULT  
!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!
et,2,solid87
et,3,surf152
mp,kxx,2,400
mp,c,2,0.39*1000
mp,dens,2,8900  
VSEL, , , ,3
VATT,       2, ,   2,0  
allsel,all  
FLST,5,4,5,ORDE,2   
FITEM,5,3   
FITEM,5,-6  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   3,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*
ESIZE,.0015  
VMESH,3                        ! MESH VOLUME  
amesh,3
amesh,4
amesh,5
amesh,6
FLST,5,4,5,ORDE,2   
FITEM,5,3   
FITEM,5,-6  
ASEL,S, , ,P51X
eplot
SFA,all,,FSIN,1,
DA,10,TEMP,30
finish
!!!!!!!!!!!!!!!!!!
/solu   
mfan,on                     ! Activate MFS analysis
mfel,1,1                    ! Field #1: rf
mfel,2,2,3                  ! Field #2: Thermal
mfor,1,2
mfti,2                      ! Time at end of analysis
mfdt,1                      ! One field loop within a stagger
mfit,2                      ! Max 5 stagger loops
mfre,all,0.5                ! Field transfer relaxation parameter
mffn,1,rf                   ! Field #1 filename
mffn,2,Thermal              ! Field #2 filename
antyp,HARMIC
eqslv,iccg,,2  
HARFRQ,FREQ  
mfcm,1                      ! Write rf analysis options
antype,static
eqslv,SPARSE,,2         
mfcm,2                      ! Write thermal analysis options
MFSURFACE,1,1,HFLU,2
solve
finish

/post1
file,thermal,rth             ! Thermal results file
set,last