Création de tenseurs ordre 4 en 3D et 2D. Il y a une classe générale Tens4_2 et Tens4_3 qui ne contient que la dimension et les paramètres ij33 et 33kl pour Tens4_2. Ces derniers éléments permettent de tenir compte des termes hors plan. Il dérive de ces classes des tenseurs isotropes Tens4_iso2 et Tens4_iso3 qui contiennent les termes en dijdkl, dildjk, dikdjl. Et des tenseur sparses Tens4_spa2 et Tens4_spa3 qui contiennent tous les termes. Outre les fonction +=, =… qui sont définies, il y a les fonctions add_dik_djl …, dérivée de jauman, Vect*Tens*Vect …
La matrice de raideur par défaut (itan=0) est maintenant très simple (si ce n'est qu'il faut tenir compte du report de pression ce qui triple le nombre de lignes). Il y a un appel au materiau de get_hijkl_mate, get _hijkl_geom puis l'opération BIj Hijkl BJl.
L'ancienne méthode peut encore être appelée par itan=2.
Il faut encore programmer les get_hijkl pour les matériau hyperélastique.