Research Article
Validation of a geometric compliance by the measurement of distances in the displacements space
SYMME Laboratory, University of Savoie,
BP 80439,
74944
AnnecyleVieux Cedex, France
^{⁎} Corresponding author: sergent@univsmb.fr
Received:
4
March
2016
Accepted:
9
December
2016
The objective of the study is to show that a oneway machining simulation can be qualitatively satisfactory to validate a range of parts machined on a manufacturing with machinetools 3 axis, even if they must meet geometric tolerances. A reflection on a method of distribution of defects is proposed. A compliance criterion based on the variance domain differences is used. Then, the main component analysis of defects is applied. The proposed model is validated by experimental analysis for which the distributions of defects are given. The results show the major character of the dimensional defect component compared to the angular components (repositioning and machining). A unidirectional approach seems so reliable when the manufactured part and its range of machining remains modest complexity.
Key words: model and measurement 3D and 1D / spread of production / clearance space and deviation space / inertia and center of gravity / noise of measurement and dispersion
© D. Duret and A. Sergent, published by EDP Sciences, 2017
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
At present, CAD tools encourage use at each stage a 3D approach (from the product design stage to final inspection). The 3D approach is relatively well controlled and well integrated, especially in design office [1] and metrology [2]. By cons, during manufacture, process control is complex [3]. In practice, it is mainly through the steering translation adjustment parameters (gauges tools, original program, instructions axes). The driving of setting parameters of rotation (parallelism, perpendicularity …) is fixed by the geometry of the process (machine tools, assembly) [4]. In particular, experimental work of Bui [5] showed that:

Quantification of angular defects requires means of complex measurement and analysis [2,6].

The dispersion of results obtained on the components of rotation is same order than the intrinsic defects obtained in production [7].
In this article, we will try to show that a 3D approach is not necessary in all cases, and that, a 1D approach is often sufficient and avoids a measure of heaviness and treatment. This proposal will be supported by a real experiment on machine tools.
2 Definition of a distance in space travel
2.1 Displacement
The objective is to compare two characteristics of the part − a target characteristic and a realization (a machining part will be called individual).
The characteristic corresponds to the relative position of two elements of references (reference systems related to surfaces or groups of surfaces). This is a displacement of solids (reference systems of associated references):

a translation. T_{O,A→B} (for target) and T_{k,A→B} (for the kth individual);

a rotation R_{O,A→B} (for target) and R_{k,A→B} (for the kth individual).
In the book [8], it is said: it will always express the relative position of two surfaces by displacement said “significant” that superimposes geometric parameters of position of the two surfaces.
If the association of reference systems to surfaces is controlled, we can then compare the two displacements (individual and target). They are defined by the passage of the R1 reference system to R2 [9]. These two coordinate systems are related by the displacement D ={ T_{O1→O2} ; R_{O1→O2} }, which is shown in Figure 1.
The position of a geometric vector can be expressed in two reference systems (Eq. (1)): (1)
If we called R matrix (Eq. (2)): (2)
We obtain the following equation: (3)
Particular case where A is equal to O1 (Eq. (4)): (4)
This transformation preserves the lengths and angles. It is an isometry.
Fig. 1
Systems of coordinate. 
2.2 Identification of moving target (definition)
In the case of two parallel planes, the choice of the coordinate system is quite arbitrary. There are three degrees of freedom for this choice: two translations in the plan and one orientation. The “Good mechanical meaning” restricts this choice. There are still many opportunities.
2.2.1 Didactics example
In this example (Fig. 2), we will confine ourselves to the position of the upper surface to the basic plane A.
The choice of R1 is pretty easy on this part because we have three references (surfaces A, B and C). The choice is unique if we take the framed sides of “40” to fix the origin. Edge effect of B and C surfaces is very important [10]. The unique definition of a coordinate system is often implicit. This may generate an error of interpretation in the declaration of conformity.
For the particular case of target example (Fig. 3), we have the following equation: (5)
Fig. 2
Drawing of definition of the used part. 
Fig. 3
Presentation of target part. 
2.3 Identification of displacement on an individual
To compare the two movements (target and individual), they must be defined in the same way. This is more difficult [7,8]. We must first define R1 on the part [using the Coordinate Measuring Machine (CMM)]. R1 is associated with the part by numerical methods (least squares type).
For the particular case of individual part (Fig. 4), we have the following equation: (7)
Fig. 4
Presentation of individual part. 
3 Differences between two geometric movements
We will make the assumption that the two movements are known by their components p. It can therefore be represented by a point in a space of the same size p (space of individuals).
3.1 Small displacement torsor
The distance between the two displacements (individual and target) is often modeled by a torsor of small displacements (Fig. 5). It is assumed that the two reference systems are aligned (R1 target and R1 individual). Geometrically, we pass from the R2 target to the R2 individual by using a torsor of small displacement (linearization to first order).
Fig. 5
Illustration of the field of geometric deviations. 
3.1.1 Educational example
For example, it holds three parameters to know the difference, namely:

