suanPan
DafaliasManzari.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 DAFALIASMANZARI_H
30#define DAFALIASMANZARI_H
31
33
35 const double shear_modulus = 125.; // shear modulus
36 const double poissons_ratio = .05; // poisson's ratio
37 const double ac = 1.25;
38 const double lc = .02;
39 const double e0 = .9;
40 const double xi = .7;
41 const double m = .01;
42 const double h0 = 7.;
43 const double h1 = .1;
44 const double ch = .9;
45 const double nb = 1.1;
46 const double a = -.7;
47 const double nd = 3.5;
48 const double zm = 4.;
49 const double cz = 6E2;
50 const double pc = -130.;
51 const double gr = .2;
52};
53
54class DafaliasManzari final : protected DataDafaliasManzari, public Material3D {
55 static constexpr unsigned max_iteration = 20u;
56 static const mat unit_dev_tensor;
57
58 const double pr = (2. + 2. * poissons_ratio) / (3. - 6. * poissons_ratio);
59 const double gi = gr * shear_modulus * abs(pc);
60
61 static constexpr uword sa = 0, si = 0, sj = 1;
62 static const span sb, sk, sl, sm;
63
64public:
66 unsigned, // tag
67 double, // g0
68 double, // nu
69 double, // ac
70 double, // lc
71 double, // e0
72 double, // xi
73 double, // m
74 double, // h0
75 double, // h1
76 double, // ch
77 double, // nb
78 double, // a
79 double, // nd
80 double, // zm
81 double, // cz
82 double, // pc
83 double, // gr
84 double = 0. // density
85 );
86
87 int initialize(const shared_ptr<DomainBase>&) override;
88
90
91 [[nodiscard]] double get_parameter(ParameterType) const override;
92
93 int update_trial_status(const vec&) override;
94
95 int clear_status() override;
96 int commit_status() override;
97 int reset_status() override;
98
99 void print() override;
100};
101
102#endif
103
ParameterType
Definition: ParameterType.h:21
The DafaliasManzari class.
Definition: DafaliasManzari.h:54
void print() override
Definition: DafaliasManzari.cpp:379
int update_trial_status(const vec &) override
Definition: DafaliasManzari.cpp:45
unique_ptr< Material > get_copy() override
Definition: DafaliasManzari.cpp:41
DafaliasManzari(unsigned, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double=0.)
Definition: DafaliasManzari.cpp:29
int clear_status() override
Definition: DafaliasManzari.cpp:355
int reset_status() override
Definition: DafaliasManzari.cpp:371
double get_parameter(ParameterType) const override
Definition: DafaliasManzari.cpp:43
int commit_status() override
Definition: DafaliasManzari.cpp:363
int initialize(const shared_ptr< DomainBase > &) override
Definition: DafaliasManzari.cpp:33
The Material3D class.
Definition: Material3D.h:37
Definition: DafaliasManzari.h:34
const double ch
Definition: DafaliasManzari.h:44
const double e0
Definition: DafaliasManzari.h:39
const double poissons_ratio
Definition: DafaliasManzari.h:36
const double a
Definition: DafaliasManzari.h:46
const double shear_modulus
Definition: DafaliasManzari.h:35
const double m
Definition: DafaliasManzari.h:41
const double h0
Definition: DafaliasManzari.h:42
const double nd
Definition: DafaliasManzari.h:47
const double pc
Definition: DafaliasManzari.h:50
const double h1
Definition: DafaliasManzari.h:43
const double zm
Definition: DafaliasManzari.h:48
const double gr
Definition: DafaliasManzari.h:51
const double ac
Definition: DafaliasManzari.h:37
const double nb
Definition: DafaliasManzari.h:45
const double xi
Definition: DafaliasManzari.h:40
const double cz
Definition: DafaliasManzari.h:49
const double lc
Definition: DafaliasManzari.h:38