From 77631c6397a840f6b75582969591c45de3aabc53 Mon Sep 17 00:00:00 2001 From: "Michael H. Scott" Date: Sun, 3 Dec 2023 09:34:10 -0800 Subject: [PATCH] Update ElasticPlaneStress.cpp --- SRC/material/nD/ElasticPlaneStress.cpp | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) diff --git a/SRC/material/nD/ElasticPlaneStress.cpp b/SRC/material/nD/ElasticPlaneStress.cpp index 572de5ad7c..371aa93c09 100644 --- a/SRC/material/nD/ElasticPlaneStress.cpp +++ b/SRC/material/nD/ElasticPlaneStress.cpp @@ -211,8 +211,14 @@ const Vector& ElasticPlaneStress :: getStrain( ) const Vector& ElasticPlaneStress :: getStress( ) { - stress_vec = this->getTangent() * strain_vec; + //stress_vec = this->getTangent() * strain_vec; + double den = 1 - nu*nu; + double G = E / (2*(1+nu)); + stress_vec(0) = E/den*(strain_vec(0) + nu*strain_vec(1)); + stress_vec(1) = E/den*(strain_vec(1) + nu*strain_vec(0)); + stress_vec(2) = G*strain_vec(2); + return stress_vec ; } @@ -228,7 +234,7 @@ const Matrix& ElasticPlaneStress :: getTangent( ) // double den = 1 - nu*nu; - double G = E / (2*(1-nu)); + double G = E / (2*(1+nu)); tangent_matrix(0,0) = E / den ; tangent_matrix(1,1) = E / den ;