Example 721: Material parameter estimation for an inflated 3D prolate spheroid using measured forces

A single prolate spheroidal element is defined with fibres oriented 65 degrees wrt the horizontal on the inside and -55 degrees wrt the horizontal on the outside. Geometry is interpolated isochorically (Focus^3.Cosh[lambda].Sinh^2[lambda]). Hydrostatic pressure is interpolated quartically in Xi3 (element based) and the pole-zero law defines the incompressible + fluid movement material behaviour. The prolate is inflated with a pressure of 2kPa on the inside face while the the pressure is held at zero on the outside face. The fibre and sheet axis coefficients of the material law are optimised wrt the difference between the calculated nodal forces and the experimentally determined forces (which in this case come from the solution of the boundary value problem in example_544).


The comfile run by this example is as follows:

#Example_721: Material parameter optimisation: deformed mesh of a single 
#inflated prolate element with rotating fibre and sheet fields; 
#estimating pole-zero coefficients

fem
fem define coordinates;r;param_opt;example           #3D prolate spheroidal coords
fem define nodes;r;;example                          #4 nodes of example_544
fem define bases;r;;example                          #trilinear + quartic pressure
fem define element;r;;example                        #single element
fem define fibre;r;;example                          #65deg inside, -55deg outside
fem define element;r;;example fibre                  #fibre elements
fem define equation;r;;example lock                  #Finite elasticity, FEM
fem define material;r;;example                       #pole-zero law
fem define initial;r;;example                        #inflation
fem define solve;r;;example                          #newton iterations
fem define nodes;r;param_opt_def;example dependent 1 #read deformed mesh
fem define optimise;r;param_opt;example              #opt for pole-zero coeffs
fem update pressure increment 1.0                #update pressure parameters
optimise                                         #optimised for material params

Additional testing commands:

#test update pressure calculation
fem list variable number 1
#test optimised parameters
fem list optim

Files used by this example are:

Name                         Modified     Size

example_721.com 10-Apr-2000 1.2k param_opt.ipbase 10-Apr-2000 4.3k param_opt.ipcoor 10-Apr-2000 707 param_opt.ipelem 10-Apr-2000 409 param_opt.ipelfb 10-Apr-2000 353 param_opt.ipequa 02-May-2004 2.1k param_opt.ipfibr 30-Jan-2001 1.4k param_opt.ipinit 12-Dec-2002 1.8k param_opt.ipmate 12-Dec-2002 5.9k param_opt.ipnode 10-Apr-2000 1.2k param_opt.ipopti 16-Sep-2001 1.3k param_opt.ippara 17-Dec-2002 5.0k param_opt.ipshee 10-Apr-2000 747 param_opt.ipsolv 16-Aug-2010 2.3k param_opt.ipsolv.old 13-Apr-2007 2.2k param_opt_def.ipinit 12-Dec-2002 6.1k param_opt_def.ipnode 10-Apr-2000 1.4k param_opt_def_old.ipinit 12-Dec-2002 6.3k param_opt_def_old.ipnode 10-Apr-2000 1.4k param_opt_def_recent.ipnode 10-Apr-2000 1.4k test_output.com 10-Apr-2000 103

Download the entire example:

Name                      Modified     Size

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

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:02:04 20161
cm-debugSuccessSat Mar 5 00:02:24 20160
mips-irix
cmSuccessSun Aug 19 01:26:14 20073
cm-debugSuccessWed Aug 15 01:24:38 20077
cm-debug-clear-mallocSuccessSat Aug 18 01:26:09 200710
cm-debug-clear-malloc7SuccessMon Aug 20 01:24:34 20079
cm64SuccessSun Aug 19 01:27:52 20074
cm64-debugSuccessTue Aug 21 01:22:18 20077
cm64-debug-clear-mallocSuccessThu Apr 1 10:40:57 20043
rs6000-aix
cmSuccessWed Mar 4 01:07:29 20091
cm-debugSuccessMon Mar 2 01:06:25 20092
cm64SuccessWed Mar 4 01:07:29 20091
cm64-debugSuccessTue Mar 3 01:11:56 20092
x86_64-linux
cmSuccessSun Mar 6 00:01:03 20160
cm-debugSuccessSat Mar 5 00:01:14 20160

Testing status by file:


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

Input last modified: Mon Aug 16 11:21:50 2010


CMISS Help / Examples / 7 / 72 / 721