background image
Code_Aster
®
Version
5.0
Titrate:
Identification of the model of Weibull
Date:
30/11/00
Author (S):
R. MASSON, W. LEFEVRE
Key:
R7.02.09-A
Page:
1/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/00/006/A
Organization (S):
EDF/RNE/MTC
Manual of Reference
R7.02 booklet: Breaking process
Document: R7.02.09
Identification of the model of Weibull
Summary
One tackles here the problem of the identification of the parameters of the model of WEIBULL on a sample of tests
representative of behavior with rupture of a fragile material (typically, ferritic steel with low
temperature). The method of regression linear and the method of the maximum of probability are both
adopted methods. One details of it the principle as well as the associated methods of resolution, resting in
two cases on an iterative process. Lastly, one shows their extension if one of the two parameters of it
model (the stress of cleavage) depends on the temperature.
background image
Code_Aster
®
Version
5.0
Titrate:
Identification of the model of Weibull
Date:
30/11/00
Author (S):
R. MASSON, W. LEFEVRE
Key:
R7.02.09-A
Page:
2/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/00/006/A
Contents
1 Introduction ............................................................................................................................................ 3
2 Recalls .................................................................................................................................................. 3
2.1 The model of WEIBULL .................................................................................................................. 3
2.2 Identification of the parameters .......................................................................................................... 3
3 Method of the linear regression .......................................................................................................... 4
3.1.1 Principle .................................................................................................................................. 4
3.1.2 Resolution .............................................................................................................................. 5
4 Method of the maximum of probability .............................................................................................. 6
4.1 Principle ............................................................................................................................................ 6
4.2 ........................................................................................................................................ Resolution 6
5 Dependence of the parameters with the temperature ................................................................................ 7
5.1 Linear regression .......................................................................................................................... 7
5.2 Maximum of probability ............................................................................................................ 8
6 Conclusion ............................................................................................................................................. 8
7 Bibliography .......................................................................................................................................... 8
background image
Code_Aster
®
Version
5.0
Titrate:
Identification of the model of Weibull
Date:
30/11/00
Author (S):
R. MASSON, W. LEFEVRE
Key:
R7.02.09-A
Page:
3/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/00/006/A
1 Introduction
When they call upon the model of WEIBULL (cf.
POST_ELEM
[U4.81.22]), the study of modeling
brittle fracture of steels in general require a preliminary identification of the parameters of
this model. In order to avoid a hard identification “with the hand” of these parameters which would require
to start again the operation repeatedly
POST_ELEM
with the option
WEIBULL
, a procedure of retiming
automatic was established in Code_Aster.
In this document, one briefly recalls the equations of the model of WEIBULL then one defines it
problem of identification posed. One then describes the principle of the two adopted methods of resolution
(linear and maximum regression of probability) by including the case where one of the two parameters of
model depends on the temperature.
2 Recalls
2.1
The model of WEIBULL
One considers a structure of behavior elastoplastic subjected to a stress
thermomechanics. It is supposed that the probability of cumulated rupture of this structure follows the law of
WEIBULL [bib1] with two parameters following:
P
F
W
W
U
m
(
)
exp
= -
-


