FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX Struct Reference

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

Inheritance diagram for FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX:

Collaboration diagram for FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX:

Public Member Functions
	OpAleLhsWithDensitySingularElement_dX_dX (const std::string row_field, const std::string col_field, boost::shared_ptr< HookeElement::DataAtIntegrationPts > &data_at_pts, boost::shared_ptr< VectorDouble > rho_at_gauss_pts, boost::shared_ptr< MatrixDouble > rho_grad_at_gauss_pts, const double rho_n, const double rho_0, boost::shared_ptr< MatrixDouble > singular_displacement)

Public Attributes
boost::shared_ptr< VectorDouble >	rhoAtGaussPtsPtr
boost::shared_ptr< MatrixDouble >	rhoGradAtGaussPtsPtr
boost::shared_ptr< MatrixDouble >	singularDisplacement
const double	rhoN
const double	rHo0
MatrixDouble	K
MatrixDouble	transK
VectorDouble	nF
boost::shared_ptr< HookeElement::DataAtIntegrationPts >	dataAtPts
VectorInt	rowIndices
VectorInt	colIndices
int	nbRows
	number of dofs on rows More...
int	nbCols
	number if dof on column More...
int	nbIntegrationPts
	number of integration points More...
bool	isDiag
	true if this block is on diagonal More...

Protected Member Functions
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
MoFEMErrorCode	iNtegrate (DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
MoFEMErrorCode	aSsemble (DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)

Detailed Description

Definition at line 664 of file CrackFrontElement.hpp.

Constructor & Destructor Documentation

OpAleLhsWithDensitySingularElement_dX_dX()

FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::OpAleLhsWithDensitySingularElement_dX_dX	(	const std::string	row_field,
		const std::string	col_field,
		boost::shared_ptr< HookeElement::DataAtIntegrationPts > &	data_at_pts,
		boost::shared_ptr< VectorDouble >	rho_at_gauss_pts,
		boost::shared_ptr< MatrixDouble >	rho_grad_at_gauss_pts,
		const double	rho_n,
		const double	rho_0,
		boost::shared_ptr< MatrixDouble >	singular_displacement
	)

Definition at line 1187 of file CrackFrontElement.cpp.

    : UserDataOperator(row_field, col_field, OPROWCOL, true),
      dataAtPts(data_at_pts), rhoAtGaussPtsPtr(rho_at_gauss_pts),
      rhoGradAtGaussPtsPtr(rho_grad_at_gauss_pts),
      singularDisplacement(singular_displacement), rhoN(rho_n), rHo0(rho_0) {}

Member Function Documentation

aSsemble()

MoFEMErrorCode FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::aSsemble ( DataForcesAndSourcesCore::EntData &  row_data, DataForcesAndSourcesCore::EntData &  col_data  )

protected

Definition at line 1309 of file CrackFrontElement.cpp.

                                               {
  MoFEMFunctionBegin;
 
  const int *row_indices = &*row_data.getIndices().data().begin();
  const int *col_indices = &*col_data.getIndices().data().begin();
 
  auto &data = *dataAtPts;
  Mat B = getFEMethod()->ksp_B != PETSC_NULL ? getFEMethod()->ksp_B
                                             : getFEMethod()->snes_B;
  // assemble local matrix
  CHKERR MatSetValues(B, nbRows, row_indices, nbCols, col_indices,
                      &*K.data().begin(), ADD_VALUES);
 
  if (!isDiag && sYmm) {
    // if not diagonal term and since global matrix is symmetric assemble
    // transpose term.
    transK.resize(K.size2(), K.size1(), false);
    noalias(transK) = trans(K);
    CHKERR MatSetValues(B, nbCols, col_indices, nbRows, row_indices,
                        &*transK.data().begin(), ADD_VALUES);
  }
  MoFEMFunctionReturn(0);
}

doWork()

MoFEMErrorCode FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, DataForcesAndSourcesCore::EntData &  row_data, DataForcesAndSourcesCore::EntData &  col_data  )

protected

Definition at line 1199 of file CrackFrontElement.cpp.

                                               {
 
  MoFEMFunctionBegin;
  nbRows = row_data.getIndices().size();
  if (!nbRows)
    MoFEMFunctionReturnHot(0);
 
  nbCols = col_data.getIndices().size();
  if (!nbCols)
    MoFEMFunctionReturnHot(0);
 
  K.resize(nbRows, nbCols, false);
  K.clear();
 
  nbIntegrationPts = getGaussPts().size2();
  if (row_side == col_side && row_type == col_type) {
    isDiag = true;
  } else {
    isDiag = false;
  }
 
  CHKERR iNtegrate(row_data, col_data);
  CHKERR aSsemble(row_data, col_data);
 
  MoFEMFunctionReturn(0);
}

iNtegrate()

MoFEMErrorCode FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::iNtegrate ( DataForcesAndSourcesCore::EntData &  row_data, DataForcesAndSourcesCore::EntData &  col_data  )

protected