w: translation following axis Oz_{1};

α: rotation about axis Ox_{1};

β: rotation about axis Oy_{1}.
The other three parameters are identically zero. A graphic illustration in three dimensions is possible.
4 Reflections on the distribution of the distance between two displacements
The objective of the analysis of displacements between “target and individual” is to determine all components of displacements (individual target) and show that some are very small compared to others. This will validate a 1D approach is generally sufficient and much simpler implementation.
4.1 Formatting data
We have n manufactured individuals. Each individual is characterized by components p (p ≤ 6).
The measurement protocol is as follows:

the measuring coordinate system is constructed from the three planes (A–B–C);

the association criterion of least squares is used to build the associated numerical element;

moving components machined upper plane (h = 30) (individuals) are calculated at the centroid of the element;

the measurement protocol of each part is defined in Table 1;

these data are grouped in an array X:

An individual e_{i} is a vector at p components (it corresponds to a row in the array). A variable x^{j} (column) is identified by the n measured values for the same parameter.

For our example, we get the following equation:
4.2 The space of individuals (displacements)
4.2.1 Center of gravity
Displacement (line) is defined by its p coordinates. It is an element of the space of individuals. The set of n individuals form a cloud of points with its center of gravity g: (11)
It will be assumed that every individual has the same weight. The center of gravity is an indicator of the rightness of the means of production or measure.
4.2.2 Choice of a metric in the space of individuals
The components of the displacement have different units (angles and lengths). The general formula of distance is of the form: (13)
The most commonly used metric is (14) where s is the experimental standard deviation.
This is equivalent to multiplying each component by 1/s (this allows to be dimension less but with the drawback of s_{i} close to zero).
For example, the distance reduced between the 1 and 2 is equal to the following equation: (16)
4.2.3 Inertia
Inertia in g is equal to the following equation: (17)
We consider that every individual has a weight p_{i}.
Either the inertia in O_{2cible} (18)
The relationship of Huygens gives the following equation: (19)
This scalar indicator brings together the global dispersion of the gap field and its position relative to the target. Unfortunately, it is rather difficult to interpret in the case of physical parameters of a different nature. For our example, we get with (D = I) and (p_{i} = 1): (20) with D = M, inertia is defined as the sum of squared distances reduced: (21)
This definition has the advantage of use values dimensionless.
5 Definition of compliance
In the conventional control frame, the p variables of the individual are defined by limits. Very often, these terminals are linked (not independent). For example, for the location of the upper surface, the link between the local variance α, β and w can be represented by a field of permissible differences [4,9].
The shape of the area deviations permitted depends strongly on the point of calculation of small displacements torsor (expressed in point O_{2cible} in Fig. 6) [11].
The series of parts will be consistent, if the domain (or ball) measured deviations is included in the permissible differences field.
For our didactics example, variables angular α and β are not a problem. However, w represents the projection of the vector on the axis Oz_{2cible} (or Oz_{1cible} for this particular case). Referring to Section 2.1, we can write: (22)
In our case, we go from R1 to R2 by a simple translation, R is the identity I matrix. (23)
In doing so the scalar product with z_{2cible}, yields well w: (24)
If the CMM, calculates the standard deviation w following this procedure, the location of the point O_{2individu} is not fundamental.
Fig. 6
The shape of the area deviations permitted. 
6 The space of variables (components of displacement)
Each variable x^{j} is a list of n numeric values. It can be represented by a vector of R^{n}.
The “vicinity” of variables will be given by the metric defined by the diagonal matrix of weights [12]. (25) with for example, pi = 1/n, so each individual has the same weight.
The scalar product of two variables x^{j} and x^{k} worth: (26)
If the variables are centered, the scalar product is the covariance s_{jk}. The magnitude of a variable is given by . The angle θ_{jk} between twocentered variables is given by the following equation: (27)
The cosine is the linear correlation coefficient.
6.1 Covariance matrix
With weights equal to 1/n, we obtain: (28)
6.2 Linear correlation coefficients
6.3 Principal components analysis
In the most general case, we saw that we could have six settings by individual (three in this example). The principal component analysis allows us to have a representation of individuals by their projection in a lower dimensional, that with the minimum of distortion of the distances between individuals.
In our example, you can search if the two angular axes α and β can be replaced by a single axis (corresponding to that of greatest slope where assembly would result in a defect of systematic orientation). This is not the case for this experiment.
7 Experimentation
The experiment is carried out following 3 steps (Tab. 1):

–
Repeat a measurement of the same part with disassembly and reassembly on the fixture between each measurement (23 times). That allows to build a cloud of points, which shows the influence of the repositioning defect on the measuring machine.

–
Measurement of a sample of parts (23 parts). The parts were products over a long period of production. A complete setting of the machine was made for each part (gauges). The series of 23 parts was machining by students. Measuring gauge tools is therefore not strictly controlled. This can highlight the sources of production defaults.

