#include "users_modules/eshelbian_plasticity/src/PlasticIncrementalOptimization.hpp"

Inheritance diagram for EshelbianPlasticity::PlasticIncrementalOptimizationProblem:

Collaboration diagram for EshelbianPlasticity::PlasticIncrementalOptimizationProblem:

Public Member Functions
	PlasticIncrementalOptimizationProblem (EshelbianCore &ep, SmartPetscObj< TS > ts)

std::string	getName () const override

MoFEMErrorCode	validate () override

DM	getControlDM () const override

MoFEMErrorCode	createControlData (SmartPetscObj< Vec > &committed_control, SmartPetscObj< Vec > &candidate_control) override

MoFEMErrorCode	configureConstraintDM (DM constraint_dm) override

MoFEMErrorCode	initialiseInequalityMultipliers (DM constraint_dm, Vec multipliers) override

MoFEMErrorCode	initialiseReferenceControl (Vec committed_control, Vec reference_control) override

MoFEMErrorCode	solveState (Vec state, EquilibriumSolveReport &report) override

MoFEMErrorCode	restoreState (Vec state) override

MoFEMErrorCode	evaluateConservativeValue (double &value) override

MoFEMErrorCode	assembleConservativeGradient (Vec smooth_gradient) override

MoFEMErrorCode	evaluateDissipation (Vec control, double &value) override

MoFEMErrorCode	assembleObjectiveGradient (Vec control, Vec smooth_gradient, Vec objective_gradient) override

MoFEMErrorCode	evaluateInequalityConstraints (Vec control, Vec constraints) override

MoFEMErrorCode	assembleInequalityJacobian (Vec control, Vec smooth_gradient, Mat jacobian) override

MoFEMErrorCode	assembleInequalityJacobianWithoutSmoothGradient (Vec control, Mat jacobian) override

MoFEMErrorCode	validateSolution (DM constraint_dm, Vec control, Vec smooth_gradient, Vec inequality_multipliers, PetscReal gradient_tolerance, PetscReal constraint_tolerance, std::string &diagnostics) override

MoFEMErrorCode	commit (Vec solution, Vec equilibrated_state, Vec committed_control, Vec candidate_control, Vec reference_control) override

const std::string &	getPlasticFlowFieldName () const

const std::string &	getPlasticKappaFieldName () const

double	getInitialYieldStress () const

double	getIsotropicHardeningModulus () const

double	getDissipationRegularizationEpsilon () const

const std::string &	getVolumeElementName () const

MoFEMErrorCode	addConstitutiveGeometryOperators (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pipeline) const

Public Member Functions inherited from EshelbianPlasticity::IncrementalOptimizationProblem
virtual	~IncrementalOptimizationProblem ()=default

Private Attributes
EshelbianCore &	ep

SmartPetscObj< TS >	timeSolver

SmartPetscObj< Vec >	fallbackJacobianSmoothGradient

double	initialYieldStress = 1

double	isotropicHardeningModulus = 0

double	dissipationRegularizationEpsilon = 0

Detailed Description

Plasticity adapter; no plastic data is present in the generic context.

Definition at line 20 of file PlasticIncrementalOptimization.hpp.

Constructor & Destructor Documentation

◆ PlasticIncrementalOptimizationProblem()

EshelbianPlasticity::PlasticIncrementalOptimizationProblem::PlasticIncrementalOptimizationProblem	(	EshelbianCore &	ep,
		SmartPetscObj< TS >	ts
	)

Definition at line 148 of file PlasticIncrementalOptimization.cpp.

150 : ep(ep), timeSolver(std::move(ts)) {}

EshelbianPlasticity::PlasticIncrementalOptimizationProblem::timeSolver

SmartPetscObj< TS > timeSolver

Definition PlasticIncrementalOptimization.hpp:78

EshelbianPlasticity::PlasticIncrementalOptimizationProblem::ep

EshelbianCore & ep

Definition PlasticIncrementalOptimization.hpp:77

Member Function Documentation

◆ addConstitutiveGeometryOperators()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::addConstitutiveGeometryOperators ( boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > & pipeline ) const

Add geometry operators for local constitutive integration.

