OGS
SmallDeformationProcess.h
Go to the documentation of this file.
1
11#pragma once
12
14#include "ProcessLib/Process.h"
16
17namespace ProcessLib
18{
19namespace LIE
20{
21namespace SmallDeformation
22{
23template <int DisplacementDim>
24class SmallDeformationProcess final : public Process
25{
26 static_assert(DisplacementDim == 2 || DisplacementDim == 3,
27 "Currently LIE::SmallDeformationProcess "
28 "supports only 2D or 3D.");
29
30public:
32 std::string name,
33 MeshLib::Mesh& mesh,
34 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&&
35 jacobian_assembler,
36 std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const&
37 parameters,
38 unsigned const integration_order,
39 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
40 process_variables,
42 SecondaryVariableCollection&& secondary_variables);
43
46 bool isLinear() const override;
48
49 void computeSecondaryVariableConcrete(double const t, double const dt,
50 std::vector<GlobalVector*> const& x,
51 GlobalVector const& x_prev,
52 int const process_id) override;
53
54private:
56
57 void constructDofTable() override;
58
60 NumLib::LocalToGlobalIndexMap const& dof_table,
61 MeshLib::Mesh const& mesh,
62 unsigned const integration_order) override;
63
64 void assembleConcreteProcess(const double t, double const dt,
65 std::vector<GlobalVector*> const& x,
66 std::vector<GlobalVector*> const& x_prev,
67 int const process_id, GlobalMatrix& M,
68 GlobalMatrix& K, GlobalVector& b) override;
69
71 const double t, double const dt, std::vector<GlobalVector*> const& x,
72 std::vector<GlobalVector*> const& x_prev, int const process_id,
73 GlobalVector& b, GlobalMatrix& Jac) override;
74
75 void preTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
76 double const t, double const dt,
77 const int process_id) override;
78
79private:
81
82 std::vector<std::unique_ptr<LocalAssemblerInterface>> _local_assemblers;
83
84 std::unique_ptr<NumLib::LocalToGlobalIndexMap>
86
87 std::vector<MeshLib::Element*> _vec_matrix_elements;
88 std::vector<int> _vec_fracture_mat_IDs;
89 std::vector<std::vector<MeshLib::Element*>> _vec_fracture_elements;
90 std::vector<std::vector<MeshLib::Element*>> _vec_fracture_matrix_elements;
91 std::vector<std::vector<MeshLib::Node*>> _vec_fracture_nodes;
92 std::vector<MeshLib::Node*> _vec_junction_nodes;
93 std::vector<std::vector<MeshLib::Element*>>
95
96 std::vector<std::unique_ptr<MeshLib::MeshSubset const>>
98 std::unique_ptr<MeshLib::MeshSubset const> _mesh_subset_junction_nodes;
99 std::unique_ptr<MeshLib::MeshSubset const> _mesh_subset_matrix_nodes;
100};
101
102extern template class SmallDeformationProcess<2>;
103extern template class SmallDeformationProcess<3>;
104
105} // namespace SmallDeformation
106} // namespace LIE
107} // namespace ProcessLib
Global vector based on Eigen vector.
Definition EigenVector.h:25
void initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
Process specific initialization called by initialize().
std::vector< std::vector< MeshLib::Element * > > _vec_junction_fracture_matrix_elements
void computeSecondaryVariableConcrete(double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id) override
std::unique_ptr< NumLib::LocalToGlobalIndexMap > _local_to_global_index_map_single_component
std::vector< std::vector< MeshLib::Node * > > _vec_fracture_nodes
std::vector< std::unique_ptr< MeshLib::MeshSubset const > > _mesh_subset_fracture_nodes
void preTimestepConcreteProcess(std::vector< GlobalVector * > const &x, double const t, double const dt, const int process_id) override
void assembleWithJacobianConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalVector &b, GlobalMatrix &Jac) override
std::vector< std::vector< MeshLib::Element * > > _vec_fracture_elements
std::unique_ptr< MeshLib::MeshSubset const > _mesh_subset_junction_nodes
SmallDeformationProcess(std::string name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters, unsigned const integration_order, std::vector< std::vector< std::reference_wrapper< ProcessVariable > > > &&process_variables, SmallDeformationProcessData< DisplacementDim > &&process_data, SecondaryVariableCollection &&secondary_variables)
void assembleConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b) override
std::unique_ptr< MeshLib::MeshSubset const > _mesh_subset_matrix_nodes
std::vector< std::vector< MeshLib::Element * > > _vec_fracture_matrix_elements
SmallDeformationProcessData< DisplacementDim > _process_data
std::vector< std::unique_ptr< LocalAssemblerInterface > > _local_assemblers
std::string const name
Definition Process.h:362
Handles configuration of several secondary variables from the project file.