一个简单2D场路耦合模型，求解时有点问题（附程序）

arthur0203069

自己建的一个2D场路耦合的例子，模型建立没有什么问题，求解时总出现错误
Error：Element 1 has negative or zero cross section specified.check real constant table11
请高手指教，是不是我的求解程序有问题？
下面是我的程序
/PREP7                  !开始前处理程序
/UNITS, SI              !选取国际单位
ET,1,53,0               ! AIR
ET,2,53,3               ! STATOR COIL(CIRCUIT-COUPLED STRANDED COIL)
ET,3,53,0               ! 铁心和衔铁

EMUNIT, MKS           !选取（米、千克、秒）电磁单位制
MP,MURX,1,1           !AIR
MP,MURX,2,1           !STATOR WINDING
mp,rsvx,2,0.0434e-6   !导体电阻率
mp,murx,3,2000

LOCAL,11,0            !选取直角坐标系统
csys,11
K, 1,-0.05,  0.000
K, 2,-0.05,  0.100
K, 3, 0.05,  0.100
K, 4, 0.05,  0.000
K, 5,-0.02,  0.010
K, 6,-0.02,  0.030
K, 7, 0.02,  0.030
K, 8, 0.02,  0.010
K, 9,-0.02,  0.050
K, 10,-0.02,  0.070
K, 11, 0.02,  0.070
K, 12, 0.02,  0.050
K, 13,-0.04,  0.030
K, 14,-0.04,  0. 075
K, 15,-0.02,  0.075
K, 16, 0.02,  0.075
K, 17, 0.04,  0.075
K, 18, 0.04,  0.030
K, 19,-0.04,  0.080
K, 20,-0.04,  0.090
K, 21, 0.04,  0.090
K, 22, 0.04,  0.080
LSTR, 1, 2                   ! L1  第一条线（Straight line）
LSTR, 2, 3
LSTR, 3, 4
LSTR, 4, 1
LSTR, 5, 6                   ! L5
LSTR, 6, 7
LSTR, 7, 8
LSTR, 8, 5
LSTR, 9, 10
LSTR, 10, 11                 ! L10
LSTR, 11, 12
LSTR, 12, 9
LSTR, 6, 13
LSTR, 13, 14
LSTR, 14, 15                 ! L15
LSTR, 15, 10
LSTR, 11, 16
LSTR, 16, 17
LSTR, 17, 18
LSTR, 18, 7                  ! L20
LSTR, 19, 20
LSTR, 20, 21
LSTR, 21, 22
LSTR, 22, 19
LSTR, 14, 19                 ! L25
LSTR, 15, 16
LSTR, 17, 22                 ! L27
LESIZE,1,,,20                !求解区边界，分得粗一些
LESIZE,2,,,20
LESIZE,3,,,20
LESIZE,4,,,20
LESIZE,5,,,5                 !绕组下部，分得稍密一些
LESIZE,6,,,10
LESIZE,7,,,5
LESIZE,8,,,10
LESIZE,9,,,5                 !绕组上部，分得稍密一些
LESIZE,10,,,20               !气隙部分要分得最密
LESIZE,11,,,5
LESIZE,12,,,10
LESIZE,13,,,5                !铁心
LESIZE,14,,,10
LESIZE,15,,,10               !气隙部分要分得最密
LESIZE,16,,,3
LESIZE,17,,,3                ! 与L16相同
LESIZE,18,,,10               ! 与L15相同
LESIZE,19,,,10               ! 与L14相同
LESIZE,20,,,5                ! 与L13相同
LESIZE,21,,,4                ! 衔铁部分要分得密
LESIZE,22,,,30
LESIZE,23,,,4
LESIZE,24,,,40
LESIZE,25,,,3
LESIZE,26,,,20
LESIZE,27,,,3
AL, 5, 6, 7, 8               ! A1 绕组下部
AL, 9, 10, 11, 12            ! A2 绕组上部

