METROLOGIA E ENSAIOS

Incerteza de Medição
Prof. Alexandre Pedott
pedott@producao.ufrgs.br

ENG09007 – 2012/2
 Incerteza da Medição

Distribuição de freqüência das

Freqüência de ocorrência

medidas

Erro Sistemático (Tendência)

Erro de Repetitividade
(Aleatoriedade)

Erro Total combinado

Xverdadeiro(VV) Xm x

ENG09007 – 2012/2
 Incerteza da Medição

 A incerteza está relacionada a um valor de medição e não

ao valor verdadeiro do mensurando. O resultado da medição
é apenas a melhor estimativa do valor verdadeiro, que na
ausência de efeitos sistemáticos, é obtido pela média
aritmética de N medições repetidas do mesmo mensurando.

 Caracteriza uma faixa de dispersão ou intervalo e não um

valor pontual. Não deve ser confundida com um “erro.

ENG09007 – 2012/2
 Incerteza da Medição

 É a dúvida remanescente associada ao resultado da

medição. Mede o grau de desconhecimento sobre aquilo
que está sendo medido.

 É um intervalo de valores que podem ser atribuídos

fundamentadamente ao mensurando, de forma
fundamentada e realista, não devendo ser entendida como
uma “faixa de segurança”.

ENG09007 – 2012/2
 Incerteza da Medição

 É a faixa de valores dentro da qual deve se situar o valor

verdadeiro do mensurando.
Associado a um
valor probabilístico
RM = (RB ± IM) unidade

 Resultado base é a estimativa do valor do mensurando que,

acredita-se, mais se aproxime do seu valor verdadeiro.

 Incerteza da medição é o tamanho da faixa simétrica, e

centrada em torno do resultado base, que delimita a faixa onde se
situam as dúvidas associadas à medição.

ENG09007 – 2012/2
 Termos e Definições

 Incerteza de medição: parâmetro associado ao resultado de

uma medição, que caracteriza a dispersão dos valores que
podem ser razoavelmente atribuídos ao mensurando.
 Incerteza padrão: incerteza do resultado de uma medição
expressa como um desvio padrão.
 Incerteza padrão combinada: incerteza padrão do resultado
de uma medição, quando este resultado é obtido por meio dos
valores de várias outras grandezas, sendo igual a raiz
quadrada positiva de uma soma de termos, que constituem as
variâncias e covariâncias destas outras grandezas, ponderadas
de acordo com quanto o resultado da medição varia com
mudanças nestas grandezas.

ENG09007 – 2012/2
 Importância

 Na calibração de equipamentos: instrumentos e padrões,

para verificar se os mesmos encontram-se dentro das
tolerâncias definidas.
 Em ensaios: para verificar se o resultado do ensaio pode
ser aprovado ou não.
 Metrologia legal: para verificar a conformidade de
resultados de medições com limites de tolerâncias legais.
 Na tomada de decisão: reduzir o risco de erro de aprovar
ou rejeitar uma amostra.
 Pode ser um diferencial competitivo entre laboratórios.

ENG09007 – 2012/2
Expressão da Incerteza de Medição

 ISO GUM – Guide to Expression of

Uncertainty in Measurement
 Monte Carlo – Simulação
 Integração Numérica

ENG09007 – 2012/2
 Seminários

Grupo Tema Fatores / Medidas

1 Ensaio Climático (Temperatura, umidade)

2 Calibração Relógio Comparador

3 Calibração Vibração

4 Ensaio Resistência Mecânica

Escolher uma característica de qualidade de um produto.

Apresentação dia: 24/10/2012

ENG09007 – 2012/2
 Expressão da Incerteza de Medição – ISO GUM
1 Definir o modelo
matemático da medição

2 Identificar as componentes
da Incerteza

3 Estimar as incertezas
padrão
Calcular os coeficientes
4 de sensibilidade
Avaliar a existência de
5 correlação
Calcular a incerteza
6 combinada

7 Calcular a incerteza
expandida

8 Corrigir o resultado da
medição
ENG09007 – 2012/2
 Planilha de Incerteza

Estimativa Distribuição de Incerteza Coeficiente de

Grandeza GDL
(média) Probabilidade Padrão Sensibilidade

X1 x1 u(x1) C1 n1

X2 x2 u(x2) C2 n2

XN xN u(xN) C3 nN

Incerteza Combinada Fator de Abrangência Incerteza Expandida

uC(y) kP U

ENG09007 – 2012/2
 Modelo Matemático
1 y = f(x1, x2, ..., xN)
Y é a grandeza do mensurando
Xi são as grandezas que influenciam a medição de Y.

