Dans le projet ICOGEN, j'ai rencontrés des problèmes de convergence récurrents avec les matériaux orthotropes élastoplastiques avec ou sans endommagement (EpIsoHOrthoHypoMaterial et DamageEpIsoHOrthoHypoMaterial). Les problèmes se manifestaient au niveau du Newton-Raphson de la correction plastique mais uniquement avec certaines courbes d'écrouissage, ce qui explique que je ne m'en étais pas rendu compte avant. J'ai modifié les matériaux orthotropes et les critères de plasticité de façon à rendre la procédure plus robuste. Le calcul est sans doute un peu plus lent mais je n'ai plus de problèmes de convergence depuis les modifications.

Il est maintenant possible de définir les axes d'orthotropie du matériau au moyen d'une fonction des coordonnées spatiales. C'est par exemple nécessaire pour des structures composites ondulées, comme dans le cas-test apps.composite.lam2DCompBase. La syntaxe est un peu lourde, en particulier parce qu'il faut définir séparément les trois composantes de deux des axes d'orthotropie.

axis1X = PythonMultiParameterFunction( [fonction x1(x,y,z)], 3 )
axis1Y = PythonMultiParameterFunction( [fonction y1(x,y,z)], 3 )
axis1Z = PythonMultiParameterFunction( [fonction z1(x,y,z)], 3 )
material.put( ORTHO_AX1_X, 1.0 )
material.put( ORTHO_AX1_Y, 1.0 )
material.put( ORTHO_AX1_Z, 1.0 )
material.depend( ORTHO_AX1_X, axis1X, FieldList( Field1D(TX,AB), Field1D(TY,AB), Field1D(TZ,AB) ) )
material.depend( ORTHO_AX1_Y, axis1Y, FieldList( Field1D(TX,AB), Field1D(TY,AB), Field1D(TZ,AB) ) )
material.depend( ORTHO_AX1_Z, axis1Z, FieldList( Field1D(TX,AB), Field1D(TY,AB), Field1D(TZ,AB) ) )

Des matériaux identifiés par l'UCL ou trouvés dans la littérature et utilisés dans le projet ICOGEN ont été ajoutés dans apps/toolbox/createCompositeMaterial.py:

• fibres carbone isotropes et orthotropes (AS4, IM7, G0926).
• résine RTM6 identifiée par l'UCL à plusieurs températures.
• pli composite unidirectionnel composé de fibres IM7 et de résine RTM6, identifié par homogénéisation.

[a]:
[r]:

[a]: mtParasolid\tests\tools\numisheet1993BM1.py
[a]: mtParasolid\tests\tools\numisheet1993BM1_full_die.xmt_txt
[a]: mtParasolid\tests\tools\numisheet1993BM1_full_holder.xmt_txt
[a]: mtParasolid\tests\tools\numisheet1993BM1_full_punch.xmt_txt
[a]: mtParasolid\tests\tools\numisheet1993BM1_quarter_die.xmt_txt
[a]: mtParasolid\tests\tools\numisheet1993BM1_quarter_holder.xmt_txt
[a]: mtParasolid\tests\tools\numisheet1993BM1_quarter_punch.xmt_txt
[a]: mtParasolid\tests\tools\numisheet2011BM2.py
[a]: mtParasolid\tests\numisheet1993BM1.py
[a]: mtParasolid\tests\numisheet2011BM2.py
[r]: mtParasolid\tests\numisheet\

— Geoffrey 2017/03/10