1025x
001935
11.02.2025

Calcul de déformation directe d’une poutre en béton armé en tenant compte du fluage et du retrait

Cet article technique traite du calcul direct de la déformation d’une poutre en béton armé en tenant compte des effets à long terme de fluage et de retrait. Le calcul direct est expliqué conformément à l’Eurocode 2 (EN 1992-1-1, 7.4.3) à l’aide d’une poutre à travée simple. Les armatures à la traction, le comportement en état fissuré en utilisant le facteur de distribution (paramètre d’endommagement) ainsi que la prise en compte du comportement au retrait et au fluage font l’objet d’une attention particulière.

Dans cet article technique, le calcul direct de la déformation d’une poutre en béton armé est effectué, en prenant également en compte les effets à long terme de fluage et de retrait. Un exemple montre comment ces effets influencent la déformation d’un élément structurel et sont intégrés dans le calcul. L’article détaille également les entrées requises pour considérer correctement tous les facteurs pertinents dans RFEM 6 et comment le facteur de distribution influence la rigidité de l’élément.

Données d’entrée

La géométrie, les armatures et les charges sont décrites par les paramètres suivants :

Système

  • Type de poutre : Poutre à travée simple
  • Portée : l = 4,210 m

Section

  • Épaisseur de la dalle : h = 20 cm
  • Largeur de la dalle b = 100 cm
  • Matériau : béton C20/25 avec Ecm = 30 000 MN/m² et B 500A
  • Armatures : As,-z,(bas) = 4,45 cm² avec 7 ∅ 9 et d1 = 30 mm
  • Hauteur utile des armatures inférieures : ddef,+z (bas) = 17 cm

Charges permanentes

  • Poids propre : gs = 0,20 m ⋅ 1m ⋅ 25 kN/m³ = 5,00 kN/m
  • Revêtement et enduit : gbp = 1,50 kN/m
  • Total : gk,total = 6,5 kN/m

Charges variables

  • Charge d’exploitation (bureau) : qb = 2,00 kN/m avec ψ2 = 0,3
  • Cloison : qt = 1,25 kN/m avec ψ2 = 1,0

Charge quasi-permanente

  • 6,5 kN/m + 0,3 ⋅ 2,00 kN/m + 1,0 ⋅ 1,25 kN/m = 8,35 kN/m

Moment fléchissant de calcul pour le calcul de la flèche

  • My,Ed,def = 8,35 kN/m ⋅ (4,21 m)² / 8 = 18,50 kNm

Valeurs initiales pour le calcul des déformations

  • Module d’élasticité moyen du béton : Ecm = 30.000 MN/m²
  • Taux d’armatures longitudinales : ρ = As / Ac = 4,45 cm² / ( 20 cm ⋅ 100 cm) = 0,223 %
  • Déformation due au retrait : εsh = -0,5 ‰
  • Coefficient de fluage : φ = 2

L'option « Propriétés dépendantes du temps avancées du béton » doit être activée dans les paramètres de section pour pouvoir définir le coefficient de fluage de manière personnalisée.

Dans l’onglet désormais disponible, cochez d’abord « Fluage » et « Retrait » pour pouvoir visualiser et modifier les « valeurs fondamentales des propriétés dépendantes du temps ». Le facteur de fluage φ a été ici défini par l’entrée de φ0, ϵcd,0 et ϵca(∞).

Fluage

Les effets de fluage sont pris en compte par une réduction du module d’élasticité Ec du béton.

Le module d’élasticité effectif Ec,eff prend en compte les effets à long terme du béton, notamment le fluage. Le fluage décrit la déformation à long terme du béton sous une charge constante. Grâce au facteur de fluage φ, le module d’élasticité Ecm (module d’élasticité moyen du béton) est réduit pour que la rigidité réelle du béton soit représentée sur une longue période. Cette valeur est utilisée dans les calculs tels que le moment d’inertie ou le rapport des rigidités.


Ec,eff = 30.000 MN/m² / ( 1 + 2 ) = 10.001,2 MN/m²

Module de cisaillement effectif du béton Gc,eff
Le module de cisaillement effectif décrit la capacité du béton à résister aux déformations de cisaillement et est déterminé par le rapport de déformation transversale à la déformation longitudinale (le coefficient de Poisson du béton). Cette valeur est particulièrement importante pour les calculs de déformations de section et les vérifications de la résistance au cisaillement.