Definition at line 176 of file PlasticIncrementalOptimization.cpp.

                                                                          {
  // Only transformed quadrature weights are needed. Keep the P0 bases
  // unchanged while evaluating vector-owned trial coefficients.
  return MoFEM::AddHOOps<SPACE_DIM, SPACE_DIM, SPACE_DIM>::add(
      pipeline, {NOSPACE}, ep.materialH1Positions);
}

◆ assembleConservativeGradient()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::assembleConservativeGradient ( Vec smooth_gradient )

overridevirtual

Assemble the reduced conservative gradient.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 299 of file PlasticIncrementalOptimization.cpp.

                         {
  return assembleReducedHelmholtzGradient(ep, timeSolver, smooth_gradient);
}

◆ assembleInequalityJacobian()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::assembleInequalityJacobian	(	Vec	control,
		Vec	smooth_gradient,
		Mat	jacobian
	)

overridevirtual

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 325 of file PlasticIncrementalOptimization.cpp.

                                                    {
  return assemblePlasticInequalityJacobian(*this, control, smooth_gradient,
                                           jacobian);
}

◆ assembleInequalityJacobianWithoutSmoothGradient()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::assembleInequalityJacobianWithoutSmoothGradient	(	Vec	control,
		Mat	jacobian
	)

overridevirtual

Assemble a finite generalized row without a usable smooth gradient.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 331 of file PlasticIncrementalOptimization.cpp.

                                                                   {
  MoFEMFunctionBegin;
  if (!fallbackJacobianSmoothGradient)
    SETERRQ(ep.mField.get_comm(), MOFEM_DATA_INCONSISTENCY,
            "Fallback plastic Jacobian scratch vector is unavailable");
  CHKERR VecZeroEntries(fallbackJacobianSmoothGradient);
  // At the exact cone apex, zero is an admissible subgradient. Away from the
  // apex, and on the rounded branch, the row depends only on the control.
  CHKERR assemblePlasticInequalityJacobian(
      *this, control, fallbackJacobianSmoothGradient, jacobian);
  MoFEMFunctionReturn(0);
}

◆ assembleObjectiveGradient()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::assembleObjectiveGradient	(	Vec	control,
		Vec	smooth_gradient,
		Vec	objective_gradient
	)

overridevirtual

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 312 of file PlasticIncrementalOptimization.cpp.

                                                              {
  return assembleIncrementalObjectiveGradient(*this, control, smooth_gradient,
                                              objective_gradient);
}

◆ commit()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::commit	(	Vec	solution,
		Vec	equilibrated_state,
		Vec	committed_control,
		Vec	candidate_control,
		Vec	reference_control
	)

overridevirtual

Commit physics history and all next-step control vectors atomically. solution and equilibrated_state are read-only accepted inputs; a successful implementation may modify only physics history and the three control outputs.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 740 of file PlasticIncrementalOptimization.cpp.

                                                  {
  return commitIncrementalOptimizationLoadStep(
      ep, solution, equilibrated_state, committed_control, candidate_control,
      reference_control);
}

◆ configureConstraintDM()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::configureConstraintDM ( DM constraint_dm )

overridevirtual

Add physics-specific rows, columns, and elements to constraintDM.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 248 of file PlasticIncrementalOptimization.cpp.

                      {
  MoFEMFunctionBegin;
  if (!constraint_dm)
    SETERRQ(ep.mField.get_comm(), MOFEM_DATA_INCONSISTENCY,
            "Plastic constraint DM is null");
  CHKERR DMMoFEMAddSubFieldRow(constraint_dm, ep.plasticKappaField,
                               ep.plasticVolumes);
  CHKERR DMMoFEMAddSubFieldCol(constraint_dm, ep.plasticFlowField,
                               ep.plasticVolumes);
  CHKERR DMMoFEMAddSubFieldCol(constraint_dm, ep.plasticKappaField,
                               ep.plasticVolumes);
  CHKERR DMMoFEMAddElement(constraint_dm, ep.elementVolumeName);
  MoFEMFunctionReturn(0);
}

◆ createControlData()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::createControlData	(	SmartPetscObj< Vec > &	committed_control,
		SmartPetscObj< Vec > &	candidate_control
	)

