!Simulation !aerosimul, 15.11.2017 !Lufttransformator zwei Wicklungen !harmonische Rechnung finish $/clear *abbr,aerosimul3,/input,aerosimul3,txt !-----Parameter----- abst=30e-3 !Abstand beider Spulen radau=120e-3 radau=80e-3 !Radius außen Betrachtungsraum strom=100 !Speisestrom frequenz=100e3 !Frequenz !---Primärspule--- radwiin1=30.5e-3 !Radius Spule prim innen radwiau1=35.5e-3 !Radius Spule prim außen howi1=12e-3 !Höhe Spule prim wizahl1=25 !Windungszahl asp1=(radwiau1-radwiin1)*howi1 !Fläche Primärwicklung strdi1=strom*wizahl1/asp1 !Stromdichte !---Sekundärspule--- radwiin2=14.5e-3 !Radius Spule sek innen radwiau2=15.5e-3 !Radius Spule sek außen howi2=7e-3 !Höhe Spule sek wizahl2=100 !Windungszahl !strdi2=strom*wizahl2/((radwiau2-radwiin2)*howi2) !Stromdichte asp2=(radwiau2-radwiin2)*howi2 !Fläche Sekudärwicklung /prep7 et,1,plane53,,,1 !Spule prim et,2,plane53,3,,1 !Spule sek et,3,plane53,,,1 !Luft mp,murx,1,1 !Spule prim !mp,rsvx,1,1/56e6 !Spule prim Kupfer mp,murx,2,1 !Spule sek mp,rsvx,2,1/56e6 !Spule sek Kupfer mp,murx,3,1 !Luft r,2,asp2,wizahl2 !real,2 !-----Geometrie----- rectng,radwiin1,radwiau1,abst/2,abst/2+howi1 !Spule prim rectng,radwiin2,radwiau2,-abst/2,-abst/2-howi2 !Spule sek pcirc,0,radau,-90,90 !Umgebungsluft aovlap,all numcmp,area !A1: obere Wicklung A2: untere Wicklung A3: Luft !/eof !-----Vernetzung----- !---Primärspule--- type,1 mat,1 esize,(radwiau1-radwiin1)/4 !Elementgröße amesh,1 !---Sekundärspule--- type,2 mat,2 real,2 esize,(radwiau2-radwiin2)/4 !Elementgröße amesh,2 !-----Luft----- type,3 mat,3 esize,howi1 !Elementgröße amesh,3 eplot !-----Randbedingung----- csys,1 !Zylinderkoordinaten nsel,s,loc,x,radau d,all,az,0 !Vektorpotential auf Null csys,0 !-----Lasten----- !---Sekundärspule--- esel,s,type,,2 nsle,s cp,1,emf,all cp,2,curr,all !---Primärspule--- esel,s,type,,1 nsle bfe,all,js,1,,,strdi1 allsel !-----Schaltung----- *get,n_max,node,,num,max !höchste Knotennummer *get,e_max,elem,,num,max !höchste Elementnummer esel,s,type,,2 !Koppelknoten nsle,s *get,n_kopp,node,,num,min abm=30e-3 !für Schaltungsdarstellung n,n_max+1,-abm,0 n,n_max+2,-abm,abm !-----Koppelelement----- !/icscale,,.05 !Darstellungsgröße et,4,circu124,5,0 !Elementvereinbarung r,4,1 !Symmetriefaktor type,4 !Elementaufruf real,4 e,n_max+1,n_max+2,n_kopp !Elementgenerierung !-----Schaltungsmasse----- d,n_max+1,volt,0 !-----Randbedingungen----- !allsel !lsel,s,ext !nsll,s,1 !d,all,az,0 !-----Lösung----- /solu allsel antype,harm harfre,frequenz solve !-----Auswertung----- /post1 set,1,1,1,0 plns,b,sum plf2d *get,u_real,elem,e_max+1,smisc,1 ! UR_R set,1,1,1,1 *get,u_imag,elem,e_max+1,smisc,1 ! UR_I u_betrag=sqrt(u_real*u_real+u_imag*u_imag)