Articulo de referencia

Proceso de renovación generalizada

En la teoría matemática de la probabilidad , un proceso de renovación generalizado (PRG) o proceso de renovación G es un proceso estocástico puntual que se utiliza para modelar ...

En la teoría matemática de la probabilidad , un proceso de renovación generalizado (PRG) o proceso de renovación G es un proceso estocástico puntual que se utiliza para modelar el comportamiento de falla/reparación de sistemas reparables en ingeniería de confiabilidad . El proceso puntual de Poisson es un caso particular de PRG.

Modelo probabilístico

Era virtual

El proceso de G-renovación es introducido por Kijima y Sumita a través de la noción de la era virtual . [ 1 ]

yi=qti{\displaystyle y_{i}=qt_{i}}
dónde:
ti{\displaystyle t_{i}}yyi{\displaystyle y_{i}}es la edad real y virtual (respectivamente) del sistema en/después de la i -ésima reparación,
q{\displaystyle q}es el factor de restauración (también conocido como factor de eficacia de la reparación),
q=0{\displaystyle q=0}, representa la condición de una reparación perfecta, donde la edad del sistema se restablece a cero después de la reparación. Esta condición corresponde al Proceso de Renovación Ordinario .
q=1{\displaystyle q=1}, representa la condición de una reparación mínima, donde la condición del sistema después de la reparación permanece igual que justo antes de la reparación. Esta condición corresponde al proceso de Poisson no homogéneo .
0<q<1{\displaystyle 0<q<1}, representa la condición de una reparación general, donde la condición del sistema está entre reparación perfecta y reparación mínima. Esta condición corresponde al Proceso de Renovación Generalizado .

Kaminskiy y Krivtsov [ 2 ] extendieron los modelos de Kijima permitiendo q  >  1, de modo que la reparación daña (envejece) el sistema en un grado mayor que el que tenía justo antes de la falla correspondiente.

ecuación de renovación G

Matemáticamente, el proceso de renovación G se cuantifica mediante la solución de la ecuación de renovación G:

