[bn_ansys]

Добрый день,

Я работаю над моделью силового кабеля ЛЭП. Однако, я получаю странные результаты даже если в модели есть только один единственный проводник. Код, который я использую, приведён ниже.

Если построить график плотности тока начиная от центра проводника, то видно, что ток в центре проводника и его окрестностях принимает отрицательные значения. Значения тока ниже, чем при постоянном токе, можно ожидать из-за поверхностного эффекта, однако, на сколько мне известно, ток не должен принмать отрицательных значений.

Может ли кто-нибудь помочь выяснить, в чем тут дело?

Заранее спасибо за Вашу помощь.


http://imgur.com/Wx9MOHZ - mesh

http://imgur.com/J0lIxbr - current density graph

/PREP7

d_al=30/1000 !diameter of aluminium conductor
I_cable=1 !current
ff=250 !frequency
R=0.0000000282 !Resistivity of aluminium

r_air=0.3 !radius of circle (air around cable) (m)
airlayer=0.005 !thin air layer around the cable

pi=acos(-1)

ET,1,PLANE53 !PLANE53, AZ DOF, (FOR AIR)
ET,2,PLANE53,1 !PLANE53, AZ VOLT DOF, (FOR CONDUCTOR)

MP,MURX,1,1 !material one is air around the cable
MP,RSVX,1,2e16
MP,MURX,2,1
MP,RSVX,2,R !material 2 is aluminium
MP,MURX,4,1 !material 4 is thin air layer around the cable
MP,RSVX,4,2e16

CYL4,1,1,r_air
CYL4,1,1,d_al/2
CYL4,1,1,(d_al/2+airlayer)

AOVLAP,ALL

*GET,MinAreaNum,AREA,,NUM,MIN
*GET,MaxAreaNum,AREA,,NUM,MAX
ASEL,S,AREA,,MinAreaNum !select the central wire (minimum area number)
AATT,2,,2,0,
ASEL,S,AREA,,(MaxAreaNum-1)
AATT,1,,1,0, !select the air around cable (maximum area number - 1)
ASEL,S,AREA,,MaxAreaNum
AATT,4,,1,0, !select the thin air layer around the cable

SMRT,1
MSHAPE,1,2D
MSHKEY,0
ASEL,S,MAT,,4
AMESH,ALL
AREFINE,ALL,,,2,
ASEL,S,MAT,,2
AMESH,ALL
AREFINE,ALL,,,3,
SMRT,2
ASEL,S,MAT,,1
AMESH,ALL

ANTYPE,HARMIC

LOCAL,83,1,1,1,0, , , ,1,1,
CSYS,83, !local coordinate system

ESEL,ALL
NSEL,S,LOC,X,r_air
D,ALL,AZ,0 !boundary conditions

ESEL,S,MAT,,2
NSLE,S
CP,1,VOLT,ALL
*GET,N1,NODE,,NUM,MIN !applying current on one node - lowest number node in selected set (material 2 - aluminium)
F,N1,AMPS,I_cable

ASEL,S,MAT,,2
ASUM
*GET,AlumArea,AREA,,AREA
*SET,J_cab,I_cable/Alumarea !calculation of conductor's area and current density (DC)

CSYS,0

ALLSEL,ALL

FINISH

/SOLU

HARFRQ,ff
SOLVE
FINISH

/POST1

SET,1,1
ETABLE,currdens,NMISC,7

PATH,JSPath,2,80,80
d_cable=1 !1st point - center of cable
dist_cable=0.005 !2nd point - distance from the cable's surface
PPATH,1,,d_cable,1,
PPATH,2,,d_cable+d_al+dist_cable,1,
PDEF,jspath,ETAB,currdens !current density plot

ALLSEL,ALL