SphericalArcLengthControl Struct Reference

Implementation of spherical arc-length method. More...

#include <users_modules/basic_finite_elements/src/ArcLengthTools.hpp>

Inheritance diagram for SphericalArcLengthControl:

Collaboration diagram for SphericalArcLengthControl:

Public Member Functions
DEPRECATED	SphericalArcLengthControl (ArcLengthCtx *arc_ptr_raw)
	SphericalArcLengthControl (boost::shared_ptr< ArcLengthCtx > &arc_ptr)
virtual	~SphericalArcLengthControl ()
MoFEMErrorCode	preProcess ()
MoFEMErrorCode	operator() ()
MoFEMErrorCode	postProcess ()
virtual double	calculateLambdaInt ()
	Calculate f_lambda(dx,lambda) More...
virtual MoFEMErrorCode	calculateDb ()
	Calculate db. More...
virtual MoFEMErrorCode	calculateDxAndDlambda (Vec x)
virtual MoFEMErrorCode	calculateInitDlambda (double *dlambda)
virtual MoFEMErrorCode	setDlambdaToX (Vec x, double dlambda)

Public Attributes
ArcLengthCtx *	arcPtrRaw

Private Attributes
boost::shared_ptr< ArcLengthCtx >	arcPtr

Detailed Description

Implementation of spherical arc-length method.

\[ \alpha \| \Delta\mathbf{x} \|^2 + \Delta\lambda^2 \beta^2 \| \mathbf{F}_\lambda \|^2 = s^2 \]

This is particular implementation of ArcLength control, i.e. spherical arc length control. If beta is set to 0 and alpha is non-zero it is cylindrical arc-length control. Works well with general problem with non-linear geometry. It not guarantee dissipative loading path in case of physical nonlinearities.

Definition at line 374 of file ArcLengthTools.hpp.

Constructor & Destructor Documentation

SphericalArcLengthControl() [1/2]

SphericalArcLengthControl::SphericalArcLengthControl

(

ArcLengthCtx *

arc_ptr_raw

)

Deprecated:

use constructor with shared_ptr

Definition at line 661 of file ArcLengthTools.cpp.

    : FEMethod(), arcPtrRaw(arc_ptr_raw) {}

SphericalArcLengthControl() [2/2]

SphericalArcLengthControl::SphericalArcLengthControl

(

boost::shared_ptr< ArcLengthCtx > &

arc_ptr

)

Definition at line 664 of file ArcLengthTools.cpp.

    : FEMethod(), arcPtrRaw(arc_ptr.get()), arcPtr(arc_ptr) {}

~SphericalArcLengthControl()

SphericalArcLengthControl::~SphericalArcLengthControl ( )

virtual

Definition at line 668 of file ArcLengthTools.cpp.

{}

Member Function Documentation

calculateDb()

MoFEMErrorCode SphericalArcLengthControl::calculateDb ( )

virtual

Calculate db.

\[ \textrm{d}\mathbf{b} = 2 \alpha \Delta\mathbf{x} \]

Reimplemented in ShellArcLength.

Definition at line 691 of file ArcLengthTools.cpp.

                                                      {
  MoFEMFunctionBegin;
  CHKERR VecCopy(arcPtrRaw->dx, arcPtrRaw->db);
  CHKERR VecScale(arcPtrRaw->db, 2 * arcPtrRaw->alpha);
  MoFEMFunctionReturn(0);
}

calculateDxAndDlambda()

MoFEMErrorCode SphericalArcLengthControl::calculateDxAndDlambda ( Vec  x )

virtual

