On the unicity of Weibull's distribution function.
G. Díaz (1), V. Martínez (2) and
P. Kittl (2)
(1) Departamento de Ingeniería de los
Materiales, IDIEM,
Facultad de Ciencias Físicas y Matemáticas,
Universidad de Chile, Casilla 1420, Santiago, Chile.
(2) Departamento de Ingeniería Mecánica,
Facultad de Ciencias Físicas y Matemáticas,
Universidad de Chile, Casilla 2777, Santiago, Chile.
|
Abstract
The present work aims to obtain a generalization of the fundamental
functional equation of the Probabilistic Strength of Materials. This was
not possible but this way has the advantage of proving that only one statistics,
the one of Weibull, is compatible with the boundary conditions reached
in the cases of fracture and yield stress in brittle and ductile materials,
respectively.
|
Introduction.
A large number of mechanical properties whose knowledge is required
for designing structures are non-deterministic by nature and they follow
a certain distribution function of probability. Fracture stresses in brittle
materials as well as yielding stresses in ductile materials obey Weibull's
distribution function [1,2]. As a matter of fact this last function has
been widely used to explain the probabilistic distribution of said stresses
in a large collection of materials and in diverse strain-states[3-8]. The
fundamentals of the Probabilistic Strength of Materials as well as the
sundry deductions of Weibull's distribution function have been discussed
in many works [2,9-12]. In a previous work, Kittl and Dìaz [12]
have discussed five different deductions of Weibull's statistics. It is
interesting to find a more general distribution function in order to see
if anomalous behavior in fracture materials can be explained through this
manner.
Generalization of the functional equation.
Consider an isotropic homogeneous body of volume V subjected to a constant
state of stress, and an arbitrary division of V into volumes V1
and V2 without common elements so that V=V1 + V2.
Let
N be the number of similarly-fabricated samples of a material and that
are to be subjected to a fracture or yielding test when the material is
brittle or ductile, respectively. Three possibilities may arise in this
connection, namely: volume V1 fails; volume V2 fails;
and there is no failure. Let n1 be the number of fractures started
within volume V1, n2 the number of fractures started
within volume V2, and n the number of tests without any fracture
at all, so that then N = n1 + n2 + n. All these tests
are made by applying, gradually, some stress s
so that 0 £s£si where
i= 1,2,3,..., N.
According to the Probabilistic Strength of Materials, the fundamental
functional equation originating the same is [2]:
|
|
(1)
|
where is the cumulative
probability of non-fracture of the body of volume V while
and are the cumulative probabilities
of non-fracture of the bodies of respective volumes V1 and V2
when the material is subjected to a constant stress s
. The functional equation (1) is based on the hypothesis of the independence
of the cumulative probabilities of non-fracture. Hence, considering the
mentioned experience with N tests, the following probabilities derive therefrom:
|
|
(2)
|
where F1 and F2 are the cumulative probabilities
of fracture of the bodies whose respective volumes are V1 and
V2. Using equation (2) allows to rewrite equation (1):
|
|
(3)
|
and here, to simplify writing, the argument of the functions of cumulative
probability of non-fracture have been omitted.
Equation (3) can be generalized as follows:
|
|
(4)
|
where F12 is the cumulative probability of fracture of body
of volume V and l is a real number. In equation
(4), to the set of cases without fracture there must be added the set of
cases with fracture in volume V1 only (none in V2),
and the set of cases with fracture in volume V2 only (none in
V1).
Using again equation (2) allows to rewrite equation (4):
|
|
(5)
|
and this expression becomes in terms of the functions of cumulative
probabilities of fracture and of non-fracture:
|
|
(6)
|
Obviously functional equation (6) yields functional equation (1) in
the particular instance of l =1.
Functional equation (6) can be converted into:
|
|
(7)
|
Unicity of cumulative probabilty of non-fracture.
Now it will be shown that the function of cumulative probability of
non-fracture is unique, using therefor the following boundary condition:
|
|
(8)
|
that is, the non-fracture probability will be approaching the value
1 in the case of some body whose very small volume V ®
0 is devoid of whatever defect liable to start some fracture. According
to equation (7), if V1=0 and V2=V:
|
|
(9)
|
and likewise if V2=0 and V1=V:
|
|
(10)
|
Consequently the function of cumulative probability of non-fracture
is unique, that is to say:
|
|
(11)
|
Hence the functional equation (7) can be written as follows:
|
|
(12)
|
Solution of the generalized functional equation.
Now functional equation (12) and, therefore, functional equation (7)
will be solved. Small volume-increments yield:
|
|
(13)
|
Considering the following Taylor-series expansion:
|
|
(14)
|
its introduction into equation (13) produces:
|
|
(15)
|
Similarly, also can be
expanded as a Taylor series:
|
|
(16)
|
Equating equations (15) and (16) and making some rearrangements yield
the following differential equation allowing to solve the functional equation
(7):
|
|
(17)
|
Putting:
|
|
(18)
|
differential equation (17) becomes:
|
|
(19)
|
whose solution is:
|
|
(20)
|
where k is a constant determined using the boundary condition .
Thus the solution to differential equation (19) is:
|
|
(21)
|
After carrying out some transformations the solution of functional equation
(7) is:
|
|
(22)
|
The expression for is
already known and can be found in [2], for example:
|
|
(23)
|
where V0 is the unit volume whereas f
(s ) is the Weibull’s specific risk function.
Hence equation (22) becomes:
|
|
(24)
|
Equation (24) complies with boundary condition .
In addition, there is also another boundary condition in the ambit of Probabilistic
Strength of Materials, namely:
|
|
(25)
|
that is, non-fracture probability will be approaching the value zero
in the case of some body whose very large volume V ®¥
always includes some defect liable to start some fracture. When V ®¥
, equation (24) is reduced to:
|
|
(26)
|
whose value is not zero except when l = 1.
Accordingly, the probability function obtained by solving functional equation
(7) is only compatible with boundary conditions set forth in equations
(8) and (25) when l = 1, so that instance l
= 0 remains to be studied below.
Extreme instances l = 0 and l
= 1 regarding functional equation (7) define two distribution functions.
Thus, if l = 1 then the following distribution
function is obtained from equation (24):
|
|
(27)
|
and like result is obtained by putting l
= 1 in equation (17). Using the cumulative probability of fracture function
the equation (27) yields:
|
|
(28)
|
which is the well-known Weibull’s distribution function [1]. On the
other hand, if l =0 then the following differential
equation is obtained from equation (17):
|
|
(29)
|
and its solution is:
|
|
(30)
|
where C1 is a constant determined using the boundary condition
defined in equation (8) and its value is 0. Hence the solution to equation
(29) is:
|
|
(31)
|
Employing the cumulative probability of fracture function and introducing
equation (23) into equation (31) yields:
|
|
(32)
|
The distribution function of cumulative probability defined by equation
(32) is a Weibull statistics when the volume of the material, or the stress
applied thereto, is approaching the value zero, as can be seen when equation
(27) is expanded in Taylor series.
In short, according to the generalization proposed herein for the functional
equation of Probabilistic Strength of Materials, only one statistics compatible
with the boundary conditions is existing, that is, the Weibull’s distribution
function when l = 1, equation (28), because
the distribution function when l =0, equation
(32), can be obtained from the above equation (28) by making (V/V0)f
(s)® 0.
Acknoledgments.
The authors are grateful to the Fondo Nacional de Desarrollo Científico
y Tecnológico, FONDECYT, for the Grant of Projects N° 1931056
and 1961105.
References.
-
Weibull, W., 1939, "A Statistical Theory of the Strength of Materials",
Ing. Vetenskaps Akad. Handl., Vol. 151, pp. 1-45.
-
Kittl, P. and Díaz, G., 1988, "Weibull’s Fracture Statistics, or
Probabilistic Strength of Materials: State of the Art", Res Mech., Vol.
24, pp. 99-207.
-
Kittl, P., Baeza, J. and Crisóstomo, J., 1984, "Fracture Statistics
of a Brittle Welding", Res Mech., Vol. 12, pp. 275-282.
-
Kittl, P. and Aldunate, R., 1983, "Compression Fracture Statistics of Compacted
Cement Cylinders", J. Mater. Sci., Vol. 18, pp. 2947-2950.
-
Díaz, G. and Morales, M., 1988, "Fracture Statistics of Torsion
in Glass Cylinders", J. Mater. Sci., Vol. 23, pp. 2444-2448.
-
Martínez, V., Kittl, P. and Díaz, G., 1992, "Fracture Statistics
of Torsion and Flexure in Glass Rectangular Bars", J. Mater. Sci., Vol.
27, pp. 1457-1463.
-
Kittl, P., Martínez, V., Díaz, G., Bölcich, J.C. and
Fernández, L., 1993, "Non-linear Deformation Probability of an ASTM-516
Steel", J. Mater. Sci. Lett., Vol. 12, pp. 823-824.
-
Weil, N. A. and Daniel, I. M., "Analysis of the Fracture in Nonuniformly
Stressed Brittle Materials", J. Am. Ceram. Soc., Vol. 47, pp. 268-274.
-
Freudenthal, A. M., 1968, Statistical Approach to Brittle Fracture, In:
Fracture, An Advanced Treatise, H. Liebowitz (ed), Academic Press, New
York, Vol. 2, Chap. 6, pp. 591-619.
-
Kittl, P., Díaz, G. and León, M., 1986, "Theoretical and
Experimental Investigations on Fracture Statistics Carried Out at the IDIEM
(Chile)", Materiales de Construcción, Vol. 36, pp. 3-19.
-
Kittl, P., 1984, "Analysis of the Weibull Distribution Function", J. Appl.
Mech., Vol. 51, pp. 221-222.
Kittl, P. and Díaz, G., 1990, "Five Deductions of Weibull’s
Distribution Function in the Probabilistic Strength of Materials", Engng.Fract.Mech.,
Vol. 36, pp. 749-762.
|
|