Deformation of a tricubic hermite wedge

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:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:07:34 2016	7
cm-debug	Success	Sat Mar 5 01:00:57 2016	47
mips-irix
cm	Success	Sun Aug 19 02:58:37 2007	83
cm-debug	Success	Wed Aug 15 04:48:40 2007	424
cm-debug-clear-malloc	Success	Sat Aug 18 05:57:43 2007	466
cm-debug-clear-malloc7	Success	Mon Aug 20 05:19:57 2007	479
cm64	Success	Sun Aug 19 03:05:15 2007	95
cm64-debug	Success	Tue Aug 21 03:43:09 2007	471
cm64-debug-clear-malloc	Success	Tue Feb 1 09:55:00 2005	189
rs6000-aix
cm	Success	Wed Mar 4 01:14:03 2009	11
cm-debug	Success	Mon Mar 2 01:41:50 2009	153
cm64	Success	Wed Mar 4 01:14:04 2009	11
cm64-debug	Success	Tue Mar 3 01:40:20 2009	139
x86_64-linux
cm	Success	Sun Mar 6 00:01:16 2016	4
cm-debug	Success	Sat Mar 5 00:04:26 2016	21

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

CMISS Help / Examples / 5 / 54 / 542