Gc,eff = 10.001,2 MN/m² / ( 2 ⋅ ( 1 + 0,2 ) ) = 4.167,180 MN/m²

Rapport des modules d’élasticité pour l’état non fissuré (charge à long terme) αe,l
Le rapport αe,l indique la rigidité relative de l’acier par rapport au béton sous charge à long terme. Es est le module d’élasticité de l’acier, Ec,l est le module d’élasticité effectif du béton dans l’état non fissuré (identique à Ec,eff). Étant donné que le béton présente une rigidité réduite due au fluage, la valeur de αe,l est plus élevée dans cet état. Ce rapport est utilisé pour le calcul du centre de gravité et des valeurs de section effectives.


αe,l = 2 ⋅ 105 MN/m² / 10.001,2 MN/m² = 20

Rapport des modules d’élasticité pour l’état non fissuré (charge à court terme) αe,I,st
Le rapport αe,I,st décrit le rapport de rigidité de l’acier au béton sous charge à court terme. Contrairement à αe,l, le module d’élasticité moyen Ecm est utilisé sans prendre en compte les effets de fluage. Cela reflète la situation de charge réelle lorsque le béton est soumis à une charge à court terme. Cette valeur est particulièrement pertinente pour la vérification des charges à court terme.


αe,I,st = 2 ⋅ 105 MN/m² / 30.000 MN/m² = 6,67

Rapport des modules d’élasticité pour l'état fissuré αe,II
Dans l’état fissuré, le béton dans la zone tendue est considéré comme non porteur. Le rapport αe,II prend cela en compte en incluant seulement le module d’élasticité effectif Ec,eff du béton. Cette valeur montre que la rigidité de l’acier est plus élevée par rapport au béton dans l’état fissuré, ce qui souligne l’importance des armatures dans de tels cas.


αe,II = 2 ⋅ 105 MN/m² / 10.001,2 MN/m² = 20,00

Paramètres géométriques non fissurés

Ladistance du centre de gravité de la section effective dans l’état non fissuré sous charge à long terme, zI, décrit la position du centre de gravité en tenant compte de la surface du béton et des armatures. L’influence des armatures est mise à l’échelle avec un facteur de conversion αe,l, qui représente le rapport entre le module d’élasticité de l’acier et le module d’élasticité effectif du béton. Ceci est particulièrement important car les charges à long terme, telles que le fluage, affaiblissent le béton. Le centre de gravité influence le calcul des moments et des déformations de la section et est donc un paramètre central pour l’analyse statique.


L’aire effective de la section dans l’état non fissuré sous charge à long terme, AI représente l’aire effective qui supporte les charges. Outre la surface du béton, l’aire d’armatures est également prise en compte, complétée par le facteur αe,l. Cela représente plus précisément la rigidité de la section. Cette valeur est cruciale pour évaluer la capacité portante et pour le calcul de la déformation de l’élément.


AI = 1.000 mm ⋅ 200 mm + 20 ⋅ ( 4,45 cm² + 0 cm² ) = 2.089,05 cm²

Le moment d’inertie effectif du centre de gravité idéal dans l’état non fissuré sous charge à long terme, II décrit la résistance de la section en flexion. Il prend en compte à la fois la surface du béton et les armatures, cette dernière générant des moments supplémentaires par rapport à sa position relative au centre de gravité. Ce moment d’inertie est un facteur central pour le calcul des déformations et montre dans quelle mesure la section peut résister aux moments fléchissants.


L’excentrement du centre de gravité idéal dans l’état non fissuré, eI indique l’écart du centre de gravité par rapport au centre géométrique de la section. Cet excentrement est important car il influence les moments qui se développent dans la section, ce qui affecte directement les déformations.


eI = 103 mm - 200 mm / 2 = 3 mm

La distance du centre de gravité de la section idéale dans l’état non fissuré sous charge à court terme, zI,st décrit la position du centre de gravité sous des charges qui ne tiennent pas compte des effets de fluage ou de retrait. Le facteur de conversion αe,I,st, utilisé dans le calcul à court terme, est donc inférieur à celui des charges à long terme. Cette distance du centre de gravité est cruciale pour la distribution des charges et la détermination des moments sous charges à court terme.


