! rk, 2006-12-10
! Holz-Stahlblechverbindung

/clear,nostart
/filname,Test,1
/title, Stahlblech mit Dicke t=10mm, bündig
!/cmap,'color','cmap'


/VIEW,1,-1
/ANG,1, 90,ZS,1 
/ANG,1,-60,YS,1
/ANG,1, 15,XS,1
/auto,1
/pnum,volu,1	! Controls entity numbering/coloring on plots



! Parameter
! ---------
! Holz
*set,LAENGE,0.165	! Angaben in SI-Einheiten
*set,BREITE,0.02
*set,HOEHE ,0.05

! Blech
*set,DICKE ,0.005
*set,TIEFE ,0.110



! Geometrie des Modells
! ---------------------
/prep7
   ! Volumen Holz
   block,0,BREITE,0,LAENGE,0,HOEHE

   ! Volumen Blech
   block,0,BREITE,0.055,LAENGE,0,DICKE
   vsbv,1,2
   block,0,BREITE,0.055,LAENGE,0,DICKE

   vglue,all
finish



! Materialeigenschaften
! ---------------------
/prep7

   ! Materialeigenschaften von Holz
   mp,kxx,1,1.7
   mp,c,1,1500
   mp,dens,1,480
   mp,emis,1,.56
   ! Materialeigenschaften von Stahl
   mp,kxx,2,40
   mp,c,2,500
   mp,dens,2,7800
   mp,emis,2,.56
finish


! Vernetzung des Modells
! ----------------------
/prep7

   et,1,solid70	
   et,2,surf152,,,,1,1
   keyopt,2,8,2		! Oberflaechenelemet mit(!) Konvektion
   keyopt,2,9,1		! Oberflaechenelemet mit(!) Strahlung
   r,2,1,5.67E-8		! Formfaktor = 1, Stefan-Boltzmann-Konstante
   mshape,0,3D
   mshkey,1

   ! Blech
   type,1
   mat ,2

   vsel,,volu,,1
   vplot
   ldiv,1,,,4
   ldiv,3,,,4
   numcmp,line

   lesize, 1,0.001
   lesize,16,0.002
   lesize,17,0.0025
   lesize,18,0.0025

   lesize,21,0.001
   lesize,20,0.002
   lesize,19,0.0025
   lesize, 3,0.0025

   lesize,11,,,2
   lesize,23,0.001

   vsweep,all
   allsel,all


   ! Holz
   type,1
   mat ,1

   vsel,,volu,,2
   vplot
   ldiv,28,,,3
   ldiv,29,,,3

   ldiv, 6,,,6
   ldiv, 8,,,6

   ldiv, 9,,,10
   ldiv,10,,,10
   numcmp,line

   lesize, 8,0.001
   lesize,39,0.002
   lesize,40,0.0025
   lesize,41,0.0025
   lesize,42,0.0025
   lesize,43,0.0025

   lesize,38,0.001
   lesize,37,0.002
   lesize,36,0.0025
   lesize,35,0.0025
   lesize,34,0.0025
   lesize, 6,0.0025

   lesize,33,0.001
   lesize,29,0.0025
   lesize,32,0.0025

   lesize,31,0.001
   lesize,28,0.0025
   lesize,30,0.0025

   lesize, 4,,,2
   lesize,26,0.0025

   lesize,22,0.001
   lesize,23,0.001
   lesize,10,0.001
   lesize,52,0.001

   lesize,49,0.001
   lesize,50,0.001
   lesize,51,0.001

   lesize,53,0.001
   lesize,54,0.001
   lesize,55,0.001

   lesize,44,0.0025
   lesize,45,0.0025
   lesize,46,0.0025
   lesize,47,0.0025
   lesize,48,0.0025
   lesize, 9,0.0025

   lesize,56,0.0025
   lesize,57,0.0025
   lesize,58,0.0025
   lesize,59,0.0025
   lesize,60,0.0025
   lesize,61,0.0025

   vsweep,all
   allsel,all   

   ! Oberflaechenelemente
   nsel,s,loc,y,.165   
   esln,s
   type,2
   real,2
   n,10000,BREITE/2,LAENGE*1.2,HOEHE/2	! Umgebungsknoten
   esurf,10000 					! GENERATE SURF152 ELEMENTS
   allsel,all 
finish


! Loesung des Gleichungssystems
/solu
   antype,trans
   autots,1
   nsubst,100,200,20
   IC,all,temp,273!tunif,273			! Assigns a uniform temperature to all nodes

   ! Startbedingungen (Lastschritt 1)
   nsel,s,loc,x,BREITE
   nsel,a,loc,y,LAENGE
   nsel,a,loc,z,HOEHE 
   cm,AUSSEN,node

   *set,ALPHA,25			! Waermeuebergangskoeffizient
   sf,all,conv,ALPHA,273	! Konvektion

   allsel,all
   d,10000,temp,273
   time,1e-6
   solve
finish