!-----------------------------------
!
! current distribution conductor
!
!
!-----------------------------------


fini
/clear

/prep7


kappacu  = 58e6     ! Leitfähigkeit Kupfer
kappaal  = 58e6      ! Leitfähigkeit Al-Legierung

! conductor Al
ET,1,PLANE53               ! Define PLANE233 as element type
KEYOPT,1,1,1                ! Use AZ+VOLT
KEYOPT,1,3,0                ! Use plane analysis option
!KEYOPT,1,7,1                ! Condense forces at the corner nodes

! conductor Cu
ET,2,PLANE53               ! Define PLANE233 as element type
KEYOPT,2,1,1                ! Use AZ+VOLT
KEYOPT,2,3,0                ! Use plane analysis option
!KEYOPT,2,7,1                ! Condense forces at the corner nodes

!air
ET,3,PLANE53               ! Define PLANE233 as element type
KEYOPT,3,1,0                ! Use AZ
KEYOPT,3,3,0                ! Use plane analysis option
!KEYOPT,3,7,1                ! Condense forces at the corner nodes

et,4,infin110,,,0


mp,murx,1,1.0        ! Leiter Al - Permeabilitätskonstante
mp,rsvx,1,1/kappaal  ! Leiter Al - spez. Widerstand
mp,murx,2,1.0        ! Leiter Cu - Permeablitätskonstante
mp,rsvx,2,1/kappacu  ! Leiter Cu - spez. Widerstand
mp,murx,3,1.0        ! Luft Dielektrizitätskonstante (ideal)



! Einheit m
ri=4.0e-3/2
ra=5.0e-3/2

r= 25e-3/2

! mesh size
esizeleiter  = 0.05e-3   ! Elementgroesse im Leiter
esizeluft    = 0.2e-3    ! Elementgroesse Luft


! ---Geometrie---


! inner conductor
cyl4,0,0,0,0,ri,90
cyl4,0,0,0,90,ri,180
cyl4,0,0,0,180,ri,270
cyl4,0,0,0,270,ri,360
numcmp,all

! outer conductor
cyl4,0,0,ri,0,ra,90
cyl4,0,0,ri,90,ra,180
cyl4,0,0,ri,180,ra,270
cyl4,0,0,ri,270,ra,360
numcmp,all

! Luft
cyl4,0,0,0,0,r,90
cyl4,0,0,0,90,r,180
cyl4,0,0,0,180,r,270
cyl4,0,0,0,270,r,360
aovlap,all
!lovlap,all
numcmp,all

!/eof


!asel,s,area,,11$aplot
/pnum,area,1
aplot



! ---Vernetzung---
! A1:blank+ A2...AN:coil AN+1:Luft

MSHAPE,0,2D
! Mapped meshing benutzen
MSHKEY,2 

csys,1
! Vernetzung innerer Leiter
esize,esizeleiter
type,1
mat,1
amesh,1,4

! Vernetzung aeusserer Leiter
esize,esizeleiter
type,2
mat,2
amesh,5
amesh,7,9

! Vernetzung Luft
esize,esizeluft
type,3
mat,3
amesh,6
amesh,10,12

nsel,s,loc,x,r
esln
emodif,all,type,4
sf,all,inf
csys,0
alls



! ---Kopplung---
asel,s,area,,1,4
asel,a,area,,5
asel,a,area,,7,9
esla,
nsle, 
*get,nnodemin,node,,num,min  ! niedrigste Knotennummer
cp,next,volt,all          ! Kopplung des Potentials (Spannung)
alls


!save,,,,model

!/eof

! ---Lösung---
stromamp = 1          ! Stromamplitude (A)
freq = 5000           ! 5khZ

/solu
antype,harmic
harfrq,freq
!hrout,off,
!outres,all,all

f,nnodemin,amps,stromamp
solve


fini
/post1

!set,1,1  ! Realteil
!plf2d,27,0,10,1 

/plopts,minm,0
!-----------------
! realteil
set,1,1  ! Realteil
esel,s,mat,,1,2
etable,jtr,nmisc,7 !Stromdichtebetrag
pletab,jtr,avg

!------------------
! imaginaerteil
set,1,1,,1 
presol,elem
esel,s,mat,,1,2
etable,jti,nmisc,7 !Stromdichtebetrag
pletab,jti,avg

!-----------------
! berechnung absolutwert stromdichte
smult,quadjr,jtr,jtr
smult,quadji,jti,jti
sadd,sumj,quadjr,quadji,1,1
sexp,absj,sumj,,0.5

pletab,absj,avg

!-----------------
! berechnung gesamtstrom
esel,s,mat,,1,2            ! Leiter + Beschichtung
sabs,0 
etable,flaeche,volu        ! Elementflaeche
smult,gesamtstrom,absj,flaeche
ssum
*get,sstrom,ssum,,item,gesamtstrom

/eof

!---------------
esel,s,mat,,1,2
path,radius,2
ppath,1,,0,0,0
ppath,2,,ra,0,0
!etable,jt1,jt,z
pdef,jt2,etab,absj
/yrange,0,1.5e5
plpath,jt2

!---------------
!ples,jt,sum
!plnsol,jt,z,1   


!-----------------
esel,s,mat,,1,2           !Leiter
sabs,0
!etable,wider,nmisc,8  
etable,stromd,nmisc,7
etable,flaeche,volu        !Elementvolumen
!smult,strom,absj,flaeche
smult,strom,stromd,flaeche
ssum
*get,sstrom,ssum,,item,strom
gesstrom=sstrom