1
éq 2.1-1
expression in which the module of WEIBULL
m
> 1 described the tail of the statistical distribution
sizes of the defects at the origin of cleavage,
U
is the stress of cleavage and
W
is the stress of
WEIBULL which depends on the history of the main stress field in the area plasticized on
structure. For example, in the case of a monotonous way of loading, it is written:
()
W
IP
m
p
p
m
V
V
=
0
.
éq 2.1-2
The summation relates to volumes of matter
V
p
plasticized,
IP
indicating the main stress
maximum in each one of these volumes (
V
0
is a volume characteristic of material).
2.2
Identification of the parameters
In a very general way, one considers an experimental base made up of tests the different ones
natures (type 1, 2,…, N), each type of test being carried out
N
J
time so that the total number
tests rises with:
NR
N
J
J
J N
=
=
=
1
.
This experimental base could for example consist of tests on axisymmetric test-tubes
notched different radii of notch led to various temperatures. Taking into account nature
random of the properties with rupture of material considered, this base constitutes only one sample.
The more important the number of these samples will be, the more it will be representative of the behavior of
material considered.
background image
Code_Aster
®
Version
5.0
Titrate:
Identification of the model of Weibull
Date:
30/11/00
Author (S):
R. MASSON, W. LEFEVRE
Key:
R7.02.09-A
Page:
4/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/00/006/A
Among the various methods of identification suggested in the literature (see for example [bib2]),
we retain two of them: method of regression linear, often used, like that of
maximum of probability recommended by the “Structural European Integrity Society (ESIS)” [bib3].
Note:
A comparative systematic study of the results given by these two methods [bib2] in
function of the number of sample taken by chance on a theoretical distribution showed
that method of the maximum of probability led to a better estimate of
parameters of the model of WEIBULL. Method of regression linear remaining
nevertheless very much used, we integrated it into our developments.
In the two adopted methods of retiming, one carries out the first calculation of the stresses of
WEIBULL with a play of parameter given (typically, m=20,
U
=3000 MPa). These NR is classified
tests using their stress of WEIBULL reached at the instant of the failure. One thus lays out
of an increasing list of stresses of WEIBULL (
W
1
,…,
W
I
,…,
wN
), such as for each
()
I
, it
a many test-tubes broken with a stress of lower WEIBULL or equalize with
W
I
is
N
W
I
(in
General
N
I
W
I
=
). Among the various possible estimators of the probability of cumulated rupture
P
fi
agent with
W
I
[bib2], we choose that generally recommended:
P
I
NR
fi
=
+
1
.
Note:
In the particular case where the stress of WEIBULL depends on the temperature, it
preceding classification must be made temperature by temperature, each temperature
agent with a different statistical law. The estimator of the probability of rupture
precedent thus becomes:
P
I
NR
fi
T
=
+
1
, if the test-tube
()
I
was broken with
temperature T, for which there was
NR
T
tests.
The two adopted methods of retiming are valid as long as [éq 2.1-1] remains true. If identification
is carried out on test results anisothermes whereas the stress of cleavage is supposed
to depend on the temperature, this condition is not checked any more (cf.
POST_ELEM
[U4.81.22]). In it
particular case, one will not be able to thus apply the developments which follow.
3
Method of the linear regression
3.1.1 Principle
The variation theory-experiment is measured by the expression:
LogLog
P
LogLog
P
fi
F
W
I
I
(
)
(
(
))
1
1
1
1
2
-
-
-


éq 3.1.1-1
(“Log” indicates the Napierian logarithm). One wants to minimize this variation compared to (m,
U
).
background image
Code_Aster
®
Version
5.0
Titrate:
Identification of the model of Weibull
Date:
30/11/00
Author (S):
R. MASSON, W. LEFEVRE
Key:
R7.02.09-A
Page:
5/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/00/006/A
3.1.2 Resolution
The method of retiming usually used is based on successive linear regressions: with
the iteration (K), values (
m
K
the U.K.
,
()
) of the module and the stress of cleavage are known. It is thus
possible, with these values, to calculate the stresses of WEIBULL
W K
I ()
at the various moments of
rupture thanks to [éq 2.1-1]. One then classifies these new stresses of WEIBULL per amplitude
increasing and one deduces the new estimates from them from the probability of rupture
P
F K
I ()
with the iteration (K).
For these values of stresses of fixed WEIBULL, the minimization of [éq 3.1.1-1] is brought back to one
simple linear regression on the scatter plot (
Log
W K
I
(
)
()
,
LogLog
P
F K
I
(
)
()
1
1
-
) since if one
defer
LogLog
P
F
(
)
1
1
-
according to
Log
W
(
)
, one obtains a line of slope m which cuts the axis
X-coordinates in (
Log
U
()
). New values (
m
K
the U.K.
+
+
1
1
,
(
)
) of these parameters are thus
data by (cancellation of the derivative partial of [éq 3.1.1-1] compared to each parameter):
m
NR
X
Y
Y
X
NR
X
X
X
K
I K
I J
J K
I K
I K
I
I K
J K
I J
I K
I
+
=
-
-
1
2
1
1
()
,
()
()
()
()
()
,
()
éq 3.1.2-1
the U.K.
I K
I
I K
I
NR
X
m
Y
(
)
()
()
exp
+
=
-










