problem with Stress

For developers writing C++, Fortran, Java, code who have questions or comments to make.

Moderators: silvia, Moderators

problem with Stress

Postby mohammmad » Thu Oct 11, 2018 1:25 am

Hi all,
I have an 2ِD element that in each section of the element there are three component: strainXX, strainYY, and strainXY. I use ND material, first plain-strain J2 Plasticity Material for steel section and second plain-strain Drucker Prager Material.
Because of the specific cross - section conditions, I set the variable (strainYY) in such a way that stressYY is 0.0. In order to set stressYY=0.0, I got the tangent stiffness matrix(getTangent) of the material each time and did this:

see th picture :
const Matrix& newK = theMat->getTangent();
strain(1) = -((strain(0)*newK(1, 0) + strain(2)*newK(1, 2)) / (newK(1, 1)));

But, the problem is when I set this new strain vector with "theMat->setTrialStrain(strain);" and after that, I call "theMat->getStress();", in the resultant stress vector stressYY (in vector is stress(1)) is not zero.
Posts: 19
Joined: Tue Aug 21, 2018 4:10 am

Return to Framework

Who is online

Users browsing this forum: No registered users and 1 guest