CompressibleNonlinearElasticitySolver< DIM > Class Template Reference

#include <CompressibleNonlinearElasticitySolver.hpp>

Inheritance diagram for CompressibleNonlinearElasticitySolver< DIM >:

Inheritance graph
[legend]
Collaboration diagram for CompressibleNonlinearElasticitySolver< DIM >:

Collaboration graph
[legend]

List of all members.

Public Member Functions

 CompressibleNonlinearElasticitySolver (QuadraticMesh< DIM > *pQuadMesh, AbstractMaterialLaw< DIM > *pMaterialLaw, c_vector< double, DIM > bodyForce, double density, std::string outputDirectory, std::vector< unsigned > &fixedNodes, std::vector< c_vector< double, DIM > > *pFixedNodeLocations=NULL)
 CompressibleNonlinearElasticitySolver (QuadraticMesh< DIM > *pQuadMesh, std::vector< AbstractMaterialLaw< DIM > * > &rMaterialLaws, c_vector< double, DIM > bodyForce, double density, std::string outputDirectory, std::vector< unsigned > &fixedNodes, std::vector< c_vector< double, DIM > > *pFixedNodeLocations=NULL)
 ~CompressibleNonlinearElasticitySolver ()

Protected Member Functions

virtual void AssembleOnElement (Element< DIM, DIM > &rElement, c_matrix< double, STENCIL_SIZE, STENCIL_SIZE > &rAElem, c_matrix< double, STENCIL_SIZE, STENCIL_SIZE > &rAElemPrecond, c_vector< double, STENCIL_SIZE > &rBElem, bool assembleResidual, bool assembleJacobian)
virtual void AssembleOnBoundaryElement (BoundaryElement< DIM-1, DIM > &rBoundaryElement, c_matrix< double, BOUNDARY_STENCIL_SIZE, BOUNDARY_STENCIL_SIZE > &rAelem, c_vector< double, BOUNDARY_STENCIL_SIZE > &rBelem, c_vector< double, DIM > &rTraction, bool assembleResidual, bool assembleJacobian)
void AssembleSystem (bool assembleResidual, bool assembleJacobian)

Protected Attributes

std::vector
< AbstractCompressibleMaterialLaw
< DIM > * > 		mMaterialLaws

Static Protected Attributes

static const size_t NUM_NODES_PER_ELEMENT = AbstractNonlinearElasticitySolver<DIM>::NUM_NODES_PER_ELEMENT
static const size_t NUM_VERTICES_PER_ELEMENT = AbstractNonlinearElasticitySolver<DIM>::NUM_VERTICES_PER_ELEMENT
static const size_t NUM_NODES_PER_BOUNDARY_ELEMENT = AbstractNonlinearElasticitySolver<DIM>::NUM_NODES_PER_BOUNDARY_ELEMENT
static const size_t STENCIL_SIZE = DIM*NUM_NODES_PER_ELEMENT
static const size_t BOUNDARY_STENCIL_SIZE = DIM*NUM_NODES_PER_BOUNDARY_ELEMENT

Friends

class TestCompressibleNonlinearElasticitySolver

Detailed Description

template<size_t DIM>
class CompressibleNonlinearElasticitySolver< DIM >

Finite elasticity solver. Solves static *compressible* nonlinear elasticity problems with arbitrary (compressible) material laws and a body force.

Uses quadratic basis functions for displacement, and is therefore outside the other assembler or solver hierarchy.

Definition at line 56 of file CompressibleNonlinearElasticitySolver.hpp.

Constructor & Destructor Documentation

template<size_t DIM>
CompressibleNonlinearElasticitySolver< DIM >::CompressibleNonlinearElasticitySolver ( QuadraticMesh< DIM > *  pQuadMesh,
AbstractMaterialLaw< DIM > *  pMaterialLaw,
c_vector< double, DIM >  bodyForce,
double  density,
std::string  outputDirectory,
std::vector< unsigned > &  fixedNodes,
std::vector< c_vector< double, DIM > > *  pFixedNodeLocations = NULL 
) [inline]

Constructor for homogeneous problems.

