Research on the Inﬂuence of Production Technologies on the Positioning Accuracy of a Robotic Arm for Low-Handling Weights

: The subject of the paper is the research of production technologies’ inﬂuence on positioning accuracy of a robotic arm. The aim was to ﬁnd out whether different production technologies (additive and conventional) and the related design differences of the robotic arm affect its operational functionality. In the research, positioning accuracy of a robotic arm formed by three partial arms was speciﬁcally investigated, while the ﬁrst partial arm, Arm I, was manufactured by two different technologies. On the robotic arm, the research was carried out in such a way that the ﬁrst partial arm, Arm I, was being continuously changed and was available for research purposes in two variants. Each of the Arm I variants was manufactured using a different technology (additive and conventional) while, at the same time, the individual variants also differed in construction. The design differences of both variants were related to the production technology used. The measurement of positioning accuracy was performed with the use of two methods. Speciﬁcally, a contact and a non-contact method were used. The contact method was implemented on a 3D-measuring machine, RAPID, and the second contactless method was performed using an inductive sensor.


Introduction
The use of robots for handling and positioning, welding, subtractive and additive production has been expanding in recent years in line with the Industry 4.0 concept [1]. A robotic arm can be a separate mechanism or part of a more complex robot. It is a type of mechanical hand that is usually programmable and has similar functions to a human hand [2].
At present, we encounter the use of various production technologies in the production of robotic arm components. CNC machining technologies are the most widespread [3]. Covaciu and Filip [4] designed and validated a robotic arm made with two CNC machines. Gordaninejad and Vaidyaraman are among the first researchers to compare the positioning accuracy of robotic arms made of conventional metallic materials and advanced composites [5].
Until recently, it was not possible to find a published use of an additive manufacturing component in an industrial robotic arm. However, the great expansion of 3D printing technology is creating increasingly suitable conditions for this [6]. A conventional 3D printer uses layering in a horizontal plane to produce a 3D printed part. Ishak and Larochelle [7] proposed the integration of existing additive manufacturing process technologies with the arm of an industrial robot to create a multilayer layered 3D printer. Their approach allows printing to be layered in multiple planes, while existing conventional 3D printers 2 of 23 are limited to one plane of the toolpath. The integration of the robotic arm and the extruder allows the movement of several planes of the tool path to be used in the production of structural parts.
Hajash et al. [8] introduced a fast liquid printing device that can freely print in any direction, rather than layer by layer, depositing liquid material in a granular gel to create 3D structures. Recently, however, several scientific papers have appeared on the additive manufacturing of robotic arm components. Mick et al. [9] proposed a prototype of a robotic arm made using 3D printing, which is also economically more advantageous compared to the price of an industrial robot thanks to conventional drives. Junia Santillo Costa et al. [10] implemented and validated a 3D printed open source robotic arm with 6 degrees of freedom made of ABS (Acrylonitrile Butadiene Styrene) material, due to higher mechanical strength. Ismail et al. [11] designed and developed a robotic arm for lifting light parts with 4 degrees of freedom. They used a 3D printing method to make the robotic arm components, which provided more accurate dimensions and time and cost savings. Wang et al. [12] proposed a special holder for the InnoMotion robotic arm, where the main components of the holder were made by 3D printing using plastic material and are fully compatible with the MR and CT system and the robotic arm.
The positioning accuracy of the robotic arm can vary widely in the workspace, which is influenced by various factors. Therefore, researchers propose various methods and methodologies to improve reliability and repeatability of positioning accuracy [13]. One of the parameters influencing the accuracy of the robotic arm is vibration. Elvira-Ortiz et al. [14] proposed a methodology to improve the estimation of kinematic parameters on industrial robots by correctly suppressing the vibrational components present on the signals obtained from the two primary sensors: the accelerometer and the gyroscope. Their results prove that the sensor fusion technique, accompanied by correct vibration suppression, provides a better estimate of the kinematic parameters than other proposed techniques. The accuracy of the robotic arm is also affected by sensing the position of each joint with a high-resolution optical coding device that cannot detect certain mechanical deformations, thus reducing the accuracy of the robot's positioning and orientation. Research in this area has been addressed by Rodriguez-Donate et al. [15], who developed an intelligent processor using Kalman filters to filter and fuse information from a network of sensors. They used two primary sensors: an optical encoder and a 3-axis accelerometer. Calibration of a robot arm is an important factor in the accuracy of robot positioning. A simple, lowcost calibration procedure that improves the surface positioning accuracy of a SCARA (Selective Compliance Assembly Robot Arm) double-arm robot was published by Joubair et al. [16]. One of the key problems in examining the positioning accuracy of robotic arms is the working temperature. Kluz et al. Analyzed the analysis of the influence of temperature on the positioning accuracy of the robot arm [17]. The results obtained were subjected to statistical analysis using the Shapiro-Wilk test, which confirmed that the three-sigma rule can be used to calculate the value of the total positioning error of the robot arm.
The use of cooperating robotic solutions is also supported by current trend of automation and data exchange in production technologies, the so-called Industry 4.0 [18]. The goal of Industry 4.0 is ultimately to achieve efficiency, reduce costs and increase productivity through integrated automation. In Industry 4.0, production systems are characterized by individualized products in conditions of highly flexible mass production. A literature review of current standards for human-robot collaboration shows they can be used in a wide range of different regimes [19]. The field of application of human robot collaborative systems includes handling, welding, assembly and automotive [20]. Nowadays human robot collaboration systems are widely used in manufacturing companies that operate according to the concept of Industry 4.0. Due to the necessity of meeting the needs of individual customers, interest in the application for knowledge transfer support is growing. Ballestar et al. [21] provided knowledge related to the interconnection of industrial robotics and productivity of work in small and medium enterprises (SMEs). Patalas-Maliszewska and Krebs [22] developed the knowledge transfer approach is to select the best character-Appl. Sci. 2021, 11, 6104 3 of 23 istics of a knowledge worker in order to achieve an improvement of effective use of the application for knowledge transfer support among employees. This approach is based on survey and data obtained from Polish production plants.
The research so far has mostly focused on positioning accuracy control, research of control systems, control methods and positioning, but we have not been able to find publications aimed at examining the relationship between positioning accuracy of the robotic arm and manufacturing technologies with which it was manufactured. It is from this knowledge that the main research intention of the presented paper is based, the aim of which is to expand the field of knowledge concerning the accuracy of robotic arm positioning in relation to used production technologies.
The subject of the paper is the research of production technologies impact on positioning accuracy of a robotic arm. The aim was to find out whether different production technologies (additive and conventional) and the related design differences of the robotic arm affect its operational functionality. In the research, positioning accuracy of a robotic arm formed by three partial arms was specifically investigated, while the first partial arm Arm I was manufactured by two different technologies. On the robotic arm, the research was carried out in such a way that the first partial arm, Arm I, was being continuously changed and was available for research purposes in two variants. Each of the Arm I variants was manufactured using a different technology (additive and conventional) while, at the same time, the individual variants also differed in construction. The design differences of both variants were related to the production technology used. The measurement of positioning accuracy was performed with the use of two methods. Specifically, a contact and a non-contact method were used. The contact method was implemented on a 3D-measuring machine, RAPID, and the second contactless method was performed using an inductive sensor. The maximum working load of the robotic arm was 2 kg, therefore the positioning accuracy was examined at three degrees of operating load equal to 0, 50 and 100% of the maximum workload.
The research result proves that the technology of robotic arm production does not have a direct influence on its positioning. This result is based on results of mathematicalstatistical analysis. However, production technology affects its design. This fact might be considered the secondary aspect, which, however, can already affect the positioning accuracy. However, the hypothesis, which emerged from the research, needs to be further investigated. Within the presented research results, inaccuracies in the positioning of the robotic arm were manifested that are to be attributed to various constructions due to the production technology used.