overridevirtual

Create the problem's committed and candidate control vectors.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 199 of file PlasticIncrementalOptimization.cpp.

                                           {
  MoFEMFunctionBegin;
  committed_control =
      createDMVector(ep.dmIncrementalOptimization, RowColData::ROW);
  if (!ep.incrementalTrialControl)
    SETERRQ(ep.mField.get_comm(), MOFEM_DATA_INCONSISTENCY,
            "Plastic trial-control vector is not available");
  candidate_control = ep.incrementalTrialControl;
  // Mesh field data store committed kappa, whereas every vector with this DM
  // layout stores fresh increments (Delta H_p, Delta kappa).
  CHKERR VecZeroEntries(committed_control);
  CHKERR VecZeroEntries(candidate_control);
  fallbackJacobianSmoothGradient = vectorDuplicate(committed_control);
  CHKERR VecZeroEntries(fallbackJacobianSmoothGradient);
  CHKERR PetscOptionsGetReal(PETSC_NULLPTR, "",
                             "-incremental_optimization_yield_stress",
                             &initialYieldStress, PETSC_NULLPTR);
  CHKERR PetscOptionsGetReal(
      PETSC_NULLPTR, "",
      "-incremental_optimization_isotropic_hardening_modulus",
      &isotropicHardeningModulus, PETSC_NULLPTR);
  CHKERR PetscOptionsGetReal(PETSC_NULLPTR, "",
                             "-incremental_optimization_dissipation_epsilon",
                             &dissipationRegularizationEpsilon, PETSC_NULLPTR);
  if (!std::isfinite(initialYieldStress) || !(initialYieldStress > 0))
    SETERRQ(ep.mField.get_comm(), MOFEM_INVALID_DATA,
            "Initial plastic yield stress must be finite and positive");
  if (!std::isfinite(isotropicHardeningModulus) ||
      isotropicHardeningModulus < 0)
    SETERRQ(ep.mField.get_comm(), MOFEM_INVALID_DATA,
            "Isotropic hardening modulus must be finite and non-negative");
  if (!std::isfinite(dissipationRegularizationEpsilon) ||
      dissipationRegularizationEpsilon < 0)
    SETERRQ(ep.mField.get_comm(), MOFEM_INVALID_DATA,
            "Plastic dissipation epsilon must be finite and non-negative");
  PlasticKappaHistory kappa_history;
  CHKERR savePlasticKappaHistory(ep, kappa_history);
  MOFEM_LOG("EP", Sev::inform)
      << "Plastic model parameters: initial yield stress "
      << initialYieldStress << ", isotropic hardening modulus "
      << isotropicHardeningModulus
      << ", dissipation epsilon " << dissipationRegularizationEpsilon
      << ", constraint scaling sqrt(reference cell volume)";
 
  MoFEMFunctionReturn(0);
}

◆ evaluateConservativeValue()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::evaluateConservativeValue ( double & value )

overridevirtual

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 291 of file PlasticIncrementalOptimization.cpp.

                   {
  MoFEMFunctionBegin;
  CHKERR evaluateEquilibratedMechanicalValue(ep, timeSolver, value);
  MoFEMFunctionReturn(0);
}

◆ evaluateDissipation()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::evaluateDissipation	(	Vec	control,
		double &	value
	)

overridevirtual

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 305 of file PlasticIncrementalOptimization.cpp.

                                                                          {
  // The generic "dissipation" slot is the complete local incremental
  // resistance: irreversible yield work plus stored hardening-energy change.
  return evaluatePlasticIncrementalResistanceValue(*this, control, value);
}

◆ evaluateInequalityConstraints()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::evaluateInequalityConstraints	(	Vec	control,
		Vec	constraints
	)

overridevirtual

Return sqrt(cell volume) times the exact or rounded epigraph gap.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 319 of file PlasticIncrementalOptimization.cpp.

                                  {
  return evaluatePlasticInequalityConstraints(*this, control, constraints);
}

◆ getControlDM()

DM EshelbianPlasticity::PlasticIncrementalOptimizationProblem::getControlDM ( ) const

overridevirtual

