!-------------------------------------------------------------!
!                                                             !
!                                                             !
! Projektname: waermeintegrierterReaktor                      !
! Dateiname:   main.txt                                       !
!                                                             !
!                                                             !
!-------------------------------------------------------------!
!-------------------------------------------------------------!
! Programmierer: 	Hendrik Dubbe                             !
! Datum:	     	28.02.2012                                !
!-------------------------------------------------------------!
! Datei-Inhalt:                                               !
! Testcase 5 - Dichtung (Punkte separiert) (1/4 Symmetrie)    ! 
!-------------------------------------------------------------!
!
! Definition des Dateinamens und des Titels 
/FILNAM,testCase4
/TITLE,testCase4

! Parameter
s    = 0.01     ! m

L_1  =  4*s     ! m
L_2  = 14*s     ! m
L_3  = 15*s     ! m
L_4  = 16*s     ! m
L_5  = 17*s     ! m
L_6  = 26*s     ! m
L_7  = 30*s     ! m

H_1  =  8*s     ! m
H_2  =  9*s     ! m
H_3  = 10*s     ! m
H_4  = 11*s     ! m
H_5  = 12*s     ! m

TOL  = 1e-5     ! m
ra   = 12*s     ! m
L    = 30*s     ! m
ri   = 8*s      ! m
F_Sp = 10       ! N
PP   = 2e5      ! N/m*m
rm = 0.5*(ra+ri)!m 
p_sp = PP*(rm*rm)/(ra*ra-ri*ri)+F_Sp/(3.14*(ra*ra-ri*ri)) ! N/m*m

/PREP7  

! Elementtyp 1
ET, 1, SOLID226, 11

! Elementtyp 2 (Kontaktelemente)
ET, 2, CONTA174, 1

! Elementtyp 3 (Target-Elemente)
ET, 3, TARGE170, 1

TREF,20
MPTEMP, 1,  20,  50, 

! Stoffeigenschaften - notwendige Kennwerte zur Ermittlung des 
! Temperaturfeldes:
!-------------------------------------------------------------!

! Dichte [kg/m^3]
MPDATA, DENS, 1, 1, 7800.0, 7788.3, 

! Wärmeleitfähigkeit [W/mK]
MPDATA, KXX, 1, 1, 15.00, 15.38, 

! spezifische Wärmekapazität [J/kgK]
MPDATA, C, 1, 1, 500, 500, 

! Werkstoffeigenschaften - notwendige Kennwerte zur Ermittlung 
! der thermischen Verformungen und Spannungen
!-------------------------------------------------------------!

! Querdehnzahl [mm/mm]
MPDATA, PRXY, 1, 1, 0.3, 0.3, 

! Wärmeausdehnungskoeffizient [1/K] 
MPDATA, ALPX, 1, 1,  0.00e-6, 16.17e-6,

! Elastizitätsmodul [N/m^2]
MPDATA, EX, 1, 1, 2.00e+11, 1.98e+11, 

! Reibungskoeffizient
MPDATA, MU, 1, 1, 0.1, 0.1

! Punkte
K,1001,  0,  0,  0, 
K,1002,  0, H_1,  0, 
K,1003,  0, H_2,  0, 
K,1004,  0, H_3,  0, 
K,1005,  0, H_4,  0, 
K,1006,  0, H_5,  0, 

K,1007, L_1, 0,  0, 
K,1008, L_1, H_1,  0, 
K,1009, L_1, H_2,  0, 
K,1010, L_1, H_3,  0, 
K,1011, L_1, H_4,  0, 
K,1012, L_1, H_5,  0, 

K,1013, L_2, H_1,  0, 
K,1014, L_2, H_2,  0, 
K,1015, L_2, H_3,  0, 
K,1016, L_2, H_4,  0, 
K,1017, L_2, H_5,  0, 

K,1018, L_3, H_1,  0, 
K,1019, L_3, H_2,  0, 
K,1020, L_3, H_3,  0, 
K,1021, L_3, H_4,  0, 
K,1022, L_3, H_5,  0,
 
