MechanicsLinearFEOperator.h

Go to the documentation of this file.

#ifndef included_AMP_MechanicsLinearFEOperator
#define included_AMP_MechanicsLinearFEOperator
 
/* AMP files */
#include "AMP/operators/libmesh/LinearFEOperator.h"
#include "AMP/operators/mechanics/MechanicsConstants.h"
#include "AMP/operators/mechanics/MechanicsLinearElement.h"
#include "AMP/operators/mechanics/MechanicsLinearFEOperatorParameters.h"
#include "AMP/operators/mechanics/MechanicsLinearUpdatedLagrangianElement.h"
#include "AMP/vectors/Variable.h"
 
#include <vector>
 
 
namespace AMP::Operator {
 
class MechanicsLinearFEOperator : public LinearFEOperator
{
public:
    explicit MechanicsLinearFEOperator( std::shared_ptr<const OperatorParameters> params );
 
    virtual ~MechanicsLinearFEOperator() {}
 
    void preAssembly( std::shared_ptr<const OperatorParameters> params ) override;
 
    void postAssembly() override;
 
    void preElementOperation( const AMP::Mesh::MeshElement & ) override;
 
    void postElementOperation() override;
 
    void printStressAndStrain( AMP::LinearAlgebra::Vector::shared_ptr disp,
                               const std::string &fname );
 
    std::shared_ptr<MechanicsMaterialModel> getMaterialModel() { return d_materialModel; }
 
protected:
    void getDofIndicesForCurrentElement();
    MechanicsLinearFEOperator( std::shared_ptr<const MechanicsLinearFEOperatorParameters>, bool );
 
protected:
    std::vector<std::vector<size_t>> d_dofIndices; 
    std::vector<std::vector<double>> d_elementStiffnessMatrix; 
    std::shared_ptr<MechanicsLinearElement> d_mechLinElem; 
    std::shared_ptr<MechanicsLinearUpdatedLagrangianElement>
        d_mechLinULElem; 
    std::shared_ptr<MechanicsMaterialModel> d_materialModel; 
    AMP::LinearAlgebra::Vector::shared_ptr d_refXYZ; 
    AMP::LinearAlgebra::Vector::shared_ptr d_dispVec;
 
    bool d_useUpdatedLagrangian;
};
} // namespace AMP::Operator
 
#endif