baiyan118 发表于 2010-8-19 10:31

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

这是用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 34:独立电流源 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
页: [1]
查看完整版本: 请教:建了一个一个发射线圈和一个接收线圈感应2D模型,结果不对