suanPan
Mindlin.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 ******************************************************************************/
30#ifndef MINDLIN_H
31#define MINDLIN_H
32
34
35class Mindlin final : public MaterialElement2D {
36 struct IntegrationPoint final {
38 const double eccentricity, factor;
42 };
43
44 vec coor;
45 mat strain_mat;
46 vector<SectionIntegrationPoint> sec_int_pt;
47 explicit IntegrationPoint(vec&&);
48 };
49
50 static constexpr unsigned p_node = 4, p_dof = 3, p_size = p_dof * p_node;
51
52 const double thickness;
53 const unsigned num_section_ip;
54
55 vector<IntegrationPoint> int_pt;
56
57 mat penalty_stiffness;
58
59public:
60 Mindlin(
61 unsigned, // element tag
62 uvec&&, // node tag
63 unsigned, // material tag
64 double, // thickness
65 unsigned = 2 // integration points along thickness
66 );
67
68 int initialize(const shared_ptr<DomainBase>&) override;
69
70 int update_status() override;
71
72 int commit_status() override;
73 int clear_status() override;
74 int reset_status() override;
75
76 vector<vec> record(OutputType) override;
77
78 void print() override;
79
80#ifdef SUANPAN_VTK
81 void Setup() override;
82 void GetData(vtkSmartPointer<vtkDoubleArray>&, OutputType) override;
83 void SetDeformation(vtkSmartPointer<vtkPoints>&, double) override;
84#endif
85};
86
87#endif
88
OutputType
Definition: OutputType.h:23
Definition: MaterialElement.h:64
A Mindlin plate class.
Definition: Mindlin.h:35
int clear_status() override
Definition: Mindlin.cpp:128
vector< vec > record(OutputType) override
Definition: Mindlin.cpp:140
int update_status() override
Definition: Mindlin.cpp:106
Mindlin(unsigned, uvec &&, unsigned, double, unsigned=2)
Definition: Mindlin.cpp:40
void print() override
Definition: Mindlin.cpp:146
int initialize(const shared_ptr< DomainBase > &) override
Definition: Mindlin.cpp:45
int reset_status() override
Definition: Mindlin.cpp:134
int commit_status() override
Definition: Mindlin.cpp:122
const double factor
Definition: Mindlin.h:38
const double eccentricity
Definition: Mindlin.h:38
unique_ptr< Material > p_material
Definition: Mindlin.h:39
SectionIntegrationPoint(double, double, unique_ptr< Material > &&)
Definition: Mindlin.cpp:26