Parameters:
pQuadMesh The quadratic mesh to solve on
pMaterialLaw A single material law to use on all elements
bodyForce The body force if constant. (If not constant, pass in a zero vector and call SetFunctionalBodyForce()
density The density (assumed constant)
outputDirectory The output directory
fixedNodes Which nodes are fixed in space (the displacement is assumed to be zero unless the next parameter is given
pFixedNodeLocations Optional new locations of the fixed nodes.

Definition at line 534 of file CompressibleNonlinearElasticitySolver.cpp.

References EXCEPTION, AbstractNonlinearElasticitySolver< DIM >::Initialise(), and CompressibleNonlinearElasticitySolver< DIM >::mMaterialLaws.

template<size_t DIM>
CompressibleNonlinearElasticitySolver< DIM >::CompressibleNonlinearElasticitySolver ( QuadraticMesh< DIM > *  pQuadMesh,
std::vector< AbstractMaterialLaw< DIM > * > &  rMaterialLaws,
c_vector< double, DIM >  bodyForce,
double  density,
std::string  outputDirectory,
std::vector< unsigned > &  fixedNodes,
std::vector< c_vector< double, DIM > > *  pFixedNodeLocations = NULL 
) [inline]

Variant constructor taking a vector of material laws for heterogeneous problems

Parameters:
pQuadMesh The quadratic mesh to solve on
rMaterialLaws Vector of material laws for each element
bodyForce The body force if constant. (If not constant, pass in a zero vector and call SetFunctionalBodyForce()
density The density (assumed constant)
outputDirectory The output directory
fixedNodes Which nodes are fixed in space (the displacement is assumed to be zero unless the next parameter is given
pFixedNodeLocations Optional new locations of the fixed nodes.

Definition at line 561 of file CompressibleNonlinearElasticitySolver.cpp.

References EXCEPTION, AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM >::GetNumElements(), AbstractNonlinearElasticitySolver< DIM >::Initialise(), and CompressibleNonlinearElasticitySolver< DIM >::mMaterialLaws.

template<size_t DIM>
CompressibleNonlinearElasticitySolver< DIM >::~CompressibleNonlinearElasticitySolver (  )  [inline]

Destructor.

Definition at line 592 of file CompressibleNonlinearElasticitySolver.cpp.

Member Function Documentation

template<size_t DIM>
void CompressibleNonlinearElasticitySolver< DIM >::AssembleOnElement ( Element< DIM, DIM > &  rElement,
c_matrix< double, STENCIL_SIZE, STENCIL_SIZE > &  rAElem,
c_matrix< double, STENCIL_SIZE, STENCIL_SIZE > &  rAElemPrecond,
c_vector< double, STENCIL_SIZE > &  rBElem,
bool  assembleResidual,
bool  assembleJacobian 
) [inline, protected, virtual]

Assemble residual or jacobian on an element, using the current solution stored in mCurrrentSolution. The ordering assumed is (in 2d) rBElem = [u0 v0 u1 v1 .. u5 v5].

Parameters:
rElement The element to assemble on.
rAElem The element's contribution to the LHS matrix is returned in this n by n matrix, where n is the no. of nodes in this element. There is no need to zero this matrix before calling.
rAElemPrecond The element's contribution to the matrix passed to PetSC in creating a preconditioner
rBElem The element's contribution to the RHS vector is returned in this vector of length n, the no. of nodes in this element. There is no need to zero this vector before calling.
assembleResidual A bool stating whether to assemble the residual vector.
assembleJacobian A bool stating whether to assemble the Jacobian matrix.

Definition at line 166 of file CompressibleNonlinearElasticitySolver.cpp.

References QuadraticBasisFunction< ELEMENT_DIM >::ComputeBasisFunctions(), LinearBasisFunction< ELEMENT_DIM >::ComputeBasisFunctions(), AbstractNonlinearElasticitySolver< DIM >::ComputeStressAndStressDerivative(), QuadraticBasisFunction< ELEMENT_DIM >::ComputeTransformedBasisFunctionDerivatives(), AbstractNonlinearElasticitySolver< DIM >::dTdE, AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetIndex(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNodeGlobalIndex(), GaussianQuadratureRule< ELEMENT_DIM >::GetNumQuadPoints(), GaussianQuadratureRule< ELEMENT_DIM >::GetWeight(), Inverse(), AbstractNonlinearElasticitySolver< DIM >::mBodyForce, AbstractNonlinearElasticitySolver< DIM >::mCurrentSolution, AbstractNonlinearElasticitySolver< DIM >::mCurrentTime, AbstractNonlinearElasticitySolver< DIM >::mDensity, CompressibleNonlinearElasticitySolver< DIM >::mMaterialLaws, AbstractNonlinearElasticitySolver< DIM >::mpBodyForceFunction, AbstractNonlinearElasticitySolver< DIM >::mpQuadMesh, AbstractNonlinearElasticitySolver< DIM >::mpQuadratureRule, AbstractNonlinearElasticitySolver< DIM >::mUsingBodyForceFunction, CompressibleNonlinearElasticitySolver< DIM >::NUM_NODES_PER_ELEMENT, CompressibleNonlinearElasticitySolver< DIM >::NUM_VERTICES_PER_ELEMENT, and GaussianQuadratureRule< ELEMENT_DIM >::rGetQuadPoint().

Referenced by CompressibleNonlinearElasticitySolver< DIM >::AssembleSystem().

template<size_t DIM>
void CompressibleNonlinearElasticitySolver< DIM >::AssembleOnBoundaryElement ( BoundaryElement< DIM-1, DIM > &  rBoundaryElement,
c_matrix< double, BOUNDARY_STENCIL_SIZE, BOUNDARY_STENCIL_SIZE > &  rAelem,
c_vector< double, BOUNDARY_STENCIL_SIZE > &  rBelem,
c_vector< double, DIM > &  rTraction,
bool  assembleResidual,
bool  assembleJacobian 
) [inline, protected, virtual]

Compute the term from the surface integral of s*phi, where s is a specified non-zero surface traction (ie Neumann boundary condition) to be added to the Rhs vector.

Parameters:
rBoundaryElement 
rAelem The element's contribution to the LHS matrix is returned in this n by n matrix, where n is the no. of nodes in this element. There is no need to zero this matrix before calling.
rBelem The element's contribution to the RHS vector is returned in this vector of length n, the no. of nodes in this element. There is no need to zero this vector before calling.
rTraction 
assembleResidual A bool stating whether to assemble the residual vector.
assembleJacobian A bool stating whether to assemble the Jacobian matrix.

Definition at line 468 of file CompressibleNonlinearElasticitySolver.cpp.

References QuadraticBasisFunction< ELEMENT_DIM >::ComputeBasisFunctions(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetIndex(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNodeGlobalIndex(), GaussianQuadratureRule< ELEMENT_DIM >::GetNumQuadPoints(), GaussianQuadratureRule< ELEMENT_DIM >::GetWeight(), AbstractNonlinearElasticitySolver< DIM >::mCurrentTime, AbstractNonlinearElasticitySolver< DIM >::mpBoundaryQuadratureRule, AbstractNonlinearElasticitySolver< DIM >::mpQuadMesh, AbstractNonlinearElasticitySolver< DIM >::mpTractionBoundaryConditionFunction, AbstractNonlinearElasticitySolver< DIM >::mUsingTractionBoundaryConditionFunction, CompressibleNonlinearElasticitySolver< DIM >::NUM_NODES_PER_BOUNDARY_ELEMENT, and GaussianQuadratureRule< ELEMENT_DIM >::rGetQuadPoint().

Referenced by CompressibleNonlinearElasticitySolver< DIM >::AssembleSystem().

template<size_t DIM>
void CompressibleNonlinearElasticitySolver< DIM >::AssembleSystem ( bool  assembleResidual,
bool  assembleJacobian 
) [inline, protected, virtual]

Assemble the residual vector (using the current solution stored in mCurrentSolution, output going to mpLinearSystem->rGetRhsVector), or Jacobian matrix (using the current solution stored in mCurrentSolution, output going to mpLinearSystem->rGetLhsMatrix).

Parameters:
assembleResidual A bool stating whether to assemble the residual vector.
assembleJacobian A bool stating whether to assemble the Jacobian matrix.

Implements AbstractNonlinearElasticitySolver< DIM >.

Definition at line 49 of file CompressibleNonlinearElasticitySolver.cpp.

References CompressibleNonlinearElasticitySolver< DIM >::AssembleOnBoundaryElement(), CompressibleNonlinearElasticitySolver< DIM >::AssembleOnElement(), CompressibleNonlinearElasticitySolver< DIM >::BOUNDARY_STENCIL_SIZE, AbstractNonlinearElasticitySolver< DIM >::FinishAssembleSystem(), PetscTools::GetMyRank(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNodeGlobalIndex(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetOwnership(), AbstractNonlinearElasticitySolver< DIM >::mBoundaryElements, AbstractNonlinearElasticitySolver< DIM >::mCurrentSolution, AbstractNonlinearElasticitySolver< DIM >::mJacobianMatrix, AbstractNonlinearElasticitySolver< DIM >::mNumDofs, AbstractNonlinearElasticitySolver< DIM >::mpQuadMesh, AbstractNonlinearElasticitySolver< DIM >::mPreconditionMatrix, AbstractNonlinearElasticitySolver< DIM >::mResidualVector, AbstractNonlinearElasticitySolver< DIM >::mSurfaceTractions, CompressibleNonlinearElasticitySolver< DIM >::NUM_NODES_PER_BOUNDARY_ELEMENT, CompressibleNonlinearElasticitySolver< DIM >::NUM_NODES_PER_ELEMENT, CompressibleNonlinearElasticitySolver< DIM >::STENCIL_SIZE, PetscMatTools::Zero(), and PetscVecTools::Zero().

Member Data Documentation

template<size_t DIM>
const size_t CompressibleNonlinearElasticitySolver< DIM >::NUM_NODES_PER_ELEMENT = AbstractNonlinearElasticitySolver<DIM>::NUM_NODES_PER_ELEMENT [static, protected]

Number of nodes per element

Reimplemented from AbstractNonlinearElasticitySolver< DIM >.

Definition at line 62 of file CompressibleNonlinearElasticitySolver.hpp.

Referenced by CompressibleNonlinearElasticitySolver< DIM >::AssembleOnElement(), and CompressibleNonlinearElasticitySolver< DIM >::AssembleSystem().

template<size_t DIM>
const size_t CompressibleNonlinearElasticitySolver< DIM >::NUM_VERTICES_PER_ELEMENT = AbstractNonlinearElasticitySolver<DIM>::NUM_VERTICES_PER_ELEMENT [static, protected]

Number of vertices per element

Reimplemented from AbstractNonlinearElasticitySolver< DIM >.

Definition at line 64 of file CompressibleNonlinearElasticitySolver.hpp.

Referenced by CompressibleNonlinearElasticitySolver< DIM >::AssembleOnElement().

template<size_t DIM>
const size_t CompressibleNonlinearElasticitySolver< DIM >::NUM_NODES_PER_BOUNDARY_ELEMENT = AbstractNonlinearElasticitySolver<DIM>::NUM_NODES_PER_BOUNDARY_ELEMENT [static, protected]

Number of nodes per boundary element

Reimplemented from AbstractNonlinearElasticitySolver< DIM >.

Definition at line 66 of file CompressibleNonlinearElasticitySolver.hpp.

Referenced by CompressibleNonlinearElasticitySolver< DIM >::AssembleOnBoundaryElement(), and CompressibleNonlinearElasticitySolver< DIM >::AssembleSystem().

template<size_t DIM>
const size_t CompressibleNonlinearElasticitySolver< DIM >::STENCIL_SIZE = DIM*NUM_NODES_PER_ELEMENT [static, protected]

Stencil size - number of unknowns per element (DIM*NUM_NODES_PER_ELEMENT displacement unknowns, no pressure unknowns

Definition at line 71 of file CompressibleNonlinearElasticitySolver.hpp.

Referenced by CompressibleNonlinearElasticitySolver< DIM >::AssembleSystem().

template<size_t DIM>
const size_t CompressibleNonlinearElasticitySolver< DIM >::BOUNDARY_STENCIL_SIZE = DIM*NUM_NODES_PER_BOUNDARY_ELEMENT [static, protected]

Boundary stencil size

Definition at line 73 of file CompressibleNonlinearElasticitySolver.hpp.

Referenced by CompressibleNonlinearElasticitySolver< DIM >::AssembleSystem().

template<size_t DIM>
std::vector<AbstractCompressibleMaterialLaw<DIM>*> CompressibleNonlinearElasticitySolver< DIM >::mMaterialLaws [protected]

The material laws for each element. This will either be of size 1 (same material law for all elements, ie homogeneous), or size num_elem.

Definition at line 80 of file CompressibleNonlinearElasticitySolver.hpp.

Referenced by CompressibleNonlinearElasticitySolver< DIM >::AssembleOnElement(), and CompressibleNonlinearElasticitySolver< DIM >::CompressibleNonlinearElasticitySolver().

The documentation for this class was generated from the following files:
Generated on Tue May 31 14:33:07 2011 for Chaste by  doxygen 1.5.5