Managing N.-R. Iterations

Mechanical iterations parameters are imposed by functions associated to metafor.getMechanicalIterationManager(), when thermal iterations parameters are imposed by functions associated to metafor.getThermalIterationManager().

  mim = metafor.getMechanicalIterationManager() 
  tim = metafor.getThermalIterationManager()
  
  mim/tim.setResidualComputationMethod(meth)    
  mim/tim.setVerbose()     
  mim/tim.setResidualTolerance(prec)
  mim/tim.setMaxNbOfIterations(itma)
  mim/tim.setNbOfIterationsForStiffUpdate(irea)
  mim/tim.setMaxNbOfLineSearchIterations(lsma)
  mim/tim.setLineSearchTolerance1(lsp1)
  mim/tim.setLineSearchTolerance2(lsp2)
  mim.setPredictorComputationMethod(predMeth)
  mim.setConstrainedDofsVAComputationMethod(VACMeth)
  mim/tim.setPreviousResidualsDecreaseCriterion()
  mim/tim.setCpuDependency(cpuDep)
  mim/tim.setForceOneIte(bool)  

where meth defines the method used to compute the $Residual$.

The residual is adimensionalized by the norm of external forces and internal reactions . First :

$$Rmin = \frac{||FreeExternalForces|| + ||CstrInternalForces||+||CstrInertialForces||}{nreac} $$

If $$ Rmin < limitNormFactor $$

$Rmin$ is set to $$ Rmin = limitNormFactor $$

and the residual is adimensionalized :

$$Residual = \frac{||FreeResidual||}{ndofs*Rmin}$$

where

$ndofs$ : total number of degrees of freedom

$nreac$ : number of fixed degrees of freedom (having an associated reaction)

The residual is adimensionalized by the square of the norm of external forces and internal reactions. First:

$$Rmin = \frac{\sqrt{||FreeExternalForces||^2 + ||CstrInternalForces||^2+||CstrInertialForces||^2}}{nreac}$$

If $$ Rmin < limitNormFactor $$ then $$ Rmin = limitNormFactor $$

and the residual is adimensionalized :

$$Residual = \frac{||FreeResidual||}{ndofs*Rmin}$$

where

$ndofs$ : total number of degrees of freedom

$nreac$ : number of fixed degrees of freedom (having an associated reaction)

The residual is adimensionalized by the norm of external and inertial forces, and internal reactions (value to use preferentially for thermal simulations. Careful, the user must define it).

If $$ ||FreeExternalForces|| < limitNormFactor $$ then $$ Residual = max(|FreeResidual|) $$ Else, $$Rmin = ||FreeExternalForces|| + ||CstrInternalForces||+||CstrInertialForces|| $$ and $$Residual = \frac{||FreeResidual||}{ndofs*Rmin}$$

This method is the same as Method1ResidualComputation() but without dividing by the number of dofs nor nreac, so

$$Rmin = ||FreeExternalForces|| + ||CstrInternalForces||+||CstrInertialForces|| $$

If $$ Rmin < limitNormFactor$$ then $$ Rmin = limitNormFactor $$

and the residual is adimensionalized :

$$Residual = \frac{||FreeResidual||}{Rmin}$$

The residual is adimensionalized by the factor adimFactor, given by the user, without dividing by the number of dofs.

$$Residual = \frac{||FreeResidual||}{adimFactor }$$

The residual is adimensionalized by the factor adimFactor, given by the user, and by the number of dofs. $$Residual = \frac{||FreeResidual||}{ndofs * adimFactor }$$

where

$ndofs$ : total number of dofs