From 8b3eecc7a3b6c48e7eb599b02e4894ec8e8493fb Mon Sep 17 00:00:00 2001
From: "vinugholap@gmail.com" <Vinugitlab$123>
Date: Wed, 17 Mar 2021 16:57:49 +0100
Subject: [PATCH] More corrections
---
NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp | 2 --
dG3D/src/dG3DTerms.cpp | 8 ++++----
2 files changed, 4 insertions(+), 6 deletions(-)
diff --git a/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp b/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp
index ec6fffe42..7fd097bf7 100644
--- a/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp
+++ b/NonLinearSolver/materialLaw/mlawLinearElecMagTherMech.cpp
@@ -394,7 +394,6 @@ void mlawLinearElecMagTherMech::constitutive(
dqdV(i) += 0.0;
djydT(i) += (_seebeck + dseebeckdT * T) * js0(i);
- //djydV(i) += 0.0;
// derivatives of sourceVectorField
dsourceVectorFielddV(i) += 0.0;
@@ -405,7 +404,6 @@ void mlawLinearElecMagTherMech::constitutive(
dfluxDdA(i,j) += 0.0;
dfluxTdA(i,j) += 0.0;
djydA(i,j) += 0.0;
- djydV(i) += _l0(i, j) * (An(j) - A0(j)) * invdh;
dsourceVectorFielddA(i, j) += 0.0;
dsourceVectorFielddB(i, j) += 0.0;
diff --git a/dG3D/src/dG3DTerms.cpp b/dG3D/src/dG3DTerms.cpp
index d2a81743a..98e1167cf 100644
--- a/dG3D/src/dG3DTerms.cpp
+++ b/dG3D/src/dG3DTerms.cpp
@@ -1859,13 +1859,13 @@ void dG3DStiffnessBulk::get(MElement *ele,int npts,IntPt *GP,fullMatrix<double>
dHdB.operator()(m, n) *
CurlcurlVals[k + numEdgeShapeFuncTotalLastCol][n]);
mStiff(j + numEdgeShapeFuncTotalLastRow + numStandardDof,
- k + numEdgeShapeFuncTotalLastCol + numStandardDof) -=
+ k + numEdgeShapeFuncTotalLastCol + numStandardDof) +=
getConstitutiveCurlEqRatio() * ratio *
(curlVals[j + numEdgeShapeFuncTotalLastRow][m] *
dSourceVectorFielddVectorPotential.operator()(m, n) *
curlVals[k + numEdgeShapeFuncTotalLastCol][n]);
mStiff(j + numEdgeShapeFuncTotalLastRow + numStandardDof,
- k + numEdgeShapeFuncTotalLastCol + numStandardDof) -=
+ k + numEdgeShapeFuncTotalLastCol + numStandardDof) +=
getConstitutiveCurlEqRatio() * ratio *
(curlVals[j + numEdgeShapeFuncTotalLastRow][m] *
dSourceVectorFielddVectorCurl.operator()(m, n) *
@@ -1994,7 +1994,7 @@ void dG3DStiffnessBulk::get(MElement *ele,int npts,IntPt *GP,fullMatrix<double>
(dVectorFielddExtraDofField.operator()(m) *
CurlcurlVals[j+numEdgeShapeFuncTotalLastRow][m] *
Vals[k+nbFFTotalLastCol]);
- mStiff(j+numEdgeShapeFuncTotalLastRow+numStandardDof,k+nbFFTotalLastCol) -= getConstitutiveCurlEqRatio() * ratio *
+ mStiff(j+numEdgeShapeFuncTotalLastRow+numStandardDof,k+nbFFTotalLastCol) += getConstitutiveCurlEqRatio() * ratio *
(dSourceVectorFielddExtraDofField.operator()(m) *
curlVals[j+numEdgeShapeFuncTotalLastRow][m] *
Vals[k+nbFFTotalLastCol]);
@@ -2004,7 +2004,7 @@ void dG3DStiffnessBulk::get(MElement *ele,int npts,IntPt *GP,fullMatrix<double>
(CurlcurlVals[j+numEdgeShapeFuncTotalLastRow][m] *
dVectorFielddGradExtraDofField.operator()(m,n) *
Grads[k+nbFFTotalLastCol][n]);
- mStiff(j+numEdgeShapeFuncTotalLastRow+numStandardDof,k+nbFFTotalLastCol) -= getConstitutiveCurlEqRatio() * ratio *
+ mStiff(j+numEdgeShapeFuncTotalLastRow+numStandardDof,k+nbFFTotalLastCol) += getConstitutiveCurlEqRatio() * ratio *
(curlVals[j+numEdgeShapeFuncTotalLastRow][m] *
dSourceVectorFielddGradExtraDofField.operator()(m,n) *
Grads[k+nbFFTotalLastCol][n]);
--
GitLab