RateDepSmearCrackModel.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 "SymmElasticityTensor.h"
 
#include "petscsys.h"
#include "petscblaslapack.h"
 
#if PETSC_VERSION_LESS_THAN(3, 5, 0)
extern "C" void FORTRAN_CALL(dgetri)(...); // matrix inversion routine from LAPACK
#endif
 
class RateDepSmearCrackModel;
 
template <>
InputParameters validParams<RateDepSmearCrackModel>();
 
class RateDepSmearCrackModel : public ConstitutiveModel
{
public:
  RateDepSmearCrackModel(const InputParameters & parameters);
 
  virtual ~RateDepSmearCrackModel();
 
protected:
  virtual void computeStress(const Elem & current_elem,
                             const SymmElasticityTensor & elasticity_tensor,
                             const SymmTensor & stress_old,
                             SymmTensor & strain_increment,
                             SymmTensor & stress_new);
 
  virtual void initQpStatefulProperties();
 
  virtual void initVariables();
 
  virtual void solve();
  virtual void postSolveVariables();
  virtual void postSolveStress();
 
  virtual void updateVariables();
 
  bool getConvergeVar();
 
  virtual void calcResidual();
 
  virtual void calcStateIncr();
 
  virtual void calcJacobian();
 
  int matrixInversion(std::vector<Real> & A, int n) const;
 
  Real _ref_damage_rate; 
  unsigned int _nstate;  
  Real _exponent;
  unsigned int _maxiter; 
  Real _tol;             
  Real _zero_tol;        
  Real _intvar_incr_tol; 
  bool _input_rndm_scale_var; 
  Real _rndm_scale_var;       
 
  MaterialProperty<std::vector<Real>> & _intvar;
  const MaterialProperty<std::vector<Real>> & _intvar_old;
 
  MaterialProperty<SymmTensor> & _stress_undamaged;
  const MaterialProperty<SymmTensor> & _stress_undamaged_old;
 
  std::vector<Real> _intvar_incr;
  std::vector<Real> _intvar_tmp, _intvar_old_tmp;
  std::vector<Real> _resid;
  std::vector<Real> _jac;
  std::vector<Real> _dvar;
 
  SymmElasticityTensor _elasticity;
  SymmTensor _stress_old, _dstrain, _stress_new;
  SymmTensor _stress0, _dstress0;
  bool _nconv;   
  bool _err_tol; 
 
private:
};