Example 724: Parameter estimation by minimising sums of squares of data errors

This example estimates a neo-hookian parameter by minimising the error between tracked data points and points predicted using finite elasticity.


The comfile run by this example is as follows:

#Example_724  

if ($OPTIMISE != 1) {
	fem define coordinate;r;uniaxial;example 
	fem define node;r;;example
	fem define base;r;;example
	fem define element;r;;example 
	fem define fibre;d
	fem define element;r;;example fibre
	fem define equation;r;;example
	fem define material;r;;example
	fem define initial;r;;example
	fem define solve;r;;example

	# this loads in the undeformed data,
	# and generate the measured data.
	fem define data;r;;example
	fem define xi;c contain
	fem solve step 1 error 1.0D-10
	fem define data;c from_xi def

	fem define optimise;r;;example

	$OPTIMISE = 1;
	optimise;
} else {
	fem solve incr 0.0 error 1.0D-10
	fem evaluate resid wrt data_fitting function deformed
}

Additional testing commands:

fem list stress;

Files used by this example are:

Name                 Modified     Size

example_724.com 17-Aug-2014 701 test_output.com 17-Aug-2014 17 uniaxial.ipbase 19-Apr-2004 3.4k uniaxial.ipcoor 19-Apr-2004 610 uniaxial.ipdata 19-Apr-2004 507 uniaxial.ipelem 19-Apr-2004 446 uniaxial.ipelfb 19-Apr-2004 275 uniaxial.ipequa 02-May-2004 2.1k uniaxial.ipinit 19-Apr-2004 1.5k uniaxial.ipmate 19-Apr-2004 2.3k uniaxial.ipnode 19-Apr-2004 1.9k uniaxial.ipopti 24-Jun-2005 1.3k uniaxial.ipsolv 16-Aug-2010 2.3k uniaxial.ipsolv.old 13-Apr-2007 2.2k

Download the entire example:

Name                      Modified     Size

examples_7_72_724.tar.gz 18-Aug-2014 7.3k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:02:05 20162
cm-debugSuccessSat Mar 5 00:05:07 20161
mips-irix
cmSuccessSun Aug 19 01:48:28 200715
cm-debugSuccessWed Aug 15 01:37:29 200731
cm-debug-clear-mallocSuccessSat Aug 18 01:50:18 200742
cm-debug-clear-malloc7SuccessMon Aug 20 01:46:23 200744
cm64SuccessSun Aug 19 01:48:16 200717
cm64-debugSuccessTue Aug 21 01:38:07 200733
rs6000-aix
cmSuccessWed Mar 4 01:08:57 20093
cm-debugSuccessMon Mar 2 01:08:54 20095
cm64SuccessWed Mar 4 01:08:57 20093
cm64-debugSuccessTue Mar 3 01:13:54 20094
x86_64-linux
cmSuccessSun Mar 6 00:01:04 20161
cm-debugSuccessSat Mar 5 00:01:15 20161

Testing status by file:


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

Input last modified: Sun Aug 17 20:33:32 2014


CMISS Help / Examples / 7 / 72 / 724