suanPan
LeeNewmarkBase.h
Go to the documentation of this file.
1/*******************************************************************************
2 * Copyright (C) 2017-2024 Theodore Chang
3 *
4 * This program is free software: you can redistribute it and/or modify
5 * it under the terms of the GNU General Public License as published by
6 * the Free Software Foundation, either version 3 of the License, or
7 * (at your option) any later version.
8 *
9 * This program is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 * GNU General Public License for more details.
13 *
14 * You should have received a copy of the GNU General Public License
15 * along with this program. If not, see <http://www.gnu.org/licenses/>.
16 ******************************************************************************/
29#ifndef LEENEWMARKBASE_H
30#define LEENEWMARKBASE_H
31
32#include <Domain/Factory.hpp>
35
36class LeeNewmarkBase : public Newmark {
37public:
38 enum class StiffnessType {
39 INITIAL,
40 CURRENT,
41 TRIAL
42 };
43
44protected:
45 const uword n_block;
46
48
49 bool first_iteration = true;
50
51 bool if_iterative = false;
52
54
55 mutable vec residual;
56
57 unique_ptr<MetaMat<double>> stiffness = nullptr;
58
59 shared_ptr<Factory<double>> factory = nullptr;
60
61 [[nodiscard]] virtual uword get_total_size() const = 0;
62
63 virtual void update_stiffness() const = 0;
64 virtual void update_residual() const = 0;
65
66 int erase_top_left_block() const;
67
68public:
69 explicit LeeNewmarkBase(unsigned, double, double, StiffnessType = StiffnessType::CURRENT);
70
71 int initialize() override;
72
73 int update_internal(const mat&) final;
74
75 int solve(mat&, const mat&) final;
76 int solve(mat&, const sp_mat&) final;
77 int solve(mat&, mat&&) final;
78 int solve(mat&, sp_mat&&) final;
79
80 vec get_force_residual() final;
81 vec get_displacement_residual() final;
82
83 void commit_status() final;
84 void clear_status() final;
85 void reset_status() final;
86};
87
88#endif
89
A LeeNewmarkBase class defines a solver using Newmark algorithm with Lee damping model.
Definition: LeeNewmarkBase.h:36
bool first_iteration
Definition: LeeNewmarkBase.h:49
const uword n_block
Definition: LeeNewmarkBase.h:45
unique_ptr< MetaMat< double > > stiffness
Definition: LeeNewmarkBase.h:57
vec current_internal
Definition: LeeNewmarkBase.h:53
virtual void update_stiffness() const =0
LeeNewmarkBase(unsigned, double, double, StiffnessType=StiffnessType::CURRENT)
Definition: LeeNewmarkBase.cpp:52
vec get_displacement_residual() final
Definition: LeeNewmarkBase.cpp:110
void clear_status() final
Definition: LeeNewmarkBase.cpp:124
StiffnessType
Definition: LeeNewmarkBase.h:38
bool if_iterative
Definition: LeeNewmarkBase.h:51
int update_internal(const mat &) final
Definition: LeeNewmarkBase.cpp:84
vec residual
Definition: LeeNewmarkBase.h:55
int initialize() override
Definition: LeeNewmarkBase.cpp:57
void commit_status() final
Definition: LeeNewmarkBase.cpp:116
int erase_top_left_block() const
Definition: LeeNewmarkBase.cpp:22
const StiffnessType stiffness_type
Definition: LeeNewmarkBase.h:47
virtual uword get_total_size() const =0
int solve(mat &, const mat &) final
Definition: LeeNewmarkBase.cpp:90
shared_ptr< Factory< double > > factory
Definition: LeeNewmarkBase.h:59
vec get_force_residual() final
Definition: LeeNewmarkBase.cpp:104
void reset_status() final
Definition: LeeNewmarkBase.cpp:132
virtual void update_residual() const =0
vec trial_internal
Definition: LeeNewmarkBase.h:53
A Newmark class defines a solver using Newmark algorithm.
Definition: Newmark.h:45