Example 5i2: Mesh with 2 collased elements and an element adjacent to a collapsed one.

Mesh with 2 collased elements and an element adjacent to a collapsed one.

The comfile run by this example is as follows:

# Example 5i2 

fem def para;r;memory
#set diag on

fem de coord;r;allow_mapping

$name = "example_5i2";

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

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

   #fem update nodes deriv $j
}
#set diag off

# the scale factors are not correct after the above so fix them 
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

fem def init;r;displace_zx_3elements

# for testing
#fem def node;w;after
#fem def node;w;dep_after dep 1

#the above 2 files are the same

fem def solve;r;cubeherm

# for testing
#fem def node;w;after2
#fem def node;w;dep_after2 dep 1

# the above 2 files are different

# kludge to fix issue above
fem def mapping;r;$name
fem update solution geom

#fem eval reac
fem solv inc 0

fem solv inc 0.0001 iter 20 error 1e-4
fem solv inc 0.0009 iter 20 error 1e-4
fem solv inc 0.009 iter 20 error 1e-4
fem solv inc 0.09 iter 20 error 1e-6
fem solv inc 0.4 iter 20 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_5i2.com               09-Mar-2007  1.9k
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  658 
cubeherm.ipsolv               16-Aug-2010  3.1k
cubeherm.ipsolv.old           13-Apr-2007  3.0k
deformed.exelem               09-Mar-2007   24k
deformed.exnode               09-Mar-2007   14k
displace_zx_3elements.ipinit  09-Mar-2007  6.4k
example_5i2.ipelem            09-Mar-2007  4.1k
example_5i2.ipmap             09-Mar-2007   19k
example_5i2.ipnode            09-Mar-2007   32k
example_5i2_exp.ipelem        09-Mar-2007  4.3k
example_5i2_exp.ipnode        09-Mar-2007   30k
initial.exelem                09-Mar-2007   21k
initial.exnode                09-Mar-2007   13k
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
partial.ipinit                09-Mar-2007   94k
view.com                      09-Mar-2007  2.3k

Download the entire example:

Name                      Modified     Size

examples_5_5i_5i2.tar.gz  17-Aug-2010   24k

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

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