BilinearMixedModeCohesiveZoneModel.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 "CohesiveZoneModelBase.h"
#include "TwoVector.h"
 
class BilinearMixedModeCohesiveZoneModel : virtual public CohesiveZoneModelBase
{
public:
  static InputParameters validParams();
 
  BilinearMixedModeCohesiveZoneModel(const InputParameters & parameters);
 
  virtual void initialize() override;
  virtual void finalize() override;
 
  // Getters for analysis output
  Real getModeMixityRatio(const Node * const node) const;
  Real getCohesiveDamage(const Node * const node) const;
  Real getLocalDisplacementNormal(const Node * const node) const;
  Real getLocalDisplacementTangential(const Node * const node) const;
  Real getLocalDisplacementTangentialOne(const Node * const node) const;
  Real getLocalDisplacementTangentialTwo(const Node * const node) const;
  Real getLocalDisplacementTangentialEffective(const Node * const node) const;
  Real getCohesiveTractionNormal(const Node * const node) const;
  Real getCohesiveTractionTangentialMagnitude(const Node * const node) const;
  Real getCohesiveTractionTangentialOne(const Node * const node) const;
  Real getCohesiveTractionTangentialTwo(const Node * const node) const;
  Real getCohesiveTractionEffective(const Node * const node) const;
 
protected:
  virtual void computeQpProperties() override;
  virtual void computeQpIProperties() override;
 
  virtual bool constrainedByOwner() const override { return false; }
 
  // @{
  // Compute CZM bilinear traction law.
  virtual void computeModeMixity(const Node * const node);
  virtual void computeCriticalDisplacementJump(const Node * const node);
  virtual void computeFinalDisplacementJump(const Node * const node);
  virtual void computeEffectiveDisplacementJump(const Node * const node);
  virtual void computeDamage(const Node * const node) override;
  // @}
 
  virtual void computeCZMTraction(const Node * const node) override;
 
  const bool _set_compressive_traction_to_zero;
 
  std::unordered_map<const DofObject *, ADReal> _dof_to_mode_mixity_ratio;
 
  std::unordered_map<const DofObject *, RealVectorValue> _dof_to_local_czm_traction;
 
  const MaterialProperty<Real> & _normal_strength;
 
  const MaterialProperty<Real> & _shear_strength;
 
  const MaterialProperty<Real> & _GI_c;
 
  const MaterialProperty<Real> & _GII_c;
 
  ADReal _normal_strength_interpolation;
 
  ADReal _shear_strength_interpolation;
 
  ADReal _GI_c_interpolation;
 
  ADReal _GII_c_interpolation;
 
  // Penalty stiffness for bilinear traction model
  const Real _penalty_stiffness_czm;
 
  enum class MixedModeCriterion
  {
    POWER_LAW,
    BK
  } _mix_mode_criterion;
 
  const Real _power_law_parameter;
 
  const Real _viscosity;
 
  const Real _regularization_alpha;
 
  // @{
  // Strength material properties at the nodes
  std::unordered_map<const DofObject *, ADReal> _dof_to_normal_strength;
  std::unordered_map<const DofObject *, ADReal> _dof_to_shear_strength;
  // @}
 
  // @{
  // Fracture properties at the nodes
  std::unordered_map<const DofObject *, ADReal> _dof_to_GI_c;
  std::unordered_map<const DofObject *, ADReal> _dof_to_GII_c;
  // @}
 
  // @{
  // The parameters in the damage evolution law: Maps node to damage
  std::unordered_map<const DofObject *, ADReal> _dof_to_delta_initial;
  std::unordered_map<const DofObject *, ADReal> _dof_to_delta_final;
  std::unordered_map<const DofObject *, ADReal> _dof_to_delta_max;
  // @}
};