Class to manage nested solution for augmented Lagrange contact. More...

#include <AugmentedLagrangianContactProblem.h>

Inheritance diagram for AugmentedLagrangianContactProblem:

Public Member Functions
	AugmentedLagrangianContactProblem (const InputParameters &params)

virtual	~AugmentedLagrangianContactProblem ()

virtual void	timestepSetup () override

virtual MooseNonlinearConvergenceReason	checkNonlinearConvergence (std::string &msg, const PetscInt it, const Real xnorm, const Real snorm, const Real fnorm, const Real rtol, const Real divtol, const Real stol, const Real abstol, const PetscInt nfuncs, const PetscInt max_funcs, const PetscBool force_iteration, const Real ref_resid, const Real div_threshold) override

Private Attributes
int	_num_lagmul_iterations

int	_max_lagmul_iters

Detailed Description

Class to manage nested solution for augmented Lagrange contact.

The AugmentedLagrangianContactProblem manages the nested solution procedure, repeating the solution until convergence has been achieved, checking for convergence, and updating the Lagrangian multipliers.

Definition at line 28 of file AugmentedLagrangianContactProblem.h.

Constructor & Destructor Documentation

◆ AugmentedLagrangianContactProblem()

AugmentedLagrangianContactProblem::AugmentedLagrangianContactProblem ( const InputParameters & params )

Definition at line 41 of file AugmentedLagrangianContactProblem.C.

   : ReferenceResidualProblem(params),
     _num_lagmul_iterations(0),
     _max_lagmul_iters(getParam<int>("maximum_lagrangian_update_iterations"))
 {
 }

◆ ~AugmentedLagrangianContactProblem()

virtual AugmentedLagrangianContactProblem::~AugmentedLagrangianContactProblem ( )

inlinevirtual

Definition at line 32 of file AugmentedLagrangianContactProblem.h.

32 {}

Member Function Documentation

◆ checkNonlinearConvergence()

MooseNonlinearConvergenceReason AugmentedLagrangianContactProblem::checkNonlinearConvergence	(	std::string &	msg,
		const PetscInt	it,
		const Real	xnorm,
		const Real	snorm,
		const Real	fnorm,
		const Real	rtol,
		const Real	divtol,
		const Real	stol,
		const Real	abstol,
		const PetscInt	nfuncs,
		const PetscInt	max_funcs,
		const PetscBool	force_iteration,
		const Real	ref_resid,
		const Real	div_threshold
	)

overridevirtual

Definition at line 56 of file AugmentedLagrangianContactProblem.C.

 {
  
   Real my_max_funcs = std::numeric_limits<int>::max();
   Real my_div_threshold = std::numeric_limits<Real>::max();
  
   MooseNonlinearConvergenceReason reason =
       ReferenceResidualProblem::checkNonlinearConvergence(msg,
                                                           it,
                                                           xnorm,
                                                           snorm,
                                                           fnorm,
                                                           rtol,
                                                           divtol,
                                                           stol,
                                                           abstol,
                                                           nfuncs,
                                                           my_max_funcs,
                                                           force_iteration,
                                                           ref_resid,
                                                           my_div_threshold);
  
   _console << "Augmented Lagrangian contact iteration " << _num_lagmul_iterations << "\n";
  
   bool _augLM_repeat_step;
  
   if (reason == MooseNonlinearConvergenceReason::CONVERGED_FNORM_ABS ||
       reason == MooseNonlinearConvergenceReason::CONVERGED_FNORM_RELATIVE ||
       reason == MooseNonlinearConvergenceReason::CONVERGED_SNORM_RELATIVE)
   {
     if (_num_lagmul_iterations < _max_lagmul_iters)
     {
       NonlinearSystemBase & nonlinear_sys = getNonlinearSystemBase();
       nonlinear_sys.update();
  
       const ConstraintWarehouse & constraints = nonlinear_sys.getConstraintWarehouse();
  
       std::map<std::pair<unsigned int, unsigned int>, PenetrationLocator *> * penetration_locators =
           NULL;
  
       bool displaced = false;
       _augLM_repeat_step = false;
       if (getDisplacedProblem() == NULL)
       {
         GeometricSearchData & geom_search_data = geomSearchData();
         penetration_locators = &geom_search_data._penetration_locators;
       }
       else
       {
         GeometricSearchData & displaced_geom_search_data = getDisplacedProblem()->geomSearchData();
         penetration_locators = &displaced_geom_search_data._penetration_locators;
         displaced = true;
       }
  
       for (const auto & it : *penetration_locators)
       {
         PenetrationLocator & pen_loc = *(it.second);
  
         BoundaryID slave_boundary = pen_loc._slave_boundary;
  
         if (constraints.hasActiveNodeFaceConstraints(slave_boundary, displaced))
         {
           const auto & ncs = constraints.getActiveNodeFaceConstraints(slave_boundary, displaced);
  
           for (const auto & nc : ncs)
           {
             if (std::dynamic_pointer_cast<MechanicalContactConstraint>(nc) == NULL)
               mooseError("AugmentedLagrangianContactProblem: dynamic cast of "
                          "MechanicalContactConstraint object failed.");
  
             if (!(std::dynamic_pointer_cast<MechanicalContactConstraint>(nc))
                      ->AugmentedLagrangianContactConverged())
             {
               (std::dynamic_pointer_cast<MechanicalContactConstraint>(nc))
                   ->updateAugmentedLagrangianMultiplier(false);
               _augLM_repeat_step = true;
               break;
             }
           }
         }
       }
  
       if (_augLM_repeat_step)
       {
         // force it to keep iterating
         reason = MooseNonlinearConvergenceReason::ITERATING;
         _console << "Augmented Lagrangian Multiplier needs updating." << std::endl;
         _num_lagmul_iterations++;
       }
       else
         _console << "Augmented Lagrangian contact constraint enforcement is satisfied."
                  << std::endl;
     }
     else
     {
       // maxed out
       _console << "Maximum Augmented Lagrangian contact iterations have been reached." << std::endl;
       reason = MooseNonlinearConvergenceReason::DIVERGED_FUNCTION_COUNT;
     }
   }
  
   return reason;
 }

◆ timestepSetup()

void AugmentedLagrangianContactProblem::timestepSetup ( )

overridevirtual

Definition at line 49 of file AugmentedLagrangianContactProblem.C.

 {
   _num_lagmul_iterations = 0;
   ReferenceResidualProblem::timestepSetup();
 }

Member Data Documentation

◆ _max_lagmul_iters

int AugmentedLagrangianContactProblem::_max_lagmul_iters

private

Definition at line 54 of file AugmentedLagrangianContactProblem.h.

Referenced by checkNonlinearConvergence().

◆ _num_lagmul_iterations

int AugmentedLagrangianContactProblem::_num_lagmul_iterations

private

Definition at line 53 of file AugmentedLagrangianContactProblem.h.

Referenced by checkNonlinearConvergence(), and timestepSetup().

The documentation for this class was generated from the following files:

contact/include/problems/AugmentedLagrangianContactProblem.h
contact/src/problems/AugmentedLagrangianContactProblem.C

Public Member Functions

Private Attributes