!copyright @ Wolfgang Schätzing und UglyJoe(2)

finish
/clear
*abbr,uglcad,/input,uglcad,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
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

/PNUM,AREA,1		!PlotCtrls>Numbering>AREA: on <OK>
			!Flächennummerierung einschalten
APlot			!U.M.>Plot>Areas

!A1: Basisplatte  
!A2: Luft





!ESize,nd


!Linieneinteilung
!Preprocessor>Meshing>SizeControls>ManualSize>Picked Lines



!************
!* 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			!U.M.>Select>Everything 
			!Selektiert alles
nsel,s,loc,x,0
!NPlot			!U.M.>Plot>Nodes
d,all,volt,0

!-----Potentialvorgabe---------
lsel,s,,,4,7		!Leiter
nsll,s,1
D,all,volt,pot2		

			!3.	Lösung
AllSel			!U.M.>Select>Everything !Selektiert alles
!-------Lösung-----
/solu			!Aktivieren des Prozessors SOLUTION
			!M.M.>Solution>
solve			!	"	Solve>Current LS

			!4.	Postprozessor
!------Auswertung--------
/post1			!Startet den allgemeinen Postprozessor POST1

			!4.1. Ausgabe der Resultate
plnsol,volt,,1		!M.M.>Gen. Post.>Plot Results>Contour Plot>Nod.Sol.>


			!4.2 	Kapazitätsbestimmung
!4.2.1.	Kapazitätsbestimmung über die Energie (Elementtabelle)

!W=C/2*U*U   !Gesamtanordnung

AllSel					!U.M.>Select>Everything !Selektiert alles
ETable,energie,sene			!M.M.>Gen. Post.>Element Table>Define Table> 
					!Add; Lab: SENE; Energy: SENE; OK; Close;
					!Füllt eine Tabelle mit Elementwerten  der "Steifigkeits"energie SENE
					!für weitere Berechnungen aus
					!Gibt den Inhalt SENE der Elementtabelle aus
!PlEtab,energie,avg  			!M.M.>Gen. Post.>Element Table>Plot Elem Table> 
					!Itlab: SENE; Avglab: Yes <OK>
					!Stellt den Inhalt SENE der Elementtabelle grafisch durch 
					!Mittelwertsbildung dar 
SSum					!M.M.>Gen. Post.>Element Table>Sum of each item> OK;
					!Ermittelt und gibt die Summe aller Eintragungen der Elementtab.aus
*get,sumenergie,ssum,,item,energie	!U.M.>Parameters>Get scalar data> 
					!results data: elem table sums; par: w; OK
					!Auslesen des Wertes von SSUM und speichern als Parameter W
!sumenergie=2*sumenergie			!Achtung! 
					!Die Gesamtenergie ist der 2-fache Wert der Halbanordnung.

!C=2*sumenergie/((pot2-pot1)*(pot2-pot1))	
!kapazität=53,43 pF/m  (53,31)
!auf die halbierte Geometrie bezogen, besser so
C=sumenergie/(pot2*pot2)