Code_Aster ®
Version
5.0
Titre :
Identification du modèle de Weibull
Date :
30/11/00
Auteur(s) :
R. MASSON, W. LEFEVRE
Clé :
R7.02.09-A
Page :
1/8
Organisme(s) : EDF/RNE/MTC
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
Document : R7.02.09

Identification du modèle de Weibull
Résumé
On aborde ici le problème de l'identification des paramètres du modèle de WEIBULL sur un échantillon d'essais
représentatif du comportement à rupture d'un matériau fragile (typiquement, acier ferritique à basse
température). La méthode de régression linéaire et la méthode du maximum de vraisemblance sont les deux
méthodes retenues. On en détaille le principe ainsi que les méthodes de résolution associées, s'appuyant dans
les deux cas sur un procédé itératif. Enfin, on montre leur extension au cas où l'un des deux paramètres de ce
modèle (la contrainte de clivage) dépend de la température.
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/00/006/A

Code_Aster ®
Version
5.0
Titre :
Identification du modèle de Weibull
Date :
30/11/00
Auteur(s) :
R. MASSON, W. LEFEVRE
Clé :
R7.02.09-A
Page :
2/8
Table des matières
1 Introduction ............................................................................................................................................ 3
2 Rappels .................................................................................................................................................. 3
2.1 Le modèle de WEIBULL.................................................................................................................. 3
2.2 Identification des paramètres .......................................................................................................... 3
3 Méthode de la régression linéaire .......................................................................................................... 4
3.1.1 Principe .................................................................................................................................. 4
3.1.2 Résolution .............................................................................................................................. 5
4 Méthode du maximum de vraisemblance .............................................................................................. 6
4.1 Principe............................................................................................................................................ 6
4.2 Résolution........................................................................................................................................ 6
5 Dépendance des paramètres avec la température................................................................................ 7
5.1 Régression linéaire.......................................................................................................................... 7
5.2 Maximum de vraisemblance............................................................................................................ 8
6 Conclusion ............................................................................................................................................. 8
7 Bibliographie .......................................................................................................................................... 8
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/00/006/A

Code_Aster ®
Version
5.0
Titre :
Identification du modèle de Weibull
Date :
30/11/00
Auteur(s) :
R. MASSON, W. LEFEVRE
Clé :
R7.02.09-A
Page :
3/8
1 Introduction
Lorsqu'elles font appel au modèle de WEIBULL (cf. POST_ELEM [U4.81.22]), les étude de modélisation
de la rupture fragile des aciers nécessitent en général une identification préalable des paramètres de
ce modèle. Afin d'éviter une laborieuse identification « à la main » de ces paramètres qui nécessiterait
de relancer itérativement l'opération POST_ELEM avec l'option WEIBULL, une procédure de recalage
automatique a été implantée dans Code_Aster.
Dans ce document, on rappelle brièvement les équations du modèle de WEIBULL puis on définit le
problème d'identification posé. On décrit alors le principe des deux méthodes de résolution retenues
(régression linéaire et maximum de vraisemblance) en incluant le cas où l'un des deux paramètres du
modèle dépend de la température.
2 Rappels
2.1
Le modèle de WEIBULL
On considère une structure de comportement élastoplastique soumise à une sollicitation
thermomécanique. On suppose que la probabilité de rupture cumulée de cette structure suit la loi de
WEIBULL [bib1] à deux paramètres suivante :
m

P
w
1


f ( w ) = - exp -

éq 2.1-1
u



expression dans laquelle le module de WEIBULL m > 1 décrit la queue de la distribution statistique
des tailles des défauts à l'origine du clivage, u est la contrainte de clivage et w est la contrainte de
WEIBULL qui dépend de l'histoire du champ de contrainte principale dans la zone plastifiée de la
structure. Par exemple, dans le cas d'un trajet de chargement monotone, elle s'écrit :
V

