FractureMechanics::MWLSApprox::OpMWLSStressAtGaussUsingPrecalulatedCoeffs Struct Reference

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

Inheritance diagram for FractureMechanics::MWLSApprox::OpMWLSStressAtGaussUsingPrecalulatedCoeffs:

Collaboration diagram for FractureMechanics::MWLSApprox::OpMWLSStressAtGaussUsingPrecalulatedCoeffs:

Public Member Functions
	OpMWLSStressAtGaussUsingPrecalulatedCoeffs (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)
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
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

Definition at line 377 of file MWLS.hpp.

Constructor & Destructor Documentation

OpMWLSStressAtGaussUsingPrecalulatedCoeffs()

FractureMechanics::MWLSApprox::OpMWLSStressAtGaussUsingPrecalulatedCoeffs::OpMWLSStressAtGaussUsingPrecalulatedCoeffs ( 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 384 of file MWLS.hpp.

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

Member Function Documentation

doMWLSWork()

MoFEMErrorCode FractureMechanics::MWLSApprox::OpMWLSStressAtGaussUsingPrecalulatedCoeffs::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 1504 of file MWLS.cpp.

                                                                      {
  MoFEMFunctionBegin;
 
  if (testing)
    cerr << "Element " << getFEEntityHandle() << endl;
 
  const bool use_global_base =
      mwlsApprox->getUseGlobalBaseAtMaterialReferenceConfiguration();
  if (use_global_base)
    SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSSIBLE_CASE,
            "Coefficients can be precalculated only for local base in "
            "reference configuration");
 
  FTensor::Index<'i', 3> i;
 
  Tag th;
  CHKERR mwlsApprox->mwlsMoab.tag_get_handle(stressTagName.c_str(), th);
 
  const int nb_gauss_pts = data.getN().size1();
  auto &inv_AB_map = mwlsApprox->invABMap.at(getFEEntityHandle());
  auto &influence_nodes_map =
      mwlsApprox->influenceNodesMap.at(getFEEntityHandle());
  auto &dm_nodes_map = mwlsApprox->dmNodesMap.at(this->getFEEntityHandle());
 
  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) {
 
    FTensor::Tensor1<double, 3> t_mat_pos;
    t_mat_pos(i) = t_mwls_material_positions(i);
 
    if (testing) {
      cerr << "material_positions " << t_mwls_material_positions << endl;
    }
 
    for (int dd : {0, 1, 2})
      mwlsApprox->approxPointCoords[dd] = t_approx_base_point(dd);
 
    mwlsApprox->influenceNodes = influence_nodes_map[gg];
    mwlsApprox->invAB = inv_AB_map[gg];
    mwlsApprox->shortenDm = dm_nodes_map[gg];
 
    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::OpMWLSStressAtGaussUsingPrecalulatedCoeffs::calculateDerivative

Definition at line 381 of file MWLS.hpp.

stressTagName

std::string FractureMechanics::MWLSApprox::OpMWLSStressAtGaussUsingPrecalulatedCoeffs::stressTagName

Definition at line 380 of file MWLS.hpp.

testing

const bool FractureMechanics::MWLSApprox::OpMWLSStressAtGaussUsingPrecalulatedCoeffs::testing

Definition at line 382 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