这是我建的.就发觉剖分的网格不是关于Y轴对称的.这个可以避免吗?
/prep7
et,1,plane121
mp,perx,1,4.4 !oil
mp,perx,2,1 !copper
mp,rsvx,1,1e12 !oil
mp,rsvx,2,0.0175 !copper
k,1,750,0
k,2,750,25000
k,3,-750,25000
k,4,-750,0
k,5,0,-1289.2881
k,6,0,-1292.5381
k,7,5.629,-1289.2881
k,8,-5.629,-1289.2881
lstr,1,2
lstr,2,3
lstr,3,4
lstr,4,8
larc,6,7,5,6.5
larc,6,8,5,6.5
lstr,7,1
lsel,all
al,all
kdele,5
gaplength=50000
k,11,25000,-1292.5381-gaplength !构建板电极
k,12,25000,-3000-gaplength
k,13,-25000,-3000-gaplength
k,14,-25000,-1292.5381-gaplength
*do,i,11,13,1
lstr,i,i+1
*enddo
lstr,14,11
lsel,s,line,,8,11
al,all
rectng,-125000,125000,-150000,150000
rectng,-30000,30000,-10000-gaplength,10000-gaplength
rectng,-10000,10000,-10000,30000
rectng,-50000,50000,-40000-gaplength,60000
allsel,all
aovlap,all
numcmp,all
allsel,all
asel,s,area,,5
aatt,1,,1
SMRTSIZE,1
mshape,1
amesh,all
asel,s,area,,3
aatt,1,,1
SMRTSIZE,4
mshape,1
amesh,all
asel,s,area,,6
aatt,1,,1
SMRTSIZE,4
mshape,1
amesh,all
asel,s,area,,4
aatt,1,,1
SMRTSIZE,6
mshape,1
amesh,all
EREF,ALL, , ,1,0,1,1
CMDELE,_Y1
MCHECK
MCHECK
/solution
allsel,all
asel,s,area,,2
lsel,s,ext
dl,all,,volt,0
lsel,s,line,,12,15
dl,all,,volt,0
asel,s,area,,1
lsel,s,ext
dl,all,,volt,35000
allsel,all
solve
finish
/POST1
plnsol,ef,sum
N=20
path,A,N+20,20,1
ppath,1,,0,-1292.5381
ppath,2,,0,-1293
ppath,3,,0,-1296
ppath,4,,0,-1300
ppath,5,,0,-1320
ppath,6,,0,-1360
ppath,7,,0,-1390
ppath,8,,0,-1420
ppath,9,,0,-1440
ppath,10,,0,-1460
ppath,11,,0,-1480
*do,i,1,8,1
ppath,i+11,,0,-1292.5381-gaplength*i/100
*enddo
*do,i,1,N,1
ppath,i+19,,0,-1292.5381-gaplength/6-gaplength*i/(1.5*N)
*enddo
ppath,N+20,,0,-1292.5381-gaplength
pdef,EFSUM,ef,sum,avg
paget,items,table
*dim,y_path,array,N+20,1
*dim,efsum_path,array,N+20,1
*do,i,1,N+20
y_path(i)=items(i,2)
efsum_path(i)=items(i,5)
*enddo
*cfopen,EF_RESULT_5_0mm_ac,txt,,
*do,i,1,N+20
f1=y_path(i)
f2=efsum_path(i)
*vwrite,f1,f2
(f12.0,f30.20)
*enddo
*cfclos
谢谢了. |