Approche numérique du calcul à la rupture :

application au dimensionnement des structures

Éléments finis pour le génie civil

Jérémy Bleyer ([email protected])

Séance 10 (2022-2023)

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 1 / 33

 Calcul à la rupture

Sommaire

1 Calcul à la rupture

2 Applications

3 Optimisation topologique

4 Formulation à deux phases et méthode des bielles-et-tirants

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 1 / 33

 Calcul à la rupture

Méthodes actuelles de dimensionnement des structures

Dimensionnement des structures aux états limites ultimes (ELU)

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 2 / 33

 Calcul à la rupture

Méthodes actuelles de dimensionnement des structures

Dimensionnement des structures aux états limites ultimes (ELU)

Limites des méthodes actuelles

réglementation (Eurocodes) : approche "semi-empirique", coefficients de sécurité,
marge de sécurité non quantifiable

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 2 / 33

 Calcul à la rupture

Méthodes actuelles de dimensionnement des structures

Dimensionnement des structures aux états limites ultimes (ELU)

Limites des méthodes actuelles

réglementation (Eurocodes) : approche "semi-empirique", coefficients de sécurité,
marge de sécurité non quantifiable
calcul complet pas-à-pas : temps de calcul très longs, convergence
en bureau d’études, 98% des calculs sont élastiques !
Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 2 / 33
 Calcul à la rupture

Besoins d’un outil de dimensionnement fiable et performant

Structures complexes, mixtes, impression 3D

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 3 / 33

 Calcul à la rupture

Besoins d’un outil de dimensionnement fiable et performant

Structures complexes, mixtes, impression 3D

situations complexes pas couvertes par les codes, impact environnemental

⇒ besoin d’outils robustes et efficaces

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 3 / 33
 Calcul à la rupture

Théorie du calcul à la rupture [J. Salençon, 1983]

Une alternative aux approches traditionnelles du calcul des ouvrages

approche directe rigoureuse : compatibilité équilibre + critère de résistance
encadrement de la charge de ruine : bornes inférieure et supérieure
visualisation des mécanismes de ruine
milieux continus généralisés (poutres, plaques, coques)

approche cinématique

charge limite

approche statique

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 4 / 33

 Calcul à la rupture

Théorie du calcul à la rupture

Chargement de référence : F => charge ultime (facteur de sécurité) Flim = λ+ F ?

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 5 / 33

 Calcul à la rupture

Théorie du calcul à la rupture

Chargement de référence : F => charge ultime (facteur de sécurité) Flim = λ+ F ?

Approche statique par l’intérieur

λ+ ≥ λs = max λ
λ,σ
σ S.A. avec λF (équilibre)
t.q.
f (σ(x)) ≤ 0 ∀x ∈ Ω (critère de résistance)

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 5 / 33

 Calcul à la rupture

Théorie du calcul à la rupture

Chargement de référence : F => charge ultime (facteur de sécurité) Flim = λ+ F ?

Approche statique par l’intérieur

λ+ ≥ λs = max λ
λ,σ
σ S.A. avec λF (équilibre)
t.q.
f (σ(x)) ≤ 0 ∀x ∈ Ω (critère de résistance)

Approche cinématique par l’extérieur

Principe des puissances virtuelles + fonction d’appui du critère de résistance :
π (d) = sup σ·d
σ tq f (σ)≤0 Z Z
λ+ ≤ λk = min Prm (U ) = π (∇s U )dΩ + π (n; [[U ]])dS
U C.A. Ω Γ
t.q. Pext (U ) = F · U = 1

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 5 / 33

 Calcul à la rupture

Théorie du calcul à la rupture

Chargement de référence : F => charge ultime (facteur de sécurité) Flim = λ+ F ?

Approche statique par l’intérieur

λ+ ≥ λs = max λ
λ,σ
σ S.A. avec λF (équilibre)
t.q.
f (σ(x)) ≤ 0 ∀x ∈ Ω (critère de résistance)

Approche cinématique par l’extérieur

Principe des puissances virtuelles + fonction d’appui du critère de résistance :
π (d) = sup σ·d
σ tq f (σ)≤0 Z Z
λ+ ≤ λk = min Prm (U ) = π (∇s U )dΩ + π (n; [[U ]])dS
U C.A. Ω Γ
t.q. Pext (U ) = F · U = 1

λs ≤ λ + ≤ λk

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 5 / 33

 Calcul à la rupture

Formulation discrète des problèmes de calcul à la rupture

