!linearer Widerstand
!als Vorübung für Varistor
!13.10.08
!Schätzing
!-----------------------------------------

finish
/clear

*abbr,R_lin,/input,R_lin.txt  

radpil=2.5e-3      !Radius Pille  													
dipil=2e-3         !Dicke Pille

		
!-------------------
/prep7                              	 	 
!---Geometrie---
!a=0.1e-3
!lang=20*a
!/triad,off
!/plopts,info,0
!Cylind,0,radpil,0,dipil
!cylind,0,a,-lang,0
!cylind,0,a,dipil,dipil+lang
!/color,volu,lgra,1
!/color,volu,dgra,2,3
!/color,area,dgra,all
!/color,area,lgra,4
!aplot

!/eof

et,1,plane67,,,1  !Rotationssysmmetrisch															

!-----Keramik-Material--------
mp,kxx,1,2.2   !	             
mp,dens,1,2400														 
mp,c,1,550     !
mp,rsvx,1,1.2    !willkürlicher Wert


!-----Geometrie-----
rectng,0,radpil,0,dipil

!-----Vernetzung-----
esize,,3
amesh,1
!-----Randbedingungen-elektrisch-----
nsel,s,loc,y,dipil
cp,1,volt,all   		
n_ein=node(0,dipil,0)   !Einspeiseknoten
nsel,s,loc,y,0
d,all,volt,0

!-----Randbedingungen-thermisch-----
nsel,s,loc,x,radpil
nsel,a,loc,y,0
nsel,a,loc,y,dipil

sf,all,conv,10,20
allsel
!ic,all,temp,20	

!/eof

!-----Speisung mit Spannungsquelle------
u1=80                  !Versorgungsspannung
r1=0.1e-3             !Vorwiderstand
et,3,circu124,0       !Widerstand
r,3,r1
abst=2e-3             !für Knotendarstellung
*get,n_max,node,,num,max
n,n_max+1,0,dipil+abst
type,3  $real,3
e,n_max+1,n_ein       !Vorwiderstand
*get,e_max,elem,,num,max
!-------------------
tunif,20
!/eof

!-----Lösung-----
allsel
/solu
zeitende=20e-3
pi=4*atan(1)
antype,trans
outres,all,all
kbc,1
schritte=20
*do,z,0,schritte
time,zeitende/schritte*z+1e-12
arg=sin(2*pi*z/schritte)
d,n_max+1,volt,u1*arg     !Einspeisung
solve																			!Lösung anfordern
*enddo
!-----Auswertung-----

/post1
!/plopts,info,0
/cont,,9,20,,33
plnsol,temp


!/eof

/post26
nsol,2,n_ein,temp,,Temperatur
esol,3,e_max,,smisc,2,Strom
!----Darstellung des Widerstandes-----
nsol,4,n_ein,volt,,Potential    !Spannung
!quot,5,4,3,,Widerst  !Spannung/Strom


/grtyp,2
/yrange,20,32,1     !Temperatur
/yrange,-1,1,2      !Strom
/yrange,-100,100,3  !Spannung
/gmarker,2,1        !Strom mit Dreiecken
/axlab,x,Zeit in s
plvar,2,3,4