XFEMAction Class Reference

#include <XFEMAction.h>

Inheritance diagram for XFEMAction:

Public Member Functions
	XFEMAction (InputParameters params)
virtual void	act ()

Protected Attributes
std::vector< UserObjectName >	_geom_cut_userobjects
std::string	_xfem_qrule
std::string	_order
std::string	_family
bool	_xfem_cut_plane
bool	_xfem_use_crack_growth_increment
Real	_xfem_crack_growth_increment
bool	_use_crack_tip_enrichment
UserObjectName	_crack_front_definition
std::vector< VariableName >	_enrich_displacements
std::vector< VariableName >	_displacements
std::vector< BoundaryName >	_cut_off_bc
Real	_cut_off_radius

Detailed Description

Definition at line 20 of file XFEMAction.h.

Constructor & Destructor Documentation

XFEMAction()

XFEMAction::XFEMAction

(

InputParameters

params

)

Definition at line 81 of file XFEMAction.C.

  : Action(params),
    _geom_cut_userobjects(getParam<std::vector<UserObjectName>>("geometric_cut_userobjects")),
    _xfem_qrule(getParam<std::string>("qrule")),
    _xfem_cut_plane(false),
    _xfem_use_crack_growth_increment(getParam<bool>("use_crack_growth_increment")),
    _xfem_crack_growth_increment(getParam<Real>("crack_growth_increment")),
    _use_crack_tip_enrichment(getParam<bool>("use_crack_tip_enrichment"))
{
  _order = "CONSTANT";
  _family = "MONOMIAL";
  if (isParamValid("output_cut_plane"))
    _xfem_cut_plane = getParam<bool>("output_cut_plane");
 
  if (_use_crack_tip_enrichment)
  {
    if (isParamValid("crack_front_definition"))
      _crack_front_definition = getParam<UserObjectName>("crack_front_definition");
    else
      mooseError("To add crack tip enrichment, crack_front_definition must be provided.");
 
    if (isParamValid("displacements"))
      _displacements = getParam<std::vector<VariableName>>("displacements");
    else
      mooseError("To add crack tip enrichment, displacements must be provided.");
 
    if (isParamValid("enrichment_displacements"))
    {
      _enrich_displacements = getParam<std::vector<VariableName>>("enrichment_displacements");
      if (_enrich_displacements.size() != 8 && _displacements.size() == 2)
        mooseError("The number of enrichment displacements should be total 8 for 2D.");
      else if (_enrich_displacements.size() != 12 && _displacements.size() == 3)
        mooseError("The number of enrichment displacements should be total 12 for 3D.");
    }
    else
      mooseError("To add crack tip enrichment, enrichment_displacements must be provided.");
 
    if (isParamValid("cut_off_boundary"))
      _cut_off_bc = getParam<std::vector<BoundaryName>>("cut_off_boundary");
    else
      mooseError("To add crack tip enrichment, cut_off_boundary must be provided.");
 
    if (isParamValid("cut_off_radius"))
      _cut_off_radius = getParam<Real>("cut_off_radius");
    else
      mooseError("To add crack tip enrichment, cut_off_radius must be provided.");
  }
}

Member Function Documentation

act()

void XFEMAction::act ( )

virtual

Definition at line 131 of file XFEMAction.C.