Definition at line 742 of file ArcLengthTools.cpp.

                                                                     {
  ArcFunctionBegin;
  // dx
  CHKERR VecCopy(x, arcPtrRaw->dx);
  CHKERR VecAXPY(arcPtrRaw->dx, -1, arcPtrRaw->x0);
  CHKERR VecGhostUpdateBegin(arcPtrRaw->dx, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(arcPtrRaw->dx, INSERT_VALUES, SCATTER_FORWARD);
  // dlambda
  if (arcPtrRaw->getPetscLocalDofIdx() != -1) {
    double *array;
    CHKERR VecGetArray(arcPtrRaw->dx, &array);
    arcPtrRaw->dLambda = array[arcPtrRaw->getPetscLocalDofIdx()];
    array[arcPtrRaw->getPetscLocalDofIdx()] = 0;
    CHKERR VecRestoreArray(arcPtrRaw->dx, &array);
  }
  CHKERR VecGhostUpdateBegin(arcPtrRaw->ghosTdLambda, INSERT_VALUES,
                             SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(arcPtrRaw->ghosTdLambda, INSERT_VALUES,
                           SCATTER_FORWARD);
  // dx2
  CHKERR VecDot(arcPtrRaw->dx, arcPtrRaw->dx, &arcPtrRaw->dx2);
  MOFEM_LOG_C("WORLD", Sev::verbose, "\tdlambda = %6.4e dx2 = %6.4e\n",
              arcPtrRaw->dLambda, arcPtrRaw->dx2);
  MoFEMFunctionReturn(0);
}

calculateInitDlambda()

MoFEMErrorCode SphericalArcLengthControl::calculateInitDlambda ( double *  dlambda )

virtual

Definition at line 769 of file ArcLengthTools.cpp.

                                                               {
  ArcFunctionBegin;
  *dlambda = std::sqrt(pow(arcPtrRaw->s, 2) /
                  (pow(arcPtrRaw->beta, 2) * arcPtrRaw->F_lambda2));
  if (!(*dlambda == *dlambda)) {
    MOFEM_LOG("WORLD", Sev::error)
        << "s " << arcPtrRaw->s << " " << arcPtrRaw->beta << " "
        << arcPtrRaw->F_lambda2;
    SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSSIBLE_CASE,
            "Increment of lambda is not a number");
  }
  MoFEMFunctionReturn(0);
}

calculateLambdaInt()

double SphericalArcLengthControl::calculateLambdaInt ( )

virtual

Calculate f_lambda(dx,lambda)

\[ f_\lambda(\Delta\mathbf{x},\lambda) = \alpha \| \Delta\mathbf{x} \|^2 + \Delta\lambda^2 \beta^2 \| \mathbf{F}_\lambda \|^2 \]

Reimplemented in ShellArcLength.

Definition at line 685 of file ArcLengthTools.cpp.

                                                     {
  return arcPtrRaw->alpha * arcPtrRaw->dx2 + pow(arcPtrRaw->dLambda, 2) *
                                                 pow(arcPtrRaw->beta, 2) *
                                                 arcPtrRaw->F_lambda2;
}

operator()()

MoFEMErrorCode SphericalArcLengthControl::operator()

(

)

Definition at line 698 of file ArcLengthTools.cpp.

                                                     {
  ArcFunctionBegin;
  switch (snes_ctx) {
  case CTX_SNESSETFUNCTION: {
    arcPtrRaw->res_lambda = calculateLambdaInt() - pow(arcPtrRaw->s, 2);
    CHKERR VecSetValue(snes_f, arcPtrRaw->getPetscGlobalDofIdx(),
                       arcPtrRaw->res_lambda, ADD_VALUES);
    MOFEM_LOG_C("WORLD", Sev::verbose, "\tres_lambda = %6.4e\n",
                arcPtrRaw->res_lambda);
  } break;
  case CTX_SNESSETJACOBIAN: {
    arcPtrRaw->dIag =
        2 * arcPtrRaw->dLambda * pow(arcPtrRaw->beta, 2) * arcPtrRaw->F_lambda2;
    CHKERR MatSetValue(snes_B, arcPtrRaw->getPetscGlobalDofIdx(),
                       arcPtrRaw->getPetscGlobalDofIdx(), 1, ADD_VALUES);
  } break;
  default:
    break;
  }
  MoFEMFunctionReturn(0);
}

postProcess()

MoFEMErrorCode SphericalArcLengthControl::postProcess

(

)

Definition at line 720 of file ArcLengthTools.cpp.

                                                      {
  ArcFunctionBegin;
  switch (snes_ctx) {
  case CTX_SNESSETFUNCTION: {
    MOFEM_LOG_C("WORLD", Sev::verbose, "\tlambda = %6.4e\n",
                arcPtrRaw->getFieldData());
  } break;
  case CTX_SNESSETJACOBIAN: {
    CHKERR VecGhostUpdateBegin(arcPtrRaw->ghostDiag, INSERT_VALUES,
                               SCATTER_FORWARD);
    CHKERR VecGhostUpdateEnd(arcPtrRaw->ghostDiag, INSERT_VALUES,
                             SCATTER_FORWARD);
    CHKERR MatAssemblyBegin(snes_B, MAT_FLUSH_ASSEMBLY);
    CHKERR MatAssemblyEnd(snes_B, MAT_FLUSH_ASSEMBLY);
    MOFEM_LOG_C("WORLD", Sev::verbose, "\tdiag = %6.4e", arcPtrRaw->dIag);
  } break;
  default:
    break;
  }
  MoFEMFunctionReturn(0);
}

preProcess()

MoFEMErrorCode SphericalArcLengthControl::preProcess

(

)

Definition at line 670 of file ArcLengthTools.cpp.

                                                     {
  MoFEMFunctionBegin;
  switch (snes_ctx) {
  case CTX_SNESSETFUNCTION: {
    CHKERR calculateDxAndDlambda(snes_x);
    CHKERR calculateDb();
  } break;
  case CTX_SNESSETJACOBIAN: {
  } break;
  default:
    break;
  }
  MoFEMFunctionReturn(0);
}

setDlambdaToX()

MoFEMErrorCode SphericalArcLengthControl::setDlambdaToX ( Vec  x, double  dlambda  )

virtual

Definition at line 783 of file ArcLengthTools.cpp.

                                                                             {
  ArcFunctionBegin;
  // check if local dof idx is non zero, i.e. that lambda is accessible from
  // this processor
  if (arcPtrRaw->getPetscLocalDofIdx() != -1) {
    double *array;
    CHKERR VecGetArray(x, &array);
    double lambda_old = array[arcPtrRaw->getPetscLocalDofIdx()];
    if (!(dlambda == dlambda)) {
      MOFEM_LOG("WORLD", Sev::error)
          << "s " << arcPtrRaw->s << " " << arcPtrRaw->beta << " "
          << arcPtrRaw->F_lambda2;
      SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSSIBLE_CASE,
              "Increment of lambda is not a number");
    }
    array[arcPtrRaw->getPetscLocalDofIdx()] = lambda_old + dlambda;
    MOFEM_LOG_C("WORLD", Sev::verbose, "\tlambda = %6.4e, %6.4e (%6.4e)\n",
                lambda_old, array[arcPtrRaw->getPetscLocalDofIdx()], dlambda);
    CHKERR VecRestoreArray(x, &array);
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

arcPtr

boost::shared_ptr<ArcLengthCtx> SphericalArcLengthControl::arcPtr

private

Definition at line 412 of file ArcLengthTools.hpp.

arcPtrRaw

ArcLengthCtx* SphericalArcLengthControl::arcPtrRaw

Definition at line 376 of file ArcLengthTools.hpp.

---

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

• users_modules/basic_finite_elements/src/ArcLengthTools.hpp
• users_modules/basic_finite_elements/src/impl/ArcLengthTools.cpp