27namespace SmallDeformation
29template <
int DisplacementDim>
33 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&& jacobian_assembler,
34 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const& parameters,
35 unsigned const integration_order,
36 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
40 :
Process(std::move(name), mesh, std::move(jacobian_assembler), parameters,
41 integration_order, std::move(process_variables),
42 std::move(secondary_variables)),
43 _process_data(std::move(process_data))
52 MeshLib::getOrCreateMeshProperty<double>(
53 const_cast<MeshLib::Mesh&
>(mesh),
"principal_stress_vector_1",
57 MeshLib::getOrCreateMeshProperty<double>(
58 const_cast<MeshLib::Mesh&
>(mesh),
"principal_stress_vector_2",
62 MeshLib::getOrCreateMeshProperty<double>(
63 const_cast<MeshLib::Mesh&
>(mesh),
"principal_stress_vector_3",
67 MeshLib::getOrCreateMeshProperty<double>(
73 DisplacementDim>::getReflectionDataForOutput(),
78template <
int DisplacementDim>
84template <
int DisplacementDim>
88 unsigned const integration_order)
96 auto add_secondary_variable = [&](std::string
const& name,
97 int const num_components,
98 auto get_ip_values_function)
100 _secondary_variables.addSecondaryVariable(
104 std::move(get_ip_values_function)));
107 ProcessLib::Reflection::addReflectedSecondaryVariables<DisplacementDim>(
109 DisplacementDim>::getReflectionDataForOutput(),
110 _secondary_variables, getExtrapolator(), _local_assemblers);
117 add_secondary_variable);
121 _process_data.solid_materials, _local_assemblers,
122 _integration_point_writer, integration_order);
130 *_local_to_global_index_map);
133template <
int DisplacementDim>
135 double const t,
double const dt, std::vector<GlobalVector*>
const& x,
136 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
139 DBUG(
"Assemble SmallDeformationProcess.");
143 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_table = {
144 _local_to_global_index_map.get()};
151 _global_output(t, process_id, M, K, b);
154template <
int DisplacementDim>
157 double const t,
double const dt, std::vector<GlobalVector*>
const& x,
158 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
161 DBUG(
"AssembleWithJacobian SmallDeformationProcess.");
165 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_table = {
166 _local_to_global_index_map.get()};
171 x_prev, process_id, M, K, b, Jac);
173 transformVariableFromGlobalVector(b, 0, *_local_to_global_index_map,
174 *_nodal_forces, std::negate<double>());
176 _global_output(t, process_id, M, K, b, &Jac);
179template <
int DisplacementDim>
181 std::vector<GlobalVector*>
const& x,
182 std::vector<GlobalVector*>
const& x_prev,
double const t,
double const dt,
183 int const process_id)
185 DBUG(
"PostTimestep SmallDeformationProcess.");
186 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
187 dof_tables.reserve(x.size());
188 std::generate_n(std::back_inserter(dof_tables), x.size(),
189 [&]() { return _local_to_global_index_map.get(); });
196 std::unique_ptr<GlobalVector> material_forces;
198 material_forces, _local_assemblers, *_local_to_global_index_map,
201 material_forces->copyValues(*_material_forces);
204template <
int DisplacementDim>
206 double const t,
double const dt, std::vector<GlobalVector*>
const& x,
209 DBUG(
"Compute the secondary variables for SmallDeformationProcess.");
210 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
211 dof_tables.reserve(x.size());
212 std::generate_n(std::back_inserter(dof_tables), x.size(),
213 [&]() { return _local_to_global_index_map.get(); });
void DBUG(fmt::format_string< Args... > fmt, Args &&... args)
Global vector based on Eigen vector.
bool isAxiallySymmetric() const
std::vector< Element * > const & getElements() const
Get the element-vector for the mesh.
Properties & getProperties()
virtual void postTimestep(std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
virtual void computeSecondaryVariable(std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id)
virtual void initialize(std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)
std::vector< std::size_t > const & getActiveElementIDs() const
std::vector< std::unique_ptr< MeshLib::IntegrationPointWriter > > _integration_point_writer
Handles configuration of several secondary variables from the project file.
void assembleWithJacobian(std::size_t const mesh_item_id, LocalAssemblerInterface &local_assembler, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, 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, GlobalMatrix &Jac)
void assemble(std::size_t const mesh_item_id, LocalAssemblerInterface &local_assembler, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const 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)
void addReflectedIntegrationPointWriters(ReflData const &reflection_data, std::vector< std::unique_ptr< MeshLib::IntegrationPointWriter > > &integration_point_writers, unsigned const integration_order, std::vector< std::unique_ptr< LocAsmIF > > const &local_assemblers)
SecondaryVariableFunctions makeExtrapolator(const unsigned num_components, NumLib::Extrapolator &extrapolator, LocalAssemblerCollection const &local_assemblers, typename NumLib::ExtrapolatableLocalAssemblerCollection< LocalAssemblerCollection >::IntegrationPointValuesMethod integration_point_values_method)
void setIPDataInitialConditions(std::vector< std::unique_ptr< MeshLib::IntegrationPointWriter > > const &_integration_point_writer, MeshLib::Properties const &mesh_properties, LocalAssemblersVector &local_assemblers)
static void executeSelectedMemberOnDereferenced(Method method, Container const &container, std::vector< std::size_t > const &active_container_ids, Args &&... args)
static void executeSelectedMemberDereferenced(Object &object, Method method, Container const &container, std::vector< std::size_t > const &active_container_ids, Args &&... args)
static void executeMemberOnDereferenced(Method method, Container const &container, Args &&... args)