suanPan
NonlinearCamClay.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 ******************************************************************************/
40#ifndef NONLINEARCAMCLAY_H
41#define NONLINEARCAMCLAY_H
42
44
46 const double elastic_modulus; // elastic modulus
47 const double poissons_ratio; // poisson's ratio
48 const double square_beta; // beta squared
49 const double m; // radius ratio
50 const double pt; // tensile yield hydrostatic stress
51};
52
54 static constexpr unsigned max_iteration = 20u;
55 static const double sqrt_three_two;
56 static const mat unit_dev_tensor;
57
58 const double shear = elastic_modulus / (2. + 2. * poissons_ratio);
59 const double bulk = elastic_modulus / (3. - 6. * poissons_ratio);
60 const double six_shear = 3. * elastic_modulus / (1. + poissons_ratio);
61 const double square_m = m * m;
62
63 [[nodiscard]] virtual double compute_a(double) const = 0;
64 [[nodiscard]] virtual double compute_da(double) const = 0;
65
66public:
68 unsigned, // tag
69 double, // elastic modulus
70 double, // poisson's ratio
71 double, // beta
72 double, // m
73 double, // pt
74 double = 0. // density
75 );
76
77 int initialize(const shared_ptr<DomainBase>&) override;
78
79 [[nodiscard]] double get_parameter(ParameterType) const override;
80
81 int update_trial_status(const vec&) override;
82
83 int clear_status() override;
84 int commit_status() override;
85 int reset_status() override;
86};
87
88#endif
89
ParameterType
Definition: ParameterType.h:21
The Material3D class.
Definition: Material3D.h:37
The NonlinearCamClay class.
Definition: NonlinearCamClay.h:53
double get_parameter(ParameterType) const override
Definition: NonlinearCamClay.cpp:36
int clear_status() override
Definition: NonlinearCamClay.cpp:114
NonlinearCamClay(unsigned, double, double, double, double, double, double=0.)
Definition: NonlinearCamClay.cpp:24
int reset_status() override
Definition: NonlinearCamClay.cpp:130
int initialize(const shared_ptr< DomainBase > &) override
Definition: NonlinearCamClay.cpp:28
int update_trial_status(const vec &) override
Definition: NonlinearCamClay.cpp:38
int commit_status() override
Definition: NonlinearCamClay.cpp:122
Definition: NonlinearCamClay.h:45
const double m
Definition: NonlinearCamClay.h:49
const double pt
Definition: NonlinearCamClay.h:50
const double elastic_modulus
Definition: NonlinearCamClay.h:46
const double square_beta
Definition: NonlinearCamClay.h:48
const double poissons_ratio
Definition: NonlinearCamClay.h:47