Parent DM containing every optimization-control field.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 156 of file PlasticIncrementalOptimization.cpp.

                                                             {
  return ep.dmIncrementalOptimization;
}

◆ getDissipationRegularizationEpsilon()

double EshelbianPlasticity::PlasticIncrementalOptimizationProblem::getDissipationRegularizationEpsilon ( ) const

inline

Definition at line 68 of file PlasticIncrementalOptimization.hpp.

                                                     {
    return dissipationRegularizationEpsilon;
  }

◆ getInitialYieldStress()

double EshelbianPlasticity::PlasticIncrementalOptimizationProblem::getInitialYieldStress ( ) const

inline

Definition at line 64 of file PlasticIncrementalOptimization.hpp.

64{ return initialYieldStress; }

◆ getIsotropicHardeningModulus()

double EshelbianPlasticity::PlasticIncrementalOptimizationProblem::getIsotropicHardeningModulus ( ) const

inline

Definition at line 65 of file PlasticIncrementalOptimization.hpp.

                                              {
    return isotropicHardeningModulus;
  }

◆ getName()

std::string EshelbianPlasticity::PlasticIncrementalOptimizationProblem::getName ( ) const

overridevirtual

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 152 of file PlasticIncrementalOptimization.cpp.

                                                               {
  return "plasticity";
}

◆ getPlasticFlowFieldName()

const std::string & EshelbianPlasticity::PlasticIncrementalOptimizationProblem::getPlasticFlowFieldName ( ) const

Plastic-flow control field; its vector values are Delta H^p.

Definition at line 161 of file PlasticIncrementalOptimization.cpp.

                                                                     {
  return ep.plasticFlowField;
}

◆ getPlasticKappaFieldName()

const std::string & EshelbianPlasticity::PlasticIncrementalOptimizationProblem::getPlasticKappaFieldName ( ) const

Scalar control/constraint field; its vector values are Delta kappa.

Definition at line 166 of file PlasticIncrementalOptimization.cpp.

                                                                      {
  return ep.plasticKappaField;
}

◆ getVolumeElementName()

const std::string & EshelbianPlasticity::PlasticIncrementalOptimizationProblem::getVolumeElementName ( ) const

Definition at line 171 of file PlasticIncrementalOptimization.cpp.

                                                                  {
  return ep.elementVolumeName;
}

◆ initialiseInequalityMultipliers()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::initialiseInequalityMultipliers	(	DM	constraint_dm,
		Vec	multipliers
	)

overridevirtual

Initialise multipliers for sqrt(cell volume)-scaled inequality rows.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 265 of file PlasticIncrementalOptimization.cpp.

                                       {
  return initialisePlasticInequalityMultipliers(*this, constraint_dm,
                                                multipliers);
}

◆ initialiseReferenceControl()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::initialiseReferenceControl	(	Vec	committed_control,
		Vec	reference_control
	)

overridevirtual

Define the zero-increment/reference point for a new load step.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 272 of file PlasticIncrementalOptimization.cpp.

                                                  {
  MoFEMFunctionBegin;
  (void)committed_control;
  // Both Delta H_p and Delta kappa are fresh increment variables. Committed
  // kappa remains in mesh field storage.
  CHKERR VecZeroEntries(reference_control);
  MoFEMFunctionReturn(0);
}

◆ restoreState()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::restoreState ( Vec state )

overridevirtual

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 287 of file PlasticIncrementalOptimization.cpp.

                                                                            {
  return setStateOnMesh(ep, state);
}

◆ solveState()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::solveState	(	Vec	state,
		EquilibriumSolveReport &	report
	)

overridevirtual

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 282 of file PlasticIncrementalOptimization.cpp.

                                               {
  return solveEquilibriumStateTS(ep, timeSolver, state, report);
}

◆ validate()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::validate ( )