{
 
  std::shared_ptr<XFEMInterface> xfem_interface = _problem->getXFEM();
  if (xfem_interface == NULL)
  {
    _pars.set<FEProblemBase *>("_fe_problem_base") = &*_problem;
    std::shared_ptr<XFEM> new_xfem(new XFEM(_pars));
    _problem->initXFEM(new_xfem);
    xfem_interface = _problem->getXFEM();
  }
 
  std::shared_ptr<XFEM> xfem = MooseSharedNamespace::dynamic_pointer_cast<XFEM>(xfem_interface);
  if (xfem == NULL)
    mooseError("dynamic cast of xfem object failed");
 
  if (_current_task == "setup_xfem")
  {
    xfem->setXFEMQRule(_xfem_qrule);
 
    xfem->setCrackGrowthMethod(_xfem_use_crack_growth_increment, _xfem_crack_growth_increment);
    xfem->setDebugOutputLevel(getParam<unsigned int>("debug_output_level"));
  }
  else if (_current_task == "add_variable" && _use_crack_tip_enrichment)
  {
    auto var_params = _factory.getValidParams("MooseVariable");
    var_params.set<MooseEnum>("family") = "LAGRANGE";
    var_params.set<MooseEnum>("order") = "FIRST";
 
    for (const auto & enrich_disp : _enrich_displacements)
      _problem->addVariable("MooseVariable", enrich_disp, var_params);
  }
  else if (_current_task == "add_kernel" && _use_crack_tip_enrichment)
  {
    for (unsigned int i = 0; i < _enrich_displacements.size(); ++i)
    {
      InputParameters params = _factory.getValidParams("CrackTipEnrichmentStressDivergenceTensors");
      params.set<NonlinearVariableName>("variable") = _enrich_displacements[i];
      params.set<unsigned int>("component") = i / 4;
      params.set<unsigned int>("enrichment_component") = i % 4;
      params.set<UserObjectName>("crack_front_definition") = _crack_front_definition;
      params.set<std::vector<VariableName>>("enrichment_displacements") = _enrich_displacements;
      params.set<std::vector<VariableName>>("displacements") = _displacements;
      _problem->addKernel(
          "CrackTipEnrichmentStressDivergenceTensors", _enrich_displacements[i], params);
    }
  }
  else if (_current_task == "add_bc" && _use_crack_tip_enrichment)
  {
    for (unsigned int i = 0; i < _enrich_displacements.size(); ++i)
    {
      InputParameters params = _factory.getValidParams("CrackTipEnrichmentCutOffBC");
      params.set<NonlinearVariableName>("variable") = _enrich_displacements[i];
      params.set<Real>("value") = 0;
      params.set<std::vector<BoundaryName>>("boundary") = _cut_off_bc;
      params.set<Real>("cut_off_radius") = _cut_off_radius;
      params.set<UserObjectName>("crack_front_definition") = _crack_front_definition;
      _problem->addBoundaryCondition(
          "CrackTipEnrichmentCutOffBC", _enrich_displacements[i], params);
    }
  }
  else if (_current_task == "add_aux_variable" && _xfem_cut_plane)
  {
    auto var_params = _factory.getValidParams("MooseVariableConstMonomial");
 
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut_origin_x", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut_origin_y", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut_origin_z", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut_normal_x", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut_normal_y", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut_normal_z", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut2_origin_x", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut2_origin_y", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut2_origin_z", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut2_normal_x", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut2_normal_y", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_cut2_normal_z", var_params);
    _problem->addAuxVariable("MooseVariableConstMonomial", "xfem_volfrac", var_params);
  }
  else if (_current_task == "add_aux_kernel" && _xfem_cut_plane)
  {
    InputParameters params = _factory.getValidParams("XFEMVolFracAux");
    params.set<ExecFlagEnum>("execute_on") = EXEC_TIMESTEP_BEGIN;
    params.set<AuxVariableName>("variable") = "xfem_volfrac";
    _problem->addAuxKernel("XFEMVolFracAux", "xfem_volfrac", params);
 
    params = _factory.getValidParams("XFEMCutPlaneAux");
    params.set<ExecFlagEnum>("execute_on") = EXEC_TIMESTEP_END;
 
    // first cut plane
    params.set<unsigned int>("plane_id") = 0;
 
    params.set<AuxVariableName>("variable") = "xfem_cut_origin_x";
    params.set<MooseEnum>("quantity") = "origin_x";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut_origin_x", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut_origin_y";
    params.set<MooseEnum>("quantity") = "origin_y";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut_origin_y", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut_origin_z";
    params.set<MooseEnum>("quantity") = "origin_z";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut_origin_z", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut_normal_x";
    params.set<MooseEnum>("quantity") = "normal_x";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut_normal_x", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut_normal_y";
    params.set<MooseEnum>("quantity") = "normal_y";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut_normal_y", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut_normal_z";
    params.set<MooseEnum>("quantity") = "normal_z";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut_normal_z", params);
 
    // second cut plane
    params.set<unsigned int>("plane_id") = 1;
 
    params.set<AuxVariableName>("variable") = "xfem_cut2_origin_x";
    params.set<MooseEnum>("quantity") = "origin_x";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut2_origin_x", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut2_origin_y";
    params.set<MooseEnum>("quantity") = "origin_y";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut2_origin_y", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut2_origin_z";
    params.set<MooseEnum>("quantity") = "origin_z";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut2_origin_z", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut2_normal_x";
    params.set<MooseEnum>("quantity") = "normal_x";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut2_normal_x", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut2_normal_y";
    params.set<MooseEnum>("quantity") = "normal_y";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut2_normal_y", params);
 
    params.set<AuxVariableName>("variable") = "xfem_cut2_normal_z";
    params.set<MooseEnum>("quantity") = "normal_z";
    _problem->addAuxKernel("XFEMCutPlaneAux", "xfem_cut2_normal_z", params);
  }
}

Member Data Documentation

_crack_front_definition

UserObjectName XFEMAction::_crack_front_definition

protected

Definition at line 36 of file XFEMAction.h.

Referenced by act(), and XFEMAction().

_cut_off_bc

std::vector<BoundaryName> XFEMAction::_cut_off_bc

protected

Definition at line 39 of file XFEMAction.h.

Referenced by act(), and XFEMAction().

_cut_off_radius

Real XFEMAction::_cut_off_radius

protected

Definition at line 40 of file XFEMAction.h.

Referenced by act(), and XFEMAction().

_displacements

std::vector<VariableName> XFEMAction::_displacements

protected

Definition at line 38 of file XFEMAction.h.

Referenced by act(), and XFEMAction().

_enrich_displacements

std::vector<VariableName> XFEMAction::_enrich_displacements

protected

Definition at line 37 of file XFEMAction.h.

Referenced by act(), and XFEMAction().

_family

std::string XFEMAction::_family

protected

Definition at line 31 of file XFEMAction.h.

Referenced by XFEMAction().

_geom_cut_userobjects

std::vector<UserObjectName> XFEMAction::_geom_cut_userobjects

protected

Definition at line 28 of file XFEMAction.h.

_order

std::string XFEMAction::_order

protected

Definition at line 30 of file XFEMAction.h.

Referenced by XFEMAction().

_use_crack_tip_enrichment

bool XFEMAction::_use_crack_tip_enrichment

protected

Definition at line 35 of file XFEMAction.h.

Referenced by act(), and XFEMAction().

_xfem_crack_growth_increment

Real XFEMAction::_xfem_crack_growth_increment

protected

Definition at line 34 of file XFEMAction.h.

Referenced by act().

_xfem_cut_plane

bool XFEMAction::_xfem_cut_plane

protected

Definition at line 32 of file XFEMAction.h.

Referenced by act(), and XFEMAction().

_xfem_qrule

std::string XFEMAction::_xfem_qrule

protected

Definition at line 29 of file XFEMAction.h.

Referenced by act().

_xfem_use_crack_growth_increment

bool XFEMAction::_xfem_use_crack_growth_increment

protected

Definition at line 33 of file XFEMAction.h.

Referenced by act().

---

The documentation for this class was generated from the following files:

• xfem/include/actions/XFEMAction.h
• xfem/src/actions/XFEMAction.C