Comment modéliser la diffusion des charges de roues à travers le remblai dans Robot Structural Analysis ?
Lors du dimensionnement d'ouvrages enterrés (cadres en béton armé, buses, ouvrages de soutènement, structures pousse-tubes), la prise en compte des charges routières ou ferroviaires en surface est une étape critique. Cependant, une roue de camion n'applique pas la même pression sur une dalle située à $0,50\text{ m}$ de profondeur que sur une dalle située à $4,00\text{ m}$.
Le sol agit comme un milieu répartiteur. Autodesk Robot Structural Analysis (RSA) possède un moteur d'analyse par Éléments Finis (EF) très puissant, mais il ne calcule pas automatiquement la diffusion d'une force à travers un massif de sol 3D si l'on modélise uniquement la structure en béton.
Cet article détaille les hypothèses théoriques, les formules de diffusion, la méthodologie de modélisation dans Robot Structural Analysis, ainsi qu'un exemple de calcul complet.
1. Les Hypothèses Théoriques et Méthodes de Diffusion
Avant d'ouvrir le logiciel, il est obligatoire de déterminer comment la charge de la roue (charge ponctuelle ou répartie sur une empreinte de pneu) se diffuse à une profondeur donnée $H$.
Méthode 1 : La diffusion forfaitaire pyramidale (Règles Eurocode / Fascicule 61)
C'est la méthode classique la plus répandue en ingénierie de structure pour les remblais homogènes et compactés. On considère que la charge se propage à l'intérieur d'une pyramide inversée selon une pente définie par l'angle $\alpha$.
- Sous chaussée souple/remblai courant : On utilise généralement une pente de 2:1 (vertical : horizontal), ce qui correspond à un angle de diffusion $\alpha$ tel que $\tan(\alpha) = 0,5$ ($\approx 26,6^\circ$) ou parfois une diffusion à $45^\circ$ ($\tan(\alpha) = 1$) selon la compacité du remblai et les spécifications contractuelles (Eurocode 1 Partie 2).
Méthode 2 : La théorie d'élasticité de Boussinesq
Pour des calculs plus rigoureux (milieu semi-infini, élastique, homogène et isotrope), la contrainte verticale $\sigma_z$ à une profondeur $H$ et une distance radiale $r$ induite par une charge ponctuelle $P$ s'exprime par la formule :
$$\sigma_z = \frac{3P}{2\pi H^2} \left[ \frac{1}{1 + \left(\frac{r}{H}\right)^2} \right]^{5/2}$$
2. Formules Mathématiques d'Empreinte Diffusée (Pente 2:1)
Pour l'implémentation dans Robot Structural Analysis, nous traduisons la charge d'une roue (définie par une force d'intensité $Q$ sur une surface rectangulaire initiale de dimensions $a_0 \times b_0$) en une charge répartie équivalente $q_{\text{diff}}$ à la profondeur du remblai $H$.
Calcul des dimensions de l'empreinte diffusée :
- Longueur diffusée : $a_H = a_0 + 2 \cdot H \cdot \tan(\alpha)$
- Largeur diffusée : $b_H = b_0 + 2 \cdot H \cdot \tan(\alpha)$
Si l'on adopte la diffusion classique Eurocode standard à $45^\circ$ ($\tan(\alpha) = 1$) :
$$a_H = a_0 + 2H$$
$$b_H = b_0 + 2H$$
$$a_H = a_0 + 2H$$
$$b_H = b_0 + 2H$$
Calcul de la pression équivalente sur la structure :
La contrainte verticale uniforme appliquée sur la dalle de l'ouvrage est :
$$q_{\text{diff}} = \frac{Q}{a_H \cdot b_H}$$
$$q_{\text{diff}} = \frac{Q}{a_H \cdot b_H}$$
Note importante : Si les empreintes de deux roues adjacentes (par exemple l'essieu d'un camion) se chevauchent à la profondeur $H$, il faut fusionner les surfaces et diviser la somme des charges par la surface totale combinée.
3. Méthodologie de Modélisation dans Robot Structural Analysis
Puisque le remblai n'est pas modélisé volumétriquement (pour éviter des temps de calcul prohibitifs et des modélisations de contact complexes), voici les deux approches recommandées pour appliquer cette charge diffuse dans Robot Structural Analysis :
Approche A : La méthode des "Lignes de Contour" (Recommandée)
Cette méthode est idéale si vous simulez des charges fixes ou si vous connaissez l'emplacement exact le plus défavorable.
- Générer le modèle de structure : Créez votre dalle ou élément de structure (panneau EF avec son épaisseur et ses appuis).
- Définir un plan de contour : Allez dans le menu supérieur :
Chargements>Définir charges...> OngletSurface. - Cliquez sur "Charge surfacique uniforme sur contour".
- Saisissez la valeur calculée de la pression $q_{\text{diff}}$ (valeur négative selon l'axe Z, ex:
-12.5 kN/m²). - Cliquez sur "Définition du contour". Choisissez une forme rectangulaire ou polygonale.
- Entrez les coordonnées spatiales correspondant à votre rectangle de diffusion ($a_H \times b_H$) centré sur la zone d'impact de l'ouvrage.
Approche B : Utilisation des "Bardages" (Cladding)
Si vous souhaitez automatiser la répartition sans modifier la structure des mailles EF de votre panneau structurel :
- Créez un Bardage (
Structure>Bardages) positionné exactement au niveau supérieur de votre remblai fictif, ou directement sur votre dalle. - Appliquez la charge de roue réelle ($a_0 \times b_0$) sur ce bardage.
- Configurez la répartition du bardage en mode bidirectionnel vers la dalle inférieure. Attention : Cette technique ne simule pas l'élargissement de l'empreinte géométrique. Elle est donc à réserver uniquement si vous réalisez manuellement l'expansion des dimensions du bardage aux cotes $a_H \times b_H$.
4. Exemple Pratique d'Application
Énoncé du problème :
Soit un cadre fermé (cadre en béton armé de type passage inférieur) sous une route.
- Épaisseur du remblai ($H$) : $1,50\text{ m}$
- Type de charge : Roue lourde du système de charge d'un convoi de type Eurocode LM1 (Tandem).
- Charge de la roue ($Q$) : $150\text{ kN}$
- Empreinte réelle au sol de la roue : $a_0 = 0,40\text{ m}$ (sens de circulation) et $b_0 = 0,40\text{ m}$.
- Hypothèse de diffusion : Pente de $45^\circ$ ($\tan(\alpha) = 1,0$).
Étape 1 : Calcul analytique de la diffusion
- Calcul des dimensions de la zone d'impact à la profondeur $1,50\text{ m}$ :
$$a_H = 0,40 + (2 \times 1,50 \times 1,0) = 0,40 + 3,00 = 3,40\text{ m}$$
$$b_H = 0,40 + (2 \times 1,50 \times 1,0) = 0,40 + 3,00 = 3,40\text{ m}$$ L'empreinte initiale de $0,40 \times 0,40\text{ m}$ s'est donc élargie pour devenir une zone de $3,40 \times 3,40\text{ m}$. - Calcul de la pression répartie uniforme ($q_{\text{diff}}$) :
$$\text{Surface diffusée } A_H = 3,40 \times 3,40 = 11,56\text{ m}^2$$
$$q_{\text{diff}} = \frac{150\text{ kN}}{11,56\text{ m}^2} \approx 12,98\text{ kN/m}^2$$
Étape 2 : Entrée des données pas à pas dans Robot Structural Analysis
Pour appliquer cette charge sur la dalle supérieure de notre ouvrage dans RSA :
- Création du cas de charge :
- Allez dans
Chargements>Cas de charge. - Créez un cas d'exploitation nommé
Exploitation_Tandem_Diffusé.
- Allez dans
- Paramétrage de la charge surfacique :
- Ouvrez la boîte de dialogue
Définir charges. - Sélectionnez l'onglet
Surfacepuis cliquez surCharge planaire uniforme sur contour. - Dans le champ PZ, saisissez
-12.98(l'unité par défaut est le $kN/m^2$).
- Ouvrez la boîte de dialogue
- Définition géométrique du contour :
- Cliquez sur le bouton
Définition du contour. - Si le centre de votre ouvrage se trouve à l'origine $(X=0, Y=0)$, les 4 sommets de votre rectangle de diffusion seront :
- Point 1 : $X = -1.70, Y = -1.70$
- Point 2 : $X = 1.70, Y = -1.70$
- Point 3 : $X = 1.70, Y = 1.70$
- Point 4 : $X = -1.70, Y = 1.70$
- Cliquez sur
Appliquer.
- Cliquez sur le bouton
- Maillage et Analyse :
- Lancez le maillage automatique du panneau (
Outils>Maillage>Génération du maillage EF). - Robot va automatiquement découper la charge sur le contour et l'affecter aux nœuds des éléments finis de la dalle.
- Lancez les calculs via l'icône Calculer (calculs statiques linéaires).
- Lancez le maillage automatique du panneau (
5. Synthèse comparative des sollicitations
Pour illustrer l'importance de ce calcul, voici l'impact sur les résultats de moments fléchissants ($M_{xx}$) au centre de la dalle d'épaisseur $30\text{ cm}$ de notre ouvrage testé sous RSA :
| Paramètre mesuré | Sans diffusion (Pression directe au sol) | Avec diffusion à $1,50\text{ m}$(Notre calcul) |
|---|---|---|
| Surface d'application | $0,40 \times 0,40\text{ m}$ | $3,40 \times 3,40\text{ m}$ |
| Intensité de la pression | $937,5\text{ kN/m}^2$ | $12,98\text{ kN/m}^2$ |
| Moment Fléchissant Max | $\approx 42,5\text{ kNm/m}$ (Pic aigu localisé) | $\approx 14,2\text{ kNm/m}$(Répartition douce) |
Le calcul de diffusion évite ainsi un surdimensionnement majeur des armatures en flexion et élimine les faux problèmes de poinçonnement sur les dalles d'ouvrages enterrés.
Conclusion et Bonnes Pratiques
Le calcul de la diffusion à travers le remblai est une étape indispensable pour obtenir des notes de calcul optimisées et conformes aux exigences réglementaires. Lors de vos modélisations sur Robot Structural Analysis, retenez ces trois règles d'or :
- Validez toujours l'angle de diffusion réglementaire avec le bureau d'études géotechniques.
- Prenez garde au recouvrement des empreintes pour les ponts-cadres étroits.
- Utilisez l'option des charges sur contour pour garder un modèle propre sans surcharger le maillage EF géométrique de la structure.
Pour affiner votre note de calcul, souhaitez-vous :
- Un script d'automatisation ou un tableau Excel de calcul des empreintes pour plusieurs profondeurs ?
- La procédure pour créer des charges roulantes dynamiques (convois mobiles) associées à cette diffusion dans RSA ?
- La prise en compte du coefficient de majoration dynamique (coefficient d'amplification) selon l'Eurocode ?

