Evaluate stress at integration points. More...

#include <users_modules/fracture_mechanics/src/MWLS.hpp>

Inheritance diagram for FractureMechanics::MWLSApprox::OpMWLSStressAtGaussPts:

Collaboration diagram for FractureMechanics::MWLSApprox::OpMWLSStressAtGaussPts:

Public Member Functions
	OpMWLSStressAtGaussPts (boost::shared_ptr< MatrixDouble > mat_pos_at_pts_ptr, boost::shared_ptr< MatrixDouble > mat_grad_pos_at_pts_ptr, boost::shared_ptr< CrackFrontElement > fe_singular_ptr, boost::shared_ptr< MWLSApprox > mwls_approx, const std::string stress_tag_name, bool calculate_derivative, bool testing=false)

virtual	~OpMWLSStressAtGaussPts ()=default

Public Member Functions inherited from FractureMechanics::MWLSApprox::OpMWLSBase< VolumeElementForcesAndSourcesCore >
	OpMWLSBase (boost::shared_ptr< MatrixDouble > mat_pos_at_pts_ptr, boost::shared_ptr< MatrixDouble > mat_grad_pos_at_pts_ptr, boost::shared_ptr< CrackFrontElement > fe_singular_ptr, boost::shared_ptr< MWLSApprox > mwls_approx, bool testing=false)

virtual	~OpMWLSBase ()=default

MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
std::string	stressTagName

const bool	calculateDerivative

const bool	testing

boost::shared_ptr< MatrixDouble >	matGradPosAtPtsPtr

Public Attributes inherited from FractureMechanics::MWLSApprox::OpMWLSBase< VolumeElementForcesAndSourcesCore >
boost::shared_ptr< MatrixDouble >	matPosAtPtsPtr

boost::shared_ptr< MatrixDouble >	matGradPosAtPtsPtr

boost::weak_ptr< CrackFrontElement >	feSingularPtr

boost::shared_ptr< MWLSApprox >	mwlsApprox

Protected Member Functions
MoFEMErrorCode	doMWLSWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)
	Do specific work in operator. More...

Detailed Description

Evaluate stress at integration points.

Definition at line 404 of file MWLS.hpp.

Constructor & Destructor Documentation

◆ OpMWLSStressAtGaussPts()

FractureMechanics::MWLSApprox::OpMWLSStressAtGaussPts::OpMWLSStressAtGaussPts	(	boost::shared_ptr< MatrixDouble >	mat_pos_at_pts_ptr,
		boost::shared_ptr< MatrixDouble >	mat_grad_pos_at_pts_ptr,
		boost::shared_ptr< CrackFrontElement >	fe_singular_ptr,
		boost::shared_ptr< MWLSApprox >	mwls_approx,
		const std::string	stress_tag_name,
		bool	calculate_derivative,
		bool	testing = `false`
	)

inline

Definition at line 410 of file MWLS.hpp.

         : OpMWLSBase<VolumeElementForcesAndSourcesCore>(
               mat_pos_at_pts_ptr, mat_grad_pos_at_pts_ptr, fe_singular_ptr,
               mwls_approx, testing),
           matGradPosAtPtsPtr(mat_pos_at_pts_ptr),
           stressTagName(stress_tag_name),
           calculateDerivative(calculate_derivative), testing(testing) {}

◆ ~OpMWLSStressAtGaussPts()

virtual FractureMechanics::MWLSApprox::OpMWLSStressAtGaussPts::~OpMWLSStressAtGaussPts ( )

virtualdefault

Member Function Documentation

◆ doMWLSWork()

MoFEMErrorCode FractureMechanics::MWLSApprox::OpMWLSStressAtGaussPts::doMWLSWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

protectedvirtual

Do specific work in operator.

Is protected to prevent direct call. This can be only called by doWork method.

Parameters

side
type
data

Returns: MoFEMErrorCode

Implements FractureMechanics::MWLSApprox::OpMWLSBase< VolumeElementForcesAndSourcesCore >.