Approche statique par l’intérieur :

λ+ = max λ
λ,σ
div σ + λf = 0
σ · n = 0 sur ∂ΩT
t.q.
[[σ]] · n = 0 à travers Γ
f ( σ ) ≤ 0 ∀x ∈ Ω

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 6 / 33

 Calcul à la rupture

Formulation discrète des problèmes de calcul à la rupture

Approche statique par l’intérieur :

λ+ ≥ λs = max λ
λ,σ
Hσ + λF = 0
t.q. Cσ = 0
f ( σ i ) ≤ 0 ∀i

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 6 / 33

 Calcul à la rupture

Formulation discrète des problèmes de calcul à la rupture

Approche statique par l’intérieur :

λ+ ≥ λs = max λ
λ,σ
Hσ + λF = 0
t.q. Cσ = 0
f ( σ i ) ≤ 0 ∀i

Le critère f (σ i ) est convexe + contraintes d’égalité linéaires

⇒ problème d’optimisation convexe

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 6 / 33

 Calcul à la rupture

Formulation discrète des problèmes de calcul à la rupture

Approche statique par l’intérieur :

λ+ ≥ λs = max λ
λ,σ
Hσ + λF = 0
t.q. Cσ = 0
f ( σ i ) ≤ 0 ∀i

Le critère f (σ i ) est convexe + contraintes d’égalité linéaires

⇒ problème d’optimisation convexe

Approche cinématique par l’extérieur :

Z Z
λ+ = min π (d)dΩ + π (n, V )dS
U ZΩ Γ
t.q. f · UdΩ = 1
Ω
d = (∇U + ∇UT )/2
V = [[U ]]

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 6 / 33

 Calcul à la rupture

Formulation discrète des problèmes de calcul à la rupture

Approche statique par l’intérieur :

λ+ ≥ λs = max λ
λ,σ
Hσ + λF = 0
t.q. Cσ = 0
f ( σ i ) ≤ 0 ∀i

Le critère f (σ i ) est convexe + contraintes d’égalité linéaires

⇒ problème d’optimisation convexe

Approche cinématique par l’extérieur :

λ+ ≤ λk = min
U,d,V
∑ ωi π (di ) + ∑ ωj0 π (n; Vj )
i j
t.q. FT U = 1
d = Bu
V = Du

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 6 / 33

 Calcul à la rupture

Formulation discrète des problèmes de calcul à la rupture

Approche statique par l’intérieur :

λ+ ≥ λs = max λ
λ,σ
Hσ + λF = 0
t.q. Cσ = 0
f ( σ i ) ≤ 0 ∀i

Le critère f (σ i ) est convexe + contraintes d’égalité linéaires

⇒ problème d’optimisation convexe

Approche cinématique par l’extérieur :

λ+ ≤ λk = min
U,d,V
∑ ωi π (di ) + ∑ ωj0 π (n; Vj )
i j
t.q. FT U = 1
d = Bu
V = Du

Les fonctions π (di ), π (n; Vj ) sont convexes + contraintes d’égalité linéaires

⇒ problème d’optimisation convexe
Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 6 / 33
 Calcul à la rupture

