!copyright @ Wolfgang Schätzing und UglyJoe(2)
!Kapazitätsermittlung über Verschiebungsfluss
!am Halbmodell

finish
/clear
*abbr,uglcad3,/input,uglcad3,txt

!-----Parameter aus der Aufgabenstellung-----
h=20e-3     !BasisPlattenDicke
w=4e-3      !BandBreite
s=12e-3     !Bänderabstand
t=0.1e-3    !BandDicke
pot2=+1     !Potenzial rechter Leiter
pot1=-1     !Potenzial linker Leiter
ladung=1e-9 !Ladung auf eine Elektrode
epsr2=9.8   !Dielektrizitätszahl
tand=1e3    !Verlustfaktor

!-----Parameter gewählt-------
b=(2*w+s)*4          !BasisPlattenBreite
b=b*2
radin=b/2+10e-3      !Radius innen für Berechnungsbereich
radau=radin+60e-3    !Radius außen für Berechnungsbereich


!------------------------------
/prep7
et,1,plane121
mp,perx,1,1      !Permittivität Luft
mp,perx,2,epsr2  !Permittivität Basisplatte
mp,lsst,2,tand



!----Geometrie--------
pcirc,0,radin,-90, 90 !Berechnungsfläche
rectng,+s/2,+s/2+w,0,t  !rechter Bandleiter
AsbA,1,2 
rectng,0,b/2,0,-h  !Basisplatte
aovlap,all
numcmp,area
aplot
!A1: Basisplatte  
!A2: Luft


!************
!* Leiter
!************
LESize,4,,,4
LESize,5,,,1
LESize,6,,,4
LESize,7,,,1



!************
!* Isolator
!************
LESize,8,,,20
LESize,9,,,5
LESize,13,,,20
LESize,11,,,5
LESize,12,,,5



!************
!* Luft
!************
LESize,1,,,20			
LESize,15,,,12
LESize,14,,,10


!-----Vernetzung-----
Mat,1
AMesh,2

Mat,2
AMesh,1,


!-----Symmetrierandbedingung---------

!Symmetriebedingung Knoten Auswahl des 0 Potentials (gesammte Symmetrieachse)

AllSel
nsel,s,loc,x,0
d,all,volt,0

!-----Ladungsvorgabe---------
lsel,s,,,4,7		!Leiter
nsll,s,1
cp,1,volt,all
*get,n1,node,,num,min
f,n1,chrg,ladung
			
AllSel			
!-------Lösung-----
/solu			
solve
!------Auswertung--------
/post1			
plnsol,volt,,1	
	
!Kapazitätsbestimmung über die Ladung
!C=Q/U  

kap_HM=Ladung/volt(n1)
!Das ist die Kapazität des Halbmodells
!beim Gesamtmodell sind zwei davon in Reihe
kap_GM=kap_HM/2
	
!kapazität=53,47 pF/m