Modelo para calibração: E = I – VV

A incerteza está associada ao mensurando (I) e ao
padrão (VV).

Modelo para ensaio: s = f(F, A) = F / A

A incerteza está associada a medição da força F e
da área da seção transversal A.
ENG09007 – 2012/2
 Componentes de Incerteza
2 São incertezas associadas a:
 Condições ambientais;
 avaliador;
 Equipamentos;
 padrões;
 Método de medição;
 amostragem.

As componentes devem aparecer no modelo.

Diagrama de Causa e Efeito das componentes.

ENG09007 – 2012/2
 Componentes de Incerteza
2 EM
Padrão de
referência
- Desvio padrão da - IM na calibração do BP
série de medições
- Efeito da interpolação - Erro do BP
- Efeito da Resolução
- Resolução
Fontes de incerteza
na calibração de
micrômetro analógico
Paralelismo -
- Afastamento da temperatura
de referência

- Diferença de temperatura
entre BP e EM
Efeitos Temperatura
geométricos

ENG09007 – 2012/2
 Fontes de Incerteza
2
Desvio padrão da média de medições repetidas.

Incerteza da calibração do padrão – incerteza herdada.

Deriva (estabilidade) dos padrões: a variação do padrão no intervalo entre
suas duas últimas calibrações.
Incerteza do fator de correção para erros sistemáticos do padrão.

Erro máximo admissível para o padrão.

Efeito da resolução finita de leitura dos equipamentos.

Efeito das condições ambientais (desvio de temperatura).

Histerese: a indicação de alguns equipamentos pode variar quando as
medições são realizadas no sentido ascendente ou descendente – média do
avanço – média do retorno.

ENG09007 – 2012/2
 Estimativas da Incerteza Padrão
3 O cálculo da incerteza segue a “Lei da propagação das
Incertezas”. A incerteza de cada componente do ensaio ou
calibração é propagada para a obtenção da incerteza final da
medição através de uma expansão de Série de Taylor, truncada
geralmente em primeira ordem.
Supõe-se, para a aplicabilidade do Teorema do Limite
Central, que a medição e a sua respectiva incerteza possuem
uma função densidade de probabilidade (FDP) normal, ou t-
Student para pequenas amostras.
Pressupõem-se observações idêntica e independentemente
distribuídas, com médias e variâncias constantes.
A estabilidade do sistema de medição deve ser assegurada.

ENG09007 – 2012/2
 Estimativas da Incerteza Padrão
3 As incertezas associadas a cada componente devem ser
determinadas.
Tipo A – A incerteza padrão é expressa como um desvio
padrão da distribuição de valores medidos.

Tipo B – A incerteza padrão é expressa como um desvio

padrão, a partir de outro meios de obtenção:
 Experiência ou conhecimento geral do comportamento e
propriedades de materiais relevantes e instrumentos;
 Dados fornecidos em certificados de calibrações e outros
certificados;
 Dados de medições anteriores;
 Valores aceitos de constantes associadas com materiais e
grandezas;
 Especificações de fabricantes.

ENG09007 – 2012/2
 Estimativas da Incerteza Padrão
3
sxi  A incerteza padrão depende
uxi   da distribuição de dados da
Divisor componente de variação.

Tipo de Distribuição de
Componente de incerteza Divisor GDL
avaliação probabilidade
Desvio padrão experimental t-student n n-1
Tipo A Desvio padrão de R&R Normal 1 n-1
Desvio padrão de CC Normal 1 n-1
Erros Retangular 3 infinitos
Resolução Retangular 2 3 infinitos
Tipo B Incertezas herdadas t-student k n eff
Outras Triangular 6 infinitos
Outras Bimodal 2 infinitos

ENG09007 – 2012/2
 Estimativas da Incerteza Padrão
3

ENG09007 – 2012/2 Fonte: Rede Metrológica RS

 Estimativas da Incerteza Padrão
3 Estimativa obtida a partir da especificação, manuais e
catálogos do fabricante, ou do certificado de calibração.

O certificado de calibração de um padrão de aço inoxidável,

de massa m = 1000,000325g, indica uma incerteza de
240mg, para um nível de confiança com k = 3.

A incerteza padrão da massa padrão é dada por:

(Distribuição Normal)

240mg
u (m)   80mg
3

ENG09007 – 2012/2
 Estimativas da Incerteza Padrão
3 Certificado de calibração de um resistor padrão RS
Valor nominal: 10W
Valor estimado: 10,000742W (23ºC)
Incerteza de medição: 129mW a nível de confiança de 99%

