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;
  // @}
};