A12-Comportement Au FEU Des Ouvrages Et Outils Disponibles PDF
A12-Comportement Au FEU Des Ouvrages Et Outils Disponibles PDF
A12-Comportement Au FEU Des Ouvrages Et Outils Disponibles PDF
Octobre 2012
1 / 122
SOMMAIRE
INTRODUCTION PN ISI ACTION N 12 ...................................................................................... 4
I. MTHODES ET OUTILS DE VRIFICATION DES STRUCTURES EN BTON SOUMISES
AU FEU. ........................................................................................................................................ 7
I.1. INTRODUCTION : .................................................................................................................. 9
I.2. MTHODES DE VRIFICATION DES STRUCTURES EN BTON EN SITUATION
DINCENDIE ................................................................................................................................ 10
I 2.1.
I 2.2.
I 2.3.
I 2.4.
INTRODUCTION ............................................................................................................ 72
III.2.
PROPRITS AUX TEMPRATURES LEVES DES ACIERS DE CONSTRUCTION
72
III 2.1 PROPRITS DE LACIER AU CARBONE ................................................................................ 72
III 2.2 PROPRITES DES ACIERS INOXYDABLES ............................................................................. 74
III.3.
CALCUL DE LECHAUFFEMENT DES STRUCTURES EN ACIER ET MIXTE ACIER
BTON 75
III 3.1 MTHODES DE CALCUL SIMPLIFIES ................................................................................... 75
2 / 122
CONCLUSION ............................................................................................................... 87
ANNEXE A
A.1
A.2
A.3
A.4
A.5
A.6
A.7
A.8
A.9
A.10
A.11
A.12
A.13
A.14
A.15
IV. OUTILS DE CALCUL DE RSISTANCE AU FEU DES STRUCTURES EN VERRE ....... 115
IV 1.
IV 1.1.
IV 1.2.
IV 1.3.
IV 2.
FEU
IV 2.1.
IV 2.2.
3 / 122
L'action 12 du PN ISI est destine tablir un tat de l'art des mthodes de vrification des structures en
situation d'incendie et d'effectuer une corrlation entre les mthodes et les ouvrages viss. Les rflexions
qui ont t menes au cours de cette action s'articulent autour des cinq points suivants :
1.
2.
Classement des mthodes recenses en fonction des types douvrages viss et des objectifs de
scurit attendus.
3.
4.
Dfinition de critres de choix des mthodes, prise en compte des systmes constructifs dans
lapprciation des critres.
5.
Proposition dune grille de correspondance entre type douvrage et type danalyse. Avantages et
inconvnients. Niveaux de simplification attendus.
Lanalyse du comportement mcanique des structures exposes au feu peut tre ralis par les trois
approches suivantes
analyse par lment, dans laquelle chacun des lments de la structure sera vrifi en le
considrant totalement spar des autres lments. Les liaisons tant remplaces par des
conditions aux limites appropries ;
analyse de parties de la structure, dans laquelle une partie de la structure sera directement prise
en compte en utilisant des conditions aux limites appropries pour reprsenter sa liaison avec le
reste de la structure ;
analyse de structure globale, dans laquelle la structure totale sera utilise dans le calcul
Par rapport ces procdures danalyse relatives au comportement mcanique des structures en situation
dincendie, les remarques suivantes peuvent tre faites :
Lanalyse par lment sapplique aux lments de structure isols (lment par lment), donc
facile utiliser, en particulier avec les mthodes de calcul simplifies et par consquent
largement utilise sous conditions dincendie normalis (par exemple : feu ISO-834) ;
Lanalyse de parties de la structure ou lanalyse globale de la structure considre au moins
plusieurs lments de la structure, si bien que leffet dinteraction entre eux ainsi que la
redistribution de charge des parties chauffes (parties affaiblies et assouplies lintrieur du
compartiment en feu) aux parties froides (parties plus rsistantes et plus raides en dehors du
compartiment en feu) peuvent tre pris en compte directement et de manire exacte. Par
consquent, lanalyse de la structure globale permet ainsi dobtenir une valuation plus raliste
du comportement au feu de la structure.
4 / 122
Le premier classement des mthodes a consist en une rpartition naturelle des mthodes selon le
matriau composant la structure. Les quatre matriaux considrs dans ce rapport sont ceux le plus
couramment rencontrs dans le domaine du btiment : le bton, le bois, lacier (et mixte acier-bton) et le
verre. Le prsent rapport est dailleurs structur en quatre parties traitant chacune dun matriau.
Bien que les mthodes utilises pour les quatre matriaux soient trs diffrentes, il existe plusieurs
champs communs. Le premier de ces champs apparat sous la forme des trois types de mthodes que lon
peut utiliser selon la prcision recherche :
Les mthodes par valeurs tabules qui indiquent directement des rsultats selon plusieurs
critres. Ces mthodes sont ncessairement enveloppe pour couvrir l'ensemble de cas
envisageables.
Les mthodes simplifies qui permettent de calculer une rsistance au feu l'aide de quelques
rgles simples. Elles sont donc plus prcises que les mthodes tabules.
Les mthodes avances qui, bass sur les proprits caractristiques des matriaux et sur la
ralit des phnomnes physiques, permettent par un calcul thermomcanique dvaluer de
manire prcise le comportement au feu des structures. Ces mthodes requirent en gnral des
outils de calcul numrique de type lments finis. A lheure actuelle, il existe de nombreux
modles numriques (valids par rapport aux rsultats dessais) qui dcrivent adquatement la
rponse mcanique des lments en acier ou mixtes acier-bton exposs au feu. Pour le bton,
la plupart des modles avancs disponibles pour le calcul mcanique ncessitent encore dtre
valids.
Les tableaux 1 et 2 rsument les diffrentes possibilits dapplication des ces trois mthodes de calcul.
Leur emploi est aujourdhui encadr par larrt du 22 mars 2004 relatif la rsistance au feu des produits,
des lments de construction et douvrage, qui mane du ministre de lintrieur.
Les mthodes de calcul simplifies (ou valeurs tabules) qui sont bases principalement sur les conditions
de lincendie normalis, permettent uniquement de vrifier la stabilit au feu des lments les plus
courants (poteau, poutre, ) vis--vis des exigences descriptives imposes en termes de rsistance au
feu dans les rglementations traitant de la scurit incendie.
Les mthodes de calculs avances sont gnralement employes en ingnierie de la scurit incendie
(prenant en compte des scnarios dincendie rels) pour permettre dvaluer de manire optimale les
mesures de protection ncessaires afin datteindre les objectifs fixs en matire de scurit (des
personnes des biens). Toutefois, leur utilisation ncessite de recourir, dans un grand nombre de
situations, un avis sur tude (devant tre ncessairement fait par un laboratoire agre).
5 / 122
Procdure de Calcul
Rglementation
descriptive
Analyse par
lment
Analyse de parties
de la structure
Analyse globale
de la structure
Rglementation
objectif
Analyse par
lment
Analyse de parties
de la structure
Analyse globale
de la structure
Calcul des
actions
mcanique et
conditions aux
limites
Slection des
actions
mcaniques
Calcul des
actions
mcanique et
conditions aux
limites
Slection des
actions
mcaniques
Valeur
tabule
Modles de
calcul
simplifis
Modles
de calcul
avancs
OUI
OUI
OUI
NON
OUI, si
disponible
OUI
NON
NON
OUI
NON
OUI, si
disponible
OUI
NON
NON
OUI
NON
NON
OUI
Tableau 1 : Domaine dapplication des procdures danalyse pour le calcul des structures
soumises au feu
Action
thermique
Approche pour la
vrification du
comportement au feu
Utilisation par
les Bureaux
d'tudes
Rglementation
descriptive
Valeurs tabules
OUI
NON
Mthodes simplifies
OUI
NON
Mthodes avances
OUI
OUI
Valeurs tabules
Non applicable
Mthodes simplifies
OUI (si
disponible)
OUI
Mthodes avances
OUI
OUI
Rglementation
objectif
6 / 122
7 / 122
Avant propos :
Pour les structures en bton, les mthodes de calcul au feu peuvent tre classes soit par degr de complexit
(mthode) soit par type douvrage (problme). La mthode la plus simple consiste utiliser des valeurs tabules mais le
champ dapplication reste limit (feu ISO et lment isol). Les mthodes simplifies permettent de couvrir dautres
types de feu et donner une analyse plus prcise que les valeurs tabules (fondes pour couvrir les cas les plus
dfavorables) mais lanalyse reste toujours limite un seul lment (poutre, poteau, dalle, voile). Pour tudier une
structure plus complexe et plus globale, on fait appel des mthodes de calculs plus sophistiques soit par la mthode
des lments finis en utilisant des lments massifs tridimensionnels (lourds en temps de calcul) soit en utilisant des
modles gnralises en simplifiant la structure sous forme dassemblage de poutres multifibres (EF 1D) et de plaques
multicouches (EF 2D).
Elments de type
poutre sans effets du
second ordre
Elment poteau ou
poutre avec effets du
second ordre
(rotule plastique,
assemblage)
Structure complte
Globale
Mthode globale
volumique
Programmes
gnraux E.F
Modle gnralis
Modle fibres
+Autres modles
gnraliss
simplifis (barres,
poutres Timoshenko ou
Bernoulli)
Modle local
Mthodes simplifies
Valeurs
tabules
Calcul de
Section (diagramme
moment/courbure)
NON
Isotherme 500(EC2
1-2)
OUI
Programmes
gnraux E.F
Modle fibres
+
Autres modles
gnraliss
simplifies (barres,
poutres Timoshenko ou
Bernoulli)
Programmes
gnraux E.F
Modle couches
Modles fibres
(dalles)
+
Autres modles
gnraliss
simplifies (plaques,
coques,)
NON
Modle couches
Modles fibres
(dalles)
NON
Modles fibres
(sous rserve)
Elments
spcifiques
(ressorts, rotules, )
Modle
adhrence
glissement
Formules simplifies
(DTU Feu-Bton)
OUI
Modle fibres +
modle couches
+ lments
spcifiques
NON
NON
NON
Programmes
gnraux E.F
Programmes
gnraux E.F
Programmes
gnraux E.F
NON
OUI
Calcul de
Section(diagramme
moment/courbure)
Isotherme 500
OUI
Mthode par zones
Annexe E (EC2 1-2)
Mthode par zone
Mthode de Faessel
Annexe B.3 (EC2 1-2)
8 / 122
OUI
I.1. Introduction :
Comme le montre le tableau prcdent, les mthodes de calcul des structures en bton en situation dincendie sont de
plus en plus nombreuses grce lavance technologique dans le domaine de linformatique et les mthodes de calcul
numriques mais aussi grce au dveloppement de modles thoriques de plus en plus sophistiqus. Il est possible
donc de dimensionner un lment en bton pour une exigence descriptive (feu normalis) en utilisant des mthodes
tabules bases sur des approches exprimentales et des hypothses gnralement conservatives qui nous placent du
ct de la scurit. Mais, il est aussi intressant de pouvoir raffiner le dimensionnement ou de llargir un feu naturel
(ingnierie de la scurit incendie) par des mthodes plus pousses se basant sur une approche thorique valide par
des expriences.
On a le choix entre diffrentes mthodes suivant le degr de prcision quon souhaite avoir et le domaine dapplicabilit
de ces dernires : une mthode simplifie est, certes, simple utiliser mais son domaine dapplication est gnralement
limit. On a recours alors des modles encore plus sophistiqus qui ont des domaines dapplication de plus en plus
larges : les mthodes numriques avances ou encore les mthodes volumiques globales se basant sur des logiciels
de calcul trs puissants.
Le but de ce chapitre est de prsenter les diffrentes mthodes et outils de calculs quon a pu trouver dans le domaine
du dimensionnement au feu des structures en bton tout en progressant des mthodes simples aux mthodes les plus
complexes et tout en prcisant chaque fois les limites et le domaine dapplicabilit. Lorsque cest possible, on fait un
classement de ces mthodes par type dlment (poutre, poteau, dalle, voile,). Cela est de moins en moins vrai avec
des mthodes de plus en plus avances qui sont applicables avec tous les types dlments.
9 / 122
fi =
Gk + fi Qk ,1
G Gk + Q ,1Qk ,1
am =
A f a
A f
si
si
ki
ki
O
fki=fyki pour les aciers de bton arm et fki=fpki pour les aciers de prcontrainte.
10 / 122
I 2.1.1. Poutres
La rsistance au feu de poutres prfabriques en bton arm et prcontraint est gnralement dtermine partir de
valeurs tabules. Pour ce faire, il faut remplir les conditions suivantes:
Les poutres sont normalement exposes sur trois faces (la face suprieure est protge par des lments de
plancher ou de toiture.)
Les tableaux sont valables pour les sections reprises la Figure I-2.
d eff
a eff = a1.85
bmin
bw
b
min
et la rgle
Pour les armatures situes dans les angles des poutres un lit d'armatures, la distance de l'axe la face latrale de
la poutre doit tre augmente de 10 mm pour les poutres de largeurs jusqu' celles dfinies dans la colonne (4) du
Tableau 1
Les ouvertures dans l'me n'ont pas d'influence sur la rsistance au feu pour autant que la section subsistante de
l'lment dans la zone de traction ne soit pas infrieure Ac = 2b min o bmin est donn dans le Tableau 1.
11 / 122
Le Tableau 1 indique les valeurs minimales des distances des axes des armatures aux parements infrieurs et latraux
de poutres sur appuis simples sans moments sur appuis, ainsi que les valeurs minimales des largeurs de poutres et de
leurs mes. Ce tableau provient de l'Eurocode 2 partie 1-2 - tableau 5.5.
Tableau1 : Dimensions et distances de l'axe des armatures au parement pour les poutres sur appuis simples
sans moment sur appuis en bton arm et prcontraint
I 2.1.2. Poteaux
La rsistance au feu des poteaux dpend de diffrents paramtres:
Taille et lancement des poteaux
Niveau de chargement
Excentricit du 1er ordre
Rsistance du bton
Armatures
Distance de l'axe des armatures la surface
De nombreux essais au feu ont t effectus sur des poteaux en bton dans le pass. Il existe maintenant des
programmes informatiques complexes qui permettent de dterminer la rsistance au feu de poteaux arms et
prcontraints, qui tiennent compte de tous les paramtres cits ci-dessus, y compris le flambement. Les rgles de
conception suivantes sont limites aux constructions contreventes.
Selon l'Eurocode 2 partie 1-2 la rsistance au feu des poteaux sous feu ISO peut tre dtermine selon deux
mthodes dites de valeurs tabules .
12 / 122
La mthode A
Elle est base sur l'quation suivante:
O :
As f yd
Ac f cd
cc est le facteur de rsistance la compression, qui tient compte des effets long terme sur la rsistance mcanique.
La valeur pour les actions normales est entre 0,8 et 1,0.
La valeur prendre en compte ici est la valeur froid soit 0,85.
En situation dincendie, la rsistance de calcul la compression du bton est
fcd,fi = cc,fi . fck/c,fi =1,00.fck/1,00 = fck comparer avec fcd = 0,85.fck/c=0.56 fck froid.
L'Eurocode 2 partie 1-2 donne au Tableau 5.2a les dimensions minimales conservatives ( =0 et 0,fi =3m ) pour la
largeur des poteaux et la distance de l'axe des armatures au parement pour des poteaux de longueur 0,fi 3m.
La mthode A est applicable pour une excentricit maximale du premier ordre entre 0,15h 0,40h.
La mthode B
Elle est galement uniquement valable pour des constructions contreventes. Elle est base sur des calculs empiriques
prenant en compte les effets du second ordre et la rduction des caractristiques mcaniques du bton et des aciers en
fonction du champ de temprature.
Elle est valable sous les conditions suivantes :
excentricit du premier ordre e/b 0,25 avec emax = 100 mm.
lancement fi = 0,fi /i 30
Elle dpend du :
niveau de chargement n = N0Ed,fi/0.7(Ac fcd+As fyd)
pourcentage d'armatures
As f yd
Ac f cd
I 2.1.3. Dalles
Lpaisseur minimale des dalles hs et la distance minimale de laxe au parement a sont donnes dans le Tableau 5.8 de
lEurocode 2 partie 1-2 :
13 / 122
I 2.1.4. Cloisons
Les cloisons ne doivent satisfaire quaux critres d'isolation thermique "I" et d'intgrit "E".
L'paisseur minimale de cloison ne peut tre infrieure aux valeurs du Tableau 3. Si on utilise des granulats calcaires,
les valeurs donnes peuvent tre rduites de 10%. Pour viter une dformation thermique excessive et les pertes
d'tanchit qui en rsultent entre le voile et la dalle, il convient que le rapport de la hauteur libre du voile lpaisseur
du voile nexcde pas 40.
Tableau 4 : Dimensions et distances de laxe des armatures au parement pour les voiles porteurs en bton
arm
14 / 122
Figure I-3 : Section rduite d'une poutre expose au feu sur trois faces
15 / 122
16 / 122
s est le coefficient partiel sur les matriaux pour les aciers, comme dfini dans lEN 1992-1-1
s,fi est le coefficient partiel sur les matriaux pour les aciers en conditions dincendie
ks() est le facteur de rduction de la rsistance de lacier
MEd est le moment appliqu pour le calcul froid selon lEN 1992-1-1
As,prov est la section dacier tendu mise en place
As,req est la section dacier tendu exige pour le calcul froid conformment lEN 1992-1-1
Le rapport As,prov / As,req doit tre limit 1.3 dans les calculs.
-a est la distance moyenne requise de laxe des aciers au parement infrieur, donne dans le Tableau 5.5, colonne 5
pour les poutres et dans le Tableau 5.8, colonne 3 pour les dalles de lEurocode2 partie 1-2.
-d est la hauteur utile de la section
Le rapport As,prov / As,req doit tre limit 1.3 dans les calculs.
MRd,fi doit tre rduit par le facteur ks(cr) ou kp(cr) lorsque la temprature des aciers suprieurs au-dessus des appuis
dpasse :
350 C pour les armatures de BA
100 C pour les armatures de BP
17 / 122
Dalles et poutres en bton arm sur appuis simples sans moments sur appuis
Aucune zone comprime nest donc soumise, dans ce cas, des hautes tempratures. La rsistance est alors
simplement une fonction de la temprature des aciers de bton arm. Le moment rsistant maximal en situation
dincendie est donn donc par :
af = As fy,T/0.85fcb
af = As fy,T/0.85fcb
18 / 122
As fy,T/fc,Tb df <0.30
Ceci revient imposer une paisseur de zone comprime rectangulaire quivalente qui ne dpasse pas les 35% de
lpaisseur effective de la section df
Cette vrification nest pas ncessaire si des armatures longitudinales ont t prvues dans la zone comprime.
19 / 122
Pour dterminer lpaisseur de la zone endommage az, on commence par calculer le facteur de rduction moyen :
Pour les poutres et les dalles, lpaisseur de la zone endommage az est calcule par :
Pour les poteaux et les voiles (pour lesquels il faut prendre les effets du second ordre), lpaisseur de la zone
endommage az est calcule par :
TM est la temprature moyenne sur la section mdiane et w est gale lpaisseur de llment sil est expos sur une
seule face et la demi-paisseur sil est expos deux faces.
Connaissant la section droite rduite, le calcul suit la mme dmarche que la mthode de lisotherme 500C.
I.2.2.2.2. Mthode de lannexe B.3 : moment ultime avec effets du second ordre
Lannexe B.3 de lEN2 -1-2 prsente une mthode de dtermination du moment ultime avec prise en compte des effets
du second ordre :
tablir le diagramme de moment-courbure pour leffort normal concomitant NEd,fi sans prise en compte des effets du
second ordre.
dterminer le moment nominal du second ordre (proportionnel la courbure avec une pente de l0/c o c et l0
dpendent des conditions dappui du poteau) :
N
20 / 122
x
l
la section mdiane :
soit e = e0 +
y
= = f
x
l
l )
Il y aura quilibre du poteau si l'on peut trouver une courbure pour laquelle :
Mtemp
Nappliqu
e=M/N
= e0 +
e=M/N
e=M/N
e = e0 + l /
e= Mtemp/Nappliqu
e0 max
Situation d'quilibre
Divergence
21 / 122
Les mthodes de calcul avances incluent un modle de rponse thermique et un modle de rponse mcanique trs
sophistiqus intgrant tout type de feu et la gomtrie globale de tout ou partie de la structure.
Eclatement
Lclatement explosif doit tre vit, ou son incidence sur les exigences de performance (R et/ou EI) doit tre prise en
compte. Lclatement explosif est improbable lorsque la teneur en eau du bton est infrieure 3 % du poids. Au
dessus de 3 %, il convient dtudier plus prcisment linfluence de la teneur en eau, du type de granulat, de la
permabilit du bton et de la vitesse dchauffement.
Pour les poutres, les dalles et les lments tendus, si la teneur en eau du bton est suprieure 3 % du poids,
linfluence de lclatement explosif sur la fonction porteuse R peut tre value en supposant une perte locale
denrobage sur une armature de bton arm ou un paquet darmatures dans la section droite puis en vrifiant la
capacit portante rduite de la section. Pour cette vrification, la temprature des autres armatures de bton arm peut
tre suppose identique celle dune section sans clatement. Cette vrification nest pas ncessaire pour les lments
structuraux dont le comportement correct vis vis de lclatement explosif a t vrifi exprimentalement ou sur
lesquels une protection complmentaire est applique et valide par des essais.
Procdure de calcul
La procdure de calcul pour dterminer la rsistance ultime dune section est dtaille sur le schma (Fig : 14)
22 / 122
Le champ de temprature donn par le calcul lments finis est transfr au programme numrique qui permet de
calculer par exemple :
Pour les dalles ou les poutres : le moment rsistant en sommant les forces dans toutes les petites sections de
bton jusqu quilibrer la somme des efforts dans les aciers (figure I-14). Ceci permet de trouver la profondeur de la
zone comprime x et par consquent le bras de levier ncessaire pour le calcul du moment ultime. La loi momentcourbure peut tre galement tablie par ce calcul de section
Pour les voiles ou les poteaux : leffort normal concomitant maximal
- Sans prise en compte des effets du second ordre : dans ce cas, un calcul de section est suffisant pour avoir la
loi moment-courbure. Leffort normal maximal est dduit connaissant lexcentricit initiale.
- Avec prise en compte des effets du second ordre : dans ce cas, la connaissance des conditions dappui et
donc de la longueur effective de llment est ncessaire pour dterminer lexcentricit finale par double intgration
de la courbure.
e0
: courbure de la section :
( y , e0 , 0 ) = e0 + 0 . y
A froid :
As s ( ( y s , e0 , 0 ) ) + b( y ) c ( ( y, e0 , 0 )) dy = N 0
As s ( ( y s , e0 , 0 ) ) y s + b( y ) c ( ( y, e0 , 0 )) y dy = M 0
Fs = As s ( ( y s , e0 , 0 ) )
A chaud :
23 / 122
As s , ( ( y s , e0, , 0, ) s ( )) + c , ( ( y, e0, , 0, ) c ( ( x, y )) dx dy = N 0
As s , ( ( y s , e0, , 0, ) s ( )) y s + c , ( ( y, e0, , 0, ) c ( ( x, y )) y dx dy = M 0
Fs , fi = As s , ( ( y s , e0, , 0, ) s ( ))
I 2.3.2. Calcul de lquilibre dun lment en bton prcontraint charg
Le calcul est ralis par lquilibre de la section, en considrant une dformation initiale des aciers de prcontrainte
froid, toutes pertes comptabilises sauf pertes lastiques :
p 0 0,9 x 0,95
f pk
Ep
A froid :
A p p ( ( y p , e0 , 0 ) + p 0 ) + b( y ) c ( ( y, e0 , 0 )) dy = N 0
A p p ( ( y p , e0 , 0 ) + p 0 ) y p + b( y ) c ( ( y, e0 , 0 )) y dy = M 0
La rsolution de ces quations permet dobtenir e0 , dformation de la section au niveau de laxe de calcul et
courbure de la section.
0 ,
Fp = Ap p ( ( y p , e0 , 0 ) + p 0 )
A chaud :
A p p , ( ( y p , e0, , 0, ) + p 0 p ( )) + c, ( ( y, e0, , 0, ) c ( ( x, y )) dx dy = N 0
A p p , ( ( y p , e0, , 0, ) + p 0 p ( )) y p + c , ( ( y, e0, , 0, ) c ( ( x, y )) y dx dy = M 0
Fp , fi = Ap p , ( ( y p , e0, , 0, ) + p 0 p ( ))
En premire approximation on peut prendre :
Fp , fi = Ap p , ( p 0 p ( ))
I 2.3.3. Etablissement du diagramme moment courbure froid ou chaud
On suit la mme procdure que dans les deux paragraphes prcdents.
Pour un effort normal N0 fix, on dtermine pour un ensemble de courbures possibles de section 0 i , les dformations
correspondantes
e0 i
par :
As s , ( ( y s , e0i , , 0i , ) s ( )) + c, ( ( y, e0i , , 0i , ) c ( ( x, y )) dx dy = N 0
24 / 122
M 0i = As s , ( ( y s , e0i , , 0i , ) s ( )) y s + c , ( ( y, e0i , , 0i , ) c ( ( x, y )) y dx dy
Le diagramme moment courbure est tabli ensuite partir de la courbe ( 0 i ,M0i)
I 2.3.4. Prise en compte des effets du second ordre pour une poutre ou une dalle
Lquilibre de la section sous un effort normal N0 effectu pour diffrentes valeurs de courbure permet dobtenir : la
courbe moment courbure (M)
A partir de la distribution des moments M(x) issue du schma statique (par ex :
x
M ( x) = p ( x l )
2
pour une
distribution uniforme de charge p) on peut attribuer pour chaque abscisse le long de llment, la courbure
(daprs le diagramme moment-courbure).
(x)
y ' ' ( x) = ( x)
y (0) = y (l ) = 0
Le calcul est alors reprendre partir de lquation dquilibre :
La distribution des moments doit alors prendre en compte linteraction flche-effort normal M(x) + y(x) N0
Le calcul est ritr plusieurs fois. Si la flche y(x) converge, leffort N0 peut tre support.
Cette mthode permet de dterminer leffort normal supportable maximal partir duquel la flche y(x) diverge au cours
des itrations.
(M)
e0 (M)
x( y ) = l 1 + y ' ( x) 2 dx
0
x(e) = e( x) dx
0
xappui =
x( y ) + x(e)
2
I 2.3.5. Prise en compte des effets du second ordre pour un poteau ou un voile
Lquilibre de la section effectu pour diffrentes valeurs de courbure permet dobtenir : la courbe moment courbure
(M).
On suppose quon a un effort normal concomitant N0 avec une excentricit initiale e0 (imperfection gomtrique ou effort
excentr) et un moment sollicitant M0.
25 / 122
la courbure
La flche y(x) est donne par double intgration des courbures partir de lquation diffrentielle :
y ' ' ( x) =
y (0) = y (l ) = 0
y ' ' ( x) =
(cas encastr)
y ( 0) = y ' ( 0) = 0
Lquilibre est itr pour dterminer la courbure (x ) sous un effort normal
(cas bi-articul) ou
N0 et un moment total : M ( x )
= M 0 + (e0 + y ( x)) N 0
(x)
).
Cette mthode permet de dterminer pour un effort normal N0 donn, lexcentricit initiale maximale e0max partir de
laquelle le calcul itratif diverge.
Elle permet galement de dterminer pour une excentricit initiale e0, leffort normal maximal supportable N0max.
Dans le cas o la teneur en eau nest pas nulle, un pic de chaleur spcifique apparat entre 100 et 115C (phnomne
dvaporation) suivi dune dcroissance linaire entre 115 et 200C.
La valeur de ce pic dpend de la teneur en eau de faon quasi-linaire :
26 / 122
La conductivit thermique c peut tre prise dans le fuseau dfini par les limites infrieure et suprieure.
La limite infrieure est donne par :
pour 140C,
pour 140 < 160C,
pour > 160C,
27 / 122
limite infrieure
1,8
140C
limite suprieure
1,6
courbe retenue
Lamda (W/m K)
1,4
1,2
160C
1
0,8
0,6
0,4
0,2
0
0
200
400
600
800
1000
1200
Temprature (C)
28 / 122
29 / 122
Dans le cas dun modle avanc avec un feu naturel et pour la branche de temprature dcroissante, la rsistance de
compression rsiduelle du bton chauff une temprature max puis refroidi la temprature ambiante de 20C peut
tre donne daprs lannexe C de lEurocode 4 partie1-2 :
La relation contrainte-dformation pour lacier en contrainte uni-axiale est donne par la figure suivante :
Les rsistances des aciers de bton arm de types 3 et 4 sont donnes par : fe=fe . s
30 / 122
O :
s = 1 pour 20C<<400C
s = 1-0.85(T-400)/180 pour 400C<<580C
s = 0.15-0.15(T-580)/170 pour 580C<<750C
Les rsistances des aciers de prcontrainte sont donnes par : fp=fp . p
O :
p = 1 pour 20C<<400C
p = 1-0.85(T-400)/180 pour 400C<<580C
p = 0.15-0.15(T-580)/170 pour 580C<<750C
31 / 122
32 / 122
I.4. Conclusion :
Le dimensionnement au feu des structures en bton peut tre ralis avec diffrentes mthodes et diffrents degrs de
complexit. Certaines sont applicables directement (valeurs tabules, mthodes simplifies) pour diffrents types
dlments (poutre, poteaux, dalles, voiles,) pour une approche descriptive (feu normalis) et dautres sont plus
complexes et ncessitent un traitement numrique ou une programmation sur un logiciel spcialis (lments finis,
volumes finis, diffrences finies,) mais permettent de couvrir un domaine plus large (feu naturel, prise en compte
dune partie ou de la totalit de la structure). Ces logiciels de calcul gnralistes ou ddis sont multiples. Toutefois, ils
nont pas pu tre encore tous explorer dans le dtail afin de voir quelles sont leurs limites (lois de comportement
utilises, rfrences de validation, notices,). Une investigation de ces logiciels pour en donner une explication plus
dtaille et faciliter leur validation (par exemple par une fiche type doutil informatique) serait dun grand intrt et
permettrait une accessibilit par un plus grand nombre dutilisateurs dans le domaine de lingnierie de la scurit
incendie :
33 / 122
34 / 122
II 1. INTRODUCTION
Le bois est le matriau de construction pour lequel les travaux relatifs au comportement mcanique douvrages en
situation dincendie sont les moins nombreux. Les efforts de recherche se sont principalement concentrs sur le
comportement de poutres et ltude de la raction au feu du matriau compte tenu de sa combustibilit.
Depuis 2001, la scurit incendie a t clairement identifie comme lun des obstacles majeurs lutilisation du bois
dans la construction [1]. A linitiative des industries du bois en Europe, un rseau europen [2] a t mis en place pour
recenser les travaux et dvelopper les travaux de recherche.
Dans un premier temps, ce rapport dresse un tat de l'art des connaissances et des mthodes utilises ou utilisables
pour le calcul de la rsistance au feu des structures en bois. Les diffrentes mthodes sont ensuite utilises sur des
exemples couvrant une partie de leur champ d'application afin de comparer les diffrents niveaux de scurit apport et
de mieux apprhender les limites de certaines d'entre elles.
= 0,60 0,10
500
300
0,40
t = k x c1, 23
avec
t, le temps en minutes ;
xc, lpaisseur de bois carbonis ;
k, une fonction dpendant de la densit, du taux dhumidit et de lessence.
Pour des feux naturels (ou feux paramtriques), les tudes ralises flux constant se traduisent par une relation de la
forme [7] :
= 8,8
avec
Ces travaux sont complts par ltude de la masse volumique, de la conductivit thermique et de la chaleur spcifique
pour le domaine de temprature en situation dincendie [8, 9]. La figure II-1 rsume les variations de ces
caractristiques ncessaires pour un calcul thermomcanique.
35 / 122
Figure II-2 : Evolution de la loi de comportement longitudinal du bois en fonction de la temprature [10]
36 / 122
Figure II-4 : Variation du module longitudinal E0, des rsistances en flexion fm et en compression longitudinale
fc,0 en fonction de la temprature
Figure II-5 : Temps la rupture en fonction de linertie des poutres en bois lamell-coll GL24 [12]
II 1.3. Assemblages
Le comportement des assemblages entre lments structuraux et sa modlisation constitue un sujet essentiel pour
lanalyse de la scurit structurale en situation dincendie. En principe, les assemblages bois doivent prsenter une
rsistance au feu quivalente ou suprieure celle des lments, les jeux entre pices infrieurs 5 mm ninfluenant
pas la carbonisation des pices bois.
Un nombre trs limit de travaux aborde le comportement au feu des assemblages bois. La figure II-6 prsente des
rsultats obtenus pour des assemblages clous, boulonns, par connecteurs mtalliques et par anneaux. Ces deux
derniers modes dassemblage assurant le transfert deffort par contact surfacique, les rsultats montrent la ncessit
dune protection adapte la tenue au feu exige. Mobilisant un contact linique au sein des pices, les organes de
type tige prsentent un meilleur comportement.
37 / 122
Figure II-6 : Comportement des assemblages bois sous feu normalis et naturel [10]
38 / 122
Figure II-7 : Comportement des assemblages clous sous feu normalis [13]
Des tudes rcentes [14, 15, 16] ont abord le comportement au feu des assemblages par broches ou boulons. La
figure II-8 prsente les principaux rsultats issus des essais raliss. La variation interne des raideurs due au transfert
thermique constitue un paramtre important.
Figure II-8 : Temps la rupture selon le taux de chargement : (a) et (b)- assemblages bois-mtal [14] et boisbois [14], (c)- assemblages bois-mtal par broches 6,3mm [15], (d)- temprature au sein des organes [14]
Pour les organes de type tige, la portance localise du bois et le moment de plastification partielle des organes sont les
caractristiques matrielles de base. A partir de ces donnes, un modle de calcul de la rsistance par plan de
cisaillement aprs un temps t dexposition thermique est driv du modle de Johansen [17,18] :
avec, pour les assemblages bois-mtal :
39 / 122
f h , 2 ( , t )
f h ,1 ( , t )
Figure II- 9 : Comparaison qualitative des flammes observes lors des essais (t 6min.) avec les rsultats de la
simulation (FDS3 model)
Lapplication de lingnierie incendie aux structures bois constitue galement un second axe de recherche important au
niveau mondial pour le bois dans la construction. Cette approche permet en effet de lever les nombreux blocages
empiriques dus la combustibilit du matriau [21].
40 / 122
Figure II- 10 : Application de lingnierie incendie dans le cas dun btiment commercial en bois
[21]
En intgrant les diffrents modles thermiques, humains et dintervention, cette tude met en exergue la problmatique
associe la formulation de la fonction dtat limite. Pour les simulations de type Monte Carlo, il est retenu une relation
fonction des efforts normaux, des caractristiques de lincendie et de la section rduite ainsi quun facteur associ
lala de modlisation.
Dans ce cas, il est noter que seul le mode de rupture associ aux contraintes normales est intgr. Selon la
gomtrie de la poutre, le comportement de ce type dlment peut tre conditionn par les tats de contraintes de
cisaillement et de traction perpendiculaire. Ce point constitue un second mode de rupture devant tre intgr avec
ventuellement une hypothse dindpendance.
R =t
0
2
+ 15 0
Les valeurs des vitesses de combustion pour les panneaux drivs du bois sont celles du tableau 1. Les valeurs des
vitesses de combustion pour les bois naturels ou lamells ou abouts sont celles du tableau 2.
Ces valeurs sont affecter dun coefficient majorateur K1 qui est appliquer dans le cas des parois horizontales ou
inclines situes au-dessus du feu, l'inverse, un coefficient minorateur K1 est appliquer aux parois horizontales
situes sous le feu.
Un coefficient majorateur K2 est appliquer si le dbut de l'exposition au feu du bois ou des panneaux drivs du bois
est suprieur 15 minutes.
Les coefficients K1 et K2 peuvent se cumuler le cas chant.
41 / 122
42 / 122
(1)
et/ou
1,1 G + 0,7 p P + Pc
(2)
avec :
p = 0 ou 1
G = Sollicitation due aux charges permanentes
P = Sollicitation due aux surcharges dexploitation
Pc = Sollicitation due aux surcharges climatiques normales
Pc = soit Nn (neige normale),
soit Vn (vent normal), soit Nn + Vn.
Les rgles BF88 donnent galement la dure de protection thermique apporte par les matriaux base de pltre. Ces
valeurs figurent dans le tableau 3.
Tableau 3 : Dure de protection thermique apporte prendre en compte dans le calcul des
structures
Ltude de la fonction sparative pour les murs ou les planchers bois est traite dans les rgles BF88 sous forme
dexemples de solutions et de dispositions constructives, pour des parois verticales, horizontales et inclines selon
diverses configurations et dures de stabilit au feu.
Pour les panneaux drivs du bois, cette vitesse dpend de lpaisseur hp et de la masse volumique :
0 , ,t = 0
450 20
k hp
43 / 122
II 2.2.1.
Elments structuraux
Pour la dure exige, la section efficace est alors dfinie en considrant une profondeur de carbonisation def sur chaque
face expose :
d ef = d char ,n + k 0 d 0
Avec :
Figure II-12 : Vitesse de carbonisation des lments protgs : a) tf tchar ; b) tf > tchar
(tf : temps de chute de la protection))
Pour une protection dpaisseur hp, le dbut de la carbonisation du bois tchar est donn par :
44 / 122
panneaux bois
plaque de pltre de type A ou H
Dans le cas dune exposition un feu naturel (ou paramtrique), les calculs de la section rduite sont bass sur une
vitesse de carbonisation char et une profondeur de carbonisation dchar. Cette mthode nest prcise que pour des
lments bois rsineux non protgs.
La rduction de la section est alors dfinie par :
Av Ai hi
o : O est le facteur douverture : O =
A
At
v
0,5
et b est labsorptivit :
b = c
Av est la surface totale des ouvertures dans les parois verticales du compartiment,
At est la surface totale des parois du compartiment,
Ai est la surface douverture de hauteur hi dans une paroi,
t0 est la priode vitesse constante de carbonisation : t 0
= 0,009
qt ,d
O
40 min ,
b
4
h
d char
d char
Sur ces bases, la rsistance des composants est justifie sur la base de la formulation suivante :
E d , fi
k fi Rk
M , fi
Avec :
Ed , fi
: Valeur de calcul de leffet des actions en situation de feu dtermine conformment la norme EN 1991-1-2,
Rk
: Facteur de conversion
: est la valeur caractristique de la capacit rsistante dun assemblage pour une temprature normale sans
45 / 122
fictives (n) et la profondeur de carbonisation fictive (dchar,n). On calcule ensuite le facteur de modification vis--vis du
feu, kmod,fi, qui, pour t 20 minutes, prend les valeurs suivantes :
pour la rsistance en flexion :
k mod, fi = 1,0
1
p
200 Ar
k mod, fi = 1,0
1 p
125 Ar
k mod, fi = 1,0
1 p
330 Ar
II 2.2.2.
Lanalyse de la fonction sparative des systmes de mur et de plancher est traite dans lannexe E de lEurocode 5.
Les exigences relatives lisolation (critre I) et lintgrit (critre E) y sont supposes satisfaites condition de
respecter les dtails structuraux prconniss. Il est a not que les rsultats des calculs prsents au paragraphe II.6.4
ne correspondent pas toujours pas aux constatations faites au cours des essais.
Cette annexe sapplique aux lments dossature bois, aux parements constitus de panneaux base de bois et aux
plaques de pltre de type A, F et H.
Pour ce type dlments, il convient de vrifier que :
tins t req
o tins est le temps correspondant latteinte de laugmentation de temprature sur la face non expose qui conduit la
perte de lexigence disolation et treq est le temps exig pour la rsistance au feu.
La valeur de tins est calcule par le cumul des contributions des couches individuelles utilises dans llment
considr :
avec tins,0,i la valeur disolation de base de la couche i en minutes donne dans le tableau suivant en fonction de
lpaisseur en millimtres du panneau :
Nature du panneau
Contreplaqu
3
( 450 kg/m )
Panneaux de particules et
3
de fibre ( 600 kg/m )
Panneautage bois
3
( 400 kg/m )
Pltre de type A, F, R et H
tins (mn)
0,95 hp
1,1 hp
0,5 hp
1,4 hp
Lorsque les cavits sont remplies disolation dpaisseur hins, les valeurs de tins correspondantes sont :
pour de la laine de roche : tins , 0,i
pour de la laine de verre :
46 / 122
Le coefficient de position kpos dpend de la position du panneau considr par rapport au feu. Ses valeurs sont donnes
dans les tableaux 4 et 5 avec les quations suivantes :
0,02h p + 0,54
1
k pos = min
(E.9)
(E.10)
Tableau 4 : Coefficients de
position kpos pour des
panneaux simple couche sur la
face expose
Tableau 5 : Coefficients de
position kpos pour les
panneaux simple couche
sur la face non expose
Pour les murs avec panneaux double couche, les valeurs de kpos sont prendre dans le tableau 6 suivant :
47 / 122
Enfin, kj est un coefficient dassemblage qui tient compte de leffet des assemblages dans les panneaux de plaques. Le
tableau 7 prsente les valeurs de kj pour les panneaux de plaques de pltre non recouverts par une volige :
II 2.2.3.
Assemblages
Pour les assemblages, laction thermique induit une rduction de raideur des organes dassemblage :
k fi = k u f
, avec f = 0,2 pour les vis et les pointes et 0,67 pour les autres organes.
Pour des pices externes en bois dpaisseur t1, la rsistance minimale des assemblages non protgs est dfinie
forfaitairement par :
Lapproche simplifie est alors base soit sur le calcul de la protection ncessaire pour la tenue au feu requise treq (avec
treq 30 mn), soit sur la mthode des efforts rduits.
t ch t req 0,5t fi , d
Ou
t ch t req 1, 2 t fi , d
O
tch : temps de dmarrage de la combustion
treq : temps exig pour la rsistance au feu standard
tfi,d : rsistance au feu intrinsque de lassemblage non protg (paragraphe prcdent) charg par une valeur de calcul
de leffet des actions FED,fi.
Dans le cas de plaques mtalliques internes, la figure II-13 prsente les configurations ne ncessitant pas de protection
spcifique des plaques.
48 / 122
= e
k t
fi
o k est
II 2.3.1 Wood-frame
La catgorie wood-frame contient les structures ou les lments de structure raliss partir de pices en bois de
faible section. En raison de cette faible section, la rsistance au feu doit tre assure par des matriaux protecteurs
(plaques de pltre, etc.). Il nest pas donn de mthodes de calcul pour cette catgorie mais des mthodes tabules
sont proposes, regroupant un grand nombre de configurations et de solutions constructives testes sous conditions de
feu normalis.
Exemple de tableau indiquant des temps de rsistance au feu pour des structures du type wood-frame selon la
norme UCL-S101 :
49 / 122
Structure
supporte
Bois massif
(largeur x
hauteur) mm x
mm
Lamell coll
(largeur x
hauteur)
mm x mm
diamtre
mm
140 x 191
130 x 190
180
89 x 140
80 x 152
89 x 140
80 x 152
Arcs supports au
niveau ou proche du
niveau du sol
140 x 140
130 x 152
Poteaux
191 x 191
175 x 190
200
140 x 241 ou
191 x 191
130 x 228 ou
175 x 190
lment de structure
Toiture seule
Poteaux
Arcs supports
des murs ou
cules
par
des
Poutres et treillis
Planchers,
planchers
toiture
plus
Poutres, poutres
treillis et arcs
eff =
1,2 n
t 0,187
avec : n = 2 si les quatre cts sont exposs et n = 1 si la face suprieure de la poutre est protge,
et kfi prend les valeurs suivantes :
50 / 122
II 5. CODES DE CALCUL
Le calcul des stuctures bois en situation dincendie est possible laide de logiciels de calcul tels que Robobat ou
AccordBat qui permettent deffectuer des calculs froid sur les sections rduites.
Des codes de calcul gnralistes par lments finis permettent galement daborder le comportement
thermomcanique des structures bois : ABAQUS, ANSYS, ASTER, CASTEM, MARC.
Ces outils sont performants mais leur utilisation exige un haut niveau de qualification et induit un cot lev dtudes.
De ce fait, ils sont dapplication limite la recherche ou pour des projets trs particuliers.
Par ailleurs, une limitation importante de leur utilisation rsulte de labsence de critre de rupture en situation dincendie.
Comme pour les autres matriaux, seul un critre cinmatique est utilisable en terme de dformation limite ou de
vitesse de dformation.
51 / 122
52 / 122
Avec :
uinst =
M l2
9,6 E d I
Ed =
et
E0,05
u inst
l
400
7400
= 5692 Mpa
1,3
On obtient alors :
uinst =
M 6000 12
= 0,03124 M 15mm
9,6 5,69 75 1503
2
g
0,3 et M = ( g + q)l
g+q
8
On effectue maintenant le calcul de la poutre au feu. On considre alors la poutre expose au feu sur trois
faces (protge de laction directe du feu sur la face suprieure).
Selon lEurocode 5, pour le calcul de cet lment, les efforts et moments Efi,d peuvent tre obtenus partir des valeurs
normales de calcul Ed par la relation :
E fi , d = Ed
Avec :
M fi , d
Md
et
M fi =
q fi , d l
8
q fi = 1 g + fi q = g + 0,5q = 0,0695kN / m
Et = 0,447
fi =0,5 selon la norme EN 1990 correspondant un btiment dhabitation ou de bureau pour le bois.
II 6.1.2.1 Mthode de la section rduite
On calcule dans un premier les dimensions de la section efficace :
d eff = d char , n + k 0 d 0
Avec : d0 = 7 mm
dchar,n = n.t
d eff = 0,8t + k0 d 0
f d , fi = k mod, fi
k fi f m , k
m , fi
Avec : kmod,fi = 1
53 / 122
kmod, fi = 1,0
1 p
200 Ar
En effectuant les calculs avec les dimensions de section rsiduelle, on obtient alors : t = 37,5 minutes
Q = 1,1g + 0,8q
Ce qui donne : 1,1 x 0,032+ 0,8 x 0,075 = 0,0952 kN/m.
On calcule ensuite la section rsiduelle avec :
0 = 0,7 mm/min (bois naturels)
k1 = 1,6 (coefficient majorateur pour les faces horizontales au-dessus du feu)
Les contraintes lissue du temps de stabilit au feu sont gales :
2,25 fois la contrainte admissible si lpaisseur rsiduelle est > 30mm
1,75 fois la contrainte admissible si lpaisseur rsiduelle est < 30 mm
Avec ici adm = 11 MPa
On obtient alors : t = 45 minutes
On effectue les mmes calculs en faisant varier la section des poutres. On obtient les rsultats suivants :
section rduite
proprits rduites
120
Rgles BF88
100
t (min)
80
60
40
20
0
50 x
150
63 x
175
75 x
150
75 x
225
S (mm)
Figure II-16 : Dure de rsistance au feu pour trois mthodes diffrentes (C24)
54 / 122
dt (min)
S(mm)
Ecart DTU / section rduite
15
dt (min)
10
5
0
-5
50 x 63 x 75 x 75 x
150 175 150 225
100
x
100
100
x
150
120
x
175
100 100
x
x
200 300
120
x
210
125
x
125
175
x
200
200
x
200
-10
-15
S (mm)
Figure II-17 : Ecarts des dures de rsistance au feu de trois mthodes de calcul (C24).
En comparant les rsultats de la mthode de la section rduite et ceux obtenus avec la mthode des proprits
rduites, on peut constater que la mthode de la section rduite est toujours celle qui donne les valeurs les plus faibles
(les carts varient entre 1 et 3,5 minutes). Elle est par consquent lgrement plus scuritaire, ce qui expliquerait
pourquoi lEurocode la prconise plutt lautre mthode. Elle est par ailleurs plus facile appliquer et son domaine
dapplication est plus large que celui de la mthode des proprits rduites.
La comparaison entre les rsultats obtenus avec la mthode de la section rduite et ceux obtenus avec les rgles BF88
montre tout dabord que la mthode de la section rduite est plus scuritaire pour les petites sections. Cette diffrence,
de lordre de 10 minutes au dpart, diminue avec laugmentation de la section des poutres jusqu sinverser pour les
plus grosses sections. Les rgles BF88 donnent alors les valeurs les plus faibles. La hauteur rsiduelle de la poutre
dcrot plus vite dans le calcul BF88 que dans le calcul de la section rduite en raison du coefficient K1 de 1,6 pour les
faces horizontales exposes au-dessus du feu. Ainsi, pour des sections importantes, la valeur de la contrainte de calcul
est plus importante dans le calcul BF88.
Un autre calcul sur des poutres en lamell-coll, donc de sections plus importantes que les poutres en bois massif,
confirme ces rsultats. Le bois tudi est de classe GL28 :
55 / 122
dt (min)
2
0
-2
90 x
225
90 x
270
90 x
315
90 x
360
115 x
315
115 x
360
115 x
405
140 x
360
140 x
405
140 x
450
140 x
495
-4
-6
-8
S (mm)
dt (min)
4
3
2
1
0
90 x
225
90 x
270
90 x
315
90 x
360
115 x
315
115 x
360
115 x
405
140 x
360
140 x
405
140 x
450
S(mm)
Figure II-18 : Ecarts des dures de rsistance au feu de trois mthodes de calcul (GL28).
rel , y =
y =
l
I
A
f c ,0, k
E0,05
avec :
- c,0,d
et qui vaut :
c , 0 ,d
k c , y f c , 0 ,d
avec :
c , 0 ,d =
N
A
f c ,0,d = k mod
f c , 0 ,k
56 / 122
kc, y =
1
k y + k y rel , y
2
o ky vaut :
c, 0,d
k c , y f c , 0 ,d
N
N
= 100 100 =
1
0,285 12,93 36850,5
En prenant
g
0,3 , il vient : g = 7,63 kN et q = 17,80 kN.
g+q
q fi = 1 g + fi q = g + 0,5q = 16,53kN
fi =0,5 selon la norme EN 1990 correspondant un btiment dhabitation ou de bureau pour le bois.
II 6.1.3.1 Mthode de la section rduite
On effectue le mme calcul de la profondeur de carbonisation que pour les poutres. En revanche, la vrification de
llment seffectue au flambement comme au paragraphe prcdent mais en ne considrant plus que la section
efficace. Comme nous avons ici une section de poteau carre, les valeurs des lancements selon les axes y et z sont
similaires.
En effectuant les calculs avec les dimensions de section efficace, on obtient alors
c , 0 ,d
k c f c , 0 ,d
linstant
t = 14 minutes.
k mod, fi = 1,0
1 p
125 Ar
En effectuant les calculs avec les dimensions de section rsiduelle, on obtient alors
c , 0 ,d
k c f c , 0 ,d
linstant
t = 15 minutes.
Q = 1,1g + 0,8q
Ce qui donne ici : 1,1 x 14,64+ 0,8 x 34,16 = 43,43 kN
On calcule llancement comme prcdemment.
Le CB71 prcise, pour la vrification des pices soumises la compression simple, de multiplier la contrainte due
leffort normal par le coefficient
K=
1
k
coefficient k est reprsent par les ordonnes de deux courbes qui ont pour expression :
57 / 122
pour 75
k = 1 0,8
100
2
k=
et
En multipliant la contrainte admissible par 2 ou 1,5 suivant la valeur de lpaisseur rsiduelle (> ou < 30 mm), on
obtient alors le temps de rupture suivant : t = 14,5 minutes.
On effectue les mmes calculs pour des poteaux de diffrentes sections. On obtient les rsultats suivants :
section rduite
proprits rduites
rgles BF88
80
70
60
t (min)
50
40
30
20
10
0
100 x 100 100 x 150 100 x 200 150 x 200 200 x 200 200 x 250 200 x 300 300 x 300
S (mm)
2
1
0
100 x 100
-1
100 x 150
100 x 200
150 x 200
200 x 200
200 x 250
200 x 300
300 x 300
-2
-3
-4
S (mm)
dt (min)
3,5
3
2,5
2
1,5
1
0,5
0
100 x
100
100 x
150
100 x
200
150 x
200 x
200
200
S (mm)
200 x
250
200 x
300
300 x
300
Figure II-20a : Ecarts des dures de rsistance au feu de trois mthodes de calcul.
Dans un premier temps, nous pouvons constater que, comme pour les poutres, la mthode des proprits rduites est
celle qui donne les rsultats les moins scuritaires. Nous observons des carts peu variables entre cette mthode est
celle de la section rduite, compris entre 1 et 4 minutes.
En comparant la mthode de la section rduite et celle du DTU bois feu, on constate une que les valeurs sont
relativement proche avec une prdominance vers des rsultats plus scuritaire dans le cas de la mthode de la section
rduite. Cependant, cette tendance semble sinverser lorsque les sections des poteaux deviennent plus importantes
58 / 122
II 6.2 Exemple dutilisation dun logiciel pour le calcul dune poutre bois
Le module de dimensionnement Bois de ROBOT Millenium permet de vrifier, dimensionner et doptimiser les profils
bois pour tous les lments de la structure conformment aux prescriptions du CB71 et de lEurocode 5. Ce logiciel
permet deffectuer des calculs froid sur la section rduite dun lment.
Nous allons ici effectuer la vrification au feu de la poutre de 75 x 150 mm tudie prcdemment en exemple.
59 / 122
II 6.3 Utilisation des mthodes avances pour le calcul dune poutre bois
Le principe est ici dutiliser les donnes de lannexe B de lEurocode 5 partie 1-2 ( Mthodes de calcul avances )
pour effectuer un calcul numrique par lments finis dune poutre bois charge uniformment. .
g = 0 + 345 log(8t + 1)
o 0 est la temprature initiale en K et t est le temps dexposition en minutes.
A un instant t, la densit de flux reue par la pice de bois vaut alors :
o Tgt est la temprature du gaz en K, T est la temprature du solide en K, h est le coefficient dchange convectif
(W/m.K), est le coefficient rsultant dmissivit et est la constante de Boltzmann (5,67.10-8 W/mK4).
60 / 122
Dans le cadre de cette tude, nous prenons un coefficient dchange convectif gal 25 W/m.K et un coefficient
rsultant dmissivit gal 0,8. Ces valeurs sont justifies exprimentalement.
Les proprits thermophysiques des bois rsineux sont donnes par lannexe B :
61 / 122
Figure II-25 : champs de temprature en Kelvin du quart de poutre aprs heure de simulation
Proprits du bois
EL (Mpa)
ET (Mpa)
G (Mpa)
11000
370
690
II 6.3.2.4 Rsultats
Compte tenu de labsence de critre de rupture en situation dincendie pour le bois, il nest pas possible de dterminer,
comme avec les autres mthodes, une dure de rsistance au feu de la poutre. Nous pouvons en revanche connatre
ltat de contraintes en tout point de la poutre jusquau temps de simulation calcul.
La figure ci-dessous montre les valeurs des contraintes dans la poutre. La figure II-26 suivante est un dtail des
contraintes au milieu de la poutre dforme avant et aprs chargement thermique.
62 / 122
Figure II-27 : Valeurs maximales des contraintes au milieu de la poutre dforme t=0 et t=1800 s
63 / 122
On peut voir apparatre les zones de fortes contraintes sur la figure de droite et constater les faibles valeurs de
contraintes sur la face infrieure, totalement carbonise. Les deux diagrammes ci-dessous reprsentent ltat de
contrainte sur la face intrieure de la poutre dans les deux configurations prcdentes.
64 / 122
180
Laine de verre
laine de roche
Cavit vide
160
140
t (min)
120
100
80
60
40
20
0
0
50
100
150
200
250
p isolant
140
Laine de verre
laine de roche
Cavit vide
120
100
t (min)
80
40
20
0
0
50
100
150
200
250
p isolant
Figure II-29 : Comparaison des dures de rsistance au feu (isolation) de diffrentes configurations de parois
On peut constater que la courbe de la laine de roche nest pas linaire. La variation des coefficients de position au-del
de 95 mm augmente de manire importante la valeur disolation. Les diffrences entre les cas o la cavit est vide et
ceux o elle est remplie de laine de roche sont assez importantes. Pour les cavits les plus larges, il est possible davoir
3 voire 4 fois la valeur disolation du mur avec cavit vide. Des essais canadiens ont montr quavec de la laine de
roche dans la cavit et une protection par plaque de pltre type F, on pouvait augmenter la rsistance au feu de 50%
mais quil tait difficile dobtenir un gain disolation suprieur. Dautres tests ont galement montr que dans la plupart
des cas, la prsence de laine de verre dans la cavit naugmentait pas la rsistance au feu du mur.
Comparons maintenant les rsultats obtenus avec une protection simple couche et ceux obtenus avec une protection
double couche. Lisolation est faite en laine de roche, les parements sont en plaques de pltre de 12 mm et nous
considrons des joints parfaits entre les panneaux de protection :
160
140
simple couche
120
Double couche
t (min)
100
80
60
40
20
0
0
50
100
150
200
250
p isolant
Figure II-30 : Comparaison des dures de rsistance au feu (isolation) de diffrentes configurations de parois
A partir denviron 120 mm disolant, on constate que les murs avec parement simple couche donnent de meilleurs
rsultats que les murs avec parement double couche. Ainsi, partir dune certaine paisseur de laine de roche, un mur
avec parement simple couche auquel on aurait rajout une couche de 12 mm de pltre perdrait de sa rsistance au feu.
Ce rsultat apparat comme erron et sa validit est compromise. Il est principalement d aux coefficients de position
des parements prcds de laine de roche qui sont constants pour les murs doubles couches et qui augmentent avec
lpaisseur de laine de roche pour les murs simple couche.
Ainsi, il apparat daprs cette tude que lannexe E de lEurocode peut donner des rsultats trs contradictoires. Le
comportement des cloisons en fonction de leur configuration peut diffrer normment.
65 / 122
25
EC5
DTU Certifi
DTU non certifi
20
35
30
25
t (min)
t (min)
15
10
20
15
10
EC5
DTU , au contact du bois
5
0
0
7
12
17
e (mm)
22
27
12
17
22
e (mm)
Figure II-31 : Comparaison des dures de rsistance au feu (isolation) de deux mthodes diffrentes
De ces deux courbes, on peut constater que les rsultats de lEurocode et du DTU sont assez proches (cart maximal
de 3 minutes). Dans la plupart des cas, et surtout pour les faibles paisseurs, lEurocode est moins scuritaire que le
DTU.
II 6.4.4 Comparaison de solutions du DTU avec les rsultats Eurocode pour les parois
verticales
On tudie les quatre cas ci-dessous dlments de mur intrieur simple. Le DTU annonce une rsistance au feu de
heure pour ces solutions.
e2 = 100 mm
e3 = 100 mm
Liaisons remplies
kj le plus dfavorable
Liaisons remplies
kj le plus dfavorable
Cas 1
57,36
38,08
65,36
46,08
Cas 2
Non applicable
Non applicable
90
61,74
Cas 3
59,43
40,46
67,43
48,46
Cas 4
Non applicable
Non applicable
95,5
60,3
Les cas avec laine de roche dans les cavits prsentent des valeurs toujours suprieures heure, les cas avec joints
parfaits entre panneaux atteignent mme 1h30. Le DTU ne prsente pas plusieurs types de liaisons entre panneaux, il
stipule uniquement que les protections thermiques doivent tre au contact de la surface du bois et ne prsenter aucun
interstice ou fente suprieur 1 mm, localement sur des espaces limits 2 mm.
Pour les cas avec laine de verre, les rsultats sont moins levs mais toujours suprieurs heure. LEurocode valide
donc tous les cas des DTU sous certaines conditions. Cependant, la diffrence entre les deux codes peut se rvler
trs importante.
II 6.4.5 Comparaison des solutions du DTU avec les rsultats Eurocode pour les parois
horizontales
On ralise la mme dmarche avec les solutions pour parois horizontales du DTU. On considre les 3 cas ci-dessous
de plancher avec cavit vide. La rsistance au feu indique dans les rgles BF88 pour cette solution est de heure.
66 / 122
Le parement suprieur est ralis avec du contreplaqu de 19 mm dpaisseur (valeur la plus faible) et la cavit a une
largeur de 45 mm (valeur minimale de lEurocode) :
kj le plus
dfavorable
Cas 1
29,79
16,08
Cas 2
27,04
16,34
Cas 3
48,25
21,91
On constate de nouveau que pour des joints parfaits, les rsultats sont nettement suprieurs ceux du DTU. Pour les
configurations de joints les plus dfavorables, on est juste au-dessus des rsultats DTU.
Dautres tudes ont t ralises avec des planchers dont la cavit est remplie disolant et les conclusions sont les
mmes que pour les solutions verticales : lEurocode semble moins scuritaire que le DTU BF88 car les valeurs quil
fournit sont la plupart du temps plus leves, sauf dans des cas particuliers de liaisons. Dans le cas o lisolant est de la
laine de roche, les valeurs fournies par lEurocode peuvent dpasser le double de celles du DTU.
II 7. Conclusion
Les mthodes de calcule de la rsistance au feu des structures en bois utilises ce jour sont essentiellement les
mthodes simplifies. Celles-ci sont dcrites dans les diffrents textes normatifs et sont utilisables par l'ensemble des
intervenants du secteur de la construction. Ces mthodes, compares des rsultats exprimentaux, offrent un niveau
de fiabilit satisfaisant pour un calcul par lments. De plus, elles peuvent tre facilement appliques grce des
logiciels de calcul de structures bois du commerce ce qui en simplifie l'accs.
La comparaison des diffrentes mthodes simplifies a permit plusieurs constatations :
Premirement, la mthode des proprits rduites de lEurocode est de toutes ces mthodes celle fournissant les
rsultats les moins scuritaires. Son champ dapplication est, de plus, trs limit (bois rsineux uniquement, tenue au
feu obligatoirement suprieure 20 minutes). Il est noter que la mthode des proprits rduites na pas t retenue
dans lannexe nationale de lEurocode 5. Il apparat ainsi prfrable dutiliser la mthode de la section rduite. Celle-ci
et les rgles BF88 donnent des rsultats sensiblement diffrents en fonction de la section des lments tudis.
Deuximement, une mthode de justification par le calcul des lments de mur et de plancher proposes dans lannexe
E de lEurocode, qui contient de nombreuses restrictions au domaine dapplication, donne parfois des valeurs qui
semblent surleves et apparat moins scuritaire que les solutions proposes par le DTU BF88.
Les mthodes avances sont en revanche mois bien dcrites dans les textes normatifs et sont donc rarement
pratiques sur des projets rels. L'utilisation correcte de telles mthodes repose encore largement sur les
connaissances de l'ingnieur charg de les mettre en uvre. Dans ce cadre, le calcul par lment est ralisable mais
prsente un intrt limit par rapport aux mthodes simplifies. Les mthodes avances prsenteront un plus grand
intrt lorsque les travaux en cours concernant le comportement des matriaux et des assemblages permettront le
calcul d'ouvrages ou de partie d'ouvrage.
67 / 122
Rfrences
[1]
[2]
[3]
R.H. White
Charring rate of composite timber products
Proceedings. Wood and fire safety; Zvolen, Slovakia, 2000.
[4]
A. Frangi, M. Fontana
Charring rates and temperature profiles of wood section
Fire and Materials, Vol 27, p. 91-102, 2003.
[5]
[6]
[7]
E. Mikkola
Charring of wood based materials
Proceedings third international symposium Fire safety science, London, 1991.
[8]
P. Clancy
A parametric study on the time-to-failure of wood framed walls in fire.
Fire Technology, p.243-269, 2002.
[9]
[10]
A. Buchanan
Fire performance of timber construction
Progress in Structural Engineering and Materials, Vol 2, p.278-289, 2000.
[11]
A. Frangi, M. Fontana
Effect of temperature on the structural behaviour of bond lines in glulam
Int. Symposium on advanced timber elements, COST E29, Florence, 2004.
[12]
E.L. Schaffer
Structural fire design: wood
Research paper FPL 450, For. Prod. Lab., Wisconsin, USA, 1984.
68 / 122
[13]
J. Noren
Load-bearing capacity of nailed joints exposed to fire.
Fire and materials, Vol.20, p. 133-143, 1996.
[14]
[15]
A. Frangi, A. Mischler
Fire tests on timber connections with dowel-type fasteners
CIB-W18, XXXVIIme Meeting, Edinburgh, UK, 2004.
[16]
[17]
D. Povel
Tragfhigkeit von Holzverbindungen mit stabfrmigen Verbindungsmitteln im Brandfall.
PhD Thesis, Technical University Berlin, 2002.
[18]
[19]
[20]
T. Korhonen, J. Hietaniemi
Fire safety of wooden faades in residential suburb multi-storey buildings
VTT working papers N32, Espoo, Finland, 109 p., 2005.
[21]
[22]
[23]
[24]
[25]
K.H. Almand, L.T. Phan, T.P. McAllister, M.A. Starnes, J.L. Gross
R&D roadmap for structural fire safety design and retrofit of structures.
Proceedings NISTIR 7133, NIST-SFPE workshop, Baltimore, USA, April 2004.
69 / 122
[26]
70 / 122
71 / 122
III.1. INTRODUCTION
Le prsent chapitre fait un tat de lart des outils (plus ou moins sophistiqus) actuellement disponibles en
ingnierie du comportement au feu pour les ouvrages structure en acier ou mixtes acier-bton.
Conductivit thermique :
La conductivit thermique de lacier est beaucoup plus importante que celle du bton. Elle peut tre dfinie partir des
relations suivantes (voir figure III-1) :
pour 20 C a < 800 C :
a = 27,3 W/mK
Chaleur spcifique :
La chaleur spcifique de lacier de construction (en J/KgK) a tendance augmenter avec la temprature. La courbe
donnant sa variation en fonction de la temprature peut tre dfinie par les relations suivantes (voir figure III-1) :
Pour :
20 C a < 600 C
900 C a 1200 C
ca = 650
Nous pouvons noter que la courbe donne par lEurocode 4 prsente un pic autour de 735C, dont le sommet denviron
2500 J/KgK napparat pas sur la figure. Ce pic tient compte du caractre endothermique du changement de phase
cristallographique qui a lieu vers cette temprature.
300
0
250
0
200
0
150
0
100
0
50
0
0
Chaleur spcifique
(J/kgK)
acier
bton
20
0
40 60 80 100 120
0
0
0
0
0
Temprature
(C)
Conductivit
6 thermique (W/mK)
0
5
0
4
acie
0
3
r
0
2
0
1
bton
0
0
20 40 60 80 100 120
0
0
0 0
0
0 0
Temprature
(C)
Masse volumique :
Lincendie a peu dinfluence sur la masse volumique de lacier. Une valeur de 7850 kg / m3 est gnralement adopte
pour tous les types dacier.
72 / 122
O l est la longueur 20 C ; l est la dilatation induite par la temprature ; a est la temprature de lacier [C].
En pratique, on nobserve pas de grandes diffrences entre la dilatation thermique des diffrents types dacier au
carbone. La courbe de dilatation montre un ralentissement, avec dans certains cas une phase passagre de contraction
pour les tempratures dans la plage 750 850C. Ce ralentissement est d au changement de phase cristallographique
(transformation allotropique). Au-del de 800 900C, la dilatation reprend plus rapidement quaux tempratures
infrieures et suit une loi linaire en fonction de la temprature. Toutefois, dans leur domaine de validit commun, les
deux normes donnent des valeurs trs proches.
L/L (x103)
20
15
10
acier
5
0
20
0
40
60
80
100
0temprature
0
0 (C)0
120
0
Contrainte
Module tangent
p,
Ea,
Ea,
p, < < y,
y, t,
t, < < u,
= u,
Paramtres
f p, - c + (b/a) a 2 - ( y, - )2
b( y, - )
fy,
f y, [1 - ( - t, )/ ( u, - t, )]
0,00
p, = fp, / Ea,
a a 2 ( y, - )2
0,5
Contrainte
f y,
0,5
f p,
y, = 0,02
t, = 0,15
u, = 0,20
2
a = ( y, - p, )( y, - p, + c/ E a, )
2
2
b = c ( y, - p, ) E a, + c
Fonctions
c=
(f
y,
- f p,
E a, = tan
( y, - p, )E a, - 2( f y, - f p, )
p,
y, =0.02
t, =0.15
u, =0.20
fy, = ky, fy
Ea, = kE, Ea
Les facteurs de rduction appliquer en fonction de la temprature ces paramtres sont prsents sur la figure III-4.
73 / 122
Facteur de rduction
(par rapport fy)
pour la limite
dlasticit efficace
Facteur de rduction
(par rapport fy)
pour la limite de
proportionnalit
Facteur de rduction
(par rapport Ea)
pour la pente du domaine
lastique linaire
ky, = fy, / fy
kp, = fp, / fy
kE, = Ea, / Ea
20 C
1,000
1,000
1,000
100 C
1,000
1,000
1,000
200 C
1,000
0,807
0,900
300 C
1,000
0,613
0,800
400 C
1,000
0,420
0,700
500 C
0,780
0,360
0,600
600 C
0,470
0,180
0,310
700 C
0,230
0,075
0,130
800 C
0,110
0,050
0,090
900 C
0,060
0,0375
0,0675
1000 C
0,040
0,0250
0,0450
1100 C
0,020
0,0125
0,0225
1200 C
0,000
0,0000
0,0000
Rsistance
% de la valeur normale
Rsistance
effective
10
Contrainte normalise
20 200
400
0.
500
6
4
0.
Module
2
lastique
0 30 60
0.
600
C
700
800
0.
90 120
Temperature (C)
Dformation ( %)
Masse volumique :
Lincendie a peu influence sur la masse volumique de lacier. Une valeur de 7850 kg / m3 est gnralement adopte
pour tous les types dacier.
74 / 122
Contrainte
Module tangent Et
c,
E
1+ a b
E (1 + a b a b b )
(1 + a b ) 2
c, < u,
f 0.2 p, - e + (d/c) c ( u , )
Paramtres
a=
b=
Fonctions
d + ( u , - )
f0.2p ()
= arctan( Ect())
(1 c , Ect , / f 0.2 p , ) E a , c ,
( E a , c, / f 0.2 p , 1) f 0.2 p ,
e
c 2 = ( u , c , ) u , c , +
Ect ,
d 2 = e ( u , c, ) Ect , + e2
e=
Contrainte
c c 2 ( u , - )
Ea , c , f 0.2 p ,
f 0.2 p , c,
fu()
( f u , f 0.2 p , ) 2
= arctan(E( ))
u()
c()
( u , c , ) Ect , 2 ( f u , f 0.2 p , )
75 / 122
Figure III-7 : Facteur de massivet Am / V pour les lments en acier non protgs
76 / 122
Am /V
h net,d t
ca a
o :
ksh
Am/V
caa
h net,d est la valeur de calcul du flux thermique net [W/m/K]. Ce terme qui reprsente laction thermique, dpend du
modle de feu utilis (conditions de feu normalis, de feu naturel) ;
Le flux thermique est exprim partir de la valeur de la temprature des gaz chauds, et est habituellement dcompos
en la somme de deux flux distincts, savoir : une composante convective et une composante radiative. En situation
dincendie, les deux modes de transmission de la chaleur (rayonnement et convection) sont cumulatifs. Le flux de
chaleur net transmis la surface de llment lors de son exposition au feu est dfini par la relation :
77 / 122
est le facteur de forme (<1). Dans de nombreux cas pratique, le facteur de forme est pris gal 1 ;
m est lmissivit de paroi. Elle dpend du type de matriau appliqu la surface. Pour lacier, m = 1 ;
f est la temprature du rayonnement proximit de llment (C) ;
m est la temprature de surface de llment (C) ;
hnet,c = c ( g m )
o
c est le coefficient de convection (W/mK). Il varie en pratique entre 25 (conditions de feu normalis ISO) 50
W/mK (condition de feu hydrocarbure) ;
g la temprature des gaz proximit de llment (C) ;
m la temprature de surface de llment (C) ;
Le flux de chaleur net sur la surface expose des lments de structure dpend donc de deux paramtres : lmissivit
de paroi et le coefficient dchange par convection hc.
Pour les sections de type I sous conditions de feu nominal, leffet dombre est dcrit de faon raisonnable en prenant :
ksh = [Am/V]box/[Am/V]
o [Am/V]box est la valeur en caisson du facteur de massivet.
Dans tous les autres cas, la valeur de ksh doit tre prise comme :
ksh = [Am/V]box/[Am/V]
a.t =
p / d p Ap
( g .t a.t )t (e / 10 1) g ,t
c a a V 1 + / 3
avec
cpp
ca a
dp
Ap
V
o
Ap/V est le facteur de massivet dun lment isol par un matriau de protection contre le feu ;
Ap est la surface concerne de la protection par unit de longueur (m) ;
V est le volume de llment par unit de longueur (m3) ;
Ca est la chaleur spcifique de lacier (J/kg/K) ;
Cp est la chaleur spcifique de la protection (J/kg/K) ;
dp est lpaisseur de protection (m) ;
.t est lintervalle de temps (s) ;
a,t est la temprature de lacier au temps t (C) ;
g,t est la temprature des gaz au temps t (C) ;
g,t est laugmentation de temprature pendant lintervalle de temps .t au temps t (C)
p est la conductivit thermique du matriau de protection (W/m/K) ;
a est la masse volumique de lacier (kg/m3) ;
p est la masse volumique du matriau de protection (kg/m3) ;
Il convient de souligner que les caractristiques des matriaux de protection sont habituellement dtermines sous
condition de feu normalis. Elles sont donc en principe limites ces conditions. En effet, aucune information
systmatique nest disponible pour dcrire lvolution de ces caractristiques lorsque les conditions de feu sont
considrablement diffrentes des conditions de feu normalis. Aussi, lapplication du concept de scurit incendie bas
sur le feu naturel pour vrifier la rsistance au feu des structures en acier protges doit tre mene avec prcaution.
En particulier, un tel calcul ne peut tre effectu que si lon dispose des donnes appropries ou que lon peut montrer
que les conditions de feu nont pas deffets significatifs sur les proprits thermiques et sur lintgrit des matriaux de
protection. En particulier, iI est important de sassurer que les dformations dun llment de structure expos au feu ne
conduiront pas une dtrioration prmature de la protection rapporte sur cet lment
78 / 122
en situation dincendie, il est suppos que les changes de chaleur se font essentiellement par convection et
rayonnement au niveau des surfaces exposes. A lintrieur des lments homognes tels que lacier, la chaleur est
uniquement transfre par conduction. Pour les matriaux poreux tels que le bton ou lorsquil existe des cavits
lintrieur des lments (vides), lchange de chaleur est en principe beaucoup plus complexe : les trois phnomnes
rayonnement, convection et conduction tant mis en cause simultanment, auquel sajoute le transferts de masse.
Toutefois, dans la pratique on peut se limiter au processus dchange prpondrant en tenant compte des processus
secondaires par lintermdiaire de facteurs de correction. Ainsi, on suppose que le bton est un matriau homogne et
que lchange thermique se fait principalement par conduction. Les changes convectifs et rayonnant qui se
manifestent dans les pores sont considrs comme des phnomnes de second ordre, intervenant dans le calcul par
une majoration approprie du coefficient de conductivit thermique. Par ailleurs, les transferts de masse sont
gnralement ngligs et on se contente de reprsenter la vaporisation de leau libre contenue dans le bton. Celle ci
est suppose uniformment rpartie dans le bton et pour obtenir numriquement les paliers de vaporisation, la chaleur
spcifique du bton peut tre artificiellement augmente laide dun pic, situ entre 100 et 200C. Enfin, le contact
entre lacier et le bton peut tre suppos parfait.
Poutres mixtes
Poteaux mixtes
Dalles
(isolation thermique)
Dalle
Bton
pour
isolation
Figure III-9 : Domaine dapplication des mthodes de calcul par valeurs tabules
Lapplication de la mthode de calcul par valeur tabule peut tre ralise dans deux situations diffrentes, lune en cas
de vrification o les dimensions des lments de structure sont dj connues et lautre en cas de prdimensionnement o, en revanche, les actions mcaniques de calcul sont dj dfinies.
79 / 122
En cas de vrification, les dimensions de la section de llment ainsi que sa capacit portante
Rd
vis vis de la
charge sont dj connues. Par contre, on doit calculer leffet des actions mcaniques en situation dincendie
E fi,d
afin
fi,t = E fi,d R d
fi, d
d et
on connat les effets des actions
partir des combinaisons de charges pour le dimensionnement
temprature normale et en situation dincendie. Dans ce cas, on peut admettre de manire scuritaire que le niveau de
=E
fi,t
fi,d
d
chargement sexprime par
. A partir de cette valeur du niveau de chargement et du degr de
rsistance au feu normalis, les dimensions minimales de la section ainsi que les dispositions constructives
correspondantes de llment peuvent tre dfinies. Enfin la section ainsi dfinie doit tre vrifie pour la condition
Rd Ed
Tirants
Poutres
(acier et mixtes)
Poteaux
Dalles mixtes
(isolation thermique et
rsistance mcanique)
80 / 122
Elments soumis aux effets combins en moment de flexion et effort axial, comme les poteaux lancs chargs avec
une importante excentricit, les poutres avec dversement latral, etc. ; pour ce type dlment, le modle de calcul
simplifi prend en compte leffet combin des diffrentes sollicitations en assemblant les deux modles de calcul cidessus correspondant la condition simple de chargement en flexion ou en compression.
Les mthodes de calcul simplifies fournissent beaucoup plus de libert dans le choix des sections transversales. Bien
que plus simples, elles comportent un nombre important de paramtres rendant le calcul manuel relativement long. En
outre, elles ne permettent que le dimensionnement en situation dincendie normalis (ISO).
, a
a
tape 1
tape 2
Temprature
acier a
tape 3
tape 4
Rreq
temps
Efi,d,t Rfi,d,t
Figure III-11 : Principe dapplication des mthodes de calcul simplifies pour les lments en acier
Note 1 : lincendie, des mthodes de calcul sont galement proposes pour lacier inoxydable. Ces mthodes sont
bases sur celles relatives aux aciers au carbone. Elles sappliquent :
Aux lments dune classe de section 1, 2, 3 ou 4 (lments prsentant des risques de voilement local importants);
Aux tirants (lments tendus) ;
Aux poutres sur appuis simples (lments flchis) ;
Aux poteaux (lments comprims et lments comprims flchis).
Les sections transversales en acier inoxydable sont classes de la mme faon que les aciers au carbone. Cest le
rapport largeur-paisseur des parois dune section transversale qui dtermine si cette paroi a une grande capacit de
rotation ou non, ou bien si elle est sensible ou non au voilement.
Note 2 : LEurocode 3 Partie 1.2 donnent galement des mthodes permettant de vrifier la rsistance au feu des
boulons soumis au cisaillement, des boulons tendus (prcontraint on non prcontraint) et des soudures.
81 / 122
Poteaux
crit
Temprature critique : cr
Facteur massivet : Am / V
Courbe ISO
a
Temprature acier a
tfi
Temps
82 / 122
- les effets des non-linarits gomtriques (grands dplacements, grandes dformations, ) amplifies par laction
de la temprature (effets des contraintes thermiques dues aux dilatations thermiques diffrentielles des
matriaux) ;
- les effets des proprits non linaires des matriaux, comprenant les effets de dchargement sur la rigidit de la
structure et les variations des proprits mcaniques avec la temprature. A tempratures leves, les relations
contraintes dformations de lacier sont bases sur un modle elliptique linaire en opposition avec la
schmatisation dans des conditions de temprature ambiante, pour laquelle un modle bilinaire suffit ;
- La dilatation thermique des matriaux ;
- Les effets des contraintes et des dformations induites la fois par l'accroissement de temprature et par des
diffrences de temprature ;
- A condition que les relations contrainte-dformation proposes dans les Eurocodes soient utilises, il n'est pas
ncessaire de prendre explicitement en considration les effets du fluage aux tempratures leves.
Pour les structures en acier et mixte acier-bton, lapplication des modles de calcul avancs dans lanalyse de
structure globale ncessite de prter une attention particulire aux points suivants :
Concernant les modles de matriaux, on doit considrer :
- la dcomposition de dformation totale en plusieurs parties aux tempratures leves ;
- le modle de matriau cinmatique pour lvolution de la temprature ;
- la rsistance de certains matriaux, par exemple le bton pendant la phase de refroidissement ;
Dcomposition de la dformation
Acier
(matriaux isotrope)
t = th + ( + c)+ r
t: dformation totale
t: dformation due dilatation thermique
: dformation due aux contraintes
r: dformation due ventuelle contrainte rsiduelle
c: dformation due au fluage
y
Section
th
t
parallle
(
1, = 0)
d
Compression
1 = (t)
1 = (t)
2 = (t+
t)
c
r
Distribution temprature
z = constante
Bton
(matriau anisotrope en
compression-traction)
d
Parallle
2, = 0)
(
d
Traction
Dformation unitaire
2 = (t+
t)
Le rgime transitoire de lchauffement de la structure lors de lincendie exige dutiliser une procdure de rsolution
incrmentale (pas pas) et itrative, prenant en compte lvolution et la distribution des tempratures chaque pas
de temps, ainsi que leurs influences sur les proprits mcaniques des matriaux : Une analyse pas pas afin
dobtenir ltat dquilibre de la structure pour diffrents instants de feu, c'est--dire pour diffrents champs de
temprature de la structure et pour chaque pas de calcul, une procdure de rsolution itrative pour rtablir lquilibre
de la structure qui se comporte de manire lasto-plastique et non-linaire ;
Les conditions aux limites existantes de la structure ;
La condition de chargement de la structure modlise en situation dincendie ;
Les modles de matriaux utiliss dans la modlisation numrique doivent tre reprsentatifs des lois de
comportement relles des matriaux aux tempratures leves ;
Dans lapplication des modles de calcul avancs aux structures en acier et mixtes acier-bton, il faut tre vigilant sur
les aspects spcifiques qui en gnral ne sont pas pris en compte dans la modlisation directe ; par exemple, la
rupture des barres darmature due une longation excessive, lcrasement et la fissuration du bton, la rsistance
des assemblages, la connexion entre lacier et le bton, etc.
Comme cela est gnralement admis en situation dincendie, dans un souci de simplification, il ny a pas dans ces
modles de calculs de rel couplage entre le comportement thermique et le comportement mcanique. Pour les besoins
de lanalyse mcanique, la distribution des tempratures au sein des lments est obtenue sparment, soit laide
dun modle thermique, soit laide des rsultats dessais.
En gnral, la vrification de la rsistance au feu des structures se limite celle des lments isols soumis lincendie
conventionnel. Cette procdure danalyse possde un avantage incontestable par sa simplicit dapplication. Toutefois,
les informations que donne ce type danalyse restent assez restreintes en ce qui concerne le comportement rel dune
structure globale soumise lincendie. Dans cette optique, lingnierie incendie soriente de plus en plus vers une
analyse du comportement global des structures. Cette dmarche reprsente en particulier les avantages suivants :
une meilleure comprhension du comportement global des structures avec la prise en compte des effets des
interactions entre les diffrents lments des structures, des transferts de charges, etc. ;
lintroduction dun chauffement localis dune structure globale au lieu dun chauffement gnralis lorsque le
compartimentage contre lincendie est prise en considration ;
ltude de la stabilit au feu dune structure globale tout en acceptant des ruines localises de la structure.
En comparaison avec les mthodes simplifies, les modles de calcul avancs, tant bass sur les proprits
caractristiques des matriaux et sur la ralit des phnomnes physiques, permettent dvaluer de manire plus
prcise le comportement au feu des structures. Ils sappliquent aux effets des incendies naturels et aux effets de
lincendie normalis. Ils peuvent tre utiliss pour lanalyse dlments de constructions isols, lanalyse de sousensembles ou encore pour lanalyse de structure globale incluant linteraction des lments entre eux. En revanche,
lanalyse avec de tels modles reprsente un calcul complexe ncessitant lemploi de logiciels sophistiqus sur
ordinateur. Outre la prise en compte des donnes gomtriques, cette modlisation ncessite de connatre, ou pour le
moins destimer les caractristiques thermo-mcaniques des matriaux concerns. Ceci peut prsenter certains
handicaps lorsque sont mis en uvre des matriaux dont les proprits en fonction de la temprature sont mal
83 / 122
connues, tels que les matriaux de protection. Cest dailleurs lune des raisons pour lesquelles il est ncessaire de
calibrer, partir dessais de rfrence, la validit des modles et donnes utiliss.
84 / 122
Pays
Description
AFCB
Luxembourg
AFCC
Luxembourg
COFIL
Canada
ELIFIR
Belgique
H-Fire
Allemagne
INSTAI
Canada
INSTCO
Canada
POTFIRE
France
RECTST
Canada
WSHAPS
Canada
85 / 122
Pays
Description
ABAQUS
USA
ALGOR
USA
ANSYS
USA
BoFire
Allemagne
CEFICOSS
Belgique
CMPST
France
COSMOS
USA
FASBUS
USA
FIRES-T3
USA
LENAS
France
NASTRAN
USA
SAFIR
Belgique
SISMEF
France
TASEF
Sude
THELMA
UK
VULCAN
UK
86 / 122
III.6. CONCLUSION
De nombreux modles de calcul ont sont actuellement disponibles pour la vrification du comportement au
feu des structures en acier ou mixte acier-bton (voir tableau 4). Avec un degr de prcision limit et
satisfaisant, le comportement mcanique des lments de structure en acier ou mixte acier-bton peut
tre valu laide de mthodes de calculs simplifies, telles que celles proposes dans les parties
feu des Eurocodes. Ces mthodes, bases essentiellement sur des conditions dincendie normalis,
permettent de vrifier rapidement et facilement la stabilit au feu des lments les plus courants (poteau,
poutre, ) vis--vis des exigences descriptives imposes en termes de rsistance au feu dans les
rglementations en scurit incendie actuelles. Le seul cas pour lequel les mthodes de calcul simplifies
peuvent tre galement utilises dans des conditions dincendie rel concerne les lments en acier avec
ou sans protection passive contre lincendie et entirement entours par le feu.
Paralllement, les dveloppements rcents dans le domaine de lingnierie de la scurit incendie offrent
maintenant la possibilit dvaluer, de manire numrique, le comportement au feu dune structure en
acier ou mixte acier bton (considre partiellement dans sa globalit) sous incendie rel ou incendie
normalis.
87 / 122
Type
Valeurs tabules
Mthodes simplifies
Modles de calcul
avancs
Valeurs tabules
Mthodes
simplifies
Modles de calcul
avancs
Elment en acier :
Tirant, Poutre,
Panne, Poteau,
Acier nu ou protg
Non
Oui Mthode de la
temprature critique,
Mthode de vrification
en rsistance
Oui Mthode de
calcul aux lments
finis
Non
Oui
(lments non
protgs)
Oui Mthode de
calcul aux lments
finis
Elment mixte :
Poutre (protge ou
non), Poteau, Dalle,
Plancher
Oui
Oui Mthode de
vrification en
rsistance
Oui Mthode de
calcul aux lments
finis
Non
Non
Oui Mthode de
calcul aux lments
finis
Portique en acier
Non
Oui
Oui Mthode de
calcul aux lments
finis
Non
Non
Oui Mthode de
calcul aux lments
finis
Portique mixte
Non
Non
Oui Mthode de
calcul aux lments
finis
Non
Non
Oui Mthode de
calcul aux lments
finis
Structure en acier
ou mixte
Non
Non
Oui Mthode de
calcul aux lments
finis
Non
Non
Oui Mthode de
calcul aux lments
finis
Elment
individuel
Partie de la
structure
Structure
globale
Tableau 4 : Mthodes de calcul disponibles pour la rsistance au feu des structures en acier et mixte acier-bton
88 / 122
Annexe A
A.1
Description :
Le programme AFCB permet de calculer les moments flchissants maximum des poutres mixtes
temprature ambiante conformment la mthode de lEUROCODE 4 Partie 1-1 (ENV 1994-1-1) et
pour les classes de rsistance au feu normalis R30, R60, R90, R120 et R180 selon lEUROCODE 4
Part 1.2 (ENV 1994-1-2).
Le logiciel a larchitecture suivante :
Donnes dentre :
Dalle : Bien que le programme ne soit pas fait pour ce genre de calcul, il ncessite des informations sur
la dalle (type, largeur efficace, paisseur) afin de pouvoir dterminer sa participation dans la rsistance
de la poutre et les rductions prendre en compte pour la situation dincendie.
Armatures : lutilisateur doit dfinir les armatures dans le bton entre les semelles du profil et les
armatures dans la dalle.
Proprits des matriaux : Dfinit les proprits mcaniques de chaque matriau : limite lastique dans le
profil en acier, rsistance caractristique du bton mesur sur prouvette cylindrique (dalle et enrobage
du profil) et la limite lastique des armature dans le profil et des maillages dans la dalle.
Facteur de scurit des matriaux : Lutilisateur peut choisir les facteurs appliquer la rsistance de
chaque matriau pour les deux cas : charge de service et charge incendie.
c)
Calcul de la rsistance en section : Dtermine les valeurs des moments rsistants positifs et ngatifs
dune section donne.
Dimensionnement sous diffrentes charges : Lutilisateur peut dfinir lui-mme ses charges. La
poutre sera dabord teste froid. Si la rsistance est insuffisante, une modification sur la section
sera demande. Si la rsistance est suffisante froid, le calcul en situation dincendie est effectu. Si
dans ce deuxime cas, la section est insuffisante, le logiciel va essayer dautres combinaisons
darmatures afin den trouver une qui conviendra. La combinaison darmatures est indique dans le
fichier rebars.reb . Lutilisateur peut le modifier.
Dimensionnement pour une rsistance en section minimale donne : similaire au b). La principale
diffrence est que dans ce cas, les valeurs de rsistance voulues dans la section ne sont pas
89 / 122
calcules en utilisant les charges mais doivent tre introduites directement. A nutiliser que si les
moments rsistants peuvent tre calculs la main ou avec un autre programme.
Rsultats :
Le logiciel peut calculer :
Pour des calculs de type b) et c), le logiciel calcule le taux de chargement et les renforts (sils sont
prvus) ;
Dtails : Le calcul complet en conditions de service et pour un feu ISO donn est rappel dans ce fichier.
Lutilisateur peut trouver toutes les donnes dentre par rapport au bord suprieur de la dalle en bton.
De plus, ce fichier texte contient aussi toutes les valeurs rduites des moments positif et ngatif ;
Graphique : Selon le type de calcul, diffrents graphes peuvent tre obtenus : dessin de la section, la
courbe du moment et des graphes pour le calcul de la rsistance en section ;
Formules utilises : La mthode de calcul est indique dans l Eurocode 4 Partie 1.1 et Partie 1.2.
Hypothses :
-
Ce programme traite les poutres sur simples appuis et les poutres continues.
Lannexe H de lENV 1994-1-1 nest pas applique.
Limitations :
-
La vrification des efforts tranchants nest pas effectue dans le programme. Cela doit tre fait
sparment.
Documentation :
Publications et exemples de validations : le logiciel est valid puisque il suit la mthode de lEurocode 4.
Ergonomie :
Interface : Windows
Donnes dentre/sortie : Lutilisateur peut imprimer les donnes de sortie sous forme rsume ou sous
forme exhaustive. Dans cette dernire forme, toutes les donnes dentre et de sortie (rsultats sur les
moments et cisaillement maximums et sur les moments rsistants en section froid et chaud) sont
imprimes ;
Graphique : Le logiciel dessine la section et trace les courbes des moments rsistants pour les moments
positifs et ngatifs froid et chaud.
A.1.3 Conclusions
Ergonomique.
A.2
90 / 122
Version : 3.05 ;
Anne : 2003 ;
Pays : Luxembourg ;
Taille : 2,5 MB ;
Contact information;
Arcelor LCS Research Centre;
66, rue de Luxembourg
L-4221 Esch-sur-Alzette
Tlphone (+352) 5313-3007
Fax (+352) 5313-3095
E-mail : europrofil.dsm@profilarbed.lu
Internet : www.sections.arcelor.com
Description :
Le logiciel AFCC permet de calculer les charge ultimes des poteaux mixtes (poteaux
partiellement enrobs de bton) temprature ambiante selon lEurocode 4 Partie 1.1 (ENV
1994-1-1) et sous feu ISO en classe R30, R60, R90 et R120, selon lEurocode 4 Part 1.2
(ENV 1994-1-2).
Matriaux : Dfinir les proprits mcaniques de chaque matriau : limite lastique dans le profil en acier,
rsistance caractristique du bton mesur sur prouvette cylindrique (dalle et enrobage du profil) et la
limite lastique des armatures dans le profil et des maillages dans la dalle.
Facteur de scurit des matriaux : Lutilisateur peut choisir les facteurs appliquer la rsistance de
chaque matriau pour les deux cas : charge de service et charge incendie
Longueur de flambement : Lutilisateur doit dfinir les longueurs de flambement pour laxe fort et laxe faible
des poteaux mixtes partiellement enrobs de bton en conditions de chargement normal et dincendie.
Excentricit : Excentricits de la charge dans laxe faible et dans laxe fort du profil.
Rsultats : Dans les cinq configurations suivantes : en condition normale dexploitation (20C) et pour les
tempratures atteintes 30, 60, 90 et 120 min (R30 R120), le logiciel calcule :
91 / 122
Dtails : Les dtails complets du calcul (longueur de flambement, charge plastique, charge critique,
lancement rduit, coefficient de flambement) en exploitation normale et pour les diffrentes classes de
rsistance au feu (R30 R120) y sont inscrits. Par ailleurs, la masse linique du poteau est donne,
incluant des informations spares pour le bton et les armatures.
Graphiques : Cette partie du logiciel montre une revue gnrale de la section droite que lutilisateur a dfini
(donnes gomtriques, position des armatures, ).
Formules utilises : la mthode de calcul est celle de lEurocode 4 Parties 1.1 et 1.2 (annexe F). Une fois la
distribution des tempratures connue aprs un temps dexposition au feu donn, le principe de la mthode
consiste dterminer la capacit portante du poteau partir de la rsistance plastique de la section mixte,
corrige par lintermdiaire dun coefficient pour prendre en compte le flambement des poteaux.
Hypothses :
-
Ce logiciel calcule la rsistance dun poteau sous des charges faiblement excentres.
Ce logiciel permet uniquement lanalyse des profils qui comportent deux plans de symtrie et dont la
section droite est constante sur la longueur.
Lannexe H of ENV 1994-1-1 nest pas applique.
Le pourcentage darmature doit respecter les rgles de lENV 1994-1-1, 4.8.3.1 (3e) et 4.8.2.5 (3) et de
lENV 1994-1-2, 4.3.6.2 (2).
Limitations :
-
Documentation :
Publications et exemples de validation : le logiciel est valid puisquil suit la mthode de lEurocode 4.
Ergonomie :
Interface : Windows
Donnes dentre/sortie : Lutilisateur peut imprimer les donnes de sortie sous forme rsume ou sous
forme exhaustive. Sous cette dernire forme, toutes les donnes dentre et de sortie (rsultats sur les
moments et cisaillement maximum et sur les moments rsistants en section froid et chaud) sont
imprims.
A.2.3 Conclusions
Ergonomique
A.3
Elefir
Nom : Elefir
Version : 2.1
Anne : 1998
Domaine dApplication : Modle de rsistance au feu des structures mtalliques Mthode simplifie
Pays : Belgique
92 / 122
Taille : 8 MB
Description :
ELEFIR est un logiciel qui permet de calculer la rsistance au feu dlments individuels en acier, de section
ouverte et chargs selon leurs axes forts.
Diffrentes formes de profils sont disponibles : HD, HE, HL, HP, IPE, UB, UC, W, L.
Deux options dexposition au feu : exposition sur trois faces ou quatre faces.
Les proprits de nombreux matriaux de protection sont disponibles : laine de roche, laine de verre, plaque
de pltre ; lutilisateur peut lui-mme dfinir de nouveaux matriaux de protection
De nombreuses courbes de monte en temprature sont disponibles : courbes ISO, courbe de feu extrieur,
courbe hydrocarbure ou toute autre courbe dfinie par lutilisateur.
Calcul du temps dexposition au feu (temps critique) ncessaire pour que llment atteigne sa temprature
critique ;
Calcul de la temprature atteinte par llment aprs un temps dexposition au feu donn ;
Calcul de la temprature critique de llment et du temps dexposition au feu correspondant pour les profils
comprims, tendus ou comprims et flchis. Pour les lments tendus et comprims, il est galement
possible de calculer la charge maximale que lon peut appliquer llment en fonction de sa temprature
critique ou du temps critique.
93 / 122
Formules utilises :
-
Hypothses :
-
Limitations :
-
Seules les sections ouvertes et comportant deux plans de symtrie peuvent tre analyses ;
Option dexposition lincendie : 3 ou 4 faces ;
Si durant le calcul les sections deviennent des sections de classe 4, le calcul sarrte. Il ne prend pas en
compte la dernire modification de lEN 1993-1-2 qui permet de considrer que la section du profil reste
la mme, quil soit en situation dincendie ou en situation normale dexploitation ;
Documentation :
Manuel de lutilisateur : non disponible, mais peu utile (prise en main facile) ;
Publications et exemples de validation : non disponible. Toutefois, le logiciel peut tre considr comme
valid puisquil suit la mthode de lEurocode 4 ;
Ergonomie :
Interface : Windows
A.3.3 Conclusions
Ergonomique
A.4
H-Fire
Nom : H-Fire ;
Version : 04.1 ;
Anne : 2004 ;
Domaine dApplication : Modle de rsistance au feu des lments mixtes Mthode simplifie ;
Pays : Allemagne ;
Taille : 12.6 MB ;
Prix : Gratuit ;
Contact : www.stahlbau.uni-hannover.de
Description:
94 / 122
H-fire permet le calcul de la rsistance mcanique de lments mixtes (dalles mixtes, poutres mixtes, poutres
mixtes avec profil en acier partiellement enrob de bton et poteaux mixtes partiellement enrobs de bton)
exposs au feu selon les mthodes de calcul simplifies de lEN 1994-1-2. Il inclut plusieurs feuilles de calcul
Excel :
La feuille Excel H-Fire-Slab1-04-1.xls permet de calculer les moments rsistants positif et ngatif des dalles
mixtes exposes au feu. La tle en acier doit tre choisie dans la base de donnes profilenew.mdb ou doit y
tre ajoute. Avant de donner les rsultats, le programme vrifie que le critre d'isolation "I" est bien
respect. Si ce nest pas le cas, le programme envoie un message d'erreur et demande l'utilisateur
dajuster ses donnes d'entre. Les sorties incluent la classe de rsistance au feu normalis demande, le
moment rsistant de calcul en situation dincendie et un rsum des donnes d'entre. Il est possible de
sauver et/ou imprimer les rsultats du calcul.
La feuille H-Fire-Beam1-04-1.xls permet de calculer le moment rsistant des poutres mixtes avec profil en
acier partiellement enrob de bton exposes au feu. Avant de donner les rsultats, le programme vrifie
que le domaine dapplication de la mthode est bien respect. Si ce nest pas le cas, le programme envoie
un message d'erreur et demande l'utilisateur dajuster ses donnes d'entre. Les sorties incluent la classe
de rsistance au feu normalis demande, le moment rsistant de calcul en situation dincendie et un
rsum des donnes d'entre. Il est possible de sauver et/ou imprimer les rsultats du calcul.
La feuille H-Fire-Beam2-04-1.xls permet de calculer le moment rsistant des poutres mixtes avec profil en
acier sans enrobage de bton exposes au feu. Lutilisateur doit indiquer si la dalle est une dalle bton ou
une dalle mixte. Les tempratures du profil en acier et la rsistance de feu demande doivent tre
spcifies. La rsistance au cisaillement et le nombre de connecteurs de cisaillement doivent tre indiqus.
Les sorties incluent les moments rsistants positif et ngatif en situation dincendie et un rsum des
donnes d'entre. Il est possible de sauver et/ou imprimer les rsultats du calcul.
La feuille H-Fire-Column1-04-1.xls permet danalyser les poteaux mixtes partiellement enrobs de bton.
Avant de donner les rsultats, le programme vrifie que le domaine dapplication de la mthode est bien
respect. Si ce nest pas le cas, le programme envoie un message d'erreur et demande l'utilisateur
dajuster ses donnes d'entre. Les sorties incluent la classe de rsistance au feu standard demande, la
rsistance au feu (charge axiale de compression) et un rsum des donnes d'entre. Il est possible de
sauver et/ou imprimer les rsultats du calcul.
Formules utilises : Les calculs sont raliss selon les mthodes de calcul simplifies de lENV 1994-1-2
(Eurocode 4), except pour les dalles pour lesquelles les calculs sont faits selon le prEN 1994-1-2 ;
Hypothses : les mmes que pour celles adoptes dans les mthodes simplifies ;
95 / 122
Documentation :
Manuels techniques : ENV 1994-1-2 (Eurocode 4), et prEN 1994-1-2 pour les dalles disponible sur :
www.stahlbau.uni-hannover.de
Ergonomie :
Donnes dentre/sortie : Le programme prsente la plupart des donnes dentre et toutes les sorties.
A.4.3 Conclusions
Ergonomie : Moyenne
A.5
Potfire
Nom : Potfire ;
Version : 1.11 ;
Anne : 2001 ;
Domaine dApplication : Modle de rsistance au feu des poteaux mixtes Mthode simplifie ;
Pays : France ;
Auteur(s) : Genevive Fouquet, George Tabet, Bin Zhao, Christophe Renaud, Julien Kruppa ;
Organisme(s) : CTICM, CIDECT ;
Matriel requis : Pentium 200 Mhz, W95, CD-Rom, and 24 MB RAM ;
Taille : 15 MB ;
Disponible sur : www.cidect.org;
Contact : www.cidect.org;
Description :
PotFire est un logiciel d'valuation de la stabilit au feu des poteaux non-protgs en
profils creux remplis de bton bas sur la mthode de lannexe G de lENV 1994-1-2. Il
permet de calculer dune part, la dure de rsistance au feu dun poteau sous charge
fixe, et dautre part, dvaluer sa capacit portante maximale pour une dure
dexposition au feu ISO donne. La dtermination de la capacit portante des poteaux
exige le calcul de la distribution et lvolution des tempratures sur la section
transversale des poteaux. Ces tempratures sont obtenues par une procdure
numrique base sur la mthode des diffrences finies.
Formules utilises : Ltat ultime des poteaux aprs un certain temps dexposition au feu correspond la
condition dgalit entre la rsistance plastique de la section mixte (dpendant de ltat des contraintes
normales) et la charge critique eulrienne du poteau, calcule selon le concept du module tangent
dEngesser, appliqu ici aux diffrents matriaux (acier du tube, bton de remplissage et acier darmatures),
eux-mmes des tempratures leves. Une valeur de la rigidit effective du poteau est calcule partir
des relations contrainte-dformation des matriaux, en incrmentant par pas successifs la dformation
96 / 122
longitudinale, suppose constante en tout point de la section. Cette formulation prend en compte
explicitement les non-linarits et les htrognits matrielles, amplifies par llvation des tempratures
Hypothses : Lutilisateur de POTFIRE doit prendre bonne note des dtails de la modlisation du pied et de
la tte dun poteau simple ou des jonctions dun poteau continu, afin de sassurer que les charges sont
appliques de la bonne manire dans le poteau et que les transferts de charges sont maintenus en cas
dincendie.
Limitations : LEurocode 4 Partie 1.2 Annexe G est limit un panel de poteaux (diamtres et longueurs). On
peut prciser encore que la mthode de lAnnexe G sapplique uniquement aux poteaux non protgs,
chargs axialement ou avec une faible excentricit de charge et exposs lincendie conventionnel.
Certaines limites dimensionnelles doivent tre respectes. Elles concernent le type de poteau, (carr ou
rond), les dimensions de la section (140 mm largeur b ou diamtre de la section transversale 400 mm), la
longueur de flambement ( 4,5 m), et le pourcentage darmatures ( 5 %). Enfin, la mthode nautorise pas
lanalyse des poteaux pour une stabilit au feu suprieure 120 minutes.
Section : trois choix : carre rectangulaire ou circulaire. Lutilisateur doit dfinir les dimensions du tube
(diamtre/ largeur et paisseur);
Matriaux : Lutilisateur doit dfinir les proprits mcaniques de chaque matriau : limite lastique du tube
en acier, rsistance la compression du bton et limite lastique des armatures ;
Longueur de flambement : lutilisateur doit dfinir la longueur de flambement du poteau pour la situation
dincendie ;
Rsultats : Charge ultime pour dure donne ou temps de rsistance au feu pour charge donne.
Documentation :
Manuel de lutilisateur : inclus dans le logiciel. Lensemble des quations utilises avec ce modle pour
dcrire le comportement thermique et structurel des poteaux est donn dans lAnnexe 2 du manuel.
Manuels techniques : Bons conseils donns dans lEurocode 4, Part 1-2 et dans le CIDECT Guide de
Conception 4 "Guide de conception des poteaux creux remplis de bton".
97 / 122
Ergonomie :
Interface : Windows.
Donnes dentre/sortie : toutes les donnes saffichent dans une fentre de contrle.
Graphique : aucun.
A.5.3 Conclusions
Ergonomique ;
A.6
COFIL
Nom : COFIL
Version : 1
Anne : 2002
Domaine dApplication : Modle de rsistance au feu des profils creux circulaire remplis de bton lger
Pays : CANADA
Organisme(s) : National Fire Laboratory Institute for Research in Construction, National Research Council of
Canada.
Matriel requis : PC
Taille : 50 Ko
Prix : N/A
Contact : Venkatesh Kodur, Ph.D., P.Eng. Research Officer Institute for Research in Construction National
Research Council of Canada Bldg. M59, 1500 Montreal Road Ottawa, ON, K1A 0R6 Canada Tel: (613) 9939729 Fax: (613) 954-0483
Description :
Donnes dentre :
Matriaux : Dfinir les proprits mcaniques de chaque matriau : limite lastique dans le profil en acier,
rsistance caractristique du bton et la limite lastique des armatures.
Longueur de flambement : lutilisateur doit dfinir les longueurs de flambement en situation dincendie ;
Formules utilises : COFIL calcule lvolution et la distribution des tempratures laide de la mthode des
diffrences finies et, en utilisant la mthode des lments finis, le temps de stabilit au feu (rsistance au
feu) des poteaux pour une exposition au feu standard nord-amricain (ou tout autre feu dfini par
lutilisateur). La rsistance au feu calcule correspond au temps o le poteau nest plus capable de supporter
la charge applique.
98 / 122
Limitations :
o
o
Documentation :
Ergonomie :
Interface : aucune
Donnes dentre/sortie : toutes les donnes saffichent dans une fentre de contrle.
Graphique : aucun
A.7
ABAQUS
Nom : Abaqus
Version : 6.4
Anne : 2003
Pays : USA ;
Langue : Anglais ;
Matriel requis :
Contact : www.abaqus.com ;
Description :
Abaqus est code de calcul gnraliste bas sur la mthode des lments finis.
ABAQUS/Standard : Ce modle est spcialement conu pour la rsolution de problmes linaires et
non-linaires de complexit variable en gomtrie 1D, 2D, Axi ou 3D. Il propose une grande varit de
procdures d'analyses dans le domaine temporel ou frquentiel.
ABAQUS/Explicit : Ce produit fournit une analyse non-linaire, transitoire et dynamique de structures
utilisant la mthode explicite d'intgration du temps. Ses puissantes capacits en contact, sa fiabilit
et son efficacit de calcul sur de gros modles rendent aussi ce produit trs efficace pour des
analyses quasi-statiques impliquant un comportement non-linaire important.
99 / 122
A.7.2 Conclusions
Mthode de calcul :
ABAQUS/Standard : Il fournit un grand nombre de procdures danalyse statique allant de la simple analyse
linaire aux analyses non-linaires (non-linarit matrielle, non linarit gomtrique, contact), rsolues
de manire fiable et efficace (Newton-Raphson, Riks). Le logiciel permet de simuler toute une varit de
phnomnes physiques tels que le transfert thermique, la diffusion de masse, les phnomnes lectriques
(pizo et thermique) et lacoustique, en plus des problmes de mcanique des structures.
ABAQUS excute des analyses de transfert thermique en rgime permanent ou transitoire et des analyses
thermo-mcaniques couples. Il traite galement les problmes d'analyse modale, de rponse transitoire ou
frquentielle et de rponse alatoire ;
ABAQUS/Explicite : il fournit des techniques bases sur les mthode des lments finis qui permettent de
simuler un grande varit de phnomnes quasi-statiques et dynamiques (plus particulirement les
problmes de crash et les problmes avec fortes discontinuits), de manire prcise, robuste et efficace. Il
ne traite pas seulement les problmes de mcanique des structures mais aussi les problmes de couplage
dynamique transitoire temprature-dplacement, dacoustique et de couplage vibro-acoustique.
o
o
o
Documentation :
Lamont, S., The Behaviour of Multi-storey Composite Steel Framed Structures in Response to
Compartment Fires, University of Edinburgh, 2001;
Design Tools of the behaviour of multi-storey steel framed buildings exposed to natural fire conditions
(CARD (2)), Technical Steel Research, Final Report EUR 20953 EN, 2004;
Calculation rules of lightweight steel sections in fire situations, technical steel research, Final report
EUR 21426, 2005;
Ergonomie :
Interface : Windows
Donnes dentre/sortie : Donne dentres prsente dans le fichier input (*.inp) et les sorties dans le fichier
(*.odb).
A.7.3 Conclusions
Pas ergonomique.
100 / 122
A.8
BoFire
Nom : BoFire
Version : 7
Anne : 2004
Pays : Allemagne
Langue : Allemand
Taille : 200 kB
Description :
BoFire est un logiciel bas sur la mthode des lments finis, qui gre les problmes non linaires (non linarit
gomtrique et matrielle) en rgime transitoire. Le programme inclut deux modules de calcul : un modle de
transfert thermique pour permettre le calcul de lvolution et de la distribution des tempratures dans les lments
de structure soumis lincendie (normalis ou non) et un modle de rponse mcanique pour simuler le
comportement au feu de la structure en prenant en compte le changement des proprits matrielles aux
tempratures leves (modle de rponse mcanique). Les proprits thermiques et mcaniques des matriaux
sont celles de lENV 1994-1-2. Les lments de structures en acier, bton et mixte acier-bton peuvent tre
analyss.
Formules utilises : Dans lanalyse mcanique, les quations d'quilibre entre efforts internes et efforts
externes sont rsolus laide dune procdure incrmentale et itrative non-linaire prenant en compte les
effets du second ordre gomtrique et les non-linarits matrielles.
Hypothses :
-
Ce programme permet lanalyse des poutres, des poteaux et des structures planes ;
Limitations :
-
Pas de structure 3D ;
Pas de cisaillement dans la section droite des lments (Hypothses de Bernoulli).
Documentation :
Manuels techniques :
o
101 / 122
Ergonomie :
Interface : Windows
Les paramtres entre/sortie sont donns sous forme de fichier texte. Lapplication HaFront permet de
gnrer le fichier dentre.
Graphique : Le logiciel inclut une librairie DISLIN qui permet de tracer les iso-tempratures en 2D ou 3D ainsi
que les diagrammes des contraintes et des dformations.
A.8.3 Conclusions
Ergonomique.
A.9
LENAS
Nom : Lenas
Version : ---
Anne : 1998
Pays : France
Organisme(s) : CTICM
Langue : Anglais
Taille : ---
Logiciel indisponible
Description :
LENAS-MT (Large Elasto-plastic Numerical Analysis of Structures-Member in Transient state) est un code de
calcul tridimensionnel, bas sur la mthode des lments finis des structures barres. Le logiciel permet de
simuler lvolution des efforts internes et des dformations des structures mtalliques tridimensionnelles
soumises des tempratures voluant avec le temps.
Formules utilises : Lenas permet deffectuer des calculs thermo-lastiques en grande dformation. Les
quations d'quilibre entre efforts internes et efforts externes sont rsolues laide dune procdure
incrmentale et itrative prenant en compte les effets du second ordre gomtrique et les non-linarits
matrielles.
Domaine dutilisation :
o
o
o
o
LENAS est applicable aux lments individuels, aux parties douvrage ainsi quaux structures
compltes. Les lments tudis peuvent tre soumis des tempratures leves homognes en
section ou prsentant des gradients, tant en section (dans les deux directions) que sur la longueur.
Il prend en compte les non-linarits gomtriques (grands dplacements) et matrielles (variation des
lois de comportement lasto-plastique avec la temprature).
LENAS permet de prendre en compte des appuis avec dilatation libre ou limite et des assemblages
rigides, semi-rigides ou articuls.
Les lments traits sont des lments de barre uniaxiaux et tridimensionnels qui possdent 2 nuds
leurs extrmits et sont initialement rectilignes. Ils sont capables de combiner les phnomnes de
flexion, de compression, de torsion et ventuellement de gauchissement.
102 / 122
Les proprits thermomcaniques de l ENV 1994-1-2 (1994) sont implmentes. Les matriaux sont
supposs isotropes avec une fonction de plasticit conforme au critre de Von Mises avec un seul
paramtre dcrouissage.
Afin de permettre le calcul des structures avec nimporte quel type dassemblage, un lment spcifique
(lment de liaison semi-rigide bas sur la mthode de condensation) a t intgr au modle.
Limitations :
Documentation :
Ergonomie :
Interface : aucune
Graphique : aucun
A.9.3 Conclusions
Peu ergonomique ;
A.10 SISMEF
A.10.1 Informations gnrales
Nom : SISMEF
Version : ---
103 / 122
Anne : 1998
Pays : France
Organisme(s) : CTICM
Langue : Francais
Taille : ---
Logiciel indisponible
Description :
SISMEF (Simulation lIncendie des Structures Mixtes par Elments Finis) est un code de calcul bidimensionnel,
bas sur la mthode des lments finis. Il permet de prvoir de faon numrique le comportement des poutres ou
des structures acier, bton ou mixtes (acier-bton) soumis ou non lincendie et ce jusqu la ruine (domaine
plastique), en considrant une connexion complte ou partielle entre la section acier et la section bton pour les
poutres et les poteaux. Il permet de traiter diffrents types d'assemblages, savoir : rigides, simples et semirigides, ainsi que les lments "barre" dans une ou deux dimensions. De plus, SISMEF dtecte les phnomnes
d'instabilit.
Formules utilises : Le modle permet de prendre en compte, dans le cadre dune formulation classique
introduite sous forme incrmentale (la modlisation repose sur les rsultats de la thorie des poutres en
flexion plane respectant lhypothse de Bernoulli), les diverses non-linarits lies au comportement des
matriaux et aux effets du second ordre gomtrique (grands dplacements). Pour les besoins de lanalyse
mcanique, la distribution des tempratures au sein des lments est obtenue, soit laide dun modle
thermique, soit laide des rsultats dessais. Le comportement mcanique de llment est alors dtermin
dans une procdure itrative, prenant en compte lvolution et la distribution des tempratures chaque pas
de temps, ainsi que leurs influences sur les proprits mcaniques des matriaux. Des algorithmes
appropris, de nature incrmentale et itrative, ont t introduits dans le modle pour rsoudre les diverses
non linarits lies au comportement des matriaux, lvolution de ce comportement avec la temprature,
au fluage et aux effets du second ordre gomtrique.
Domaine dutilisation :
SISMEF est applicable aux lments individuels, aux parties douvrage ainsi quaux structures
compltes. Les lments tudis peuvent tre soumis des tempratures leves homognes en
section ou prsentant des gradients, tant en section (dans les deux directions) que sur la longueur.
Le modle permet dtudier des structures dont les lments sont soumis diffrents types de
sollicitations, combinant la flexion et la compression.
SISMEF permet de prendre en compte des appuis avec dilatation libre ou limite et des assemblages
rigides, semi-rigides ou articuls.
Le modle intgre les non linarits lies la plasticit des matriaux, amplifies par laction de la
temprature, incluant les effets du chargement et du dchargement sur la rigidit des lments de
structure.
Le modle intgre les effets des contraintes et des dformations thermiques dues llvation des
tempratures et aux gradients de temprature.
Le modle permet de prendre en compte les effets du second ordre gomtrique, permettant ainsi de
prvoir (en tenant compte des imperfections initiales) les risques dinstabilit par flambement des
poteaux.
Les effets des contraintes rsiduelles et thermiques sont traits explicitement. En prsence dun
champ de temprature, la dformation longitudinale est donne par la somme de quatre termes
indpendants lis respectivement aux contraintes normales, aux dilations thermiques, aux contraintes
rsiduelles et au fluage.
Prise en compte des effets du glissement linterface acier-bton et des effets du second ordre
gomtrique, en prsence dun champ de tempratures voluant de manire quelconque.
Limitations :
104 / 122
o Les effets de torsion ne sont pas considrs et les effets du cisaillement sur lnergie de dformation
de llment sont ngligs.
o Seules les structures planes peuvent tre tudies.
Documentation :
Manuels techniques : SISMEF, logiciel de simulation du comportement mcanique des structures mixtes
(acier-bton) soumises un chauffement - Prsentation et justifications , CTICM, Mai 1998.
Btiment : PARIS : Maison de la Radio (partie IGH et partie ERP) ; Structures en bton arm : poutres,
poteaux, portiques ; Type de feu : feu naturel ; Objectif : Vrification de la rsistance au feu.
Ouvrage : Couverture du Bd priphrique (alle couverte) - Paris Porte des Lilas ; Structures en bton
arm et bton prcontraint : poutres, poteaux, portiques ; Structures mixtes (acier bton arm) :
poutres, portiques ; Type de feu : feu HCM ; Objectif : Vrification de la rsistance au feu.
Nombreuses utilisations pour la vrification au feu de poteaux mixtes en profil creux rempli de bton et
de poutres mixtes.
Le modle mcanique a t confront des essais de poteaux mixtes publis dans la littrature.
Ergonomie :
Interface : aucune
Graphique : aucun
A.10.3 Conclusions
Documentation : dtaille ;
Peu ergonomique ;
A.11 SAFIR
A.11.1 Information gnrale
Nom : SAFIR
Version : SAFIR98a
Anne : 1998
Pays : Belgique
105 / 122
Langue : Anglais
Prix de base : 10 000 EUROS. 40% de rduction pour version rduite (analyse thermique seule ou calcul 2D)
;
Contact : jm.franssen@ULG.AC.BE
Quelques fichiers d'exemples, ainsi que le pr-processeur et le post-processeur graphiques sont disponibles pour
le grand public. La version complte du logiciel est indisponible pour le grand public. Une version de
dmonstration est en libre accs afin de pouvoir valuer le logiciel, pour la recherche ou lducation.
Description :
SAFIR est un logiciel de simulation numrique dvelopp pour lanalyse des structures en acier, bton ou mixte
acier-bton soumises des tempratures leves. Le logiciel bas sur lancien programme par lments finis
CEFICOSS peut tre utilis pour tudier le comportement au feu de structures 2D ou 3D. Il permet de calculer la
distribution et l'volution des tempratures dans la structure exposes au feu et ensuite de simuler le
comportement au feu de la structure (effort, contraintes, dplacements, moments, etc.) jusqu' la ruine. Dans un
souci de simplification, il ny a pas de rel couplage entre le comportement thermique et le comportement
mcanique des structures. Lhistoire des tempratures est calcule indpendamment des charges appliques et
de ltat des contraintes existant dans la structure.
Il rsout les problmes concernant les lments structuraux isols, des sous-ensembles ou des structures
entires.
106 / 122
SAFIR est un programme par lments finis dvelopp spcifiquement pour l'analyse thermique et
structurale des lments de structure exposs au feu. Les pr et post-processeurs rendent le programme
efficace et facile utiliser. SAFIR propose diffrents types d'lment et modles de matriaux pour
diffrentes idalisations et procdures de calcul. Il offre de nombreuses possibilits avec des lments de
type poutre, coque, plaque, solide 2D et 3D et un maillage qui peut tre cartsien ou cylindrique. Le logiciel
inclut deux modules de calcul : un module thermique et un module mcanique.
L'analyse thermique permet dobtenir les champs de tempratures dans des lments de construction. La
prsence ventuelle de matriaux isolants, d'humidit et lvolution des proprits thermiques avec la
temprature peuvent tre considres dans cette analyse. L'exposition au feu normalis ou au feu naturel
dfini par lutilisateur peut tre indique. Le calcul de la distribution des tempratures est transitoire et non
linaire. Diffrents types dlment finis sont disponibles : lments surfaciques 2D (isoparamtriques
linaires 3 ou 4 nuds) ou volumiques 3D ( 6 ou 8 nuds). Les lois de comportement des matriaux
(pour lacier, le bton, le pltre et tous autres matriaux de protection) aux tempratures leves sont bases
sur celles proposes dans les Eurocodes. Elles sont fortement non-linaires et voluent galement en
fonction de la temprature. Les calculs peuvent galement tre faits avec tout autre matriau de proprits
thermiques connues. Le logiciel permet de prendre en compte lvaporation (mais sans migration et sans
transfert de masse) de leau et les phnomnes de rayonnement dans les cavits internes (2D). A lheure
actuelle, la prise en compte des vides internes nest possible que dans le cas des structures planes et les
cavits convexes.
Le champ de temprature issu de l'analyse thermique est utilis comme entre pour l'analyse structurale.
SAFIR peut analyser les sections planes, les structures 2D et 3D, et les lments prcontraints. Les
structures sont discrtises en lments barres de treillis, poutres ( 3 nuds et 7 degrs de liberts,
hypothse de Bernoulli) coques (lments plans triangulaires 3 nuds) ou solides (avec prise en compte
des grands dplacements. Les effets des contraintes thermiques et des empchements de dilatation sont
pris en considration ainsi que les proprits matrielles non-linaires dpendantes de la temprature. De
nombreuses proprits matrielles sont introduites dans le programme SAFIR pour diffrents types
danalyse. Pour les simulations sous temprature ambiante seule, les matriaux valides peuvent avoir des
proprits lastiques, bilinaires ou de Ramberg-Osgood. Pour les tempratures leves, les proprits de
lEurocode sont donnes pour l'acier de construction, lacier des armatures et lacier de prcontrainte. Les
modles de lEurocode et de Schneider sont utiliss pour le bton calcaire et siliceux. Des matriaux valides
tels que l'isolation et le pltre peuvent tre employs pour des structures non porteuses. Diffrentes relations
contrainte-dformation prenant implicitement en compte le fluage thermique sont tablies dans le
programme.
Documentation :
Guide dutilisation : Nwosu, D., Kodur, V., Franssen, JM and Hum, J., User Manual for SAFIR: A Computer
Program for Analysis of Structures at Elevated Temperature Conditions, NRC-CNRC, Int. Report 782, Oct
1999.
o
o
o
A comparison between five structural fire code applied to steel elements, J.M. FRANSSEN, J.B.
SCHLEICH, L.G. CAGOT, D.TALAMONA, B. ZHAO, L. TWILT and K. BOTH, Fourth international
symposium on fire safety science, Ottawa (1994), 1125-1136.
Evaluation of the thermal part of the code SAFIR by comparison with the code Tasef, D. Pintea and
th
J.M. FRANSSEN, Proc. 8 Int. Conf. on Steel structures, Vol.2, (1197), 636-643.
Behaviour of steel frames under fire conditions, Nwosu, D. and Kodur, V.K.R (1998), Canadian journal
of civil Engineering, 26, 156-167.
Appui la normalisation - Volet II. Partie 2: Validation de Programmes Numriques Modlisant les
Structures en cas d'incendie, J.-M. Franssen and J. Unanua, Research Report, Universit de Lige,
2000.
Franssen J.M., Schleich J.B.,.Cajot L.G., Talamona P., Zhao B., Twilt L., Both K., A compararison
between five structural fire model applied to steel elements, Fire Safety Science Proceeding of the
fourth international sumposium, p1125-1136.
Le programme a t valid et employ dans plusieurs tudes de cas effectues dans le pass.
107 / 122
Ergonomie :
Interface : Windows
Graphique : ---
A.11.3 Conclusions
Ergonomique
A.12 VULCAN
A.12.1 Informations gnrales
Nom : Vulcan
Version : ---
Anne : ---
Auteur(s) : ---
Langue : Anglais
Contact : info@vulcan-solutions.com
Description :
Vulcan est un logiciel par lments finis dvelopp l'universit de Sheffield pour lanalyse du comportement
global des btiments (en acier, en bton ou mixte acier-bton) en situation dincendie.
Vulcan permet de simuler le comportement global 3D des btiments charpente mtallique, en bton ou
mixte acier-bton en situation dincendie.
L'analyse considre l'action globale de la charpente et inclut les non-linarits gomtriques et matrielles.
Le logiciel permet dintroduire des assemblages semi-rigides (qui se dgradent avec la temprature) et
permet de prendre en compte une interaction partielle entre le profil acier et la dalle de bton.
108 / 122
Les distributions de temprature sur les lments de structure peuvent tre non-uniformes, permettant de
prendre en compte leffet des dilatations thermiques diffrentielles et une distribution htrogne des
proprits lastiques et non lastiques sur la section transversale des lments.
Un btiment en charpente mixte acier-bton est modlis comme un assemblage dlments de poutre, de
ressorts, de connecteurs en cisaillement et d'lments de coque. Les lments de poutre sont des lments
liniques 3 noeuds. La section transversale de chaque lment est divise en un certain nombre de
segments pour permettre des variations bidimensionnelles de la distribution des tempratures, des
contraintes et des dformations sur la section transversale des lments. Les assemblages entre lments
mtalliques sont modliss avec des lments de type ressort deux noeuds (avec les mmes degrs de
libert nodaux que les lments de poutre) et de longueur nulle. L'lment de connecteur en cisaillement ( 2
noeuds et 3 degrs de libert par nuds : 1 translation et deux rotations) et de longueur nulle permet de
modliser une interaction complte, partielle ou nulle l'interface entre la dalle de bton et le profil en acier.
Les dalles sont modlises partir dlments de coque multicouche et orthotropique, bass sur la thorie
de Mindlin/Reissner. Chaque couche peut avoir une temprature et des proprits matrielles diffrentes.
Le logiciel permet de simuler le comportement en grand dplacement des dalles de plancher exposes au
feu, y compris le flambement d aux empchements de dilatation thermique et leffet membrane.
Documentation :
o
o
o
Huang, Z.; Burgess, I.W.; Plank, R.J., Non-linear modelling of three full-scale structural fire test;
nd
Proceedings of the first international workshop Structures in Fire, 2 edition; Copenhagen; Denmark;
Editor: J.-M. Franssen, pp.53-70, 2000;
Huang, Z.; Burgess, I.W.; Plank, R.J., Non-linear structural modelling of a fire test subjected to high
restraint; in Fire Safety Journal, Elsevier, Vo.36, No.8, pp.795-814, 2001;
Cai, J.; Burgess, I.; Plank, R.J., A generalised steel/reinforced concrete beam-column element model for
fire conditions, in Engineering Structures, Elsevier, Vol.25, pp.817-833, 2003;
SAAB, H.A., NETHERCOT, D.A., 1991, Modelling of Steel Frame Behaviour under Fire Conditions,
Engineering Structures, v.13, n.4, pp. 371-382.
Le programme a t vrifi avec les rsultats exprimentaux des essais au feu de Cardington.
Ergonomie :
Interface : Windows
A.12.3 Conclusions
Documentation : dtaille ;
Ergonomique ;
A.13 ANSYS
A.13.1 Informations gnrales
Nom : ANSYS
Version : 10
Anne : 2006
Pays : USA
109 / 122
Auteur(s) : ---
Langue : Anglais
Prix : ---
Contact : www.ansys.com
Description :
ANSYS est un logiciel de simulation bas sur la mthode des lments finis. Ses domaines d'application sont :
lanalyse structurale (phnomnes quasi statiques ou dynamiques), la thermique, la mcanique des fluides,
llectromagntisme, lacoustique (propagation dondes) et lanalyse en champ coupls : Thermo-mcanique,
Thermo-lectrique, Magnto-thermique, pizolectrique, Electro-magnto-thermo mcanique.
De nombreux modules ANSYS peuvent tre utilises pour le calcul lincendie :
ANSYS Multiphysics : Ce module est le plus complet. Il permet les tudes linaires et non linaires, statiques et
dynamiques, en mcanique et thermique, mcanique des fluides, lectromagntisme, couplages multiphysiques,
etc.
ANSYS Mechanical : Ce module permet les analyses linaires et non linaires, statiques et dynamiques, en
mcanique et thermique.
ANSYS Structural : Ce module permet les analyses linaires et non linaires, statiques et dynamiques en
mcanique uniquement.
ANSYS LS-DYNA : Ce module permet les analyses en dynamique rapide (crash, emboutissage, etc.). Le solveur
explicite LS-DYNA est intgr dans l'environnement pr/post processeur d'ANSYS.
Mcanique des solides : .Phnomne statique ou transitoire. Domaines linaires et non linaires :
gomtrique (grandes dformations, grands dplacements, ), matriaux (la plasticit, hyper-lasticit ou le
fluage), les problmes de contact et lments non linaires. Processus de calcul itratif et incrmental
.Technique du line search . Diffrents types de chargements mcaniques (charges ponctuelles (forces et
moments) charges rparties, pressions, chargements de volume, dplacements imposs ou bloqusforces
suiveuses et non suiveuses).
110 / 122
Elments disponibles: diffrents type dlments disponibles : Elments de barre, poutre et tuyaux ( pipe )
avec catalogue de section standard; Elment de structure 2D (plaques, membranes et coques) et 3D (solides
4, 5, 6 jusqu 20 nuds) ; Formulation linaire et quadratique ; Elments thermiques (Masse, 1D, 2D et
3D) ; Ressorts, masses, amortisseurs, Cbles, Elments de contact (de type gap ou de type surface de
contact 2D et 3D);
Documentation :
Calculation rules of lightweight steel sections in fire situations, technical steel research, Final report
EUR 21426, 2005
Demonstration of real fire tests in car parks and high buildings, final report of ECSC project-EUR
20466 EN 2002, CTICM, Brussels, December 2003
Design Tools of the behaviour of multi-storey steel framed buildings exposed to natural fire conditions
(CARD (2)), Technical Steel Research, Final Report EUR 20953 EN, 2004;
Competitive steel buildings through natural fire safety concept, final report of CEC agreement 7210
SA125, British Steel, March 1999;
o
o
o
o
Btiment : Entrepts de stockage automatis Savoye Logistics, classs sous la rubrique 1510 des
ICPE ; Type de structure : Racks autoporteurs en acier ; Scnarios de feu : feu rel ; Objectif : Analyse
du mode ruine et vrification des critres de ruine vers lextrieur et de ruine en chane.
Btiment : BRON - Universit de Lyon II ERP Type de structure : structure extrieure avec poteaux
et poutres en acier ; Scnarios de feu : tude sous feu rel ; Objectif : Vrification de la stabilit au feu.
Btiment : Chtenay Malabry (rgion parisienne) - Facult de Pharmacie ERP ; Type de structure :
poteaux en acier extrieurs, poutres et plancher mixtes ; Scnarios de feu : feu rel ; Objectif :
Vrification de la stabilit au feu.
Btiment : Parking ouvert - Aix-les-Milles ; Type de structure : poteaux en acier, poutres et plancher
mixtes ; Scnarios de feu : feu rel ; Objectif : vrification de la stabilit au feu.
Btiment : NANTES Palais de justice ; Type de structure : Charpente mtallique ; Scnarios de feu :
feu localis ; Objectif : vrification de la stabilit au feu.
Ergonomie :
Interface : Windows
111 / 122
A.13.3 Conclusions
Ergonomique ;
A.14 DIANA
A.14.1 Informations gnrales
Nom : DIANA
Version : 8.1
Anne : 2002
Pays : Pays-Bas
Auteur(s) : ---
Organisme(s) : TNO
Langue : Anglais
Contact : www.tnodiana.com
Description :
DIANA est un logiciel par lments finis appropri pour la simulation des rponses thermique et mcanique des
structures en acier, en bton et mixte acier-bton en traitant les non-linarits gomtrique et matrielle, les
effets dynamiques et les problmes dpendant du temps et de la temprature. Le programme propose un large
ventail de modles de matriaux pour l'analyse du comportement non linaire du bton qui prsente des
phnomnes de cisaillement au niveau des fissures et des joints, dcrasement, de fissuration et en fissures et
joints, des techniques spciales pour modliser les renforts et les cbles prcontraints, la dtermination et
l'intgration du fluage et schage et solutions avances pour l'analyse des btons durcissant. Le programme
fournit un ventail de modles pour la simulation du comportement non linaire et dpendant du temps des sols,
tels que les sables, les argiles, les vases, des schistes et les roches pour des applications aussi diverses que des
fondations, des remblais, des tunnels, des excavations, des mines et des barrages. Le programme permet
l'analyse des problmes de flambement, de post-flambement et deffondrement des structures mtalliques, de
fatigue dans les structures soudes et rupture.
Mcanique des solides : Domaines linaires et non linaires : gomtrique (grandes dformations, grands
dplacements), matriaux (la plasticit, hyper-lasticit ou le fluage). Processus de calcul itratif et
incrmental avec Formulation Lagrangienne Totale et Lagrangienne actualise ; diffrents schmas
dintgration : Newton-Raphson, Newton-Raphson modifi, Secant stiffness , Constant stiffness
Technique du arc length.
Dynamique : Le logiciel permet deffectuer des analyses modale, harmonique, transitoire linaire et non
linaire, spectrale, ;
Thermique : conduction, convection, radiation, rgime stationnaire ou transitoire, analyse linaire ou non
linaire ; chargements et conditions aux limites : Tempratures imposes, Parois adiabatiques, Flux de
chaleur, Source de chaleur,
112 / 122
Elments disponibles: Diffrents type dlments disponibles : Elments barre de treillis, poutre, plaque,
coques, volumique, axisymtrique, lments de renforts (barre et treillis), lment dcoulement; Formulation
linaire, quadratique ; Des lments spcifiques tels que des lments de renforts noys, des lments
d'interface et de contact, des lments d'coulement et de mlange, des ressorts et des masses sont
galement disponibles ;
Documentation :
Rfrences : Manuel technique de DIANA sous fichiers html (Getting started, Analysis Procedures, Material
Library, Pre and Post-processing, Element Library, Analysis Examples, Cumulative Index).
Design Tools of the behaviour of multi-storey steel framed buildings exposed to natural fire conditions
(CARD (2)), Technical Steel Research, Final Report EUR 20953 EN, 2004;
Franssen J.M., Schleich J.B.,.Cajot L.G., Talamona P., Zhao B., Twilt L., Both K., A compararison
between five structural fire model applied to steel elements, Fire Safety Science Proceeding of the
fourth international sumposium, p1125-1136.
Ergonomie :
Interface : Windows
A.14.3 Conclusions
Ergonomique ;
A.15 TASEF
A.15.1 Informations gnrales
Nom : TASEF
Version : ---
Anne : 2001
Pays : Sude
Auteur(s) : Wickstrm, U.
Langue : Anglais
113 / 122
Description :
TASEF est un logiciel de transfert thermique bas sur la mthode des lments finis permettant de calculer un
champ de temprature bidimensionnel dans les lments de structures exposs au feu.
Le champ de temprature bidimensionnel est bas sur une discrtisation rgulire base de mailles
rectangulaires de la section de llment.
TASEF permet danalyser les lments constitus de diffrents matriaux : chaque maille de la discrtisation
pouvant se voir attribuer un type de matriau diffrent acier, bton ou matriau isolant avec des proprits
voluant en fonction de la temprature.
TASEF permet de prendre en compte la prsence de cavits internes lintrieur desquelles des transferts
thermiques ont lieu par rayonnement et par convection.
TASEF permet de modliser les changes radiatif et convectif avec le milieu ambiant. Les conditions
imposes aux limites peuvent tre, soit des tempratures, soit des flux de chaleur. Des parois adiabatiques
peuvent galement tre prise en compte.
Documentation :
Manuel de lutilisateur : Sterner, E. and Wickstrom, U., TASEF Temperature Analysis of Structures
Exposed to Fire- Users manual. Swedish National Testing and Research Institute (SP), SP Report 1990:05,
Boras 1990;
Manuels techniques :
o Wickstrom, U., TASEF 2 A computer Program for Temperature Analysis of Structures Exposed to
Fire. Lund Institute of Technology, Sweden, Report No. 79-2, Lund 1979.
o Wickstrm, U. and Tuovinen, H., Calculation of Fluid Temperatures in Circular Tubes Using TubeTASEF. SP Swedish Testing and Research Institute, SP REPORT 1997:29, Bors 1997.
o
o
Wickstrm, U. and Plsson, J. Scheme for Verification of Computer Codes for Calculating
Temperature in Fire Exposed Structures, SP Swedish Testing and Research Institute, SP REPORT
1999:36, Bors 1999.
Thomas, G. C., Fire Resistance of Light Timber Framed Walls and Floors, Fire Engineering Research
Report 97/7 (Chapter 3), School of Engineering, University of Canterbury, New Zealand, 1997.
Gerlich, J. T. C., Collier, P. C. R. and Buchanan, A. H., Design of Light Steel-framed Walls for Fire
Resistance, Fire and Materials, vol. 20, pp 79 96, (1996).
Ergonomie :
Interface : aucune
Graphique : aucun
A.15.3 Conclusions
Peu ergonomique ;
114 / 122
115 / 122
Depuis une vingtaine dannes, il y a une tendance considrer le verre comme un matriau de structure, au mme titre ou en
substitution de matriaux plus couramment utiliss pour de telles applications : bois, acier, bton, pierre Cest un nouveau
dveloppement de ce matriau conduisant par exemple lutilisation de poteaux ou de poutres en verre et ouvrant ainsi de
nouvelles possibilits aux architectes et ingnieurs afin de concevoir des structures transparentes. Mais cette demande
grandissante du verre en tant que matriau de structure confronte les concepteurs un problme : la plupart des standards
nationaux de construction ne couvre pas les structures en verre porteur et il nexiste actuellement pas de rgles de calcul pour
la vrification de la rsistance au feu des structures en verre. Une mthode de dimensionnement froid est prsente dans le
CPT du GS2 publi dans le cahier du CSTB n 3448 : Dalles de plancher et marches descalier en verre conditions
gnrales de conception, fabrication et mise en uvre .
IV 1.
BIBLIOGRAPHIE
Cette premire partie a pour objectif de recenser les connaissances relatives au matriau verre en tant qulment de structure
et au comportement au feu de ce mme matriau.
Figure IV-1 : Variation du volume spcifique avec la temprature pour deux vitesses de refroidissement
Quand le verre est chauff jusqu la temprature T2 (point A de la figure IV-1) en dessous de Tg (temprature de
transformation), la structure du verre continue voluer vers un tat stable : cet tat correspond au verre liquide (point B de la
figure IV-1). Cette volution est appele relaxation structurelle directe reprsente par le chemin AB.
116 / 122
Cette volution est trs importante si le verre se trouve une temprature proche de Tg. Par contre, si la temprature est
largement en dessous de Tg, ce phnomne devient tellement lent quil nest pas observable. Si le verre est chauff brutalement
jusqu T1 au-del de Tg, le volume du verre volue vers une valeur largement en dessous de la droite liquide (point C de la
figure IV-1). Lvolution suit le prolongement de la droite verre . Cette configuration instable volue vers un tat stable : le
verre liquide (point D de la figure IV-1). Cette volution est appele relaxation structurelle inverse.
Quand sa structure est stabilise, le verre est un matriau thermorhologiquement simple. Un matriau viscolastique est dit
thermorhologiquement simple si un changement de temprature peut tre pris en compte par un changement de lchelle de
temps ; le matriau prsente les mmes comportements diffrentes tempratures et seule la vitesse des mcanismes
microstructuraux de rorganisation (temps moyen de relaxation par exemple) varie avec la temprature.
= 0,975 + 8,58.10 4 T
avec T en C.
Conductivit
(W/m.K)
1.8
1.7
1.6
1.5
1.4
1.3
1.2
1.1
1
0.9
Temprature (C)
0
200
400
600
800
1000
1,8.10 7
T2
C p , s = 893 + 0,4T
C p ,l = 1433 + 6,5.10 3 T
Chaleur
spcifique
(J/kg.K)
1500
1450
1400
1350
1300
1250
1200
1150
1100
1050
1000
273
Temprature (K)
473
673
873
1073
117 / 122
IV 2.
Il nexiste actuellement pas de rgles pour la vrification de la rsistance au feu des structures en verre. Les mthodes utilises
sont donc essentiellement exprimentales. De nombreux essais de rsistance au feu ont t raliss sur des cloisons
composes de panneaux de verre.
Il existe cependant des exemples de calculs thermomcaniques concernant le verre. Ces calculs, effectus par lments finis,
sont destins simuler le comportement du verre au cours de lopration de trempe (phnomne faisant intervenir la fois le
temps et la temprature). Lobjectif de la simulation de la trempe thermique du verre est de dterminer les contraintes
transitoires et rsiduelles dues ce traitement. Ce traitement thermique provoque en effet lapparition de contraintes au cours
du refroidissement, en raison des forts gradients thermiques dans llment, et de contraintes rsiduelles car le verre a un
comportement viscolastique variable avec la temprature. Le calcul des contraintes transitoires permet dvaluer les risques
de rupture du verre en cours de trempe et de pouvoir ventuellement quantifier lvolution de dfauts prexistants. Les
contraintes rsiduelles donnent ltat de prcontrainte du verre la fin du traitement. La connaissance de ces contraintes est
indispensable pour prdire la rsistance du verre tremp.
Un modle de comportement thermomcanique incluant les diffrents phnomnes de relaxation des contraintes fait
actuellement rfrence : le modle de Narayanaswamy. Le verre est dcrit en tant que matriau viscolastique. De plus, son
comportement variable avec le temps est aussi fonction de la temprature. Ce modle prend donc en compte la relaxation des
contraintes mais aussi la relaxation structurelle (influence de ltat de la structure du verre sur son comportement).
M v (t ) =
V (t ) V , 2
V 0, 2 V , 2
T f T2
T1 T2
La dpendance en temprature de la fonction rponse est obtenue au moyen du temps rduit . MV(t) peut tre considre
comme la fonction de relaxation volumique qui permet de dfinir les variations de la temprature fictive :
T f (t ) = T (t ) M V [ (t ) (t ' )]
0
dT (t ' )
dt '
dt '
(t , T ) =
0
ref
(T , t ' )
ln( (T )) =
H
R
1
1
ref T
H = nergie dactivation ;
R = constante des gaz parfaits.
118 / 122
de
La figure ci-dessous reprsente lvolution des proprits (proprits dpendantes de la temprature) du verre en fonction de la
temprature :
Figure IV-4 : Variation dune proprit p du verre avec la temprature pour une vitesse de chargement thermique
spcifie
Lquation gnrale pour la dtermination de la proprit p une temprature quelconque T de la zone de transition est la
suivante :
o T0 est la temprature initiale (temprature ambiante) et et sont respectivement les pentes des droites AB et CD.
Cette quation montre que la valeur de p recherche peut tre obtenue en considrant le verre comme liquide au-dessus de la
temprature fictive Tf et comme un solide froid en dessous de cette temprature.
Le modle de Narayanaswamy permet de prdire correctement les valeurs des contraintes transitoires et rsiduelles et permet
de connatre ltat de contrainte final en tout point du composant en verre tremp. Par ailleurs, ce modle prend en compte les
mmes proprits du verre que celles ncessaires un calcul de rsistance au feu et pourrait ainsi permettre, aprs
implmentation dans un code aux lments finis, deffectuer un tel calcul.
d
d
= C ( 1 2 ) + 2
sin cot
dz
dz
avec tel que :
d
d
=
cos
dz
dz
1 et 2 sont les contraintes principales secondaires et est langle des directions principales secondaires avec laxe du
polariseur.
Grce ce matriel, il est possible de suivre les isochromes au cours des essais de chargement. Les images obtenues servent
de moyen de validation du calcul numrique aux lments finis, ralis pour simuler ces essais.
Ainsi, en ralisant des images photolastiques sur des structures existantes, on peut, grce un programme danalyse
dimages, raliser une analyse inverse de ltat de contraintes dans le matriau. Cette mthode permet de vrifier in situ ltat
de la structure en verre, notamment pour chercher savoir si la surface du verre est dcomprime ou non.
119 / 122
Rfrences
[1]
H. Carr
Etude du comportement la rupture dun matriau fragile prcontraint : le
Thse soutenue le 11 dcembre 1996
[2]
[3]
[4]
B. Beer
Structural glass engineering a review of project specific testing
Glass Processing Days 2005
[5]
O.S. Narayanaswamy
A model of structural relaxation in glass
Journal of the American Ceramic Society, vol.54, p. 491-498, 1971
[6]
O.S. Narayanaswamy
Evolution of glass tempering models
Glass Processing Days, 18-21 juin 2001
[7]
L. Daudeville, F. Bernard, R. Gy
La birfringence du verre : une proprit remarquable pour lingnieur
Photonique (numro spcial Photomcanique), p. 26-29
[8]
L. Daudeville, H. Carr
Dtermination de la rsistance du verre tremp
Verre, vol.3, p. 54-64
[9]
L. Daudeville, F. Bernard, R. Gy
Une mesure des contraintes pour le contrle en service des structures en verre
Revue I2M, vol.5, p. 1-16, mars 2005
verre tremp
120 / 122
121 / 122
Type
Valeurs
tabules
Mthodes
simplifies
Modles
de calcul
avancs
Non
Oui Mthode
de la
temprature
critique,
Mthode de
vrification en
rsistance
Elments en
bton :
Poutre,
Poteau,
Dalle, Mur,
Elments
mixtes :
Poutre
(protge ou
non),
Poteau,
Dalle,
Plancher,
Mthodes
simplifies
Modles
de calcul
avancs
Oui
Mthode
E.F.
Non
Oui
(lments
non
protgs)
Oui
Mthode
E.F.
Oui
Oui Mthode
de lisotherme
500C, Mthode
par zones
Mthode de
vrification en
rsistance
Oui
Mthode
E.F.
Non
Non
Oui
Mthode
E.F.
Oui
Oui Mthode
de vrification
en rsistance
Oui
Mthode
E.F.
Non
Non
Oui
Mthode
E.F.
Elments
Bois :
Poutre,
Poteau,
Non
Oui Mthode
de la section
rduite
Mthode des
proprits
rduites
Oui
Mthode
E.F.
Non
Non
Oui
Mthode
E.F.
Elments
Verre
Non
Non
Non
Non
Non
Non
Partie de
la
structure
Portique en
acier, bton
au mixte
acier-bton
Non
Oui (si
disponible)
Oui
Mthode
E.F.
Non
Non
Oui
Mthode
E.F.
Structure
globale
Structure en
acier, en
bton, en
bois ou
mixte acierbton
Non
Non
Oui
Mthode
E.F.
Non
Non
Oui
Mthode
E.F.
Elment en
acier :
Tirant,
Poutre,
Panne,
Poteau,
Acier nu ou
protg,
Elment
individuel
122 / 122