Example 7183: Estimating material parameters of in-vivo human skin

This example is a testing version of example 7182, using a reduced numbers of data sets and bi-linear basis functions. This example executes in a few minutes, while the full examples 7181 and 7182 requires hours to complete

The comfile run by this example is as follows:

#This example is a scaled down testing vesion of example 7182. Two material parameters are estimated using two deformed states and #bi-linear basis functions
# If the example path is not set, default to current directory
if (!defined $example) {
    $example = "./";
}
# Drop off the trailing / in the example path
$chopped = chop $example;
if ($chopped ne "/") {
    $example .= $chopped;
}
#**************************************************************
@DATASETS = (12,14)								#Deformed states

fem define parameters;r;multiaxes;example
	
fem define bases;r;bases_line;example						#Bi-linear
fem define nodes;r;coordinates_line_undeformed;example
fem define elements;r;elements_line;example					#Circular mesh with 192 elements

fem define fibre;r;fibre;example
fem define elements;r;fibreelem_line;example fibre

fem define equation;r;equation;example						#Incompressible membrane (uses plain-stress)
fem define material;r;tongfung;example						#Tong and Fungs skin function

fem group nodes 1..32 as BOUNDARY

fem define data;r;datapoints_undeformed;example					#Reads the undeformed coordinates of the datapoints
foreach $index (@DATASETS[1..$#DATASETS]){
	fem define data;r;datapoints_undeformed;example append
}

fem define xi;c closest								#Defines the xi positions of the datapoints

fem define data;r;datapoints_${DATASETS[0]};example as IMAGE_${DATASETS[0]}	#Reads the deformed coordinates of the datapoints
foreach $index (@DATASETS[1..$#DATASETS]){
	fem define data;r;datapoints_${index};example append as IMAGE_${index}
}
symlink "$example/evalresid_line.com","evalresid_line.com"
fem define initial;r;initskinsolve_line;example					#Defines a force boundary problem
fem define solve;r;solve;example
fem define optimise;r;tongfung_line;example					#Estimate twe exponential material parameters

fem define nodes;r;displacements_line_$DATASETS[-1];example dep 2 nc 1
fem define nodes;r;force_line_$DATASETS[-1];example dep 2 nc 2			#Reads the experimentaly applied forces
fem define nodes;r;force_line_residual;example dep 1 nc 2			#Reads the estimated residual forces

fem solve increment 0.1 iterate 100
for $M (1..3)
{
	fem solve increment 0.3 iterate 100					#Solves the force boundary
}

optimise									#The optimiser calles the eval_resid_line.com file

fem list data errors deformed
fem list material constitutive

Additional testing commands:

fem list material constitutive

Files used by this example are:

Name                                Modified     Size

example_7183.com                    24-May-2004  2.3k
bases_line.ipbase                   03-May-2004  1.1k
coordinates_line_undeformed.ipnode  03-May-2004   31k
datapoints_12.ipdata                03-May-2004   37k
datapoints_14.ipdata                03-May-2004   37k
datapoints_undeformed.ipdata        03-May-2004   37k
displacements_line_12.ipnode        03-May-2004   31k
displacements_line_14.ipnode        03-May-2004   31k
elements_line.ipelem                03-May-2004   51k
equation.ipequa                     03-May-2004  1.8k
evalresid_line.com                  03-May-2004  744 
fibre.ipfibr                        03-May-2004   53k
fibreelem_line.ipelfb               03-May-2004   79k
force_line_12.ipnode                03-May-2004   31k
force_line_14.ipnode                03-May-2004   31k
force_line_residual.ipnode          03-May-2004   31k
initskin_line.ipinit                03-May-2004  1.3k
initskinsolve_line.ipinit           03-May-2004  1.2k
multiaxes.ippara                    14-Apr-2008  5.8k
skinforce_line.com                  03-May-2004  1.1k
solve.ipsolv                        16-Aug-2010  2.3k
solve.ipsolv.old                    13-Apr-2007  2.2k
test_output.com                     03-May-2004   31 
tongfung.ipmate                     24-May-2004  5.2k
tongfung_line.ipopti                24-May-2004  1.2k

Download the entire example:

Name                           Modified     Size

examples_7_71_718_7183.tar.gz  17-Aug-2010   91k

Testing status by version:

Testing status by file:

• cmiss_test.out
  

  i686-linux
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  mips-irix
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
  Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
  Success	cm64:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug:	ndiff test: no significant differences with generic answer.
  rs6000-aix
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.
  Success	cm64:	ndiff test: no significant differences with generic answer.
  Success	cm64-debug:	ndiff test: no significant differences with generic answer.
  x86_64-linux
  Success	cm:	ndiff test: no significant differences with generic answer.
  Success	cm-debug:	ndiff test: no significant differences with generic answer.

  
  

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

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