suanPan
DKT4.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 DKT4_H
31
#define DKT4_H
32
33
#include <
Element/MaterialElement.h
>
34
35
class
DKT4
final :
public
MaterialElement2D
{
36
struct
IntegrationPoint final {
37
struct
SectionIntegrationPoint
final {
38
const
double
eccentricity
,
factor
;
39
unique_ptr<Material>
p_material
;
40
SectionIntegrationPoint
(
double
,
double
,
unique_ptr<Material>
&&);
41
SectionIntegrationPoint
(
const
SectionIntegrationPoint
&);
42
SectionIntegrationPoint
(
SectionIntegrationPoint
&&) noexcept = default;
43
SectionIntegrationPoint
& operator=(const
SectionIntegrationPoint
&) = delete;
44
SectionIntegrationPoint
& operator=(
SectionIntegrationPoint
&&) noexcept = delete;
45
~
SectionIntegrationPoint
() = default;
46
};
47
48
vec coor;
49
mat strain_mat;
50
vector
<
SectionIntegrationPoint
> sec_int_pt;
51
explicit IntegrationPoint(vec&&);
52
};
53
54
static constexpr
unsigned
p_node = 4, p_dof = 3, p_size = p_dof * p_node;
55
56
const
double
thickness;
57
const
unsigned
num_section_ip;
58
59
vector
<IntegrationPoint> int_pt;
60
61
static field<mat> form_transform(const mat&);
62
63
public:
64
DKT4
(
65
unsigned
,
// element tag
66
uvec&&,
// node tag
67
unsigned
,
// material tag
68
double
,
// thickness
69
unsigned
= 3
// integration points along thickness
70
);
71
72
int
initialize
(const shared_ptr<
DomainBase
>&) override;
73
74
int
update_status
() override;
75
int
clear_status
() override;
76
int
commit_status
() override;
77
int
reset_status
() override;
78
79
vector
<vec>
record
(
OutputType
) override;
80
81
void
print
() override;
82
83
#ifdef SUANPAN_VTK
84
void
Setup()
override
;
85
void
GetData(vtkSmartPointer<vtkDoubleArray>&,
OutputType
)
override
;
86
void
SetDeformation(vtkSmartPointer<vtkPoints>&,
double
)
override
;
87
#endif
88
};
89
90
#endif
91
MaterialElement.h
OutputType
OutputType
Definition:
OutputType.h:23
DKT4
A DKT4 plate class.
Definition:
DKT4.h:35
DKT4::update_status
int update_status() override
Definition:
DKT4.cpp:153
DKT4::record
vector< vec > record(OutputType) override
Definition:
DKT4.cpp:188
DKT4::print
void print() override
Definition:
DKT4.cpp:194
DKT4::clear_status
int clear_status() override
Definition:
DKT4.cpp:170
DKT4::reset_status
int reset_status() override
Definition:
DKT4.cpp:182
DKT4::initialize
int initialize(const shared_ptr< DomainBase > &) override
Definition:
DKT4.cpp:99
DKT4::commit_status
int commit_status() override
Definition:
DKT4.cpp:176
DomainBase
The DomainBase class is a template.
Definition:
DomainBase.h:104
MaterialElement2D
Definition:
MaterialElement.h:64
unique_ptr< Material >
suanpan::vector
std::vector< T > vector
Definition:
container.h:53
DKT4::IntegrationPoint::SectionIntegrationPoint
Definition:
DKT4.h:37
DKT4::IntegrationPoint::SectionIntegrationPoint::factor
const double factor
Definition:
DKT4.h:38
DKT4::IntegrationPoint::SectionIntegrationPoint::SectionIntegrationPoint
SectionIntegrationPoint(double, double, unique_ptr< Material > &&)
Definition:
DKT4.cpp:24
DKT4::IntegrationPoint::SectionIntegrationPoint::p_material
unique_ptr< Material > p_material
Definition:
DKT4.h:39
DKT4::IntegrationPoint::SectionIntegrationPoint::SectionIntegrationPoint
SectionIntegrationPoint(SectionIntegrationPoint &&) noexcept=default
DKT4::IntegrationPoint::SectionIntegrationPoint::eccentricity
const double eccentricity
Definition:
DKT4.h:38
Element
Plate
DKT4.h
Generated by
1.9.4