Example 541: Deformation of a trilinear wedge

Example of explicit pressure boundary conditions applied to a trilinear wedge in 3 dimensional rectangular cartesian coordinates.


The comfile run by this example is as follows:

# This example solves a finite elasticity deformation of a wedge 
# loaded by pressure on one face. The symmetry of the problem 
# likens itself to a square-cornered cylinder with a pressure bc 
# on the 'inside' face
#    
#      y
#      ^
#      |
#    2 |\
#      | \
#      |  \     with unit depth
#    1 |\  \
#      | \  \
#      |  \  \
#      --------------> x  
#         1  2
#
# where the x-face (1,0) -> (2,0) was constrained to slide along the
# x-axis, the y-face (0,1) -> (0,2) was constrained to slide 
# along y and one z-face is completely fixed


#
# Set up local paths
#

$OUTPUT_DIR = "output" # local directory for problem output


#
# Set up problem
#
fem
fem define parameters;r;pressure_inflation;example
fem define coordinates;r;3d_Standard_RC;example reg 1
fem define bases;r;3Linear_3x3x3Gauss;example reg 1          # Geometry
fem define ;add bases;r;2Linear_3x3Gauss;example reg 1       # Faces
fem define ;add bases;r;5Auxiliary_3x3x3Gauss;example reg 1  # Pressure
fem define nodes;r;wedge;example reg 1 
fem define elements;r;wedge;example reg 1
fem define fibres;r;wedge;example reg 1
fem define elements;r;wedge;example fibre reg 1
fem define equation;r;Incompressible_Finite_Elasticity;example reg 1
fem define material;r;Orthotropic_Incompressible;example reg 1
fem define active;r;active0_00;example reg 1
fem define initial;r;pressure_inflation;example reg 1
fem define solve;r;newton;example reg 1

#
# Inital geometry
#
fem export nodes;wedge;$OUTPUT_DIR as wedge
fem export elem;wedge;$OUTPUT_DIR as wedge

#
#Solve 
#
set out;steps_inflate on

$NAME1 = "wedge_def_00"
$NAME2 = "_active_0000"

$ITERS = 5
$TOL = 0.0001

#
#increase pressure incrementally
#
for $PRESS ( 0..5 ) 
{

  if( $PRESS < 1 )
  {
    $FILENAME = "${NAME1}${PRESS}$NAME2"
    $INCREM = 0.0
  } else {
    $FILENAME = "${NAME1}${PRESS}$NAME2"
    $INCREM = 0.2
  }

  list assign 
  fem solve increment $INCREM iter $ITERS error $TOL
  fem export nodes;$FILENAME;$OUTPUT_DIR field as wedge
  fem export elements;$FILENAME;$OUTPUT_DIR field as wedge
  fem define initial;w;$FILENAME;$OUTPUT_DIR

}


set out off

fem define window
fem draw lines
fem draw lines deformed dotted 

Additional testing commands:

fem list node solution




Files used by this example are:

Name                                     Modified     Size

example_541.com 09-Jul-2001 2.1k 2Linear_3x3Gauss.ipbase 09-Jul-2001 1.1k 3Linear_3x3x3Gauss.ipbase 09-Jul-2001 1.4k 3d_Standard_RC.ipcoor 10-Apr-2000 606 5Auxiliary_3x3x3Gauss.ipbase 10-Apr-2000 1.9k Incompressible_Finite_Elasticity.ipequa 02-May-2004 2.1k Orthotropic_Incompressible.ipmate 05-Dec-2002 6.2k active0_00.ipacti 03-Mar-2004 777 newton.ipsolv 16-Aug-2010 2.8k newton.ipsolv.old 13-Apr-2007 2.7k pressure_inflation.ipinit 05-Dec-2002 1.9k pressure_inflation.ippara 12-Nov-2002 5.9k steps_inflate.out 09-Jul-2001 22k test_example.out 10-Apr-2000 26k test_output.com 10-Apr-2000 26 view.com 20-Nov-2001 1.9k wedge.ipelem 10-Apr-2000 460 wedge.ipelfb 10-Apr-2000 232 wedge.ipfibr 30-Jan-2001 1.0k wedge.ipnode 10-Apr-2000 2.2k

Download the entire example:

Name                      Modified     Size

examples_5_54_541.tar.gz 17-Aug-2010 14k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:02:01 20160
cm-debugSuccessSat Mar 5 00:05:04 20162
mips-irix
cmSuccessSun Aug 19 01:34:30 20078
cm-debugSuccessWed Aug 15 01:30:47 200718
cm-debug-clear-mallocSuccessSat Aug 18 01:38:17 200722
cm-debug-clear-malloc7SuccessMon Aug 20 01:34:46 200722
cm64SuccessSun Aug 19 01:36:05 20078
cm64-debugSuccessTue Aug 21 01:30:37 200719
cm64-debug-clear-mallocSuccessTue Feb 1 09:45:41 20058
rs6000-aix
cmSuccessWed Mar 4 01:08:05 20091
cm-debugSuccessMon Mar 2 01:08:24 20095
cm64SuccessWed Mar 4 01:08:05 20092
cm64-debugSuccessTue Mar 3 01:13:53 20095
x86_64-linux
cmSuccessSun Mar 6 00:01:03 20160
cm-debugSuccessSat Mar 5 00:01:14 20161

Testing status by file:


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

Input last modified: Mon Aug 16 11:19:15 2010


CMISS Help / Examples / 5 / 54 / 541