p m
p
w = m ( I )
.
éq 2.1-2
V
p
0
p
La sommation porte sur les volumes de matière Vp plastifiés, I désignant la contrainte principale
maximale dans chacun de ces volumes (V0 est un volume caractéristique du matériau).
2.2
Identification des paramètres
De façon très générale, on considère une base expérimentale constituée d'essais de différentes
natures (type 1, 2, ..., n), chaque type d'essai étant réalisé n j fois de sorte que le nombre total
d'essais s'élève à :
j=n
N = n j .
j=1
Cette base expérimentale pourrait par exemple être constituée d'essais sur éprouvettes axisymétriques
entaillées de rayons d'entaille différents menés à différentes températures. Compte - tenu de la nature
aléatoire des propriétés à rupture du matériau considéré, cette base ne constitue qu'un échantillon.
Plus le nombre de ces échantillons sera important, plus il sera représentatif du comportement du
matériau considéré.
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/00/006/A

Code_Aster ®
Version
5.0
Titre :
Identification du modèle de Weibull
Date :
30/11/00
Auteur(s) :
R. MASSON, W. LEFEVRE
Clé :
R7.02.09-A
Page :
4/8
Parmi les différentes méthodes d'identification proposées dans la littérature (voir par exemple [bib2]),
nous en retenons deux : la méthode de régression linéaire, souvent utilisée, ainsi que celle du
maximum de vraisemblance préconisée par l' « European Structural Integrity Society (ESIS) » [bib3].
Remarque :
Une étude systématique comparative des résultats donnés par ces deux méthodes [bib2] en
fonction du nombre d'échantillon prélevé aléatoirement sur une distribution théorique a montré
que la méthode du maximum de vraisemblance conduit à une meilleure estimation des
paramètres du modèle de WEIBULL. La méthode de régression linéaire demeurant
néanmoins très utilisée, nous l'avons intégrée à nos développements.

Dans les deux méthodes de recalage retenues, on effectue un premier calcul des contraintes de
WEIBULL avec un jeu de paramètre donné (typiquement, m=20, =3000 MPa). On classe ces N
u
essais à l'aide de leur contrainte de WEIBULL atteinte au moment de la rupture. On dispose donc
d'une liste croissante de contraintes de WEIBULL ( 1
i
N
w , ..., w , ..., w ), telle que pour chaque (i) , le
nombre d'éprouvettes rompues avec une contrainte de WEIBULL inférieure ou égale à i
i
w est nw (en
général ni = i
i
w
). Parmi les différents estimateurs possibles de la probabilité de rupture cumulée Pf
i
correspondant à i
i
w [bib2], nous choisissons celui généralement recommandé : Pf =
.
N + 1
Remarque :
Dans le cas particulier où la contrainte de WEIBULL dépend de la température, le
classement précédent doit se faire température par température, chaque température
correspondant à une loi statistique différente. L'estimateur de la probabilité de rupture

i
précédent devient donc : Pif =
, si l'éprouvette (i) a été cassée à la
NT + 1
température T, pour laquelle il y a eu NT essais.
Les deux méthodes de recalage retenues sont valables tant que [éq 2.1-1] reste vraie. Si l'identification
est menée sur des résultats d'essais anisothermes alors que la contrainte de clivage est supposée
dépendre de la température, cette condition n'est plus vérifiée (cf. POST_ELEM [U4.81.22]). Dans ce
cas particulier, on ne pourra donc pas appliquer les développements qui suivent.
3
Méthode de la régression linéaire
3.1.1 Principe
L'écart théorie-expérience est mesuré par l'expression :
2

1
1

LogLog(
)

- LogLog(
)

éq 3.1.1-1
1

- Pi
1 - P
i

i
f
f (W )
(« Log » désigne le logarithme népérien). On veut minimiser cet écart par rapport à (m, u ).
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/00/006/A

Code_Aster ®
Version
5.0
Titre :
Identification du modèle de Weibull
Date :
30/11/00
Auteur(s) :
R. MASSON, W. LEFEVRE
Clé :
R7.02.09-A
Page :
5/8
3.1.2 Résolution
La méthode de recalage habituellement utilisée s'appuie sur des régressions linéaires successives : à
l'itération (k), les valeurs ( mk ,u(k) ) du module et de la contrainte de clivage sont connus. Il est donc
possible, avec ces valeurs, de calculer les contraintes de WEIBULL iW(k) aux différents instants de
rupture grâce à [éq 2.1-1]. On classe alors ces nouvelles contraintes de WEIBULL par amplitude
croissante et on en déduit les nouvelles estimations des probabilité de rupture Pif (k) à l'itération (k).
Pour ces valeurs de contraintes de WEIBULL fixées, la minimisation de [éq 3.1.1-1] se ramène à une
1
simple régression linéaire sur le nuage de points ( Log
i
(W k )
( ) , LogLog(
) ) puisque si on
1- Pif (k)
1
reporte LogLog(
)
( ) , on obtient une droite de pente m qui coupe l'axe
1-
en fonction de Log
P
W
f
des abscisses en ( Log( ) ). Les nouvelles valeurs ( m
,
) de ces paramètres sont donc
u
k +1
u(k +1)
données par (annulation des dérivées partielles de [éq 3.1.1-1] par rapport à chaque paramètre) :
1 X Y
( )
( ) - Y
X
N
i k
j k
i(k )
i(k )
i, j
i
mk+1 = 1
éq 3.1.2-1
X X
( )
( ) - X
2
N
i k
j k
i(k )
i, j
i
1
1




