On the unicity of Weibull's distribution function
Volumen Actual  Ciencia Abierta
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.
 

  1. Weibull, W., 1939, "A Statistical Theory of the Strength of Materials", Ing. Vetenskaps Akad. Handl., Vol. 151, pp. 1-45.
  2. 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.
  3. Kittl, P., Baeza, J. and Crisóstomo, J., 1984, "Fracture Statistics of a Brittle Welding", Res Mech., Vol. 12, pp. 275-282.
  4. Kittl, P. and Aldunate, R., 1983, "Compression Fracture Statistics of Compacted Cement Cylinders", J. Mater. Sci., Vol. 18, pp. 2947-2950.
  5. Díaz, G. and Morales, M., 1988, "Fracture Statistics of Torsion in Glass Cylinders", J. Mater. Sci., Vol. 23, pp. 2444-2448.
  6. 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.
  7. 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.
  8. 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.
  9. 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.
  10. 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.
  11. Kittl, P., 1984, "Analysis of the Weibull Distribution Function", J. Appl. Mech., Vol. 51, pp. 221-222.

  12. 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.
     
índice

página anterior

página siguiente

Ciencia Abierta

 

Ciencia Abierta



Posicionamiento Web Edreams

Web Counter