ReturnMappingModel.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
#include "ConstitutiveModel.h"
#include "SingleVariableReturnMappingSolution.h"
 
class ReturnMappingModel;
 
template <>
InputParameters validParams<ReturnMappingModel>();
 
class ReturnMappingModel : public ConstitutiveModel, public SingleVariableReturnMappingSolution
{
public:
  ReturnMappingModel(const InputParameters & parameters,
                     const std::string inelastic_strain_name = "");
  virtual ~ReturnMappingModel() {}
 
  virtual void initQpStatefulProperties() override;
 
  virtual void computeStress(const Elem & current_elem,
                             const SymmElasticityTensor & elasticityTensor,
                             const SymmTensor & stress_old,
                             SymmTensor & strain_increment,
                             SymmTensor & stress_new) override;
 
  void computeStress(const Elem & current_elem,
                     const SymmElasticityTensor & elasticityTensor,
                     const SymmTensor & stress_old,
                     SymmTensor & strain_increment,
                     SymmTensor & stress_new,
                     SymmTensor & inelastic_strain_increment);
 
  virtual Real computeReferenceResidual(const Real effective_trial_stress,
                                        const Real scalar) override;
 
  Real computeTimeStepLimit();
 
  virtual Real minimumPermissibleValue(const Real /*effective_trial_stress*/) const override
  {
    return 0.0;
  }
 
  void outputIterationSummary(std::stringstream * iter_output,
                              const unsigned int total_it) override;
 
protected:
  virtual void computeStressInitialize(Real /*effectiveTrialStress*/,
                                       const SymmElasticityTensor & /*elasticityTensor*/)
  {
  }
 
  virtual void computeStressFinalize(const SymmTensor & /*inelasticStrainIncrement*/) {}
 
  Real _effective_strain_increment;
 
  Real _three_shear_modulus;
 
  MaterialProperty<Real> & _effective_inelastic_strain;
  const MaterialProperty<Real> & _effective_inelastic_strain_old;
  Real _max_inelastic_increment;
  const bool _compute_matl_timestep_limit;
  MaterialProperty<Real> * _matl_timestep_limit;
};