ReleaseNodalForce.h

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
 
#include "MathLib/LinAlg/GlobalMatrixVectorTypes.h"
#include "NumLib/DOF/LocalToGlobalIndexMap.h"
#include "ProcessLib/BoundaryConditionAndSourceTerm/BoundaryCondition.h"
 
namespace ParameterLib
{
template <typename T>
struct Parameter;
}  // namespace ParameterLib
 
namespace MeshLib
{
class Mesh;
}
 
namespace ProcessLib
{
 
class ReleaseNodalForce final : public BoundaryCondition
{
public:
    explicit ReleaseNodalForce(
        int const variable_id,
        MeshLib::Mesh const& boundary_mesh,
        std::unique_ptr<NumLib::LocalToGlobalIndexMap>& dof_table,
        ParameterLib::Parameter<double> const& time_decay_parameter);
 
    void set(GlobalVector const* r_neq);
 
    void applyNaturalBC(const double t, std::vector<GlobalVector*> const& x,
                        int const process_id, GlobalMatrix* K, GlobalVector& b,
                        GlobalMatrix* Jac) override;
 
private:
    int const variable_id_;
 
    MeshLib::Mesh const& boundary_mesh_;
 
    std::unique_ptr<NumLib::LocalToGlobalIndexMap> const dof_table_;
 
    ParameterLib::Parameter<double> const& time_decay_parameter_;
 
    std::vector<double> initial_release_nodal_force_;
 
    std::vector<GlobalIndexType> global_indices_;
    std::vector<MeshLib::Node const*> boundary_nodes_;
};
 
}  // namespace ProcessLib