L’aire effective de la section dans l’état non fissuré sous charge à court terme, AI,st est similaire à l’aire AI, mais ajustée par le facteur de conversion αe,I,st, qui ne tient pas compte des effets à long terme. Cela conduit à une surface réduite et a un impact sur le calcul de la capacité portante sous charges à court terme.


AI,st = 1000 mm ⋅ 200 mm + 6,67 ⋅ ( 4,45 cm² + 0 cm² ) = 2.029,69 cm²

Le moment d’inertie effectif du centre de gravité idéal dans l’état non fissuré sous charge à court terme, II,st représente la résistance de la section en flexion sans influence des effets à long terme. Les aires d’armatures et les surfaces en béton ainsi que leurs distances par rapport au centre de gravité sont prises en compte, ce qui est essentiel pour le calcul des déformations sous charges à court terme.


Paramètres géométriques fissurés

La distance du centre de gravité de la section idéale dans l’état fissuré, zII prend en compte la capacité portante modifiée de la section, car la zone tendue du béton ne supporte plus de charge après fissuration. La position du centre de gravité est recalculée en considérant uniquement la zone comprimée du béton et les armatures. Ce paramètre est crucial pour l’analyse de la section après fissuration et influence la capacité portante et la déformation.


L’aire effective de la section dans l’état fissuré, AII représente la surface restante après fissuration. Seule la zone de compression du béton et l’aire des armature sont prises en compte, ce qui réduit considérablement la rigidité de la section. Cette valeur est cruciale pour l’analyse de la capacité portante des sections fissurées.


AII = 1.000 mm ⋅ 46,8 mm + 20 ⋅ ( 4,45 cm² + 0 ) = 557,41 cm²

Le moment d’inertie effectif du centre de gravité idéal dans l’état fissuré, III décrit la résistance en flexion après fissuration. Comme la zone tendue n’est plus porteuse, le moment d’inertie est considérablement réduit. Cette valeur est un facteur essentiel pour le calcul des déformations et l’évaluation de la capacité portante des sections fissurées.


L’excentrement du centre de gravité idéal dans l’état fissuré, eII décrit le déplacement du centre de gravité dû à la fissuration. Ce déplacement influence les moments et les déformations de la section et est donc un paramètre important pour l’analyse statique.


eII = 46,8 mm - 200 mm / 2 = -53,2 mm

Retrait

L’effort normal dû au retrait, Nsh se développe parce que les armatures ne participent pas à la déformation due au retrait du béton et donc absorbe des forces. Ces forces résultent de l’interaction entre l’effort de traction du béton et la réaction des armatures. La valeur calculée montre à quel point les armature sont sollicitées par le retrait. À ce stade, la déformation due au retrait définie indirectement par l’utilisateur εsh = -0,5 ‰ est utilisée.


Nsh = -2 ⋅ 105 MN/m² ⋅ ( -0,000,5 ) ⋅ ( 4,45 cm² + 0,00 ) = 44,532 kN

L’excentrement de la force due au retrait par rapport au centre de la section idéale à l'état non fissuré, esh,I décrit la position de la force résultante de retrait par rapport au centre de la section. Une plus grande excentricité mène à des moments plus élevés et à de plus grandes déformations.


esh,I = ( 4,45 cm² ⋅ 170 mm + 0 ) / ( 4,45 cm² + 0 ) - 103 mm = 67 mm

Le moment de retrait pour l’état non fissuré, Msh,I résulte de la force due au retrait Nsh et de l’excentrement esh,I. Il décrit comment la force de retrait génère un moment sur la section. Ce moment influence considérablement les déformations et les tensions dans la section et doit être pris en compte dans la conception.


Msh,I = 44,532 kN ⋅ 67 mm = 2,98 kNm

La courbure pour l’état non fissuré ksh,I indique comment le moment de retrait agit par rapport à l’effort normal et à l’excentrement. Elle montre comment la distribution de la force due au retrait et la position du centre de gravité influencent les déformations de l’élément. Cette valeur est essentielle pour décrire complètement les déformations de la section dues au retrait.


ksh,I = ( 2,98 kNm + 18,5 kNm - 0 ) / ( 18,50 kNm - 0 ) = 1,161

