! Zwei sich beruehrende Wuerfel, die anfangs ein unterschiedliches Temp.-Niveau besitzen,
! was sich aber dann 5 min ausgleichen soll.
! Alle Temperaturen jeden Zeitschrittes sollen von jedem Knoten gespeichert werden.


  fini
  /clear
  
  t0=1e-3
  tend=5*60
  nstep=60
  dt=(tend-t0)/(nstep)


! ### PREPROCESSOR ################################
 
 /prep7

! --- geometrie ----------------------------------

! Wuerfel A:
  x0=0
  x1=1
  y0=0
  y1=1
  z0=0
  z1=1
  
  block,x0,x1,y0,y1,z0,z1

! Wuerfel B:
  x2=1
  x3=2
  y2=1
  y3=0
  z2=1
  z3=0
  
  block,x2,x3,y2,y3,z2,z3
 

! --- Elementtyp + Vernetzen ------------------------
  et,1,70
  mp,ex,1,210e3
  mp,kxx,1,300
  mp,dens,1,7.850
  mp,c,1,500

  esize,z1/10
  mshape,1
  vmesh,all

  /pnum,kpoi,0
  /pnum,line,0
  /pnum,area,1
  /pnum,volu,0


! --- Kontakt: area 6 + area 11: --------------------
   
  et,2,170
  et,3,174

  r,3,,,1,0.1,0
  keyopt,2,2,0
  keyopt,2,3,0
  keyopt,2,4,0
  keyopt,2,5,1

  keyopt,3,1,2
  keyopt,3,2,2
  keyopt,3,4,2
  keyopt,3,5,0
  keyopt,3,7,0
  keyopt,3,8,0
  keyopt,3,9,1
  keyopt,3,10,2
  keyopt,3,11,0
  keyopt,3,12,6

!Target area
  asel,s,,,6
  nsla,s,1
  esln,s,0
  type,2
  esurf
  
!Contact area
  asel,s,,,11
  nsla,s,1
  esln,0
  type,3
  esurf
  
  alls

! --- Selektieren der Knoten: ----------------------

   nsel,s,loc,x,x0,x1,z1/10
   cm,volume_1,node
   alls

   nsel,s,loc,x,x2,x3,z1/10
   cm,volume_2,node
   alls
  

! ### SOLUTION #####################################

   /solu
  
! Anfangstemperaturen

 	  antype,trans,new
  	  timint,off,therm
  	  d,volume_1,temp,20+273.15
  	  d,volume_2,temp,80+273.15
 	  time,t0   
 	  kbc,1
 	  alls
 	  solve

 
!*** Berechnung der Temperatur von t0 bis tend
 
 *get,ndnbr,node,0,count		! Initialisierung Speichervariable
 *dim,tpvec,array,ndnbr,1		! für individ. Knotentemperatur
 *dim,knoten,array,ndnbr,1		! Knotenvektor
 *dim,tvec,array,nstep+1,1		! Zeitvektor
 *dim,tpvert,array,ndnbr,nstep+1	! Vektor Temperaturverteilung

     *do,i,1,nstep+1,1
        *if,i,eq,1,then,
           tvec(i,1)=t0
             *elseif,i,ge,2,then,
           	tvec(i,1)=tvec(1,1)+(i-1)*dt
            *endif
    *enddo
   

     *do,j,1,nstep,1

	  timint,on,therm   
   	  time,tvec(j+1)		! Zeit aus Zeitvektor festgelegt
   	  deltim,dt
   	  autotst,off
   	  outres,all,all
  	  solve

     *enddo


 ! Berechnungsschleifen
   
	*do,lj,1,ndnbr,1
	  nsel,s,,,lj
          *get,tpval,node,lj,temp	! Ablesen der Temp aus der DB für Knotn lj
          *set,tpvec(lj,1),tpval+273.15	! Wertzuweisung
          *set,knoten(lj,1),lj
	  alls
	*enddo


     *do,j,1,nstep,1			! diskrete Zeiten auswählen
      *do,k,1,ndnbr,1			! alle Knoten durchlaufen
        nsel,s,,,k
       
          tpvert(0,j,1)=tvec(j,1)
          tpvert(k,0,1)=knoten(k,1)
          tpvert(j,k,1)=tpvec(j,k)
            
        alls
      *enddo
     *enddo


! ### POSTPROCESSOR ################################

   /post1
   plnsol,temp