ComputeSmearedCrackingStress.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 "ComputeMultipleInelasticStress.h"
#include "SmearedCrackSofteningBase.h"
#include "Function.h"

class ComputeSmearedCrackingStress : public ComputeMultipleInelasticStress
{
public:
  static InputParameters validParams();

  ComputeSmearedCrackingStress(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(RealVectorValue & strain_in_crack_dir);

  void updateStressTensorForCracking(RankTwoTensor & tensor, const RealVectorValue & sigma);

  bool previouslyCracked();


  const VariableValue & _cracking_stress;

  std::vector<unsigned int> _prescribed_crack_directions;

  const unsigned int _max_cracks;

  const Real _cracking_neg_fraction;

  const Real _shear_retention_factor;

  const Real _max_stress_correction;

  const enum class CrackedElasticityType { DIAGONAL, FULL } _cracked_elasticity_type;

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

  MaterialProperty<RealVectorValue> & _crack_flags;

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

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

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

  RankFourTensor _local_elasticity_tensor;

  std::vector<SmearedCrackSofteningBase *> _softening_models;

  std::vector<Real> _local_elastic_vector;
};