请教：建了一个一个发射线圈和一个接收线圈感应2D模型，结果不对

baiyan118

这是用53和circuit124建立的场路耦合模型，模式是一个发射线圈和一个接收线圈的二维模型，采用电流源激励 想得到磁力线图但结果不对，实在找不出是哪里的错误，把命令流帖上来，劳烦有经验的朋友帮帮忙，给我提点建议。
    FINI
    /CLEAR,START
    /prep7
    /title, TWO-HARM
   
    ET,1,PLANE53,0
    ET,2,PLANE53,3
   
    mp,murx,1,1            !空气材料特性
    mp,murx,2,1500        
    mp,murx,3,1            !线圈材料特性
    mp,rsvx,3,1.7e-8
      
    rectng,0,57e-2,0,48.2e-2
    rectng,25.1E-2,26.1e-2,0,480.2e-3
    rectng,28.6E-2,29.6e-2,0,480.2e-3
     
    ASBA,1,2
    NUMCMP,AREA  !压缩面标号
    ASBA,2,1
    NUMCMP,AREA
   
    rectng,25.3E-2,25.9e-2,5E-3,7.5e-2
    rectng,25.3E-2,25.9e-2,30.8E-2,37.8e-2
      
    rectng,28.8E-2,29.4e-2,30.8E-2,37.8e-2
    rectng,28.8E-2,29.4e-2,5E-3,7.5e-2
      
    rectng,-3E-2,60e-2,-5E-2,50.2e-2
   
    NUMCMP,AREA
    AOVLAP,all !（面搭接）
    AGLUE, all
    /PNUM,KP,0  
    /PNUM,LINE,0
    /PNUM,AREA,1
    /PNUM,VOLU,0
    /PNUM,NODE,0
    /PNUM,TABN,0
    /PNUM,SVAL,0
    /NUMBEA,0
    NUMCMP,AREA
   
   
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!!!!!!!!!!!!!!!!!!!!参数赋值
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    N1=200  !发射线圈匝数
    N2=104  !接收线圈匝数  
   
    ASUM, DEFAULT           ! 计算面的各个状态
    ASEL,S,AREA, ,      3
    *GET,A,AREA,1,AREA      !将A3的面积求出赋值给变量A
    ALLSEL,ALL
   
    ASEL,S,AREA, ,      2
    *GET,B,AREA,1,AREA      !将A2的面积求出赋值给变量B
    ALLSEL,ALL
   
    leng_1=56.8e-3            !线圈在z方向的长度
    leng_2=51e-3            !线圈在z方向的长度
    dianzu=17E-9              !电阻率
    R1=92.4                   !发射线圈电阻
    R2=0.5                    !接收线圈电阻
    R3=1000000                !负载电阻
   
    asel,s,area,,3        ! primary  
    aatt,3,2,2  
    asel,s,area,,4        ! primary  
    aatt,3,4,2  
   
    asel,s,area,,2        ! secondary  
    aatt,3,1,2  
    asel,s,area,,5        ! secondary  
    aatt,3,3,2  
      
    asel,s,area,,1        ! window  
    aatt,1,1,1  
   
    asel,s,area,,6       ! window  
    aatt,1,1,1  
    alls
   
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!!!!!!!!!!!!!!!!!!!!网格智能划分
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    asel,s,area,,2,5
    MSHAPE,0,2D   !MSHAPE,网格形状（0：四边形，六面体 1：三角形 四面体），维数
    MSHKEY,1      !MSHKEY,网格划分方法（0：自由网格划分 1：使用扫掠 2：尽可能用扫掠,不可以的情况用自由网格
                  !智能划分在扫掠网格中不起作用
    AMESH,ALL
    alls
   
    AREFINE,2,5,,2,0,1,1
    alls
   
    SMRTSIZE,1
    asel,s,area,,1
    asel,A,area,,6
    MSHAPE,1,2D   !MSHAPE,网格形状（0：四边形，六面体 1：三角形 四面体），维数
    MSHKEY,0
    AMESH,ALL
    alls
   
     *SET,c_fac,.653184
    r,1,B,N2,leng_2,1,c_fac   !!!R,NUM,AREA,TURN,LENG,电流方向，槽满率
    r,3,B,N2,leng_2,-1,c_fac  !!!R,NUM,AREA,TURN,LENG,电流方向，槽满率
    r,2,A,N1,leng_1,1,c_fac   !!!R,NUM,AREA,TURN,LENG,电流方向，槽满率
    r,4,A,N1,leng_1,-1,c_fac  !!!R,NUM,AREA,TURN,LENG,电流方向，槽满率

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  
!!!!!定义电路模型及其与有限元模型的连接
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    allsel
    numcmp,node  !压缩节点数
      c=0
     
    R,5,0,1,0.05,0.001,0.001,0.248,
    !R,5,0,1,1000
    allsel
    *Get,Nmax,Node,0,Num,Max  
    N,Nmax+1,68e-2,80e-3,0   ! 定义节点I，节点的号，是默认最大加1，节点位置定位
    N,Nmax+2,68e-2,20e-3,0    !J
    N,Nmax+3,68e-2,50e-3,0  !k
  
    ET,3,CIRCU124,4,0    !定义type 3  4:独立电流源 0:直流交流电压
    TYPE,3  
    REAL,5  
    MAT,1  
    E,Nmax+1,Nmax+2,Nmax+3
   
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!!!!!!!!!!!!!线圈1
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    R,6,1,  
    allsel
    asel,s,area,,3
    nsla,,1                  ! 选中属于所选面的所有节点
    *get,nn1,node,,num,min   !获得所有选中节点的最小节点号
    allsel
    n,Nmax+4,83e-2,80e-3,0
    RMOD,6,15,0,1
    ET,4,CIRCU124,5,0    !定义type5
    TYPE,4  
    REAL,6  
    MAT,1
    E,Nmax+1,Nmax+4,nn1
   
    R,7,1,  
    allsel
    asel,s,area,,4
    nsla,,1                  ! 选中属于所选面的所有节点
    *get,nn2,node,,num,min   !获得所有选中节点的最小节点号
    allsel
    n,Nmax+5,83e-2,20e-3,0
    RMOD,7,15,0,2
    ET,5,CIRCU124,5,0    !定义type5
    TYPE,5  
    REAL,7  
    MAT,1
    E,Nmax+4,Nmax+5,nn2
   
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!!!!!!!!!!!!!!!!!!!!!!!发射线圈电阻
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    r,8,R1
    RMOD,8,15,0,1
    ET,6,CIRCU124,0,0  !定义type9，因为是circu124，0，0所以系统可以识别是电阻，在real中就可以填写电阻的常数了  r,11,2
    TYPE,6
    REAL,8
    MAT,1   
    E,Nmax+5,Nmax+2
   
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!!!!!!!!!!!!!!!!!接收线圈
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  
   
   
    R,9,1,  
    allsel
    asel,s,area,,2
    nsla,,1                  ! 选中属于所选面的所有节点
    *get,nn3,node,,num,min  !获得所有选中节点的最小节点号
    allsel
    n,Nmax+6,83e-2,0,0
    n,Nmax+7,83e-2,-60e-3,0
    RMOD,9,15,0,1
    ET,7,CIRCU124,5,0     !定义type5
    TYPE,7
    REAL,9
    MAT,1  
    E,Nmax+6,Nmax+7,nn3
   
   
    R,10,1,  
    allsel
    asel,s,area,,5
    nsla,,1                  ! 选中属于所选面的所有节点
    *get,nn4,node,,num,min  !获得所有选中节点的最小节点号
    allsel
    n,Nmax+8,68e-2,-60e-3,0
    RMOD,10,15,0,2
    ET,8,CIRCU124,5,0      !定义type5
    TYPE,8
    REAL,10
    MAT,1
    E,Nmax+7,Nmax+8,nn4
   
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!!!!!!!!!!!!!!!!!!!!!!!!接收线圈电阻
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    r,11,R2
    allsel
    n,Nmax+9,68e-2,0,0
    RMOD,11,15,0,2
    ET,9,CIRCU124,0,0  !定义type9，因为是circu124，0，0所以系统可以识别是电阻，在real中就可以填写电阻的常数了  r,11,2
    TYPE,9
    REAL,11
    MAT,1   
    E,Nmax+6,Nmax+9
   
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!!!!!!!!!!!!!!!!!!!!!!!!电压零点约束
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    allSEL,ALL
    NSEL,S,NODE,,Nmax+2
    NSEL,A,NODE,,Nmax+8
    D,ALL,VOLT,0            !电压零点约束
    allSEL,ALL
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!!!!!!!!!!!!!!!!!!!!!!!耦合自由度
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    asel,s,area,,2 ! 选接收线圈的左边耦合
    nsla, ,1
    cp,1,curr,all   
    cp,2,emf,all
    allSEL,ALL
   
    asel,s,area,,3  ! 选发射线圈的左边耦合
    nsla, ,1
    cp,3,curr,all
    cp,4,emf,all
    allSEL,ALL
   
    asel,s,area,,4  ! 选接收线圈的右边耦合
    nsla, ,1
    cp,5,curr,all   
    cp,6,emf,all
    allSEL,ALL
   
    asel,s,area,,5  ! 选发射线圈的右边边耦合
    nsla, ,1
    cp,7,curr,all   
    cp,8,emf,all
    allSEL,ALL
   
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!!!!!!!!!!!!!!!!边界条件磁通量平行边界条件
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    alls
    esel,,ename,,53
    nsle
    nsel,r,ext  
    d,all,az
    alls
    FINISH
   
     /solu
     antype,harmic !谐波(交流)求解
     eqslv,sparse
    !csys,1
     HARFRQ,26000             !设置分析频率
     solve
     save
     FINISH
     
     /post1
      set,1
      plf2d,27,0,10,1
     FINISH