CALCULO DE CONFIABILIDAD
La Distribucin de Weibull es por dems el modelo
estadstico ms popular del mundo para el uso de los datos
de vida.
Walodi Weibull naci el 18 de junio de 1887 .Su familia vino
originalmente de Schleswig-Holstein, en ese momento
estrechamente conectado con Dinamarca.
l era un alfrez, alumno en el Guardia de la Costa Sueco
Real en 1907, Capitn en 1916, y Mayor en 1940. Tom cursos
en el Instituto Real de Tecnologa donde ms tarde se hizo
profesor a tiempo completo (1924) y graduado ese mismo ao
1924. Su doctorado es de la Universidad de Uppsala en 1932.
Trabaj en las industrias suecas y alemanas como inventor
(bolas y rolineras de rodillo, y el martillo elctrico), y como
ingeniero consultor.
l public muchos trabajos sobre resistencia de materiales,
fatiga, ruptura en los slidos, rodamientos, y claro, la
distribucin de Weibull.
La familia matemtica del Valor Extremo de las
distribuciones, a las que ambos tipos Gumbel y Weibull
pertenece, es muy aplicable a los materiales, estructuras y
sistemas biolgicos porque tienen una creciente proporcin de
fracasos y pueden describir el uso de los procesos.
La Distribucin de Weibull se public por primera vez en
1939, hace alrededor de 60 (+10) aos y ha demostrado ser
invaluable para el anlisis de los datos de vida en lo
aeroespacial, automotor, potencia elctrica, potencia nuclear,
dental,
la
electrnica,
en
cada
industria.
(VER
[email protected] , bibliografa completa)
Wallodi Weibull 1887-1979
Photo by Sam C. Saunders
HISTORIAL DE FALLA:
El Historial de Fallas o el Histrico de Funcionamiento, de un material, permite determinar los
TEF (Tiempos Entre Fallas) o las duraciones de vida de los componentes, as como la
frecuencias acumuladas de fallas y repuestos utilizados, entre otros. (columna 7). (Ver Tabla 1)
HISTORIAL DE FALLAS
Columnas
2
PLANTA:
Cocimiento
5
EQUIPO:
Dosificador Equipo No. 086
FECHA DE
FECHA DE
TEF
TEF
TPR
OBSERVACIONES:
OCURRENCIA
ARRANQUE
(D)
(HR)
(HR)
DESCRIPCION DE LA FALLA
08/04/1996
08/04/1996
15/05/1996
15/05/199637 888 10.0 Cambiodepiones,salidainspectorvisual
26/06/1996
26/06/199641 984 10.0 Cambiochumaceras,entradainspectorvisual
05/07/1996
05/07/19969 216 7.0 Cambiodeengranes,salidadelalavadora
15/07/1996
15/07/199610 240 6.0 Cambiocadena7-1/2,delamesadecarga
27/07/1996
27/07/199612 288 6.0 Cambiochumacera.Lavadoradegaveras
16/08/1996
16/08/199619 456 4.0 Cadena3-1/4,deladesempacadora
14/09/1996
14/09/199628 672 8.0 Cambiodechumaceras,entradainspectorvisual
30/09/1996
30/09/199616 384 8.0 Cambiopionesdobles,Motor-reductor
10
12/12/1996
13/12/199672 1,728 12.0 Eslabonesdoblados,cadena7-1/2,entradallenadora
11
30/12/1996
30/12/199617 408 11.0 Cadenade3-1/4,entradallenadora
12
20/01/1997
20/01/199720 480 4.0 Cambiopionesdobles,transporadorbotellasllenas
13
09/02/1997
09/02/199719 456 8.0 Pionesyejes,transportadorsalidaempacadora
12.0 Cadena,salidadelaempacadora
Tabla 1
Explicacin de las siete (7)
columnas:
1.- Nmero de Registros, en nuestro caso son 13 registros.
2.- Fechas de la Ocurrencias de Falla, se inicia con el primer registro en la operacin del
equipo,
es decir, primer falla, y as sucesivamente, hasta la 13 falla.
3.- Fechas de Arranque, se inicia con el primer registro en la operacin del equipo, es decir,
arranque despus de la primer falla, y as sucesivamente, hasta el 13 arranque.
4.- Clculo de los das entre fallas, la fecha de ocurrencia de falla siguiente menos la anterior,
en
nuestro primer caso, 15/05/1996 08/05/96, resta del rengln 2 menos el 3, y as
sucesivamente, resultando 37 das en total, dato colocado en el rengln del segundo
registro. LOS CALCULOS SON HECHOS EN BASE A 30 DIAS POR MES.
5.-Los clculos normalmente son hechos en horas, columna 5, por ello, la columna 4 se
multiplica
por 24horas/da, es decir, 37 das* 24 horas/da = 888 horas, y as sucesivamente.
6.-Horas que dur cada reparacin a partir del primer registro de falla, para fines prcticos,
segn las Hojas de Trabajo para los clculos de Confiabilidad y Mantenibilidad, Tablas 2 y
3,
se tomarn solamente los primeros doce (12) registros de Tiempos Para Reparar (TPR), no
tomando
el ltimo
registro en cuenta, es decir, el nmero 13.
7.-En
este rengln
se especifican:
a) Descripcin de la falla
b) Repuestos utilizados
c) Mano de obra directa, indirecta y especializada.
d) Otras observaciones pertinentes a la falla, es nuestro caso la descripcin ha sido
muy escueta
HOJA DE TRABAJO PARA EL CALCULO DE CONFIABILIDAD (TABLA 2)
1. Esta columna presenta los nmeros de los
registros de falla de la tabla N 1, menos uno.
2. Esta columna presenta los valores de los (TEF),
calculados en la Tabla N 1, en la columna 6, de
menor a mayor, es decir en orden ascendente.
3. En esta columna se calcularn las Probabilidades
de Falla Estadstico (Pf) desde el registro 1 al 12,
con la siguiente formula Pf = i/(N+1), i =1,2,3,
.,N; siendo N el valor del ltimo registro (12),
este valor deber permanecer como una
constante de 4 cifras decimales Pf1 = 1/(12+1)
= 0.0769, se colocar en la columna con dos (2)
cifras decimales redondeadas. Pf1 =0.08; Pf2 =
0.015 (calculado como Pf2 = 2Pf1 = 2* 0.0769 =
0.1538 a ser redondeado a dos cifras decimales,
etc., y as sucesivamente hasta Pf12 = 12Pf1 =
4. 0.92
Esta columna presenta el clculo de la
Probabilidad de Supervivencia Estadstica Ps ,
con la formula de Ps1 = 1 Pf1, es decir, los
valores invertidos de Pf del rengln 12 al 1.
5. Esta columna presenta el clculo de la
Probabilidad de Supervivencia Estadstica Ps ,
que es el promedio de los Ps, cuando los
Tiempos entre falla son iguales, esto para
facilitar el graficado de los datos en el papel LogLog de la Distribucin de Weibull. En la tabla 2 se
repiten los valores de TEF, aproximadamente,
384 y 408, Ps=(0.69+0.62)/2, resultando
Ps= 0.655, en el caso de los tres datos
ligeramente diferentes (6,7 y 8) su promedio
Ps ser 0.46
Cocimiento
PLANTA:
EQUIPO:
TEF (HR)
PF
PS
PS
1216.0
0.08
0.92
2240.0
0.15
0.85
3288.0
0.23
0.77
4384.0
0.31
0.69
5408.0
0.38
0.62
0.66
6456.0
0.46
0.54
7456.0
0.54
0.46
0.46
8480.0
0.62
0.38
9672.0
0.69
0.31
10888.0
0.77
0.23
11984.0
0.85
0.15
121,728.0
0.92
0.08
7,200.0
TPEFa =
TEF
N
600.00 HR
Clculo del TPEFa ( Tiempo Promedio Entre Fallas
Analtico), es igual a la sumatoria de los TEF (7,200
Hrs)entreelnmeroderegistrosN(12)delatabla2,
iguala600.00Hrs.
DISTRIBUCINDEWEIBULLCALCULODECONFIABILIDAD
A.- AJUSTE GRAFICO LA DETERMINACION DE LOS
PARAMETROS:
La determinacin de tres los parmetros de Weibull,
permitir ajustar la Ley Probabilstica en la Distribucin de
WeibulldelaTabla2.Estoselograrllevandolospuntos(TEF,
Ps)o(TEF,Ps),sobreunGrficodeDistribucindeWeibull.
UnavezobtenidoslosvaloresdeTEFyPsoPs,estosse
graficarn en el Papel Logartmico de la Distribucin de
Weibull.
A partir de estos puntos se trazar un recta ajustada,
aproximada,lneaverde.
1. Paralela a esta lnea se trazar una lnea a partir del
puntoPivote(+),hastacruzarlaescaladeK,lneanegra
yobteneraselvalordelParmetrodeFormaK.
2. Apartirdelpercentil37%setrazaunalneaverticalhasta
llegar a la lnea ajustada a los puntos, a partir de este
cruce se traza una lnea horizontal, hacia la izquierda y
con ello se obtendr el valor de la Edad Caracterstica
paraFallaroParmetrodePosicinV.
1000
10
20
V=
640
30
40
10
80
85
90
93
95
96
96.5
98
99
100
93
95
96
96.5
98
99
100
T
E
F
100
50
20
30
40
50
80
85
90
3. Apartirdeestosvaloresydelvalorde ,delatablakvs
,obtenemoselTPEFgrfico,esdecir:
Obtenemos de la Grfica de Distribucin de
Weibulllossiguientesvalores:
K=1.88V=640hr
Probabilidad % Ps o
Ps
0.511.52.0
B.-CALCULOGRAFICODELTIEMPOPROMEDIOENTREFALLASY
CALCULODELAPROBABILIDADDESUPERVIVENCIAPs(t)OCONFIABILIDAD.
DelaTabladeKVs ,obtenemoselvalorde
iguala0.888ydelasiguienteformulaobtenemos
elvalordelTiempoPromedioEntreFallagrfico:
TPEFg = V *
TPEFg =
568.32 HR
Una vez calculado el valor grfico del TPEF,
procedemosacalcularelporcentajedeerrorentreel
valoranalticoyelgrfico,esdecir:
% Error = (TPEFg - TPEFa ) / TPEFa *100
%Error = (568.32 - 600.00)/600.00 * 100 = - 5.28%
Valor aceptable, ya que este clculo permite una
toleranciadeerrorde+/-20%.
Para el clculo de la Confiabilidad o Probabilidad de
Supervivenciautilizaremoslasiguienteecuacin:
Actual:ParaunTPEFa=600hr
PS (t ) =
e
t
V
Ps (TPEFa ) =
1
600hr 1.88
(
)
640hr
*100 = 41%
Paraobtenerun80%deeficiencia,deberemos
tener un prximo Mantenimiento Preventivo a
las290hr,esdecir:
Ps (TPEF) =
1
290hr 1.88
(
)
640
hr
*100 = 80%
"K"SELEEENLAGRAFICADEWEIBULL(TEFvs%Ps)
K
0.00
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
0.09
0.5
2.000
1.930
1.865
1.806
1.752
1.702
1.657
1.614
1.575
1.538
0.6
1.505
1.473
1.444
1.416
1.390
1.366
1.344
1.322
1.302
1.284
0.7
1.266
1.249
1.233
1.218
1.204
1.191
1.178
1.166
1.154
1.143
0.8
1.133
1.123
1.114
1.105
1.096
1.088
1.080
1.073
1.066
1.059
0.9
1.052
1.046
1.040
1.034
1.029
1.023
1.018
1.013
1.009
1.004
1.0
1.000
0.996
0.992
0.988
0.984
0.981
0.977
0.974
0.971
0.968
1.1
0.965
0.962
0.959
0.957
0.954
0.952
0.949
0.947
0.945
0.943
1.2
0.941
0.939
0.937
0.935
0.933
0.931
0.930
0.928
0.927
0.925
1.3
0.924
0.922
0.921
0.919
0.918
0.917
0.916
0.915
0.914
0.912
1.4
0.911
0.910
0.909
0.909
0.908
0.907
0.906
0.905
0.904
0.903
1.5
0.903
0.902
0.901
0.901
0.900
0.899
0.898
0.898
0.898
0.897
1.6
0.897
0.896
0.896
0.895
0.895
0.894
0.894
0.893
0.893
0.893
1.7
0.892
0.892
0.892
0.891
0.891
0.891
0.890
0.890
0.890
0.890
1.8
0.889
0.889
0.889
0.889
0.888
0.888
0.888
0.888
0.888
0.888
1.9
0.887
0.887
0.887
0.887
0.887
0.887
0.887
0.886
0.886
0.886
2.0
0.886
0.886
0.886
0.886
0.886
0.886
0.886
0.886
0.886
0.886
2.1
0.886
0.886
0.886
0.886
0.886
0.886
0.886
0.886
0.886
0.886
2.2
0.886
0.886
0.886
0.886
0.886
0.886
0.886
0.886
0.886
0.886
2.3
0.886
0.886
0.887
0.886
0.886
0.886
0.886
0.886
0.886
0.886
2.4
0.886
0.887
0.886
0.887
0.887
0.887
0.887
0.887
0.887
0.887
2.5
0.887
0.887
0.886
0.888
0.888
0.888
0.888
0.888
0.888
0.888
2.6
0.888
0.888
0.889
0.889
0.889
0.889
0.889
0.889
0.889
0.889
Resumiendo: Elparmetrode formaK permiteadaptar
lasratasdefallar(t),alasdiferentesfasesdelavidade
un sistema o de un componente, Curva de la Baera.
Igualmente, puede servir de indicador para un
diagnstico,yaquerepresentaelmecanismodefalla.En
nuestro caso K= 1.88, nos dice que el equipo esta en
etapadedesgasteporfatiga.
CONFIABILIDADYMANTENIBILIDAD.
ExtensindelconceptodeConfiabilidadaloscasosdesistemasreparables.Nocinparalela,poniendoen
evidencialasimilituddeclculosentrelaMantenibilidadylaConfiabilidad.
Paraelcasodeun
SistemaReparable.
(t)RatadeReparacin
r(t)RatadeFalla
MANTENIBILIDAD
CONFIABILIDAD
Probabilidad de duracin de
una buena reparacin
Probabilidad de un buen
funcionamiento.
TPEF
TPPR
TiempoPromedioEntreFallaso
debuenfuncionamiento.
TiempoPromedioParaRepararo
tiempotcnicoparareparar.
DISPONIBILIDAD
TPEF
D=
TPEF + TPPR
Probabilidad asegurar un
sistema requerido.
Losmodelosmatemticosestudiadosapropsitodelaconfiabilidadsern,porlotanto,utilizadosparala
cuantificacindelostiemposrelativosalaMantenibilidad.
HOJA DE TRABAJO PARA EL CALCULO DE MANTENIBILIDAD (TABLA 3)
1. Esta columna presenta los nmeros de los registros
de falla de la tabla N 1, menos uno, el ltimo TPR.
2. Esta columna presenta los valores de los Tiempos
Para Reparar (TPR), relacionados en la Tabla N 1,
de menor a mayor, es decir en orden ascendente.
HOJA DE
TRABAJO
TPR (HR)
3. En esta columna se calcularn las Probabilidades
de Falla Estadstico (Pf) desde el registro 1 al 12,
con la siguiente formula Pf = i/(N+1), i =1,2,3,
.,N; siendo N el valor del ltimo registro (12), este
valor deber permanecer como una constante de 4
cifras decimales Pf1 = 1/(12+1) = 0.0769, se
colocar en la columna con dos (2) cifras decimales
redondeadas. Pf1 =0.08; Pf2 = 0.015 (calculado
como Pf2 = 2Pf1 = 2* 0.0769 = 0.1538 a ser
redondeado a dos cifras decimales, etc., y as
sucesivamente hasta Pf12 = 12Pf1 = 0.92
4. Esta columna presenta el clculo de la Probabilidad
de Falla Estadstica Pf , que es el promedio de los
Pf, cuando los Tiempos Para Reparar son iguales,
esto para facilitar el graficado de los datos en el
papel Semi-Log de la Distribucin de Gumbel Tipo I.
En la tabla 3 se repiten los valores para TPR
renglones 1 y 2, Pf=(0.08+0.15)/2, resultando Pf= 0.12,
en el caso de los dos renglones 3 y 4 su promedio
Pf ser 0.27, en los renglones 6 y 7 su promedio
ser Pf 0.50, de los renglones 8 y 9 su promedio
0.65 , finalmente los renglones 11 y 12 ser 0.88
CALCULO DE MANTENIBILIDAD
PF
PF
1 4.0
0.08
2 4.0
0.15
0.12
3 6.0
0.23
4 6.0
0.31
0.27
5 7.0
0.38
6 8.0
0.46
7 8.0
0.54
0.50
8 10.0
0.62
9 10.0
0.69
0.65
10 11.0
0.77
11 12.0
0.85
12 12.0
0.92
0.88
98.0
TPPRa =
TPR
N
= 8.17hr
ClculodelTPPRa(TiempoPromedioParaReparar
Analtico), es igual a la sumatoria de los TPR (98.0
Hr)entreelnmeroderegistrosN(12)delatabla3,
iguala8.17hr.
DISTRIBUCINDEGUMBELTIPO1CALCULOMANTENIBILIDAD
A.- AJUSTE GRAFICO LA DETERMINACION DE LOS
PARAMETROS:
La determinacin de los parmetros de Gumbel,
permitirajustarlaLeyProbabilsticaenlaDistribucinde
GumbeldelaTabla3.Estoselograrllevandolospuntos
(TPR, Pf) o (TPR, Pf), sobre un Grfico de Distribucin
deGumbelTipo1.
Apartirdeestospuntossetrazarunrectaajustada,
aproximada,lneaverde.
Apartirdelpercentil37%(correspondeaL=0)se
trazaunalneaverticalhastallegaralalneaajustadaa
los puntos, a partir de este cruce se traza una lnea
horizontal,hacialaizquierdayconelloseobtendrel
valor de la Edad Caracterstica para Reparar o
Parmetro de Posicin , Valor que seala el tiempo
idealpararealizarunareparacinalequipo.
Obtenemos de la Grfica de Distribucin de
GumbelTipo1,lossiguientesvalores:
L1=0L2=2
=6.8hrP1=6.8hrP2=12.5hr
P2=12.5
TPR(hr)
Una vez obtenidos los valores de TPR y Pf o Pf,
estos se graficarn en el Papel Semi-Logartmico de la
DistribucindeGumbelTipo1.
20
=6.8
P1=
9
8
7
6
5
4
3
2
1
10
9
8
7
6
5
4
3
2
1
10
20
30
40
50
70
75
80
85
90
93
95
96
96.5
98
99
100
-2
4
-1
5
CALCULO DEL TIEMPO PROMEDIO PARA REPARAR GRAFICO TPPRg:
Con los datos obtenidos, procederemos a
reemplazarlosenlasecuacionescorrespondientespara
elclculodelTPPRgysuporcentajedeerror:
1. ClculodelapendientePdelarectaajustada:
P=
P P2 P1 12.5 6.8
=
=
L L2 L1
20
P = 2.85hr
2.
Calculada la pendiente se calcular el valor del
Parmetro de Forma a, igual al inverso de la
pendienteP.
1
a = = 0.35hr 1
P
3. AplicaremoslaformuladeTPPRg:
0.5778
a
0.5778
TPPRg = 6.8hr +
= 8.45hr
0.35hr 1
TPPRg = +
4. Calcularemos el porcentaje de Error, el cual
deberestarentreel+/-10%:
TPPRg TPPRa
% Error =
*100
TPPRa
8.45hr 8.17 hr
*100
8.17hr
Error = 3.43%
% Error =
A continuacin calcularemos los valores de la
Mantenibilidad(ProbabilidaddeFalla):
PF (T < t ) =
ee
a ( t )
MantenibilidadActual:
PF (T < TPPRa ) =
ParaTPPRa(Actual)
*100
1
e
e 0.35 ( 8.176.8 )
PF(T)=
*100
54%
Para una mantenibilidad ptima, se requerir de
un %Pf alrededor del percentil 37%, por tanteo
obtendremos un 38% para ello utilizaremos un
valordeTPPRligeramentemayorque :
PF (T < 6.9) =
ParaTPPR(Optimo)
e 0.35 ( 6.96.8 )
PF(T)=
*100
Formulaparalacalculadora:
e(e( a * (t ))) *100
38%
CALCULO DE LA DISPONIBILIDAD ACTUAL, SEGN:
DISPONIBILIDAD
TPEF
D=
TPEF + TPPR
Probabilidad asegurar un
sistema requerido.
600hr
D=
*100
(600 + 8.17)hr
D = 98.66%
Estevalornosdaunainformacinuntantoerrnea,dado
queelTiempoPromedioParaRepararesmuchomenorque
elTiempoPromedioEntreFallas,porellomeinclinomsal
ClculodelaEfectividadGlobaldelequipo,esdecir:
Paraobtenerundatorealdelaefectividadtotaldelequipo,
sehandecalculartresvariablesseparadas:
1.LaDisponibilidad
2.LaEficienciadeDesempeo,y
3.LaTasadeCalidaddelosproductos.
EFICIENCIA DE DESEMPEO (Prdidas de Velocidad)
Estavariableesiguala:
%Efd=TasadeoperacinrealXTasanetadeoperacin
Ejemplo:
Unamquinaestprevistaparaobtener300piezasen
una hora, pero realmente se estn obteniendo 200
piezasporhora.
Tasa de operacin real : Tor
Tiempo de ciclo terico sersuproduccinterica,
serde300piezas/60min=5piezas/min,suciclo
tericoserigualaTt=1/5=0.2min/pieza.
Tiempo de Ciclo real, sera su produccin real
200 piezas/60 min= 3.33 piezas/min , es decir,
Tr=1/3.333=0.3min/pieza,finalmente:
Tiempo _ de _ Ciclo _ Terico
*100
Tiempo _ de _ ciclo _ real
0.2
Tor =
*100 = 66.67%
0.3
%Tor =
Tasa de operacin neta: Ton
La Tasa de operacin neta nos da una idea de
cmo se mantiene una velocidad de operacin en
unperododado.
%Ton =
Tiempo _ de _ proceso _ real
*100
Tiempo _ total _ de _ operacin
Sabemos que se estn produciendo 5 piezas por
minuto, si el equipo ha estado produciendo en un
tiempo real de operacin de 300 minutos, restando
tiempodelasparadasimprevistasde110min,node
reparaciones,tendremos:
Cantidadprocesada=5piezas/min*300min=1000
piezas.
Efectividad Global del equipo ser:
Eg = Disponibilidad x Eficiencia de Desempeo x
Tasa de Calidad
E g = D * E fd* Tc
Eg = 0.987 * 0.488 * 0.97 *100
Eg = 46.73%
Si tomamos en cuenta el tiempo de operacin total
delequipode410min,tendremos:
Estosignifica que el equiposeest utilizando muypor
debajodesuscondicionesidealesqueson:
Tiempo del proceso real =Cantidadprocesadax
tiempodecicloreal.
Disponibilidad:
mayordeun90%
EficienciadeDesempeo mayordeun95%
TasadeCalidad:
mayordeun99%
1,000 piezas * 0.3 min/ pieza
*100
410 min
Ton = 73.2%
Ton =
%Efd=TasadeoperacinrealXTasanetadeoperacin
Efd=66.7x73.2=48.8%
Tasa de Calidad (prdidas por defectos):
Supondremosunatasadecalidaddelosproductos:
Tc=97%
Finalmente, la efectividad global del equipo
debersermayorde:
Eg = 0.90 * 0.95 * 0.99 *100
Eg = 85%
Enelejemploanalizadosenotaqueelvalormsalejadode
las condiciones ideales es el referente a la Eficiencia de
Desempeo(enlasprdidasporvelocidaddeproduccin),
esalldondesedebenconcentrartodoslosesfuerzosenla
bsqueda de solucionespara mejorar laEfectividadGlobal
delEquipo.