suanPan
DuncanSelig.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 ******************************************************************************/
31#ifndef DUNCANSELIG_H
32#define DUNCANSELIG_H
33
35
36class DuncanSelig final : public Material2D {
37 static constexpr unsigned max_iteration = 20u;
38 static constexpr double min_ratio = 1.;
39
40 static double dev(const vec&);
41 static rowvec3 der_dev(const vec&);
42 static mat compute_stiffness(double, double);
43
44 using ds_moduli = std::tuple<double, double, rowvec3, rowvec3>;
45
46 double p_atm = 14.7;
47 double ref_elastic = 400. * p_atm, n = .6;
48 double ref_bulk = 300. * p_atm, m = .2;
49 double ini_phi = .7, ten_fold_phi_diff = .1, r_f = .7, cohesion = .5;
50
51 [[nodiscard]] std::tuple<double, double> compute_elastic(double) const;
52 [[nodiscard]] std::tuple<double, double> compute_bulk(double) const;
53 [[nodiscard]] ds_moduli compute_elastic_moduli();
54 [[nodiscard]] ds_moduli compute_plastic_moduli();
55
56 int project_onto_surface(double&);
57 int local_update(const vec&, const vec&, bool);
58
59public:
61 unsigned, // tag
62 const vec&, // parameters
63 double = 0. // density
64 );
65
66 int initialize(const shared_ptr<DomainBase>&) override;
67
69
70 int update_trial_status(const vec&) override;
71
72 int clear_status() override;
73 int commit_status() override;
74 int reset_status() override;
75
76 void print() override;
77};
78
79#endif
80
A DuncanSelig material class.
Definition: DuncanSelig.h:36
int clear_status() override
Definition: DuncanSelig.cpp:322
DuncanSelig(unsigned, const vec &, double=0.)
Definition: DuncanSelig.cpp:265
int commit_status() override
Definition: DuncanSelig.cpp:330
int update_trial_status(const vec &) override
Definition: DuncanSelig.cpp:289
void print() override
Definition: DuncanSelig.cpp:346
unique_ptr< Material > get_copy() override
Definition: DuncanSelig.cpp:287
int initialize(const shared_ptr< DomainBase > &) override
Definition: DuncanSelig.cpp:277
int reset_status() override
Definition: DuncanSelig.cpp:338
The Material2D class defines a isotropic elastic material for plane stress and plane strain problems.
Definition: Material2D.h:43