BiLinearMixedModeTraction.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 "CZMComputeLocalTractionTotalBase.h"
 
class BiLinearMixedModeTraction : public CZMComputeLocalTractionTotalBase
{
public:
  static InputParameters validParams();
  BiLinearMixedModeTraction(const InputParameters & parameters);
 
protected:
  virtual void initQpStatefulProperties() override;
 
  virtual void computeInterfaceTractionAndDerivatives() override;
 
  virtual RealVectorValue computeTraction();
 
  virtual RankTwoTensor computeTractionDerivatives();
 
  // The damage evolution law
  virtual void computeDamage();
 
  const Real _K;
 
  MaterialProperty<Real> & _d;
  const MaterialProperty<Real> & _d_old;
 
  const MaterialProperty<RealVectorValue> & _interface_displacement_jump_old;
 
  // @{
  // The parameters in the damage evolution law
  MaterialProperty<Real> & _delta_init;
  MaterialProperty<Real> & _delta_final;
  MaterialProperty<Real> & _delta_m;
  // @}
 
  const MaterialProperty<Real> & _GI_c;
 
  const MaterialProperty<Real> & _GII_c;
 
  const MaterialProperty<Real> & _N;
 
  const MaterialProperty<Real> & _S;
 
  const Real _eta;
 
  MaterialProperty<Real> & _beta;
 
  const Real _viscosity;
 
  enum class MixedModeCriterion
  {
    POWER_LAW,
    BK
  } _criterion;
 
private:
  // Compute mode mixity ratio
  void computeModeMixity();
 
  // Compute relative displacement jump at damage initiation
  void computeCriticalDisplacementJump();
 
  // Compute relative displacement jump at full degradation
  void computeFinalDisplacementJump();
 
  // Compute mixed mode relative displacement
  void computeEffectiveDisplacementJump();
 
  // @{
  // Parameters to improve numerical convergence
  const bool _lag_mode_mixity;
  const bool _lag_disp_jump;
  const Real _alpha;
  // @}
 
  // @{
  // Intermediate derivatives used in the chain rule for computing the derivative of damage w.r.t.
  // displacement jumps
  RealVectorValue _dbeta_ddelta;
  RealVectorValue _ddelta_init_ddelta;
  RealVectorValue _ddelta_final_ddelta;
  RealVectorValue _ddelta_m_ddelta;
  RealVectorValue _dd_ddelta;
  // @}
};