ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain Struct Reference

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

Inheritance diagram for ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain:

Collaboration diagram for ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain:

Public Member Functions
	TangentWithMeshSmoothingFrontConstrain (MoFEM::Interface &m_field, string lambda_field_name, int verbose=0)
virtual	~TangentWithMeshSmoothingFrontConstrain ()
PetscErrorCode	preProcess ()
PetscErrorCode	operator() ()
PetscErrorCode	postProcess ()
Public Member Functions inherited from ObosleteUsersModules::ConstrainConstantAarea
	ConstrainConstantAarea (MoFEM::Interface &m_field, Mat _C, Mat _Q, string lambda_field_name, int _verbose=0)
virtual	~ConstrainConstantAarea ()=default
template<typename DOFS >
auto	loIt (DOFS dofs, const FieldBitNumber bit, const EntityHandle ent)
template<typename DOFS >
auto	hiIt (DOFS dofs, const FieldBitNumber bit, const EntityHandle ent)
PetscErrorCode	getData (bool is_that_C_otherwise_dC, bool trans)
	get face data, indices and coors and nodal values More...
PetscErrorCode	calcDirevatives (double *diffNTRI, double *dofs_X, double *dofs_iX, double *C, double *iC, double *T, double *iT)
	calculate direvatives More...
PetscErrorCode	preProcess ()
PetscErrorCode	operator() ()
PetscErrorCode	postProcess ()

Public Attributes
Vec	frontF
Vec	tangentFrontF
const double	eps
bool	ownVectors
Tag	thFrontTangent
Public Attributes inherited from ObosleteUsersModules::ConstrainConstantAarea
MoFEM::Interface &	mField
moab::Interface &	moab
Mat	C
Mat	Q
string	lambdaFieldName
int	verbose
MatrixDouble	diffNTRI
const double *	G_TRI_W
vector< DofIdx >	DirichletBC
VectorInt	dispDofsColIdx
VectorInt	dispDofsRowIdx
VectorInt	localDispDofsRowIdx
VectorInt	lambdaDofsRowIdx
VectorInt	lambdaDofsColIdx
ublas::vector< double, ublas::bounded_array< double, 9 > >	coords
ublas::vector< double, ublas::bounded_array< double, 9 > >	dofs_X
ublas::vector< double, ublas::bounded_array< double, 3 > >	lambda
EntityHandle	face
Range	problemTets
map< DofIdx, Vec >	mapV

Detailed Description

Constrains is vector_tangent_to_front*nodal_element_quality_at_crack_front = 0, multiplying that constrain by Legrange multipliers, we get

1) vector_tangent_to_front*nodal_element_quality_at_crack_front=0 for each lambda, 2) lambda*vector_tangent_to_front*direvatives_of_quality = 0

Those two equations need to be linearsed and assembled into system of equations

Definition at line 583 of file ComplexConstArea.hpp.

Constructor & Destructor Documentation

TangentWithMeshSmoothingFrontConstrain()

ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::TangentWithMeshSmoothingFrontConstrain ( MoFEM::Interface &  m_field, string  lambda_field_name, int  verbose = 0  )

inline

Definition at line 590 of file ComplexConstArea.hpp.

      : ConstrainConstantAarea(m_field, PETSC_NULL, PETSC_NULL,
                               lambda_field_name, verbose),
        frontF(PETSC_NULL), tangentFrontF(PETSC_NULL), eps(1e-10),
        ownVectors(false) {}

~TangentWithMeshSmoothingFrontConstrain()

virtual ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::~TangentWithMeshSmoothingFrontConstrain ( )

inlinevirtual

Definition at line 598 of file ComplexConstArea.hpp.

                                                    {
    if (ownVectors) {
      CHKERR VecDestroy(&frontF);
      CHKERRABORT(PETSC_COMM_WORLD, ierr);
      CHKERR VecDestroy(&tangentFrontF);
      CHKERRABORT(PETSC_COMM_WORLD, ierr);
    }
  }

Member Function Documentation

operator()()

PetscErrorCode ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::operator() ( )

inline