u k
exp
X
+1 =
i k - iY
(
)
( )
(k

,
éq 3.1.2-2
N

m
)
i
i

1
avec X
= Log i
i(k )
(W(k) ) et Y
= LogLog
i(k )
(
) .
1 - Pif (k)
On répète ces itérations tant que la différence entre les jeux de paramètre obtenus aux itérations (k) et
(k+1) est significative (typiquement, cinq itérations). La mesure de cet écart est donnée
m


1 - m

u k +1 - u k
par : Max
k +
k

,
(
)
( )
m

k

.
u(k

)


Remarque :
Si m est fixé, u(k+ )1 est toujours donnée par [éq 3.1.2-2]. En revanche, si u est fixé,
X
Y
i(k) i(k)
m
i
k +1 n'est plus donné par [eq 3.1.2-1] mais : mk +1 =
.
X

2
( ) - log( )
X
i k
u i(k)
i
i
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/00/006/A

Code_Aster ®
Version
5.0
Titre :
Identification du modèle de Weibull
Date :
30/11/00
Auteur(s) :
R. MASSON, W. LEFEVRE
Clé :
R7.02.09-A
Page :
6/8
4
Méthode du maximum de vraisemblance
4.1 Principe
Notons p f (w ) la densité de probabilité associée à la probabilité de rupture cumulée Pf (w ) :
m-1
m
m

p
w
w


f ( w ) =


exp -

u u
u



La quantité p ( )d
f
w
W est égale à la probabilité de casser une éprouvette soumise à une
sollicitation correspondant à une contrainte de WEIBULL comprise dans l'intervalle [W ,W + dW ].
La probabilité pour que toutes les éprouvettes de la base aient cassé s'élève donc à :
p(m, ) d = p
i
( ) d
u
w
f
W
w ,
éq 4.1-1
i
p étant la fonction de vraisemblance. La méthode du maximum de vraisemblance consiste alors à
choisir les paramètres du modèle de WEIBULL de façon à ce que la fonction de vraisemblance définie
par [éq 4.1-1] (en pratique plutôt son logarithme népérien) soit maximale.
4.2 Résolution
On utilise à nouveau un procédé itératif. Là encore, à l'itération (k), ( m
i
k ,u(k ) ) ainsi que les W (k )
sont connus. Pour ces valeurs de contraintes de WEIBULL fixées, la maximisation de Log( p) conduit
à un nouveau couple ( mk+1,u(k+1) ) donné par :
i= N
i
mk
(
) +1 Log
i
(
)

N
W (k )
W (k
i= N
)
f (m
)
1
1 =
+
Log
i
(
)

- N i=
k +
= 0 éq 4.2-1
m
W (k )
i= N
k +1
i=1
i
mk
(
+

1
W(k )
)
i=1
i= N
1

m
=

( i
)mk
k
+
+
1
1
.
éq 4.2-2
1

( k + )
W ( k )
N i=1
A chaque pas, la résolution de [éq 4.2-1] peut être réalisée à l'aide de la méthode de Newton, le
gradient de f (m) étant donné par :

i= N
i= N
i= N



2

i
m
( ) Log2 i
i
m
i
m
( ) ( ) - ( ) Log i

W
W
W
W
( W )
df
1 =1
=1
=1

(m) = - N
i
i
i

+
.
dm
m2
i= N

2

i
m

(

W )


i=1


Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/00/006/A

