ISSN: 2319-9873

All submissions of the EM system will be redirected to **Online Manuscript Submission System**. Authors are requested to submit articles directly to **Online Manuscript Submission System** of respective journal.

**Indrit Vozga ^{1}^{*} and Jorgaq Kacani^{2}**

Mechanical Engineering, Production and Management Department, Albania

- *Corresponding Author:
- Indrit Vozga

Mechanical Engineering Faculty

Production and Management Department

Albania.

**E-mail:**indritvozga@gmail.com

**Received date:** 10/10/2018; **Accepted date:** 20/12/2018; **Published date:** 26/12/2018

**Visit for more related articles at** Research & Reviews: Journal of Engineering and Technology

This document is intended to describe in detail the technical actions that should be taken by authorized personnel in order to assess the uncertainty of measurement of compressive resistance. We try to fulfill the specific requirements of SSH ISO/IEC 17025 standard 5.4.6 paragraph "Assessment of uncertainty of measurements"

Ceramic, Masonry, Resistance

**Measurement Procedure and Calculation**

Working standard SSH 548/1-87. Determination of strength limit in compression

**Calculation**

Rc = F/S (1)

where F-maximum force applied from compression till failure the S-base surface of the sample that is subjected to this force.

**Identification of Uncertainty Sources**

Since it is not possible to measure the true value of each measurement, it is important to report not only the measurement value but also the uncertainty of this value together with the level of reliability. All possible factual sources that really influence the assessment of measurement uncertainty should be evaluated and continuous measures should be taken to eliminate them as much as possible [1]. Based on the testing method and the real conditions of its implementation, we identify all the components of the uncertainty by making a metrological and statistically valid estimation.

The main factors influencing the measurement process are specified as follows:

Specimen geometry

Specimen dimensions

Planarity of the specimen surface

Angles between faces (plans) that form the sample

• Curing conditions that have to do with temperature and relative humidity

• The conductivity of test mode

• The speed of application of load force of the compressor

• Temperature and air humidity condition during test performance

• Necessary corrections to be made regarding the reading of force transmitted by the testing machine according to calibration modalities

**Measurement Model and Quantitative Assessment of Uncertainty Sources**

Measurement model is called the function: Y=f(x1, x2, x3 ... xn) and x1, x2, x3 ... xn are the input values. Typical input values are those derived from measurement process, those that are reported from calibration certificates of used instruments as well as influenced values that are environmental variables such as air temperature, air humidity, and air pressure.

During this test, we have the uncertainty of category A and category B.

Category A. When value Xi can be estimated directly from the laboratory by repurposing a metering process under controlled conditions, resulting in a series of values.

Statistical theory tells us how to use the collected information. Practically because of cost and time, repeated measurements of a magnitude form a limited set of values. The variable that rationalizes exactly the limited series values is studied by the Student:

Category B. the estimation of uncertainty is based on a series of repeated observations. This kind of evaluation is based on scientific methods

These evaluations are carried out in a manner different from that based on a series of repeated observations. The minimum information situation is presented by an interval defined by two values of ximax and ximin outside which the possibilities are not to be found in the range of the magnitude, while within the interval all values have the same probability called quadratic width equal to ximax - ximin.

In this case, it can be attributed as an estimation of xi the average value of the interval:

Below is given the expression of composed uncertainty (absolute value). In frequent applications for simplicity, we assume the measurement function with no more than 3 input values.

The function Y=A/B in this case we have a ratio, ie composed uncertainty in absolute value:

Below are the uncertainties associated with the flatness of the sample surface and the angles between its faces. From the equation:

Considering factors of influence, the following model of measurement is derived:

Where: F is the value of the maximum force transmitted by the compressor, which corresponds to the load of the sample failure

ΔF1 Correction determined by the compressor setting

ΔF2 Correction determined by the application speed of the load

ΔF3 Correction determined by the conditions of aging

ΔF4 Correction determined by nonlinearity

Ssample base surface over which the fore is applied. This is determined by the equation:

In which ℓ1 and ℓ2 are the long and the width of the sample surface respectively and Δℓ1 and Δℓ2 are the corrections determined by relative uncertainty related to their measurement.

In the following equations, the two terms in denominator and numerator of formulas are considered statistically independent.

**Calculation of Composed Uncertainty**

Absolute composed uncertainty may be derived from the equation:

**Defining the Best Evaluation of F and its Uncertainties**

The distribution of the values of the measurement object, ie the resistance to compression, is evaluated by performing
a series of repeated measurements** Table 1**.