Logiciels
OptumG2 (http://www.optumce.com)

OptumG3 (version 3D)

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 7 / 33

 Calcul à la rupture

Logiciels
DS-Steel pour assemblages métalliques : société Strains (http://www.strains.fr/)

En développement : applications béton (structures massives), maçonnerie...

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 8 / 33

 Applications

Sommaire

1 Calcul à la rupture

2 Applications

3 Optimisation topologique

4 Formulation à deux phases et méthode des bielles-et-tirants

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 8 / 33

 Applications

Vérification d’assemblages métalliques 3D [Thèse de C. El Boustani]

Platine de base pour poteau avec boulons

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 9 / 33

 Applications

Vérification d’assemblages métalliques 3D [Thèse de C. El Boustani]

Dimensionnement à l’Eurocode 3 : vérification manuelle de 20 mécanismes de ruine (cf.

méthode des lignes de rupture)

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 10 / 33

 Applications

Résistance de la platine [Thèse de C. El Boustani]

Approches cinématique : éléments finis 3D + contact unilatéral platine/boulons

El Boustani, Bleyer, Arquier, Ferradi, Sab, Eng. Str., in press

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 11 / 33
 Applications

Structures massives en béton armé [Thèse de H. Vincent]

approche "manuelle" : modèles bielles et tirants ⇒ très difficile en 3D

vérification automatique des structures massives en béton armé
aide à l’optimisation du ferraillage

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 12 / 33

 Applications

Critère de résistance 3D pour le béton

béton = critère de résistance de Rankine
(
Σ  ft I
−fc ≤ σI,II,III ≤ ft ⇔
− fc I  Σ

Vincent, Arquier, Bleyer, de Buhan, Int. J. Num. Anal. Meth. Geomech., 2018
Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 13 / 33
 Applications

Résistance au poinçonnement et à la flexion d’une fondation

circulaire

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 14 / 33

 Applications

Résistance au poinçonnement et à la flexion d’une fondation

circulaire

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 14 / 33

 Applications

Application au viaduc d’Echinghen

Modélisation d’un bloc déviateur de voussoir précontraint (viaduc d’Echinghen)

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 15 / 33

 Applications

Application au viaduc d’Echinghen

structure en béton armé et câbles de précontrainte

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 15 / 33

 Applications

Application au viaduc d’Echinghen

Vérification de résistance (approche statique) et post-traitement sous forme d’un schéma

bielles et tirants 3D

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 15 / 33

 Optimisation topologique

Sommaire

1 Calcul à la rupture

2 Applications

3 Optimisation topologique

4 Formulation à deux phases et méthode des bielles-et-tirants

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 15 / 33

 Optimisation topologique

Optimisation topologique

Trouver Ω ⊆ D qui minimise de la souplesse à volume donné :

[TopOpt in Python, DTU]

Formulation 1 :
Z
min T · u dS
Z
Ω,u ∂Ω min σ : C : σ dx
Ω,σ Ω
s.t. σ = C : ε[u] dans Ω ⇔
σ SA dans Ω
σ SA
|Ω| ≤ η |D|
|Ω| ≤ η |D|

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 16 / 33

 Optimisation topologique

Optimisation topologique

Trouver Ω ⊆ D qui minimise de la souplesse à volume donné :

[TopOpt in Python, DTU]

Formulation 2 :
Ω remplacé par une fonction binaire ρ(x) = 0 ou 1 :
Z
min T · u dS
Z
ρ,u ∂D min ρ−1 σ : C−1 : σ dx
s.t. σ = ρC : ε[u] dans D ⇔
ρ,σ D
R SA dans D
σ
R SA
σ
D ρ dx ≤ η |D|
D ρ dx ≤ η |D| ρ(x) ∈ {0, 1}
ρ(x) ∈ {0, 1}

⇒ la contrainte binaire rend le problème très dur à résoudre en pratique

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 16 / 33

 Optimisation topologique

Optimisation topologique

Trouver Ω ⊆ D qui minimise de la souplesse à volume donné :

[TopOpt in Python, DTU]

Formulation 3 :
Ω remplacé par une fonction continue ρ(x) ∈ [0; 1] :
Z
min ρ−1 σ : C−1 : σ dx
ρ,σ D
R SA dans D
σ
D ρ dx ≤ η |D|
0≤ρ≤1

⇒ problème convexe ! facile mais produit des niveaux de gris 0 < ρ < 1

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 16 / 33

 Optimisation topologique

Optimisation topologique

Trouver Ω ⊆ D qui minimise de la souplesse à volume donné :

[TopOpt in Python, DTU]

Formulation 4 :
densité intermédiaires pénalisées, classiquement(SIMP) : C(ρ) = ρp C, p>1
Z
min ρ−p σ : C−1 : σ dx
ρ,σ D
R SA dans D
σ
D ρ dx ≤ η |D|
0≤ρ≤1

⇒ problème non-convexe, procédure itérative

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 16 / 33

 Optimisation topologique

Exemple
Exemple FEniCS (https://comet-fenics.readthedocs.io)

Améliorations
filtrage pour éliminer la dépendence au maillage, etc.
autres fonctions coût
contraintes supplémentaires : limite d’élasticité

Autre question
Trouver Ω ⊆ D qui maximise la résistance à volume donné ?

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 17 / 33

 Optimisation topologique

Maximisation de la capacité portante [Thèse de L. Mourad]

Objectif : trouver la structure Ω possédant une charge limite maximale pour un volume
donné

λ+ (η ) = max λ
λ,σ,Ω
s.t. σ en éq. avec λF dans Ω
σ∈G dans Ω
|Ω| ≤ η |D|

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 18 / 33

 Optimisation topologique

Maximisation de la capacité portante [Thèse de L. Mourad]

Objectif : trouver la structure Ω possédant une charge limite maximale pour un volume
donné

λ+ (η ) = max λ
λ,σ,Ω
s.t. σ en éq. avec λF dans D
σ∈G dans Ω
σ=0 dans D \ Ω
|Ω| ≤ η |D|
extension par σ = 0 en dehors de Ω

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 18 / 33

 Optimisation topologique

Maximisation de la capacité portante [Thèse de L. Mourad]

Objectif : trouver la structure Ω possédant une charge limite maximale pour un volume
donné

λ+ (η ) = max λ
λ,σ,ρ
s.t. σ en éq. avec λF dans D
R ∈ ρG
σ dans D
D ρ dx ≤ η |D|
ρ ∈ {0; 1}
fonction caractéristique ρ = 0 ou ρ = 1 de Ω

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 18 / 33

 Optimisation topologique

Maximisation de la capacité portante [Thèse de L. Mourad]

Objectif : trouver la structure Ω possédant une charge limite maximale pour un volume
donné

λ+ (η ) = max λ
λ,σ,ρ
s.t. σ en éq. avec λF dans D
(LOAD-MAX)
R ∈ ρG
σ dans D
D ρ dx ≤ η |D|
ρ ∈ [0; 1]
convexification du problème

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 18 / 33

 Optimisation topologique

Maximisation de la capacité portante [Thèse de L. Mourad]

Objectif : trouver la structure Ω possédant une charge limite maximale pour un volume
donné

λ+ (η ) = max λ
λ,σ,ρ
s.t. σ en éq. avec λF dans D
(LOAD-MAX)
R ∈ ρG
σ dans D
D ρ dx ≤ η |D|
ρ ∈ [0; 1]
convexification du problème

Propriétés :
la contrainte σ ∈ ρG est convexe
problème d’analyse limite avec un champ scalaire ρ supplémentaire
λ+ (η = 1) : charge limite du domaine D

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 18 / 33

 Optimisation topologique

Formulation alternative : minimisation du volume

Objectif : trouver la structure Ω de volume de matière minimal supportant un

chargement donné
1
Z
η − (λ) = min ρ dx
ρ,σ,Ω |D| D
s.t. σ en éq. avec λF dans D (VOL-MIN)
σ ∈ ρG dans D
ρ ∈ [0; 1]

convexification du problème

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 19 / 33

 Optimisation topologique

Choix du critère de résistance

(
+ft eα ⊗ eα
acier : von Mises, béton : Rankine, on veut en fait σ ≈
−fc eα ⊗ eα

Proposition : critère L1 −Rankine

afin de favoriser les efforts uni-axiaux (modèle "treillis")
 : 
σJ σJ
σ ∈ GL1-Rankine ⇔ ∑ max − ; ≤1
J =I,II,III
fc ft

(a) résistances symétriques ft = fc (b) résistances asymétriques fc /ft = 5

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 20 / 33

 Optimisation topologique

Poutre MBB

Résultats
(VOL-MIN) = (LOAD-MAX)
On économise jusqu’à 50% de matière

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 21 / 33

 Optimisation topologique

Poutre MBB (après pénalisation)

pmax = 3, 20 itérations

(a) L1 -Rankine (b) Rankine (c) von Mises

(LOAD-MAX) pénalisé

(a) L1 -Rankine (b) Rankine (c) von Mises

(VOL-MIN) pénalisé

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 22 / 33

 Optimisation topologique

Matériaux à résistances asymétriques

1 0.5

(a) fc /ft = 10 (b) fc /ft = 0.1

Champs de contraintes principales du matériau : contraintes de traction et contraintes

de compression

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 23 / 33

 Optimisation topologique

Matériau sans résistance à la traction

ex : structures en pierres ou maçonnerie

1.0

1.5
1.0

s 1.4 s
1.8 1.8

(a) Avant pénalisation s=0 (b) Après pénalisation

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 24 / 33

 Optimisation topologique

Matériau sans résistance à la traction

ex : structures en pierres ou maçonnerie

1.0

1.5
1.0

s 1.4 s
1.8 1.8

(a) Façade de la passion s=0 (b) Après pénalisation

forme d’arc !
Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 24 / 33
 Optimisation topologique

Matériau sans résistance à la traction

ex : structures en pierres ou maçonnerie

1.0

1.5
1.0

s 1.4 s
1.8 1.8

(a) Avant pénalisation s = 0.3 (b) Après pénalisation

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 24 / 33

 Optimisation topologique

Matériau sans résistance à la traction

ex : structures en pierres ou maçonnerie

1.0

1.5
1.0

s 1.4 s
1.8 1.8

(a) Avant pénalisation s = 0.3 (b) Après pénalisation

création d’arc-boutants !

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 24 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Sommaire

1 Calcul à la rupture

2 Applications

3 Optimisation topologique

4 Formulation à deux phases et méthode des bielles-et-tirants

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 24 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Optimisation de forme à deux phases

on souhaite optimiser deux phases de façon indépendante

(+ vide)

Critère de résistance

1 2
σ = σ + σ

1 2 1 1
σ ∈ G(ρ1 , ρ2 ) ⇔ ∃σ , σ s.t. σ ∈ ρ1 G
 2
σ ∈ ρ 2 G2


avec ρ1 + ρ2 ≤ 1

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 25 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Formulation du problème

Maximisation de capacité portante Minimisation de volume

λ+ = 1
Z
max λ η − = min c(ρ1 , ρ2 ) dx
λ,σ 1 ,σ 2 ,ρ 1 ,ρ2
σ ,σ ,ρ1 ,ρ2 |D| D
1 2
s.t. (σ 1 + σ 2 ) en éq.avec λF
σ 1 ∈ ρ1 G1 s.t. (σ 1 + σ 2 ) en éq.avec λF
2 2 σ 1 ∈ ρ 1 G1
R ∈ ρ2 G
σ
σ 2 ∈ ρ 2 G2
D c(ρ1 , ρ2 ) dx ≤ η |D|
0 ≤ ρ1 ≤ 1 0 ≤ ρ1 ≤ 1
0 ≤ ρ2 ≤ 1 0 ≤ ρ2 ≤ 1
ρ1 + ρ2 ≤ 1 ρ1 + ρ2 ≤ 1
(BIMAT-LOAD-MAX) (BIMAT-VOL-MIN)

c(ρ1 , ρ2 ) = ρ1 + ρ2 : fonction de coût de matière

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 26 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Poutre MBB

critère L1 -Rankine en compression pure et traction pure

(a) η = 0.05 (b) η = 0.1

(c) η = 0.2
(d) η = 0.3

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 27 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Poutre MBB

Phase en traction : renforts selon différentes orientations

(a) Orientations selon 0◦ (b) Orientations selon 0◦ and ±30◦

(c) Orientations selon 0◦ et ±45◦ (d) Orientations selon 0◦ et 90◦

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 28 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Structure complexe

Figure – Poutre-voile d’épaisseur 25 cm étudié

Figure – Modèle BT proposé dans [(Muttoni,
dans [(Muttoni, 2015)]
2015)]

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 29 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Structure complexe

Figure – BIMAT-VOL-MIN

Figure – Modèle BT proposé dans [(Muttoni,

2015)]

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 29 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Structure complexe

Figure – BIMAT-VOL-MIN, orthotrope

Figure – Modèle BT avec renforcements
orthogonaux proposé dans [(Muttoni, 2015)]

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 29 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Exemple du corbeau

Figure – Structure du corbeau

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 30 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Exemple du corbeau

Figure – Modèle BT proposé par [(Bosc, 2008)]

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 30 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Exemple du corbeau

(a) Aciers isotropes (b) Aciers orthogonaux

Figure – Modèles optimisés avec (BIMAT-VOL-MIN)

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 31 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Exemple du corbeau

équilibre vérifié dans les trois cas

différence de l’ordre de -3% avec la méthode B&T

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 31 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Exemple du poutre-poteau

Figure – La géométrie du poutre-poteau étudié

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 32 / 33

 Formulation à deux phases et méthode des bielles-et-tirants

Exemple du poutre-poteau

Figure – Modèles optimisés par (BIMAT-VOL-MIN)

Figure – Schéma illustrant le dispositif de barres relevées

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 32 / 33
 Formulation à deux phases et méthode des bielles-et-tirants

Extensions

milieux généralisés : plaques et coques (multicouches), milieux continus à rotation

interne (Cosserat)
conditions d’incendie : dégradation des propriétés de résistance + changements de
géométrie (CSTB)
optimisation robuste : prise en compte d’incertitudes ⇒ charge limite la plus
défavorable

nominal
robust optimization

Éléments finis pour le génie civil Calcul à la rupture Séance 10 (2022-2023) 33 / 33