From a78b935faf7aa6b82b410c580de21fc5dd08778f Mon Sep 17 00:00:00 2001 From: Antoine Jacquey Date: Wed, 13 Mar 2019 17:02:32 +0100 Subject: [PATCH] Removed initializeDamageParameters as it is not needed anymore --- include/materials/LynxDamageDeformation.h | 1 - src/materials/LynxDamageDeformation.C | 24 ----------------------- 2 files changed, 25 deletions(-) diff --git a/include/materials/LynxDamageDeformation.h b/include/materials/LynxDamageDeformation.h index 7d446fd..16a884e 100644 --- a/include/materials/LynxDamageDeformation.h +++ b/include/materials/LynxDamageDeformation.h @@ -46,7 +46,6 @@ protected: virtual Real convexPlasticIncrement(Real & vol_plastic_incr, Real & eqv_plastic_incr); virtual void computeDamageProperties(const Real & pressure, const Real & eqv_stress); virtual void updateDamageParameters(); - virtual void initializeDamageParameters(); virtual void updateDamageConvexParameters(const Real & pressure, const Real & eqv_stress); virtual Real convexReferencePressure(); virtual Real dConvexPlasticYield2(const Real & rho); diff --git a/src/materials/LynxDamageDeformation.C b/src/materials/LynxDamageDeformation.C index 77d434d..1a5d6de 100644 --- a/src/materials/LynxDamageDeformation.C +++ b/src/materials/LynxDamageDeformation.C @@ -112,9 +112,6 @@ LynxDamageDeformation::LynxDamageDeformation(const InputParameters & parameters) void LynxDamageDeformation::initializeQpDeformation() { - // Need to get _xi0 - initializeDamageParameters(); - // Initialize yield derivative _dyield_dp_tr = 0.0; _dyield_dq_tr = 0.0; @@ -415,27 +412,6 @@ LynxDamageDeformation::updateDamageParameters() averageProperty(_one_on_damage_eta)); } -void -LynxDamageDeformation::initializeDamageParameters() -{ - Real sin_phi = std::sin(averageProperty(_friction_angle) * libMesh::pi / 180.0); - Real xi0 = -std::sqrt(3.0) / - std::sqrt(1.0 + 1.5 * Utility::pow<2>(averageProperty(_bulk_modulus) / - averageProperty(_shear_modulus) * sin_phi)); - Real porous_coeff = - averageProperty(_porous_coeff) + averageProperty(_porous_coeff_linear) * _porosity[_qp]; - Real p_cr = _K[_qp] * (std::sqrt(3.0) + xi0) / (3.0 * porous_coeff); - Real k0 = std::sqrt(3.0) * averageProperty(_cohesion) * - std::cos(averageProperty(_friction_angle) * libMesh::pi / 180.0); - - _damage_plasticity->fill(xi0, - averageProperty(_damage_modulus), - p_cr, - k0, - averageProperty(_one_on_plastic_eta), - averageProperty(_one_on_damage_eta)); -} - void LynxDamageDeformation::updateDamageConvexParameters(const Real & pressure, const Real & eqv_stress) { -- GitLab