A incerteza padrão da massa padrão é dada por:

129mW 129mW
u ( RS )    50mW
k 2,58

ENG09007 – 2012/2
 Estimativas da Incerteza Padrão
3 Estimativa feita a partir da amplitude de variação: limites
inferior e superior – limites simétricos (-a, +a).

-a 0 +a
O manual do fabricante estabelece que o valor do coeficiente
linear de expansão térmica de um bloco padrão de aço é

O erro máximo é: Distribuição Retangular

2 106 C 1
u ( S )   1,2 106 C 1
3
ENG09007 – 2012/2
 Componentes de Incerteza - Exercício
Exercício:
Considere a calibração de um voltímetro digital por comparação.
Código: VO341 / Série: 006-C.
Descrição: voltímetro digital.
Menor Div: 1mV
Unidade: mV
Faixa de leitura: 0 a 200mV
Temperatura ambiente: 20 ±3ºC.

Padrão de referência: multímetro digital HP - 3458 .

Resolução: 0,01mV
Incerteza expandida U = ± 0,001% com k = 2 (em relação ao fundo de escala).
Drift (instabilidade) = ± 0,002mV
Faixa de leitura: 0 a 200mV .
Resultados: A bancada foi ajustada com o voltímetro digital (a ser calibrado) e as
leituras foram realizadas com o padrão.
Tabela de valores medidos dada.
ENG09007 – 2012/2
 Componentes de Incerteza - Calibração

Exercício:
Ponto
Fonte Estimativa Tipo Distribuição Divisor Incerteza GL

ENG09007 – 2012/2
 Coeficientes de Sensibilidade
4 Os coeficientes de sensibilidade são fatores de conversão de
unidades de medida.
y
Ci 
xi
Regra 1: Se o modelo matemático da medição for uma soma de
N variáveis, os coeficientes de sensibilidade serão todos iguais a
um. Caso haja subtração, o sinal do coeficiente será negativo.

Exemplo:

Se o modelo é dado por E = I – VV,

então: CI = 1 e CVV = -1

ENG09007 – 2012/2
 Coeficientes de Sensibilidade
4
Regra 2: Se o modelo matemático da medição for apenas um
produto de N variáveis, os coeficientes de sensibilidade serão
iguais a ( y / xi ).
Note que aqui não é utilizado o valor de incerteza da variável,
mas sim, a melhor estimativa do valor da variável em si.
A divisão segue a mesma regra.

Exemplo:

F s 1 s 1
s CF   CA   F 2
A F A A A

ENG09007 – 2012/2
 Correlações entre as Componentes
5
A correlação existe quando duas grandezas de entrada, Xi e Xj, apresentam
uma relação de dependência entre elas ou com uma terceira grandeza de
entrada comum a ambas.
Exemplo: quando duas grandezas de entrada são medidas com um mesmo
equipamento.
O coeficiente de correlação, r(xi,xj), mede o grau de correlação linear entre
duas variáveis (varia de -1 a 1).
No exemplo dado, r(xi,xj) = 1.

s( xi , x j )
r ( xi , x j ) 
s( xi ) s( x j )
ENG09007 – 2012/2
 Correlações entre as Componentes
5
Efeitos de correlação podem reduzir a incerteza combinada,
se r(xi,xj) < 0.
Efeitos de correlação podem aumentar a incerteza combinada,
se r(xi,xj) > 0.

s( xi , x j )
r ( xi , x j ) 
s( xi ) s( x j )

Onde s(xi) é o desvio padrão associado a xi, s(xj) é o desvio

padrão associado a xj e s(xi,xj) é a covariância associada a xi e xj,
obtidos através dos dados das n medições repetidas de xi e xj.

ENG09007 – 2012/2
 Exemplo

Calibração de um voltímetro digital (Por comparação)

Descrição do voltímetro a calibrar
- Resolução: 0,1mV
- Faixa de leitura: 0 a 100mV
Descrição do padrão de referência
• Incerteza Expandida: ±0,001% (k=2) - Certificado
• Resolução: 0,01mV
• Drift (instabilidade): 0,002mV - Certificado
• Faixa de Leitura: 0 a 100mV

ENG09007 – 2012/2
 Valores Medidos

Leituras
VVC
(mV) 1 2 3 4 5

25 25,09 25,11 25,08 25,07 25,09

50 50,07 50,05 50,06 50,04 50,06

75 75,11 75,12 75,1 75,11 75,13

100 100,15 100,14 100,17 100,16 100,14

ENG09007 – 2012/2
 Fontes de Incerteza
