Example 5521: Passive inflation and active contraction of a trilinear prolate spheroid

A single prolate is passively inflated to 4 kPa then the made to actively contract.

The comfile run by this example is as follows:

#Example_5521  Inflation and active contraction of an orthotropic prolate element

fem
fem define param;r;prol;example                 # 3D prolate spheroidal
fem define coordinate;r;profib;example          # 4 nodes of previous example_5512
fem define node;r;profib;example                # trilinear + quadratic hydr. pressure
fem define base;r;profib;example                # single axisymmetric element
fem define element;r;profib;example             # fibre angles of previous example_5512
fem define fibre;r;profib;example               # fibre elements
fem define element;r;profib;example fibre
fem define equation;r;profib;example            # 3D, finite elasticity
# 
# INFLATION
# 
fem define material;r;profib;example            # incomp; pole-zero law
fem define initial;r;profib;example             # Internal pressure of 4kPa        
fem define solve;r;profib;example               # Newton iterations
fem solve increment 0.2 iter 8                  # Gradually increment Pressure
fem solve increment 0.3 iter 11 
fem solve increment 0.3 iter 8
fem solve increment 0.1 iter 6
fem eval reac                                   # Inflated solution is converged OK! 
#
fem define window 
fem draw lines 
fem draw lines dot def 
refresh graphics
# 
# EJECTION 
# 
fem define material;r;profib_active;example     # incomp; pole-zero law + active fibre stress
fem define active;r;profib;example              # steady state calcium/length/tension
fem define solve;r;profib;example               # Newton iterations
fem solve increment 0                           # check current solution is converged
fem draw;add lines deformed dotted rgb=blue
refresh graphics


fem define active;r;profib_active;example       # Increase activation level  
fem solve increment 0 iterate 10                # Solve with new activation level but keep press const.
fem define initial;w;profib_endeject            # Save deformed configurations for end ejection 
fem draw;add lines deformed dotted rgb=blue
refresh graphics

Additional testing commands:

fem list initial 

Files used by this example are:

Name                    Modified     Size

example_5521.com        10-Apr-2000  2.0k
profib.ipacti           03-Mar-2004  715 
profib.ipbase           10-Apr-2000  4.2k
profib.ipcoor           10-Apr-2000  707 
profib.ipelem           10-Apr-2000  409 
profib.ipelfb           10-Apr-2000  349 
profib.ipelfd           10-Apr-2000  276 
profib.ipequa           02-May-2004  2.0k
profib.ipfibr           30-Jan-2001  1.3k
profib.ipfiel           10-Apr-2000  601 
profib.ipfit            13-Apr-2007  1.5k
profib.ipinit           12-Dec-2002  1.7k
profib.ipmate           12-Dec-2002  5.7k
profib.ipnode           10-Apr-2000  1.2k
profib.ipsolv           16-Aug-2010  2.3k
profib.ipsolv.old       13-Apr-2007  2.1k
profib_active.ipacti    03-Mar-2004  717 
profib_active.ipmate    12-Dec-2002  6.2k
profib_endeject.ipinit  12-Dec-2002  6.1k
prol.ippara             12-Nov-2002  5.9k
test_output.com         10-Apr-2000   18 

Download the entire example:

Name                           Modified     Size

examples_5_55_552_5521.tar.gz  17-Aug-2010   11k

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
  Failure	cm:	ndiff test: significant differences with generic answer.
  Failure	cm-debug:	ndiff test: significant differences with generic answer.
  mips-irix
  Failure	cm:	ndiff test: significant differences with generic answer.
  Failure	cm-debug:	ndiff test: significant differences with generic answer.
  Failure	cm-debug-clear-malloc:	ndiff test: significant differences with generic answer.
  Failure	cm-debug-clear-malloc7:	ndiff test: significant differences with generic answer.
  Failure	cm64:	ndiff test: significant differences with generic answer.
  Failure	cm64-debug:	ndiff test: significant differences with generic answer.
  Failure	cm64-debug-clear-malloc:	ndiff test: significant differences with generic answer.
  rs6000-aix
  Failure	cm:	ndiff test: significant differences with generic answer.
  Failure	cm-debug:	ndiff test: significant differences with generic answer.
  Failure	cm64:	ndiff test: significant differences with generic answer.
  Failure	cm64-debug:	ndiff test: significant differences with generic answer.
  x86_64-linux
  Failure	cm:	ndiff test: significant differences with generic answer.
  Failure	cm-debug:	ndiff test: significant differences with generic answer.

  
  

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

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