Code_Aster ®
Version
5.0
Titre :
Identification du modèle de Weibull
Date :
30/11/00
Auteur(s) :
R. MASSON, W. LEFEVRE
Clé :
R7.02.09-A
Page :
7/8
Remarque :
Si m est fixé,
est donné par [4.2-2]. En revanche, si est fixé, m
n'est plus
u(k + )
1
u
k +1
solution de [4.2-1] mais de :
N
i
i
i= N


f m
=
+
Log
W ( k )
W ( k )
mk
(
)
(
)(1
1
0
1

- (
) + ) =
k +
m


.
k
i
+
=
1
1
u
u
Cette équation peut être à nouveau résolue à l'aide de la méthode de Newton, le gradient
étant à présent donné par :

df
N
i
i
i= N

(m
W
m
) = -
2
2 -
(
) Log
W
(
)

dm
m

.
i=1
u
u
5
Dépendance des paramètres avec la température
Si on souhaite caler indépendamment les deux paramètres température par température, il suffit de
décomposer la base d'essais en autant de sous - bases par température et d'appliquer à chacune de
ces sous-bases les méthodes précédentes. Si, en revanche, on souhaite seulement faire varier la
contrainte de clivage avec la température, on procède de la façon suivante.
u
5.1 Régression
linéaire
L'estimation des probabilités de rupture étant à présent effectuée température par température
(cf. remarque du [§2.2]), il suffit de caler la contrainte de clivage sur chaque nuage de points associé
aux différentes températures (T). L'équation [éq 3.1.2-2] devient donc :
1
1





u k 1 = exp
X
+
i k - iY
(
)
( )
(k
N
)


T
m
i T

i T


( NT désignant le nombre d'essais pour la sous - base correspondant à la température (T)), le module
de WEIBULL étant donné par :



1 X Y
( )
( )
Y X
N
i k
j k
i(k )
i(k
-
)
T
T i T
, j T


i
mk+1 =
.



1 X X
( )
( )
X 2
N
i k
j k
i(k
-
)
T
T i T
, j T


i
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/00/006/A

Code_Aster ®
Version
5.0
Titre :
Identification du modèle de Weibull
Date :
30/11/00
Auteur(s) :
R. MASSON, W. LEFEVRE
Clé :
R7.02.09-A
Page :
8/8
5.2
Maximum de vraisemblance
La contrainte de clivage est donnée pour chaque température (T) considérée par :

1
i
m
m
+1
u k
(T)
1
= +1
(W(k (T k
k
(
)
))
+
,
N
)
T i T

m étant solution de :
k +1
i
mk
(
) +1 Log
i
(
)

N
W (k )
W (k
i= N
)
f (m
)
1 =
+
Log
i
(
)

-
N i T

k +

= 0.
m
W (k )
T
i
m
k
k
+
+
1
1
i=1
T
(
W(k ))
i T

6 Conclusion
La commande RECA_WEIBULL du Code_Aster permet d'effectuer le calage des paramètres du modèle
de WEIBULL [U4.82.06].
L'utilisateur donne en entrée de cette commande les concepts résultats associés aux différents calculs
non linéaires effectués. L'éventuelle dépendance de la contrainte de clivage avec la température est
implicitement spécifiée lorsque des températures différentes sont associées à chacun de ces concepts
résultats (si toutes ces températures sont identiques ou si elles ne sont pas spécifiées, il n'y a pas de
dépendance avec la température de ce paramètre).
L'utilisateur peut effectuer ce recalage par la méthode du maximum de vraisemblance (METHODE :
`MAXI_VRAI') ou celle de la régression linéaire (METHODE : `REGR_LIN').
Les grandeurs déterminées par la commande RECA_WEIBULL sont reportées dans une table dans
laquelle on trouve la valeur des paramètres identifiés, les probabilités de rupture estimées à partir des
résultats expérimentaux ainsi que les probabilités de rupture théorique calculées avec les paramètres
identifiés.
7 Bibliographie
[1]
F. BEREMIN, « A local criterion for cleavage fracture of a nuclear pressure vessel steel »,
Metall. Trans. 14A, pp 2277-2287, 1981.
[2]
A. KHALILI, K. KROMP, « Statistical properties of weibull estimators », Journal of Material
Science, 26, pp 6741-6752, 1991.
[3]
ESIS, TC 1.1 on « Local Approach », Procedure to measure and calculate local approach
criteria using notched tensile specimens », P6, 1998.
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/00/006/A

Document Outline