Example 3141: Steady state heat distibution in a plate calculated using the laplace equation

This example uses the laplace equation to calculate the temperature distribution inside a plate using FEM. This problem uses bilinear bases functions.


The comfile run by this example is as follows:

#Example_314 Steady state temperature distribution in a plate (see example_111).
#This example is used to test the convergence of the Laplace or Poisson equation
#The basis function can be checked with either Bilinear or Bicubic Hermite    
# 
fem                               #Set the FEM environment
fem define parameters;r;laplin1;example  #Reads parameters needed to define the plate

fem define coordinates;r;;example #Reads in the coordinates used for the plate
fem define nodes;r;;example       #Reads in the 15 nodes to define the plate
fem define bases;r;;example       #Defines a bilinear basis function
fem define elements;r;;example    #Reads in eight elements from file
fem define wind;c                 #Calculates the graphic window size
fem draw lines                    #Makes the line segments visible
fem draw nodes                    #Makes the node numbers visible
fem define equation;r;;example    #Defines a STATIC analysis of LAPLACE's 
#                                 !equation (STANDARD) using GALERKIN 
#                                 !FINITE ELEMENTS with LINEAR ANALYSIS
fem define material;r;;example    #Define material properties (none)
fem define initial;r;;example     #Defines initial conditions
#                                 !Three sides are held at zero and the forth
#                                 !side has a heat function of Sin(Pi Y/L)
fem define solve;r;;example       #Defines solver options, default all settings
fem list node solution            #Lists the initial solutions
fem sol                           #Solves the boundary value problem
fem list node solution            #Lists the nodal solutions
fem list node flux                #Lists the nodal (heat) fluxes
fem draw field                    #This will draw colour temperature field
fem draw contours label off       #Draws 10 contours step size 0.1 deg

Additional testing commands:

fem list node solution

Files used by this example are:

Name              Modified     Size

example_3141.com 20-Nov-2001 1.8k laplin1.ipbase 10-Apr-2000 1.1k laplin1.ipcoor 10-Apr-2000 860 laplin1.ipelem 10-Apr-2000 2.2k laplin1.ipequa 26-May-2003 1.3k laplin1.ipinit 10-Apr-2000 720 laplin1.ipmate 10-Apr-2000 54 laplin1.ipnode 10-Apr-2000 2.7k laplin1.ippara 12-Nov-2002 5.9k laplin1.ipsolv 13-Apr-2007 1.1k laplin1.opvari 10-Apr-2000 2.2k test_output.com 10-Apr-2000 23

Download the entire example:

Name                           Modified     Size

examples_3_31_314_3141.tar.gz 14-Apr-2007 7.3k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:43 20160
cm-debugSuccessSat Mar 5 00:02:00 20160
mips-irix
cmSuccessSun Aug 19 01:24:31 20073
cm-debugSuccessWed Aug 15 01:20:43 20075
cm-debug-clear-mallocSuccessSat Aug 18 01:30:15 200712
cm-debug-clear-malloc7SuccessMon Aug 20 01:25:53 200712
cm64SuccessSun Aug 19 01:24:35 20073
cm64-debugSuccessTue Aug 21 01:19:48 20075
cm64-debug-clear-mallocSuccessThu Apr 1 10:12:31 20043
rs6000-aix
cmSuccessWed Mar 4 01:06:57 20091
cm-debugSuccessMon Mar 2 01:06:54 20091
cm64SuccessWed Mar 4 01:06:57 20091
cm64-debugSuccessTue Mar 3 01:11:38 20091
x86_64-linux
cmSuccessSun Mar 6 00:01:01 20160
cm-debugSuccessSat Mar 5 00:01:12 20160

Testing status by file:


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

Input last modified: Fri Apr 13 10:31:38 2007


CMISS Help / Examples / 3 / 31 / 314 / 3141