Distribuição de Weibull no Excel: forma, escala e tempo até a falha
Resumindo: quanto tempo uma peça dura antes de falhar. Um número define a vida típica; o outro diz se as falhas vêm cedo, ao acaso ou conforme as peças se desgastam.
A distribuição de Weibull descreve por quanto tempo algo dura antes de falhar. O parâmetro de escala define a vida típica e o de forma define como as falhas se distribuem no tempo: abaixo de 1, elas vêm cedo; em 1, vêm a uma taxa constante (a distribuição exponencial); acima de 1, as peças se desgastam. O Excel tem WEIBULL.DIST (DIST.WEIBULL no Excel em português), mas não a inversa, então um sorteio é =scale * (-LN(1 - RAND()))^(1 / shape). Uma bomba com forma 1,5 e escala de 20.000 horas dura 18.055 horas em média, e 25,2% das bombas falham em um ano de funcionamento contínuo.
Forma e escala
A escala (η, eta) é a vida característica: qualquer que seja a forma, 63,2% das unidades já falharam ao chegar a ela. A forma (β, beta) diz como a taxa de falha, a chance de uma unidade que ainda funciona falhar na hora seguinte, muda com a idade:
- Abaixo de 1: a taxa de falha cai com a idade. As unidades fracas falham cedo (mortalidade infantil) e as sobreviventes duram.
- Exatamente 1: a taxa de falha é constante. A idade não importa; é a distribuição exponencial.
- Acima de 1: a taxa de falha aumenta com a idade, conforme as peças se desgastam. Com forma 1,5, uma bomba com 20.000 horas tem 2,0 vezes mais chance de falhar na hora seguinte do que uma com 5.000 horas. Por volta de 3 a 4, a vida é quase simétrica em torno da escala.
Como ler o gráfico: cada linha é um valor de forma, com a mesma escala. Onde a linha está alta, muitas unidades falham naquela idade: com forma 0,8, a maioria das falhas vem logo no início; com forma 3, elas se concentram em torno da vida típica.
A distribuição de Weibull no Excel
WEIBULL.DIST(x, shape, scale, TRUE) dá a proporção de unidades que falharam até a idade x, então =WEIBULL.DIST(8760, 1.5, 20000, TRUE) devolve 25,2%, a chance de uma bomba falhar em um ano (8.760 horas). O Excel chama a forma de alfa e a escala de beta, enquanto engenheiros costumam escrever β para a forma e η para a escala; confira a ordem. Para uma vida aleatória, transforme um número aleatório u entre 0 e 1 em uma idade: a vida é escala × (−ln(1 − u))1/forma, que em uma célula fica =F10*(-LN(1-RAND()))^(1/D10).
Como ler o gráfico: cada barra conta quantos dos 10.000 valores simulados caíram naquele trecho, e a linha é a curva teórica. Quanto mais alta a barra, mais vezes os valores caem ali. As linhas pontilhadas marcam a vida B10 e a média.
| Célula | Fórmula | Valor |
|---|---|---|
| Número aleatório u (I10) | =RAND() | |
| Sorteio (J10) | =F10*(-LN(1-I10))^(1/D10) | |
| Sorteio em uma célula | =F10*(-LN(1-RAND()))^(1/D10) | |
| Média (K10) | =F10*EXP(GAMMALN(1+1/D10)) | 18.055 |
| P10 (L10) | =F10*(-LN(1-0.1))^(1/D10) | 4.462 |
| P90 (M10) | =F10*(-LN(1-0.9))^(1/D10) | 34.874 |
Conferir os números: valores exatos contra 10.000 sorteios do Excel e do xellstorm
| Exato | Excel, RAND() | xellstorm | |
|---|---|---|---|
| Média | 18.055 | 18.229 | 18.055 |
| Desvio padrão | 12.259 | 12.465 | 12.262 |
| P10 | 4.462 | 4.443 | 4.462 |
| P50 (mediana) | 15.664 | 15.628 | 15.665 |
| P90 | 34.874 | 35.730 | 34.870 |
O Excel sorteia com RAND(), totalmente ao acaso; o xellstorm distribui os sorteios de modo uniforme pela faixa (amostragem por Hipercubo Latino), por isso as estatísticas dele ficam mais perto dos valores exatos.
A vida B10, a idade até a qual 10% das unidades já falharam, é o P10: 4.462 horas para esta bomba. A vida mediana é 15.664 horas e a média, 18.055; a média fica acima da mediana porque algumas bombas duram muito mais que a maioria.
A distribuição exponencial: uma Weibull com forma 1
Quando falhas ou chegadas ocorrem ao acaso, a uma taxa média constante, o tempo entre elas é exponencial: uma quebra a cada 500 horas em média equivale a uma taxa de 0,002 por hora. Ela tem um único parâmetro, a média (ou o inverso dela, a taxa), e não tem memória, então o tempo que já se esperou não diz nada sobre o tempo que falta: uma máquina que rodou 1.000 horas sem quebrar tem a mesma chance, 71,5%, de passar pela próxima semana (168 horas) que uma máquina nova. A maioria dos intervalos é curta: 63,2% deles são menores que a média, e a mediana é de apenas 347 horas.
| Célula | Fórmula | Valor |
|---|---|---|
| Número aleatório u (I11) | =RAND() | |
| Sorteio (J11) | =-D11*LN(1-I11) | |
| Sorteio em uma célula | =-D11*LN(1-RAND()) | |
| Média (K11) | =D11 | 500 |
| P10 (L11) | =-D11*LN(1-0.1) | 53 |
| P90 (M11) | =-D11*LN(1-0.9) | 1.151 |
Conferir os números: valores exatos contra 10.000 sorteios do Excel e do xellstorm
| Exato | Excel, RAND() | xellstorm | |
|---|---|---|---|
| Média | 500 | 497 | 500 |
| Desvio padrão | 500 | 499 | 500 |
| P10 | 53 | 54 | 53 |
| P50 (mediana) | 347 | 347 | 347 |
| P90 | 1.151 | 1.132 | 1.151 |
O Excel sorteia com RAND(), totalmente ao acaso; o xellstorm distribui os sorteios de modo uniforme pela faixa (amostragem por Hipercubo Latino), por isso as estatísticas dele ficam mais perto dos valores exatos.
No Excel, EXPON.DIST(x, rate, TRUE) recebe a taxa, 1 / média, e não a média. Um sorteio é =-D11*LN(1-RAND()). O tempo total de vários desses eventos em sequência segue uma distribuição gama, e o número de eventos em um período fixo, uma distribuição de Poisson.
Quando usar Weibull ou exponencial
- Weibull para a vida de componentes e equipamentos, fadiga e qualquer tempo até a falha em que a idade importa. Com dados de falha, ajuste a distribuição. O exemplo de confiabilidade de um sistema em série combina quatro componentes Weibull, entre eles esta bomba.
- Exponencial para tempos entre eventos aleatórios: quebras de um sistema com muitas peças, chegadas de clientes, chamadas. Não a use para desgaste, em que ela subestima as falhas de unidades antigas.
- Lognormal para tempos de reparo e outras durações que raramente são curtas; costuma ser ajustada junto com a Weibull aos mesmos dados.
Como usar as distribuições de Weibull e exponencial no xellstorm
- Transforme a célula em entrada e, na etapa Distribuições, escolha Weibull e informe forma e scale; loc é um período sem falhas, em geral 0. Para uma exponencial, escolha Exponencial e informe a média como scale.
- Com tempos de falha, selecione a entrada e use “A partir de dados” no painel lateral: ele ajusta Weibull, exponencial, gama e lognormal (todas começando em zero), entre outras formas, e ordena os ajustes, para você ver qual deles os dados sustentam.
- Transforme vidas na resposta de que você precisa com fórmulas da planilha: um sistema que para na primeira falha é o
MINdas vidas dos componentes, eIF(life<8760,1,0)conta as falhas em um ano.
Modelos de suplementos são convertidos a partir de RiskWeibull, VoseWeibull e PsiWeibull, e de RiskExpon, VoseExpon e PsiExponential.
Perguntas
Como gerar números aleatórios de Weibull no Excel?
O Excel não tem WEIBULL.INV, então inverta a distribuição você mesmo: =scale*(-LN(1-RAND()))^(1/shape). Para forma 1,5 e escala 20.000, fica =20000*(-LN(1-RAND()))^(1/1.5). Os sorteios têm exatamente a distribuição de Weibull, como mostra a tabela acima.
O que significam a forma e a escala de uma distribuição de Weibull?
A escala é a vida característica, a idade até a qual 63,2% das unidades já falharam. A forma descreve a taxa de falha ao longo do tempo: abaixo de 1 ela cai (falhas precoces), em 1 é constante (falhas aleatórias) e acima de 1 sobe (desgaste).
O que é a vida B10?
A vida B10 é a idade até a qual 10% das unidades já falharam: o P10 da distribuição da vida, escala × (−ln 0,9)1/forma. Para uma bomba com forma 1,5 e escala de 20.000 horas, são 4.462 horas.
Como a distribuição exponencial se relaciona com a de Weibull?
Uma distribuição exponencial é uma Weibull com forma 1, com a média como escala. Ela descreve vidas que não dependem da idade: a chance de falhar na próxima hora é a mesma para uma unidade nova e para uma antiga.
Veja também
- Guia
Qual distribuição de probabilidade usar? Guia para modelos de Monte Carlo
Encontre a distribuição de probabilidade certa em duas perguntas: PERT, triangular, normal, lognormal, Weibull, Poisson, binomial, beta e outras, com a finalidade de cada uma e a fórmula dela no Excel. - Confiabilidade de equipamentos
Qual a chance de um sistema em série falhar em um ano?
Quatro componentes em série com tempos de vida Weibull: 79,0% das unidades falham em um ano (8.760 h) e a vida média é de 5.712 h. Simulado x exato. - Guia
Distribuição gama no Excel: tempos de espera, totais e GAMMA.INV
A distribuição gama descreve grandezas positivas e assimétricas à direita, como o tempo total de várias etapas. Forma e escala, GAMMA.INV e GAMMA.DIST no Excel, e gama x lognormal.
O xellstorm é uma ferramenta de simulação de Monte Carlo para modelos do Excel, que roda no navegador: sem suplemento, e o arquivo nunca sai do seu computador.