请明白ansys电磁学方面的朋友帮我看看,我这个问题怎么回事

浏览:68250

请懂得ansys电磁学方面的朋友帮我看看,为什么我计算结果都是零。
目的是研究通电螺旋线圈的电磁场分布,以及对磁场中工件的电磁力作用,线圈中的电流是通过电路-电磁耦合加载的。
采用三维模型,我的命令流是参考类似的二维模型的命令流写的,二维的那个可以计算出结果,但是我的模型计算结果都是零,如下图:请问怎么回事?
命令流有点长,如下,希望高手能指点一下我,非常感谢啊。
1){R$B0N_087MO%9~4WL0FF.jpg

命令流如下:
FINISH
/CLEAR
/FILNAME,COUPLE,1
/TITLE,The Analysis of Muti-physics Couple
KEYW,MAGNOD,1
/PREP7
ET,1,SOLID97,1,0 !定义工件的单元类型
ET,2,SOLID97,3 !定义线圈(有限元区域)的单元类型
ET,3,SOLID97,0 !定义近场的单元类型
ET,4,INFIN111,1,0 !定义远场的单元类型
ET,5,CIRCU124,0 !定义电阻的单元类型
ET,6,CIRCU124,1 !定义电感的单元类型
ET,7,CIRCU124,2 !定义电容的单元类型
ET,8,CIRCU124,5 !定义线圈(电路)的单元类型
ET,9,MESH200,6

!设置模型几何尺寸
X1=500e-3 !设置工件尺寸及位置
Y1=500e-3
t=14e-3 !设置工件厚度

X2=3e-3 !设置线圈横截面尺寸
Z2=10e-3
r0=0.02 !线圈内径
g=5e-3 !线圈下表面到工件上表面距离
h=Z2/2+g !线圈中心线到工件上表面距离

X_2=250E-3 !设置线圈位置
Y_2=300E-3
Z_2=t+h
R3=500E-3 !设置空气场区域半径
rm=0.125
Z_4=t+g
Z_5=t+g+Z2
!定义实常数
R,1,0.02 !回路电阻多大
R,2,6E-6 !回路电感、起始电感电流
R,3,0.001,500 !电容、起始电容电压
LOCAL,12,1
ESYS,12
EMODIF
R,4,X2*Z2,21,7.92,0,1,0,1,0.7875 !绞线圈的横截面积、线圈总匝数、模型导体体积、表示电流流向的单位矢量、对称系数、填充因子
R,5,X2*Z2,21,7.92,1,1,0.7875 !电路中线圈的横截面积、线圈总匝数、模型导体体积、表示电流流向的单位矢量、对称系数、填充因子
EMUNIT,MKS
MP,MURX,1,200 !工件的相对磁导率
MP,RSVX,1,1.2E-7 !工件的电阻率
MP,MURX,2,1 !线圈的相对磁导率
MP,RSVX,2,1.72E-8 !线圈的电阻率
MP,MURX,3,1 !空气的相对磁导率
MP,DENS,1,7850 !工件密度
MP,KXX,1,14.7 !工件导热系数
MP,EX,1,2.02E11 !弹性模量
MP,ALPX,1,1.16E-5 !线膨胀系数
MP,NUXY,1,0.28 !泊松比
TB,MISO,1,1,2 !多线性等向强化
TBPT,DEFI,116.337E-5,235e6 !应变、应力
N,1,1.5*R3,0,-0.5*R3 !设置电路单元的节点位置
N,2,1.5*R3,0,0.5*R3
N,3,2.5*R3,0,0.5*R3
N,4,2.5*R3,0,-0.5*R3
!建立有限元区域的几何模型
BLOCK,0,X1,0,Y1,0,t !工件模型
WPOFF,X_2,Y_2,0
CYLIND,0,R3,0,Z_5 !中间近场模型
BLOCK,-rm,rm,-rm,rm,Z_4,Z_5

VSBV,2,1,SEPO,DELETE,KEEP !中间近场减去工件的布尔运算
NUMCMP,VOLU
VSBV,3,2,SEPO,DELETE,DELETE !中间近场减去线圈区域的布尔运算
NUMCMP,VOLU
NUMCMP,LINE

LOCAL,11,1,X_2,Y_2,Z_2 !线圈模型
j=100
*DO,theta,0,7560,20
r=r0+1.0583e-5*theta
K,j,r,theta
j=j+1
*ENDDO
*do,j,100,472,5
bsplin,j,j+1,j+2,j+3,j+4,j+5
*enddo
NUMCMP,LINE
csys,0
LSEL,S,LINE,,43,117
CM,L1,LINE
WPCSYS,-1
WPOFF,X_2,Y_2,Z_2
wprota,,90
RECTNG,-X2/2,X2/2,-Z2/2,Z2/2
NUMCMP,AREA
WPCSYS,-1
AGEN,1,11,,,r0,0,0,,,1
VDRAG,11,,,,,,L1
VSEL,S,VOLU,,3,77
VADD,ALL
NUMCMP,VOLU
WPOFF,X_2,Y_2,Z_5
WPROTA,90
WPROTA,,90
SPHERE,,R3,0,180 !上半球型近场区域
WPCSYS,-1
WPOFF,X_2,Y_2,0
WPROTA,90
WPROTA,,-90
SPHERE,,R3,180 !下半球型近场区域
NUMCMP,VOLU
APLOT
SAVE,Emage_3D_geom.db !保存几何模型文件
SAVE
!网格划分
!设置网格划分的单元尺寸
ESIZE_1=0.05 !工件网格单元尺寸
ESIZE_2=0.01 !线圈网格单元尺寸
ESIZE_3=0.1 !近场单元尺寸
VSEL,S,VOLU,,3 !设置线圈属性
VATT,2,4,2
!划分线圈网格
ASEL,S,AREA,,11
AATT,2,4,9
LSLA,S
LESIZE,ALL,,,6
AMESH,11
ASEL,S,LOC,Z,Z_2-Z2/2,Z_2+Z2/2
LSLA,ALL
ESIZE,ESIZE_2
VSWEEP,3,11
VSEL,S,VOLU,,4,5 !切分球体,以便划分近场空气网格
VSBW,ALL
WPROTA,,,90
VSBW,ALL
NUMCMP,VOLU
VSEL,S,VOLU,,2 !设置中间近场属性
VATT,3,,3
VSEL,S,VOLU,,4,11 !设置上下球形近场属性
VATT,3,,3
VSEL,S,VOLU,,2 !划分中间近场网格
ESIZE,ESIZE_3
MSHAPE,1,3D
MSHKEY,0
VMESH,2
AREFINE,7,10,1,1,,SMOOTH !局部网格细化
AREFINE,12,,1,1,,SMOOTH
ASEL,S,AREA,,13,17
AREFINE,13,17,1,1,,SMOOTH
VSEL,S,VOLU,,4,11 !划分球形近场网格
ESIZE,ESIZE_3
MSHAPE,0,3D
MSHKEY,1
VMESH,4,11

ASEL,S,AREA,,66 !生成上半球远场单元
ASEL,A,AREA,,76
ASEL,A,AREA,,91
ASEL,A,AREA,,101
CM,A_AIR2,AREA
CMSEL,S,A_AIR2
EXTOPT,ESIZE,1 !生成一层单元
TYPE,4 !定义单元类型
MAT,3 !定义材料
LOCAL,13,2,X_2,Y_2,Z_5
VEXT,A_AIR2,,,0.01 !延伸一层单元成为远场单元
ASEL,S,AREA,,111 !生成下半球远场单元
ASEL,A,AREA,,121
ASEL,A,AREA,,136
ASEL,A,AREA,,146
CM,A_AIR3,AREA
CMSEL,S,A_AIR3
EXTOPT,ESIZE,1 !生成一层单元
TYPE,4 !定义单元类型
MAT,3 !定义材料
LOCAL,13,2,X_2,Y_2,0
VEXT,A_AIR3,,,0.01 !延伸一层单元成为远场单元
ASEL,S,AREA,,18 !生成中间远场单元
ASEL,A,AREA,,19
CM,A_AIR1,AREA
CMSEL,S,A_AIR1
EXTOPT,ESIZE,1 !生成一层单元
TYPE,4 !定义单元类型
MAT,3 !定义材料
LOCAL,13,1,X_2,Y_2,Z_2+Z2/2
VEXT,A_AIR1,,,0.01 !延伸一层单元成为远场单元
NUMCMP,ALL
VSEL,S,VOLU,,1 !设置工件属性
VATT,1,,1
VSEL,S,VOLU,,1 !划分工件网格
LESIZE,ALL,ESIZE_1
MSHAPE,0,3D
MSHKEY,1
VMESH,1
NUMCMP,AREA
NUMCMP,VOLU
SAVE,Emage_3D_mesh.db !保存网格单元数据到文件
SAVE
ESEL,S,MAT,,2
*GET,NE,NODE,,NUM,MAX !画电路图
TYPE,7 !电容器
REAL,3
E,2,1
TYPE,5 !电阻
REAL,1
E,2,3
TYPE,8 !电路线圈
REAL,5
E,3,4,19729
TYPE,6 !电感
REAL,2
E,4,1
ESEL,S,MAT,,2 !耦合线圈上的电流和电动势降自由度
NSLE,S
CP,1,CURR,ALL
CP,2,EMF,ALL
ESEL,ALL
NSEL,ALL
PHYSICS,WRITE,emag
PHYSICS,CLEAR
FINISH
/SOLU
PHYSICS,READ,emag
CSYS,0
N1=NODE(0,0,0)
D,N1,VOLT,0
ALLSEL
NSEL,ALL
D,1,VOLT,0 !电路中1号节点接地
ESEL,ALL
EPLOT
ASEL,S,AREA,,86
ASEL,A,AREA,,166
ASEL,A,AREA,,181
ASEL,A,AREA,,196
ASEL,A,AREA,,206
ASEL,A,AREA,,226
ASEL,A,AREA,,241
ASEL,A,AREA,,256
ASEL,A,AREA,,266
ASEL,A,AREA,,291
NSLA,S,1
SF,ALL,INF !无穷远标志
CSYS,2
D,ALL,AZ,0
ALLSEL
ESEL,S,MAT,,1
CM,ARM,ELEM
FMAGBC,'ARM' !施加力标志
ANTYPE,TRANS !瞬态分析
OUTRES,ALL,ALL
TIME,5E-2
AUTOTS
DELTIM,1E-5
SOLVE
FINISH
暂时只需要到求解这里吧,就是想知道为什么结果都为零,是单元问题还是边界条件问题,还是电路-电磁没有耦合上?
请大神指点,这个命令流从建模到计算一共也就几分钟,请高手运行一下帮我看看问题到底在哪?小弟感激不尽啊
邀请回答 我来回答

当前暂无回答

回答可获赠 200金币

没解决?试试专家一对一服务

换一批
    App下载
    技术邻APP
    工程师必备
    • 项目客服
    • 培训客服
    • 平台客服

    TOP