MechanicsLinearElement.h

Go to the documentation of this file.

#ifndef included_AMP_MechanicsLinearElement
#define included_AMP_MechanicsLinearElement
 
#include "AMP/operators/mechanics/MechanicsElement.h"
 
#include <memory>
#include <vector>
 
 
namespace AMP::Operator {
 
class MechanicsLinearElement : public MechanicsElement
{
public:
    explicit MechanicsLinearElement( std::shared_ptr<const ElementOperationParameters> params )
        : MechanicsElement( params ), d_elementStiffnessMatrix( nullptr )
    {
        d_JxW  = &( d_fe->get_JxW() );
        d_dphi = &( d_fe->get_dphi() );
        d_xyz  = &( d_fe->get_xyz() );
        AMP_INSIST( ( d_useJaumannRate == false ),
                    "Jaumann rate with small strain does not make any sense." );
    }
 
    virtual ~MechanicsLinearElement() {}
 
    void setElementStiffnessMatrix( std::vector<std::vector<double>> &elementStiffnessMatrix )
    {
        d_elementStiffnessMatrix = &( elementStiffnessMatrix );
    }
 
    void apply() override
    {
        if ( d_useReducedIntegration ) {
            apply_Reduced();
        } else {
            apply_Normal();
        }
    }
 
    void printStressAndStrain( FILE *fp, const std::vector<double> &dispVec );
 
    void computeStressAndStrain( const std::vector<double> &dispVec,
                                 std::vector<double> &stressVec,
                                 std::vector<double> &strainVec );
 
protected:
    void apply_Normal();
 
    void apply_Reduced();
 
    const std::vector<libMesh::Real> *d_JxW; 
    const std::vector<std::vector<libMesh::RealGradient>>
        *d_dphi; 
    const std::vector<libMesh::Point>
        *d_xyz; 
    std::vector<std::vector<double>> *d_elementStiffnessMatrix; 
private:
};
} // namespace AMP::Operator
 
#endif