From ccbff6fd0c76ec2a3ba422f2cff02330a9acebf8 Mon Sep 17 00:00:00 2001 From: Prithivirajan Veerappan Date: Thu, 23 Feb 2023 15:19:09 -0700 Subject: [PATCH] Updated comments per review comments, and modified dilatancy factor name in the tests ref# 184 --- doc/content/bib/blackbear.bib | 4 ++-- .../materials/ComputeMultipleInelasticDamageStress.md | 8 ++++++-- src/materials/DamagePlasticityStressUpdate.C | 2 +- test/tests/damage_plasticity_model/tests | 2 +- 4 files changed, 10 insertions(+), 6 deletions(-) diff --git a/doc/content/bib/blackbear.bib b/doc/content/bib/blackbear.bib index f5fa6ac0..cb2523db 100644 --- a/doc/content/bib/blackbear.bib +++ b/doc/content/bib/blackbear.bib @@ -256,11 +256,11 @@ @article{lee1998plastic publisher={American Society of Civil Engineers} } -@book{lee1996theory, +@phdthesis{lee1996theory, title={Theory and implementation of plastic-damage model for concrete structures under cyclic and dynamic loading}, author={Lee, Jeeho}, year={1996}, - publisher={University of California, Berkeley} + school={University of California, Berkeley} } @article{wilkins2020method, diff --git a/doc/content/source/materials/ComputeMultipleInelasticDamageStress.md b/doc/content/source/materials/ComputeMultipleInelasticDamageStress.md index 6fef69c6..6a8cb3ec 100644 --- a/doc/content/source/materials/ComputeMultipleInelasticDamageStress.md +++ b/doc/content/source/materials/ComputeMultipleInelasticDamageStress.md @@ -4,8 +4,12 @@ ## Description -This class `ComputeMultipleInelasticDamageStress` computes the stress with the damage obtained from - the previous time step. This is done mainly to avoid convergence problems +`ComputeMultipleInelasticDamageStress` is a class based on `ComputeMultipleInelasticStress` +that is designed specifically to be used with [DamagePlasticityStressUpdate](DamagePlasticityStressUpdate.md). +The main difference between this class and the standard `ComputeMultipleInelasticStress` class is that +it makes some modifications to both the current and old stress to account for the damage in a way that +is appropriate for the damaged plasticity model. It divides the old stress passed to the plasticity model +by $(1-D)$, and multiplies the new stress by $(1-D)$, where $D$ is the damage index. This ComputeMultipleInelasticStress is to be used with (/DamagePlasticityStressUpdate.md). ## Example Input Files diff --git a/src/materials/DamagePlasticityStressUpdate.C b/src/materials/DamagePlasticityStressUpdate.C index bd06e05a..d9006449 100644 --- a/src/materials/DamagePlasticityStressUpdate.C +++ b/src/materials/DamagePlasticityStressUpdate.C @@ -538,7 +538,7 @@ DamagePlasticityStressUpdate::ft(const std::vector & intnl) const if (_zt > sqrtPhi_t / _at) return _ft0 * std::pow(_zt - sqrtPhi_t / _at, (1. - _dt_bt)) * sqrtPhi_t; else - return _ft0 * 1.E-6; + return _ft0 * 1.E-6; // A very small number (instead of zero) is used for end of softening } Real diff --git a/test/tests/damage_plasticity_model/tests b/test/tests/damage_plasticity_model/tests index 5ac05b4a..3a9d23e2 100644 --- a/test/tests/damage_plasticity_model/tests +++ b/test/tests/damage_plasticity_model/tests @@ -41,7 +41,7 @@ cli_args = 'BCs/x_r/function=-2E-5*x*t Executioner/end_time=500 Executioner/dt=10 - Materials/damage_plasticity_model/factor_controlling_dilatancy=0.30 + Materials/damage_plasticity_model/dilatancy_factor=0.30 Outputs/file_base=dilatancy_out' detail = ' when factor_controlling_dilatancy is changed.' []