ReleaseNodalForce.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#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