1. Modelo Matemático
Desvio (d) = leitura do Voltímetro – leitura do Padrão
2. Fontes de
Grandeza Tipo
Incerteza
Repetitividade (Δt) A
Resolução do Voltímetro B
Resolução do Padrão B
Incerteza do Padrão B
Drift do Padrão B

d  t  Res(pad)  Res(volt)  U(pad)  Drift

ENG09007 – 2012/2
 Incertezas Padrão

Grandeza Tipo Distribuição Divisor

Repetitividade (Δt) A Normal 1

Resolução do Voltímetro B Retangular 2 3

Resolução do Padrão B Retangular 2 3
Incerteza do Padrão B t-student k=2

Drift do Padrão B Retangular 2 3

ENG09007 – 2012/2
 Estimativas da Incerteza Padrão

sxi  A incerteza padrão depende

uxi   da distribuição de dados da
Divisor componente de variação.

Tipo de Distribuição de
Componente de incerteza Divisor GDL
avaliação probabilidade
Desvio padrão experimental t-student n n-1
Tipo A Desvio padrão de R&R Normal 1 n-1
Desvio padrão de CC Normal 1 n-1
Erros Retangular 3 infinitos
Resolução Retangular 2 3 infinitos
Tipo B Incertezas herdadas t-student k n eff
Outras Triangular 6 infinitos
Outras Bimodal 2 infinitos

ENG09007 – 2012/2
 Incertezas Padrão

Desvio padrão de R&R

s
(Repetitividade) SX 
≠ n
Desvio padrão experimental (S).

Considerando o valor de referência para 50mV

s 0,0114
t    0,0051
n 5
t
u t   0,0051
1
ENG09007 – 2012/2
 Incertezas Padrão

Grandeza Distrib Divisor Incerteza Padrão

Repetitividade t
(Δt) (A)
Normal 1 u t   0,0051
1
Resolução do 0,1
Retangular 2 3 uRe sV   0,0289
Voltímetro (B) 2 3
Resolução do 0,01
Retangular 2 3 uRe sP   0,0029
Padrão (B) 2 3
Incerteza do U 100  0,001%
t-student k=2   0,0005
Padrão (B) k 2
Drift do Padrão 2  0,002
Retangular 2 3 uDrift   0,0011
(B) 2 3
ENG09007 – 2012/2
 Coeficientes de Sensibilidade

d  t  Res(pad)  Res(volt)  U(pad)  Drift

d
C t  1
 t
Grandeza Coeficientes
Repetitividade (Δt) 1
Resolução do Voltímetro 1
Resolução do Padrão 1
Incerteza do Padrão 1
Drift do Padrão 1
ENG09007 – 2012/2
 Correlações entre as Componentes
5
No exemplo dado as variáveis são consideradas como
independentes.
Então o coeficiente de correlação é nulo.

s( xi , x j )
r ( xi , x j )  0
s( xi ) s( x j )

s( xi , x j ) 
1 N
 X i  X Yi  Y 
N  1 i 1

ENG09007 – 2012/2
 Incerteza Combinada
6
É necessário combinar as componentes de incerteza da
mesma maneira a fim de prover um único valor de incerteza.
Quando não houver correlação entre as grandezas de
entrada, a incerteza combinada é dada pela multiplicação da
estimativa de cada incerteza padrão pelo seu respectivo
coeficiente de sensibilidade.
2
n
 f  2
u ( y )     u ( xi )
2

i 1  xi 
c

Z Z Z
Z  ( X 1 )  (
2
X 2 )  ...  (
2
X n ) 2

X 1 X 2 X n
ENG09007 – 2012/2
 Incerteza Combinada
6
Quando houver correlação entre as grandezas de entrada, a
incerteza combinada é dada por:

N 1 N
u   C u xi   2  CiC j uxi ux j r xi , x j 
N
2 2 2
C i
i 1 i 1 j i 1

ENG09007 – 2012/2
 Incerteza Combinada
6
No exemplo:

uC  0,0051  0,0289  0,0029  0,0005  0,0011

2 2 2 2 2

uC  0,00087  0,0295

ENG09007 – 2012/2
 Incerteza de Expandida
7
Através do Teorema do Limite Central, assume-se que a
distribuição de probabilidade da grandeza do mensurando (Y) será
normal, para um grau de liberdade neff.
Pressupostos:
As contribuições das incertezas combinadas tem a mesma ordem
de grandeza.
As distribuições de probabilidade associadas às contribuições de
incerteza se assemelhem da distribuição normal.

ENG09007 – 2012/2
 Incerteza de Expandida
7