Definition at line 648 of file ComplexConstArea.hpp.

                              {
    MoFEMFunctionBegin;
    // get elemet hanlde
    EntityHandle face = numeredEntFiniteElementPtr->getEnt();
    // get adjacent tets
    Range tet;
    CHKERR mField.get_moab().get_adjacencies(&face, 1, 3, false, tet);
    // get data and indices on dofs
    CHKERR getData(true, true);
 
    // Calculate tangent vector constrains
    ublas::vector<double, ublas::bounded_array<double, 9>> tangent_constrains(
        9);
    CHKERR calcDirevatives(&*diffNTRI.data().begin(), &*dofs_X.data().begin(),
                           NULL, NULL, NULL,
                           &*tangent_constrains.data().begin(), NULL);
    // take in account face orientation in respect crack surface
    Tag th_interface_side;
    CHKERR moab.tag_get_handle("INTERFACE_SIDE", th_interface_side);
    int side;
    CHKERR moab.tag_get_data(th_interface_side, &face, 1, &side);
    if (side == 1) {
      tangent_constrains *= -1;
    }
    // get smoothing element forces
    ublas::vector<double, ublas::bounded_array<double, 9>>
        f_front_mesh_smoothing(9);
    ublas::noalias(f_front_mesh_smoothing) = ublas::zero_vector<double>(9);
    double *f_front_mesh_array;
    CHKERR VecGetArray(frontF, &f_front_mesh_array);
    for (int nn = 0; nn < 3; nn++) {
      for (int dd = 0; dd < 3; dd++) {
        if (localDispDofsRowIdx[3 * nn + dd] == -1)
          continue;
        f_front_mesh_smoothing[3 * nn + dd] =
            f_front_mesh_array[localDispDofsRowIdx[3 * nn + dd]];
      }
    }
    CHKERR VecRestoreArray(frontF, &f_front_mesh_array);
    // calculate tangent matrix
    if (snes_ctx == CTX_SNESSETJACOBIAN) {
      // get radius
      double center[3];
      tricircumcenter3d_tp(&coords.data()[0], &coords.data()[3],
                           &coords.data()[6], center, NULL, NULL);
      cblas_daxpy(3, -1, &coords.data()[0], 1, center, 1);
      double r = cblas_dnrm2(3, center, 1);
      for (int nn = 0; nn < 3; nn++) {
        for (int dd = 0; dd < 3; dd++) {
          // ---> calculate tangent starts here
          ublas::vector<double, ublas::bounded_array<double, 9>> idofs_X(9, 0);
          idofs_X[nn * 3 + dd] = r * eps;
          ublas::vector<double, ublas::bounded_array<double, 9>>
              direvative_of_constrains(9);
          // use complex derivative
          CHKERR calcDirevatives(&*diffNTRI.data().begin(),
                                 &*dofs_X.data().begin(),
                                 &*idofs_X.data().begin(), NULL, NULL, NULL,
                                 &*direvative_of_constrains.data().begin());
          if (side == 1) {
            direvative_of_constrains /= -r * eps;
          } else {
            direvative_of_constrains /= +r * eps;
          }
          // dg -> C*q_quality
          double g[3] = {0, 0, 0};
          for (int nnn = 0; nnn < 3; nnn++) {
            if (lambdaDofsRowIdx[nnn] == -1)
              continue;
            for (int ddd = 0; ddd < 3; ddd++) {
              g[nnn] += direvative_of_constrains[nnn * 3 + ddd] *
                        f_front_mesh_smoothing[3 * nnn + ddd];
            }
          }
          for (int nnn = 0; nnn < 3; nnn++) {
            if (lambdaDofsRowIdx[nnn] == -1)
              continue;
            CHKERR MatSetValues(snes_B, 1, &lambdaDofsRowIdx[nnn], 1,
                                &dispDofsColIdx[3 * nn + dd], &g[nnn],
                                ADD_VALUES);
          }
          // dCT_lambda
          for (int nnn = 0; nnn < 3; nnn++) {
            for (int ddd = 0; ddd < 3; ddd++) {
              direvative_of_constrains[nnn * 3 + ddd] *= lambda[nnn];
            }
          }
          for (int nnn = 0; nnn < 3; nnn++) {
            if (lambdaDofsRowIdx[nnn] == -1)
              continue;
            CHKERR MatSetValues(snes_B, 3, &dispDofsRowIdx[3 * nnn], 1,
                                &dispDofsColIdx[3 * nn + dd],
                                &direvative_of_constrains[3 * nnn], ADD_VALUES);
          }
          // ---> calculate tangent end here
        }
      }
    }
    switch (snes_ctx) {
    case CTX_SNESSETFUNCTION: {
      ublas::vector<double, ublas::bounded_array<double, 3>> g(3);
      for (int nn = 0; nn < 3; nn++) {
        g[nn] = cblas_ddot(3, &tangent_constrains[3 * nn], 1,
                           &f_front_mesh_smoothing[3 * nn], 1);
      }
      // cerr << "g : " << g << endl;
      CHKERR VecSetValues(snes_f, 3, &*lambdaDofsRowIdx.data().begin(),
                          &*g.data().begin(), ADD_VALUES);
      CHKERR VecSetValues(tangentFrontF, 9, &dispDofsRowIdx[0],
                          &*tangent_constrains.data().begin(), ADD_VALUES);
      ublas::vector<double, ublas::bounded_array<double, 9>> f(9);
      for (int nn = 0; nn < 3; nn++) {
        for (int dd = 0; dd < 3; dd++) {
          f[nn * 3 + dd] = lambda[nn] * tangent_constrains[3 * nn + dd];
        }
      }
      // cerr << "f : " << f << endl;
      CHKERR VecSetValues(snes_f, 9, &dispDofsRowIdx[0], &*f.data().begin(),
                          ADD_VALUES);
      /*//TAG - only for one proc analysis
      for(int nn = 0;nn<3;nn++) {
      EntityHandle ent = lambda_dofs_row_ents[nn];
      if(ent == no_handle) continue;
      double *t;
      CHKERR mField.get_moab().tag_get_by_ptr(thFrontTangent,&ent,1,(const void
    **)&t); cblas_daxpy(3,+1,&tangent_constrains[3*nn],1,t,1);
    }*/
    } break;
    case CTX_SNESSETJACOBIAN: {
      for (int nn = 0; nn < 3; nn++) {
        int lambda_dof_idx = lambdaDofsColIdx[nn];
        CHKERR MatSetValues(snes_B, 3, &dispDofsRowIdx[3 * nn], 1,
                            &lambda_dof_idx, &tangent_constrains[3 * nn],
                            ADD_VALUES);
      }
    } break;
    default:
      break;
    }
    MoFEMFunctionReturn(0);
  }

