suanPan
NonlinearGurson.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 NONLINEARGURSON_H
30#define NONLINEARGURSON_H
31
33
35 const double elastic_modulus; // elastic modulus
36 const double poissons_ratio; // poisson's ratio
37
38 const double q1 = 1., q2 = 1., fn = 0., sn = 1., en = 0.;
39};
40
42 static constexpr unsigned max_iteration = 20u;
43 static const double sqrt_three_two;
44 static const mat unit_dev_tensor;
45
46 const double six_shear = 3. * elastic_modulus / (1. + poissons_ratio); // six shear modulus
47 const double bulk = elastic_modulus / (3. - 6. * poissons_ratio); // bulk modulus
48
49 const double para_a = 3. * bulk * q1 * q2;
50 const double para_b = fn / sn / sqrt(2. * datum::pi);
51
52 [[nodiscard]] virtual vec compute_hardening(double) const = 0;
53
54public:
56 unsigned, // tag
57 double, // elastic modulus
58 double, // poisson's ratio
59 double, // q1
60 double, // q2
61 double, // fn
62 double, // sn
63 double, // en
64 double = 0. // density
65 );
66
67 int initialize(const shared_ptr<DomainBase>&) override;
68
69 [[nodiscard]] double get_parameter(ParameterType) const override;
70
71 int update_trial_status(const vec&) override;
72
73 int clear_status() override;
74 int commit_status() override;
75 int reset_status() override;
76
77 vector<vec> record(OutputType) override;
78
79 void print() override;
80};
81
82#endif
83
OutputType
Definition: OutputType.h:23
ParameterType
Definition: ParameterType.h:21
The Material3D class.
Definition: Material3D.h:37
The NonlinearGurson class.
Definition: NonlinearGurson.h:41
int reset_status() override
Definition: NonlinearGurson.cpp:155
void print() override
Definition: NonlinearGurson.cpp:169
NonlinearGurson(unsigned, double, double, double, double, double, double, double, double=0.)
Definition: NonlinearGurson.cpp:25
int clear_status() override
Definition: NonlinearGurson.cpp:139
int update_trial_status(const vec &) override
Definition: NonlinearGurson.cpp:39
int commit_status() override
Definition: NonlinearGurson.cpp:147
int initialize(const shared_ptr< DomainBase > &) override
Definition: NonlinearGurson.cpp:29
double get_parameter(ParameterType) const override
Definition: NonlinearGurson.cpp:37
vector< vec > record(OutputType) override
Definition: NonlinearGurson.cpp:163
Definition: NonlinearGurson.h:34
const double en
Definition: NonlinearGurson.h:38
const double q1
Definition: NonlinearGurson.h:38
const double poissons_ratio
Definition: NonlinearGurson.h:36
const double fn
Definition: NonlinearGurson.h:38
const double sn
Definition: NonlinearGurson.h:38
const double q2
Definition: NonlinearGurson.h:38
const double elastic_modulus
Definition: NonlinearGurson.h:35