/com Simulation 2D de l'effet de peau et de proximité : simulation avec 2 conducteurs et de l'air autour


/prep7

emunit,mks


/com entrees de programme pour la simulation
freq=1e4

et,1,plane13,6			! air
mp,murx,1,1
mp,rsvx,1,1e8


et,2,plane13,6			! materiau conducteur
mp,murx,2,1
mp,rsvx,2,1.72e-8


/com declaration des dimensions
r1=0.001	! rayon des conducteurs
x=7.5e-3	! separation des cercles representant les conducteurs
x1=15e-3
y1=15e-3


/com placement des Kpoints
k,5,x1,y1		! cercle representant l air
k,6,-x1,y1
k,7,-x1,-y1
k,8,x1,-y1


local,11,1,0,0,0,,0,,1,1	! coordonnees cylindriques
/com trace des aires

cyl4,-x/2,0,0,0,r1,360
cyl4,+x/2,0,0,0,r1,360
cyl4,-x-x/2,0,0,0,r1,360
cyl4,x+x/2,0,0,0,r1,360
a,5,6,7,8			! a5

aovlap,all			! l aire 5 devient l aire 6

asel,s,,,1			! materiau conducteur
asel,a,,,2
asel,a,,,3
asel,a,,,4
aatt,2,,2

asel,s,,,6			! materiau air
aatt,1,,1


ksel,all
lsel,all
asel,all

smrtsize,2			! precision du maillage
amesh,all			! maillage de toutes les aires

asel,s,,,1,,,1
CP,1,volt,all
I=1				! courant applique
asum
*get,area,area,1,area
*get,N1,node,,num,min
F,N1,amps,I

asel,s,,,2,,,1
CP,2,volt,all
I=1
asum
*get,area,area,2,area
*get,N2,node,,num,min
F,N2,amps,I

asel,s,,,3,,,1
CP,3,volt,all
I=1
asum
*get,area,area,3,area
*get,N3,node,,num,min
F,N3,amps,I

asel,s,,,4,,,1
CP,4,volt,all
asum
*get,area,area,4,area
*get,N4,node,,num,min
F,N4,amps,I

allsel,all


nsel,s,ext
d,all,az,0

nsel,all
esel,all
finish

/solu

antype,harmic
harfreq,freq
solve

/post1
set,1,1


PLVECT,B, , , ,VECT,ELEM,ON,0   

asel,s,mat,,2		! materiau des conducteurs
esla
nsle
eplo

PLESOL,JT,SUM,0,	! affichage de la densite de courant