#include <XFEMCrackGrowthIncrement2DCut.h>

Public Member Functions
	XFEMCrackGrowthIncrement2DCut (Real x0, Real y0, Real x1, Real y1, Real t0, Real t1)

virtual bool	cutElementByCrackGrowthIncrement (const Elem *elem, std::vector< CutEdgeForCrackGrowthIncr > &cut_edges, Real time)

Real	cutCompletionFraction (Real time)

Protected Attributes
const std::pair< Real, Real >	_time_range

Private Attributes
const std::pair< Point, Point >	_cut_line_endpoints

Detailed Description

Definition at line 26 of file XFEMCrackGrowthIncrement2DCut.h.

Constructor & Destructor Documentation

◆ XFEMCrackGrowthIncrement2DCut()

XFEMCrackGrowthIncrement2DCut::XFEMCrackGrowthIncrement2DCut	(	Real	x0,
		Real	y0,
		Real	x1,
		Real	y1,
		Real	t0,
		Real	t1
	)

Definition at line 16 of file XFEMCrackGrowthIncrement2DCut.C.

   : _time_range(std::make_pair(t0, t1)),
     _cut_line_endpoints(std::make_pair(Point(x0, y0, 0.0), Point(x1, y1, 0.0)))
 {
 }

Member Function Documentation

◆ cutCompletionFraction()

Real XFEMCrackGrowthIncrement2DCut::cutCompletionFraction ( Real time )

Definition at line 24 of file XFEMCrackGrowthIncrement2DCut.C.

Referenced by cutElementByCrackGrowthIncrement().

 {
   Real fraction = 0.0;
   if (time >= _time_range.first)
   {
     if (time >= _time_range.second)
       fraction = 1.0;
     else
       fraction = (time - _time_range.first) / (_time_range.second - _time_range.first);
   }
   return fraction;
 }

◆ cutElementByCrackGrowthIncrement()

bool XFEMCrackGrowthIncrement2DCut::cutElementByCrackGrowthIncrement	(	const Elem *	elem,
		std::vector< CutEdgeForCrackGrowthIncr > &	cut_edges,
		Real	time
	)

virtual

Definition at line 38 of file XFEMCrackGrowthIncrement2DCut.C.

Referenced by XFEM::markCutEdgesByState().

 {
   bool cut_elem = false;
 
   Real fraction = cutCompletionFraction(time);
 
   if (fraction > 0.0)
   {
     unsigned int n_sides = elem->n_sides();
 
     for (unsigned int i = 0; i < n_sides; ++i)
     {
       // This returns the lowest-order type of side, which should always
       // be an EDGE2 here because this class is for 2D only.
       std::unique_ptr<const Elem> curr_side = elem->side_ptr(i);
       if (curr_side->type() != EDGE2)
         mooseError("In cutElementByGeometry element side must be EDGE2, but type is: ",
                    libMesh::Utility::enum_to_string(curr_side->type()),
                    " base element type is: ",
                    libMesh::Utility::enum_to_string(elem->type()));
 
       const Node * node1 = curr_side->node_ptr(0);
       const Node * node2 = curr_side->node_ptr(1);
       Real seg_int_frac = 0.0;
 
       if (Xfem::intersectSegmentWithCutLine(
               *node1, *node2, _cut_line_endpoints, fraction, seg_int_frac))
       {
         cut_elem = true;
         CutEdgeForCrackGrowthIncr mycut;
         mycut._id1 = node1->id();
         mycut._id2 = node2->id();
         mycut._distance = seg_int_frac;
         mycut._host_side_id = i;
         cut_edges.push_back(mycut);
       }
     }
   }
   return cut_elem;
 }