1
1
1
,
éq 3.1.2-2
with
X
Log
I K
W K
I
()
()
(
)
=
and
Y
LogLog
P
I K
F K
I
()
()
(
)
=
-
1
1
.
One repeats these iterations as long as the difference between the plays of parameter obtained with the iterations (K) and
(k+1) is significant (typically, five iterations). The measurement of this variation is given
by:
Max m
m
m
K
K
K
the U.K.
the U.K.
the U.K.
+
+
-
-




1
1
,
(
)
()
()
.
Note:
If
m
is fixed,
the U.K.
(
)
+
1
is always given by [éq 3.1.2-2]. On the other hand, if
U
is fixed,
m
K
+
1
is not given any more by [eq 3.1.2-1] but:
m
X
Y
X
X
K
I K
I
I K
I K
I
U
I K
I
+
=
-
1
2
() ()
()
()
log (
)
.
background image
Code_Aster
®
Version
5.0
Titrate:
Identification of the model of Weibull
Date:
30/11/00
Author (S):
R. MASSON, W. LEFEVRE
Key:
R7.02.09-A
Page:
6/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/00/006/A
4
Method of the maximum of probability
4.1 Principle
Let us note
p
F
W
(
)
density of probability associated with the probability of cumulated rupture
P
F
W
(
)
:
p
m
F
W
U
W
U
m
W
U
m
(
)
exp
=




-


-
1
Quantity
p
D
F
W
W
(
)
is equal to the probability of breaking a test-tube subjected to one
stress corresponding to a stress of WEIBULL included/understood in the interval
[
]
W
W
W
D
,
+
.
The probability so that all the test-tubes of the base broke thus raises with:
p m
D
p
D
U
W
F
W
I
I
W
(,
)
(
)
=
,
éq 4.1-1
p being related to probability. The method of the maximum of probability consists then with
to choose the parameters of the model of WEIBULL so that the function of definite probability
by [éq 4.1-1] (in practice rather its Napierian logarithm) that is to say maximum.
4.2 Resolution
An iterative process again is used. There still, with the iteration (K), (
m
K
the U.K.
,
()
) as well as
W K
I ()
are known. For these values of stresses of fixed WEIBULL, the maximization of
Log p
()
conduit
with a new couple (
m
K
the U.K.
+
+
1
1
,
(
)
) given by:
F m
NR
m
Log
NR
Log
K
K
W K
I
W K
I
m
W K
I
I
I NR
W K
I
m
I
I NR
I
I NR
K
K
(
)
(
)
(
)
(
)
(
)
()
()
()
()
+
+
=
=
=
=
=
=
=
+
-
=
+
+
1
1
1
1
1
1
1
0
éq 4.2-1
(
)
()
(
)
K
W K
I
m
I
I NR
m
NR
K
K
+
=
=
=
+
+
1
1
1
1
1
.
éq 4.2-2
With each pitch, the resolution of [éq 4.2-1] can be carried out using the method of Newton, it
gradient of
F m
()
being given by:
df
DM m
NR m
Log
Log
W
I
m
W
I
I
I NR
W
I
m
I
I NR
W
I
m
W
I
I
I NR
W
I
m
I
I NR
()
(
)
(
)
(
)
(
)
(
)
(
)
= -
+






-














