Homec4science
Diffusion akantu (features/52-anisotropic-and-at1-phase-field-models) Browse

rAKA/src/model/phase_field/phasefields/phasefield_linear_inline_impl.hh
a4e99c2ef93bfeatures/52-anisotropic…

Last Change
Raw File

phasefield_linear_inline_impl.hh
View Options

#include "phasefield_linear.hh"
namespace akantu {
inline void PhaseFieldLinear::computeDissipatedEnergyOnQuad(
const Real & dam, const Real & dam_prev, const Vector<Real> & grad_d,
Real & edis, Real & g_c_quad) {
edis = 0.;
for (auto i : arange(spatial_dimension)) {
edis += 3. / 8 * g_c_quad * this->l0 * grad_d[i] * grad_d[i];
}
edis += 3 * g_c_quad * dam * dam / (8 * this->l0);
Real damage_penalized = std::min(0., dam - dam_prev);
edis += 0.5 * this->gamma * damage_penalized * damage_penalized;
}
/* -------------------------------------------------------------------------- */
inline void PhaseFieldLinear::computeDrivingForceOnQuad(
const Real & phi_quad, Real & driving_force_quad, const Real & g_c_quad) {
driving_force_quad = phi_quad - 3 * g_c_quad / (16 * this->l0);
}
/* -------------------------------------------------------------------------- */
inline void PhaseFieldLinear::computeDamageEnergyDensityOnQuad(
const Real & phi_quad, Real & dam_energy_quad, const Real & g_c_quad) {
dam_energy_quad = 2 * phi_quad;
}
/* -------------------------------------------------------------------------- */
inline void PhaseFieldLinear::computePhiOnQuad(const Matrix<Real> & strain_quad,
Real & phi_quad) {
Real trace = strain_quad.trace();
Real trace_plus = std::max(Real(0.), trace);
Matrix<Real> strain_dev(dim, dim);
Matrix<Real> strain_tmp = Matrix<Real>::Zero(dim, dim);
strain_tmp.topLeftCorner(spatial_dimension, spatial_dimension) = strain_quad;
strain_dev = strain_tmp -
trace / Real(dim) *
Matrix<Real>::Identity(dim, dim);
Real kpa = this->lambda + 2. * this->mu / Real(dim);
phi_quad = 0.5 * kpa * trace_plus * trace_plus +
this->mu * strain_dev.doubleDot(strain_dev);
}
/* -------------------------------------------------------------------------- */
inline void
PhaseFieldLinear::computePhiIsotropicOnQuad(const Matrix<Real> & strain_quad,
Real & phi_quad) {
Real trace = strain_quad.trace();
phi_quad = 0.5 * this->lambda * trace * trace +
this->mu * strain_quad.doubleDot(strain_quad);
}
} // namespace akantu
c4science · Help