suanPan
F31.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 ******************************************************************************/
34#ifndef F31_H
35#define F31_H
36
39
40class F31 final : public SectionElement3D {
41 struct IntegrationPoint final {
42 double coor, weight;
43 unique_ptr<Section> b_section;
44 mat strain_mat; // axial strain, curvature about the z-axis (major), curvature about the y-axis (minor)
45 IntegrationPoint(double, double, unique_ptr<Section>&&);
46 };
47
48 static constexpr unsigned b_node = 2, b_dof = 6, b_size = b_dof * b_node;
49
50 static const span b_span;
51
52 const unsigned int_pt_num, orientation_tag;
53
54 const double length = 0., torsion_stiff = 0.;
55
56 vector<IntegrationPoint> int_pt;
57
58 unique_ptr<Orientation> b_trans;
59
60 mat initial_local_flexibility;
61 mat current_local_flexibility, trial_local_flexibility;
62 vec current_local_deformation, trial_local_deformation;
63 vec current_local_resistance, trial_local_resistance;
64
65public:
66 F31(
67 unsigned, // tag
68 uvec&&, // node tag
69 unsigned, // section tag
70 unsigned, // orientation tag
71 unsigned = 6, // integration points
72 bool = false // nonlinear geometry switch
73 );
74
75 int initialize(const shared_ptr<DomainBase>&) override;
76
77 int update_status() override;
78
79 int clear_status() override;
80 int commit_status() override;
81 int reset_status() override;
82
83 vector<vec> record(OutputType) override;
84
85 void print() override;
86
87#ifdef SUANPAN_VTK
88 void Setup() override;
89 void GetData(vtkSmartPointer<vtkDoubleArray>&, OutputType) override;
90 void SetDeformation(vtkSmartPointer<vtkPoints>&, double) override;
91#endif
92};
93
94#endif
95
OutputType
Definition: OutputType.h:23
The F31 class.
Definition: F31.h:40
int clear_status() override
Definition: F31.cpp:126
void print() override
Definition: F31.cpp:168
int update_status() override
Definition: F31.cpp:91
int commit_status() override
Definition: F31.cpp:137
F31(unsigned, uvec &&, unsigned, unsigned, unsigned=6, bool=false)
Definition: F31.cpp:36
vector< vec > record(OutputType) override
Definition: F31.cpp:159
int initialize(const shared_ptr< DomainBase > &) override
Definition: F31.cpp:41
int reset_status() override
Definition: F31.cpp:148
Definition: SectionElement.h:77