急求耦合高手相助:电磁与流场耦合!!!!

mao8479549

我做的电磁与流体耦合,电磁的体载荷加到流场中.第一步求解后进如流场求解的时候它不运算,出现的错误是:
Elements reference a solid material number(mat>1).Yet no solid material properties were specified .Affected elements are:
  

我的命令流如下,希望高手指点啊,不然完全做不下去了,急等着毕业啊:
/TITLE,magics and fluids
/prep7  
RECTNG,0,0.12,0,0.293,  
RECTNG,0,0.065,0.04,0.253,  
RECTNG,0.075,0.092,0.04,0.055,  
RECTNG,0.075,0.092,0.058,0.073,
RECTNG,0.075,0.092,0.076,0.091,
RECTNG,0.075,0.092,0.094,0.109,
RECTNG,0.075,0.092,0.112,0.127,
RECTNG,0.075,0.092,0.13,0.145,  
RECTNG,0.075,0.092,0.148,0.163,
RECTNG,0.075,0.092,0.166,0.181,
RECTNG,0.075,0.092,0.184,0.199,
RECTNG,0.075,0.092,0.202,0.217,
RECTNG,0.075,0.092,0.22,0.235,  
RECTNG,0.075,0.092,0.238,0.253,  建立几何模型
FLST,2,14,5,ORDE,2  
FITEM,2,1   
FITEM,2,-14
AOVLAP,P51X
!*  
ET,1,PLANE13
!*  
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,5,0
!*  
!*  
MPTEMP,,,,,,,,  定义材料属性
MPTEMP,1,0  
MPDATA,DENS,1,,2740
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,1,,8
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,RSVX,1,,2.1e-7   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,2,,1
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,RSVX,2,,1.75e-8  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,3,,1
CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,12,5,ORDE,2  
FITEM,5,3   
FITEM,5,-14
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,      15
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       3, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,2,1,5,ORDE,1   
FITEM,2,2   
AESIZE,P51X,0.008,  
MSHAPE,0,2D
MSHKEY,1
!*  
CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
APLOT   
MSHAPE,1,2D
MSHKEY,0
!*  
FLST,5,13,5,ORDE,2  
FITEM,5,3   
FITEM,5,-15
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   划分网格
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FLST,2,3,4,ORDE,2   
FITEM,2,1   
FITEM,2,-3  
SFL,P51X,INF
/REPLOT,RESIZE  
FLST,2,3,4,ORDE,3   
FITEM,2,8   
FITEM,2,57  
FITEM,2,-58
DL,P51X, ,ASYM  
!*  
ANTYPE,3
FINISH  
/SOL
HARFRQ,0,2000,  
NSUBST,100,
KBC,0   
!*  
FINISH  
/PREP7  
FLST,2,12,5,ORDE,2  
FITEM,2,3   
FITEM,2,-14
!*  
BFA,P51X,JS, , ,890000,0
PHYSICS,WRITE,magics,magics, ,  写电磁场物理文件

PHYSICS,CLEAR   

et,1,141                       定义流场单元
asel,s,,,3,15,1
esla,s,1
type,0                          （因为线圈部分和空气部分流场计算忽略，所以定义为0）
allsMPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2740
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,VISC,1,,0.00412       定义流场材料属性
FLST,5,2,4,ORDE,2   
FITEM,5,6   
FITEM,5,8   
CM,_Y,LINE  
LSEL,R, , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
CMDELE,_Y   
!*  
!*  
!*  
DL,_Y1, ,VX,0,1    施加x轴方向的速度边界为0
!*  
CMDELE,_Y1  
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,5   
FITEM,5,7   
CM,_Y,LINE  
LSEL,R, , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
CMDELE,_Y   
!*  
!*  
!*  
DL,_Y1, ,VY,0,1   施加y轴方向的速度边界为0
!*  
CMDELE,_Y1  
!*  
FLDATA1,SOLU,TRAN,1
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,,Transient Analysis,1  
*SET,_z4,-1
*SET,_z14,1
!*  
FLDATA4,TIME,STEP,-1
FLDATA4,TIME,ISTEP,0.1,
FLDATA4,TIME,NUMB,100,  
FLDATA4,TIME,TEND,1.0e06
FLDATA4,TIME,GLOB,200,  
FLDATA4,TIME,VX,0.01,   
FLDATA4,TIME,VY,0.01,   
FLDATA4,TIME,VZ,0.01,   
FLDATA4,TIME,PRES,1e-006,   
FLDATA4,TIME,TEMP,1e-006,   
FLDATA4,TIME,ENKE,0.01,
FLDATA4,TIME,ENDS,0.01,
FLDATA4A,STEP,OVER,0,   
FLDATA4,TIME,OVER,0
FLDATA4A,STEP,APPE,0,   
FLDATA4,TIME,APPE,1.0e6
FLDATA4A,STEP,SUMF,10,  
FLDATA4,TIME,SUMF,1.0e6
FLDATA4,TIME,BC,0   
FLDATA4,TIME,TEND,1000000,  
!*  
FLDATA12,PROP,DENS,0
FLDATA13,VARY,DENS,0
FLDATA12,PROP,VISC,0
FLDATA13,VARY,VISC,0
FLDATA12,PROP,COND,0
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,0
FLDATA13,VARY,SPHT,0
!*  
FLDATA7,PROT,DENS,CONSTANT  
FLDATA8,NOMI,DENS,2740,
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,CONSTANT  
FLDATA8,NOMI,VISC,0.00412,  
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
!*  
ACEL,0,9.8,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,
!*  
FLDATA24,TURB,ININ,0.01,
FLDATA24,TURB,INSF,0.01,
FLDATA24,TURB,RATI,1000,
FLDATA24,TURB,SCTT,0.85,
FLDATA24,TURB,SCTM,1,   
!*  
FLDATA24,TURB,WALL,VAND
FLDATA24,TURB,TRAN,11.5,
FLDATA24,TURB,KAPP,0.4,
FLDATA24,TURB,EWLL,9,   
FLDATA24,TURB,VAND,26,  
FLDATA24,TURB,KS,0,
FLDATA24,TURB,CKS,0.5,  
!*  
FLDATA24,TURB,BUC3,1,   
FLDATA24,TURB,BUC4,0,   
FLDATA24,TURB,BETA,0,   
!*  
PHYSICS,WRITE,fluids,fluids, ,      写流场的物理文件
finish  

/SOLU   
!*  
PHYSICS,READ,MAGICS                  读电磁物理文件并计算电磁场
/STATUS,SOLU
SOLVE   
finish  

/SOLU   
!*  
PHYSICS,READ,FLUIDS                  读流场物理文件
/assign,esave,fluids,esav            定义重启动文件  
/assign,emat,fluids,emat
LDREAD,FORC,,, ,0,'001','rmg',' '    加载电磁场的体载荷
SOLVE

mao8479549

过路的兄台,如有知晓的,可以指教 以下啊.感激!!!!!!!

yulanqiuhui

请问你是做哪一方面的啊，可以交流一下吗，留个联系方式吧：）

mao8479549

我是做电磁和流体耦合的.我的qq:37824451
欢迎互相交流学习

北极熊甲

...For example, a fluid may have magnetic characteristics. Any region that will be a fluid region at any time must have a material number of 1. If you cannot do this, you will have to modify the appropriate elements between performing the different physics solutions.

mao8479549

谢谢楼主的提点，我在定义材料号的时候把流体的材料号统一的定义为1号材料。线圈和空气分别为2号和3号。不管我怎么改，这个问题还是出现。我试了很多种方法，都解决不了。可否请楼主明示我的错误？小弟驽钝。请楼主赐教