!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Durchfuehrung einer thermischen Simulation mit anschließender struktur-     !!!
!!! mechanischen Berechnung der Tragfaehigkeit                                  !!!
!!! Brandsimulation                                                             !!!
!!! Jan  17/10/2013                                                             !!!
!!! PLANE55-ELEMENTE                                                            !!!
!!! SURF151-ELEMENTE                                                            !!!
!!! EINHEITEN mm, s, t, N                                                       !!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


!capone

!LEEREN DER DATENBASIS

  finish
  /clear,start



!DEFINITION DER PARAMETER
   LNGA=4800         !TRAEGERLAENGE
   LNGB=290-14       
   LNGC=300          !STEGBREITE
   LNGD=14           !FLANSCHDICKE
   LNGE=8.5          !STEGDICKE
   LS=100            !Streckenlast
   LAST=LS*4800      !EINHEITSLAST - RESULTIERENDE AUS DER FLAECHENLAST

   N_EL_L=50         !ANZAHL DER ELEMTE ÜBER DIE TRAEGERLAENGE
   N_EL_F=5          !ANZAHL DER ELEMTE ÜBER DEN FLANSCH
   N_EL_S=6          !ANZAHL DER ELEMTE ÜBER DEN STEG


   N_EW=2            !ANZAHL DER EIGENWERTE

!Materialparameter

   QDehn=0.3         !QUERDEHNZAHL               EINHEIT PRÜFEN
   RO1=7.9E3         !Dichte Stahl

! Temperaturdaten
  

!-----------------------------------------------------------------------------------

! Preprozessor starten
  /prep7

!UMSTELLUNG VON KELVIN ZU CELSIUS
  TOFFST,274                 !Sprich eine Erhoehung der Temperatur um 274Grad

! DEFINITION DER KEYPOINTS

  k,1,0,0,0
  k,2,-LNGC/2,0,0
  k,3,LNGC/2,0,0
  k,4,0,LNGB/2,0,0
  k,5,0,LNGB,0
  k,6,-LNGC/2,LNGB,0
  k,7,LNGC/2,LNGB,0

  k,11,0,0,LNGA 
  k,12,-LNGC/2,0,LNGA 
  k,13,LNGC/2,0,LNGA 
  k,14,0,LNGB/2,LNGA 
  k,15,0,LNGB,LNGA 
  k,16,-LNGC/2,LNGB,LNGA 
  k,17,LNGC/2,LNGB,LNGA 

! Definition der Linien 
  
 l,1,2     !Untergurt Stabanfang 
 l,1,3     !Untergurt Stabanfang
 l,1,4     !Steg      Stabanfang
 l,4,5     !Steg      Stabanfang
 l,5,6     !Obergurt  Stabanfang 
 l,5,7     !Obergurt  Stabanfang

 l,11,12   !Untergurt Stabende 
 l,11,13   !Untergurt Stabende
 l,11,14   !Steg      Stabende
 l,14,15   !Steg      Stabende
 l,15,16   !Obergurt  Stabende 
 l,15,17   !Obergurt  Stabende

 l,1,11    
 l,2,12   
 l,3,13
 l,4,14
 l,5,15
 l,6,16
 l,7,17

! Unterteilen der Linien
 
  lsel,s,line,,9,10,1
  lsel,a,line,,3,4,1  
  lesize,all,,,N_EL_S

  lsel,all
  lsel,s,line,,1,2,1
  lsel,a,line,,5,8,1
  lsel,a,line,,11,12,1  
  lesize,all,,,N_EL_F

  lsel,all
  lsel,s,line,,13,19,1
  lesize,all,,,N_EL_L
  lsel,all
  
! Definition der Flaechen
 
 al,11,18,5,17
 al,12,17,6,19
 al,10,17,4,16
 al,9,16,3,13
 al,8,13,2,15
 al,7,13,1,14

 aglue,all          ! FLAECHEN VERKLEBEN
   


