This example uses calculates the transient temperature distribution inside a plate using FEM. An online analytic solution is available for the problem.
#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
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
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