Télécharger le numéro complet

l.IAILLEPréparation du mail lageMAPLANmodule planHATRID I uncAPomodule tridi..mensionnel II module structure/ rrsT i \vérifi cation\ du maillageImauvais maillageI Dessins des résultatsI A. TEST, impression, II perforationIGROUPE ical cu] dugroupe iI Rêsultats imprimés,/ pe"forés ou sur' support magnêtiqueITRACE ireprésentationgraphiqueI Rêsul tats sous forme tf sraphique IFig. 1 Structure du programme ROSALTEpour tout champ de déplacements virtuels admissible(c'est-à-dire vérifiant les conditions de déplacementimposées) ôU, avec les champs de déformations ôe etô0 associés, et pour toute distribution d'incrément decharge ôH admissible (c'est-à-dire vérifiant les conditionsaux limites imposées à la charge hydraulique H),avec le champ de dérivées partielles ôg associé.Après discrétisation de I'espace et du temps (par laméthode de Galerkin), I'analyse du problème débouchesur une procédure de résolution itérative sur labase de l'équation matricielleGo,'Vr*at : Lo, * Gi.. V,avec les notations suivantes :R-CGar:-6rr - 2at KLat :.r*lJ^t :-v -T1-; Rz_6r)cî^12 tF(t) +2F(t+At)l+ ta(t) + 2e(t+ at)lv,*a, : till i îlil,V, : riltil(3)R matrice de rigidité du squelette du sol,C matrice de couplage entre l'écoulement et lesdéformations du squelette,At incrément de temps,F matrice-colonne des forces nodales imposées(forces volumiques et forces de surface confondues),O matrice-colonne des débits nodaux imposés.Cette équation matricielle se résout .. facilement,,, à lacondition d'inverser la matrice Go,, qui est symétriqueet tridiagonalisable par blocs. Dans le programmeRosALlE-Groupe 9, I'inversion est réalisée par laméthode directe de Choleski par bande. La résolutiondonne les déplacements U et les charges H au tempst + at, qui permettent de calculer les autres inconnues,et en particulier les contraintes (les conventions designe du calcul sont celles de la mécanique des milieuxcontinus et non celles de la mécanique des sols).ll faut noter que I'inversion de la matrice Go., qu iconditionne pour une grande part le coût total ducalcul, doit être effectuée chaque fois que I'on modifieI'incrément de tennps at. En pratique, les calculsdoivent pour cette raison n'utiliser qu'un nombre limitéde valeurs de At dans chaque cas traité.L'introduction de l'élastoplasticité dans la procédurede calcul de la consolidation a été effectuée sousforme incrémentale, ên su pposant que la fonctionF(o', k) : 0 est à la fois le critère de plasticité et lepotentiel plastique (c'est-à-dire qu'on lui applique leprincipe de normalité), o' désignant l'état des contrainteseffectives et k un paramètre d'écrouissage.REVUE FRANçA|SE DE GEOTECHNTOUE NUMERO 19 40