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:

StatusTestedReal time (s)
i686-linux
cmFailureSun Mar 6 00:13:50 20163
last breakTue Aug 12 00:37:00 20146
last successSun Nov 2 01:12:00 20037
cm-debugFailureSun Mar 6 00:20:46 20165
last breakTue Aug 12 00:43:00 201410
last successTue Nov 4 01:15:00 200311
mips-irix
cmFailureTue Aug 21 02:05:54 200724
last breakThu Apr 13 01:40:00 200617
last successSun Nov 2 01:37:00 200315
cm-debugFailureTue Aug 21 03:26:07 200753
last breakWed Apr 12 01:55:00 200646
last successTue Nov 4 01:47:00 200343
cm-debug-clear-mallocFailureTue Aug 21 04:04:32 200754
last breakSat Apr 15 01:51:00 200681
last successSun Nov 2 22:37:00 200349
cm-debug-clear-malloc7FailureTue Aug 21 04:15:25 200754
last breakMon Apr 17 01:51:00 200655
cm64FailureTue Aug 21 02:13:56 200719
last breakThu Apr 13 01:42:00 200616
cm64-debugFailureTue Aug 21 04:04:31 200749
last breakTue Apr 11 01:48:00 200649
cm64-debug-clear-mallocFailureTue Aug 21 05:36:02 200765
last breakThu Jun 15 01:37:00 200680
rs6000-aix
cmFailureWed Mar 4 01:12:48 20092
last breakFri Nov 3 02:12:00 20064
last successSun Nov 30 02:25:00 20033
cm-debugFailureWed Mar 4 01:23:47 200915
last breakThu Sep 13 02:09:00 200714
last successTue Nov 4 02:19:00 200314
cm64FailureWed Mar 4 01:12:34 20093
last breakFri Nov 3 02:13:00 20063
last successTue Feb 18 20:21:00 20034
cm64-debugFailureWed Mar 4 01:23:31 200914
last breakThu Sep 13 02:09:00 200713
x86_64-linux
cmFailureSun Mar 6 00:01:24 20161
last breakTue Aug 12 00:02:00 20142
cm-debugFailureSun Mar 6 00:01:53 20162
last breakTue Aug 12 00:40:00 20143

Testing status by file:


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

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


CMISS Help / Examples / 5 / 55 / 552 / 5521