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:

Testing status by file:

• cmiss_test.log
  

  i686-linux
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  mips-irix
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  Success	cm-debug-clear-malloc:	cmiss_test.log.retain.
  Success	cm-debug-clear-malloc7:	cmiss_test.log.retain.
  Success	cm64:	cmiss_test.log.retain.
  Success	cm64-debug:	cmiss_test.log.retain.
  Success	cm64-debug-clear-malloc:	cmiss_test.log.retain.
  rs6000-aix
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.
  Success	cm64:	cmiss_test.log.retain.
  Success	cm64-debug:	cmiss_test.log.retain.
  x86_64-linux
  Success	cm:	cmiss_test.log.retain.
  Success	cm-debug:	cmiss_test.log.retain.

  
  
• cmiss_test.out
  

  i686-linux
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  mips-irix
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
  Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
  Success	cm64:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
  rs6000-aix
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  Success	cm64:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug:	ndiff test: no significant differences with generic answer.
  x86_64-linux
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.

  
  

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

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