Example 331: Transient heat flow in a plate

This example uses calculates the transient temperature distribution inside a plate using FEM. An online analytic solution is available for the problem.


The comfile run by this example is as follows:

#Example_331 Transient heat flow in a plate                          
                                                                    
fem                            #set the FEM environment               
fem define nodes;r;trplate;example #6 nodes to describe a plate
fem define window;c            #calculate the window size             
fem define bases;r;;example        #read the bilinear basis function
fem define elements;r;;example     #defines two elements
fem draw lines                 #make the line segments visible        
fem draw elements              #make the element numbers visible       
fem draw nodes                 #make the node numbers visible          
fem define equation;r;;example      #we are dealing with a TIME DEPENDENT
#                              !  ADVECTION-DIFFUSION equation.  Use 
#                              !  default options for other choices
fem define material;r;;example      #see manual for advection-diffusion 
#                              !  equation.  If in doubt, entering ? at
#                              !  the question prompt will give further
#                              !  explanation.  We wish to set k 
#                              !  (storage coefficient) and each Di
#                              !  (diffusivity) to 1 and the rest to 0 
fem define initial;r;;example       #set the boundary conditions as in 
#                              !  Example_311 (with zero as the initial 
#                              !  condition)                           
fem define solve;r;;example         #use a linear time integration algorithm
#                              !  with automatic stepping and an 
#                              !  initial time increment of 0.001s. Use
#                              !  default options for other choices  
fem list equation              #list the equation parameters      
fem list materials             #list the material parameters           
fem list initial               #list the initial and boundary 
#                              !  conditions                          
fem list solve                 #list the solution parameters.  The time
#                              !  integration is (by default) from t=0 
#                              !  to t=1s. The time increment is 
#                              !  initially 0.0001s  but this increases
#                              !  when the solution is changing slowly
fem solve                      #performs the time integration          
fem display history            #shows the temperature time history at 
#                              !  all nodes                           
#fem display section at 1,2,3  !shows the changing temperature profile 
#                              !yet to be updated
fem draw contour dependent value 0..10:1  #draws contours 0 to 10 in 
#                              !  steps of 1 degree                    
fem draw field                 #this will give a shading of the         
#                              !  temperature field                  

Files used by this example are:

Name                Modified     Size

example_331.com 13-Feb-2001 3.0k analytic_331.c 10-Apr-2000 5.7k analytic_331.cgi 10-Apr-2000 34k analytic_331.html 10-Apr-2000 1.2k makefile 10-Apr-2000 427 plate.binhis 10-Apr-2000 129k plate.com 10-Apr-2000 145 plate.ipbase 10-Apr-2000 1.1k plate.ipelem 10-Apr-2000 15k plate.ipequa 26-May-2003 1.2k plate.ipequa.old 18-Jun-2002 1.2k plate.iphist 10-Apr-2000 64k plate.ipinit 10-Apr-2000 1.4k plate.ipmate 07-Apr-2003 2.6k plate.ipmate.old 10-Apr-2000 2.0k plate.ipmesh 10-Apr-2000 743 plate.ipnode 10-Apr-2000 12k plate.ipsolv 13-Apr-2007 1.7k trplate.history 10-Apr-2000 128k trplate.ipbase 10-Apr-2000 1.1k trplate.ipelem 10-Apr-2000 611 trplate.ipequa 26-May-2003 1.2k trplate.ipequa.old 18-Jun-2002 1.2k trplate.iphist 13-Feb-2001 219k trplate.ipinit 13-Feb-2001 1.1k trplate.ipmate 07-Apr-2003 2.5k trplate.ipmate.old 10-Apr-2000 1.9k trplate.ipnode 10-Apr-2000 1.3k trplate.ipsolv 13-Apr-2007 2.8k

Download the entire example:

Name                      Modified     Size

examples_3_33_331.tar.gz 14-Apr-2007 118k

Html last generated: Sun Mar 6 05:50:17 2016

Input last modified: Fri Apr 13 10:32:22 2007


CMISS Help / Examples / 3 / 33 / 331