CappedWeakInclinedPlaneStressUpdate.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 "CappedWeakPlaneStressUpdate.h"
 
class CappedWeakInclinedPlaneStressUpdate : public CappedWeakPlaneStressUpdate
{
public:
  static InputParameters validParams();
 
  CappedWeakInclinedPlaneStressUpdate(const InputParameters & parameters);
 
  bool requiresIsotropicTensor() override { return false; }
 
protected:
  virtual void initQpStatefulProperties() override;
 
  RealVectorValue _n_input;
 
  MaterialProperty<RealVectorValue> & _n;
 
  const MaterialProperty<RealVectorValue> & _n_old;
 
  RealTensorValue _rot_n_to_z;
 
  RealTensorValue _rot_z_to_n;
 
  RankTwoTensor _rotated_trial;
 
  RankFourTensor _rotated_Eijkl;
 
  virtual void initializeReturnProcess() override;
  virtual void finalizeReturnProcess(const RankTwoTensor & rotation_increment) override;
 
  virtual void preReturnMap(Real p_trial,
                            Real q_trial,
                            const RankTwoTensor & stress_trial,
                            const std::vector<Real> & intnl_old,
                            const std::vector<Real> & yf,
                            const RankFourTensor & Eijkl) override;
 
  virtual void computePQ(const RankTwoTensor & stress, Real & p, Real & q) const override;
 
  virtual void setEppEqq(const RankFourTensor & Eijkl, Real & Epp, Real & Eqq) const override;
 
  virtual void setStressAfterReturn(const RankTwoTensor & stress_trial,
                                    Real p_ok,
                                    Real q_ok,
                                    Real gaE,
                                    const std::vector<Real> & intnl,
                                    const yieldAndFlow & smoothed_q,
                                    const RankFourTensor & Eijkl,
                                    RankTwoTensor & stress) const override;
 
  virtual void consistentTangentOperator(const RankTwoTensor & stress_trial,
                                         Real p_trial,
                                         Real q_trial,
                                         const RankTwoTensor & stress,
                                         Real p,
                                         Real q,
                                         Real gaE,
                                         const yieldAndFlow & smoothed_q,
                                         const RankFourTensor & Eijkl,
                                         bool compute_full_tangent_operator,
                                         RankFourTensor & cto) const override;
 
  virtual RankTwoTensor dpdstress(const RankTwoTensor & stress) const override;
 
  virtual RankTwoTensor dqdstress(const RankTwoTensor & stress) const override;
 
  virtual RankFourTensor d2qdstress2(const RankTwoTensor & stress) const override;
};