28namespace LargeDeformation
30template <
int DisplacementDim>
34 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&& jacobian_assembler,
35 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const& parameters,
36 unsigned const integration_order,
37 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
41 :
Process(std::move(name), mesh, std::move(jacobian_assembler), parameters,
42 integration_order, std::move(process_variables),
43 std::move(secondary_variables)),
44 _process_data(std::move(process_data))
50 MeshLib::getOrCreateMeshProperty<double>(
51 const_cast<MeshLib::Mesh&
>(mesh),
"principal_stress_vector_1",
55 MeshLib::getOrCreateMeshProperty<double>(
56 const_cast<MeshLib::Mesh&
>(mesh),
"principal_stress_vector_2",
60 MeshLib::getOrCreateMeshProperty<double>(
61 const_cast<MeshLib::Mesh&
>(mesh),
"principal_stress_vector_3",
65 MeshLib::getOrCreateMeshProperty<double>(
71 DisplacementDim>::getReflectionDataForOutput(),
76template <
int DisplacementDim>
82template <
int DisplacementDim>
86 unsigned const integration_order)
94 auto add_secondary_variable = [&](std::string
const& name,
95 int const num_components,
96 auto get_ip_values_function)
98 _secondary_variables.addSecondaryVariable(
102 std::move(get_ip_values_function)));
105 ProcessLib::Reflection::addReflectedSecondaryVariables<DisplacementDim>(
107 DisplacementDim>::getReflectionDataForOutput(),
108 _secondary_variables, getExtrapolator(), _local_assemblers);
115 add_secondary_variable);
119 _process_data.solid_materials, _local_assemblers,
120 _integration_point_writer, integration_order);
128 *_local_to_global_index_map);
131template <
int DisplacementDim>
133 double const t,
double const dt, std::vector<GlobalVector*>
const& x,
134 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
137 DBUG(
"Assemble LargeDeformationProcess.");
141 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_table = {
142 _local_to_global_index_map.get()};
149 _global_output(t, process_id, M, K, b);
152template <
int DisplacementDim>
155 double const t,
double const dt, std::vector<GlobalVector*>
const& x,
156 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
159 DBUG(
"AssembleWithJacobian LargeDeformationProcess.");
163 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_table = {
164 _local_to_global_index_map.get()};
169 x_prev, process_id, M, K, b, Jac);
171 transformVariableFromGlobalVector(b, 0, *_local_to_global_index_map,
172 *_nodal_forces, std::negate<double>());
174 _global_output(t, process_id, M, K, b, &Jac);
177template <
int DisplacementDim>
179 std::vector<GlobalVector*>
const& x,
180 std::vector<GlobalVector*>
const& x_prev,
double const t,
double const dt,
181 int const process_id)
183 DBUG(
"PostTimestep LargeDeformationProcess.");
184 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
185 dof_tables.reserve(x.size());
186 std::generate_n(std::back_inserter(dof_tables), x.size(),
187 [&]() { return _local_to_global_index_map.get(); });
195template <
int DisplacementDim>
197 double const t,
double const dt, std::vector<GlobalVector*>
const& x,
200 DBUG(
"Compute the secondary variables for LargeDeformationProcess.");
201 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
202 dof_tables.reserve(x.size());
203 std::generate_n(std::back_inserter(dof_tables), x.size(),
204 [&]() { 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)