–
The noise of measurement and the dispersion of recovery has been characterized by the following:

–
The spread of production has been characterized by the following:
8 Conclusion
At present, production simulations favor a 3D approach [6,7]. But in a typical production system, angular instructions are rarely changed. Indeed, they are fit for mounting or machine tool (geometry). A good quality of the product fixation associated with very good geometry of machine tools take allows to control angular defects.
To obtain the geometrical conformity of a product, we need to verify that the displacement components defined previously remain below a threshold.
Otherwise, we can intervene only on some parameters of the means of production (origins programs, gauges …).
This means that other shortcomings noted only depend on the intrinsic quality of the means of production.
In this study, we are able to quantify the components of distance that cannot be changed by the control of the means of production.
The inertia indicators defined above enable to compare and classify different means of production.
From our results, we believe that a 1D simulation is often sufficient to validate the feasibility of a range of production. Figure 7 shows that the main defect found on the batch of manufactured parts corresponds to the component along the Zdirection. The components of angular defect are much lower [13,14].
This is identical to the fixture of the part of Figure 8. Main machining defect correction can be directly accomplished on machine tools using slight modification (gauges, correctives, origin of the frame). Angular defect correction requires significant changes (change of fixture or machining).
In summary, a simulation of production validated in a 1D approach has a high probability to see actual production to meet customer requirements.
Illustration of protocol measurement.
Components of center of gravity of noise of measurement.
Components of inertia of noise of measurement.
Components of center of gravity of spread of production.
Components of inertia of spread of production.
Fig. 7
Spread of the point cloud of defects in production. 
Fig. 8
Spread of the point cloud of noise of measurement and dispersion of recovery. 
References
 A. Desrochers, A. Rivière, A matrix approach to the representation of tolerance zones and clearances, Int. J. Adv. Manuf. Technol. 13, 630 (1997) (In the text)
 P. Bourdet, L. Mathieu, C. Lartigue, A. Ballu, The concept of the small displacement torsor in metrology, Ser. Adv. Math. Appl. Sci. 40, 110 (1996) (In the text)
 M. Kamalinejad, F. Vignat, F. Villeneuve, Simulation of the geometrical defects of manufacturing, Int. J. Adv. Manuf. Technol. 45, 631 (2009) (In the text)
 S. Tichadou, M.K. Nejad, F. Vignat, O. Legoff, 3D manufacturing dispersions: two experimental applications, in 10th CIRP International Seminar on Computer Aided Tolerancing, Paris (2007) (In the text)
 M.H. Bui, A. Sergent, F. Villeneuve, Model of manufactured part and Monte Carlo simulation with experimental data, in 12th CIRP Conference on Computer Aided Tolerancing Manufacturing Tolerance Analysis (2010), Vol. 22 (In the text)
 B.F. Alistair, Least squares bestfit geometric elements, NPL Report DITC Affiliation (Nationale Physical Laboratory, 1989) (In the text)
 M.H. Bui, Identification et simulation des incertitudes de fabrication, PhD from the University of Grenoble, 2011 (In the text)
 A. Rivière, A. Clément, M. Temmerman, Cotation Tridimensionnelle des Systèmes Mécaniques (Pyc livres, Paris, 1994) (In the text)
 A. Clément, Une science nommée technologie, Technol. Form. 108, 32 (2003) (In the text)
 A. Duret, Capabilité et pilotage d'un processus d'usinage, Technol. Form. 115, 31 (2004) (In the text)
 A. Giordano, Clearance space and deviation space. Application to threedimensional chain of dimensions and positions, in Proc. of the 3rd CIRP Seminars on ComputerAided Tolerancing, 1993 (1993) (In the text)
 G. Saporta, Probabilités, analyses des données et statistiques (Editions TECHNIP, 2006) (In the text)
 A. Sergent, Model and optimisation of a machining process thanks to the concept of elementary geometric minimum datum (EGMD), in 14th International CIRP, 2004 (2004) (In the text)
 M. Hien Bui, F. Villeneuve, A. Sergent, Manufacturing tolerance analysis based on the model of manufactured part and experimental data, Proc. Inst. Mech. Eng. B: J. Eng. Manuf. 5, 690 (2012) (In the text)
Cite this article as: Daniel Duret, Alain Sergent, Validation of a geometric compliance by the measurement of distances in the displacements space, Int. J. Metrol. Qual. Eng. 8, 12 (2017)
All Tables
All Figures
Fig. 1
Systems of coordinate. 

In the text 
Fig. 2
Drawing of definition of the used part. 

In the text 
Fig. 3
Presentation of target part. 

In the text 
Fig. 4
Presentation of individual part. 

In the text 
Fig. 5
Illustration of the field of geometric deviations. 

In the text 
Fig. 6
The shape of the area deviations permitted. 

In the text 
Fig. 7
Spread of the point cloud of defects in production. 

In the text 
Fig. 8
Spread of the point cloud of noise of measurement and dispersion of recovery. 

In the text 