17 Mai 2005	By LPX

Néant

Ajout d'une fonction à Matr3 (pour assembler les matrices jacobiennes 2D dans une 3D):

OEMTK::Matr3 operator+=(const OEMTK::Matr2 &)

Split Thermomech : (suite)

Nettoyage de la classe Tm2VolumeMetaElement de tout ce qui est reutilisable a partir de VolumeMetaElement (toDofset, fillLocel, npgt, mechanoThermalSfIntegrator,..)

Nettoyage des initialisations des ShapesFunctions - MtSfIntegrator:

- Ajout dans les éléments de fonctions donnant un objet NbIntegrationPoints pour les différents cas (avec utilisations de modèles par defaut a partir des object OptimalGpTemplate)
- Déplacement de l'initialisation des MtSfIntegratorTemplate des classes dérivées (Mt...Integrator : classes supprimées) au profit des nouvelles classes OptimalGp... ceci afin de ne pas introduire de nouveaux copier-coller
- Suppression des fonctions initialiseShapeFunction des éléments
- Suppression des fonctions calculant le nombre de points de Gauss total des éléments vers classes OptimalGp
- Suppression des valeurs du nombre de points de Gauss des VolumeElShcuts

Reste a faire : propager ces modifs vers les elements boundaries & mechanism

Nettoyages des calculs des jacobiens & determinants :

- Remontée du remplissage des tableaux de jacobiens & detJ dans les classes MtSfIntegrator & MtSfIntegrationPoint
- Gestion des Jacobiens 2DEpe/2DAxi/3D au niveau de la classe MtSfIntegrationPoint
- Gestions des detJ <= 0 et rayon < 0 au niveau de MtSfIntegrationPoint
- Suppression des fonctions handleNegativeJacobian & handleNegativeRadius (dans volumeMetaElement && MechanicalVolumeIntegrationMethod)
- Le remplissage du flag Metafor de gestion de l'erreur se faisant au moment ou l'on "catch" l'exception

Matrices de raideurs tangentes numériques :

Observation : les tests de monosThermoMeca a matrice de raideur numérique ont le double d'itérations thermique que ceux a matrice analytique.

Le dper à 1.0e-8 uniforme ne donne pas satisfaction dans tous les cas : on observe une "precision" de la matrice de raideur à
- 1.0e-4 si dper = 1.0e-8
- 1.0e-8 si dper = 1.0e-4

sans oublier la non relation au volume, ...

faut prevoir qqch de mieux

Added Files :

mtShapeFunction/OptimalGpHexa.h/cpp
mtShapeFunction/OptimalGpLine.h/cpp
mtShapeFunction/OptimalGpQuad.cpp
mtShapeFunction/OptimalGpTemplate.h/cpp
mtShapeFunction/OptimalGpTetra.h/cpp
mtShapeFunction/OptimalGpTriangle.h/cpp

Removed Files :

mtIntegrator/MtHexaIntegrator.h/cpp
mtIntegrator/MtLineIntegrator.h/cpp
mtIntegrator/MtQuadIntegrator.h/cpp
mtIntegrator/MtTetraIntegrator.h/cpp
mtIntegrator/MtTriangleIntegrator.h/cpp

mtElements/volumes/Tm2VolumeMetaElement_mec.inl
mtElements/volumes/Tm2VolumeMetaElement_ther.inl

Added tests :

apps/monosThermoMeca(2)/thElastEas2dEpeNum.py
apps/monosThermoMeca(2)/thElastEas3dNum.py

	created :12 mai 2005		modified : 17 Mai 2005
	contact :L.Papeleux@ulg.ac.be