Definition at line 1229 of file CrackFrontElement.cpp.

                                               {
  MoFEMFunctionBegin;
  auto get_tensor2 = [](MatrixDouble &m, const int r, const int c) {
    return FTensor::Tensor2<FTensor::PackPtr<double *, 3>, 3, 3>(
        &m(r + 0, c + 0), &m(r + 0, c + 1), &m(r + 0, c + 2), &m(r + 1, c + 0),
        &m(r + 1, c + 1), &m(r + 1, c + 2), &m(r + 2, c + 0), &m(r + 2, c + 1),
        &m(r + 2, c + 2));
  };
 
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
  FTensor::Index<'l', 3> l;
 
  double vol = getVolume();
 
  auto t_w = getFTensor0IntegrationWeight();
 
  auto t_row_diff_base = row_data.getFTensor1DiffN<3>();
  const int row_nb_base_fun = row_data.getN().size2();
 
  auto rho = getFTensor0FromVec(*rhoAtGaussPtsPtr);
  auto t_grad_rho = getFTensor1FromMat<3>(*rhoGradAtGaussPtsPtr);
 
  auto t_eshelby_stress =
      getFTensor2FromMat<3, 3>(*dataAtPts->eshelbyStressMat);
  auto t_invH = getFTensor2FromMat<3, 3>(*dataAtPts->invHMat);
  auto &det_H = *dataAtPts->detHVec;
 
  auto t_initial_singular_displacement =
      getFTensor1FromMat<3>(*singularDisplacement);
 
  for (int gg = 0; gg != nbIntegrationPts; ++gg) {
 
    double a = t_w * vol * det_H[gg];
    const double stress_dho_coef = (rhoN / rho);
 
    int rr = 0;
    for (; rr != nbRows / 3; ++rr) {
 
      auto t_m = get_tensor2(K, 3 * rr, 0);
 
      FTensor::Tensor1<double, 3> t_row_diff_base_pulled;
      t_row_diff_base_pulled(i) = t_row_diff_base(j) * t_invH(j, i);
 
      FTensor::Tensor1<double, 3> t_row_stress;
      t_row_stress(i) = a * t_row_diff_base_pulled(j) * t_eshelby_stress(i, j);
      FTensor::Tensor2<double, 3, 3> t_a;
      t_a(i, k) = t_row_stress(i) * stress_dho_coef * t_grad_rho(k);
 
      auto t_col_diff_base = col_data.getFTensor1DiffN<3>(gg, 0);
      auto t_col_base = col_data.getFTensor0N(gg, 0);
      for (int cc = 0; cc != nbCols / 3; ++cc) {
 
        t_m(i, k) +=
            t_a(i, k) * t_col_diff_base(l) * t_initial_singular_displacement(l);
 
        ++t_col_base;
        ++t_col_diff_base;
        ++t_m;
      }
      ++t_row_diff_base;
    }
 
    for (; rr != row_nb_base_fun; ++rr)
      ++t_row_diff_base;
 
    ++t_initial_singular_displacement;
    ++t_w;
    ++t_eshelby_stress;
    ++t_invH;
    ++rho;
    ++t_grad_rho;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

colIndices

VectorInt FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::colIndices

Definition at line 682 of file CrackFrontElement.hpp.

dataAtPts

boost::shared_ptr<HookeElement::DataAtIntegrationPts> FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::dataAtPts

Definition at line 679 of file CrackFrontElement.hpp.

isDiag

bool FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::isDiag

true if this block is on diagonal

Definition at line 687 of file CrackFrontElement.hpp.

K

MatrixDouble FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::K

Definition at line 675 of file CrackFrontElement.hpp.

nbCols

int FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::nbCols

number if dof on column

Definition at line 685 of file CrackFrontElement.hpp.

nbIntegrationPts

int FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::nbIntegrationPts

number of integration points

Definition at line 686 of file CrackFrontElement.hpp.

nbRows

int FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::nbRows

number of dofs on rows

Definition at line 684 of file CrackFrontElement.hpp.

nF

VectorDouble FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::nF

Definition at line 677 of file CrackFrontElement.hpp.

rHo0

const double FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::rHo0

Definition at line 672 of file CrackFrontElement.hpp.

rhoAtGaussPtsPtr

boost::shared_ptr<VectorDouble> FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::rhoAtGaussPtsPtr

Definition at line 667 of file CrackFrontElement.hpp.

rhoGradAtGaussPtsPtr

boost::shared_ptr<MatrixDouble> FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::rhoGradAtGaussPtsPtr

Definition at line 668 of file CrackFrontElement.hpp.

rhoN

const double FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::rhoN

Definition at line 671 of file CrackFrontElement.hpp.

rowIndices

VectorInt FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::rowIndices

Definition at line 681 of file CrackFrontElement.hpp.

singularDisplacement

boost::shared_ptr<MatrixDouble> FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::singularDisplacement

Definition at line 669 of file CrackFrontElement.hpp.

transK

MatrixDouble FractureMechanics::OpAleLhsWithDensitySingularElement_dX_dX::transK

Definition at line 676 of file CrackFrontElement.hpp.

---

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

• users_modules/fracture_mechanics/src/CrackFrontElement.hpp
• users_modules/fracture_mechanics/src/impl/CrackFrontElement.cpp