K,1023, L_4, H_1,  0, 
K,1024, L_4, H_2,  0, 
K,1025, L_4, H_3,  0, 
K,1026, L_4, H_4,  0, 
K,1027, L_4, H_5,  0, 

K,1028, L_5, H_1,  0, 
K,1029, L_5, H_2,  0, 
K,1030, L_5, H_3,  0, 
K,1031, L_5, H_4,  0, 
K,1032, L_5, H_5,  0,
 
K,1033, L_6, 0,  0, 
K,1034, L_6, H_1,  0, 
K,1035, L_6, H_2,  0, 
K,1036, L_6, H_3,  0, 
K,1037, L_6, H_4,  0, 
K,1038, L_6, H_5,  0,
 
K,1039, L_7, 0,  0, 
K,1040, L_7, H_1,  0, 
K,1041, L_7, H_2,  0, 
K,1042, L_7, H_3,  0, 
K,1043, L_7, H_4,  0, 
K,1044, L_7, H_5,  0, 

K,1045, L_4, H_3,  0, 
K,1046, L_5, H_3,  0, 
K,1047, L_5, H_4,  0, 
K,1048, L_4, H_4,  0, 

K,1049, L_3, H_4,  0, 
K,1050, L_4, H_4,  0, 
K,1051, L_5, H_4,  0, 

! Flächen
A,1001,1007,1008,1002,     ! A1
A,1002,1008,1009,1003,     ! A2
A,1003,1009,1010,1004,     ! A3
A,1004,1010,1011,1005,     ! A4
A,1005,1011,1012,1006,     ! A5

A,1008,1013,1014,1009,     ! A6
A,1009,1014,1015,1010,     ! A7
A,1010,1015,1016,1011,     ! A8
A,1011,1016,1017,1012,     ! A9

A,1015,1020,1021,1016,     ! A10
A,1018,1023,1024,1019,     ! A11
A,1020,1025,1026,1021,     ! A12
A,1049,1050,1027,1022,     ! A13
A,1023,1028,1029,1024,     ! A14

! Dichtung
A,1045,1046,1047,1048,     ! A15

A,1026,1031,1032,1027,     ! A16

A,1028,1034,1035,1029,     ! A17
A,1029,1035,1036,1030,     ! A18
A,1030,1036,1037,1051,     ! A19
A,1051,1037,1038,1032,     ! A20

A,1033,1039,1040,1034,     ! A21
A,1034,1040,1041,1035,     ! A22
A,1035,1041,1042,1036,     ! A23
A,1036,1042,1043,1037,     ! A24
A,1037,1043,1044,1038,     ! A25

! Linienteilung
LSEL, S, LINE, ,   1, , 1, 0
LSEL, A, LINE, ,   3, , 1, 0
LSEL, A, LINE, ,   6, , 1, 0
LSEL, A, LINE, ,   9, , 1, 0
LSEL, A, LINE, ,  12, , 1, 0
LSEL, A, LINE, ,  15, , 1, 0

LSEL, A, LINE, ,  63, , 1, 0
LSEL, A, LINE, ,  65, , 1, 0
LSEL, A, LINE, ,  68, , 1, 0
LSEL, A, LINE, ,  70, , 1, 0
LSEL, A, LINE, ,  72, , 1, 0
LSEL, A, LINE, ,  74, , 1, 0
LESIZE,ALL,,,4,1,0

LSEL, S, LINE, ,   2, , 1, 0
LSEL, A, LINE, ,   4, , 1, 0
LSEL, A, LINE, ,  64, , 1, 0
LSEL, A, LINE, ,  66, , 1, 0
LESIZE,ALL,,,8,1,0

LSEL, S, LINE, ,  17, , 1, 0
LSEL, A, LINE, ,  19, , 1, 0
LSEL, A, LINE, ,  21, , 1, 0
LSEL, A, LINE, ,  23, , 1, 0
LSEL, A, LINE, ,  25, , 1, 0
LESIZE,ALL,,,10,1,0

