Example 723: Nonlinear material parameter optimisation with foward solving

This example demonstrates nonlinear optimisation of material parameters by solving the forward problem

The comfile run by this example is as follows:

#Example 723
# This example demonstrates nonlinear optimisation 
# of material parameters by solving the forward problem
#
# This command file is used by the optimiser to solve the 
# forward problem for each objective evaluation


$TRUE = 1;  # Initialise logicals
$FALSE = 0;

if($RESOLVE ne "TRUE") 
{ 
  # Don't need to redo these commands when 
  # resolving during optimisation
  fem define parameters;r;shear;example;
  fem define coordinates;r;shear;example;

  fem def base;r;linear;example;
  fem def ;add base;r;lagr;example;
}

fem def node;r;shear_specimen_ref;example;
fem def elem;r;shear_ref;example;
fem def fibre;r;shear_ref;example;
fem def elem;r;shear_ref;example fibre;

fem group nodes 1,3,5,7 as LEFT;
fem group nodes 2,4,6,8 as RIGHT;
fem group nodes 1,2,5,6 as BOTTOM;
fem group nodes 3,4,7,8 as TOP;
fem group nodes 1,2,3,4 as BACK;
fem group nodes 5,6,7,8 as FRONT;

# Refinement cannot be done during in 
# the optimisation run. Refine mesh
# and write it out before the optimisation
# run e.g.
#
# (new nodes and elements will automatically
#  be added to the correct groups)
#
#fem refine xi 3;
#fem def node;w;shear_specimen_ref;
#fem def elem;w;shear_ref;
#fem def fibre;w;shear_ref; 
#fem def elem;w;shear_ref fibre;

#fem export nodes;shear;example as cube;
#fem export element;shear;example as cube;

fem def equa;r;shear;example;

fem def mate;r;guiccione_can_mid;example;

if($RESOLVE eq "TRUE") 
{ 
  # Need to update the material parameters 
  # to those being evaluated by the optimiser
  fem update optimisation
}

fem def init;r;shear_NF;example;

fem def solve;r;shear;example;

fem solve increment 0 error 0.01;

fem solve increment 0.2 iter 30 error 0.01;
if($RESOLVE eq "TRUE") 
{
  # Evaluate residuals & load residual vector - 10%
  fem evaluate reaction nodes FRONT direction 1 residual 1
  fem evaluate reaction nodes FRONT direction 2 residual 2
  fem evaluate reaction nodes FRONT direction 3 residual 3
}


fem solve increment 0.2 iter 30 error 0.01;
if($RESOLVE eq "TRUE") 
{
  # Evaluate residuals & load residual vector - 20% shear
  fem evaluate reaction nodes FRONT direction 1 residual 4
  fem evaluate reaction nodes FRONT direction 2 residual 5
  fem evaluate reaction nodes FRONT direction 3 residual 6
}

fem solve increment 0.2 iter 30 error 0.01;
if($RESOLVE eq "TRUE") 
{
  # Evaluate residuals & load residual vector - 30% shear
  fem evaluate reaction nodes FRONT direction 1 residual 7
  fem evaluate reaction nodes FRONT direction 2 residual 8
  fem evaluate reaction nodes FRONT direction 3 residual 9
}

fem solve increment 0.2 iter 30 error 0.01;
if($RESOLVE eq "TRUE") 
{
  # Evaluate residuals & load residual vector - 40% shear
  fem evaluate reaction nodes FRONT direction 1 residual 10
  fem evaluate reaction nodes FRONT direction 2 residual 11
  fem evaluate reaction nodes FRONT direction 3 residual 12
}

fem solve increment 0.2 iter 30 error 0.01;
if($RESOLVE eq "TRUE") 
{
  # Evaluate residuals & load residual vector - 50% shear
  fem evaluate reaction nodes FRONT direction 1 residual 13
  fem evaluate reaction nodes FRONT direction 2 residual 14
  fem evaluate reaction nodes FRONT direction 3 residual 15
}


if($RESOLVE eq "TRUE") 
{
  $ITERATION=$ITERATION + 1;
  print "\033[0;30;43m                    Iteration = ${ITERATION}                  \033[0m\n"; 
} 

if($RESOLVE ne "TRUE") 
{
  # Load the experimental data, force components in the 
  # x, y and z directions are in indecies 1, 2 and 3 
  # respectively
  fem define data;r;shear_NF;example;
  #
  # Create a vector of the experimental reactions
  # for use by the optimiser 
  fem update residual data point 1 index 1 residual 1
  fem update residual data point 1 index 2 residual 2
  fem update residual data point 1 index 3 residual 3
  #
  fem update residual data point 2 index 1 residual 4
  fem update residual data point 2 index 2 residual 5
  fem update residual data point 2 index 3 residual 6
  #
  fem update residual data point 3 index 1 residual 7
  fem update residual data point 3 index 2 residual 8
  fem update residual data point 3 index 3 residual 9
  #
  fem update residual data point 4 index 1 residual 10
  fem update residual data point 4 index 2 residual 11
  fem update residual data point 4 index 3 residual 12
  #
  fem update residual data point 5 index 1 residual 13
  fem update residual data point 5 index 2 residual 14
  fem update residual data point 5 index 3 residual 15

  $ITERATION = 1;
  $RESOLVE = "TRUE";
  fem define optimise;r;shear;example;           
  optimise; 
}


#fem export nodes;shear_def;example field as cube;
#fem export element;shear_def;example field as cube;

Additional testing commands:

#Testing example 723
FEM update material optimise;
FEM li material constitutive;

Files used by this example are:

Name                       Modified     Size

example_723.com            24-Oct-2003  4.6k
guiccione_can_mid.ipmate   13-Oct-2003  2.8k
junk                       05-Jun-2003   63k
lagr.ipbase                05-Jun-2003  1.1k
linear.ipbase              05-Jun-2003  2.4k
shear.exelem               05-Jun-2003  6.8k
shear.exnode               05-Jun-2003  2.6k
shear.ipbase               05-Jun-2003  5.6k
shear.ipcoor               05-Jun-2003  570 
shear.ipequa               02-May-2004  2.0k
shear.ipopti               05-Jul-2005  1.5k
shear.ippara               05-Jun-2003  5.9k
shear.ipsolv               16-Aug-2010  2.7k
shear.ipsolv.old           13-Apr-2007  2.5k
shear_NF.ipdata            05-Jun-2003  117 
shear_NF.ipinit            05-Jun-2003  1.7k
shear_def.exelem           05-Jun-2003  6.8k
shear_def.exnode           05-Jun-2003  2.6k
shear_def.opresi           05-Jun-2003    0 
shear_ref.ipelem           05-Jun-2003  802 
shear_ref.ipelfb           05-Jun-2003  688 
shear_ref.ipfibr           05-Jun-2003  2.8k
shear_specimen_ref.ipnode  05-Jun-2003  3.0k
test_output.com            05-Jun-2003   80 

Download the entire example:

Name                      Modified     Size

examples_7_72_723.tar.gz  17-Aug-2010   16k

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:31 2016

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