Definition at line 1582 of file MWLS.cpp.

                                                                       {
   MoFEMFunctionBegin;
  
   if (testing)
     cerr << "Element " << getFEEntityHandle() << endl;
  
   FTensor::Index<'i', 3> i;
  
   const int nb_gauss_pts = data.getN().size1();
   const bool use_global_base =
       mwlsApprox->getUseGlobalBaseAtMaterialReferenceConfiguration();
   Tag th;
   CHKERR mwlsApprox->mwlsMoab.tag_get_handle(stressTagName.c_str(), th);
  
   MatrixDouble &stress = mwlsApprox->stressAtGaussPts;
   stress.resize(6, nb_gauss_pts, false);
   MatrixDouble &diff_stress = mwlsApprox->diffStressAtGaussPts;
   diff_stress.resize(18, nb_gauss_pts, false);
  
   auto t_mwls_material_positions =
       getFTensor1FromMat<3>(mwlsApprox->mwlsMaterialPositions);
   auto t_approx_base_point = getFTensor1FromMat<3>(mwlsApprox->approxBasePoint);
   for (int gg = 0; gg != nb_gauss_pts; ++gg) {
  
     auto t_material_positions = getFTensor1FromMat<3>(*matPosAtPtsPtr);
  
     FTensor::Tensor1<double, 3> t_mat_pos;
     t_mat_pos(i) = t_mwls_material_positions(i);
     FTensor::Tensor1<double, 3> t_approx_base_pos;
     t_approx_base_pos(i) = t_approx_base_point(i);
  
     if (testing) {
       cerr << "material_positions " << t_mwls_material_positions << endl;
       cerr << "approx_base_point " << t_approx_base_point << endl;
     }
     CHKERR mwlsApprox->getInfluenceNodes(&t_approx_base_pos(0));
     CHKERR mwlsApprox->calculateApproxCoeff(true, use_global_base);
     CHKERR mwlsApprox->evaluateApproxFun(&t_mat_pos(0));
     CHKERR mwlsApprox->evaluateFirstDiffApproxFun(&t_mat_pos(0),
                                                   use_global_base);
  
     CHKERR mwlsApprox->getTagData(th);
     const auto &vals = mwlsApprox->getDataApprox();
  
     for (int ii = 0; ii != 6; ++ii)
       stress(ii, gg) = vals[ii];
     if (calculateDerivative) {
       const auto &diff_vals = mwlsApprox->getDiffDataApprox();
       for (int ii = 0; ii != 6; ++ii) {
         for (int jj = 0; jj != 3; ++jj)
           diff_stress(jj * 6 + ii, gg) = diff_vals(jj, ii);
       }
     }
     ++t_mwls_material_positions;
     ++t_approx_base_point;
   }
  
   if (testing) {
     cerr << "stress " << stress << endl;
     cerr << endl;
   }
  
   MoFEMFunctionReturn(0);
 }

Member Data Documentation

◆ calculateDerivative

const bool FractureMechanics::MWLSApprox::OpMWLSStressAtGaussPts::calculateDerivative

Definition at line 407 of file MWLS.hpp.

◆ matGradPosAtPtsPtr

boost::shared_ptr<MatrixDouble> FractureMechanics::MWLSApprox::OpMWLSStressAtGaussPts::matGradPosAtPtsPtr

Definition at line 409 of file MWLS.hpp.

◆ stressTagName

std::string FractureMechanics::MWLSApprox::OpMWLSStressAtGaussPts::stressTagName

Definition at line 406 of file MWLS.hpp.

◆ testing

const bool FractureMechanics::MWLSApprox::OpMWLSStressAtGaussPts::testing

Definition at line 408 of file MWLS.hpp.

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

users_modules/fracture_mechanics/src/MWLS.hpp
users_modules/fracture_mechanics/src/impl/MWLS.cpp

Public Member Functions

Public Attributes

Protected Member Functions