LSEL, S, LINE, ,  51, , 1, 0
LSEL, A, LINE, ,  53, , 1, 0
LSEL, A, LINE, ,  55, , 1, 0
LSEL, A, LINE, ,  58, , 1, 0
LSEL, A, LINE, ,  61, , 1, 0
LESIZE,ALL,,,9,1,0

LSEL,ALL
LESIZE,ALL,,,1,1,0


! Rotationskörper erzeugen
ASEL, S, AREA, ,   1 , , 1, 0
ASEL, A, AREA, ,   2 , , 1, 0
ASEL, A, AREA, ,   3 , , 1, 0
ASEL, A, AREA, ,   4 , , 1, 0
ASEL, A, AREA, ,   5 , , 1, 0
ASEL, A, AREA, ,   6 , , 1, 0
ASEL, A, AREA, ,   7 , , 1, 0
ASEL, A, AREA, ,   8 , , 1, 0
ASEL, A, AREA, ,   9 , , 1, 0
ASEL, A, AREA, ,   10 , , 1, 0
ASEL, A, AREA, ,   12 , , 1, 0
VROTAT,  ALL,  ,  ,  ,  ,  , 1001, 1039, -90, 1

ASEL, S, AREA, ,   15 , , 1, 0
VROTAT,  ALL,  ,  ,  ,  ,  , 1001, 1039, -90, 1

ASEL, S, AREA, ,   11 , , 1, 0
ASEL, A, AREA, ,   13 , , 1, 0
ASEL, A, AREA, ,   14 , , 1, 0
ASEL, A, AREA, ,   16 , , 1, 0
ASEL, A, AREA, ,   17 , , 1, 0
ASEL, A, AREA, ,   18 , , 1, 0
ASEL, A, AREA, ,   19 , , 1, 0
ASEL, A, AREA, ,   20 , , 1, 0
ASEL, A, AREA, ,   21 , , 1, 0
ASEL, A, AREA, ,   22 , , 1, 0
ASEL, A, AREA, ,   23 , , 1, 0
ASEL, A, AREA, ,   24 , , 1, 0
ASEL, A, AREA, ,   25 , , 1, 0
VROTAT,  ALL,  ,  ,  ,  ,  , 1001, 1039, -90, 1

! Linienteilung Umfang
LSEL,ALL
LSEL, U, LOC, Z, 0-TOL,0+TOL
LSEL, U, LOC, Y, 0-TOL,0+TOL
LESIZE,ALL,,,8,1,0

! Meshen
TYPE,1
MAT,1
ALLSEL,ALL
VSEL, S, VOLU, ,   1 , , 1, 0
VSEL, A, VOLU, ,   2 , , 1, 0
VSEL, A, VOLU, ,   3 , , 1, 0
VSEL, A, VOLU, ,   4 , , 1, 0
VSEL, A, VOLU, ,   5 , , 1, 0
VSEL, A, VOLU, ,   6 , , 1, 0
VSEL, A, VOLU, ,   7 , , 1, 0
VSEL, A, VOLU, ,   8 , , 1, 0
VSEL, A, VOLU, ,   9 , , 1, 0
VSEL, A, VOLU, ,  10 , , 1, 0
VSEL, A, VOLU, ,  11 , , 1, 0
VSWEEP,ALL

VSEL, S, VOLU, ,  12 , , 1, 0
VSWEEP,ALL

VSEL, S, VOLU, ,  13 , , 1, 0
VSEL, A, VOLU, ,  14 , , 1, 0
VSEL, A, VOLU, ,  15 , , 1, 0
VSEL, A, VOLU, ,  16 , , 1, 0
VSEL, A, VOLU, ,  17 , , 1, 0
VSEL, A, VOLU, ,  18 , , 1, 0
VSEL, A, VOLU, ,  19 , , 1, 0
VSEL, A, VOLU, ,  20 , , 1, 0
VSEL, A, VOLU, ,  21 , , 1, 0
VSEL, A, VOLU, ,  22 , , 1, 0
VSEL, A, VOLU, ,  23 , , 1, 0
VSEL, A, VOLU, ,  24 , , 1, 0
VSEL, A, VOLU, ,  25 , , 1, 0
VSWEEP,ALL


