37 static constexpr unsigned max_iteration = 20u;
38 static constexpr double min_ratio = 1.;
40 static double dev(
const vec&);
41 static rowvec3 der_dev(
const vec&);
42 static mat compute_stiffness(
double,
double);
44 using ds_moduli = std::tuple<double, double, rowvec3, rowvec3>;
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;
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();
56 int project_onto_surface(
double&);
57 int local_update(
const vec&,
const vec&,
bool);
66 int initialize(
const shared_ptr<DomainBase>&)
override;
76 void print()
override;
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