FractureMechanics::AnalyticalForces Struct Reference

#include <users_modules/fracture_mechanics/AnalyticalFun.hpp>

Inheritance diagram for FractureMechanics::AnalyticalForces:

Collaboration diagram for FractureMechanics::AnalyticalForces:

Public Member Functions
PetscErrorCode	getForce (const EntityHandle ent, const VectorDouble3 &coords, const VectorDouble3 &normal, VectorDouble3 &forces)
Public Member Functions inherited from NeumannForcesSurface::MethodForAnalyticalForce
virtual	~MethodForAnalyticalForce ()=default
Public Member Functions inherited from FractureMechanics::AnalyticalOptions
PetscErrorCode	getOptions ()
	AnalyticalOptions ()

Additional Inherited Members
Public Attributes inherited from FractureMechanics::AnalyticalOptions
double	aLpha
bool	optionsInitialised

Detailed Description

Definition at line 51 of file AnalyticalFun.hpp.

Member Function Documentation

getForce()

PetscErrorCode FractureMechanics::AnalyticalForces::getForce ( const EntityHandle  ent, const VectorDouble3 &  coords, const VectorDouble3 &  normal, VectorDouble3 &  force  )

inlinevirtual

User implemented analytical force

Parameters

coords	coordinates of integration point
normal	normal at integration point
force	returned force

Returns

error code

Reimplemented from NeumannForcesSurface::MethodForAnalyticalForce.

Definition at line 54 of file AnalyticalFun.hpp.

                                                                              {
 
    MoFEMFunctionBeginHot;
 
    // initailise options
    if (!optionsInitialised) {
      ierr = getOptions();
      CHKERRQ(ierr);
    }
 
    const double x = coords[0];
    const double y = coords[1];
    const double z = coords[2];
    double &fx = forces[0];
    double &fy = forces[1];
    double &fz = forces[2];
 
    const double n0 = normal[0] / sqrt(pow(normal[0], 2) + pow(normal[1], 2) +
                                       pow(normal[2], 2));
    const double n1 = normal[1] / sqrt(pow(normal[0], 2) + pow(normal[1], 2) +
                                       pow(normal[2], 2));
    const double n2 = normal[2] / sqrt(pow(normal[0], 2) + pow(normal[1], 2) +
                                       pow(normal[2], 2));
 
    const double pi = M_PI;
    const double alpha = aLpha;
    const double theta = atan2(z, x * sin(alpha) + y * cos(alpha));
 
    const double rxy =
        sqrt(pow(x * sin(alpha) + y * cos(alpha), 2) + pow(z, 2));
    const double sref1 = 1 / sqrt(2 * pi * rxy) * cos(theta / 2.) *
                         (1. - sin(theta / 2.) * sin(3. * theta / 2.));
    const double sref2 = 1 / sqrt(2 * pi * rxy) * cos(theta / 2.) *
                         (1. + sin(theta / 2.) * sin(3. * theta / 2.));
    const double sref3 = 1 / sqrt(2 * pi * rxy) * sin(theta / 2.) *
                         cos(theta / 2.) * cos(3. * theta / 2.);
 
    const double scal1 = sref1 * pow(sin(alpha), 2);
    const double scal2 = sref1 * pow(cos(alpha), 2);
    const double scal3 = sref2;
    const double scal4 = sref1 * sin(alpha) * cos(alpha);
    const double scal5 = sref3 * sin(alpha);
    const double scal6 = sref3 * cos(alpha);
 
    // This is just example
    fx = scal1 * n0 + scal4 * n1 + scal5 * n2;
    fy = scal4 * n0 + scal2 * n1 + scal6 * n2;
    fz = scal5 * n0 + scal6 * n1 + scal3 * n2;
 
    fx *= -1;
    fy *= -1;
    fz *= -1;
 
    MoFEMFunctionReturnHot(0);
  }

---

The documentation for this struct was generated from the following file:

• users_modules/fracture_mechanics/AnalyticalFun.hpp