From 974113c9b870f2aa5f80f543a3c6c12e2437d898 Mon Sep 17 00:00:00 2001 From: FLE_Knight <ujwalkishore.jinaga@mail.polimi.it> Date: Fri, 12 May 2023 12:59:55 +0200 Subject: [PATCH] commented the numerical Derivatives for mechSrc --- NonLinearSolver/materialLaw/mlawNonLinearTVP.cpp | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/NonLinearSolver/materialLaw/mlawNonLinearTVP.cpp b/NonLinearSolver/materialLaw/mlawNonLinearTVP.cpp index 88b4e0ff3..47241a444 100644 --- a/NonLinearSolver/materialLaw/mlawNonLinearTVP.cpp +++ b/NonLinearSolver/materialLaw/mlawNonLinearTVP.cpp @@ -1394,7 +1394,7 @@ void mlawNonLinearTVP::plasticCorrector_TVP(const STensor3& F0, const STensor3& corKir_e = q1->_kirchhoff; const STensor3& KS = q1->_kirchhoff; - const STensor3& DcorKirDT_TVEcorrection = q1->_DcorKirDT; // Get the corrected value directly from TVE + // const STensor3& DcorKirDT_TVEcorrection = q1->_DcorKirDT; // Get the corrected value directly from TVE getModifiedMandel(Ce, q0, q1); // update Mandel const STensor3& MS = q1->_ModMandel; @@ -2288,14 +2288,15 @@ void mlawNonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow(const STensor3& dDgammaDGamma = kk*(A+Gamma*dAdGamma); // mistake in the paper (VD 2016) this->getYieldCoefficientDerivatives(q1,q1->_nup,Da); - + + /* // DEBUG FLE double a0 = a(0); double a1 = a(1); double a2 = a(2); double da0 = Da(0); double da1 = Da(1); - double da2 = Da(2); + double da2 = Da(2);*/ dfdDgamma = Da(2)*pow(PhiEq,_n) - Da(1)*ptilde -Da(0); dfdDgamma += a(2)*_n*pow(PhiEq,(_n-1.))*He - a(1)*Hp; @@ -3311,7 +3312,7 @@ void mlawNonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow(const STensor3& // mechSourceTVP derivatives - +/* // Numerical Derivatives double T_plus(0.), Wm_TVE_plus(0.), Wm_TVP_plus(0.), dWm_TVEdT_plus(0.), dWm_TVPdT_plus(0.), DKDTsum_plus(0.), DGDTsum_plus(0.); static STensor3 F_plus, Kepr_plus, Mepr_plus, Ke_plus, Me_plus, N_plus, P_plus, DphiPDF_plus, dWm_TVEdF_plus, dWm_TVPdF_plus; @@ -3342,7 +3343,7 @@ void mlawNonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow(const STensor3& dWm_TVEdF_plus(i,j) = (Wm_TVE_plus - Wm_TVE)/_perturbationfactor; dWm_TVPdF_plus(i,j) = (Wm_TVP_plus - Wm_TVP)/_perturbationfactor; } - } + }*/ // conversion tensor for mechSourceTVE Terms static STensor43 DCeDFe, DEeDFe; -- GitLab