Example 5i1: Mesh with 2 collapsed elements

Mesh with 2 collapsed elements and no elements adjacent to the collapsed elements.

The comfile run by this example is as follows:

#Example 5i1

# Finite elasticity using a mesh with elements having collapsed nodes

fem def para;r;memory        

fem de coord;r;allow_mapping

$name = "example_5i1";

fem define bas;r;TriCubicHerm
fem def;add base;r;BiCubicHerm
fem def;add base;r;TriLinearPress
fem def;add base;r;BiLinear

fem define node;r;$name
fem define elem;r;$name

fem def fibre;d
fem def elem;d          fibre

for $j (1..3){
   #This option will make all the lines straight
   fem update nodes deriv $j versions individual

   # this option will make the lines curved
   #fem update nodes deriv $j

# not sure why this is needed, but the scale factors might
# be wrong if it isn't included
fem update node scale_factor

# Export initial mesh
$ref = "initial";

fem export nodes;$ref as $ref;
fem export elements;$ref as $ref;

fem define node;w;$name."_exp"
fem define elem;w;$name."_exp"

# Refine doesn't work because the mappings are not updated.
#fem refine xi 1
#fem refine xi 2
#fem refine xi 3

fem def equa;r;mechanics          lock
fem def mat;r;mooney_rivlin

fem group element all external s3=0 as BOTTOM;
fem group nodes xi3=0 element BOTTOM as BOTTOMNODES;

fem group element all external s3=1 as TOP;
fem group nodes xi3=1 element TOP as TOPNODES;

fem list nodes groups

# define boundary conditions, in this example we displacements
# defined

fem def init;r;displace_zx

fem def solve;r;cubeherm

# After loading the ipsolv file ZP may not have the correct
# values. This doesn't matter here because the affected DOFs
# are fixed by the boundary conditions.

#explicitly setting the mapping fixes the above

fem def mapping;r;$name
fem update solution geometry

fem solv inc 0 error 1e-4

fem solv inc 0.0001 iter 5 error 1e-4
fem solv inc 0.0009 iter 5 error 1e-4
fem solv inc 0.009 iter 5 error 1e-4
fem solv inc 0.09 iter 5 error 1e-6
fem solv inc 0.4 iter 5 error 1e-8
fem solv inc 0.5 iter 30 error 1e-10

$name = "deformed";

fem export nodes;$name as $name field offset 0
fem export elements;$name as $name field offset_elem 0

#view in cmgui if desired
system "cmgui view.com"

Files used by this example are:

Name                    Modified     Size

example_5i1.com 09-Mar-2007 2.2k BiCubicHerm.ipbase 09-Mar-2007 1.6k BiLinear.ipbase 09-Mar-2007 1.2k TriCubicHerm.ipbase 09-Mar-2007 1.9k TriLinear.ipbase 09-Mar-2007 1.4k TriLinearPress.ipbase 09-Mar-2007 1.9k allow_mapping.ipcoor 09-Mar-2007 688 cubeherm.ipsolv 16-Aug-2010 3.1k cubeherm.ipsolv.old 13-Apr-2007 3.0k deformed.exelem 09-Mar-2007 16k deformed.exnode 09-Mar-2007 11k demo1_exp.ipelem 09-Mar-2007 3.5k demo1_exp.ipnode 09-Mar-2007 23k displace_zx.ipinit 09-Mar-2007 7.0k example_5i1.ipelem 09-Mar-2007 3.3k example_5i1.ipmap 09-Mar-2007 19k example_5i1.ipnode 09-Mar-2007 24k example_5i1_exp.ipelem 09-Mar-2007 3.5k example_5i1_exp.ipnode 09-Mar-2007 23k initial.exelem 09-Mar-2007 14k initial.exnode 09-Mar-2007 10k mechanics.ipequa 09-Mar-2007 2.1k memory.ippara 09-Mar-2007 5.9k mooney_rivlin.ipmate 09-Mar-2007 2.3k newton.ipsolv 16-Aug-2010 3.1k newton.ipsolv.old 13-Apr-2007 3.0k view.com 09-Mar-2007 2.2k

Download the entire example:

Name                      Modified     Size

examples_5_5i_5i1.tar.gz 17-Aug-2010 17k

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

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

CMISS Help / Examples / 5 / 5i / 5i1