ADComputeSmearedCrackingStress.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 "ColumnMajorMatrix.h"
#include "ADComputeMultipleInelasticStress.h"
#include "ADSmearedCrackSofteningBase.h"
#include "Function.h"

class ADComputeSmearedCrackingStress : public ADComputeMultipleInelasticStress
{
public:
  static InputParameters validParams();

  ADComputeSmearedCrackingStress(const InputParameters & parameters);

  virtual void initialSetup() override;
  virtual void initQpStatefulProperties() override;
  virtual void computeQpStress() override;

protected:
  void updateLocalElasticityTensor();

  virtual void updateCrackingStateAndStress();

  virtual unsigned int getNumKnownCrackDirs() const;

  void computeCrackStrainAndOrientation(ADRealVectorValue & strain_in_crack_dir);

  void updateStressTensorForCracking(ADRankTwoTensor & tensor, const ADRealVectorValue & sigma);

  bool previouslyCracked();


  const ADVariableValue & _cracking_stress;

  std::vector<unsigned int> _prescribed_crack_directions;

  const unsigned int _max_cracks;

  const ADReal _cracking_neg_fraction;

  const ADReal _shear_retention_factor;

  const ADReal _max_stress_correction;

  const enum class CrackedElasticityType { DIAGONAL, FULL } _cracked_elasticity_type;

  ADMaterialProperty<RealVectorValue> & _crack_damage;
  const MaterialProperty<RealVectorValue> & _crack_damage_old;

  ADMaterialProperty<RealVectorValue> & _crack_flags;

  ADMaterialProperty<RankTwoTensor> & _crack_rotation;
  const MaterialProperty<RankTwoTensor> & _crack_rotation_old;

  ADMaterialProperty<RealVectorValue> & _crack_initiation_strain;
  const MaterialProperty<RealVectorValue> & _crack_initiation_strain_old;

  ADMaterialProperty<RealVectorValue> & _crack_max_strain;
  const MaterialProperty<RealVectorValue> & _crack_max_strain_old;

  ADRankFourTensor _local_elasticity_tensor;

  std::vector<ADSmearedCrackSofteningBase *> _softening_models;

  std::vector<ADReal> _local_elastic_vector;
};