!Durchschnittstemperatur
!
!13.10.08
!Schätzing
!-----------------------------------------

finish
/clear

*abbr,d_temp,/input,d_temp.txt  

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

		
!-------------------
/prep7                              	 	 
!---Geometrie---

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,1000,20
allsel
ic,all,temp,20	


!-----Speisung mit Spannungsquelle------
u1=12                 !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
!-------------------


!-----Lösung-----
allsel
/solu
antype,trans
outres,all,all
autots,on
nsubst,10,20,5
zeitende=20
kbc,1
d,n_max+1,volt,u1       !Einspeisung
time,zeitende
solve																			!Lösung anfordern
!-----Auswertung-----

/post1
esel,s,type,,1
plnsol,temp



 
!----Durchschnittstemperatur-----------
etable,temperatur,temp
plet,temperatur    !nur zur Kontrolle
etable,volumen,volu
!plet,volumen
smult,tempvol,temperatur,volumen
!plet,tempvol
ssum
*get,sumvol,ssum,,item,volumen
*get,sumtevol,ssum,,item,tempvol
Durchtemp=sumtevol/sumvol