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


CMISS Help / Examples / 4 / 48