Dear FEA way, thanks for your reply. I have check my matrix assembly and element formulation and they look ok!
I am able to fix my code to some extent but have still some issue. What I mean is that when I consider one material (e.g a fluid) for my whole domain (9 elements), my plain strain...
Why in the case of plane strain results from FEM code in Matlab and Abaqus don't match whereas the plane stress results do?
I am trying to simulate a simple pure shear problem while considering a viscous fluid with my own FEM code and I compare it with Abaqus. I use 4 node bilinear quadrilateral...