Tesis Analisis de Weibull
Tesis Analisis de Weibull
Tesis Analisis de Weibull
PRESENTADO POR.
CARLOS AMILCAR ZAMORA LOPEZ
ASESOR
ING. ARTURO GAMARRA CHINCHAY
CALLAO PERU
2012
DEDICATORIA:
A mi Madre Graciela, Esposa Roxana,
mis hijos Alondra y Leonardo, por ser el
motor de mi vida.
INDICE
1. INTRODUCCION Pg.
2. OBJETIVOS
2.1 OBJETIVO GENERAL.
2.2 OBJETIVO ESPECIFICO.
3. ATLAS COPCO
3.1 VISION.
3.2 VALORES FUNDAMENTALES.
3.3 ATLAS COPCO EN EL MUNDO.
3.4 PRINCIPALES EVENTOS HISTORICOS.
3.5 ORGANIZACIN MUNDIAL DE ATLAS COPCO.
5. EVALUACION TECNICA-ECONOMICA.
5.1 CASO ROTURA DE CABLE DE AVANCE EQUIPOS DE
PERFORACION.
5.1.1 HORAS TRABAJADAS DE LA FLOTA EQUIPOS DE
PERFORACIN.
5.1.2 COSTO DE MANTENIMIENTO GLOBAL PARA UNA
POLTICA CORRECTIVA.
5.1.3 COSTO DE MANTENIMIENTO GLOBAL PARA UNA
POLTICA PREVENTIVA.
5.1.4 MTBF
5.1.5 MTBI (MEAN TIME BETWEEN INTERVENTIONS) TIEMPO
MEDIO ENTRE INTERVENCIONES.
5.1.6 COSTO GLOBAL HORARIO PARA UNA POLTICA
CORRECTIVO.
5.1.7 COSTO GLOBAL HORARIO PARA UNA POLTICA
PREVENTIVA.
5.1.8 GRAFICO DE CONFIABILIDAD VS COSTOS.
5.1.9 CALCULO AHORRO USANDO POLTICA PREVENTIVA.
6. CONCLUSIONES Y RECOMENDACIONES.
6.1 CASO ROTURA DE CABLE DE AVANCE EQUIPOS DE
PERFORACION.
6.1.1 CONCLUSIONES.
6.1.2 RECOMENDACIONES.
6.2 CASO ROTURA PERNOS CILINDRO VOLTEO SCOOPTRAM.
6.2.1 CONCLUSIONES.
6.2.2 RECOMENDACIONES.
7. INFORMES TECNICOS.
7.1 ANALISIS CAUSA RAIZ (FALLA BOMBA AGUA BOOMER 282).
7.2 FALLA EN SISTEMA DE DIRECCION ST1030 NUMERO
INTERNO S-604, CON N/S AVO 10 X 083.
7.3 ANALISIS COSTO BENEFICIO LABIO CUCHARA VS UAS EN
SCOOPTRAM ST-1030.
7.4 ROTURA CARDAN Y TRANSVERTER SCOPTRAM ST-7 S-406.
7.5 ANALISIS JACK NIFE EN EQUIPO DE PERFORACION S1D J-
124.
7.6 ANALISIS CONFIABILIDAD PEDAL ACELERACION EN PRUEBA
EN SCOOPTRAM ST-7.
7.7 FALLA EN MOTOR DIESEL CAMION MT2010 AVO 11X005 D-
207.
8. BIBLIOGRAFIA.
9. APENDICES Y ANEXOS.
9.1 INTRODUCCION A LA INGENIERIA DE CONFIABILIDAD.
9.2 FUNDAMENTO ESTADISTICO.
9.3 DATOS Y TIPOS DE DATOS.
9.4 LIMITES DE CONFIANZA.
9.5 LA DISTRIBUCION DE WEIBULL.
1. INTRODUCCION.
1. OBJETIVO GENERAL
3.1 VISION.
La visin del Grupo Atlas Copco es ser la Primera Opcin y la Primera
Eleccin para todos nuestros clientes.
Sistema de Avance
de un Equipo de
Perforacin.
Se obtiene la ecuacin:
t
F (t ) 1 e (2)
Tomando el logaritmo natural de ambos lados de la ecuacin:
t
ln1 F (t )
t
ln ln1 F (t ) ln
o:
ln ln1 F (t ) ln ( ) ln(t ) (3)
Donde:
N : Es el numero de datos o fallas.
: Este factor nos indica si la falla es prematura, aleatoria o
por desgaste.
: Este factor nos indica la vida caracterstica.
: Este factor indica cuanto tiempo el componente puede
trabajar sin falla, es decir con una confiabilidad de 100%.
MTBF : Es el tiempo promedio entre fallas.
4.4.4 TABLA Y GRAFICA DE CONFIABILIDAD (ROTURA DE
CABLES DE AVANCE).
Donde:
N : Nro. de datos (fallas y censurados).
: Indica si la falla es prematura, aleatoria o por desgaste.
: Indica la vida caracterstica.
: Indica el periodo que el componente puede trabajar sin
falla.
5. EVALUACION TECNICA-ECONOMICA.
5.1 CASO ROTURA DE CABLE DE AVANCE EQUIPOS DE
PERFORACION.
Cgc = $ 7 495.70
Cgc = $ 243.40
5.1.4 MTBF
Para la distribucin Weibull de tres parmetros, el MTBF se
calcula as:
MTBF 242,5
Ts Ts 1 1
0
R(t )dt inc ,
Para Ts <
MTBI (Ts) = Ts
Para Ts >
(T ) 1 1
MTBI * DISTR.GAMMA(( s ) ; ;1;1) * EXP (GAMMA.LN ( ))
CostoTotalCorrectivo
C gc
MTBF
$7495.70
C gc
242.50
$
C gc 30,91
hr
C gp xR (Ts ) C gc xF (Ts )
Cg
MTBI (Ts )
Con esta ecuacin tenemos lo siguiente:
El costo mas
bajo lo
tenemos con
Aplicar una
una poltica
poltica
preventiva en
preventiva con
50 horas, a
una frecuencia
11,19 $/hr
de 125 horas
tambin tiene
sentido, pues
nuestro costo
horario seria
15,33 $/hr.
Muy por
debajo de los
30,91 $/hr que
cuesta una
poltica
correctiva.
5.1.8 GRAFICO DE CONFIABILIDAD VS COSTOS.
En esta zona
podemos aplicar
nuestro plan
preventivo, a bajo
costo y alta
confiabilidad.
5.1.9 CALCULO AHORRO USANDO POLTICA PREVENTIVA.
En esta seccin vamos a calcular el ahorro usando poltica
preventiva, para eso vamos a comparar el costo horario global de
usar poltica correctiva ($30,91 $/hr) menos el costo global para
una poltica preventiva, dependiendo de la frecuencia que
seleccionemos, para entender mejor se adjunta la siguiente tabla.
Zona de mximo
ahorro y alta
confiabilidad.
Cgc = $ 52 944,66
Cgp = $ 25 976,16
5.2.4 MTBF.
CostoTotalCorrectivo
C gc
MTBF
$52 944,66
C gc
2791,70
$
C gc 18,97
hr
Por lo tanto el costo horario para una poltica correctiva de cambio
de pernos de cilindro de volteo es de 18,97 $/hr.
5.2.7 COSTO GLOBAL HORARIO PARA UNA POLTICA
PREVENTIVA.
C gp xR (Ts ) C gc xF (Ts )
Cg
MTBI (Ts )
1. CONCLUSIONES.
2. RECOMENDACIONES.
1. CONCLUSIONES.
Para esta falla de rotura de pernos de cilindro de volteo, no
podemos aplicar mantenimiento preventivo basado en el uso
porque el factor (1,14) se encuentra muy cerca de 1. Es mejor
modelarlo como falla aleatoria.
Comunicar al fabricante Atlas Copco del desgaste prematuro,
pues el MTBF (2791,70 horas) es menor al MTBF que
recomienda fabrica 6000 horas.
Despus de este informe Fabrica realiz una campaa de cambio
de cilindro de volteo, con la tapa mejorada, esperemos mejoras
en el MTBF.
2. RECOMENDACIONES.
Realizar un plan de inspecciones semanales para los pernos del
cilindro de volteo, porque la falla es aleatoria (puede suceder en
cualquier momento), las inspecciones nos ayudara a ver el estado
del cilindro hidrulico.
Redisear cilindro de volteo, por parte de fbrica Atlas Copco
Suecia (est en ejecucin cambio de cilindros de volteo).
7. INFORMES TECNICOS.
En esta seccin vamos a presentar los informes ms
resaltantes, que han originado cambios en la Gestin de Mantenimiento
durante mi experiencia profesional.
QUE OCURRIO?
La Bomba de agua, presenta sonido extrao y no genera presin.
CUANDO OCURRIO?
El da 09-01-2011 a las 15:15 horas.
DONDE OCURRIO?
En el NV 920 BY 02 ACC 12.
FRECUENCIA.
02 en el ao.
IMPACTO.
Perdida de produccin, 123.10 TON / HORA.
Prdida de Bomba de Agua.
Menor productividad del equipo.
Menores indicadores de gestin de mantenimiento.
ANALISIS CAUSA RAIZ BOMER 282 BAJA PRESION AGUA PERFORACION.
BOOMER 282 BAJA
PRESION AGUA
PERFORACION.
CAMBIAR POR
BOMBA
DISEADA PARA
3500 RPM, 440
VOLTIOS Y 60 HZ.
RECOMENDACIONES Y CONCLUSIONES.
Cambiar de diseo de la bomba FLY actual.
Si no es posible cambiar de diseo la bomba, regresar con la marca
anterior, para disminuir la frecuencia de fallas de 03 @ 01 por cada ao.
El cliente por este tipo de bomba est perdiendo en productividad $117
094.06 cada ao.
68
7.2 FALLA EN SISTEMA DE DIRECCION ST1030 NUMERO
INTERNO S-604, CON N/S AVO 10 X 083.
ANTECEDENTES:
13-08-2010, Se entrega el equipo con los parmetros generales
del equipo segn los STD de fbrica.
En diciembre 2010, Personal tcnico de Lima, Eliot Vargas
sugiere cambiar el Block de vlvulas quedando la maquina
inoperativa del 25-12-2010 @ 04-01-2011.
04-01-2011, El taller de Lima presta 01 Block de direccin,
dejando el equipo operativo, sin embargo el operador segua
reportando fallas en la direccin.
20 y 21-01-2011, Se realiza MP de 2000 horas, en esta parada de
equipo se cambia el Block de vlvulas completo y nuevo, se
entrega la maquina operativa.
03-02-2011, El operador Mamani reporta que la direccin sigue
con los mismos problemas a pesar que se haba cambiado el
Block de vlvulas nuevo, personalmente ingreso a mina con los
tcnicos Abraham Sarzo, Ivn Quiquia y Carlos Cortabrazo,
encontramos al equipo trabajando en el NV1020, prob la
maquina en mnimo y a su rpm mximo, efectivamente cuando
69
estaba girando a la derecha se perdi la direccin a medio
recorrido de los cilindros, se probo tambin la maquina sin carga y
presentaba el mismo problema, en ese instante decid parar el
equipo realizamos pruebas que indica fabrica y todos los
resultados hidrulicos indican que la maquina esta dentro del
STD, sin embargo el problema de direccin persiste.
04 y 05-02-2011, En mina se encontraban dos tcnicos de Lima
(Manuel Cisneros y Jhon Llamoca) entregando equipos (ST1030 y
MT2010), nos dan el apoyo realizando algunas pruebas, al inicio
tenan la percepcin que era un problema del operador, sin
embargo cuando empezaron a realizar mas pruebas con nuestro
tcnico Javier Camarena la falla se presento, verificaron
presiones de todo el sistema hidrulico encontrndose dentro del
STD, sacaron los cilindros hidrulicos y las lainas pero el
problema persiste.
06-02-2011, La maquina continua inoperativa, estamos revisando
la articulacin central.
A continuacin detallo los parmetros que se tomaron cuando se
instalo el Block de vlvulas que enviaron de Lima, este trabajo se
realizo el 04-01-2011 con los tcnicos Abraham Sarzo y Aquilino
Aguedo.
70
71
CAUSA FISICA DE LA FALLA
CONCLUSIONES.
La presin del sistema de direccin esta dentro del STD 3250 PSI, el
caudal de las bombas esta tambin dentro del STD, inclusive muy por
encima de lo que dice fabrica, por ejemplo segn fabrica a 750 rpm el
caudal debe ser de 33.86 LPM, con nuestro caudalmetro da 40 LPM.
Finalmente solicitamos el apoyo de tcnicos especialistas en
SCOOPTRAM ST1030 a nuestra central en Lima.
72
7.3 ANALISIS COSTO BENEFICIO LABIO CUCHARA VS UAS EN
SCOOPTRAM ST-1030.
ANTECEDENTES
EL 01-06-2011 en coordinacin con la Jefatura de Mantenimiento
Equipo Pesado (Ing. Lucio Lara) se decide instalar un juego de uas al
scooptram S-605 con la finalidad de realizar un estudio y ver si era ms
rentable que los labios de cuchara.
73
3. Se est considerando el precio del metal a $180.00 TON, este precio
es variable de acuerdo al metal que procesan las compaas mineras
y al mercado externo.
4. La produccin de la maquina se considera variable, en el desarrollo
del estudio presentamos una grafica de Beneficio VS la produccin.
5. Se ha considerado una distribucin normal para los datos
estadsticos, porque las fallas por desgaste se aproximan mucho a
esta distribucin, sin embargo para una mayor precisin se puede
usar la distribucin de Weibull, pero necesitamos ms datos (min 20).
6. No se ha trabajado con los promedios (media) porque para una
distribucin normal significa el 50% de probabilidad que esto pase,
por eso es que se ha ajustado los valores menores al promedio para
trabajar con una mayor confiabilidad.
7. Se ha considerado monto fijo del soldador $ 2300.00 mensual y un
solo soldador para realizar el trabajo.
74
CASO 01: CAMBIO DE LABIO CUCHARA
Segn datos de cambio de labio de cuchara tenemos los
siguientes datos:
75
Vamos a calcular la probabilidad de que los equipos trabajen 400
horas por mes.
76
CASO 01: CAMBIO DE LABIO CUCHARA
Segn datos de cambio de labio de cuchara el cambio de labio de
cuchara nos tomo los siguientes tiempos (hrs).
77
COSTO MATERIALES (REPUESTOS E INSUMOS)
78
CASO 02: CAMBIO DE UAS DE CUCHARA
Segn nuestra data estadstica tenemos los siguientes datos.
79
2. Costo del mineral producido, para nuestro caso estamos
considerando que la Tonelada de Zinc tiene un valor de $ 180.00
(ciento ochenta dlares americanos)
Por lo tanto podemos decir que el impacto a la produccin es.
80
Se concluye que los beneficios estn relacionados directamente a
la produccin, eso significa que a mayor produccin de la maquina
(TON/HR) mayor es el beneficio, y nuestra ecuacin de beneficio es la
siguiente:
B=6244.20P-11 332.82
Esta ecuacin ser graficada a continuacin.
81
RECOMENDACIONES Y CONCLUSIONES.
82
7.4 ROTURA CARDAN Y TRANSVERTER SCOPTRAM ST-7 S-406.
ANTECEDENTES
25-10-2011 A las 2:45 p.m. reportan rotura de cardan, el tcnico Juan
Solano al inspeccionar la maquina observa que el cardan de salida del
transverter esta empotrado en el transverter, intenta retirar el cardan y
en su evaluacin final observa que el cardan y transverter estn
averiados, presentan rotura.
26-10-2011 @ 27-10-2011 No se puede trasladar equipo a taller por falta
de apoyo de scooptram.
83
28-10-2011 Se remolca equipo a superficie y se desmonta
TRANSVERTER, realizado por Javier Camarena, Martin lvarez y Ral
Almonacid.
RECOMENDACIONES
84
SEGUNDO, SE
ROMPIO ESTE
PERNO
PRIMERO, SE
ROMPIO ESTE
PERNO
TERCERO, SE ROMPIO
YOKE DEL CARDAN N/P
5572 2651 00
CUARTO, LA
ROTURA DEL
CARDAN OCASIONO
LA ROTURA DEL
TRANSVERTER
85
ZONA DE LA FALLA
86
7.5 ANALISIS JACK NIFE EN EQUIPO DE PERFORACION S1D J-
124.
ANTECEDENTES
En el ao 2011 la flota de equipos de perforacin ha sido la ms
desfavorable con los indicadores de gestin (disponibilidad mecnica,
utilizacin, MTBF, MTTR) por eso es que se est realizando un estudio
de las fallas que se han presentado durante el ao 2011 y presentar un
plan de accin para las mas criticas.
87
A continuacin Presentamos el Grafico de Dispersin de Tiempos
(Jack Knife).
Cuadrante ms
Crtico por prdida
de Disponibilidad,
Confiabilidad y
Mantenibilidad.
88
Antecedentes:
Si observamos el detalle de las fallas del motor diesel se observa que en
Mayo el equipo estuvo inoperativo por 194.75 horas y a fines de
Octubre @ primer da de Noviembre estuvo inoperativo 265.50 horas.
Antecedentes:
Este equipo est trabajando en la zona baja en el NV920 Bypass 02,
zona con problemas de alimentacin de energa.
89
Si observamos el detalle de las fallas en abril 2010 el equipo estuvo
inoperativo 70 horas por cambio de motor elctrico y en julio estuvo
inoperativo 28 horas por cambio de motor elctrico.
Antecedentes:
Esta falla crnica nos est haciendo perder confiabilidad y disponibilidad
(la mantenibilidad est debajo del promedio), sin embargo como la
mantenibilidad es baja nos hemos acostumbrado a pensar que la falla
est controlada, pero no es as, porque si vemos el detalle de las fallas,
se observa un buen porcentaje de cambio de mangueras hacindonos
perder confiabilidad.
90
Capacitacin a los operadores en el manejo defensivo (evitar
accidentes).
Antecedentes:
Esta falla tambin podramos considerarlo como crnica, porque est
bajo control, sin embargo nos est quitando mantenibilidad, si
observamos el detalle de las fallas se observa que las mangueras se
cambian por rotura y se han cambiado en promedio 3.50 horas, si el
taller de mangueras estara en mina se cambiaria como mximo en 2
horas, mejorando nuestra mantenibilidad, por ende la disponibilidad
mecnica.
91
7.6 ANALISIS CONFIABILIDAD PEDAL ACELERACION EN
PRUEBA EN SCOOPTRAM ST-7.
SENSOR PEDAL
ACELERACION N/P
5580 0065 81
92
pedales de aceleracin no han fallado, por lo tanto simularemos que han
fallado para poder calcular su curva de confiabilidad.
93
Con estos datos podemos calcular los parmetros de forma y el
tiempo medio entre fallas.
94
A continuacin presentamos los parmetros de forma y el tiempo
medio entre fallas, para el nuevo pedal de aceleracin.
95
pedal nuevo para 800 horas de operacin nos da una confiabilidad de
42% es decir hay 58% probabilidad de falla.
2. Y as sucesivamente podemos observar cada punto de la grafica.
< 1, falla
prematura
96
A continuacin calcularemos la perdida que le ocasionaremos a
nuestro cliente, si continuamos usando el nuevo pedal de aceleracin,
salvo se corriga algunas problemas tcnicos, para mejorar la
confiabilidad de nuestro nuevo pedal de aceleracin, necesitamos dos
datos para poder calcular la perdida.
97
Ahora estamos en capacidad de calcular los costos horarios para
cada caso.
RECOMENDACIONES Y CONCLUSIONES.
98
7.7 FALLA EN MOTOR DIESEL CAMION MT2010 AVO 11X005 D-
207.
Por medio del presente se informa de LA FALLA DEL MOTOR
DIESEL montado en el CAMION MT2010 D-207.
ANTECEDENTES
99
25-12-2011 en el turno noche los tcnicos Juan Sosa (especialista
motores Deutz), Marco ngel (Capacitado en motores Cummins)
y Roger Guivin, adems se entrego programa INSITE para una
mejor evaluacin, detectando baja presin de aceite con el programa
INSITE, intentan arrancar el motor para ver qu valor de presin de
aceite est registrando pero el motor no arranca y se detecta un ruido
extrao en el turbocompresor, al inspeccionar el turbocompresor se
detecta que la turbina esta averiada, se deja el equipo inoperativo por
falta de turbocompresor y baja presin aceite motor diesel (alarma
roja).
26-12-2011 CIA VOLCAN nos entrego 01 turbocompresor usado de
ST1030 (usan el mismo modelo de motor QSL9), los tcnicos Ral
Almonacid y Vctor Mercado montan el turbocompresor, al intentar
dar arranque el motor se observa dos advertencias 01 mbar y 01
roja, se analiza nuevamente con el Insite dndonos cdigo activo
0143 (baja presin aceite motor), se intenta resetear cdigo de falla
pero no se puede permanece activo, posteriormente se verifico filtro
de aceite, vlvula termosttica de aceite y by pass del filtro
encontrndose en buen estado, se verifico ingreso de aceite al filtro y
se observo que no tiene aceite, se comunica a la Jefatura de
Mantenimiento que la falla del motor se ha complicado y es necesario
remolcar a superficie, paralelo a ello se gestiona subida de tcnicos
especialistas en motores CUMMINS.
27-12-2011 Los tcnicos Vctor Mercado y Aquilino Aguedo remolcan
el equipo al taller de superficie.
28-12-2011 A las 2:30 p.m. llegan los tcnicos de MITSUI Oscar
Nolasco Corpus y Carlos Torres Martnez, se comunica los trabajos y
evaluacin que hemos realizado, a final de guardia me comunican
que destaparan toda la sincronizacin para poder sacar la bomba de
aceite, porque el filtro de aceite estaba sin aceite.
28-12-2011 A las 7:00 p.m. se consulta al jefe de guardia Ismael
Untiveros por los operadores, donde me informa que el da 25-12-
2011 dio orden al operador Abel Inga que trabaje con el camin de
100
4:00 @ 12:00, el operador Inga no trabajo hasta las 11:15 a.m.
cuando decide prender el camin y sacarlo a superficie, a las 13:00
horas el operador entrante Percy Huaroc arranca el motor diesel y al
observa que la alarma roja no se apagaba decide apagar el motor,
despus de unos minutos intenta nuevamente arrancar el motor
diesel y observa que la alarma sigue prendida reportando a la central
de emergencias para un auxilio mecnico por parte de Atlas Copco.
29-12-2011 Personal de MITSUI da con la falla, al desmontar la
bomba de aceite se detecta que esta averiada, el eje de la bomba se
desplazo ocasionando que no exista transferencia de movimiento del
cigeal hacia la bomba, se observa que no hay juego axial en el
cigeal y partculas metlicas en el aceite, por eso recomiendan
demonstrar el motor diesel.
RECOMENDACIONES Y CONCLUSIONES.
101
Eje de engranaje de
bomba de aceite fuera
de su lugar.
102
8. BIBLIOGRAFIA.
103
9. APENDICES Y ANEXOS.
9.1 INTRODUCCION A LA INGENIERIA DE CONFIABILIDAD.
104
raro encontrar problemas con definiciones de falla o diferencias de
unidad de tiempo que invalidan completamente los resultados de
pruebas y anlisis caros y que consumen bastante tiempo.
9.1.2 ESTIMACION
105
tomamos 1,000 bolas, podramos obtener los siguientes resultados
para el nmero de bolas negras: 44.5% y 49.5% en cada ensayo. En
este caso, notamos que nuestra estimacin para el porcentaje de
bolas negras tiene un rango mas angosto, 5%. Usando esto,
podemos ver que cuanto ms grande es el tamao de la muestra,
ms angosto es el rango de la estimacin y, probablemente, ms fina
es la estimacin del verdadero valor.
106
menos confiables e inadecuados programas de prueba pueden
originar productos con fallas ocultas que se presentarn en el campo.
Un ahorro rpido en el corto plazo por el uso de componentes ms
baratos o por tamao de muestras de prueba pequeos normalmente
producirn costos ms altos a largo plazo en forma de costos de
garanta o prdida de confianza del cliente. Debe lograrse un
equilibrio apropiado entre la Confiabilidad, la satisfaccin del cliente,
tiempo en el mercado, ventas y rasgos. En la figura 1.1 se ilustra este
concepto. El polgono en la izquierda representa un proyecto
apropiadamente equilibrado. El polgono en la derecha representa un
proyecto en el que se ha sacrificado la Confiabilidad y satisfaccin
del cliente por causa de las ventas y tiempo en el mercado.
http://www.weibull.com/articles/rellntro/index.htm
107
9.1.3.2 RAZONES CLAVES DE LA INGENIERIA DE
CONFIABILIDAD.
108
Ingeniera de Confiabilidad, xito que puede verse en la seguridad
del viaje areo-comercial de hoy.
109
proceso de produccin, a travs de su vida de uso prctica, con el
apoyo de mantenimiento y disponibilidad. La Ingeniera de
Confiabilidad cubre:
La Confiabilidad.
La Mantenibilidad.
La Disponibilidad.
Las tres reas pueden cuantificarse numricamente con el uso de
los principios y anlisis de datos de vida de la Ingeniera de
Confiabilidad.
Esta curva se traza con la vida del producto en el eje X y con la tasa
de fallas en el eje Y. La vida puede estar en minutos, horas, aos,
ciclos, actuaciones o cualquier otra unidad cuantificable de tiempo o
uso. La tasa de falla esta dada como las fallas entre las unidades
sobrevivientes por unidad de tiempo. Como puede verse de esta
grfico, muchos productos empezarn sus vidas con un tasa de falla
alta (debido a los defectos de fabricacin, pobre habilidad del
personal que lo utiliza, control de calidad pobre de partes entrantes,
etc.) y exhibe una tasa de falla decreciente. La tasa de falla
normalmente se estabiliza a una tasa aproximadamente constante
en la regin de vida til dnde las fallas observadas son fallas de
oportunidad. Cuando los productos experimentan ms uso y
desgaste, la tasa de falta empieza a subir ya que la poblacin
empieza a experimentar fallas relacionadas al desgaste. En el caso
de mortalidad humana, la tasa de mortalidad (la tasa de falla), es
ms alta durante el primer ao de vida, luego cae a un nivel
110
constante bajo durante nuestra adolescencia y la vida adulta
temprana y luego se eleva cuando progresamos en aos.
Trabajo Inicial
111
implementarse calculando la Confiabilidad ptima de tal producto. La
figura 1.3 describe tal escenario. El costo del producto total es la
suma del costo de produccin y diseo as como los costos post-
despacho. Puede verse que a un nivel de Confiabilidad ptimo, el
costo del producto total es mnimo. El nivel de Confiabilidad ptimo
es el que coincide con el costo total mnimo durante toda la vida del
producto.
112
Gua con respecto a las decisiones de accin correctiva para
minimizar las fallas y reducir los tiempos de mantenimiento y
reparacin que se eliminarn antes y despus del diseo.
Ayuda proporcionando pautas sobre prcticas de control de
calidad.
Optimizacin de la meta de Confiabilidad que debe disearse
en los productos y sistemas para un costo total mnimo de
propiedad, operacin y mantenimiento durante la vida.
Habilidad para dirigir estudios de intercambio entre parmetros
como la Confiabilidad, Mantenibilidad. Disponibilidad, Costo,
Peso, Volumen, Operatividad, Serviciabilidad y Seguridad
obtener el diseo ptimo.
Reduccin de costos de garanta o, por el mismo costo,
aumento en el periodo y monto de garanta.
Establecimiento de pautas por evaluar a los proveedores
desde el punto de vista de la Confiabilidad del producto.
Promocin de ventas en base a los ndices de Confiabilidad y
mtrica a travs de los departamentos de ventas y
comercializacin.
Aumento de la satisfaccin del cliente y un aumento de ventas
como resultado de la satisfaccin del cliente.
Aumento de las ganancias o, para la misma ganancia,
provisin de ms productos y sistemas confiables.
Promocin de imagen positiva y reputacin de la compaa.
113
sobre las fallas. Es importante tener un sistema de anlisis til,
de interpretacin y de realimentacin en todas las reas de la
compaa que se tratan con el producto desde su nacimiento
hasta su muerte
114
9.1.4 LA CONFIABILIDAD Y EL CONTROL DE LA CALIDAD.
115
9.2 FUNDAMENTO ESTADISTICO.
Variables aleatorias
En general, la mayora de los problemas en ingeniera de
confiabilidad tratan con medidas cuantitativas, como el tiempo para la
falla de un componente, o si el componente falla o no falla. Para juzgar si
un componente es defectuoso o no, slo hay dos resultados posibles.
Podemos denotar una variable aleatoria X como representante de estos
posibles resultados (es decir defectuoso o no defectuoso). En este caso,
X es una variable aleatoria que slo puede asumir estos valores.
116
9.2.1.1 DENSIDAD DE PROBABILIDAD Y LAS FUNCIONES DE
DENSIDAD ACUMULADAS.
Definiciones
Si X es una variable aleatoria contina, entonces la funcin de
densidad de probabilidad, pdf, de X, es una funcin f(x) tal que
para dos nmeros, a y b (a < = b) se representa con:
b
P(a X b) f ( x)dx and f ( x) for all x
a
(1)
Es decir, la probabilidad que X asuma un valor en el
intervalo [a, b] es el rea bajo la funcin de densidad.
117
Note que los lmites de integracin dependen del dominio
de f(x). Por ejemplo, para todas las distribuciones consideradas,
este dominio sera 0, , , o y, En el caso de
x
F ( x) f ( s)ds
(3)
f ( x)dx 1
118
Dnde es la media y es la desviacin estndar. La
distribucin normal es una distribucin de dos parmetros, es
decir con dos parmetros y .
1 t ' '
1 2 '
f (t ) e
t . 2
(4)
119
probabilidad que el tiempo para la falla est en la regin de 0 (o y)
y t. As de la Ecuacin (4):
t
F (t ) Q(t ) f ( s)ds
o,
Q(t ) R(t ) 1
R(t ) 1 Q(t )
t
R(t ) 1 f ( s)ds
,
R(t ) t f ( s)ds
Recprocamente
d ( R(t ))
f (t )
dt
120
MISSION: Mission 1 Mission 2
TIME:
time = t time = t
RESULT: SUCCESS UNKNOWN
R(T t )
R(t / T )
R(T )
f (t )
R (t )
La funcin tasa de falta tiene unidades de fallas por unidad
de tiempo entre sobrevivientes, por ejemplo 1 falla por mes.
121
9.2.1.4 FUNCION DE VIDA MEDIA.
m t * f (t )dt (5)
0,
122
9.2.1.6 FUNCION DE MODO
d f (t )
0 (7)
dt
Para una distribucin continua, la moda es el valor de la
variable que corresponde a la densidad de probabilidad mxima
(el valor dnde la pdf tiene su valor mximo).
9.2.2 DISTRIBUCIONES
123
Una de las distribuciones ms simples y ms comnmente
utilizadas (y a menudo errneamente utilizada debido a su
simplicidad), es la distribucin exponencial.
La pdf de la distribucin exponencial est matemticamente
definida como:
fi (t ) e t
1 t 2
1
f (t ) e 2
2
124
preocupe sobre cmo se estiman estos parmetros, esto se
discutir en las siguientes secciones). Por ahora, asuma que
podemos asignar valores a estos parmetros basados en los
datos.
fi (t ) e t
As, la funcin de Confiabilidad exponencial puede
derivarse como:
t s
R(t ) 1 o e ds
1 1 e * t
e *t
f (t )
(t )
R(t )
e (t )
e (t )
t f (t ) dt
0
125
t e t dt
0
Parmetro de Escala
126
Parmetro de Forma
El parmetro de forma, como implica el nombre, ayuda a
definir la forma de una distribucin. Algunas distribuciones, como
la exponencial o normal, no tiene un parmetro de forma ya que
tienen una forma predefinida que no cambia. En el caso de la
distribucin normal, la forma es siempre la forma de una
campana. El efecto del parmetro de forma en una distribucin se
refleja en las formas de la pdf, la funcin de confiabilidad y la
funcin tasa de fallas.
Parmetro de Ubicacin
El parmetro de ubicacin se usa para cambiar la posicin
de una distribucin en una u otra direccin a lo largo del eje
horizontal. El parmetro de ubicacin, normalmente se representa
como y, y define la situacin del origen de una distribucin y
puede ser positivo o negativo. En trminos de las distribuciones
de tiempo de vida, el parmetro de ubicacin representa un
cambio de tiempo.
.
Esto significa que la inclusin de un parmetro de ubicacin
para una distribucin cuyo dominio normalmente es 0,
cambiar al dominio y, , dnde y puede ser positivo o
127
se sienten incmodas con el concepto de un parmetro de
ubicacin negativo, el cual establece que el estado de falla ocurre
tericamente antes del tiempo cero. Realmente, el clculo de un
parmetro de ubicacin negativo es indicativo de fallas inmviles
(fallas que ocurren antes de que un producto se use por primera
vez) o de problemas con la fabricacin, empaquetado o durante el
proceso de envo. Se debe prestar ms atencin al concepto del
parmetro de ubicacin en las siguientes discusiones de la
distribucin exponencial y Weibull que son las distribuciones de
tiempo de vida que ms emplean el parmetro de ubicacin.
Distribucin de Weibull.
La distribucin Weibull es una distribucin de Confiabilidad
de propsito general utilizada para modelar la resistencia del
material, el tiempo para la falla de componentes electrnicos y
mecnicos, equipo o sistemas. En su caso ms general, la pdf
Weibull de tres parmetros est definido por:
1 t
t
f (t ) e
+
Con parmetros , y y, dnde = parmetro de forma,
= parmetro de escala y y = parmetro de ubicacin.
128
Si el parmetro de ubicacin y, se asume cero, la
distribucin se convierte en una distribucin Weibull de dos
parmetros o:
1 t
t
f (t ) e
Distribucin Exponencial
La distribucin exponencial normalmente se usa para
componentes o sistemas que exhiben una tasa de falla constante
y est definido en su caso ms general por:
f (t ) e (t )
(8)
(Tambin conocido como exponencial de dos parmetros),
con dos parmetros, a saber y y.
Si el parmetro de ubicacin y, se asume cero, la
distribucin se vuelve exponencial de un parmetro o,
f (t ) e t
Se presentan la distribucin exponencial y sus
caractersticas en ms detalle en el captulo correspondiente.
129
Distribucin normal
La distribucin normal se utiliza para el anlisis de
Confiabilidad general, el tiempo para la falla de componentes
electrnicos y mecnicos simples, equipo o sistemas.
La pdf de la distribucin normal esta da por:
1 t 2
1
f (t ) e 2
2
f (t ) 0, t , 0
(9)
Donde = media de los tiempos normales para falla, =
desviacin estndar de los tiempos para la falla.
Se presenta la. distribucin normal y sus caractersticas en
ms detalle en el captulo correspondiente.
130
del modo de falla bajo investigacin, principalmente a travs de
datos histricos o fsicos de falla.
S t t
fs (t ) pi i i 1e ( i ) i
i 1 i i
131
Distribucin Log normal
La distribucin log normal se usa para el anlisis de
Confiabilidad general, ciclos para la falla en la fatiga, resistencia
de materiales y en el diseo probabilstico de cargas variables.
f (t ) O, t O, T ' O
T ' 1n (t )
(10)
132
parmetros, con parmetros , y . La pdf de la distribucin esta
dada por,
ln(t )
1 *
ln( t ) e
* 2
i 1
f (t ) * * e if 0
* t T 1
2
1 ln(t ) 2
1 2
e if 0
t * 2
133
Distribucin Gamma
z ln(t )
(11)
Distribucin Logstica
ez
f (t )
(1 e z ) 2
134
t
z
0
(12)
Donde:
parmetro de ubicacin (tambin denotado como F),
parmetro de escala.
Distribucin Loglogistic.
ez
f (t )
t (1 e z ) 2
T ' '
z
f (t ) 0, t 0, T 0,
T ' ln (t )
(13)
Donde,
' = parmetro de escala, T = parmetro de forma
135
Distribucin Gumbel
La distribucin Gumbel, tambin llamada distribucin del valor
extremo ms pequeo (SEV o del tipo 1), es apropiada para modelar
esfuerzos que a veces se sesgan a la izquierda (pocas unidades
dbiles fallan a baja tensin, mientras que las restantes fallan a
tensiones ms altas). La distribucin Gumbel tambin podra ser
apropiada para modelar la vida de productos que experimentan un
desgaste muy rpido luego de alcanzar una cierta edad.
t
z
f (T ) 0, 0
(14)
Trazado de probabilidad
El mtodo menos intensivo matemticamente hablando para la
estimacin de parmetros es el mtodo de trazado de
probabilidad. Como implica el trmino, el trazando de probabilidad
136
involucra plotear los datos de probabilidad en un papel construido
especialmente para este fin. Este mtodo se realiza a mano.
T
Q(T ) 1 e
T
ln(1 Q (T ) ln e
T
ln(1 Q(T )
T
ln( ln(1 Q(T ) ) ln
1
ln ln ln (T ) 1n( )
1 Q (T )
(15)
Reemplazando y por:
1
y ln ln
1 Q (T )
y:
137
x ln(7' )
La ecuacin puede volverse a escribir como,
y x ln( )
138
probabilidad incluyen un indicador de la pendiente para simplificar
este clculo). El valor de la pendiente es el parmetro .
As a T :
Q(T ) 1 e
1 e 1
0.632
63.2%
139
necesitamos obtener el porcentaje cumulativo de falla para cada
tiempo de falla. En este ejemplo, y para 10 horas, el porcentaje
cumulativo de falla es 25%, para 20 horas es 50%, y as
sucesivamente. Este es un simple mtodo que ilustra la idea. El
problema con este mtodo simple es el hecho que el punto 100%
no est definido en la mayora de ploteos de probabilidad, as que
debe usarse una alternativa y un enfoque ms consistente. El
mtodo ampliamente utilizado para determinar este valor es el
mtodo del rango medio para cada falla. Este mtodo se discute
luego.
El rango medio
Se usa el rango medio para obtener una estimacin de la no
confiabilidad, Q(Tj), para cada falla. Es el valor que da la
verdadera probabilidad de falla, Q(Tj), que debe tener al fallar la j-
esima falla de una muestra de N unidades a un 50% nivel de
confianza. Esto significa esencialmente que sta es nuestra
estimacin ptima para la no confiabilidad. La mitad del tiempo el
verdadero valor ser mayor que el 50% estimado de confianza, la
otra la mitad del tiempo, el verdadero valor ser menor. Esta
estimacin esta basada en la solucin de la ecuacin binomial.
kN Z k (1 Z ) N k
N
P
k j
140
El rango medio se obtiene resolviendo esta ecuacin para Z a P =
0.50.
kN Z k
N
0.50 (1 Z ) N k
k j
(16)
141
Esta aproximacin del rango medio tambin es conocida como la
aproximacin de Benard.
Teora fundamental
142
mnimos cuadrados, o regresin lineal, porque la regresin se
realiza en los valores del rango, ms especficamente, los valores
del rango medio (representado en el eje y).
Rango de regresin en Y
Asuma que un juego de pares de datos (x1,y1), (x2,y2),..., (XN,YN) se
obtienen y trazan, y que los valores X son conocidos exactamente.
Entonces, segn el principio de los ,mnimos cuadrados, que
minimiza la distancia vertical entre los datos y la lnea recta
trazada, la mejor lnea recta para estos datos es la lnea
recta y a b
x (donde los smbolos introducidos ( ) indican que
ese valor es una estimacin) tal que:
N N
(a b xi y2 ) 2 min( a, b) (a bxi yi ) 2
i 1 i 1
Y donde a y b son las estimaciones de los mnimos cuadrados de
143
N N
yi xi _ _
a i 1 b i 1 y b x
N N
N N
N xi yi
xi yi i 1 i 1
N
b i 1
2
N
xi
N
2 i 1
xi N
i 1
Rango de regresin en X
b
lnea recta y a y tal que:
N N
(a b yi xi ) 2
min( a, b) (a byi xi ) 2 +
i 1 i 1
b
De nuevo, a son las estimaciones de los mnimos cuadrados
144
N N
xi yi _ _
a i 1 b i 1 x b y
N N
(19)
N N
N xi yi
N i , yi i 1 i 1
N
b i 1
2
N
yi
N
i 1
yi N
2
i 1
(20)
Las relaciones correspondientes por determinar los parmetros
para las distribuciones especficas (es decir Weibull, exponencial,
etc.), se presenta en los captulos correspondientes.
El Coeficiente de correlacin
El coeficiente de correlacin es una medida de qu tan bien
la regresin lineal carga los datos y normalmente se denota por .
En el caso de anlisis de datos de vida, es una medida de la
fuerza de la relacin lineal (correlacin) entre el rango medio y los
datos. El coeficiente de correlacin de la poblacin est definido
como sigue:
xy
x y
Dnde xy covarianza de x e y, x desviacin estndar de x, y
y = desviacin estndar de y.
dado por,
145
N N
N xi y i
xi yi i 1 i 1
N
i 1
N
2
N
2
xi yi
N N
x 2 i 1 y 2 i 1
i N i N
i 1 i 1
146
datos completos, censurados a la izquierda, censurados a la derecha
(o suspendidos) y datos en intervalo.
Teora fundamental
La idea bsica detrs del MLE es obtener valores de los parmetros
ms probables, para una distribucin dada que describan mejor los
datos.
f ( x;1 , 2 , ...., k )
147
x1, x2 , ... x R que corresponde en el caso de anlisis de datos de vida
a tiempos de falla. La funcin de probabilidad esta da por:
R
L(1 , 2 , .... k x1 , x2 , ...., x R ) L f ( xi ,1 , 2 , ... k )
i 1
i = 1, 2, , R
R
ln L ln f ( xi ; 1, 2 , ...., k )
i 1
0, j 1,2, ..., k
j
148
Comentarios sobre el Mtodo MLE
149
estn irregularmente distribuidos. Tambin puede proporcionar
estimaciones con una o ninguna falta observada, lo que el rango de
regresin no puede hacer.
150
9.3 DATOS Y TIPOS DE DATOS.
151
datos de vida en dos categoras: completos (est disponible toda la
informacin) o censurados (alguna informacin no est). Este captulo
detalla los mtodos de clasificacin de datos.
Los datos completos son mucho ms fciles de trabajar que con los
datos censurados. Por ejemplo, sera mucho ms difcil calcular el
promedio de las notas de los estudiantes si el conjunto de datos no
estuviera completo, es decir el promedio de las notas 30, 80, 60, 90,
95, tres notas mayor que 50, una menor de 70 y una nota que est
entre 60 y 80. (nota: = 70.8 para este conjunto de datos.)
152
Censura a la derecha (Suspendido)
Censura en intervalo
153
Censura a la izquierda
Datos agrupados
154
Dependiendo del mtodo de anlisis escogido, es decir de regresin
o de mxima probabilidad, Weibull++ trata de agrupar los datos
diferentes.
155
anlisis de los artculos suspendidos tambin sera considerado
cuando:
156
de las suspensiones. Si consideramos estas cinco unidades, se
utilizar la siguiente metodologa:
157
Una vez que se ha establecido el nmero de orden medio
para cada falla, obtenemos las posiciones del rango medio para
estas fallas a su nmero de orden medio. Especficamente,
obtenemos el rango medio del orden de los nmeros 1, 2.25 y
4.125 de la muestra de 5, como se dada a continuacin.
Una vez que se han obtenido los valores del rango medio,
el anlisis de la probabilidad trazada es idntico a lo presentado
anteriormente. Como podra haber notado, esta metodologa es
bastante laboriosa. Se han desarrollado otras tcnicas y atajos
para reestructurar este procedimiento. Ms detalles sobre este
mtodo, ver Kececioglu [20].
158
Para los datos dados sobre los resultados es como sigue:
Los parmetros estimados usando el mtodo descrito son los
mismo para ambos casos (1 y 2):
0.81
11,417 hr
1.33
6,900 hr
21,348 hr
159
9.4 LIMITES DE CONFIANZA.
160
Si vuelve a poner las diez bolas en la piscina y repite de
nuevo el proceso, podra conseguir seis bolas negras, lo que cambia
su estimacin a 60% de bolas negras. Cul de las dos es correcta?
Ambas estimaciones son correctas! Repitiendo este experimento
una y otra vez, podra averiguar que esta estimacin normalmente
est entre X1% y X2%, y puede asignar un porcentaje al nmero de
veces que su estimacin cae entre estos lmites. Por ejemplo, nota
que el 90% de veces esta estimacin est entre X1% y X2%
De regreso a la confiabilidad
161
ligeramente diferentes para la distribucin cada vez, y as los
resultados de confiabilidad ligeramente diferentes. Sin embargo,
empleando lmites de confianza, obtenemos un rango dentro del cual
estn los valores ms probables de confiabilidad a un cierto
porcentaje del tiempo. Esto nos ayuda a calibrar la utilidad de los
datos y la exactitud de las estimaciones resultantes. Adems,
siempre es til recordar que cada parmetro es una estimacin del
verdadero parmetro, el cual es desconocido para nosotros. Este
rango de valores crebles se llama intervalo de confianza.
Lmites bilaterales
Limites unilaterales
162
Recprocamente, un lmite unilateral inferior define un punto que un
porcentaje especifico de la poblacin es mayor.
Por ejemplo, si X es un lmite unilateral superior al 95%,
esto implicara que el 95% de la poblacin es menor que X. Si X es
un lmite unilateral inferior al 95%, esto indicara que el 95% de la
poblacin es mayor que X.
163
4.5 LA DISTRIBUCION DE WEIBULL.
donde,
f (T ) 0, T 0 or , 0, 0,
y,
parmetro de escala,
= parmetro de forma (o de inclinacin)
Y parmetro de ubicacin
164
4.1.3 Distribucin Webull de un parmetro
La Media o MTTF
La media de la pdf de Weibull, T , (tambin llamada MTTF o MTBF
por algunos autores) esta da por:
_
1
T * 1
(3)
1 1
donde 1 es la funcin gamma evaluada al valor de 1 .
La funcin gamma est definida como:
(n) e x x n 1dx
0
Esta funcin esta considerada dentro de Weibull++ para calcular los
valores de (n) a cualquier valor de n.
165
Para el caso de dos parmetros, la ecuacin (3) puede reducirse a:
1
T * 1
Note que algunas personas asumen errneamente que es igual al
MTBF o MHF. Esto es verdad para el caso de =1 desde que
1
1 (2) 1 1.
1
La Mediana
La mediana, T, esta da por:
1
T (ln 2)
(4)
La Moda
La moda, ~
T, se da por:
1
~T 1 1
(5)
La Desviacin estndar
La desviacin estndar, T , est da por:
2
2 1
T * 1 1
166
Funcin Confiabilidad Weibull
La ecuacin para la funcin de densidad cumulativa Weibull
de tres parmetros, cdf, esta da por:
T
F (T ) 1 e
T
R(T ) e
(6)
o:
T t T
R(t T ) e
167
confiabilidad de una nueva misin basada en el hecho que la unidad
o unidades ya acumularon T horas de operacin con xito.
TR * ln R(TR )
1
(7)
1
f (T ) T
(T )
R(T )
168
sustituciones apropiadas para obtener las otras formas, como la
forma de dos parmetros, dnde Y = 0, o de la forma de un
parmetro dnde = C = constante.
Para 0 1 : :
o Como T 0 (o y), f (T )
o Como T , f (T ) 0
169
o f(T) disminuye montonamente y es convexo cuando T aumenta
ms all del valor de y.
o La moda no existe.
Para 1 :
convexo.
170
o Para 1, R(T ) disminuye monotonicamente pero menos aguda
171
Para 1, (T ) , (T) tiene un valor constante de 1 / o:
1
(T )
Esto lo hace conveniente para representar la tasa de falta en el
periodo de vida til de las unidades.
172
4.5.3.2 EFECTOS CARACTERISTICOS DEL PARAMETRO DE
ESCALA, , PARA LA DISTRIBUCION WEIBULL.
173
.0 Y proporciona una estimacin del tiempo para la falla ms
temprano de tales unidades.
o El perodo de vida 0 a +y es un periodo de operacin libre de
fallas de tales unidades.
o El parmetro y puede asumir todo los valores y proporciona una
estimacin del tiempo ms temprano que una falla puede
observarse. Un y negativo puede indicar que las fallas han
ocurrido antes de comenzar la prueba, puede ser durante la
produccin, en el almacenamiento, en trnsito, durante el
chequeo previo al inicio de la misin, o antes del uso real, o y
tiene las mismas unidades que T, como horas, kilmetros,
ciclos, actuaciones, etc.
Ejemplo 1
Asuma que a seis unidades idnticas se les est haciendo
una prueba de confiabilidad a la misma aplicacin y niveles de
tensin de funcionamiento. Todas estas unidades fallan durante la
prueba despus de operar el siguiente nmero de horas, T: 93,
174
34, 16, 120, 53 y 75. Estime los valores de los parmetros para la
distribucin de Weibull de dos parmetros y determine la
confiabilidad de las unidades a las 15 horas de operacin.
i 0.3
MR% * 100
N 0.4
kN (MR) k (1 MR) N k
N
0.50 50%
k i
175
En un papel de probabilidad Weibull, trace los tiempos y
sus correspondientes rangos medios. Una muestra de un papel
de probabilidad Weibull se da en la figura 4-7 y el grfico de los
datos del ejemplo en la figura 4 8.
Q(T ) 1 e 1 e 1 0.632 63.2%)
176
Esto tambin puede obtenerse analticamente de la funcin
confiabilidad Weibull ya que las estimaciones de los dos
parmetros es conocida o:
15 15 1.4
R(t 15) e e 76
90.2%
177
o Repita hasta que los datos ploteados estn en una lnea recta
aceptable.
o El valor de y es el valor del substractor (positivo o negativo) que
pone los puntos en una lnea recta aceptable.
Ejemplo 2
Asuma que se estn realizando un test de prueba de
confiabilidad de seis unidades idnticas bajo las mismas
tensiones y condiciones. Todas las unidades se prueban hasta la
falla y se han registrado los siguientes tiempos de falla: 48, 66, 85,
107, 125 y 152 horas. Encuentre los parmetros de la distribucin
Weibull de tres parmetros usando el mtodo de ploteo de la
probabilidad.
178
4.5.4.2 USO DEL RANGO DE REGRESION EN Y PARA
CALCULAR LOS PARAMETROS DE LA
DISTRIBUCION DE WEIBULL.
T
F (T ) 1 e
(8)
Tomando el logaritmo natural de ambos lados de la ecuacin:
T
ln 1 F (T )
T
ln ln 1 F (T ) ln
o:
ln ln 1 F (T ) ln ( ) ln(T )
(9)
179
y:
b
(12)
Lo que resulta en a ecuacin lineal:
y a bx
yi xi __ __
i 1 i 1
a b y bx
N N
(13)
N N
N xi yi
xi y i i 1
N
i 1
i 1
b 2
N
xi
N
2 i 1
xi N
i 1
(14)
Una vez que a y b son obtenidos, entonces se pueden obtener
fcilmente y de las ecuaciones (11) y (12).
180
El coeficiente de correlacin.
XY
X y
y = desviacin estndar de y.
( xi x) ( yi y)
i 1
N _ N _
( xi x) 2 * ( yi y) 2
i 1 i 1
(15)
Ejemplo 3
181
6 6 6
(ln Ti ) yi ( (ln Ti )( yi ) / 6
i 1 i 1 i 1
b 6 6
(ln Ti ) 2 ( (ln Ti ) 2 / 6
i 1 i 1
__ __ yi ln Ti
i 1
a y b T b i 1
N N
(3.0070 ) 23.9068
a (1,4301) 6.19935
6 6
y la ecuacin (11):
a (6.19935)
e b
= e 1.4301
o:
76,318 hr
0.9956
182
El ejemplo anterior puede resolverse utilizando software. Si desea, la
pdf de Weibull que representan estos datos puede escribirse como:
1 1
T
f (T ) e
o:
1.4302
1
0.4302
1.4302 T
f (T ) e 70317
76 .317 76 .317
183
N N
__ __ xi y2
i 1 i 1
a x b y b
N N
Y
N N
N xi y i
xi yi i 1 i 1
N
i 1
b 2
N
yi
N
i 1
x1 N
2
i 1
donde:
yi ln ln 1 F (Ti
Y:
xi ln (Ti )
y los valores F(Ti) se obtienen de nuevo de los rangos medios.
_ _
Una vez que s obtienen a y b , resuelva la ecuacin (16) para y, la que
corresponde:
a 1
y
x
b b
Resolviendo para los parmetros de las ecuaciones (11) y (12)
conseguimos
a
a
ln( )
b
(17)
1
b
b
(18)
184
El coeficiente de correlacin se evala como antes usando la ecuacin
(15).
Ejemplo 4
Repita el Ejemplo 1 usando la regresin en x.
i 1
__ __ ln Ti yi
i 1
a x b y b i 1
6 6
O:
23.9068 (3.0070 )
a (0.6931 ) 4.318
6 6
185
Por consiguiente, de la ecuacin (18):
1 1
1.4428
0.6931
b
Y de la ecuacin (17)
a 1
*
4.3318
*
1
e b
e 0.6931 1.4428
760811 hr
0.9956
186