!Wicklung schaltungsgespeist
!3D-Modell
!-------------------------------------------------
finish  $/clear
*abbr,pro_4_4,/input,pro_4_4.inp
!------PARAMETERFESTLEGUNG------------------------
!------Geometrie----------------------------------
wiho=50e-3		!Wicklungs-Höhe
wiradin=39e-3		!Wicklungs-Radius innen
wiradau=40e-3		!Wicklungs-Radius außen
widi=wiradau-wiradin	!Wicklungs-Dicke
wianz=40		!Windungszahl
radau=200e-3		!Radius außen
!------Schaltungsparameter------------------------
kapazit=500e-6		!Kapazität
spannung=-50		!Ladespannung
widerst=3e-3		!Vorwiderstand
indukt=30e-9		!Kreisinduktivität
!------Sonstige Kenngrößen------------------------
wkappa=56e6		!spez. el. Leitfähigkeit der Wicklungen
pi=4*atan(1)            
mw=4e-3			!Maschenweite
abm=radau		!Schaltungsabmessungen
wifla=wiho*widi		!Wicklungs-Fläche 
wivol=wiho*(pi*wiradau*wiradau-pi*wiradin*wiradin)		
!------Elementetypen------------------------------
/prep7			!Aufruf des Preprozessors
et,1,solid97,0		!3D-Element->static domain
et,3,solid97,3		!3D-Element->circuit coupled stranded coil
et,4,infin111,1		!Randelement
et,7,circu124,2		!Kondensator
et,8,circu124,0		!Widerstand
et,9,circu124,5		!Koppelelement->3D/stranded coil
et,10,circu124,1	!Induktivität
!------Materialeigenschaften------------------------
mp,murx,1,1       	!rel. Permeabilität (Luft)
mp,murx,2,1		!rel. Permeabilität (Wicklung)
mp,rsvx,2,1/wkappa	!spez. el. Leitfähigkeit (Wicklung)
!------Real-Konstanten---------------------------------
r,1,wifla,wianz,wivol,0,1,0
r,7,kapazit,spannung
r,8,widerst
r,9,1
r,10,indukt
!------Geometrie------------------
!------Wicklung---------------
wplane,0,5000e-5,0
cylind,wiradin,wiradau,-wiho/2,wiho/2,-90,90
cylind,wiradin,wiradau,-wiho/2,wiho/2,90,270
vglue,all              !V1,3
numcmp,volu,all        !V1,2
wplane,0,0,0

!------Luft--------
!sphere,radau-.1e-3,radau,-90,90    !V3
!/eof
sphere,,radau,-90,90   !V3
sphere,,radau,90,270   !V4
vglue,3,4              !V3,5
numcmp,volu            !V3,4
vovlap,all             !V1,2,5,6
numcmp,volu            !V1,2,3,4
!V1,2:Wicklung   V3,4:Umgebung
!-------Vernetzung der Wicklung-------------
esize,mw*2  		!Maschenweite
mshkey,1		!(0)freie (1)mapped (2)mapped/freie Vernetzung
mshape,2,3D		!Raumvernetzung mit tetraederförmigen Elementen
type,3	$mat,2	$real,1	!Wicklung
vmesh,1,2
!------Vernetzung der Umgebung-------------
esize
mshkey,0		
mshape,1,3D		
type,1	$mat,1	
vmesh,3,4	
!/eof
!------Randelememte------------------------
csys,2
type,4
esize,,1
asel,s,ext		
vext,all,,,mw
nsel,s,ext
sf,all,inf
csys,0
!------Freiheitsgradfestlegung für Wicklung-----
esel,s,type,,3
nsle
!csys,1
cp,1,curr,all
cp,2,emf,all
*get,knr1,node,0,num,max
!-------Elementkoordinaten--------
!esel,s,type,,3
!esel,all
!nsle,s
local,11,1
emod,all,esys,11
csys,0
!------Speisende Schaltung-----------------------
allsel
numcmp,node
*get,knr,node,,num,max
n,knr+1,-2.5*abm,0,0
n,knr+2,-2.5*abm,abm,0
n,knr+3,-2*abm,abm,0
n,knr+4,-1.5*abm,abm,0
n,knr+5,-1.5*abm,0,0
type,7	$real,7	  $e,knr+1,knr+2	!Kondensator
type,10	$real,10  $e,knr+2,knr+3	!Induktivität
type,8	$real,8	  $e,knr+3,knr+4	!Vorwiderstand
type,9	$real,9	  $e,knr+4,knr+5,knr1	!Koppelelement
eplot
!------Bezugspotential-------------------------
nsel,s,node,,knr+1
nsel,a,node,,knr+5
d,all,volt,0
allsel
eplot

!------LÖSUNG----------------------------------
/solu			
allsel			
antype,trans	
outres,all,all	
kbc,1		
endzeit=0.6e-3	!Berechnungszeit
endzeit=2.4e-3
schritte=2
schritte=24
nsub,schritte
time,endzeit
solve
!/eof
!------Zeitverläufe von Strom und Spannung---------------------------
/post26			
allsel
esel,s,type,,8		!Vorwiderstand
*get,wider,elem,,num,min
esol,2,wider,,smisc,2,Strom
plvar,2
/eof
!-------------------
/post1

!/psymb,cs,0
esel,s,type,,3
plve,jt
!------------------
/post1
vsel,s,,,4
eslv
/plopts,minm,0
/triad,off
plns,a,sum

vsel,s,,,1
vsel,a,,,2
eslv