L’excentrement de la force due au retrait par rapport au centre de la section idéale dans l’état fissuré, esh,II décrit la position de la force résultante de retrait par rapport au centre de la section dans l’état fissuré. Les moments de l’aire des armatures, As,def, +z,(bas) et As,def, -z,(haut), sont considérés en fonction de leur position, def, +z,(bas) et def, -z,(haut), et divisés par l’aire totale des armatures. Le résultat est soustrait de la distance au centre de la section fissurée, zII. Cet excentrement influence le moment de retrait, car une excentricité plus grande génère un moment plus important.


esh,II = ( 4,45 cm² ⋅ 170 mm + 0 ) / ( 4,45 cm² + 0 ) – 46,8 mm = 123,2 mm

Le moment fléchissant dû à l’effort normal Nsh pour l’état fissuré, Msh,II est obtenu en multipliant la force de retrait Nsh par l’excentrement précédemment calculé esh,II. Ce moment décrit la sollicitation supplémentaire en flexion que la force de retrait induit sur la section. Surtout dans l’état fissuré, où la zone tendue du béton ne supporte plus de charges, cette valeur est pertinente.


Msh,II = 44,532 kN ⋅ 123,2 mm = 5,48 kNm

La courbure pour l’état fissuré, ksh,II indique dans quelle mesure la déformation de la section est influencée par le moment de retrait et les autres forces agissantes. Le moment de retrait Msh,II, le moment fléchissant existant My,Ed,def, ainsi que l’effort normal NEd et son excentrement eII sont pris en compte. Le calcul met en relation le moment résultant avec le moment sans retrait, fournissant ainsi une mesure de l’influence de la force de retrait.


ksh,II = ( 5,48 kNm + 18,50 kNm - 0 ) / ( 18,50 kNm - 0 ) = 1,296

Déformation de la section

La déformation de la section décrit la courbure d’un élément causée par des influences externes en tenant compte de ses paramètres de matériau et d’état.

Le calcul de la déformation de la section dans l’état non fissuré, κI décrit la courbure de la section causée par le moment de retrait et les propriétés élastiques du matériau. Le moment de retrait My,Ed,def est pris en compte, ainsi que l’effort normal NEd et son excentrement eI. Ces quantités sont multipliées par le facteur ksh,I, qui décrit l’influence du moment de retrait dans l’état non fissuré. Au dénominateur figurent le module d’élasticité effectif du béton, Ec,eff, et le moment d’inertie de la section non fissurée, II, qui déterminent la rigidité de la section.


κI = ( 1,161 ⋅ ( 18,50 kNm - 0) ) / ( 10.001,2 MN/m² ⋅ 70.844,30 cm⁴ )
= 3 mrad/m

Le calcul de la déformation de la section dans l’état fissuré, κII montre la courbure de la section après fissuration, en tenant compte du moment de retrait et de la capacité portante réduite de la section fissurée. Le moment de retrait My,Ed,def, l’effort normal NEd et son excentrement eII sont multipliés par le facteur ksh,II, qui décrit l’influence du moment de retrait dans l’état fissuré. Au dénominateur figurent le module d’élasticité effectif du béton, Ec,eff, et le moment d’inertie réduit de la section fissurée, III, qui reflètent la moindre rigidité de la section. La déformation de la section dans l’état fissuré est nettement plus grande que dans l’état non fissuré, car la rigidité de la section fissurée est réduite.


κII = ( 1,296 ⋅ ( 18,50 kNm - 0 ) ) / ( 10.001,2 MN/m² ⋅ 16.933,50 cm⁴ ) = 14,2 mrad/m

État final

L’état final décrit les contraintes maximales qui peuvent se produire dans la section non fissurée tant sous charges à long terme que sous charges à court terme, pour garantir la capacité portante et la fonctionnalité de l’élément.

La contrainte maximale dans l’état non fissuré sous charge à long terme, σmax,It décrit la contrainte la plus élevée qui peut se produire dans la section non fissurée en raison des charges à long terme. Elle se compose de deux contributions :

  • la contribution des efforts normaux NEd et Nsh
  • la contribution des moments fléchissant My,Ed,def, Msh,I et celui résultant de l’excentrement (zI - h/2) de la force normale NEd.

La deuxième contribution est amplifiée par le moment d’inertie II et la distance (h - zI).


