!von Elektriker1
!bearbeitet Schätzing, 18.7.2010

fini
/cle
*abbr,elektriker,/input,elektriker,txt


a=0.15   !Breite
b=0.1    !Länge
eps1=50
eps2=1
eps2=0.001

eps0=8.856e-12
kappa2=1e-7      !angenommen für etwa gleiche Größe von R und C
freq=50
u_speis=20000
pi=4*atan(1)
!------analytische Rechnung-----
R2=b/(kappa2*a*1)
C2=eps0*eps2*a*1/b
C1=eps0*eps1*a*1/b
XC1=-1/(2*pi*freq*C1)
Z=sqrt(R2**2+XC1**2)
wi_a=atan(XC1/R2)*180/pi
i_b_a=u_speis/z

!/eof

/prep7

!-----------------Elementtyp und Materialdaten---------------

et,2,plane121

!Hauptisolation_Grundwickel_Spaltglimmer_Nutkeil
mp,perx,1,eps1
!IPS
mp,perx,2,eps2
mp,rsvx,2,1/kappa2  !technisch nicht sinnvoll
!Blechpaket

rect,0,a,0,b
rect,0,a,b,2*b
aglue,all

asel,s,area,,1
aatt,1
asel,s,area,,3
aatt,2
allsel
esize,0.01
mshkey,1
mshape,0,2d
amesh,all




!lsel,s,line,,1
!nsll,s,1
nsel,s,loc,y,2*b
cp,1,volt,all
*get,n2,node,,num,min
d,n2,volt,u_speis
alls


!lsel,s,line,,7
!nsll,s,1
nsel,s,loc,y,0
d,all,volt,0
allsel


!-------------------------
/solu
antype,harm
harfrq,freq
solve
!----------------------
/post1
!--------real-------
set,1,1

!plns,volt
!/eof

esel,s,mat,,2
!-----------------
u_oben_real=volt(node(0,2*b,0))
u_mitte_real=volt(node(0,b,0))
u_unten_real=volt(node(0,9,0))

ples,jt,y
nsel,s,loc,y,2*b
esln,s,0
*get,e_mess,elem,,num,min
*get,j_real,elem,e_mess,jt,y
!da alle Elemente gleich groß sind
!und alle gleiche Stromdichte haben:
i_real=j_real*a*1


!--------imaginär---
allsel
set,1,1,,1
!-----------------
u_oben_imag=volt(node(0,2*b,0))
u_mitte_imag=volt(node(0,b,0))
u_unten_imag=volt(node(0,0,0))

esel,s,mat,,2
ples,jt,y
nsel,s,loc,y,2*b
esln,s,0
*get,e_mess,elem,,num,min
*get,j_imag,elem,e_mess,jt,y
i_real=i_real
i_imag=j_imag*a*1
i_betrag=sqrt(i_real**2+i_imag**2)

!------------Winkel-------
wi=atan(i_imag/i_real)*180/pi



/eof




set,1,1
nsel,s,node,,3
fsum


set,1,1,,1
nsel,s,node,,3
fsum