O intervalo de ± um desvio
padrão ao redor da estimativa do
mensurando corresponde a uma
probabilidade de abrangência de
aproximadamente 68%.
Para aumentar tal probabilidade
de abrangência, deve-se
multiplicar a incerteza
combinada pelo fator de
abrangência k.

ENG09007 – 2012/2
 Incerteza de Expandida
7 O resultado da combinação é denominado de incerteza
expandida U, dada por:

U  kPuC

Incerteza Combinada Fator de Abrangência Incerteza Expandida

uC(y) kP U

Onde kP é definido para uma determinada probabilidade de

abrangência. Para uma distribuição normal e uma probabilidade
de abrangência de 95,45%, k = 2.

ENG09007 – 2012/2
 Fator de Abrangência
7 Em estudos com poucas repetições (N<30) é recomendável
assumir a distribuição da amostra como t-student.
Neste caso, kP dependerá também do número efetivo de graus
de liberdade neff.
O número de graus de liberdade é diretamente proporcional a
confiabilidade da contribuição para a incerteza.
O Método ISO GUM recomenda o cálculo de neff através da
equação de Welch-Satterthwaite:
4
u ( y)
n eff  N 4 C
u ( xi )
i ni
ENG09007 – 2012/2
 Fator de Abrangência
7 neff k95,45% neff k95,45%
1 13,97 15 2,18

Em geral, neff não será 2 4,53 16 2,17

inteiro. 3 3,31 17 2,16

4 2,87 18 2,15
Deve-se arredondar ou 5 2,65 19 2,14
truncar o número para o 6 2,52 20 2,13
valor inteiro mais 7 2,43 25 2,11
próximo. 8 2,37 30 2,09
Localizar neff na tabela e 9 2,32 35 2,07
determinar kP. 10 2,28 40 2,06
11 2,25 45 2,06
12 2,23 50 2,05
13 2,21 100 2,025
14 2,20 acima de 100 2,00
ENG09007 – 2012/2
 Fator de Abrangência
7 No exemplo do voltímetro, o cálculo de neff será:
Incerteza
Grandeza GDL
Padrão
Repetitividade (Δt) (A) 0,0509 N-1=4

Resolução do Voltímetro (B) 0,0289 ∞

Resolução do Padrão (B) 0,0029 ∞

Incerteza do Padrão (B) 0,0005 ∞

Drift do Padrão (B) 0,0011 ∞

(0,0295) 4
n eff 
(0,0051) 0,0289 0,0029 0,0005 0,0011
4 4 4 4 4
   
4    
ENG09007 – 2012/2
 Fator de Abrangência
7 No exemplo do voltímetro, o cálculo de neff será:
(0,0295) 4
n eff  4
 4471
(0,0051)
4
Localizar neff na tabela e determinar kP. kP  2

As contribuições das incertezas devem ter a mesma ordem

de grandeza.

ENG09007 – 2012/2
 Fator de Abrangência
7 Finalmente, a incerteza expandida será:

U  k PuC  2  0,0295  0,0590

8 O resultado da incerteza expandida deve ser expresso no

máximo com dois algarismos significativos.
O valor do resultado da medição deve ser arredondado
para o mesmo número de casas decimais do valor da
incerteza expandida.

RM = 50,056 ± 0,059mV

ENG09007 – 2012/2
 Planilha de Cálculo da Incerteza da Medição
Unidade
Características do Equipamento da Características do Padrão
Grandeza
Valor
Resolução:
Nominal Certificado U k Veff
Faixa de Indicação: 1

100

Leituras Desvio
Referência Média Tendência
em () Padrão
Leitura 1 Leitura 2 Leitura 3 Leitura 4 Leitura 5 Geral
1
100

Cálculo da Incerteza da Medição

Incerteza
Componentes da Incerteza da Medição Incerteza
k Padrão
Refer. padrão neff
(95,45%) Expandida
Repetitividade Fonte 2 Fonte 3 Fonte 4 ... Combinada
U

100

ENG09007 – 2012/2
 Atividade 2

Calibração de um Ohmímetro digital (Por comparação)

Descrição do ohmímetro a calibrar
- Resolução: 0,01W
- Faixa de leitura: 0 a 10000W
Descrição do padrão de referência
• Incerteza Expandida: ±0,0001% (k=2) - Certificado
• Resolução: 0,001W
• Drift (instabilidade): 0,002W - Certificado
• Faixa de Leitura: 0 a 10000W

ENG09007 – 2012/2
METROLOGIA E ENSAIOS
Incerteza de Medição
Prof. Alexandre Pedott
pedott@producao.ufrgs.br

ENG09007 – 2012/2