suanPan
ArmstrongFrederick.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 ******************************************************************************/
36#ifndef ARMSTRONGFREDERICK_H
37#define ARMSTRONGFREDERICK_H
38
40
42 const double elastic_modulus; // elastic modulus
43 const double poissons_ratio; // poisson's ratio
44 const double yield; // yield stress
45 const double saturated;
46 const double hardening;
47 const double m;
48 const vec a, b;
49};
50
51class ArmstrongFrederick final : protected DataArmstrongFrederick, public Material3D {
52 static constexpr unsigned max_iteration = 20u;
53 static const double root_three_two;
54 static const mat unit_dev_tensor;
55
56 const unsigned size = static_cast<unsigned>(a.size());
57
58 const double shear = elastic_modulus / (2. + 2. * poissons_ratio); // shear modulus
59 const double double_shear = 2. * shear;
60 const double three_shear = 3. * shear;
61 const double root_six_shear = sqrt(6.) * shear;
62
63public:
65 unsigned, // tag
66 double, // elastic modulus
67 double, // poissons ratio
68 double, // yield stress
69 double, // saturated stress
70 double, // linear hardening modulus
71 double, // m
72 vec&&, // a
73 vec&&, // b
74 double = 0. // density
75 );
76
77 int initialize(const shared_ptr<DomainBase>&) override;
78
80
81 [[nodiscard]] double get_parameter(ParameterType) const override;
82
83 int update_trial_status(const vec&) override;
84
85 int clear_status() override;
86 int commit_status() override;
87 int reset_status() override;
88
89 void print() override;
90};
91
92#endif
93
ParameterType
Definition: ParameterType.h:21
The ArmstrongFrederick class defines a nonlinear hardening material with mixed hardening (isotropic a...
Definition: ArmstrongFrederick.h:51
int reset_status() override
Definition: ArmstrongFrederick.cpp:137
int initialize(const shared_ptr< DomainBase > &) override
Definition: ArmstrongFrederick.cpp:29
void print() override
Definition: ArmstrongFrederick.cpp:145
unique_ptr< Material > get_copy() override
Definition: ArmstrongFrederick.cpp:37
double get_parameter(ParameterType) const override
Definition: ArmstrongFrederick.cpp:39
int clear_status() override
Definition: ArmstrongFrederick.cpp:121
int commit_status() override
Definition: ArmstrongFrederick.cpp:129
ArmstrongFrederick(unsigned, double, double, double, double, double, double, vec &&, vec &&, double=0.)
Definition: ArmstrongFrederick.cpp:25
int update_trial_status(const vec &) override
Definition: ArmstrongFrederick.cpp:41
The Material3D class.
Definition: Material3D.h:37
Definition: ArmstrongFrederick.h:41
const double saturated
Definition: ArmstrongFrederick.h:45
const double hardening
Definition: ArmstrongFrederick.h:46
const double poissons_ratio
Definition: ArmstrongFrederick.h:43
const double m
Definition: ArmstrongFrederick.h:47
const vec a
Definition: ArmstrongFrederick.h:48
const double yield
Definition: ArmstrongFrederick.h:44
const double elastic_modulus
Definition: ArmstrongFrederick.h:42
const vec b
Definition: ArmstrongFrederick.h:48