! Elemente abrundet darstellen
/efacet,4

! Temperatur
TUNIF,20

ASLV, S, 1 

! Kontakte definieren
REAL,3
TYPE,2  ! CONTA174
R,,14,25000        					 ! Kontaktkonduktivitaet - Wert fuer TCC
									 ! Realkonstante 14 - CONTA Elemente
ASEL, S, AREA, , 68, ,               ! Auswahl Fläche Contact-Element
NSLA, S, 1                           ! Auswahl der zu den Flächen gehörigen Knoten
ESURF, ALL, , 
TYPE,3  ! TARGE170
ASEL, S, AREA, , 101, ,               ! Auswahl Fläche Target-Element
NSLA, S, 1                           ! Auswahl der zu den Flächen gehörigen Knoten
ESURF, ALL, , 

REAL,4
TYPE,2  ! CONTA174
R,,14,25000        					 ! Kontaktkonduktivitaet - Wert fuer TCC
									 ! Realkonstante 14 - CONTA Elemente
ASEL, S, AREA, , 69, ,               ! Auswahl Fläche Contact-Element
NSLA, S, 1                           ! Auswahl der zu den Flächen gehörigen Knoten
ESURF, ALL, , 
TYPE,3  ! TARGE170
ASEL, S, AREA, , 86, ,               ! Auswahl Fläche Target-Element
NSLA, S, 1                           ! Auswahl der zu den Flächen gehörigen Knoten
ESURF, ALL, , 

REAL,5
TYPE,2  ! CONTA174
R,,14,25000        					 ! Kontaktkonduktivitaet - Wert fuer TCC
									 ! Realkonstante 14 - CONTA Elemente
ASEL, S, AREA, , 77, ,               ! Auswahl Fläche Contact-Element
NSLA, S, 1                           ! Auswahl der zu den Flächen gehörigen Knoten
ESURF, ALL, , 
TYPE,3  ! TARGE170
ASEL, S, AREA, , 65, ,               ! Auswahl Fläche Target-Element
NSLA, S, 1                           ! Auswahl der zu den Flächen gehörigen Knoten
ESURF, ALL, , 

REAL,6
R,,14,25000        					 ! Kontaktkonduktivitaet - Wert fuer TCC
									 ! Realkonstante 14 - CONTA Elemente
TYPE,2  ! CONTA174
ASEL, S, AREA, , 70, ,               ! Auswahl Fläche Contact-Element
NSLA, S, 1                           ! Auswahl der zu den Flächen gehörigen Knoten
ESURF, ALL, , 
TYPE,3  ! TARGE170
ASEL, S, AREA, , 64, ,               ! Auswahl Fläche Target-Element
NSLA, S, 1                           ! Auswahl der zu den Flächen gehörigen Knoten
ESURF, ALL, , 

R,,14,25000        					 ! Kontaktkonduktivitaet - Wert fuer TCC
									 ! Realkonstante 14 - CONTA Elemente

ALLSEL,ALL

! Symmetriebedingung
csys,0                              ! kartesisches Koordinatensystem
NSEL,S,LOC,Z,0-TOL,0+TOL            ! alle Knoten der Ebene z = 0 selektieren
DSYM,SYMM,Z,                        ! Symmetriebedingung setzen
NSEL,S,LOC,Y,0-TOL,0+TOL            ! alle Knoten der Ebene y = 0 selektieren
DSYM,SYMM,Y,                        ! Symmetriebedingung setzen

! Festlager im Ursprung
NSEL,S,LOC,X,0-TOL,0+TOL
NSEL,R,LOC,Y,0-TOL,0+TOL
NSEL,R,LOC,Z,0-TOL,0+TOL
D,ALL,UX,0                          ! Festlagerung in x-Richtung setzen
D,ALL,UY,0                          ! Festlagerung in y-Richtung setzen
D,ALL,UZ,0                          ! Festlagerung in z-Richtung setzen

! Innendruck
ASEL, S, AREA, ,  26 , , 1, 0
ASEL, A, AREA, ,  46 , , 1, 0
ASEL, A, AREA, ,  47 , , 1, 0
ASEL, A, AREA, ,  50 , , 1, 0
ASEL, A, AREA, ,  59 , , 1, 0