W(t)=0t(gramo(τ0)+0τw(incógnita)(gramo(τincógnitaincógnita)dincógnita)dτ{\displaystyle W(t)=\int _{0}^{t}(g(\tau \mid 0)+\int _{0}^{\tau }w(x)\cdot (g(\tau -x\mid x)\,dx)\,d\tau }
dónde,
gramo(τincógnita)=(t+qincógnita,θ)1F(qincógnita,θ){\displaystyle g(\tau \mid x)={\frac {(t+qx,\theta )}{1-F(qx,\theta )}}}t,incógnita0{\displaystyle t,x\geq 0}
w(incógnita)=dW(incógnita)dincógnita{\displaystyle w(x)={\frac {dW(x)}{dx}}}
f ( t ) es la función de densidad de probabilidad (PDF) de la distribución subyacente del tiempo de falla,
F ( t ) es la función de distribución acumulativa (CDF) de la distribución del tiempo de falla subyacente,
q es el factor de restauración,
θ{\displaystyle {\theta }}es el vector de parámetros de la distribución subyacente del tiempo de falla.

No es posible obtener una solución analítica para la ecuación de renovación G. Además, las aproximaciones numéricas son difíciles de obtener debido a la serie infinita recurrente. Kaminiskiy y Krivtsov desarrollaron un método basado en el método de Monte Carlo para resolver la ecuación de renovación G. [ 2 ] [ 3 ]

Estimación estadística

El proceso de renovación G solo alcanzó popularidad práctica en la ingeniería de confiabilidad después de que se dispuso de métodos para estimar sus parámetros.

Enfoque de Monte Carlo

La estimación LSQ no lineal del proceso de renovación G fue propuesta por primera vez por Kaminskiy y Krivtsov. [ 2 ] Un tiempo entre llegadas aleatorio de un proceso de renovación G parametrizado viene dado por:

incógnitai=F1(1Ui[1F(qSi1)])qSi1{\displaystyle X_{i}=F^{-1}(1-U_{i}[1-F(qS_{i-1})])-qS_{i-1}}
dónde,
Si1{\displaystyle S_{i-1}}es la edad real acumulada antes del i -ésimo intervalo entre llegadas,
Ui{\displaystyle U_{i}}es una variable aleatoria con distribución uniforme ,
F{\displaystyle F}es la función de distribución acumulada (CDF) de la distribución subyacente del tiempo de falla.

La solución de Monte Carlo fue posteriormente mejorada [ 4 ] e implementada como un recurso web. [ 5 ]

enfoque de máxima verosimilitud

Los procedimientos de máxima verosimilitud fueron posteriormente analizados por Yañez et al. [ 6 ] y Mettas & Zhao [ 7 ] . La estimación del factor de restauración G-renovación fue abordada en detalle por Kahle & Love [ 8 ] .

Método de regularización para la estimación de parámetros GRP

La estimación de los parámetros del proceso de renovación G es un problema inverso mal condicionado y, por lo tanto, la solución puede no ser única y es sensible a los datos de entrada. Krivtsov y Yevkin [ 9 ] [ 10 ] sugirieron primero estimar los parámetros de distribución subyacentes utilizando únicamente el tiempo hasta el primer fallo. Luego, los parámetros obtenidos se utilizan como valores iniciales para el segundo paso, donde todos los parámetros del modelo (incluidos los factores de restauración) se estiman simultáneamente. Este enfoque permite, por un lado, evitar soluciones irrelevantes (máximos o mínimos locales erróneos de la función objetivo) y, por otro lado, mejorar la velocidad de cálculo, ya que el número de iteraciones depende significativamente de los valores iniciales seleccionados.

Limitaciones

Una limitación del Proceso de Renovación Generalizado es que no puede tener en cuenta las reparaciones "mejores que nuevas". [ 11 ] Se ha desarrollado el proceso de renovación G1 , que aplica el factor de restauración al parámetro de vida de una distribución de ubicación-escala para poder tener en cuenta las reparaciones "mejores que nuevas" además de otros tipos de reparación.

Referencias

  1. Kijima, Masaaki; Sumita, Ushio (1986). "Una generalización útil de la teoría de renovación: procesos de conteo gobernados por incrementos markovianos no negativos" . Journal of Applied Probability . 23 (1). Applied Probability Trust: 71– 88. doi : 10.2307/3214117 . JSTOR 3214117. S2CID 222275620 .  
  2. 1 2 3 Kaminskiy, MP; Krivtsov, VV (1998). "Un enfoque de Monte Carlo para el análisis de la fiabilidad de sistemas reparables". Evaluación y gestión probabilística de la seguridad . Londres: Springer–Verlag. págs. 1063–1068 . 
  3. Krivtsov, VV (2000). Modelado y estimación del proceso generalizado de renovación en el análisis de confiabilidad de sistemas reparables (Tesis doctoral). Universidad de Maryland, College Park, ISBN/ISSN: 0599725877.
  4. Yevkin, A. (2011). «Enfoque de Monte Carlo para la evaluación de la disponibilidad y la intensidad de fallos bajo el modelo de proceso de renovación G». En Berenguer, Christophe; Grall, Antoine; Guedes Soares, Carlos (eds.). Avances en seguridad, fiabilidad y gestión de riesgos . Londres: CRC Press. pp. 1015–1020 . doi : 10.1201/b11939 . ISBN  9780429217265.
  5. Yevkin, A. "Calculadora del proceso de renovación G" . Consultado el 13 de mayo de 2021 .
  6. Yañez, M.; Joglar, F.; Modarres, M. (agosto de 2002). "Proceso de renovación generalizado para el análisis de sistemas reparables con experiencia limitada en fallas" . Reliability Engineering & System Safety . 77 (2): 167– 180. doi : 10.1016/S0951-8320(02)00044-3 .
  7. Mettas, A.; Zhao, W. (24 de enero de 2005). Modelado y análisis de sistemas reparables con reparación general . Simposio anual sobre fiabilidad y mantenibilidad 2005. Alexandria, VA.
  8. Kahle, W.; Love, C. (2003). «Modelado de la influencia de las acciones de mantenimiento». Métodos matemáticos y estadísticos en fiabilidad . Serie sobre calidad, fiabilidad y estadística de ingeniería. Vol. 7. págs. 387–399 . doi : 10.1142/9789812795250_0025 . ISBN   978-981-238-321-1.
  9. Krivtsov, VV; Yevkin, O. (julio de 2013). "Estimación de los parámetros del proceso de renovación G como un problema inverso mal condicionado" . Reliability Engineering & System Safety . 115 : 10–18 . doi : 10.1016/j.ress.2013.02.005 .
  10. Krivtsov, Vasiliy; Yevkin, Alex (2017). Técnicas de regularización para la predicción de fallas recurrentes bajo modelos de Kijima . Simposio anual de confiabilidad y mantenibilidad 2017. Orlando, FL.
  11. ^ Kaminskiy, diputado; Krivtsov, VV (junio de 2010). "Proceso de renovación G1 como modelo de sistema reparable". arXiv : 1006.3718 [ estad.ME ].