La contrainte maximale dans l’état non fissuré sous charge à court terme, σmax,st donne la plus grande contrainte dans la section sous charges à court terme. Contrairement aux charges à long terme, seuls les effort normaux NEd et My,Ed,def sont pris en compte ici, car il n'y a pas de charges dues au retrait.


La contrainte maximale à l'état non fissuré, σmax est la plus grande des deux valeurs de contrainte sous charges à long terme et charges à court terme. Elle garantit que la sollicitation maximale possible de la section est prise en compte.


σmax = max( 3,155 MN/m²; 2,689 MN/m² )
= 3,155 MN/m²

Le coefficient de distribution (paramètre d’endommagement), ζd décrit la transition entre le comportement de la section dans l’état non fissuré et dans l’état fissuré. Il est calculé par le rapport de la résistance à la traction caractéristique du béton, fctm, à la contrainte maximale, σmax. La non-linéarité est prise en compte par la relation exponentielle.


ζd = 1 – 0,5 ⋅ (2,200 MN/m² / 3,155 MN/m²)²
= 0,757 ≤ 1
avec :
β = 1,0 (charge à court terme)
β = 0,5 (charge à long terme ou plusieurs cycles de sollicitations répétées)
Si le coefficient de distribution ζd = 1, l'élément est complètement à l'état fissuré. Si ζd est égal à 0, le béton est entièrement non fissuré.

Informations

Le calcul direct de la déformation est fortement dépendant du coefficient de distribution. Pour en savoir plus, consultez cet article technique.

Pour le calcul du coefficient de distribution ζd, il est important de savoir quelle option a été choisie pour la détection de l’état fissuré. L'option « État fissuré calculé à partir de la charge associée » signifie que l’état fissuré (coefficient de distribution ζd) est calculé uniquement à partir de la charge actuelle (combinaison de charges) — comme dans cet exemple. Les autres options sont détaillées dans le manuel.

La courbure de la section, κf est calculée par interpolation entre l’état fissuré (κII) et l'état non fissuré (κI), pondérée par le coefficient de distribution, ζd. Cela permet une description réaliste du comportement de courbure dans l’état de transition.


κf = 0,757 ⋅ 14,2 mrad/m + (1 – 0,757) ⋅ 3 mrad/m
= 11,5 mrad/m

L’aire de section idéale, Af décrit la transition entre l’aire de la section non fissurée, AI, et l’aire de la section fissurée, AII. La pondération est également faite par le coefficient de distribution, ζd.


Le moment d’inertie idéal, Iy,f décrit le moment de la section en tenant compte du coefficient de distribution, ζd, ainsi que des moments d'inertie dans l’état non fissuré, II, et à l'état fissuré, III. Des facteurs supplémentaires comme ksh,II et ksh,I prennent en compte les effets du retrait dans chaque état.


L'excentrement du centre de gravité, ef décrit la position du centre de gravité résultant de la section, basée sur la transition entre l'état non fissuré et l'état fissuré. Elle prend en compte le coefficient de distribution, ζd, ainsi que les modules d'élasticité, Ec,eff, et les moments d'inertie, II et III.


Le moment d'inertie idéal par rapport au centre géométrique de la section, Iy,0,f prend en compte, en plus du moment d’inertie idéal, Iy,f, et de l’aire de section idéale, Af, le déplacement du centre de gravité causé par l’excentrement, ef. Ce déplacement est pris en compte via la contribution de Steiner de Af.


Iy,0,f = 16.145,50 cm⁴ + 678,30 cm² ⋅ ( -49,2 mm )²
= 32.538,80 cm⁴

Rigidités finales

Les rigidités finales d’un élément décrivent sa résistance aux déformations et aux rotations sous différentes types de charges. Elles comprennent à la fois les rigidités axiales, les rigidités en flexion, ainsi que les rigidités en torsions et les rigidités en cisaillement. Ces valeurs servent de base pour l’analyse du comportement porteur et de la fonctionnalité d’un élément.

La rigidité tangentielle de membrane, EAf décrit la rigidité axiale de la section en prenant en compte le module d’élasticité effectif du béton, Ec,eff, et l’aire de section efficace, Af.


EAf = 10.001,2 MN/m² ⋅ 678,30 cm²
= 678.387 kN

La rigidité tangentielle en flexion, EIy,0,f décrit la résistance de la section en flexion autour du centre de gravité idéal. Elle est déterminée par le module d’élasticité effectif du béton, Ec,eff, et le moment d’inertie idéal, Iy,0,f.