=
=
=
=
=
=
=
=
1
2
2
1
1
1
2
1
2
.
background image
Code_Aster
®
Version
5.0
Titrate:
Identification of the model of Weibull
Date:
30/11/00
Author (S):
R. MASSON, W. LEFEVRE
Key:
R7.02.09-A
Page:
7/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/00/006/A
Note:
If
m
is fixed,
the U.K.
(
)
+
1
is given by [4.2-2]. On the other hand, if
U
is fixed,
m
K
+
1
is not any more
solution of [4.2-1] but of:
F m
NR
m
Log
K
K
W K
I
U
W K
I
U
m
I
I NR
K
(
)
(
) (
(
)
)
()
()
+
+
=
=
=
+
-
=
+
1
1
1
1
0
1
.
This equation can be again solved using the method of Newton, the gradient
being now given by:
df
DM m
NR
m
Log
W
I
U
m
W
I
U
I
I NR
()
(
)
(
)
= -
-
=
=
2
2
1
.
5
Dependence of the parameters with the temperature
If one wishes to fix independently the two parameters temperature by temperature, it is enough to
to break up the base of tests into as much of under - bases by temperature and to apply to each one of
these subbases preceding methods. If, on the other hand, one only wishes to vary
stress of cleavage
U
with the temperature, one proceeds in the following way.
5.1 Regression
linear
The estimate of the probabilities of rupture being now carried out temperature by temperature
(cf notices [§2.2]), it is enough to fix the stress of cleavage on each associated scatter plot
at the various temperatures (T). The equation [éq 3.1.2-2] thus becomes:
the U.K.
T
I K
I T
I K
I T
NR
X
m
Y
(
)
()
()
exp
+
=
-










1
1
1
(
NR
T
indicating the number of tests for the subbase corresponding to the temperature (T)), the module
WEIBULL being given by:
m
NR
X
Y
Y
X
NR
X
X
X
K
T
I K
J K
I T J T
T
I K
I K
I
T
I K
J K
I T J T
T
I K
I
+
=




-




-
1
2
1
1
()
()
,
()
()
()
()
,
()
.
background image
Code_Aster
®
Version
5.0
Titrate:
Identification of the model of Weibull
Date:
30/11/00
Author (S):
R. MASSON, W. LEFEVRE
Key:
R7.02.09-A
Page:
8/8
Manual of Reference
R7.02 booklet: Breaking process
HI-75/00/006/A
5.2
Maximum of probability
The stress of cleavage is given for each temperature (T) considered by:
the U.K.
T
W K
I
m
I T
m
T
NR
T
K
K
(
)
()
()
(
())
+
=
+
+
1
1
1
1
,
m
K
+
1
being solution of:
F m
NR
m
Log
NR
Log
K
K
W K
I
T
W K
I
m
W K
I
I T
W K
I
m
I T
T
I
I NR
K
K
(
)
(
)
(
)
(
)
(
)
()
()
()
()
+
+
=
=
=
+
-
=
+
+
1
1
1
1
1
0
.
6 Conclusion
The control
RECA_WEIBULL
of Code_Aster allows to carry out the chock of the parameters of the model
WEIBULL [U4.82.06].
The user gives in input of this control the concepts results associated with various calculations
nonlinear carried out. The possible dependence of the stress of cleavage with the temperature is
implicitly specified when different temperatures are associated each one of these concepts
results (if all these temperatures are identical or if they are not specified, there is not
dependence with the temperature of this parameter).
The user can carry out this retiming by the method of the maximum of probability (
METHOD:
“MAXI_VRAI”
) or that of the linear regression (
METHOD: “REGR_LIN”
).
Sizes determined by the control
RECA_WEIBULL
are deferred in a table in
which one finds the value of the identified parameters, the probabilities of rupture estimated from
experimental results as well as the probabilities of theoretical rupture calculated with the parameters
identified.
7 Bibliography
[1]
F. BEREMIN, “A local criterion for cleavage fracture off has nuclear presses vessel steel”,
Metall. Trans. 14A, p 2277-2287, 1981.
[2]
A. KHALILI, K. KROMP, “Statistical properties off weibull estimators”, Newspaper off Material
Science, 26, p 6741-6752, 1991.
[3]
ESIS, TC 1.1 one “Local Approach”, Procedure to local measure and calculate approach
criteria using notched tensile specimens “, P6, 1998.