MechanicsNonlinearElement.h

Go to the documentation of this file.

#ifndef included_AMP_MechanicsNonlinearElement
#define included_AMP_MechanicsNonlinearElement
 
#include "AMP/operators/mechanics/MechanicsConstants.h"
#include "AMP/operators/mechanics/MechanicsElement.h"
 
#include <memory>
#include <vector>
 
 
namespace AMP::Operator {
 
 
class MechanicsNonlinearElement : public MechanicsElement
{
public:
    enum MaterialUpdateType { RESET, JACOBIAN };
 
    explicit MechanicsNonlinearElement( std::shared_ptr<const ElementOperationParameters> params )
        : MechanicsElement( params ), d_elementOutputVector( nullptr )
    {
        d_JxW = &( d_fe->get_JxW() );
 
        d_dphi = &( d_fe->get_dphi() );
 
        d_phi = &( d_fe->get_phi() );
 
        d_xyz = &( d_fe->get_xyz() );
 
        AMP_INSIST( ( d_useJaumannRate == false ),
                    "Jaumann rate with small strain does not make any sense." );
    }
 
    virtual ~MechanicsNonlinearElement() {}
 
    void setElementVectors( const std::vector<std::vector<double>> &elementInputVectors,
                            std::vector<double> &elementOutputVector )
    {
        d_elementInputVectors = elementInputVectors;
        d_elementOutputVector = &( elementOutputVector );
    }
 
    void apply() override
    {
        if ( d_useReducedIntegration ) {
            apply_Reduced();
        } else {
            apply_Normal();
        }
    }
 
    void initMaterialModel( const std::vector<double> &initTempVector );
 
    void updateMaterialModel( MaterialUpdateType type,
                              const std::vector<std::vector<double>> &elementInputVectors );
 
    void printStressAndStrain( FILE *fp, const std::vector<std::vector<double>> &inputVec );
 
    void computeStressAndStrain( const std::vector<std::vector<double>> &inputVec,
                                 std::vector<double> &stressVec,
                                 std::vector<double> &strainVec );
 
protected:
    void materialModelPreNonlinearElementOperation( MaterialUpdateType );
 
    void materialModelPreNonlinearGaussPointOperation( MaterialUpdateType );
 
    void materialModelNonlinearGaussPointOperation( MaterialUpdateType,
                                                    const std::vector<std::vector<double>> & );
 
    void materialModelPostNonlinearGaussPointOperation( MaterialUpdateType );
 
    void materialModelPostNonlinearElementOperation( MaterialUpdateType );
 
    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<std::vector<libMesh::Real>>
        *d_phi; 
    const std::vector<libMesh::Point>
        *d_xyz; 
    std::vector<std::vector<double>> d_elementInputVectors; 
    std::vector<double> *d_elementOutputVector; 
private:
};
 
 
} // namespace AMP::Operator
 
#endif