Example 542: Deformation of a tricubic hermite wedge

Example of explicit pressure boundary conditions applied to a Tricubic Hermite 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
#
#
# This example has been set up for various combinations of basis functions.
# To run problem for a particular combination uncomment the appropriate section
# in the basis and boundary conditions definitions
#
#
#
# Set up local paths
#

# Command updated by fixcom.sh on Wed Aug 23 17:48:57 NZT 2000
# Old command: assign variable OUTPUT_DIR "output" # local directory for problem output
$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
#
# Basis definitions
#
# CubicHermite-2Linear
#fem define bases;r;CubicHermite_2Linear_3x3x3Gauss_ArithMeanSE;example reg 1    ! Geometry
#fem define ;add bases;r;Linear_CubicHermite_3x3Gauss_ArithMeanSE;example reg 1  ! Faces
#fem define ;add bases;r;5Auxiliary_3x3x3Gauss;example reg 1                     ! Pressure
#fem define ;add bases;r;CubicHermite_Linear_3x3Gauss_ArithMeanSE;example reg 1  ! Faces
#fem define ;add bases;r;2Linear_3x3Gauss;example reg 1                          ! Faces
# OR
# 2CubicHermite-Linear constant pressure
#fem define bases;r;2CubicHermite-Linear_3x3x3Gauss_ArithMeanSE;example reg 1    ! Geometry
#fem define ;add bases;r;2CubicHermite_3x3Gauss_ArithMeanSE;example reg 1        ! Faces
#fem define ;add bases;r;5Auxiliary_3x3x3Gauss;example reg 1                     ! Pressure
#fem define ;add bases;r;CubicHermite_Linear_3x3Gauss_ArithMeanSE;example reg 1  ! Faces
#fem define ;add bases;r;Linear_CubicHermite_3x3Gauss_ArithMeanSE ;examplereg 1  ! Faces
# OR
# 2CubicHermite-Linear 3linear pressure
#fem define bases;r;2CubicHermite-Linear_3x3x3Gauss_ArithMeanSE;example reg 1    ! Geometry
#fem define ;add bases;r;2CubicHermite_3x3Gauss_ArithMeanSE;example reg 1        ! Faces
#fem define ;add bases;r;3Linear-2PressAuxXi3_3x3x3Gauss;example reg 1           ! Pressure
#fem define ;add bases;r;2Linear_3x3Gauss;example reg 1                          ! Faces
#fem define ;add bases;r;CubicHermite_Linear_3x3Gauss_ArithMeanSE;example reg 1  ! Faces
#fem define ;add bases;r;Linear_CubicHermite_3x3Gauss_ArithMeanSE;example reg 1  ! Faces
# OR 
# 3CubicHermite constant pressure
#fem define bases;r;3CubicHermite_3x3x3Gauss_ArithMeanSE;example reg 1           ! Geometry
#fem define ;add bases;r;2CubicHermite_3x3Gauss_ArithMeanSE;example reg 1        ! Faces
#fem define ;add bases;r;5Auxiliary_3x3x3Gauss;example reg 1                     ! Pressure
# OR
# 3CubicHermite 3linear pressure
fem define bases;r;3CubicHermite_3x3x3Gauss_ArithMeanSE;example reg 1            # Geometry
fem define ;add bases;r;2CubicHermite_3x3Gauss_ArithMeanSE;example reg 1         # Faces
fem define ;add bases;r;3Linear-2PressAuxXi3_3x3x3Gauss;example reg 1            # Pressure
fem define ;add bases;r;2Linear_3x3Gauss;example reg 1                           # Faces
#
# Geometry
#
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
#
# Method
#
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
#
# Boundary Conditions
#
# CubicHermite-2Linear
#fem define initial;r;pressure_CubicHermite-2Linear;example reg 1            -
# OR                                                                          | BC's on derivs 
# 2CubicHermite-Linear constant pressure                                      | in z dirn
#fem define initial;r;pressure_2CubicHermite-Linear_constpres;example reg 1   | not quite
# OR                                                                          | correct
# 2CubicHermite-Linear 3linear pressure                                       |
#fem define initial;r;pressure_2CubicHermite-Linear_3linpres;example reg 1   -
# OR
# 3CubicHermite constant pressure
#fem define initial;r;pressure_3CubicHermite_constpres;example reg 1
# OR
# 3CubicHermite 3linear pressure
fem define initial;r;pressure_3CubicHermite_3linpres;example reg 1
#
#
#
#Solve 
#
fem define solve;r;newton;example reg 1
set out;steps_inflate on
$NAME1 = "wedge_def_00"
$NAME2 = "_active_0000"
$ITERS = 40
$TOL = 0.001

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

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

  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_542.com 20-Jul-2001 5.2k 2CubicHermite-Linear_3x3x3Gauss_ArithMeanSE.ipbase 10-Apr-2000 1.8k 2CubicHermite_3x3Gauss_ArithMeanSE.ipbase 10-Apr-2000 1.5k 2Linear_3x3Gauss.ipbase 10-Apr-2000 1.1k 3CubicHermite_3x3x3Gauss_ArithMeanSE.ipbase 10-Apr-2000 1.9k 3Linear-2PressAuxXi3_3x3x3Gauss.ipbase 10-Apr-2000 1.9k 3Linear_3x3x3Gauss.ipbase 10-Apr-2000 1.4k 3d_Standard_RC.ipcoor 10-Apr-2000 606 5Auxiliary_3x3x3Gauss.ipbase 10-Apr-2000 1.9k CubicHermite_2Linear_3x3x3Gauss_ArithMeanSE.ipbase 10-Apr-2000 1.7k CubicHermite_Linear_3x3Gauss_ArithMeanSE.ipbase 10-Apr-2000 1.5k Incompressible_Finite_Elasticity.ipequa 02-May-2004 2.1k Linear_CubicHermite_3x3Gauss_ArithMeanSE.ipbase 10-Apr-2000 1.5k Orthotropic_Incompressible.ipmate 05-Dec-2002 6.2k active0_00.ipacti 03-Mar-2004 777 newton.ipsolv 16-Aug-2010 2.7k newton.ipsolv.old 13-Apr-2007 2.5k pressure_2CubicHermite-Linear_3linpres.ipinit 05-Dec-2002 4.0k pressure_2CubicHermite-Linear_constpres.ipinit 05-Dec-2002 4.4k pressure_3CubicHermite_3linpres.ipinit 17-Dec-2002 6.8k pressure_3CubicHermite_constpres.ipinit 05-Dec-2002 6.7k pressure_3Linear.ipinit 05-Dec-2002 2.3k pressure_CubicHermite-2Linear.ipinit 05-Dec-2002 2.5k pressure_inflation.ipinit 05-Dec-2002 1.9k pressure_inflation.ippara 12-Nov-2002 5.9k steps_inflate.out 17-Dec-2002 93k test_output.com 10-Apr-2000 26 view.com 20-Nov-2001 1.9k wedge.ipelem 10-Apr-2000 1.4k wedge.ipelfb 10-Apr-2000 630 wedge.ipfibr 30-Jan-2001 9.0k wedge.ipnode 10-Apr-2000 30k wedge.jpg 10-Apr-2000 20k

Download the entire example:

Name                      Modified     Size

examples_5_54_542.tar.gz 17-Aug-2010 42k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:07:34 20167
cm-debugSuccessSat Mar 5 01:00:57 201647
mips-irix
cmSuccessSun Aug 19 02:58:37 200783
cm-debugSuccessWed Aug 15 04:48:40 2007424
cm-debug-clear-mallocSuccessSat Aug 18 05:57:43 2007466
cm-debug-clear-malloc7SuccessMon Aug 20 05:19:57 2007479
cm64SuccessSun Aug 19 03:05:15 200795
cm64-debugSuccessTue Aug 21 03:43:09 2007471
cm64-debug-clear-mallocSuccessTue Feb 1 09:55:00 2005189
rs6000-aix
cmSuccessWed Mar 4 01:14:03 200911
cm-debugSuccessMon Mar 2 01:41:50 2009153
cm64SuccessWed Mar 4 01:14:04 200911
cm64-debugSuccessTue Mar 3 01:40:20 2009139
x86_64-linux
cmSuccessSun Mar 6 00:01:16 20164
cm-debugSuccessSat Mar 5 00:04:26 201621

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 / 542