Example 48: 3D static linear elasticity with graphical o/p.

This example explains the use of CMISS in static 3D linear elasticity with graphical output of stresses.

The comfile run by this example is as follows:

fem de para;r;plate;example

#
#define geometry of plate
#
fem de coo 3,1                 #3D rect. cart. coords
fem de node;r;plate;example    #16 nodes
fem de base;r;;example         #1 trilinear basis & 1 bilinear basis
fem de elem;r;;example         #3 elements
fem de fibr;r;;example         #zero fibre angle
fem de elem;r;;example fibre   #fibre angle uses trilinear basis
fem group nodes 1,2,9,10 as bottom #these nodes get tied down
fem de wind
fem dr lines                   #draw mesh
#
#define equations, material props & b.c.s
#
fem de equa;r;;example         #3D linear elasticity
fem de mate;r;;example         #Young's mod si 10GPa. Poisson's ration is 0.3
fem de init;r;;example #tie down the bottom and load nodes 6 & 14 in y-dir.n
fem de solv;r;;example
fem solve
fem list node solution         #solution variables are displacements
fem dr lines deformed rgb=blue #draw deformed mesh
#
#fit principal stresses
#
fem de fiel;r;;example         #define field variable for fitting stresses
fem de elem;r;;example field
fem up gauss stress princ      #puts 3 principal stresses in YG array
fem li gauss yg elem 1         #lists princ stresses for element 1
fem de fit;r;;example gauss    #fits gauss pt data
fem fit
fem list node                  #fitted princ stresses are listed in field vars
fem dr fiel deformed           #draw principal stress as colour plot
fem dr reac scale 20           #draw reaction vectors
#
#export mesh and solution to cmgui
#
fem ex node geom as original_plate
fem ex elem geom
#

Files used by this example are:

Name            Modified     Size

example_48.com  20-Mar-2002  1.5k
plate.exelem    10-Apr-2000  7.0k
plate.exnode    10-Apr-2000  2.8k
plate.ipbase    10-Apr-2000  2.4k
plate.ipelem    10-Apr-2000  1.0k
plate.ipelfb    10-Apr-2000  526 
plate.ipelfd    10-Apr-2000  529 
plate.ipequa    26-May-2003  1.9k
plate.ipfibr    30-Jan-2001  1.6k
plate.ipfiel    10-Apr-2000  1.6k
plate.ipfit     13-Apr-2007  3.0k
plate.ipinit    10-Apr-2000  886 
plate.ipmate    10-Apr-2000  976 
plate.ipnode    10-Apr-2000  3.1k
plate.ippara    12-Nov-2002  5.9k
plate.ipsolv    13-Apr-2007  1.1k

Download the entire example:

Name                  Modified     Size

examples_4_48.tar.gz  14-Apr-2007  7.8k

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

Input last modified: Fri Apr 13 10:33:20 2007