From 3734d0fb077f3328cafb3086f6c68c3013efb806 Mon Sep 17 00:00:00 2001
From: FLE_Knight <ujwalkishore.jinaga@uliege.be>
Date: Wed, 21 May 2025 17:41:22 +0200
Subject: [PATCH] [MINOR PATCH] Bug patched in TMDMN rotated tangents.
---
dG3D/src/dG3DMaterialLaw.cpp | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/dG3D/src/dG3DMaterialLaw.cpp b/dG3D/src/dG3DMaterialLaw.cpp
index c60fb150a..277a12626 100644
--- a/dG3D/src/dG3DMaterialLaw.cpp
+++ b/dG3D/src/dG3DMaterialLaw.cpp
@@ -6192,7 +6192,7 @@ void StochTMDMNDG3DMaterialLaw::stress(IPVariable*ipv, const IPVariable*ipvprev,
for(int k=0; k<3; k++)
for(int l=0; l<3; l++)
for(int p=0; p<3; p++){
- dUdF(m,n,k,l) += I4_first_type(m,p,k,l)*R_def(n,p);
+ dUdF(m,n,k,l) += R_def(p,m)*I4_first_type(p,n,k,l); // R_def_T is used here
for(int q=0; q<3; q++)
dUrot_dU(m,n,k,l) += R(m,p)*I4_symm(p,q,k,l)*RT(q,n);
}
@@ -6223,7 +6223,7 @@ void StochTMDMNDG3DMaterialLaw::stress(IPVariable*ipv, const IPVariable*ipvprev,
// dPdF
static STensor43 dSdF, dSrot_dF, L_rot;
static STensor33 dqdF_rot;
- static STensor3 dqdgradT_rot, dwdf_rot, dmechSourcedf_rot;
+ static STensor3 dwdf_rot, dmechSourcedf_rot;
STensorOperation::zero(L_rot);
STensorOperation::zero(dqdF_rot);
@@ -6271,7 +6271,7 @@ void StochTMDMNDG3DMaterialLaw::stress(IPVariable*ipv, const IPVariable*ipvprev,
STensorOperation::multSTensor3STensor43(dmechSourcedf,dFrot_dF,dmechSourcedf_rot);
// Reassign
- L = L_rot; dqdF = dqdF_rot; dqdgradT = dqdgradT_rot; dwdf = dwdf_rot; dmechSourcedf = dmechSourcedf_rot;
+ L = L_rot; dqdF = dqdF_rot; dwdf = dwdf_rot; dmechSourcedf = dmechSourcedf_rot;
}
STensorOperation::multSTensor3(F_unrot,S,P);
--
GitLab