Materials and Methods
The robot developed with the most used angular structure with rotational movements in 3 axes, Figure 1, was used for the experiment.
The robot is in shape of a human arm with swivel joints. The working space consists of spherical areas. Such robots are suitable for a wide range of activities involving the use of three rotary motion axes. The robot itself consists of a robot base and three arms. The base of the robot is usually anchored horizontally, Arm I is mounted on it and rotates around the vertical axis Z. The remaining two axes of rotation are horizontal and parallel to each other. They consist of Arm II and Arm III. Arm III can work in proximity to the Z axis. The individual arms of the robot are connected to each other by gear mechanisms driven by servo drives. As the arms move, their elastic deformation occurs, which also affects positioning accuracy of the end robotic arm with the manipulation effector. Positioning accuracy research was carried out on two identical versions of Arm I, which differed in production technology used. The robot is in shape of a human arm with swivel joints. The working space consists of spherical areas. Such robots are suitable for a wide range of activities involving the use of three rotary motion axes. The robot itself consists of a robot base and three arms. The base of the robot is usually anchored horizontally, Arm I is mounted on it and rotates around the vertical axis Z. The remaining two axes of rotation are horizontal and parallel to each other. They consist of Arm II and Arm III. Arm III can work in proximity to the Z axis. The individual arms of the robot are connected to each other by gear mechanisms driven by servo drives. As the arms move, their elastic deformation occurs, which also affects positioning accuracy of the end robotic arm with the manipulation effector. Positioning accuracy research was carried out on two identical versions of Arm I, which differed in production technology used.

