LinearViscoelasticityBase.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 "ComputeElasticityTensorBase.h"
#include "RankFourTensor.h"
#include "RankTwoTensor.h"
 
class LinearViscoelasticityBase;
 
template <>
InputParameters validParams<LinearViscoelasticityBase>();
 
class LinearViscoelasticityBase : public ComputeElasticityTensorBase
{
public:
  enum class IntegrationRule
  {
    BackwardEuler,
    MidPoint,
    Newmark,
    Zienkiewicz,
  };
 
  static InputParameters validParams();
 
  LinearViscoelasticityBase(const InputParameters & parameters);
 
  void recomputeQpApparentProperties(unsigned int qp);
 
protected:
  virtual void initQpStatefulProperties() override;
  virtual void computeQpElasticityTensor() final;
 
  virtual void computeQpViscoelasticProperties() = 0;
 
  virtual void computeQpViscoelasticPropertiesInv();
 
  virtual void computeQpApparentElasticityTensors() = 0;
 
  virtual void computeQpApparentCreepStrain() = 0;
 
  virtual void updateQpViscousStrains() = 0;
 
  void declareViscoelasticProperties();
 
  Real computeTheta(Real dt, Real viscosity) const;
 
  IntegrationRule _integration_rule;
  Real _theta;
 
  MaterialProperty<RankFourTensor> & _apparent_elasticity_tensor;
  MaterialProperty<RankFourTensor> & _apparent_elasticity_tensor_inv;
 
  //  MaterialProperty<RankFourTensor> & _instantaneous_elasticity_tensor;
  MaterialProperty<RankFourTensor> & _elasticity_tensor_inv;
 
  bool _need_viscoelastic_properties_inverse;
  bool _has_longterm_dashpot;
  unsigned int _components;
 
  MaterialProperty<RankFourTensor> & _first_elasticity_tensor;
  MaterialProperty<RankFourTensor> * _first_elasticity_tensor_inv;
 
  std::vector<MaterialProperty<RankFourTensor> *> _springs_elasticity_tensors;
  std::vector<MaterialProperty<RankFourTensor> *> _springs_elasticity_tensors_inv;
  std::vector<const MaterialProperty<RankFourTensor> *> _springs_elasticity_tensors_inv_old;
 
  std::vector<MaterialProperty<Real> *> _dashpot_viscosities;
  std::vector<const MaterialProperty<Real> *> _dashpot_viscosities_old;
 
  std::vector<MaterialProperty<RankTwoTensor> *> _viscous_strains;
  std::vector<const MaterialProperty<RankTwoTensor> *> _viscous_strains_old;
 
  MaterialProperty<RankTwoTensor> & _apparent_creep_strain;
  const MaterialProperty<RankTwoTensor> & _apparent_creep_strain_old;
 
  const MaterialProperty<RankTwoTensor> & _elastic_strain_old;
  const MaterialProperty<RankTwoTensor> & _creep_strain_old;
 
  bool _has_driving_eigenstrain;
  std::string _driving_eigenstrain_name;
  const MaterialProperty<RankTwoTensor> * _driving_eigenstrain;
  const MaterialProperty<RankTwoTensor> * _driving_eigenstrain_old;
 
  bool _force_recompute_properties;
 
  bool & _step_zero;
};