flotran耦合电磁力到底该怎么进行耦合？为什么我总是做不对

自然而然

先对电磁场进行分析，电磁力结果如下图：

自然而然

用flotran141单元进行流场分析，耦合节点电磁力，结果如下图：

自然而然

检查耦合的电磁力载荷，发现耦合的根本就不对，磁场结果电磁力分布均匀，外部大 内部小，耦合的却是不均匀，而且内部大 外部小，而且数值相差太大。网格划分和模型都是一样的啊，怎么会这样呢？和网格划分有关系吗？
为了划分的规律一点，我把模型分成了两部分，划分结果如下图：

自然而然

用flotran做硫场分析时，耦合电磁力应该怎么耦合？是不是apply>fluid/cfd>forces>forces>from mag analysis，然后指向电磁场分析的结果，选择载荷步和实部结果就可以了？？耦合的时候应该注意什么？

自然而然

这是方坯电磁力计算的结果

自然而然

这是flotran耦合电磁力的流场结果。我不知道我到底哪里出错了。

自然而然

这是命令流，我是用gui方式做的，这是整理的log文件。我对命令流还不熟，望各位战友见谅。

/PREP7  

ET,1,FLUID141   

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,8000
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,VISC,1,,0.007

RECTNG,-0.03,0.03,-0.03,0.03,   
TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
!*  
ESIZE,0.001,0,  
MSHAPE,0,2D
MSHKEY,1
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FLDATA1,SOLU,TRAN,0
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,0
FLDATA1,SOLU,TURB,1
FLDATA1,SOLU,COMP,0
FLDATA1,SOLU,VOF,0  
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,0
FLDATA1,SOLU,ALE,0  
FLDATA1,SOLU,RDSF,1
!*  
/COM,,Steady State Analysis,0   
FLDATA2,ITER,EXEC,10000,
FLDATA2,ITER,OVER,1000,
FLDATA2,ITER,APPE,1000,
FLDATA3,TERM,VX,0.0001,
FLDATA3,TERM,VY,0.0001,
FLDATA3,TERM,VZ,0.0001,
FLDATA3,TERM,PRES,1e-008,   
FLDATA3,TERM,TEMP,1e-008,   
FLDATA3,TERM,ENKE,0.0001,   
FLDATA3,TERM,ENDS,0.0001,   
FLDATA5,OUTP,SUMF,10,   
!*  
FLDATA12,PROP,DENS,16   
FLDATA13,VARY,DENS,0
FLDATA12,PROP,VISC,16   
FLDATA13,VARY,VISC,0
FLDATA12,PROP,COND,0
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,0
FLDATA13,VARY,SPHT,0
!*  
FLDATA7,PROT,DENS,TABLE
FLDATA8,NOMI,DENS,-1
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,TABLE
FLDATA8,NOMI,VISC,-1
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS  
FLDATA7,PROT,COND,CONSTANT  
FLDATA8,NOMI,COND,-1,   
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,CONSTANT  
FLDATA8,NOMI,SPHT,-1,   
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
!*  
FLDATA26,STAB,MOME,0.5,
FLDATA26,STAB,PRES,0.5,
FLDATA26,STAB,TEMP,0.5,
FLDATA26,STAB,TURB,1,   
FLDATA26,STAB,VISC,0,   
!*  
FLDATA24,TURB,MODL,1
FLDATA24,TURB,CMU,0.09,
FLDATA24,TURB,C1,1.44,  
FLDATA24,TURB,C2,1.92,  
FLDATA24,TURB,SCTK,1,   
FLDATA24,TURB,SCTD,1.3,
!*  
LDREAD,FORC,1,2, ,2,'emag','rst','..\2d2\'  
FLST,5,4,4,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,LINE  
LSEL,R, , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
CMDELE,_Y   

DL,_Y1, ,VX,0,1
DL,_Y1, ,VY,0,1
DL,_Y1, ,VZ,0,1

CMDELE,_Y1  

FINISH

自然而然

谢谢，我耦合的方法对不对？有没有什么其他的注意事项？

JCHLI2000

你不是做的三维的吗？怎么流场是二维的？
感觉你的网格划分不太合理。

自然而然

三维的结果不好，我想先做二维的看一下。圆坯的网格可能不好，这我也想到了。所以我之后又做了个方坯的，这回网格划分应该问题了，但是耦合的结果还是不对。