CappedMohrCoulombCosseratStressUpdate.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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 "CappedMohrCoulombStressUpdate.h"
 
class CappedMohrCoulombCosseratStressUpdate : public CappedMohrCoulombStressUpdate
{
public:
  static InputParameters validParams();
 
  CappedMohrCoulombCosseratStressUpdate(const InputParameters & parameters);
 
  bool requiresIsotropicTensor() override { return false; }
 
protected:
  const Real _host_young;
 
  const Real _host_poisson;
 
  const Real _host_E0011;
 
  const Real _host_E0000;
 
  virtual void setStressAfterReturnV(const RankTwoTensor & stress_trial,
                                     const std::vector<Real> & stress_params,
                                     Real gaE,
                                     const std::vector<Real> & intnl,
                                     const yieldAndFlow & smoothed_q,
                                     const RankFourTensor & Eijkl,
                                     RankTwoTensor & stress) const override;
 
  virtual void preReturnMapV(const std::vector<Real> & trial_stress_params,
                             const RankTwoTensor & stress_trial,
                             const std::vector<Real> & intnl_old,
                             const std::vector<Real> & yf,
                             const RankFourTensor & Eijkl) override;
 
  void setEffectiveElasticity(const RankFourTensor & Eijkl) override;
 
  virtual void consistentTangentOperatorV(const RankTwoTensor & stress_trial,
                                          const std::vector<Real> & trial_stress_params,
                                          const RankTwoTensor & stress,
                                          const std::vector<Real> & stress_params,
                                          Real gaE,
                                          const yieldAndFlow & smoothed_q,
                                          const RankFourTensor & Eijkl,
                                          bool compute_full_tangent_operator,
                                          const std::vector<std::vector<Real>> & dvar_dtrial,
                                          RankFourTensor & cto) override;
 
private:
  std::vector<Real> _eigvals_scratch; // used as temporary storage to avoid heap allocation
};