! Definition der Elemente

  et,1,shell131
  et,2,surf152
  keyopt,2,5,1

 ! R,2,1.0,5.67e-8

  r,1,LNGD
  r,2,LNGD
  r,3,LNGE


!----------------------------------------------------
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!NUR PROVISORISCH ZUR BESSEREN UEBERSICHT!!!!!!!!

! Materialeigenschaften, Elementtyp  deklarieren   
 
!WAERMELEITFAEHIGKEIT
  MP,KXX,1,60
 

!WAERMEKAPAZITAET          

 MP,C,1,500


!WAERMEAUSDEHNUNGSKOEFIZIENT
 MP,THSX,1,0.01
 MP,THSY,1,0.01
 MP,THSZ,1,0.01
  
!E-MODUL                        
 
MP,EX,1,21000
MP,PRXY,1,0.3              !QUERDEHNUNG IM ELASTISCHEN BEREICH 
  
  
           



!DICHTE
  mp,dens,1,RO1

!EMISSIONSVERMOEGEN
  mp,emis,1,.70

!-------------------------------------------------------------------------
!VERNETZUNG DES MODELLS

  mat,1
  type,1

  real,1
  asel,s,area,,1
  asel,a,area,,2
  amesh,all
  asel,all


  real,2
  asel,s,area,,5,6,1
  amesh,all
  asel,all


  real,3
  asel,s,area,,3,4,1
  amesh,all
  asel,all



!-------------------------------------------------------------------------
!Temperaturbelastung aufbringen

 type,2
 real,2

 /pnum,area,1
 !aplot
  asel,s,,,1,6
  asll,s,1
  n,5555,0,LNGB/2,LNGA/2

  esurf,5555
  alls

  /pnum,type,1
  eplot

  /shrink,.2     !verkleinerte Darstellung der Elemente um zu Überprüfen ob 
                 !die SURF152 Elemente auf die Oberflaeche des Bauteils aufgetragen wurden 
  /shrink,1      !Standarddarstellung wieder festlegen



!------------------------------------------------------------------------------------------

! Preprozessor beenden
    finish


! LINEARE STATISCHE BERECHNUNG 
! Lösungsteil starten
    /solu


!System lösen
     solve
!------------------------------------------------------------------------------------------

 antype,trans             !Transiente Berechnungsmethode auswaehlen
 outres,all,all
 autots,on
 tunif,273                !Starttemperatur vor dem ersten LOADSTEP festlegen 
 


!------------------------------------------------------------------------------------------
!Konvektion fuer den Rand angeben und den Waermeuebergangkoefizienten fuer festlegen

  asel,s,,,1,6
  asll,s,1

 alpha=250
 sfe,all,1,conv,1,alpha       !film coefficients
 sfe,all,1,conv,2,273         !Bulk temperatur
 d,all,tbot,273
 ic,all,tbot,273
 cm,E_BOTTOM,elem
 alls                         !Selektierung des Randes wird aufgehoben

  asel,s,,,1,6
  asll,s,1

 alpha=250
 sfe,all,1,conv,1,alpha       !film coefficients
 sfe,all,1,conv,2,273         !Bulk temperatur
 d,all,ttop,273
 ic,all,ttop,273
 cm,E_TOP,elem



 d,5555,TEMP,273              !Temperatur der Strahlungsumgebung generieren
 TIME,1e-6

!System lösen
 solve

!LOADSTEP 2
 
 sfe,all,1,conv,1,alpha
 sfe,E_TOP,1,conv,2,273+800
 sfe,E_BOTTOM,1,conv,2,273+800
 
 d,5555,Temp,273+800
 kbc,0
 TIME,10*60
 NSUBST,60,60,20
 
!System lösen
 solve


!------------------------------------------------------------------------------------------
!Postprozessor starten
/post1

 plns,temp                  !Grapische Darstellung der Temperaturverteilung