EIy,0,f = 10.001,2 MN/m² ⋅ 32.538,80 cm⁴
= 3.254,28 kNm²

La rigidité tangentielle en flexion, EIz,0,f décrit la résistance de la section en flexion autour de l’axe local z. Elle est définie par le module d’élasticité effectif du béton, Ec,eff, et le moment d’inertie autour de l’axe z, Iz.


EIz,0,f = 10.001,2 MN/m² ⋅ 1.666.670 cm⁴
= 166.687 kNm²

Le facteur r décrit la réduction de la rigidité en cisaillement basée sur le rapport des moments d’inertie idéaux If et II.


r = 16.145,50 cm⁴/70.844,30 cm⁴
= 0,228

La rigidité en cisaillement par rapport à l’axe y, GAy,f prend en compte le module de cisaillement effectif du béton, Gc,eff, l’aire de la section, Ac,y, ainsi que le facteur de réduction r.


GAy,f = 4.167,18 MN/m² ⋅ 1.666,67 cm² ⋅ 0,228
= 158284 kN

La rigidité en cisaillement par rapport à l’axe z, GAz,f est calculée de manière analogue par rapport à l’axe y.


GAz,f = 4.167,18 MN/m² ⋅ 1.666,67 cm² ⋅ 0,228
= 158.284 kN

La rigidité en torsion, GIT,f correspond dans le cas considéré à la rigidité en torsion dans l’état non fissuré, GIT,I.


GIT,f = 7.770 kNm²

L’élément de rigidité excentré, ESy décrit la sollicitation supplémentaire de la section causée par l’excentrement ef. Il est calculé par la rigidité axiale EAf et l’excentrement ef.


ESy = 678.387 kN ⋅ ( -49,2 mm )
= -33.350,20 kNm

Flèche

Pour garantir la fonctionnalité, la déformation réelle est comparée aux limites admissibles. La flèche totale est corrigée par la contre-flèche et contrôlée par rapport aux limites données.

Dans le calcul de la flèche, la combinaison de charges principale, sans effets dépendants du temps tels que le fluage et le retrait (court terme), est prise en compte, tandis que les combinaisons de charges associées sont toujours calculées avec les propriétés dépendantes du temps (long terme). S’il y a plus d’une charge associée, la flèche de la charge ayant la plus grande valeur est prise en compte.

La flèche limite en direction z, uz,lim est calculée avec la longueur de référence en direction z, Lz,ref, et le critère de flèche limite, Lz,ref/uz,lim.


uz,lim = 4,210 m / 250 = 16,8 mm

La flèche en direction z, uz est le résultat de la différence entre la flèche totale, uz,ges, et le contre flèche à la position x, uz,c.


uz = 19,4 mm - 0 = 19,4 mm

Vérification

η = max( 19,4 mm / 16,8 mm; 0,0 mm / 16,8 mm ) = 1,155 Puisque η = 1,155 > 1, la flèche admissible est dépassée !

Conclusion

Le calcul de déformation selon les méthodes d’approximations spécifiées dans les normes, tel que lealcul de déformation conformément à l’EN 1992-1-1, 7.4.3, est réalisé à l’aide des rigidités effectives qui sont calculées dans les éléments finis en fonction de l’état limite (fissuré ou non fissuré). Ces rigidités effectives servent de base pour le calcul ultérieur de la déformation de l’élément par une autre analyse aux éléments finis.

Pour la détermination des rigidités effectives, la section de béton armé est considérée, le section de béton armé étant classé comme « fissuré » ou « non fissuré » en fonction des efforts internes déterminés pour l’état limite de service. La participation du béton entre les fissures est prise en compte par un coefficient de distribution, par exemple selon l’équation 7.19 (EN 1992-1-1). Le comportement de matériau du béton est considéré comme linéaire-élastique jusqu'à la résistance en traction du béton, ce qui suffit pour l’ELS.

La prise en compte des effets à long terme de fluage et de retrait est effectuée lors de la détermination des rigidités effectives au niveau de la section de l’élément, afin de garantir une représentation réaliste des déformations sous charges à long terme.


Auteur

Madame Dannwerth assiste les utilisateurs dans le support client et se consacre au développement dans le domaine de la géotechnique.



;