Análisis de
confiabilidad y
mantenibilidad
ME57A ./me57a/confiabilidad.ppt
Rodrigo Pascual J.
16 de agosto de 2007
14 de marzo de 2007
23-25,30 de agosto de 2006
19-31 de agosto de 2005/1,6 de abril de 2005
1 de septiembre de 2004/5 de abril de 2004
Dr. Rodrigo Pascual
4 de septiembre de 2003/16 de abril de 2003 Curso Mantención de Maquinas
Al terminar de leer este capitulo,
podrás:
Definir y calcular:
confiabilidad, Evaluación
tasa de fallas Síntesis
Análisis
Aplicación
Comprensión
Conocimiento
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Motivación
Diseño de un programa eficiente de
mantenimiento
comprensión de los fenómenos de falla
Fallas
aleatorias,
Estudiaremos
Ingeniería de la confiabilidad
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Objetivos
Reducir costo global
controlar y mejorar la confiabilidad,
Definir programas
Preventivos
predictivos
Reemplazo de equipos
Agrupamiento de intervenciones
overhaul
Indicadores
…
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Dificultades
incertidumbre
efectos de cambios en el comportamiento
Escasez de datos
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Tipos de intervenciones
Tipo
Correctivas
Preventivas
M.P. Sistemático
centrado en la condición
Calidad
Perfectas
Como nuevo
Mínimas
Como antes
Imperfectas
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Distribuciones
eventos discretos
(número de fallas..)
eventos continuos
cantidades físicas
masa, temperatura,..
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Distribuciones
1 parámetro
Exponencial
Poisson
2 parámetros
Normal
3 parámetros
Weibull
mas
Bi-Weibull
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Distribución exponencial
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Tasa de fallas
l
tasa de fallas por unidad de tiempo
1/l
tiempo medio entre fallas (MTBF)
Para una pob. exponencial
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Distribución de Poisson
Describe
numero de eventos aleatorios
Parámetro
promedio de eventos/intervalo de tiempo
l
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Poisson
probabilidad de que k eventos ocurran en el
intervalo T:
m=lT m
Valor esperado de x t
T
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Probabilidad acumulada
En excel,
POISSON(x;media;acumulado)
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
Horas de seguimiento
Equipo Stand-by Operación
1 2800 400
2 300 2500
3 2700 700
4 2500 800
5 0 3100
total fallas 7 12
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
Calcule las tasa de fallas de los equipos.
total horas standby = 2800 + 300 + 2700 + 2500 + 0
= 8300
nro de fallas en modo standby = 7
7 fallas/hora op.
tasa de fallas en modo standby =
8300
= 8:4E ¡ 4 fallas/ hora
Horas de seguimiento
Equipo total horas
Stand-by operación l= 400 + 2500 + 700 + 800 + 3100
Operación Op.
1 2800 400 Stand-by
2 300 2500 = 7500
3 2700 700
nro 4de fallas2500en operación
800 = 12
5 0 3100 12 Tiempo
tasa de fallas en operación =
total fallas 7 12 7500
= 1:6E ¡ 3 fallas/ hora
Dr. Rodrigo Pascual 7 + 12
tasa de falla global = de Maquinas
Curso Mantención
8300 + 7500
nro de fallas en modo standby = 7
7
tasa de fallas en modo standby =
8300
Solución = 8:4E ¡ 4 fallas/ hora
total horas operación = 400 + 2500 + 700 + 800 + 3100
= 7500
nro de fallas en operación = 12
12
tasa de fallas en operación =
7500 fallas/hora op.
= 1:6E ¡ 3 fallas/ hora
Horas de seguimiento 7 + 12
Equipo tasa de falla
Stand-by global
Operación =
1 2800 400 8300 + 7500
2 300 2500 = 1:2E ¡ 3 fallas/ hora
3 2700 700
4 2500 800 8:4E ¡ 4 + 1:6E ¡ 3
tasa
5 de falla0 promedio
3100 =
2
total fallas 7 12 = 1:2E ¡ 3 fallas/ hora
ser vación 2 Notese que la tasa de falla
Dr. Rodrigo global y la tasa de falla promedio
Pascual
Curso Mantención de Maquinas
son necesariamente identicas. Si se está interesado en usar un solo valor, la
l agregado
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Distribución gaussiana
Densidad de probabilidad
m , media
, varianza
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Probabilidad acumulada
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Obs
F(x)
se evalúa con cambio de variable
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejercicio
probabilidad de que
una máquina no falle durante un día
si
Promedio: 10 fallas/semana laboral (5 días)?
l=10 fallas/semana laboral
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Solución
m=lT
T:1 día
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejercicio
tasa de falla: 1/10000 fallas/ut;
probabilidad de que falle entre las 200 y 300
ut?
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Solución
0,0001
l=10-4 fallas/ut
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejercicio
Decisión:
¿Que modelo de resistencia comprar?
Criterio:
Duración de una resistencia
P([420, 720])>= 95%
Proveedor propone modelo:
Media
600 horas
desviación standard
120 horas
¿Recomienda el modelo?
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Solución
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
= 0. Sea N ¡ n el numero de component es que
Confiabilidad
n(t)
R(t) =
N
st r ibución acumulada de fallas
ción acumulada de falla F (t) se de…ne como la
falle en el int ervalo (0; t). Ent onces:
Edad
R(t) + F (t) = 1
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Obs
Disponiendo de la tasa de fallas
Numero esperado de fallas en un cierto
intervalo
Costos de falla asociados
Costos de intervención esperados
presupuesto
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Observaciones
La probabilidad de que falle en un
intervalo infinito es 1
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Observaciones
estimando l(t)
confiabilidad del equipo
plazos
de mantenimiento
preventivo
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Métodos para estimar R a partir de
historial de fallas
Estimar
directamente de los datos de falla,
R,l
Sin censuras (suspensiones)
l(t) o:preventiva
x:correctiva
o o
x
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Tiempo calendario
en el instante ti
quedan
n - i unidades i= 1 2 3 4
operando.
Edad
x x x x
tiempo
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Estimación inicial
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Por lo tanto,
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Obs
O sea,
posibilidad
nula de
que hayan unidades
operando para t > tn.
poco probable
que muestra incluya
el tiempo de
supervivencia más
largo,
se subestima R
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Además,
Es razonable que
las primeras y las F=i/n
ultimas observaciones,
tengan la misma
distancia con respecto
al 0% y 100% de
posibilidad
simetría
25% 100%
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Estimación de F (II)
-rangos medios-
F=i/n
F=1/n 25% 100%
F=1/(n+1) 20% 80%
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
III. rangos de la mediana
Importancia de las colas
F=i/n 25% 100% Rangos medianos aprox.
F=i/(n+1) 20% 80%
tabla 16% 84%
F=(i-0.3)/(n+0.4) 16% 84%
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Tabla de rangos medianos
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo, n=8
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Desviaciones c/r a método de la
mediana (tabla)
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
tasa de falla definida por tramos
al
¸ (t)
0
0 Dr. Rodrigo Pascual
Curso Mantencióntde tiempo
b Maquinas tb tu
Confiabilidad
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Modelo de Dhillon
0.9
0.8
0.7
0.6
l (t)
0.5
0.4
0.3
0.2
0.1
0 2 4 6 8 10
Dr. Rodrigo Pascual tiempo
Curso Mantención de Maquinas
MTBF si estrategia correctiva
x x x x
TBF0 TBF1 TBF2 TBF3
tiempo calendario
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
componente
confiabilidad
linealmente decreciente
La confiabilidad es 1
R(0)=1 R
R(10000)=0
Calcule su MTBF. 10000 t
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
tasa de fallas constante
-pob. exponencial-
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
Si
l= 2E-6 fallas/hora,
Calcule
confiabilidad a las 500 horas de operación
MTBF
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
a las 500 horas:
R(500) = exp(2E-6 * 500) = 0.999
MTBF = 1/2E-6 = 50000 horas
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Desgaste mecánico
Se Curso
evalúa numéricamente
Dr. Rodrigo Pascual
Mantención de Maquinas
Ejemplo
Si
b = 2E-6 fallas/h^2
a = 1E-7 fallas/h,
Calcule MTBF
edad
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
en Maple:
> MTBF:=int(exp(-(.5*1e-7*t^2+2e-6*t)),t=0..infinity);
MTBF = 3943 horas
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Mantenibilidad
Tasa de reparación
Análogo
Tasa de fallas:
Numero esperado de reparaciones por unidad
de tiempo
-reparaciones/unidad de tiempo-
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Com p. A T TTR
Ejemplo fa lla 1
fa lla 2
40.1
83
1.5
3.8
Com p. B T TTR
fa lla 1 41.4 1.3
Com p. C T TTR
fa lla 1 40.6 1.1
fa lla 2 82 1.5
TTR t im e t o re p a ir
T t ie m p o in ic io fa lla
H o riz o n t e 100
u n id a d e s h o ra s
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
Com p. A T TTR
fa lla 1 40.1 1.5
fa lla 2 83 3.8
Com p. B T TTR
fa lla 1 41.4 1.3
Com p. C T TTR
fa lla 1 40.6 1.1
fa lla 2 82 1.5
TTR t im e t o re p a ir
T t ie m p o in ic io fa lla
H o riz o n t e 100 Dr. Rodrigo Pascual
u n id a d e s h o ra s Curso Mantención de Maquinas
Distribución de Weibull
en estudios de confiabilidad
sistemas mecánicos.
Ventajas
muy flexible,
adaptable a una variedad de observaciones
experimentales.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ley de Weibull, 3 parámetros
l(t)=
adimensional
ut
ut
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Tasa de fallas de Weibull
h=2 g=0
b=3
l(t)=
l(t)
b=1
b=.5
Edad del componente
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Modelo de Weibull
h=2
1.6
1.4
1.2 b=3
3 parámetros 1
0.8
l
0.6 b=1
0.4
0.2 b=.5
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
tiempo
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Weibull
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Función Gamma
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Obs
caso g = 0, b = 1
exponencial
l= 1/h
Para b > 3
distribución normal.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Si g=0
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Weibull
Estimación grafica
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Estimación de F(t)
método de rangos medianos si la población
es pequeña:
F(i) = (i – 0.3)/(n + 0.4)
método de rangos medios:
F(i) = i/(n + 1)
i :indice de la observación (ti < ti+1).
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Aplicación practica
1. Obtener n observaciones, ordenar
2. Estimar la función de distribución F(t)
3. Calcular pares (X,Y)
4. Graficar
5. Ajustar la mejor recta
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Weibull
Si g0
cambio de variable,
t’=t-g
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Curva de Weibull para g > 0
Zona
imposible ·
log t
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo, población de
componentes
Un grupo de rodamientos tuvieron las
siguientes duraciones:
801 312 402 205 671 1150 940 495 570
Se desea conocer la confiabilidad para una
vida de 600 horas y el MTTF.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Hace en excel
En matlab
» t=[205 312 402 495 570 671 801 940 1150];
» F=.1:.1:.9;
» X=log(t);
» Y=log(log(1./(1-F)));
» P=polyfit(X,Y,1);
» beta=P(1)
beta=1.79
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
En matlab
» eta=exp(P(2)/(-P(1)))
eta =
715.97
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
1
0.5
0
log(log(1/(1-F)))
-0.5
-1
-1.5
-2
-2.5
5.2 5.4 5.6 5.8 6 6.2 6.4 6.6 6.8 7 7.2
log t
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Matlab (2)
weibplot([205 312 402 495 570 671 801 940 1150])
Weibull Probability Plot
0.96
0.90
0.75
0.50
Probability
0.25
0.10
0.05
3
10
Data
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
En maple
>MTBF=716*GAMMA(1+1/1.79)
Entonces MTTF = 636.9 horas
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo g0 , vida de 0.11. MODELOS DE CONFIABILIDAD
componentes i
1
2
V i da
2175
2800
F (i )(%)
5
10
3 3300 15
4 3800 20
5 4250 25
6 4650 30
7 5250 35
8 5840 40
9 6300 45
10 6700 50
11 7150 55
12 7800 60
13 8500 65
14 9200 70
15 10500 75
16 11000 80
17 12600 85
18 14000 90
19 15800 95
./me57a/datos/datos-conf1.xls Tabla 2: Dat os del ejercicio
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ajuste para g=0
g=0
1.5
0.5
log(log(1/(1-F)))
-0.5
-1
-1.5
-2
-2.5
-3
7.5 8 8.5 9 9.5 10
log t
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Estudio
0.3
·
0.28
0.26
norma del vector residuo
0.24
0.22
mínimo
|Ax-B| 0.2
0.18
0.16
0.14
Obs: gamma puede
0.12
800 900 1000 1100
g
1200 1300 1400 Ser negativo
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ajuste óptimo ✔
g=1280
1.5
g=0 1
1.5
1 0.5
0.5
0
log(log(1/(1-F)))
0
log(log(1/(1-F)))
-0.5
-0.5
-1
-1
-1.5 -1.5
-2
-2
-2.5
-2.5
-3
7.5 8 8.5 9 9.5 10
log t -3
6.5 7 7.5 8 8.5 9 9.5 10
log t
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Verificación de modelos
al imponer una ley
se incurre en algún error,
queremos que el riesgo sea menor:
definimos el nivel de confianza
probabilidad de que el modelo sea erróneo.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Test 2
Condición
al menos n = 50 observaciones.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Pasos
1. Se agrupan las observaciones.
Debe haber al menos 5 observaciones en cada
grupo.
Los intervalos para definir los grupos no son
necesariamente de la misma longitud.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Test 2
El test se basa en
diferencias entre
nro de observaciones en cada grupo y
nro pronosticado
por la ley seleccionada.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Test 2
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Test 2
E tiene una distribución 2 con grados de
libertad:
=r-k-1
donde
k = 1 para la ley exponencial,
k = 2 para la ley normal,
k = 3 para la ley de Weibull
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Test 2
hipótesis
“observaciones siguen modelo propuesto”
es rechazada si:
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
24
i T BF (horas) ni
Supóngase que para un grupo de 1 0-500 7
equipos similares se han 2 500-1000 8
observado los siguientes TBF: 3 1000-1500 9
4 1500-2000 10
5 2000-2500 12
6 2500-3000
P 8
54
Tabla 3: Grupos de…nidos para el
Dr. Rodrigo Pascual n=
Obser vación 5 nde¢p
Curso Mantención i es el numero de onservac
Maquinas
Ejemplo
Hipótesis
leyexponencial.
tasa de fallas l= 1/1600 fallas/hora.
Se desea realizar un test
confianza 95%
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
Modelo propuesto
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
La probabilidad de que una observación caiga
en los grupos definidos en la tabla es
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Test de aceptación
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
n = 54
= 6 - 1 - 1 = 4
tabla 2 entrega
(4;0.95) = 9.49 Ver tabla
en Matlab
>> chi2inv(0.95,4)
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
2
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
E > 2(4,0.95)
Se rechaza hipótesis
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Test Kolmogorov-Smirnov (KS)
se puede aplicar para cualquier numero de
observaciones n.
si n es grande
mejor agrupar las observaciones y
usar el test 2 .
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Test KS
compara
la función de distribución experimental con
El modelo propuesto
se usan
valores absolutos de las diferencias entre punto
y punto.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Kolmogorov-Smirnov
Sea la distribución muestreada y F(t)
la distribución propuesta.
La discrepancia para ti es:
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Kolmogorov-Smirnov
Puede demostrarse que la distribución de
Dn = max(Dni )
depende solo de n; y se puede escribir
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
Evidencia
TBF (días):
23,16,56,71,4,25,51,30
Modelo
Gaussiana
media 34
desviación standard 22,
Test con = 5%?
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Probabilidad
según población hipotética
P(t < 4) = P ((4 - 34)/22)= 0.086
En Excel
=DISTR.NORM.ESTAND((4-34)/22)
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Kolmogorov-Smirnov
Dn = max(Dni )
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Kolmogorov-Smirnov
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Kolmogorov-Smirnov
Según la tabla
Dn = 0.127
Dn, para n = 8, = 0.05 es
D8,0.05 = 0.457
se acepta la hipótesis.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Kolmogorov-Smirnov
1.4
1.2
Probabilidad acumulada de falla
0.8
0.6
0.4
0.2
-0.2
-0.4
0 10 20 30 Dr. Rodrigo
40 Pascual 50 60 70 80
Edad (ut)
Curso Mantención de Maquinas
Vida remanente esperada
Confiabilidad condicional
1
0.9
t
0.8
0.7
R t0
0.6
0.5
0.4
0.3
0 0.2 0.4 0.6 0.8 1
Dr. Rodrigo Pascual
Edad
Curso Mantención de Maquinas
Vida remanente esperada
Confiabilidad
condicional
probabilidad que el
componente
sobreviva t ut mas
dado que
ha sobrevivido t0 ut
desde que estaba
como nuevo
Dr. Rodrigo Pascual
ut:unidad de tiempo Curso Mantención de Maquinas
Vida remanente esperada
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Weibull de 2 parámetros
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Ejemplo
b=1/2,1,2,3
h= 100 ut
t0=3/2MTTF(b,h)
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Maple
>beta:=0.5;eta:=100;t0:=1.5*eta*GAMMA(1+1/beta);
> f0:=beta/eta*((t+t0)/eta)^(beta-1)*exp(-((t+t0)/eta)^beta)/exp(-(t0/eta)^beta);
>MRL:=int(t*f0,t=0...infinity);
> plot(exp(-((t+t0)/eta)^beta)/exp(-t0/eta)^beta),t=0..2*t0);
Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Comentarios
R(t)
MTBF
l(t)
Datos históricos Análisis de
MTTR
b,h Weibull
(sin suspensiones) confiabilidad g vida asegurada/
predesgaste
vida remanente
✔
esperada
confiabilidad
condicional
Edad Dr. Rodrigo Pascual
Curso Mantención de Maquinas
Bibliografía
R. Pascual. El Arte de Mantener (draft).
Universidad de Chile, 2007
P. Lyonnet. Maintenance Planning,
Methods and Mathematics. Chapman &
Hall, 1991.
Dr. Rodrigo Pascual
Curso Mantención de Maquinas