Arm I-Made by Additive Manufacturing (AdM)
Production of Arm I was realized by a 3D printer Xline 2000R [23], and the material used was AlSi10Mg ( Figure 2

Arm I-Made by Additive Manufacturing (AdM)
Production of Arm I was realized by a 3D printer Xline 2000R [23], and the material used was AlSi10Mg ( Figure 2).
Technical specification: • layer height 0.06 mm, • support structure layer height 0.12 mm, • total time of finishing and removing of support structure, finishing and sandblasting 32 h,  The robot is in shape of a human arm with swivel joints. The working space consists of spherical areas. Such robots are suitable for a wide range of activities involving the use of three rotary motion axes. The robot itself consists of a robot base and three arms. The base of the robot is usually anchored horizontally, Arm I is mounted on it and rotates around the vertical axis Z. The remaining two axes of rotation are horizontal and parallel to each other. They consist of Arm II and Arm III. Arm III can work in proximity to the Z axis. The individual arms of the robot are connected to each other by gear mechanisms driven by servo drives. As the arms move, their elastic deformation occurs, which also affects positioning accuracy of the end robotic arm with the manipulation effector. Positioning accuracy research was carried out on two identical versions of Arm I, which differed in production technology used.

Arm I-Made by Additive Manufacturing (AdM)
Production of Arm I was realized by a 3D printer Xline 2000R [23], and the material used was AlSi10Mg ( Figure 2). Technical specification:   Y, Z, X in Figure 2 represent the Cartesian coordinate system, i.e., the orientation of Arm I for the individual measurement methods, which is identical with the manipulation of the robotic arm with the load in operation.

Arm I-Made by CNC (Computer Numerical Control) Milling (CvM)
Production of Arm I was realized by CNC (computer numerical control) machine Pinnacle VMC 650S, used material is AlMg4.5Mn, DIN 1732, with dimensions 150 mm × 150 mm − 450 mm ( Figure 3). Its chemical composition is in Table 1 and physicalmechanical properties in Table 2.

Contact Measuring Method (CoM) of Robotic Arm's Position
The contact measuring method (CoM) on a 3D-measuring machine can be implemented by a stylus with a ruby ball manually or by programming the stylus for repeated measurements. In our case, stylus programming by the learning method was used. Prior to measurement itself, the coordinate system of the 3D-measuring machine and the coordinate system of the measured-scanned body were identified. An NC (numeric control) scanning program was created in the control system of the 3D-measuring machine RAPID THOME Präzision [24], which ensured scanning of positions during repeated measurements. Stylus with a 4 mm diameter ruby ball was used for scanning as in Figure 4. The defined scanning sensitivity of the 3D-measuring machine 0.001 mm was achieved by observing the design and operating conditions. A holder was mounted on the robotic arm, which allowed mounting and weight exchange, and at the same time a part of it was machined so that reading of three perpendicular planes is possible, which after alignment with the reference coordinate system 3D-measuring machine were used to sense the position.

Contact Measuring Method (CoM) of Robotic Arm's Position
The contact measuring method (CoM) on a 3D-measuring machine can be implemented by a stylus with a ruby ball manually or by programming the stylus for repeated measurements. In our case, stylus programming by the learning method was used. Prior to measurement itself, the coordinate system of the 3D-measuring machine and the coordinate system of the measured-scanned body were identified. An NC (numeric control) scanning program was created in the control system of the 3D-measuring machine RAPID THOME Präzision [24], which ensured scanning of positions during repeated measurements. Stylus with a 4 mm diameter ruby ball was used for scanning as in Figure 4. The defined scanning sensitivity of the 3D-measuring machine 0.001 mm was achieved by observing the design and operating conditions. A holder was mounted on the robotic arm, which allowed mounting and weight exchange, and at the same time a part of it was machined so that reading of three perpendicular planes is possible, which after alignment with the reference coordinate system 3D-measuring machine were used to sense the position.
The measuring chain for CoM is shown in Figure 5. After setting the reference coordinate system of 3D-measuring machine (it consists in setting a fixed zero point on the base frame from which all measured values are read), the program for positioning the robotic arm is initialized and then adjusted to the position with maximum reach. The 3D measuring machine is then initialized. The stylus with ruby ball automatically senses the end position of the robotic arm [mm] and loads it into the PC application Metrolog XG. In this application, the end position of the robotic arm is evaluated. The measuring chain for CoM is shown in Figure 5. After setting the reference coordinate system of 3D-measuring machine (it consists in setting a fixed zero point on the base frame from which all measured values are read), the program for positioning the robotic arm is initialized and then adjusted to the position with maximum reach. The 3D measuring machine is then initialized. The stylus with ruby ball automatically senses the end position of the robotic arm [mm] and loads it into the PC application Metrolog XG. In this application, the end position of the robotic arm is evaluated.

Non-Contact Measuring Method (NcM) of Robotic Arm's Position Deviation
In the second method of measuring position deviation, NcM was used using the proximity sensor MTN/EP080 Probe [25] in Figure 6. The defined sensing sensitivity of the 0.01 mm proximity sensor was achieved by adhering to the design and operating conditions.   The measuring chain for CoM is shown in Figure 5. After setting the reference coordinate system of 3D-measuring machine (it consists in setting a fixed zero point on the base frame from which all measured values are read), the program for positioning the robotic arm is initialized and then adjusted to the position with maximum reach. The 3D measuring machine is then initialized. The stylus with ruby ball automatically senses the end position of the robotic arm [mm] and loads it into the PC application Metrolog XG. In this application, the end position of the robotic arm is evaluated.

Non-Contact Measuring Method (NcM) of Robotic Arm's Position Deviation
In the second method of measuring position deviation, NcM was used using the proximity sensor MTN/EP080 Probe [25] in Figure 6. The defined sensing sensitivity of the 0.01 mm proximity sensor was achieved by adhering to the design and operating conditions.
Setting the reference coordinate system of a 3D-measuring machine

Non-Contact Measuring Method (NcM) of Robotic Arm's Position Deviation
In the second method of measuring position deviation, NcM was used using the proximity sensor MTN/EP080 Probe [25] in Figure 6. The defined sensing sensitivity of the 0.01 mm proximity sensor was achieved by adhering to the design and operating conditions.
A holder with a replaceable weight machined in three planes perpendicular to each other was reused on the robotic arm. A stable stand made of non-magnetic material with a non-conductive proximity sensor holder was mounted on the base frame of a 3D-measuring machine. The holder with a replaceable weight is used (after setting the switching position of the proximity sensor) to sense the incremental position deviation in the direction of the X, Y and Z axes in micrometers.
In the measuring chain (Figure 7), the power supply to the proximity sensor is initialized first. Subsequently, the program for positioning the robotic arm is activated. This is followed by positioning the robotic arm to the position with maximum reach. After reaching this position, a graphical recording of the maximum reach of robotic arm measurement will be performed. The incremental value of position deviation for the maximum reach of the robotic arm in µm is read from the graphic record in the LabVIEW software. Appl. Sci. 2021, 11, x FOR PEER REVIEW 7 of 23 Figure 6. Non-contact measuring method (NcM) of robotic arm's position deviation using a proximity sensor in the direction of the X, Y, Z axis.
In the second method of measuring position deviation, NcM was used using the proximity sensor MTN/EP080 Probe [25] in Figure 6. The defined sensing sensitivity of the 0.01 mm proximity sensor was achieved by adhering to the design and operating conditions.
A holder with a replaceable weight machined in three planes perpendicular to each other was reused on the robotic arm. A stable stand made of non-magnetic material with a non-conductive proximity sensor holder was mounted on the base frame of a 3D-measuring machine. The holder with a replaceable weight is used (after setting the switching position of the proximity sensor) to sense the incremental position deviation in the direction of the X, Y and Z axes in micrometers.
In the measuring chain (Figure 7), the power supply to the proximity sensor is initialized first. Subsequently, the program for positioning the robotic arm is activated. This is followed by positioning the robotic arm to the position with maximum reach. After reaching this position, a graphical recording of the maximum reach of robotic arm measurement will be performed. The incremental value of position deviation for the maximum reach of the robotic arm in µ m is read from the graphic record in the LabVIEW software.   Robotic arm with Arm I-AdM and Arm I-CvM ( Figure 8) are shown for NcM at maximum reach. For this method of measurement, the base frame of the 3D measuring machine was used only to fix the position of the robotic arm while maintaining a constant position of the base frame of the robotic arm. At the same time, the base frame of the 3D measuring machine was used to place and fix the proximity sensor holder. Robotic arm with Arm I-AdM and Arm I-CvM ( Figure 8) are shown for NcM at maximum reach. For this method of measurement, the base frame of the 3D measuring machine was used only to fix the position of the robotic arm while maintaining a constant position of the base frame of the robotic arm. At the same time, the base frame of the 3D measuring machine was used to place and fix the proximity sensor holder.

Measurement Procedure and Evaluation of Deviations in the Position of the Robotic Arm
The robotic arm was at a maximum reach of 609 mm, according to Figure 1, programmed by the learning method for Arm I-AdM and Arm I-CvM. For CoM and NcM, the same program was used in the robotic arm control system. At the beginning of the measurement, the CoM of robotic arm's position was selected, and n measurements were performed with Arm I-AdM in individual directions of the X; Y; Z coordinate system. NcM of robotic arm's position deviation followed. After its completion, the disassembly of Arm I-AdM and assembly of Arm I-CvM followed. With this arm, n measurements were performed using NcM of robotic arm's position deviation in each direction of the X; Y; Z coordinate system. Then the CoM of robotic arm's position followed. The position of the proximity sensor was always adjusted to each scanned position of the robotic arm in the corresponding coordinate system. The data obtained were recorded and statistically processed for n measurements, for each axis of the coordinate system, for CoM, NcM, AdM and CvM, while the individual measurements were marked with the order index i.

Measurement Procedure and Evaluation of Deviations in the Position of the Robotic Arm
The robotic arm was at a maximum reach of 609 mm, according to Figure 1, programmed by the learning method for Arm I-AdM and Arm I-CvM. For CoM and NcM, the same program was used in the robotic arm control system. At the beginning of the measurement, the CoM of robotic arm's position was selected, and n measurements were performed with Arm I-AdM in individual directions of the X; Y; Z coordinate system. NcM of robotic arm's position deviation followed. After its completion, the disassembly of Arm I-AdM and assembly of Arm I-CvM followed. With this arm, n measurements were performed using NcM of robotic arm's position deviation in each direction of the X; Y; Z coordinate system. Then the CoM of robotic arm's position followed. The position of the proximity sensor was always adjusted to each scanned position of the robotic arm in the corresponding coordinate system. The data obtained were recorded and statistically processed for n measurements, for each axis of the coordinate system, for CoM, NcM, AdM and CvM, while the individual measurements were marked with the order index i. The proposed procedure for measuring and evaluation of the robotic arm's position deviation is shown in Figure 9.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 23 The proposed procedure for measuring and evaluation of the robotic arm's position deviation is shown in Figure 9.

Ishikawa Diagram
The generated Ishikawa diagram ( Figure 10) defines, investigates and detects the effects of several influences and causes, which result in the variability of the robotic arm position. The expected sources of variability in deviations of the robotic arm position that affect it are: The source of variability, which is the subject of the experiment, is the accuracy of the position of the robotic arm with respect to the used production technology. However, in addition to expected resources, other factors also affect the accuracy of the robotic arm's

Ishikawa Diagram
The generated Ishikawa diagram ( Figure 10) defines, investigates and detects the effects of several influences and causes, which result in the variability of the robotic arm position. The expected sources of variability in deviations of the robotic arm position that affect it are: The source of variability, which is the subject of the experiment, is the accuracy of the position of the robotic arm with respect to the used production technology. However, in addition to expected resources, other factors also affect the accuracy of the robotic arm's position. The causes are divided into categories and represent potential sources of variability for measured position deviations in the robotic arm.  Measuring apparatus is a factor influencing the complexity of measurement, the methodology design, requiring the experience of the implementer of measurements. It also affects the time aspect.
Robotic arm construction influences the production, operation, method of measurement, the possibility of using the apparatus, the choice of measuring members, the choice of the material of the robotic arm, while directly affecting its mechanical properties.
The human factor is the sum of characteristics of the person performing the measurement. In it, physical properties can be included, e.g., the promptness and measurement speed. Furthermore, experience that can be not only theoretical but also practical. Experience, within the researched issue, means practical experience in the field of measurement, metrology, creation of measurement methodology, planning of experiments with the aim of their harmonization for the measurement of required quantities.
Position is a factor reflecting the location of the robotic arm. It includes external influences that can be mechanical or physical. Their effect on the position of the robotic arm can be in terms of duration: temporary or permanent. In terms of severity, it can be minimal, neutral or serious.
Course of measurement is a factor that affects the method of measurement, the accuracy of results, their reliability, the equipment used and the complexity of the implementation.
Other sources of variability are undesirable in the experiment and require correction using the following methods: • elimination-the conditions of the experiment will ensure that this source of variability does not occur in the experiment at all; Measuring apparatus is a factor influencing the complexity of measurement, the methodology design, requiring the experience of the implementer of measurements. It also affects the time aspect.
Robotic arm construction influences the production, operation, method of measurement, the possibility of using the apparatus, the choice of measuring members, the choice of the material of the robotic arm, while directly affecting its mechanical properties.
The human factor is the sum of characteristics of the person performing the measurement. In it, physical properties can be included, e.g., the promptness and measurement speed. Furthermore, experience that can be not only theoretical but also practical. Experience, within the researched issue, means practical experience in the field of measurement, metrology, creation of measurement methodology, planning of experiments with the aim of their harmonization for the measurement of required quantities.
Position is a factor reflecting the location of the robotic arm. It includes external influences that can be mechanical or physical. Their effect on the position of the robotic arm can be in terms of duration: temporary or permanent. In terms of severity, it can be minimal, neutral or serious.
Course of measurement is a factor that affects the method of measurement, the accuracy of results, their reliability, the equipment used and the complexity of the implementation.
Other sources of variability are undesirable in the experiment and require correction using the following methods: • elimination-the conditions of the experiment will ensure that this source of variability does not occur in the experiment at all; • minimization-targeted reduction of variability so that the rest will be part of the experimental error; • part of the experimental error-we know about this source, it is impossible to treat, thus will be reflected in a random error in the calculations.
As part of ensuring the conditions for carrying out the experiment, each known source of variability was corrected accordingly.
The measuring apparatus for measuring a robotic arm's positioning accuracy by NcM, using a proximity sensor, consisted of a voltage source, connecting electrical cables, an A/D (an analog-to-digital) converter and a proximity sensor. The measuring apparatus for measuring of robotic arm's positioning accuracy by CoM consisted of a 3D-measuring machine, RAPID, and a transducer with a complete connection of sensors.

Statistical Evaluation of Measurement Results
Statistical evaluation of measurement results was performed for both Arm I-AdM and Arm I-CvM with the use of CoM and NcM. In each case, 11 measurements described in chapter 2.5 were performed, while each group of 11 measurements formed one set. The results were processed according to Figure 9 as follows: a.
Verification that the selection of measured values comes from a population with a normal distribution.
• Check outlier method of comparing the distance of the minimum and maximum from the first and third quartiles.

•
The following descriptive statistics are used to describe the measurement groups for each setting. Position statistics-average, median, minimum, maximum. The median is the mean value in the data sorted by size, and together with the minimum and maximum values, they give us a view of the data on how densely they are distributed. • Statistics on the variability of results-range, standard deviation. While range is the difference between the minimum and maximum value and only expresses the width of the data occurrence, standard deviation refers to the data density around the average.

•
Other statistics show the shape of the frequency curve, skewness about the location of the furthest sampling value from the arithmetic mean, and kurtosis about the density of data around the mean. c.
Comparison of measurement pairs.
• For Arm I-AdM, Arm I-CvM, where the other mode settings were the same, the Two Sample t-test was used for the null hypothesis (difference in means is equal to 0) µ_1 = µ_2. To calculate the test statistic t, it was necessary to verify whether var.equal applies. This was preceded by an F test to compare two variances with the ratio of variances equal to 1.

Validation of Results by Computer Simulations
During the research implementation, verification of measured values was also carried out. FEM (finite element method) was used to verify them. The results were compared with the results of experimental calculations. As an example, we present the Z-axis deformation magnitudes in AdM_NcM determined by experimental measurements, the values of which are given in Table 3. Then, the calculation of robotic arm load was performed using FEM, and the magnitude of the deformation was determined. An example of the result of the calculation is shown in Figure 11.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 12 of 23 compared with the results of experimental calculations. As an example, we present the Zaxis deformation magnitudes in AdM_NcM determined by experimental measurements, the values of which are given in Table 3. Then, the calculation of robotic arm load was performed using FEM, and the magnitude of the deformation was determined. An example of the result of the calculation is shown in Figure 11. The comparison of experimentally measured value of deformation and with FEM calculation shows the difference of 8.5 %.
Verification was performed in the same way for the CvM_NcM arm. The magnitude of the measured deformations is given in Table 4.   The comparison of experimentally measured value of deformation and with FEM calculation shows the difference of 8.5 %.
Verification was performed in the same way for the CvM_NcM arm. The magnitude of the measured deformations is given in Table 4.  By comparing the experimentally measured value of deformation and FEM calculation, it is possible to state a difference of 1.4%.

Verification That the Selection of Measured Values Comes from a Population with a Normal Distribution
In this part, it was necessary to verify the hypothesis of a normal population distribution. The hypothesis of a normal population distribution (Table 5) was not rejected, i.e., the measurements were performed correctly, and the results of measurements can be further used for statistical evaluation using parametric tests.
Tables 5 and 6 present the results for groups of measurements without outliers. This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

Verification That the Selection of Measured Values Comes from a Population with a Normal Distribution
In this part, it was necessary to verify the hypothesis of a normal population distribution. The hypothesis of a normal population distribution (Table 5) was not rejected, i.e., the measurements were performed correctly, and the results of measurements can be further used for statistical evaluation using parametric tests.  Tables 5 and 6 present the results for groups of measurements without outliers. This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn. Since the comparison using non-parametric methods is not possible for the measurement results, a comparison for measurements in the direction of the Y axis was not performed for the aforementioned sets, they are pairs: CoM_0_Yy, CoM_1_Yy, CoM_2_Yy. Statistical conclusions show for practice that in the measurements in the Y-axis direction, there was a dimensional anomaly during the exchange of the arms, caused by noncompliance with assembly procedures and prescribed tightening torques of screws securing the arm.
The hypothesis of a normal population distribution for all NcM measurements in Table 6 was not rejected, so the measurements were performed correctly, and the measurement results can be further used for statistical evaluation using parametric tests. When replacing the arms and ensuring the correct assembly procedures and the prescribed tightening torques of the screws, measurements were achieved that also met the statistical requirements.

Descriptive Statistics
In Tables 5 and 6, for every measurement the number of n values is determined which were left out after excluding the outliners. The descriptive statistics used for the robotic arm position -were: mean, median, minimum, maximum. Additional descriptive statistics of robotic arm position variability were used in the article: range, standard deviation, and descriptive statistics of shape: skewness and kurtosis. Descriptive statistics are calculated separately for CoM and NcM due to different measurement methodology.

CoM Descriptive Statistics
Descriptive one-dimensional statistics of position and variability CoM of Arm I-AdM and Arm I-CvM are in Table 7. As follows from the CoM principle (Section 2.3), the results in Table 7 are given in mm.
It is clear from Table 7   For measurements in the direction of the X axis ( Figure 13), for Arm I-AdM and Arm I-CvM, a significant decrease in the median position of the arm showed when changing the load from 0 kg to 1 kg and a slight decrease in median position is still visible when changing the load from 1 kg to 2 kg. The variability of measurements in the direction of the X axis evaluated by range and standard deviation sd is highest for the measurement AdM_CoM_1_Xx. With both types of arm, the robot arm decreases in the X-axis direction as the load increases.   Table 7    For the measurement in the direction of the Y axis (Figure 14), for Arm I-AdM, there is a significant increase in the median position when changing the load from 0 kg to 1 kg and a slight decrease in the median position when changing the load from 1 kg to 2 kg. Arm I-CvM shows a significant decrease in the median position when changing the load from 0 kg to 1 kg, and a slight decrease in the median position when changing the load from 1 kg to 2 kg. The highest range is for measurements in the direction of the Y axis for CvM_CoM_2_Yy, namely range 0.44 mm and standard deviation sd of 0.14 mm. For measurements in the direction of the X axis (Figure 13), for Arm I-AdM and Arm I-CvM, a significant decrease in the median position of the arm showed when changing the load from 0 kg to 1 kg and a slight decrease in median position is still visible when changing the load from 1 kg to 2 kg. The variability of measurements in the direction of the X axis evaluated by range and standard deviation sd is highest for the measurement AdM_CoM_1_Xx. With both types of arm, the robot arm decreases in the X-axis direction as the load increases. For the measurement in the direction of the Y axis (Figure 14), for Arm I-AdM, there is a significant increase in the median position when changing the load from 0 kg to 1 kg and a slight decrease in the median position when changing the load from 1 kg to 2 kg. Arm I-CvM shows a significant decrease in the median position when changing the load from 0 kg to 1 kg, and a slight decrease in the median position when changing the load from 1 kg to 2 kg. The highest range is for measurements in the direction of the Y axis for CvM_CoM_2_Yy, namely range 0.44 mm and standard deviation sd of 0.14 mm. For measurements in the direction of the Z axis (Figure 15), for Arm I-AdM, there is a significant increase in the median position when changing the load from 0 kg to 1 kg and slightly decreasing in the median position when changing the load from 1 kg to 2 kg. Arm I-CvM shows a slight decrease in the median position when changing the load from 0 kg to 1 kg and a slight decrease in the median position when changing the load from 1 kg to 2 kg. The range in the direction of the Z axis is the largest for measurement AdM_CoM_2_Zz, namely the range up to 0.71 mm and the standard deviation sd of 0.20 mm. Figures 13-15 show that both Arm I-AdM and Arm I-CvM, without load, will reach a different position than when measuring with load. The resulting increase in the position of the robotic arm in the Z axis direction using Arm I-AdM with respect to Arm I-CvM, is caused by springing of the Arm I-AdM structure. For measurements in the direction of the X axis (Figure 13), for Arm I-AdM and Arm I-CvM, a significant decrease in the median position of the arm showed when changing the load from 0 kg to 1 kg and a slight decrease in median position is still visible when changing the load from 1 kg to 2 kg. The variability of measurements in the direction of the X axis evaluated by range and standard deviation sd is highest for the measurement AdM_CoM_1_Xx. With both types of arm, the robot arm decreases in the X-axis direction as the load increases. For the measurement in the direction of the Y axis (Figure 14), for Arm I-AdM, there is a significant increase in the median position when changing the load from 0 kg to 1 kg and a slight decrease in the median position when changing the load from 1 kg to 2 kg. Arm I-CvM shows a significant decrease in the median position when changing the load from 0 kg to 1 kg, and a slight decrease in the median position when changing the load from 1 kg to 2 kg. The highest range is for measurements in the direction of the Y axis for CvM_CoM_2_Yy, namely range 0.44 mm and standard deviation sd of 0.14 mm.

NcM Descriptive Statistics
Descriptive one-dimensional NcM position and variability statistics for Arm I-AdM and Arm I-CvM are in Table 8. As follows from the CoM principle (Section 2.3), the results in Table 5 are given in mm. From Table 8 and Figure 16, for Arm I-AdM measurements in the direction of the X axis, it is possible to observe a higher value of mean and lower range and standard deviation sd when the load is increased to 1 kg; and range and standard deviation sd increased when the load was increased to 2 kg and the median shifted. For Arm I-CvM in the direction of the X axis, it is possible to observe a decrease in the median value as well as the mean when the load increases, and at the same time a decrease of range and standard deviation sd. This follows from the above behavior of the position of the robotic arm Arm I-AdM and Arm I-CvM and is due to the rigidity of the sensor holder. From Table 8 and Figure 17, for Arm I-AdM, measurements in the Y axis direction, an increase in the median value can be observed as the load increases and the range and standard deviation sd change slightly. For Arm I-CvM in the direction of the Y axis, it is From Table 8 and Figure 17, for Arm I-AdM, measurements in the Y axis direction, an increase in the median value can be observed as the load increases and the range and standard deviation sd change slightly. For Arm I-CvM in the direction of the Y axis, it is possible to observe an increase in the median and the mean with increasing load, while the range and standard deviation sd is higher at higher loads. Failure to meet the condition of normal distribution of measured values is also confirmed by anomalies in Figure 17, at the same time, the stress caused by the construction of Arm I-AdM needs to be considered. From Table 8 and Figure 17, for Arm I-AdM, measurements in the Y axis direction, an increase in the median value can be observed as the load increases and the range and standard deviation sd change slightly. For Arm I-CvM in the direction of the Y axis, it is possible to observe an increase in the median and the mean with increasing load, while the range and standard deviation sd is higher at higher loads. Failure to meet the condition of normal distribution of measured values is also confirmed by anomalies in Figure  17, at the same time, the stress caused by the construction of Arm I-AdM needs to be considered.
From Table 8 and Figure 18, for Arm I-AdM measurements in the direction of the Z axis, it can be observed that the median value always decreased and the range changed slightly due to the increase in load. For Arm I-CvM in the direction of the Z axis, it is possible to observe that the median value and the mean oscillate with increasing load and From Table 8 and Figure 18, for Arm I-AdM measurements in the direction of the Z axis, it can be observed that the median value always decreased and the range changed slightly due to the increase in load. For Arm I-CvM in the direction of the Z axis, it is possible to observe that the median value and the mean oscillate with increasing load and the range is not always the same. The expected course, met by Arm I-AdM, is caused by the interruption of the Arm I-AdM suspension, and at the same time suspension of the sensor holder. The increase in Arm I-CvM at 2 kg load was again caused only by the springing of the sensor holder.
Presented results indicate a statistical significance of differences in positions of the robotic arm endpoint for Arm I-AdM and Arm I-CvM, in position and variability. Since values in the table and the graph are in micrometers, these differences are negligible from a practical point of view. For further investigation, it would be appropriate to plan experimental measurements to better examine whether the differences are random or caused by a factor that has not yet been considered. This was affected by: arm flexibility, loading time, movement time of the robotic arm, sensitivity of the proximity sensor and rigidity of the sensor holder.
Due to the fact that the production documentation of Arm I-AdM and Arm I-CvM contained dimensional tolerances in micrometers and was observed, the position shift could occur during the mutual exchange, disassembly and assembly of Arm I-AdM for Arm I-CvM.
the range is not always the same. The expected course, met by Arm I-AdM, is caused by the interruption of the Arm I-AdM suspension, and at the same time suspension of the sensor holder. The increase in Arm I-CvM at 2 kg load was again caused only by the springing of the sensor holder. Presented results indicate a statistical significance of differences in positions of the robotic arm endpoint for Arm I-AdM and Arm I-CvM, in position and variability. Since values in the table and the graph are in micrometers, these differences are negligible from a practical point of view. For further investigation, it would be appropriate to plan experimental measurements to better examine whether the differences are random or caused by a factor that has not yet been considered. This was affected by: arm flexibility, loading time, movement time of the robotic arm, sensitivity of the proximity sensor and rigidity of the sensor holder.
Due to the fact that the production documentation of Arm I-AdM and Arm I-CvM contained dimensional tolerances in micrometers and was observed, the position shift could occur during the mutual exchange, disassembly and assembly of Arm I-AdM for Arm I-CvM. Table 9 shows the results of the F test and the t-test for CoM. It is clear from Table 9 that no match of the mean values is confirmed for any pair. In four cases, for the CoM_0_Zz, CoM_1_Zz, CoM_2_Xx and CoM_2_Yy modes, the agreement of the variances is confirmed.   Table 9 shows the results of the F test and the t-test for CoM. It is clear from Table 9 that no match of the mean values is confirmed for any pair. In four cases, for the CoM_0_Zz, CoM_1_Zz, CoM_2_Xx and CoM_2_Yy modes, the agreement of the variances is confirmed.  Table 10 shows the results of the F test and the t-test for NcM. It is clear from Table 10 that no agreement of the mean values is confirmed for any pair. In six cases, for NcM_0_Y, NcM_0_Z, NcM_1_Y, NcM_2_Y, NcM_2_Z, and NcM_2_X modes, the variance agreement is confirmed. The test results in Tables 9 and 10 confirm that Arm I-AdM and Arm I-CvM are not interchangeable in this study.

Conclusions
The proposed methodology for measuring of robotic arm positioning accuracy was verified on its construction consisting of three arms Arms I, II, III, while Arm I was manufactured by two different technologies AdM and CvM. A series of repeated measurements was performed for each robotic arm configuration with two different measurement methods, CoM and NcM. The results of measurements were verified by statistical methods, based on which unsatisfactory values of measurements were excluded from the evaluation.
The use of proposed methodology is not only in the field of metrology and testing, but also to verify the interchangeability of components in construction of robotic systems. The methodology provides a detailed view of a dimensional chain quality of the robotic arm structure and determines conditions for maintaining accuracy during disassembly and reassembly of individual components of robotic system structures.
Statistical evaluation of the results verified that the data obtained were measured correctly and have a normal distribution.
For CoM with the same load of Arm I-AdM and Arm I-CvM: • homoscedasticity was confirmed for: CoM_0_Zz, CoM_1_Zz, CoM_2_Xx and CoM_2_Yy; • the conformity of the mean values has not been confirmed.
This means that for interchangeable Arm I it is necessary to modify the design of ArM I-AdM.
For NcM with the same load of Arm I-AdM and Arm I-CvM: • homoscedasticity was confirmed for: NcM_0_Y, NcM_1_Y, NcM_2_Y, NcM_0_Z, NcM_2_Z; • the conformity of the mean values for NcM_2_X has been confirmed.
This means that for mutually interchangeable Arm I-AdM, Arm I-CvM, to increase its rigidity, the structural design of the sensor holder must be modified.
Considering the results obtained, the paper follows up the work of Zhang and Wei [26]. The authors addressed the accuracy of robotic arm positioning in their work, too. However, they did it from the point of view of its control system. A similar issue is also covered in the work of Clitan and Ionut [27]. The authors also focus on accuracy of robotic arm positioning. Specifically, they dealt in detail with its payload. Another paper focusing on this issue is by Visan et al. [28]. Even in this case, the problem of robotic arm positioning accuracy was solved in detail. In this case, however, the research was focused on stepper motors that implement the positioning of the robotic arm. It follows from the above examples that the presented research brings a new area of research to the issue of positioning the robotic arm, which needs to be further addressed in order to further expand the knowledge that can be transferred to the application.
Joubair et al. [16] is another paper that dealt with positioning accuracy. The authors took a closer look at the calibration procedure, where a simple low-cost procedure improves the accuracy of surface positioning. This is the issue similar to one presented in the article. The difference is mainly that the authors used a two-arm robot, compared to our researched single-arm robotic arm.
The following contributions also deal with the calibration of the robotic arm. In our paper, the proximity sensor MTN/EP080 Probe was used for NcM. Aoyagi et al. [29] used a laser-tracking system to calibrate the kinematic parameters of the robotic arm, which achieves high positioning accuracy using a genetic algorithm.
Švaco et al. [30] used a contactless method to perform measurements of calibration points in space, using a stereovision system attached to the robot arm. This method is very similar to the NcM for which the MTN/EP080 Probe proximity sensor was used. Points (represented as spheres localized by a stereo system) are projected by the authors as circles in two planes of image capturing, regardless of the angle of view. The positioning error after calibration has been reduced to 1.29 mm.
Peng et al. [31] proposed their own method of geometric parameters' calibration of a kinematic model of a robotic arm, based on monocular vision. Similar to our paper, they used the contactless method, while in measuring accuracy of the robotic arm positioning in our paper, the CoM was also used, in which it was necessary to set the reference-calibration Cartesian coordinate system. Peng et al., to determine the kinematic parameters, first used the classical Denavit-Hartenberg (D-H) modeling method. Subsequently, they implemented nonlinear optimization and parameter compensation. Their method improved the absolute accuracy of positioning the end of the robotic arm while being universal and effective, similar to the CoM method in our paper.
Due to confirmed differences in positions for the same measured axis in different types of arm, we assume that Arm I-AdM, Arm I-CvM or an arm made by different technology may have dimensional deviations. From previous observations and calculations, for the research of other factors influencing the localization of the robotic arm position, these deviations in dimensions need to be eliminated.
Based on the results, we can conclude that production technology does not affect the positioning accuracy of the robotic arm, but the design of Arm I-AdM needs to be changed regardless of the operating load.
To further investigate, it will be appropriate to plan experimental measurements to better examine whether differences in a position of the robot arm are random or caused by loads or other factors that have not been taken into account (arm flexibility, load time, transport time, programmed trajectory of the robot arm etc.). The answer to the question will be known after analysis of accuracy of measured positions for Arm I-AdM and Arm I-CvM, depending on the measurement directions (X, Y, Z) and the load.
In the future, to increase the accuracy of measurement process, it will be necessary to use the proximity sensor MTN/EP080 Probe with higher sensitivity and to increase rigidity of the proximity sensor holder. Furthermore, we will need to ensure the same position of the sensor holder for Arm I-AdM, Arm I-CvM on the robotic arm.

Conflicts of Interest:
The authors declare that they have no conflict of interest. "made by additive manufacturing, CoM measurement, load 0 kg, measuring axis X" CvM_CoM_2_Zz "made by CNC milling, CoM measurement method, load 2 kg, measuring axis Z" AdM_NcM_0_X "made by additive manufacturing, NoC measurement method, load 0 kg, measuring axis X" CvM_NcM_2_Z "made by CNC milling, NoC measurement method, load 2 kg, measuring axis Z"