element finis avec 2 elements triangle

[diallokoula]

Bonjour,

j'utilise un code element finis triangle quadratique mais il se trouve que les deplacements(meme les deformations et contraintes) ne sont pas bonnes comparés à Abaqus.
Pour comparer les deux j'ai fait un essai de traction uniaxiale avec une pression imposé sur l'une des faces (par rapport à x) et le face opposé j'ai mis une condition limite en deplacement ( U1=0 ) suivant x. j'ai encastré un point de cette face pour des conditions de symétrie.

maintenant j'aimerai comparé les deplacements et deformations avec les resultats theoriques sur un rectangle avec 2 elements.
j'utilise l'hypothese de contraintes planes.

quelqu'un aurait il une idée pour le calcul en element finis d'un rectangle avec 2 elements finis et une interpolation quadratique.

je demarre en EF et j'ai du mal.

Merci