MODELADO DE LA CONFIABILIDAD

ALEXÁNDER REÁTIGA VILLAMIZAR1

1 Departamento de Ciencias Básicas

Universidad Pontificia Bolivariana - Bucaramanga

2023 I semestre
Viernes-Sábado
Agosto 25 y 26 / 23

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 1 / 65
 LA DISTRIBUCIÓN DE WEIBUL

La distribución de Weibull de dos parámetros es ampliamente la

distribución más usada para el análisis de datos de vida.
Se dice que una variable aleatoria continua X tiene una
DISTRIBUCIÓN DE WEIBULL con parámetros β y η, mayores que
cero, y se escribe X ∼ Weibull(β, η) si su función de densidad de
probabilidad PDF esta dada por:
 β
t β−1 − ηt
f (t) = β β e .
η
Donde:
t es el tiempo como variable.
η es la vida caracterı́stica o el parámetro de escala.
β es el parámetro de forma.
MAESTRÍA EN CONFIABILIDAD Y MANT
(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 2 / 65
 ALGUNAS DISTRIBUCIONES DENTRO DE
WEIBULL

Algunas distribuciones conocidas, son casos particulares de la distribución

de Weibull.
1 Para η = 1 y β = 3, 44 la PDF de Weibull es una buena aproximación
a la distribución normal.
1 (x−µ)2
f (x) = √ e − 2σ2 .
σ 2π
Conocida como la campana de Gauss.
2 La función de distribución de probabilidad CDF de la distribución
normal es:  x − µ i
1h
F (x) = 1 + erf √ .
2 σ 2

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 3 / 65
 ALGUNAS DISTRIBUCIONES DENTRO DE
WEIBULL

1 Para η = 1 y β = 1, se obtiene la distribución exponencial con

parámetro λ = 1.

f (x) = λ e −λx.
2 La función de distribución de probablidad de la exponencial es:

F (x) = 1 − e −λx .

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 4 / 65
 FUNCIONES EQUIVALENTES

Otra forma de escribir la función de densidad es:

 α
 α−1 − λx
f (x) = λα λx e α > 0, λ > 0.

Esta equivalencia se logra haciendo:

1
α=β y λ= .
η

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 5 / 65
 DISTRIBUCIÒN DE WEIBULL Y SU FUNCIÒN
DE DISTRIBUCIÒN DE PROBABILIDAD
ACUMULADA

Si f (t) es la función de distribución de probabilidad (función de densidad

de probabilidad) de una variable aleatoria continua X , entonces su
FUNCIÓN DE DISTRIBUCIÓN DE PROBABILIDAD ACUMULATIVA
CDF es:
Z t
F (t) = P(X ≤ t) = f (x) dx.
0
Y Z b
P(a ≤ X ≤ b) = f (x) dx.
a

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 6 / 65
 DISTRIBUCIÒN DE WEIBULL Y SU FUNCIÒN
DE DISTRIBUCIÒN DE PROBABILIDAD
ACUMULADA
Para el caso de la distribución de Weibull, se tiene que:
 x β
t
x β−1 −
  Z
η


F t =P X ≤t = β e dx
0 ηβ
 t β
 x β − .
η


Haciendo u = , se obtiene que F t = 1 − e
η
En la forma alternativa:
 α

  Z x  α−1 − λy
F x =P X ≤x = λα λy e dy
0
 α
Se logra el mismo resultado, haciendo u = λy .
MAESTRÍA EN CONFIABILIDAD Y MANT
(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 7 / 65
 DISTRIBUCIÒN DE WEIBULL Y SU FUNCIÒN
DE DISTRIBUCIÒN DE PROBABILIDAD
ACUMULADA

En los anteriores resultados, tomando el lı́mite cunado t −→ ∞ se muestra

que F (t) es una función de densidad de probabilidad .
Ya que una función de distribución de probabilidad debe cumplir:
1 f (x) ≥ 0, para toda x.
2 Z ∞
f (x) dx = 1.
−∞
3
Z b
P(a ≤ x ≤ b) = f (x) dx.
a

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 8 / 65
 FUNCIÓN GAMMA - TRANSFORMADA DE
LAPLACE
La transformada de Laplace de una función f (t), se define como:
Z ∞
L{f (t)} = e −st f (t) dt, s > 0.
0
La función Gamma se define como:
Z ∞
Γ(α) = t α−1 e −t dt, α > 0.
0

Se puede verificar que:

1

Γ(α + 1) = αΓ(α).
2
Γ(α + 1)
L{t α } = .
s α+1
MAESTRÍA EN CONFIABILIDAD Y MANT
(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 9 / 65
 ESPERANZA Y VARIANZA

Se puede verificar que la esperanza matemática y la varianza de una

variable aleatoria X con distribución de Weibull, son:
1
1 h 1i h 1i
E (X ) = Γ 1+ = ηΓ 1+ .
λ α β
2
1h  2 2
 1 i
Var (X ) = Γ 1 + − Γ 1 + .
λ2 α α

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 10 / 65
 ESPERANZA Y VARIANZA

Z ∞
E (X ) = xf (x) dx.
−∞
En este caso obtendrı́amos:
Z ∞
α
E (X ) = x λα (λx)α−1 e −(λx) . dx
0
Z ∞
α
= α (λx)α e −(λx) . dx
0
1
Haciendo la sustitución: u = (λx)α se obtiene que: 1
λα u α −1 du = dx.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 11 / 65
 ESPERANZA Y VARIANZA

La integral se convierte en:

Z ∞
1
−u 1 −1
αue u α du.
0 λα
Z ∞
1 1
−u 1  1
= u α e du = Γ 1 + .
λ 0 λ α

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 12 / 65
 VIDA MEDIA µ = MTTF

El tiempo medio para fallos se puede calcular a partir de función de

densidad de probabilidad.
Z ∞
MTTF = xf (x) dt.
0
Luego:
1  1
MTTF = Γ 1+ = E (X ).
λ α

 1
MTTF = η Γ 1 + = E (X ).
β

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 13 / 65
 ESPERANZA Y VARIANZA

La varianza de la variable aleatoria continua X , se puede conseguir como:

VAR(x) = σ 2 = E (X 2 ) − (EX )2 ,

donde Z
2
E (X ) = x 2 f (x) dx,
R
a esto se le conoce como el segundo momento de la variable aleatoria
continua.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 14 / 65
 ESPERANZA Y VARIANZA

Se puede verificar que:

1  2 1 2 1
VAR(X ) = Γ 1 + − Γ 1 + .
λ2 α λ2 α
Calculando el segundo momento se tiene:
Z Z ∞  α−1 α
2
E (X ) = 2
x f (x) dx = x 2 λα λx e −(λx) dx.
R 0
 α  α−1
Haciendo u = λx , se obtiene: du = λα λx dx y
1 2
u α = x 2 . Y ası́ se llega al resultado requerido.
λ2

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 15 / 65
 DATOS DE ENTRADA EN EL A-W

En el A-W se usan parámetros continuos como:

1 Tiempos hasta el fallo.
2 Ciclos hasta el fallo.
3 Esfuerzo mecánico.
Se trabaja con una muestra aleatoria de elementos (equipos)
operando bajo condiciones de ensayo o en servicio, para estimar
comportamiento de fiabilidad en la población de la cual se extrajo la
muestra.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 16 / 65
 DATOS DE ENTRADA EN EL A-W

Los datos analizados deben ser independientes e idénticamente

distribuidos. Se abordan métodos numéricos y gráficos para representar los
datos, para hacer la bondad del ajuste, para estimar los parámetros y para
mostrar los lı́mites de confianza.
Se interpretan gráficos en términos del riesgo como función del tiempo,
modos de fallo, poblaciones débiles, tiempo hasta el fallo o duración
mı́nima.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 17 / 65
 TÉRMINOS Y DEFINICIONES

TRUNCAMIENTO O CENSURA: Terminación de un ensayo después

de una duración dada o de un número dado de fallos.Un ensayo terminado
con elementos que no han fallado, se llama ensayo truncado, y los datos de
estos tiempos de ensayo usados se llaman datos censurados.

ELEMENTO SUSPENDIDO: Elemento sobre el cual se ha acortado el

ensayo sin encontrar fallo a considerar. El elemento puede no haber fallado
o fallar de modo distinto al estudiado.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 18 / 65
 TÉRMINOS Y DEFINICIONES

ENSAYO DE VIDA: Ensayo realizado para estimar la durabilidad de un

producto.

ELEMENTO NO REPARABLE: Elemento que después de un fallo, no

puede volver a desempeñar las funciones requeridas.

TIEMPO DE FUNCIONAMIENTO: Intervalo de tiempo en el cual el

elemento está en estado de funcionamiento.

FALLO A CONSIDERAR: Fallo que deberı́a incluirse a la hora de

interpretar los resultados de los ensayos.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 19 / 65
 TÉRMINOS Y DEFINICIONES

ENSAYO DE FIABILIDAD: experimento realizado para medir,

cuantificar o clasificar una medida o caracterı́stica de fiabilidad de un
elemento.

ELEMENTO REPARABLE: Elemento que después del fallo, puede

volver a desempeñar las funciones requeridas.

TIEMPO DE FUNCIONAMIENTO HASTA EL FALLO: Tiempos de

funcionamiento acumulados desde su primera utilización o desde su
restauración hasta la aparición de un fallo.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 20 / 65
 TÉRMINOS Y DEFINICIONES

TIEMPO ENTRE FALLOS: Tiempo transcurrido entre dos fallos

consecutivos.

VIDA B, PERCENTILES L: Tiempo en el cual una cantidad dada de

elementos ha fallado.
B10 Es el tiempo hasta el cual el 10 % de los elementos han fallado. B50
Se le llama tiempo mediano hasta el fallo. Tiempo en el cual ha fallado el
50 % de los elementos.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 21 / 65
 LA DISTRIBUCIÓN DE WEIBUL

Se dice que una variable aleatoria continua X tiene una DISTRIBUCIÓN

DE WEIBULL con parámetros β y η, mayores que cero, y se escribe
X ∼ Weibull(β, η) si su función de densidad esta dada por:
 β
t β−1 − t
η
f (t) = β β e .
η
Donde:
t es el tiempo como variable.
η es la vida caracterı́stica o el parámetro de escala.
β es el parámetro de forma.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 22 / 65
 ALGUNAS CURVAS DE LA FAMILIA

Algunas curvas de esta familia de funciones para diferentes valores del

parámetro β.

El valor de β = 3, 44 y η = 1, se parece a la
distribución normal.
MAESTRÍA EN CONFIABILIDAD Y MANT
(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 23 / 65
 ALGUNAS CURVAS DE LA FAMILIA

Algunas curvas de esta familia de funciones para diferentes valores del

parámetro η.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 24 / 65
 FUNCIÓN DE DISTRIBUCIÓN ACUMULATIVA
DE WEIBULL

La función de distribución acumulativa de Weibull es:

β
− t
F (t) = 1 − e η .

El parámetro de forma β, indica la variación de la tasa de fallo con el

tiempo.
Se incluyen ejemplos de:
1 Mortalidad infantil.
2 Fallos aleatorios.
3 Fallos por desgaste.
La variable t es genérica y puede representar tiempo, distancia, número de
ciclos o aplicaciones de esfuerzo mecánico.
η es el tiempo en que se espera que falle el 63, 2 %
MAESTRÍA EN CONFIABILIDAD Y MANT
(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 25 / 65
 FUNCIÓN DE DISTRIBUCIÓN ACUMULATIVA
DE WEIBULL

La función de distribución acumulativa de Weibull es:

β
− ηt
F (t) = 1 − e .

Tambien llamada la LA DESCONFIABILIDAD.

La función
β
t
−
R(t) = 1 − F (t) = e η .
Es llamada LA CONFIABILIDAD.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 26 / 65
 TASA DE FALLOS - FUNCIÓN DE RIESGO

La tasa de fallos λ(t) o función de riesgo h(t), de la distribución de

Weibull es:
f (t) f (t) t β−1
λ(t) = h(t) = = =β β .
1 − F (t) R(t) η

1 β < 1 indica una tasa de fallos decreciente.

2 β = 1,la distribución de Weibull se convierte en la distribución

exponencial y la tasa de fallos es constante, igual al recı́proco de η.

3 β > 1, indica una tasa de fallos creciente.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 27 / 65
 FUNCIÓN DE RIESGO ACUMULADA

La función de riesgo acumulada:

Z t Z t
f (s)
H(t) = h(s) ds = ds =
0 0 1 − F (s)
t  t β
-ln[1 − F (s)] = − ln[1 − F (t)] = .
0 η

Otra alternativa:
t t
s β−1
Z Z  t β
H(t) = λ(s) ds = β ds = .
0 0 ηβ η

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 28 / 65
 DISTRIBUCIÓN DE WEIBULL DE TRES
PARÁSMETROS

La distribución de Weibull de 3 parámetros se define como:

(t − t0 )β−1 −( t−t 0 )β
f (t) = β β
e η .
η
La función de confiabilidad R(t) de Weibull se determina por la siguiente
expresión:
Z ∞ Z ∞
(s − t0 )β−1 −( s−t0 )β
R(t) = f (s) ds = β β
e η ds
t t η

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 29 / 65
 DISTRIBUCIÓN DE WEIBULL DE TRES
PARÁMETROS

De donde se obtiene que:

t−t0 β
−( )
R(t) = e η .

Por tanto la función de distribución acumulativa (de probabilidad para el

fallo) es:
t−t0 β
−( )
F (t) = 1 − e η .

El parámetro t0 se conoce como: tiempo libre de fallos, Parámetro de

localización o vida mı́nima. Este parámetro genera un corrimiento de la
gráfica hacı́a la derecha, de tal modo que la gráfica comienza en t0 ,
tiempo en que ocurre el primer fallo.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 30 / 65
 CONSIDERACIONES SOBRE LOS DATOS -
TIPOS DE DATOS

Los datos de vida están asociados a elementos que ”envejecen” hasta el

fallo.

Los datos de Weinbull, también pueden describir datos del material cuando
el envejecimiento se debe a esfuerzo, fuerza o temperatura.

La ”edad” puede referirse a: tiempo de operación, paradas y arranques,

aterrizajes, despegues, ciclos, kilómetros y muchos otros parámetros
continuos.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 31 / 65
 CONSIDERACIONES DE LOS DATOS - TIEMPO
HASTA EL PRIMER FALLO

La variable ”tiempo” en Weibull, habitualmente se interpreta como una

medida de la vida consumida.

OTRAS CONSIDERACIONES:
1 Tiempo hasta el primer fallo de un elemento reparable.
2 Tiempo hasta el fallo de un elemento no reparable.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 32 / 65
 TAMAÑO DE LA MUESTRA

La incertidumbre relativa a la estimación del parámetro de Weibull, está

relacionada con el tamaño de la muestra.

Puede estimarse con dos fallos, pero se generarı́a una incertidumbre

excesiva.

Sin importar el tamaño de la muestra, deben calcularse y representarse los

lı́mites de confianza para evaluar la incertidumbre de las estimaciones.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 33 / 65
 DATOS CENSURADOS O TRUNCADOS Y
SUSPENDIDOS

Es importante incluir datos de vida de los elementos que no han fallado o

lo han hecho de forma diferente al fallo en estudio (datos censurados o
truncados o suspendidos).

Los datos se llaman completos cuando se han observado los tiempos hasta
el fallo de todos los elementos.

Un elemento ensayado, que no ha fallado o haya fallado de forma diferente

al estudiado, se llama un elemento censurado o una suspensión; las hay
tempranas (suspensiones antes del primer fallo) y tardı́as (se suspende el
ensayo después del último fallo) y entre fallos, se llaman aleatorias o
progresivas.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 34 / 65
 MÉTODOS GRÁFICOS Y BONDAD DE AJUSTE

Consiste en representar los datos en el papel probabilı́stico de Weibull.Se

ajusta una lı́nea a los datos se interpreta el gráfico y se estiman los
parámetros. Esto se logra transformando la ecuación de Weibull en una
expresión lineal.

Primero se hace una ordenación. Los tiempos hasta el fallo, se representan

en el eje de coordenadas X. En el eje Y se representan los rangos medianos.

Para muestras mayores de 30, el rango mediano es lo mismo que el

porcentaje de fallos.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 35 / 65
 MÉTODOS GRÁFICOS Y BONDAD DE AJUSTE

Para muestras pequeñas el rango mediano se calcula de la fómuna de

Benard:
i − 0,3
Fi = %.
N + 0,4
Donde N es el tamaño de la muestra e i es la posición de representación
de los datos.

Para N > 30, se puede utilizar:

i
Fi = × 100 %.
N

Si los datos representados siguen una tendencia lineal, se puede trazar una
recta de regresión.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 36 / 65
 MÉTODOS GRÁFICOS Y BONDAD DE AJUSTE

De la lectura del gráfico, se pueden obtener los parámetros. β es la

pendiente de la recta de regresión.

η es tiempo en el cual el 63, 2 % de los elementos ha fallado, es llamada

B63,2 .

Esta técnica es llamada la regresión de los rangos medianos (MRR).

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 37 / 65
 CÓMO CONSTRUIR UN GRÁFICO DE
PROBABILIDAD

Se deben seguir los siguientes pasos: ORDENACIÓN.

1 Se deben ordenar los datos de tiempo hasta el fallo de menor a
mayor.

2 Las representaciones para los tiempos t, se ponen en las abscisas.

3 F (t) en valores porcentuales, se acomodan en las ordenadas.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 38 / 65
 GRÁFICO DE PROBABILIDAD DE WEIBULL

Después de transformar los datos, se puede construir el gráfico utilizando

tres métodos diferentes.
1 Papel probabilı́stico de Weibull.

2 Programa informático de hojas de cálculo.

3 Sofware comercial.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 39 / 65
 EJEMPLO ECUADOR

Ejemplo norma - INCERTIDUMBRE.

Analizar estudio de caso en Ecuador:

ESTUDIO DE FIABILIDAD, MANTENIBILIDAD Y
DISPONIBILIDAD APLICADO A GRUPOS ELECTRÓGENOS
PRIME.

AUTORES:
César Marcelo Gallegos Londoño.
Mayra Alexándra Viscaı́no Cuzco.
Sergio Raúl Villacrés Parra.

Ecuador 2020.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 40 / 65
 TRATAMIENTO PARA SUSPENSIONES O
DATOS CENSURADOS

Los datos sobre los elementos censurados o suspendidos, no se pueden

ignorar. Los tiempos de estos elementos deben incluirse dentro del análisis.

En el tratamiento de estos datos, se utiliza la siguiente fórmula para

ajustar los rangos. Se deben incluir columnas para rango invertido y rango
ajustado. Ver ejemplo.

Rango invertido × Rango ajustado previo + N + 1
Rango ajustado =

Rango invertido + 1

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 41 / 65
 EJEMPLO CON DATOS CENSURADOS O
SUSPENDIDOS

Ver ejemplo en excell.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 42 / 65
 PASOS PARA REPRESENTAR DATOS CON
SUSPENSIONES

Estos son los pasos para representar datos con suspensiones.

1 Ordenar tiempos de fallos y
suspensiones de primero a
último. 1 Representar los tiempos en (x) y
2 Calcular rangos ajustados para los rangos medios en (y).
fallos, las suspensiones no se 2 Estimar η leyendo la vida B63,2 .
representan. 3 Estimar β.
3 Calcular rangos medios con
Benard.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 43 / 65
 COMPROBACIÓN DEL AJUSTE

Si los datos representados se agrupan alrededor de una recta en el gráfico

de probabilidad, es suficiente evidencia que los datos se encuentran
representados por la distribución en cuestión.

Muestras pequeñas dificultad la bondad del ajuste.

En esta norma se utiliza el cuadrado del ı́ndice de correlación. El cual se

calcula de la siguiente manera:

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 44 / 65
 ÍNDICE DE CORRELA-
CIÓN2 =COEFICIENTE DE DETERMINACIÓN

Donde: xi : Tiempos de fallo , yi : Rangos medianos. x̄, ȳ son las

medias y N es el tamaño de la muestra.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 45 / 65
 GRÁFICO DE RIESGO

Se comienza con la estimación de la tasa de fallo o función de riesgo al

fallo.

f (t)
λ(t) = h(t) =  .
1 − F (t)
La función de riesgo acumulada,
Z t
 
H(t) = h(t) dt = − ln 1 − F (t)
0

se estima como la suma acumulada de funciones de riesgo al fallo

estimadas.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 46 / 65
 GRÁFICA DE RIESGO

Para la distribución de Weibull,

 t β
H(t) = .
η

Desde la cual se puede obtener una relación lineal para ln H(t) en función
de ln(t) :

ln[H(t)] = β ln(t) − β ln(η).

Ası́: el papel de riesgo de Weibull, es por tanto papel doble logarı́tmico
(ln − ln). La pendiente de la lı́nea ajustada a los datos es β, y t = η
cuando H(t) = 1.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 47 / 65
 CÓMO UTILIZAR EL PAPEL WEIBULL?

Se deben graficar los valores de F (t) en el eje vertical Y con sus

respectivos tiempos en el eje X .

La lı́nea de regresión, permite encontrar los parámetros β como la

pendiente de la recta y η con el corte en el eje vertical.
Las componentes del papel Weibull son:

1 Número de prueba. 7 β parámetro de forma.

2 Elemento, máquina o 8 η parámetro de escala o vida
componente a evaluar. caracterı́stica.
3 Fecha.
4 Tipo de prueba.
9 γ parámetro de posición.
5 Promedio de tiempo de fallas. 10 Punto de referencia para hallar
6 Tamaño de la muestra. Pµ y β.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 48 / 65
 CÓMO UTILIZAR EL PAPEL DE WEIBULL

1 Se grafican los tiempos vs los porcentajes acumulados de falla.

2 Se traza una lı́nea, tratando que sea lo mejor aproximada a los puntos.
3 Desde el punto de referencia, se traza una perpendicular a la recta de
regresión.
4 Esta lı́nea permite encontrar los valores β y Pµ . La media para la falla
µ.
5 Desde el eje vertical, pasando por Pµ estimado, se traza una
horizontal a la lı́nea de regresión y por este corte una vertical para
encontrar (MTBF).
6 Desde el eje vertical, pasando por ”η estimador”se traza un
perpendicular a la linea de regresión y por este corte, una vertical al
eje del tiempo para encontrar la vida caracterı́stica.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 49 / 65
 CÓMO UTILIZAR EL PAPEL WEIBULL?

Ejemplo de 8 observaciones: El rango medio se obtiene haciendo:

i − 0, 3
MR = .
n + 0, 4

i falla Rango medio F(t)

1 1717 0.08333333 8.33333
2 2082 0.20238095 20.2381
3 2263 0.32142857 32.1429
4 2751 0.44047619 44.0476
5 3065 0.55952382 55.9524
6 3945 0.67857143 67.8571
7 5093 0.79761905 79.7619
8 7062 0.91666667 91.6667

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 50 / 65
 GRÁFICO DE PROBABILIDAD DE WEIBULL
LA CURVA DE LA BAÑERA: Muestra la rfelación entre el parámetro de
forma β de la distribución de Weibull y la función de riesgo a lo largo de la
vida de un elemento.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 51 / 65
 LA CURVA DE LA BAÑERA

MORTALIDAD INFANTIL: β < 1,

1 Dispositivos defectuosos o instalados indebidamente.
2 Malos diseños del fabricante.
3 Materiales no adecuados en la fabricación.
4 Esta tasa de fallos disminuye en el tiempo hasta alcanzar un valor
constante.
FALLOS ALEATORIOS: β = 1.
1 Tasa de fallos constante en el tiempo.
2 Los fallos se deben a condiciones normales de funcionamiento.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 52 / 65
 LA CURVA DE LA BAÑERA

FALLOS DE DESGASTE: β > 1.

1 Cuando el componente supera la vida útil prevista.
2 Fallos por envejecimiento.
3 Fallos de degradación por desgaste.
4 Se evidencia una aceleración en la tasa de fallos.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 53 / 65
 MÉTODOS COMPUTACIONALES Y DE BONDAD
DE AJUSTE

El método de estimación por máxima verosimilitud (EMV) es un método

computacional para el caso de muestras grandes.

Trabaja datos con mecanismos complicados de censura y suspensiones.

Primero se hace el contraste de bondad del ajuste, para comprobar la

suposición de Weibull. Si la hipótesis no es rechazada, se procede a la
EMV.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 54 / 65
 SUPOSICIONES Y CONDICIONES

1 Una muestra de n elementos no reparables provienen de la misma

población.
2 El entorno de ensayo debe ser el mismo para todos los elementos
sometidos.
3 Los elementos fallados no son reemplazados.
4 El ensayo termina en un tiempo T y hay r elementos que han fallado.
5 El tiempo T cumple que tr ≤ T .
6 Hay r tiempos hasta el fallo y cumplen: ti ≤ ti+1 , i = 1, 2, 3, ...r .

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 55 / 65
 CONTRASTE DE BONDAD DEL AJUSTE

1 Se ordenan todos los r tiempos hasta el fallo de forma ascendente.

2 Se toman los logaritmos neperianos de estos tiempos;

ln(t1 ) = x1 , ln(t2 ) = x2 , . . . , ln(tr ) = xr . x1 ≤ x2 ≤ . . . ≤ xr .

3 Calcular las cantidades li para todos los valores de r .

xi+1 − xi
li = h i
ln ln( 4(n−i−1)+3
4n+1 )/ ln( 4(n−i)+3
4n+1 )

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 56 / 65
 CONTRASTE DE BONDAD DEL AJUSTE

Calcular la cantidad H dada por:

r −1
X li
r −1
⌊ 2 ⌋
i=⌊ 2r ⌋+1
H= r .
⌊2⌋
X li
⌊ 2r ⌋
i=1

Rechazar la hipótesis de que los datos proceden de una distribución

 de 
(r −1)
Weibull al nivel de significación de γ del 100 % si H ≤ Fγ 2⌊ 2 ⌋, 2⌊ 2r ⌋
y no continuar con el análisis.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 57 / 65
 ESTIMACIONES PUNTUALES DE β y η

La EMV de β, o β̂ es aquella que satisface la ecuación:

 r 
X β
β
 ti ln(ti ) + (n − 1)T ln(T )  r
 i=1 1  1X

r − − ln(ti ) = 0.
 X β
β  r
ti + (n − r )T β i=1
 
i=1

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 58 / 65
 ESTIMACIONES PUNTUALES DE β y η

La η̂ se encuentra utilizando a β̂ según la ecuación:

( " r #) 1
β̂
1 X β̂ β̂
η̂ = ti + (n − r )T .
r
i=1

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 59 / 65
 ESTIMACIÓN PUNTUAL DEL TIEMPO MEDIO
HASTA EL FALLO

Los valores de β̂ y η̂ se utilizan para la estimación puntual del tiempo

medio hasta el fallo, de acuerdo a la siguiente ecuación:
 1
m̂ = η̂ Γ 1 + .
β̂

Donde Γ(z) es la función gamma de (z), la cual se puede buscar en tablas,

si el valor no aparece en estas, se puede hacer interpolación lineal.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 60 / 65
 ESTIMACIÓN PUNTUAL DEL PERCENTIL (10)

La estimación puntual del B10 se hace mediante la relación:

1
1
β̂
B̂10 = η̂ ln .
0,9


La fiabilidad en el instante t, t ≤ T . se consigue con la relación:

1
t β̂
−
R̂(t) = e η̂ .

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 61 / 65
 EL MÉTODO DE LA MUERTE SÚBITA

1 Se requieren de subgrupos pequeños de elementos en prueba, entre 3

y 8 elementos por subgrupo.

2 Se pone en marcha el ensayo hasta que se produzca el primer fallo, en

este caso se tiene, un fallo y el resto serán suspensiones, en cada
subgrupo.

3 10 subgrupos de 4 elementos, arrojarán 10 fallos y 4 suspensiones.

4 Se obtiene mucha incertidumbre en los datos, a cambio de una

reducción en el tiempo de ensayo.

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 62 / 65
 ALGUNAS GRÁFICAS DE LA DISTRIBUCIÓN

Se presentarán algunas curvas de la distribución de Weibull

PRESENTACIÓN DE FIGURAS
Se pueden importar figuras
Segunda columna
Comenzamos con las imágenes:
Texto en columnas.
1 Desarrollo más rápido,

Figura: El nombre de la figura

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 63 / 65
 REFERENCIAS

[1] EN 61165:2006
Application of Markov techniques. 2006.
[2] Kirkwood, James R. Markov processes. Advances in Applied
Mathematics. CRC Press, 2015
[3] Callejas-Cuervo, Mauro, Alarcón-Aldana, Andrea C. y Álvarez-Carreño,
Ana M. Modelos de calidad del software, un estado del arte. 2017
[4] Fuertes-Castro, José L. Calidad de Software. 2008

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 64 / 65
 Gracias

MAESTRÍA EN CONFIABILIDAD Y MANT

(UNIVERSIDAD PONTIFICIA BOLIVARIANA -MODELADO
BUCARAMANGA)
DE LA CONFIABILIDAD 65 / 65