!Vernetzung eines Dichtungsringes
!für Atzenpower
!Schätzing, 15.6.2009

finish
/clear

*abbr,netz3,/input,vernetzung3,inp

rad=625e-6
breitu=50e-6   !Breite unten
dick=20e-6     !Dicke
breito=30e-6   !Breite oben
radrund=5e-6   !Radius Rundung
!--------------------------------------
/prep7
et,1,mesh200,6   !Hilfselemente für Flächenvernetzung
et,2,mesh200,10  !Volumenelemente: muss konkretisiert werden
!---halbe-Querschnittsfläche-----
local,11,0,rad   !Lokales Koordinatensystem
wpcsys,,11
k,1,breitu/2,0,0
k,2,breito/2,dick,0
k,3,0,dick,0,0
k,4,0,0,0
!--------Begrenzungslinien---
l,1,2
l,2,3
l,3,4
l,4,1
!----Abrundung---------
lfillt,1,2,radrund
!--------Keypoints für Flächenunterteilung--
k,7,breito/2-radrund/2,dick-radrund/2
k,8,0,dick-radrund/2
k,9,breitu/2-radrund/2,0

!-------Flächengenerierung-----
/pnum,kpoi,1
kplot
l,7,5
l,6,7
/pnum,line,1
lplot
al,6,5,7  !A1
/pnum,kpoi,1
kplot
a,7,6,3,8,  !A2
a,1,5,7,9   !A3
a,9,7,8,4   !A4
btol,1e-6   !wegen der kleinen Abmessungen
aglue,all
/pnum,area,1
aplot

!------Vernetzung-Fläche------
mshkey,1
mshape,0,2D
type,1
esize,dick/8
amesh,1
amesh,2
amesh,3
amesh,4

!-----Symmetrie nutzen----
arsym,x,all
nummrg,kpoi,1e-7
nummrg,node,1e-7
!/eof
!------Rotation-um-Z-Achse-------
wpcsys,0
csys,0
type,2
esize,,60
k,20,0,0,0   !für Rotation
k,21,0,10e-6,0
vrotat,all,,,,,,20,21,90
aclear,all   !Löschen der Hilfselemente
eplot