空气面网格划分不了
!左右两侧子导线均镜像得到
!/bach,list
finish
/clear
/title, Transmission Line
/filename,capacitance,1
keyw,magelc,1 !指定电场分布
/UNITS,MPA
/prep7
!V1=1050000
V1a=577350
V1b=-288675
V1c=-288675
V0=0
a=15 !子导线内径（mm）
b=10 !地线1的内径
c=10.025 !地线2的内径
x1=16000
y1=28000
x21=1000 !导线组2相对于导线组1的x坐标差值
y21=20500 !导线组2相对于导线组1的y坐标差值
x31=-1000 !导线组3相对于导线组1的x坐标差值
y31=41000 !导线组3相对于导线组1的y坐标差值
xd=21000
yd=83000
d=2*yd !空气区域外径
r0=2*d !远场单元外径
/pnum,area,1 !打开面区域编号
et,1,121 !8节点二维静电场单元
et,2,121 !空气
et,3,110,1,1 !8节点二维远场单元
emunit,epzr0,8.854e-6 !以
mksv单位制设置自由空间介电常数
mp,perx,1,1 !设置相对介电常数
mp,perx,2,0.8
mp,perx,3,0.6
local,11,1,x1,y1,0 !自定义11号柱坐标系
n,1,522.6,22.5
n,2,522.6,67.5
nwpave,1 !移动坐标系到n1节点处
cyl4,0,0,a,0
nwpave,2 !移动坐标系到n2节点处
cyl4,0,0,a,0
local,12,0,x1,y1,0 !自定义12号笛卡尔坐标系
arsym,x,all !以y-z平面为镜像对称面
arsym,y,all !以y-z平面为镜像对称面
agen,2,1,8,1,x21,y21 !复制A相8根导线，
agen,2,1,8,1,x31,y31 !复制A相8根导线
csys,0
wpcsys,-1,0 !回到原始坐标系
cyl4,xd,yd,b,0 !创建地线
allsel,all
arsym,x,all !以y-z平面为镜像对称面
nummrg, node !合并重合的节点
nummr,kpoi !合并重合的关键点
cyl4, 0,0,r0,0,,180 !远场单元内径
cyl4, 0,0,2*r0,0,,180 !远场单元外径
aovlap,all
numcmp,area
asel,s,area,,1,24,1! 选择导线的面
，
从面 1到面3，
以间隔1选择
aatt,1,1,1
asel,s, area,,25! 选择地线的面
aatt,2,1,1
asel,s,area,,26,49,1!选择导线的面
aatt,1,1,1
asel,s, area,,50! 选择地线的面
aatt,2,1,1
asel,s, area,,52 ! 选择空气的面
aatt,3,1,2
asel,s,area,,51! 选择远场的面
aatt,3,1,3
asel,all
csys,1
lsel,s,loc,x,1.5*r0
lsel,a,loc,y,1.5*r0 !选择远场单元径向两条线
lesize,all,,,1
esize,,32
mshape,0,2d
mshkey,1
amesh,51
asel,all
mshape,1,2d
mshkey,0
amesh,52
asel,all
smrtsiz,4
mshape,1,2d
mshkey,0
amesh,1,50,1
allsel,all !选择所有实体
nsel,s,loc,x,2*r0 !选择远场外边界节点
sf,all,inf !设置远场标志
nsel,s,loc,x,2*r0 !选择远场外边界节点
D,ALL,VOLT,V0
*do,i,1,24,1
Asel,s,area,,i
nsla,s,0
cp,i,volt,all
*enddo
allsel,all
*do,i,26,49,1
Asel,s,area,,i
nsla,s,0
cp,i,volt,all
*enddo
allsel,all
*do,i,1,8,1
Asel,s,area,,i
nsla,s,0
D,all,volt,V1a
*enddo
allsel,all
*do,i,42,49,1
Asel,s,area,,i
nsla,s,0
D,all,volt,V1a
*enddo
allsel,all
*do,i, 9,16,1
Asel,s,area,,i
nsla,s,0
D,all,volt,V1b
*enddo
allsel,all
*do,i, 34,41,1
Asel,s,area,,i
nsla,s,0
D,all,volt,V1b
*enddo
allsel,all
*do,i,17,24,1
Asel,s,area,,i
nsla,s,0
D,all,volt,V1c
*enddo
allsel,all !选择所有实体
*do,i,26,33,1
Asel,s,area,,i
nsla,s,0
D,all,volt,V1c
*enddo
allsel,all
Asel,s,area,,25,50,25
nsla,s,0
D,ALL,VOLT,V0
csys,0 !激活全局笛卡尔坐标系
wpcsys,-1,0 !回到原始坐标系
nsel,s,loc,y,0 !选择地边界节点
D,ALL,VOLT,V0
Csys,1
nsel,s,loc,x, r0
nsla,s,0
D,ALL,VOLT,V0
Nsel,all
Finish
/solution
Antype,static
ALLSEL,all
EQSLV,
ICCG, ,0, ! 使用 ICCG 求解
solve
finish
/post1 ! 通用后处理
Etable,efx,ef,x !存储电场力
Etable,efy,ef,y
/number,1
PLNSOL,volt ！显示电位等值云图，线图PlotCtrls->DeviceOptions->Vector mode(wireframe)
PLVECT,efx,efy
PLNSOL, EF,SUM, 0
PLNSOL, EF,X, 0
PLNSOL, EF,y, 0
