suanPan
CDP.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 ******************************************************************************/
45#ifndef CDP_H
46#define CDP_H
47
48#include "NonlinearCDP.h"
49
50class CDP final : public NonlinearCDP {
51 const double a_t, cb_t, f_t;
52 const double a_c, cb_c, f_c;
53
54 [[nodiscard]] podarray<double> compute_tension_backbone(double) const override;
55 [[nodiscard]] podarray<double> compute_compression_backbone(double) const override;
56
57public:
58 explicit CDP(
59 unsigned = 0, // tag
60 double = 3E4, // elastic modulus
61 double = .2, // poissons ratio
62 double = 3., // crack stress (+)
63 double = 30., // crush stress (-)
64 double = 1E-3, // normalized crack energy (+)
65 double = 1E-1, // normalized crush energy (+)
66 double = .8, // hardening after crack stress a_t
67 double = 4., // hardening after crush stress a_c
68 double = .6, // reference damage factor at half crack stress
69 double = .6, // reference damage factor at crush stress
70 double = .2, // dilatancy parameter
71 double = 1.16, // biaxial compression strength ratio
72 double = .5, // stiffness recovery
73 double = 2400E-12 // density
74 );
75
77};
78
79#endif
80
The CDP class.
Definition: CDP.h:50
CDP(unsigned=0, double=3E4, double=.2, double=3., double=30., double=1E-3, double=1E-1, double=.8, double=4., double=.6, double=.6, double=.2, double=1.16, double=.5, double=2400E-12)
Definition: CDP.cpp:56
unique_ptr< Material > get_copy() override
Definition: CDP.cpp:80
The NonlinearCDP class.
Definition: NonlinearCDP.h:54