我做的电磁与流体耦合,电磁的体载荷加到流场中.第一步求解后进如流场求解的时候它不运算,出现的错误是:
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 |