LSEL,S,,,13,16,1,1           ! 选一组线 13,14,15,16
LSEL,A,,,17,20,1             ! 加选17,18,19,20
LSEL,A,,,6,12,3              ! 加选6,9,12
LSEL,A,,,11                  ! 加选11
AL,ALL                       ! 由所选的所有线段组成面A3
! 产生衔铁和气隙部分面积
ALLSEL                      ! 选取所有点线面
AL,21,22,23,24              ! A4 衔铁
AL,15,25,24,27,18,26        ! A5 主气隙
AL,10,16,26,17              ! A6 线圈与衔铁间气隙
AL,1,2,3,4                  ! A7  4条外边界线构成
AOVLAP,ALL                  ! 所有面重叠后产生一个新面A8
                            ! A7 面不再存在
allsel
NUMCMP,ALL

asel,s,,,5,7,1,1
AATT,1,,1
asel,s,,,1,,,1
AATT,2,11,2
asum
*get,aa,area,,area
r,11,aa,150,,1,0.8   !下线圈单元的实常数
asel,s,,,2,,,1
AATT,2,11,2
asum
*get,aa,area,,area
r,13,aa,150,,-1,0.8  !上线圈单元的实常数

asel,s,,,3,4,1,1
AATT,3,,3

ALLSEL
AMESH,ALL
finish

/prep7
local,11,0
csys,11

et,10,124,4,1 !独立电压源(用于建立单相电源)
r,10,0,115,60 !电压115V,相位为0,f=60Hz

et,11,124,5   !电路耦合绞线圈,用于线圈
r,11
r,13
et,12,124,0   !电阻,用于端环阻抗
r,12,1e-10

!(线圈自由度耦合)
asel,s,,,1,,,1
nsla,,1
cp,1,curr,all
cp,2,emf,all
asel,s,,,2,,,1
nsla,,1
cp,3,curr,all
cp,4,emf,all

allsel
numcmp,node
*get,nmax,node,,num,max
nmax1=nmax+1
nmax2=nmax+2
nmax3=nmax+3
nmax4=nmax+4
nmax11=nmax+5
nmax22=nmax+6
nmax33=nmax+7
nmax44=nmax+8
nmax5=nmax+9
asel,s,,,1,,,1
nsla,,1
*get,nn1,node,,num,min
xx=nx(nn1)
yy=ny(nn1)
n,nmax1,xx,yy,0.05
n,nmax11,xx,yy,0.10
n,nmax2,xx,yy,-0.05
n,nmax22,xx,yy,-0.10

asel,s,,,2,,,1
nsla,,1
*get,nn2,node,,num,min
xx=nx(nn2)
yy=ny(nn2)
n,nmax3,xx,yy,0.05
n,nmax33,xx,yy,0.10
n,nmax4,xx,yy,-0.05
n,nmax44,xx,yy,-0.10
xx=nx(nmax22)
yy1=ny(nmax22)
yy2=ny(nmax44)
yy=0.5*(yy1+yy2)
n,nmax5,xx,yy,0.10
allsel
type,11
real,11
e,nmax2,nmax1,nn1
real,13
e,nmax3,nmax4,nn2

type,12
real,12
e,nmax1,nmax11
e,nmax2,nmax22
e,nmax3,nmax33
e,nmax4,nmax44
e,nmax22,nmax44

type,10
real,10
e,nmax33,nmax11,nmax5

allsel
save
finish

/solu
antype,trans
trnopt,full
allsel
tm_start=1e-5
tm_end=30
tm_incr=1.5

*do,tm,tm_start,tm_end,tm_incr
time,tm
csys,11
LSEL,S,,,1,4,1,1            !选择L1,L2,L3,L4
NSLL,S                      !选择线上的所有节点
D,ALL,AZ,0                  !对边界节点赋值AZ=0
allsel
d,nmax11,volt,0  !接地
allsel
neqit,50
cnvtol,a,,1e-5
solve
finish

arthur0203069

怎么没人回啊，版主还有各位高手帮帮忙啊

arthur0203069

好像明白一点了
r,11 和r,13定义了两次，把后面那次的定义去掉就行了
求解过程没加*enddo
现在能求解了，后处理还有点不明白

liubinlw

楼主的例子做出的电流波形是这样子的，看看电流的幅值好像很大的样子，
单位应该是安培，这个结果好像不对的样子呢。