overridevirtual

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 184 of file PlasticIncrementalOptimization.cpp.

                                                               {
  MoFEMFunctionBegin;
  if (!ep.plasticVolume || ep.interfaceCrack)
    SETERRQ(ep.mField.get_comm(), MOFEM_NOT_IMPLEMENTED,
            "Plastic incremental optimization requires -plastic_volume 1 "
            "and -cohesive_interface_on 0");
  if (!ep.dmIncrementalOptimization)
    SETERRQ(ep.mField.get_comm(), MOFEM_DATA_INCONSISTENCY,
            "Plastic incremental-optimization control DM is not available");
  if (!timeSolver)
    SETERRQ(ep.mField.get_comm(), MOFEM_DATA_INCONSISTENCY,
            "Plastic incremental optimization requires a state TS");
  MoFEMFunctionReturn(0);
}

◆ validateSolution()

MoFEMErrorCode EshelbianPlasticity::PlasticIncrementalOptimizationProblem::validateSolution	(	DM	constraint_dm,
		Vec	control,
		Vec	smooth_gradient,
		Vec	inequality_multipliers,
		PetscReal	gradient_tolerance,
		PetscReal	constraint_tolerance,
		std::string &	diagnostics
	)

overridevirtual

Recover physical multipliers and check unscaled plastic feasibility.

Implements EshelbianPlasticity::IncrementalOptimizationProblem.

Definition at line 731 of file PlasticIncrementalOptimization.cpp.

                                                            {
  return validatePlasticIncrementalSolution(
      *this, constraint_dm, control, smooth_gradient, inequality_multipliers,
      gradient_tolerance, constraint_tolerance, diagnostics);
}

Member Data Documentation

◆ dissipationRegularizationEpsilon

double EshelbianPlasticity::PlasticIncrementalOptimizationProblem::dissipationRegularizationEpsilon = 0

private

Definition at line 83 of file PlasticIncrementalOptimization.hpp.

◆ ep

EshelbianCore& EshelbianPlasticity::PlasticIncrementalOptimizationProblem::ep

private

Definition at line 77 of file PlasticIncrementalOptimization.hpp.

◆ fallbackJacobianSmoothGradient

SmartPetscObj<Vec> EshelbianPlasticity::PlasticIncrementalOptimizationProblem::fallbackJacobianSmoothGradient

private

Zero force selects a finite exact-cone row when no smooth force exists.

Definition at line 80 of file PlasticIncrementalOptimization.hpp.

◆ initialYieldStress

double EshelbianPlasticity::PlasticIncrementalOptimizationProblem::initialYieldStress = 1

private

Definition at line 81 of file PlasticIncrementalOptimization.hpp.

◆ isotropicHardeningModulus

double EshelbianPlasticity::PlasticIncrementalOptimizationProblem::isotropicHardeningModulus = 0

private

Definition at line 82 of file PlasticIncrementalOptimization.hpp.

◆ timeSolver

SmartPetscObj<TS> EshelbianPlasticity::PlasticIncrementalOptimizationProblem::timeSolver

private

Definition at line 78 of file PlasticIncrementalOptimization.hpp.

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

users_modules/eshelbian_plasticity/src/PlasticIncrementalOptimization.hpp
users_modules/eshelbian_plasticity/src/impl/PlasticIncrementalOptimization.cpp

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PlasticIncrementalOptimizationProblem()

Member Function Documentation

◆ addConstitutiveGeometryOperators()

◆ assembleConservativeGradient()

◆ assembleInequalityJacobian()

◆ assembleInequalityJacobianWithoutSmoothGradient()

◆ assembleObjectiveGradient()

◆ commit()

◆ configureConstraintDM()

◆ createControlData()

◆ evaluateConservativeValue()

◆ evaluateDissipation()

◆ evaluateInequalityConstraints()

◆ getControlDM()

◆ getDissipationRegularizationEpsilon()

◆ getInitialYieldStress()

◆ getIsotropicHardeningModulus()

◆ getName()

◆ getPlasticFlowFieldName()

◆ getPlasticKappaFieldName()

◆ getVolumeElementName()

◆ initialiseInequalityMultipliers()

◆ initialiseReferenceControl()

◆ restoreState()

◆ solveState()

◆ validate()

◆ validateSolution()

Member Data Documentation

◆ dissipationRegularizationEpsilon

◆ ep

◆ fallbackJacobianSmoothGradient

◆ initialYieldStress

◆ isotropicHardeningModulus

◆ timeSolver