**Table 1**. Distribution of the values of the measurement object.

Test number (n) | Fi (N) |
---|---|

1 | 135592.8 |

2 | 140865.2 |

3 | 125048 |

4 | 135592.8 |

5 | 126948 |

6 | 135593.4 |

7 | 128079 |

8 | 140762.4 |

9 | 138487 |

10 | 126072 |

133304.1 | |

s = u(Fi) | 6178.3 |

U() = | s/1955.2 |

CV = s/ | 0.046 |

v = n-1 | 10-1=9 |

**Determination of ΔF1 Value and Its Uncertainty U(ΔF1)**

Correction obtained from reading the F value is 0.

In the case of evaluating relative uncertainty U(ΔF1) with B method, all values have the same probability called foursquare [2]. (In this case, it is not given the level of confidence). Uncertainty is given for a=1 kN=1000 N. This is the scale division of the testing machine.

(In our case the testing machine is calibrated and we obtain the specific value from the calibration certificate. According to calibration certificate no F-355/2009 with measuring range 0-500 kN and with 1 kN scale division, the expanded uncertainty is U=0.25 kN so the standard uncertainty is u=U/2 for 95% distribution level and k=2 covering factor.

**Determination of ΔF2 Value and its Uncertainty U(ΔF2)**

Value of ΔF2 = 0. As the technical documents highlight is difficult to evaluate quantitatively the effect of load application speed over compression resistance of a sample. For simplicity from technical documents that uncertainty values u(ΔF2) to be 2% of the mean value of F.

**Determination of ΔF3 Value and its Uncertainty U(ΔF3)**

Correction mean value ΔF3 = 0. As the technical documents highlight is difficult to evaluate quantitatively the effect of curing conditions over compression resistance of a sample, it is supposed for simplicity from technical documents that uncertainty values:

**Determination of ΔF4 Value and its Uncertainty U(ΔF4 )**

Correction mean value is zero. As the technical documents highlight is difficult to evaluate quantitatively the effect of the planar condition over compression resistance of a sample, it is supposed for simplicity from technical documents that uncertainty values to be 1.5% of the mean value

**Determination of U(ΔF5) Value and its Uncertainty**

Correction mean value is zero. Test methods specify conditions according to the plane surface angle between faces. From technical documents uncertainty values U(ΔF5) to be 0.1% of F

**Determination of U(ΔF6) Value and its Uncertainty**

Correction mean value is zero. Test methods specify conditions according to the positioning of plates of load appliances in the center, and this uncertainty belongs to operators since it is difficult to evaluate quantitatively this effect it is supposed for simplicity from technical documents that uncertainty values to be 0.5% of F.

**Determining the Surface of the Sample and Its Uncertainty**

The sample surface is determined by the equation:

In this case, we have : so we have obtained the input sizes composed uncertainty in absolute value:

So, in our case we can write:

(17)

where: Standard deviation or (standard uncertainty) The absolute standard deviation of the mean� The degree of freedom number

We evaluate the standard uncertainty U(Δℓ1) and U(Δℓ2) and corrections from readings Δℓ1 and Δℓ2. The corrections
obtained from Δℓ1 and Δℓ2 readings are estimated zero **Table 2.**

**Table 2.** Measurement results for ℓ1 and ℓ2.

n | ℓ1 (mm) | n | ℓ2 (mm) |
---|---|---|---|

1 | 248 | 1 | 248 |

2 | 250 | 2 | 250 |

3 | 249 | 3 | 249 |

4 | 248 | 4 | 249 |

5 | 248 | 5 | 248 |

6 | 248 | 6 | 248 |

7 | 247 | 7 | 247 |

8 | 248 | 8 | 246 |

9 | 248 | 9 | 245 |

10 | 248 | 10 | 246 |

248.3 | 246.3 | ||

S = u() | 0.625 | 0.76 | |

U() | 0.198 | 0.241 | |

v = n - 1 | 10 – 1 = 9 | 10 – 1 = 9 |

Relative uncertainty type U(Δℓ1) and U(Δℓ2) are estimated to have a value of 0.05 mm.

In conclusion, we have:

(19)

(20)

**The Best Estimation of Measurement and Calculation of Composed Uncertainty**

We replace the founded values:

For calculation of composed uncertainty we apply the following equation:

Even though composed type uncertainty u(y) might be enough to characterize a measurement, in many applications we prefer to determine a wider interval U(y) round the result y in order to fall a greater part of values [3].

This broader interval is the extended uncertainty U(y)=k . u(y)

where: k is found in variables tp of Student given in the relevant table.

To choose the value of k in Student table we should decide the level of probability (generally 95%) and calculate the number of effective free degrees that are attributed to u(y). Such a calculation can be made with the equation below:

Weich-Satterthwaite:

From which if

In our case we are calculating Veff with equation 1:

Veff = 10.65

From Student table, we have k=2.18

We find the covering factor k from the table below. This table is based in a t distribution of Student for a credibility level
of 95.45% **Table 3**.

**Table 3.** T distribution of Student for a credibility level of 95.45%.

Veff |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 10 | 20 | 50 | ∞ |

k |
13.97 | 4.53 | 3.31 | 2.87 | 2.65 | 2.52 | 2.43 | 2.37 | 2.28 | 2.13 | 2.05 | 2.00 |

U(y)=k . u(y)

Considering the approaches that are not mentioned, it is considered when Veff=l0 can be substituted the exact value
of tp with a covering coefficient k=2. This is due to the fact that the calculation of uncertainty constitutes a significant
difficulty because the measurement object is influenced by both types A uncertainty and type B uncertainty. So k=2 for
n=10 **(Table 4)**.

**Table 4.** Calculation of uncertainty constitutes.

n |
6 | 8 | 10 | 12 | 16 | 20 | 25 | ≥ 30 |

k |
2.33 | 2.19 | 2.10 | 2.05 | 1.98 | 1.93 | 1.88 | 1.64 |

This mode is used for ease of calculation.

k=2.1 for n=l0 U=2.l . 0.07=0.147=0,15 N

This is the characteristic resistance to compression coupled with its uncertainty **Table 5**.

**Table 5.** Characteristic resistance to compression coupled with its uncertainty.

Xi | Xi value | Uncertainty type u(xi) | Distribution probability | Uncertainty contribution |
---|---|---|---|---|

F | =133304 | u(F)=6178.3 | normal | U(Y)=1955 |

ΔF1 | ΔF1=0 | U(ΔF1)=577 | rectangular | U(ΔF1)=577 |

ΔF2 | ΔF2=0 | U(ΔF2)=2% =2660 |
technical assessor | U(ΔF2)=2666 |

ΔF3 | ΔF3=0 | U(ΔF3)=1.5% =2000 |
technical assessor | U(ΔF3)=2000 |

ΔF4 | U(ΔF4)=1.5% =2000 |
technical assessor | U(ΔF4)=2000 | |

ΔF5 | U(ΔF5)=0.1% =133 |
technical assessor | U(ΔF5)=133 | |

ΔF6 | U(ΔF6)=0.5% =667 |
technical assessor | U(ΔF6)=667 | |

ℓ1 | =248.3 | U(ℓ1)=0.63 | normal | U()=0,198 |

ℓ2 | =246.3 | U(ℓ2)=0.76 | normal | U()=0,241 |

Δℓ1 | Δℓ1=0 | U(Δℓ1)=0.05 | technical assessor | U(Δℓ1)=0.05 |

Δℓ2=0 | U(Δℓ1)=0.05 | technical assessor | U(Δℓ2)=0.05 | |

Δℓ2 | ||||

In percentage the uncertainty is: 0.15/2.18 × 100 = 6.88%

From the values of the table we build the graph in Microsoft Excel and look which has the greatest contribution [4] **Figure 1**. Those that have a very small contribution are neglected. (Uncertainties that are 1/5 of the main uncertainty are
neglected) **Figure 2**.

For the implementation of this procedure, the measurement uncertainty assessment register is used in the SC-Rr management system [5].

The uncertainty of the sampling process is not taken into account and it should be made clear that the result and the uncertainty associated are applied only to the samples being tested and not applied to any party from which the sample may have been taken)

The extended reported uncertainty is based on a standard multiplicity uncertainty with a coverage coefficient k=2.1 which for a margin of Veff =10.7 degrees of effective freedom provides a 95% confidence level of approximately.

- https://ec.europa.eu/jrc/en/eurl/ecvam
- http://www.nab.lt/Files/EA%204-02%20Matavimo%20neapibr%C4%97%C5%BEties%20i%C5%A1rai%C5%A1ka%20kalibravime.pdf
- https://european-accreditation.org/wp-content/uploads/2018/10/ea-4-16-g-rev00-december-2003-rev.pdf
- Curtin WG, et al. Structural masonry designer’s manual, Blackwell, Oxford. 1995.
- Saileysh SS. Bibliography on masonry structures-technical report.