CrackTipEnrichmentStressDivergenceTensors.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 "ALEKernel.h"
#include "RankTwoTensor.h"
#include "RankFourTensor.h"
#include "CrackFrontDefinition.h"
#include "EnrichmentFunctionCalculation.h"

// Forward Declarations

class CrackTipEnrichmentStressDivergenceTensors : public ALEKernel,
                                                  public EnrichmentFunctionCalculation
{
public:
  static InputParameters validParams();

  CrackTipEnrichmentStressDivergenceTensors(const InputParameters & parameters);

protected:
  virtual Real computeQpResidual() override;
  virtual Real computeQpJacobian() override;
  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;

  const std::string _base_name;

  const MaterialProperty<RankTwoTensor> & _stress;
  const MaterialProperty<RankFourTensor> & _Jacobian_mult;

  const MaterialProperty<RankTwoTensor> * _deformation_gradient;
  const MaterialProperty<RankTwoTensor> * _deformation_gradient_old;
  const MaterialProperty<RankTwoTensor> * _rotation_increment;

  const unsigned int _component;
  const unsigned int _enrichment_component;

  unsigned int _nenrich_disp;
  std::vector<unsigned int> _enrich_disp_var;

  unsigned int _ndisp;
  std::vector<unsigned int> _disp_var;

private:
  std::vector<Real> _B;
  std::vector<RealVectorValue> _dBX;
  std::vector<RealVectorValue> _dBx;
  std::vector<Real> _BI;
  std::vector<Real> _BJ;
};