postProcess()

PetscErrorCode ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::postProcess ( )

inline

Definition at line 790 of file ComplexConstArea.hpp.

                               {
    MoFEMFunctionBegin;
    switch (snes_ctx) {
    case CTX_SNESSETFUNCTION: {
      CHKERR VecAssemblyBegin(tangentFrontF);
      CHKERR VecAssemblyEnd(tangentFrontF);
      CHKERR VecGhostUpdateBegin(tangentFrontF, ADD_VALUES, SCATTER_REVERSE);
      CHKERR VecGhostUpdateEnd(tangentFrontF, ADD_VALUES, SCATTER_REVERSE);
      CHKERR VecGhostUpdateBegin(tangentFrontF, INSERT_VALUES, SCATTER_FORWARD);
      CHKERR VecGhostUpdateEnd(tangentFrontF, INSERT_VALUES, SCATTER_FORWARD);
    } break;
    case CTX_SNESSETJACOBIAN: {
      CHKERR MatAssemblyBegin(snes_B, MAT_FLUSH_ASSEMBLY);
      CHKERR MatAssemblyEnd(snes_B, MAT_FLUSH_ASSEMBLY);
    } break;
    default: {
    } break;
    }
    MoFEMFunctionReturn(0);
  }

preProcess()

PetscErrorCode ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::preProcess ( )

inline

Definition at line 608 of file ComplexConstArea.hpp.

                              {
    MoFEMFunctionBegin;
 
    switch (ts_ctx) {
    case CTX_TSSETIFUNCTION: {
      snes_ctx = CTX_SNESSETFUNCTION;
      snes_f = ts_F;
      break;
    }
    case CTX_TSSETIJACOBIAN: {
      snes_ctx = CTX_SNESSETJACOBIAN;
      snes_B = ts_B;
      break;
    }
    default:
      break;
    }
 
    CHKERR ConstrainConstantAarea::preProcess();
 
    switch (snes_ctx) {
    case CTX_SNESSETFUNCTION: {
      CHKERR VecAssemblyBegin(snes_f);
      CHKERR VecAssemblyEnd(snes_f);
      CHKERR VecZeroEntries(tangentFrontF);
      CHKERR VecGhostUpdateBegin(tangentFrontF, INSERT_VALUES, SCATTER_FORWARD);
      CHKERR VecGhostUpdateEnd(tangentFrontF, INSERT_VALUES, SCATTER_FORWARD);
      CHKERR VecAssemblyBegin(tangentFrontF);
      CHKERR VecAssemblyEnd(tangentFrontF);
    } break;
    case CTX_SNESSETJACOBIAN: {
      CHKERR MatAssemblyBegin(snes_B, MAT_FLUSH_ASSEMBLY);
      CHKERR MatAssemblyEnd(snes_B, MAT_FLUSH_ASSEMBLY);
    } break;
    default: {
    } break;
    }
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

eps

const double ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::eps

Definition at line 587 of file ComplexConstArea.hpp.

frontF

Vec ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::frontF

Definition at line 585 of file ComplexConstArea.hpp.

ownVectors

bool ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::ownVectors

Definition at line 588 of file ComplexConstArea.hpp.

tangentFrontF

Vec ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::tangentFrontF

Definition at line 586 of file ComplexConstArea.hpp.

thFrontTangent

Tag ObosleteUsersModules::TangentWithMeshSmoothingFrontConstrain::thFrontTangent

Definition at line 607 of file ComplexConstArea.hpp.

---

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

• users_modules/fracture_mechanics/src/ComplexConstArea.hpp