From 02cb60e71e73cc56cfd91c01bd0dddc979b9f674 Mon Sep 17 00:00:00 2001 From: nagelt Date: Sat, 19 Oct 2024 14:53:56 +0200 Subject: [PATCH] MC pc_char to avoid non-zero stiffness --- MaterialLib/SolidModels/MFront/ModCamClay_semiExpl.mfront | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/MaterialLib/SolidModels/MFront/ModCamClay_semiExpl.mfront b/MaterialLib/SolidModels/MFront/ModCamClay_semiExpl.mfront index 7dc9595e305..7d5224901ba 100644 --- a/MaterialLib/SolidModels/MFront/ModCamClay_semiExpl.mfront +++ b/MaterialLib/SolidModels/MFront/ModCamClay_semiExpl.mfront @@ -28,7 +28,7 @@ @Theta 1.0; // time integration scheme @Epsilon 1e-14; // tolerance of local stress integration algorithm -@MaximumNumberOfIterations 20; // for local local stress integration algorithm +@MaximumNumberOfIterations 100; // for local local stress integration algorithm @ModellingHypotheses{".+"}; // supporting all stress and strain states @Algorithm NewtonRaphson; //_NumericalJacobian;_LevenbergMarquardt @@ -95,13 +95,13 @@ v.setEntryName("VolumeRatio"); // Total volume per solid volume = inv(1 - poros sig0 = sig; // compute elastic stiffness (constant during time step) + pc_min = 0.5e-8 * pc_char; const auto p = -trace(sig) / 3; - const auto K = v0 / ka * p; + const auto K = v0 / ka * std::max(p,pc_char); const auto E = 3.0 * K * (1.0 - 2*nu); young = E; rpc = pc / young; - pc_min = 0.5e-8 * pc_char; rpc_min = pc_min / young; // stress derivative