ASEL, A, AREA, ,  109 , , 1, 0
ASEL, A, AREA, ,   91 , , 1, 0
ASEL, A, AREA, ,   82 , , 1, 0
ASEL, A, AREA, ,   72 , , 1, 0
ASEL, A, AREA, ,   75 , , 1, 0
ASEL, A, AREA, ,   74 , , 1, 0
ASEL, A, AREA, ,   84 , , 1, 0
ASEL, A, AREA, ,   97 , , 1, 0
ASEL, A, AREA, ,  63 , , 1, 0
ASEL, A, AREA, ,  67 , , 1, 0
SFA, ALL, 1, PRES, PP

! Einspannung über Druck (Keramikrohre)
ASEL, S, AREA, ,  44 , , 1, 0
ASEL, A, AREA, ,  40 , , 1, 0
ASEL, A, AREA, ,  36 , , 1, 0
ASEL, A, AREA, ,  32 , , 1, 0

ASEL, A, AREA, ,  111 , , 1, 0
ASEL, A, AREA, ,  114 , , 1, 0
ASEL, A, AREA, ,  117 , , 1, 0
ASEL, A, AREA, ,  120 , , 1, 0
SFA, ALL, 1, PRES, p_sp

! Thermische RB außen
ASEL, S, AREA, ,  28 , , 1, 0
ASEL, A, AREA, ,  32 , , 1, 0
ASEL, A, AREA, ,  36 , , 1, 0
ASEL, A, AREA, ,  40 , , 1, 0
ASEL, A, AREA, ,  44 , , 1, 0

ASEL, A, AREA, ,  43 , , 1, 0
ASEL, A, AREA, ,  57 , , 1, 0
ASEL, A, AREA, ,  56 , , 1, 0
ASEL, A, AREA, ,  61 , , 1, 0
ASEL, A, AREA, ,  80 , , 1, 0
ASEL, A, AREA, ,  79 , , 1, 0
ASEL, A, AREA, ,  88 , , 1, 0
ASEL, A, AREA, ,  104 , , 1, 0
ASEL, A, AREA, ,  121 , , 1, 0
ASEL, A, AREA, ,  120 , , 1, 0
ASEL, A, AREA, ,  117 , , 1, 0
ASEL, A, AREA, ,  114 , , 1, 0
ASEL, A, AREA, ,  111 , , 1, 0
ASEL, A, AREA, ,  107 , , 1, 0
SFA, ALL, 1, CONV, 10,  20
SFA, ALL, 1, RAD,  0.95 , 20

! Thermische RB innen
ASEL, S, AREA, ,  26 , , 1, 0
ASEL, A, AREA, ,  46 , , 1, 0
ASEL, A, AREA, ,  47 , , 1, 0
ASEL, A, AREA, ,  50 , , 1, 0
ASEL, A, AREA, ,  59 , , 1, 0

ASEL, A, AREA, ,  109 , , 1, 0
ASEL, A, AREA, ,   91 , , 1, 0
ASEL, A, AREA, ,   82 , , 1, 0
ASEL, A, AREA, ,   72 , , 1, 0
ASEL, A, AREA, ,   75 , , 1, 0
ASEL, A, AREA, ,   74 , , 1, 0
ASEL, A, AREA, ,   84 , , 1, 0
ASEL, A, AREA, ,   97 , , 1, 0
ASEL, A, AREA, ,  63 , , 1, 0
ASEL, A, AREA, ,  67 , , 1, 0
SFA, ALL, 1, CONV, 20,  40
SFA, ALL, 1, RAD,  0.95 , 40


! Alles selektieren
ALLSEL,ALL
!
! Kräfte des Solid-Modells auf das FE-Modell übertragen
SBCTRA
!
!BCSOPTION,,MINIMUM
!
! Speichern
SAVE


/solu
! Neue statische Solution
ANTYPE,STATIC,NEW		
!
! Lösen	
SOLVE
!
! Speichern
SAVE
!