Preprint
Article

Evaluating the Accuracy of the Reverse Engineering Process of Worn Non-Standard Spur Gears Using Conventional Techniques and Measuring Instruments—Pilot Studies

Altmetrics

Downloads

77

Views

27

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

26 April 2024

Posted:

26 April 2024

You are already at the latest version

Alerts
Abstract
The authors undertook pilot studies to assess the accuracy of the reverse engineering process of worn nonstandard spur gears (non-standardized modulus value, unusual profile/pressure angle, high or low teeth), using conventional techniques and measuring instruments. Eight gears were tested, the module of which ranges from 1.020 to 4.98 mm and the number of teeth from 13 to 42. The key parameter, the basic pitch, was selected to estimate the accuracy of the process. The goal is to determine the value of the profile angle. Eleven models were proposed to estimate the nominal tolerance field, using various types of random data distribution. The tested gears were made in IT grade: 6, 7, 8 and 9 according to DIN 3961. Vernier disc micrometers with various measurement ranges were used for research. It has been shown that the nominal module does not have to be treated as a random variable in the population. EIG (Equation of Identity) was developed, allowing conversion of any gear with specific values of geometric parameters into an identical gear with alternative values of these parameters. The directions and degrees of shifts of the estimated tolerance fields relative to the nominal field were determined. The most effective estimating model was selected in the context of the percentage filling of the nominal tolerance field, taking into account the symmetric Student-Fisher distribution with a confidence level of 60%. Further directions of the research activities were described.
Keywords: 
Subject: Engineering  -   Mechanical Engineering

1. Introduction

Although RE is used in many branches of modern industry, in gear engineering it particularly concerns the preparation of a technical drawing of a physically existing gear to be measured. On the basis of the technical documentation prepared in this way, it is possible to reproduce this part as a spare part (Figure 1).
The need to make a spare part for the gear wheel on your own may be caused by the lack of availability of the original spare part in the manufacturer's service box.
C M M Coordinate Measuring Machine r z Base diameter deviation m m
E I G Equation of Identical Gears t 0 Single pitch deviation m m
L S L Lower Specification Limit m m t p Chordal base pitch deviation m m
M P E Maximum Permissible Error m m Σ t 0 Total cumulative pitch deviation m m
N S Nominal Size m m
P C R Process Capability Ratio ( C p ; C p k ) c Radial clearance m m
R E Reverse Engineering c * Radial clearance factor
S - F Student–Fisher dispersion d Pitch diameter m m
U S L Upper Specification Limit m m d a Tip diameter m m
d b Base diameter m m
d e s i g Used as a subscript means: estimated d f Root diameter m m
m a x Used as a subscript means: maximal value d x Any given diameter m m
m i n Used as a subscript means: minimal value h Whole depth m m
t r u e Used as a subscript means: original h a Addendum m m
h f Dedendum m m
α Profile angle ° M n Span measurement m m
α ¯ Profile angle in radians r a d m Module m m
α x Profile angle at a given diameter ° N Number of measurements
α x ¯ Profile angle at a given diameter in radians r a d n Number of teeth spanned
φ Tolerance increment between gear tooth qualities p b Base pitch m m
σ Standard deviation m m S S–F deviation m m
s Tooth thickness at pitch diameter m m
e Total profile deviation m m s x Tooth thickness at a given diameter m m
g Tooth thickness deviation m m T Tolerance zone m m
M Span measurement deviation m m X ¯ Average m m
o Pitch error m m x Shift factor
r Runout error m m y Whole depth factor
Although this practice is known and widely used all over the world, there are no publications in the available scientific databases and literature that indicate the possibility of unambiguously reproducing the geometry of a spur gear with straight teeth using only conventional techniques and measuring instruments. In such a process, the aim is to reduce the cost and time associated with measurements. In this work, however, the authors went a step further and touched on the topic of RE for non-standard spur gears. The key issue for this topic is tolerance assignment in RE. The authors of Krogstie et al. [1] pay particular attention to the great importance of tolerance engineering on the way to increasing the efficiency of designing and manufacturing mechanical parts.
In the work of M.A. Sáenz-Nuño et al. [3] presented a useful approach to setting tolerances for any mechanical part within RE, practical for industry. Measurements of the parts were made by sufficiently qualified students. Each geometric parameter of the parts tested was measured N times, where the number N varies depending on the part being measured, one of the three: for PR 06 N 16   ; 32 , for PR 08 N 5   ; 12 and for PR 05 N 4   ; 7 . This is an important suggestion in the context of selecting the type of data dispersion in a statistical study. The four key parameters for developing the results in the cited work were: (i) Average, (ii) Median, (iii) Disper and (iii) Range. In Section 3.2. IT specification of the cited work, the authors conclude that it is more effective to use the median as the nominal dimension rather than the average. In their analysis, the authors assigned a symmetrical tolerance JS/js according to ISO 286 to each measurement. In the case of gear wheels, there is a specific analogy because the deviations of the geometric parameters of gear wheels are also symmetric. However, they themselves state in their publication that the scope of the proposed strategy is to be a simple and quick IT estimation procedure, which is why complex estimations of the normal distribution of measurements were not addressed. The classical principle is also preserved in the work cited: M P E 10 % · I T . In the case of this work, this may not necessarily be maintained due to the much higher precision of gears used in the aviation industry compared to the relatively large error limits of disc micrometers.
The works of George J. Kaisarlis et al. [4], [5], [6] describe an extensive approach to the allocation of tolerances for test parts (hole and shaft) within RE. They presented a method based on five steps: Step (a) consisting of determining the direction of analysis of measurements of diameter, geometric deviations, and surface roughness and determining the measurement uncertainty in relation to the ISO 286 (Fundamental Deviation) standard and/or the ISO 2768 (General Tolerances) standard, Step (b) consisting in designation Sets of candidate IT grades, Step (c) consisting in designation Sets of suggested fits, Step (d) in which the set of suggested fits is limited to obtain Sets of preferred fits and Step(s) of selecting the best tolerance match relative to production costs. This work contains many mathematical dependencies that determine the zero-one type of behavior (meeting one condition determines a further action, meeting an alternative condition determines a different action). However, you must pay attention to the following facts: (i) gears are tolerated according to other specialized standards regarding the specific method of dimensioning and tolerating the macro and microgeometry of gears, (ii) gears cannot be referred to as a male or female member of the assembly against which fits are determined, (iii) the issue of the possibility of using statistical techniques to determine the dispersion of measurement data should also be raised, and (iv) the cited work includes measurements using CMM. Therefore, it is a need to formulate an adequate tolerance allocation model for gears, using conventional techniques and measuring instruments.
Another publication by the same authors George J. Kaisarlis et al. [7] talks about the use of Knowledge-Based System (KBS) for tolerance allocation in RE also in relation to the use of CMM in measurements. However, the cited work contains a specific algorithm for assigning tolerances for the tested parts within RE (Figure 3 on page 531), which is the inspiration for the authors of this work on the way to creating a tolerance allocation model for gears within RE. An aspect contained in the cited work, which the authors considered potentially troublesome in the study of gears, is the use of increments of 0.5 and 0.1 (left side of the algorithm, Figure 3 on page 531 in the cited work) with respect to the NS. Even minor manipulations of the base pitch value of the order of thousandths of a millimeter significantly affect the accuracy of recreating the geometry of the tested wheel. Moreover, conventional measuring instruments have much larger error limits than CMMs.
In [8] by Zhaohui Geng et al. scientists rightly point out that since all parts were received from the original supplier or dismantled from the machines, they have certainly passed quality control by the manufacturer. Therefore, the areas defined by the upper and lower limits of the dimensional values tested should be at least within the original tolerance specifications. However, these types of estimates are conservative because the estimated tolerance areas may be too narrow, which could unnecessarily require high-precision RE manufacturing processes. Four standardized boundary estimators were used to examine a part with nine symmetrically tolerated dimensions and five geometric tolerances: (i) ASME Estimator, (ii) ISO Range Estimator, (iii) ISO 6σ Estimator and (iv) ISO PCR Estimator. As in one of the previously cited works, here the measurements were made using a CMM. By which the researchers indicate that the measurement error of a coordinate measuring machine can be approximated by the DGPPR (datum-guided partial Procrustes registration) method proposed by Geng and Bidanda [9], which is much smaller than the variability between parts. Therefore, measurement error can be ignored in the study of tolerance estimation of original manufacturing processes. The presentation of the estimated tolerance fields in the RE process is provided in Table 1 and Table 2 of the cited work. Table 2 calculates the ratios of the estimated tolerance fields to the original tolerance fields included in the manufacturing drawing of the tested part. Values are given as percentages. The estimator ratio is not determined ("–„ in the table) if both extreme values of the estimated tolerance field are not within the range of the original tolerance field. For example, the ISO PCR Estimator (Ratio) for Size 1 in Table 2 of the cited work with this, there is an error there. The value should be 0.5%, not 0.05%. Evidence: for Size 3 = 0.550 ± 0.001 (according to Figure 9) the original tolerance zone is equal T t r u e = 0.002 . The ISO Range Estimator shows an estimated tolerance zone of T d e s i g = 0.55000 0.54999 = 0.00001 (according to Table 1). And so:
I S O   R a n g e   E s t i m a t o r   R a t i o = 0.00001 0.00200 · 100 % = 0.5 %
However, both the value of 0.05% and the corrected value of 0.5% indicate a very narrow overlap (interference) of the tolerance fields, which would require increasing the machining accuracy of the re-manufactured part by several classes. In the work cited, the scientists assumed that the data distribution was normal and symmetric, and the average value was approximately equal to the unknown nominal value, which is why they used the coefficient C p , unlike the coefficient C p k . They also adopted a standard coefficient value of 1.33. In the cited work, the scientists presented a wide range of discussions and conclusions related to the estimators of the ASME [10] and ISO [11] standards. The authors will apply the same selected estimators here to verify the applicability to non-standard gears.
The aim of this work is to initially check whether it is possible to correctly recreate the geometry of a worn non-standard spur gear with straight teeth within RE, using only conventional techniques and measuring instruments, commonly available in most design offices and laboratory workshops. The most effective model (or a narrowed group of several effective models) should be selected among the available statistical models that allow the estimation of the tolerance zone.
This work consists of six sections. Previous Section 1. Introduction. This is followed by Section 2. Materials and Methods, which presents: the tested gears, the conventional techniques and measuring instruments used, the comparative model used and the statistical models used to develop the measurement results. Next, Section 3. Results and Discussion, in which the authors present the results of the estimated tolerance fields and develop a method for assessing the accuracy of the presented process. Finally, Section 4. Conclusions, which summarizes the whole and presents plans for further research.

2. Materials and Methods

2.1. Tested Gear Wheels

The study examined worn non-standard spur gears obtained from the aviation industry under an agreement. Due to legal restrictions, it is not possible to show the original technical drawing. The photos of the wheels were also cropped so that it was not possible to graphically associate the tested wheels with any commercial product (Figure 2).
The word non-standard here refers to different aspects: normal, high, or low teeth; a relatively small profile shift factor; unusual value of the profile angle; custom module value; value of the radial clearance coefficient exceeding the commonly accepted range (Table 1). Wheels are worn to varying degrees, but none of them has noticeable chipped teeth (without scoring or pitting).

2.2. Conventional Measuring Instruments Used

A regular electronic caliper was used for the preliminary measurements to identify the modulus value. For further measurements, Vernier disc micrometers with various measurement ranges were used, depending on the size of the gear wheel: 25-50 mm, 50-75 mm and 75-100 mm. All of the above-mentioned devices belong to the group of conventional, commonly used measuring devices. It is about reducing costs and time in the reverse engineering process, as well as the practicality of this process in the everyday design approach.

2.3. Conventional Measurement Techniques Used

Regarding the initial measurements with a standard caliper to determine the modulus value, the following measurements were taken:
  • measurement of the tip diameter and measurement of the root diameter in the case of gears with an even number of teeth,
  • measurement of the distance AB from the tip of the tooth to the bottom of the opposite root in the case of gears with an odd number of teeth.
These are classic methods of recreating the geometry used in the case of cylindrical gears. In the second case, the following formula was necessary:
A B ¯ = d 2 + h a + d 2 h f
In the case of determining the base pitch using a disc micrometer, span measurements were made across n and n+1 teeth, previously determining the optimal number of teeth using the formula:
n = 1 π z + 2 x cos α 2 z 2 z inv α 2 x tan α + 0.5
It was also checked whether for gears with a relatively small number of teeth and a large module it is possible to measure across n+1 teeth by checking the maximum number of teeth for span measurement from the formula:
n m a x = 1 π d a 2 d b 2 m cos α z inv α 2 x tan α + 0.5
In future research, this will be related to the need to use a workshop microscope, which will change the MPE of the device and consequently the final measurement uncertainty.
The measurement methodology and the metrological steps taken as part of it coincide with those presented by K. Konecki et al. in [12].
Measurements were made in each gear from randomly selected 10 teeth, in a random cross-sectional plane across the width of the gear. Total 8 · 10 · 2 = 160 measurements.

2.4. The model Used Compares the Nominal Gear with the Recreated Gear

In order to estimate the accuracy of reverse engineering, a key geometric parameter of the gear wheel, which is the base pitch, was selected to determine the correctness of the geometry reconstruction. This position of the involute, which unwinds from the base circle on which the tested base pitch is measured, is crucial for the quality of reproducing the geometry of the gear wheel. The original values of this parameter were compared with its recreated values. Recreating the value of the shift factor, which plays an equally important role in the context of tooth thickness, is left for a future article.

2.5. The Model Used for the Statistical Analysis of The results

The following models were used in the work to estimate the tolerance field of the basic pitch.
Model 0 is implicitly the nominal dimension and its deviations (original tolerance zone).
Model I assumes a normal symmetric distribution. The ISO estimator will be here ± 3 σ , also called estimator 6 σ . As NS is used Average:
p b d e s i g n = X ¯ = 1 N · i = 1 N X i
Model II assumes a symmetrical normal distribution. The ISO estimator here will be PCR. As NS is used Average. The estimated tolerance zone is calculated according to the formula:
T d e s i g = 6 · σ · C p
where the C p coefficient is equal to 1.33.
Model III is the Disper parameter calculated according to the formula:
D i s p e r = 2 · max M a x A v e r a g e , A v e r a g e M i n
The half Disper is the upper and lower deviations of the base pitch. Average is used as NS.
Model IV i Model V are the Range parameters calculated according to the formula:
R = M a x M i n
Half-Range is the upper and lower deviation of the base pitch. As NS it is used for models: Median:
p b d e s i g n = M e d = X i m a x + X i m i n 2
and Average.
From Model VI to Model XI it is the adoption of a symmetric S-F distribution. The ISO estimator will be ± t N ,   α σ r . The confidence level was assumed appropriately for the models: α = 99 % , α = 95 % , α = 90 % , α = 80 % , α = 60 % and α = 40 % . Average is used as NS.
The standard deviation was calculated in various ways. When a symmetric Gaussian distribution was adopted, the standard deviation for the sample was calculated according to the formula:
σ = 1 N · i = 1 N X i X ¯ 2
When a symmetric S-F distribution was adopted, the standard deviation of the sample was calculated according to the formulas:
σ t S = σ r · t N ,    α
S = 1 N 1 · i = 1 N X i X ¯ 2
σ r = S N
The value of the Student's t-coefficient depends on the adopted confidence level.
The data necessary to present the results were calculated according to the following formulas:
T t r u e = U S L L S L
E s t i m a t o r   R a t i o = T d e s i g T t r u e · 100 %
The measurement uncertainty of the disc micrometer should be added to each estimated zone; however, in order to determine the value of the base pitch, measurements were made twice: first across n+1 teeth and then across n teeth. Each of these measurements separately is burdened with the MPE of the micrometer. The uncertainty of a micrometer measurement is the geometric sum of such components [13] as:
  • micrometer screw and nut error F = ± 3   μ m ,
  • flatness error of the anvil and spindle T p = ± 0.9   μ m ,
  • parallelism error of the anvil and spindle measuring surfaces T r = 2   μ m ,
  • perpendicularity error of the spindle measuring surfaces relative to its axis T v = 1   μ m ,
  • error of the lower measurement range f A = 2   μ m ,
  • bow deformation error for P = 10   N 2   μ m ,
  • error related to the elastic deformation of the measured object 1   μ m ,
  • error related to immobilizing the spindle with a clamp 2   μ m ,
  • roughness error of the measurement surfaces 0.08   μ m ,
  • parallax error 1   μ m ,
  • temperature error (for steel and T = ± 5 ° C negligible),
  • micrometer setting error 1   μ m .
In this study, it was assumed that the limit error for a single measurement with a disc micrometer to be equal to δ = ± 4   μ m for the measurement range 25 ÷ 50   m m , δ = ± 4   μ m for the measurement range 50 ÷ 75   m m and δ = ± 5   μ m for the measuring range 75 ÷ 100   m m , for the device used during the tests according to DIN 863 and DIN 863-1, factory marked 40 SM [14], [15].
All guidelines for the course of calculations, based on the above-mentioned principles, were taken from the monograph by Jezierski J. [16], which is a classic bibliographic position for Polish scientists, metrologists and designers.

3. Results and Discussion

3.1. Module as a Random Variable in the Population

During the research, a puzzling fact for the authors was whether the reproduced value of the module should be treated as a random variable. Isn't it the case that the tested modules are a sample from the population of possible to reproduce the values of the nominal modules? The authors tried to answer the question of the relationship between the original gear wheel with a non-standard module and the reconstructed gear wheel with a different, but still non-standard module, but very similar in value. It is important that we are talking about nominal geometries here without taking into account real and observed dimensions. We are only talking about the tolerable dimensions according to ISO and the generally tolerable dimensions included in the technical drawing of the original gear wheel and the recreated gear wheel. The authors wanted to answer whether it is necessary at all (and if so to what extent) to get as close as possible to the original value of the nominal gear module when reconstructing the geometry of this gear.
It was found that for two wheels to be the same, they must have identical tooth thickness at any diameter. The base diameter of both wheels must be the same. Taking the above into account, after several transformations, an equation with two unknowns was finally obtained, allowing the selection of a pair of parameters module-shift factor, which can be used alternatively to the originally assumed/calculated ones. The diameter of the base is treated here as a known value.
First, the formula s = π 2 + 2 · x · tan α · m is inserted into the expression s x = d x · s d + i n v α i n v α x :
s x = d x · π 2 + 2 · x · tan α · m z · m + i n v   α i n v   α x
Then two steps were performed: Inserted into the above i n v α = tan α α ¯ and i n v α x = tan α x α ¯ x and then α = cos 1 d b z · m and α x = cos 1 d b d x :
s x = d x · π 2 + 2 · x · tan cos 1 d b z · m z + tan cos 1 d b z · m cos 1 d b z · m tan cos 1 d b d x cos 1 d b d x
After transformation and simplification, the following equation was obtained, called the Equation of Identical Gears (EIG) [17], [18]:
x = z · s x d x π 2 · z tan cos 1 d b z · m cos 1 d b z · m + tan cos 1 d b d x cos 1 d b d x 2 · tan cos 1 d b z · m
This proves that you should not worry about the fact that the recreated value of the module may most likely be completely different than the original value, because in the context of the correctness of the geometry reconstruction process, it has no mathematical significance. Therefore, the module cannot be used as a parameter to compare the original wheel with the recreated wheel. An explanation of this aspect is presented in Table 2, which shows that the use of alternative parameters does not change either the base pitch or the tooth thickness measured along the arc of any chosen circle.
Table 2. Alternative sets of geometric parameters of the same gear (example data).
Table 2. Alternative sets of geometric parameters of the same gear (example data).
m α y x c* pb [mm] dx [mm] sx [mm]
1.98 18.344° 1.028 –0.0876 0.200 5.904 113 1.485
1.97 17.445° 1.034 +0.0507 0.165
1.96 16.490° 1.032 +0.1981 0.180
1.95 15.466° 1.020 +0.3576 0.216
An interesting issue worth paying attention to in future research is the fact that, depending on the set (Table 2), the same gear is shifted negatively ( x < 0 ), and sometimes positively ( x > 0 ). Therefore, if there was a need for a direct tabular comparison of the basic geometric parameters of the original wheel (which would constitute the first row of the table) with the parameters of the recreated wheel (which would constitute the third row of the table), one more middle row of the table should be added with the values of the parameters of the recreated wheel converted thanks to EIG. , bringing this wheel to a common module with the original wheel. Table 3 explains this.
The situation is different from a technological perspective if the wheel being recreated is to be manufactured using hobbing processing. Previously, two Greek scientists, Christos Spitas and Vasilios Spitas [19], [20] showed that in such a case the fillet changes slightly, thus causing a change in the thickness of the tooth at the base. This, in turn, directly affects the bending strength of the tooth.
However, the question must be answered: to what extent can the module be changed in relation to the nominal tested gear so that the reconstructed geometry makes mathematical sense and practical usefulness in terms of the values of alternative geometric parameters? Such a study will be necessary due to geometric and technological limitations in converting a given set of geometric parameters into an alternative set. Figure 3 shows two graphs for an example spur gear with straight teeth and parameters: z = 18 , m = 5 , x = 0 , y = 1 , c * = 0.2 , α = 20 ° .
Figure 3. Preliminary study of the properties of EIG for an example gear.
Figure 3. Preliminary study of the properties of EIG for an example gear.
Preprints 104939 g003
It is clearly visible from the graph (a) that for module 4.70 (i.e. m = 0.30   m m ) a wheel geometry is achieved, which is impossible from the design point of view (for example, due to the negative whole depth factor y = 3.604 ). This variant is close to the left asymptote of the graph, for which α = 0 ° . However, on the right side, where the graph has no limits on the argument domain, theoretically the value of the module can be increased indefinitely. However, this graph is not coupled with the limit value of the shift factor at which tooth undercutting occurs at the base. With reference to the previous conclusion, each set of geometric parameters of the gear wheel has its alternative version in an infinite number, but nevertheless within a range limited by the value of the limit shift factor and of the limit whole depth factor. This requires further research in the future. However, during reverse engineering, the authors do not expect to obtain a recreated modulus value that differs significantly from the nominal value.
The equation discussed here (EIG) can be adapted to gear tests using the finite element method in the strength simulation of gear drives. [21]. Due to this, it will be possible to study the impact of minor manipulation of the modulus value on the strength properties of the gear tooth (this is directly related to the change in the shape of the fillet at the base of the tooth).

3.2. True Deviations and the IT Grade of the Tested Gear Wheel

Figure 4 shows a diagram for selecting the gear accuracy class of the produced depending on the industry and the application area. It is expected that each of the wheels tested is made in grades 5 to 10.
The tested gears could have been manufactured several decades ago. It is therefore impossible to reconstruct the information according to which standard the wheel designer originally calculated the values of permissible deviations. When it comes to the availability of manufacturing deviations in the drawings of the tested wheels, the situation is different. Some drawings clearly show the selected class and according to which standard, while others lack such information instead of providing deviations of selected parameters. To have an increased level of confidence in the authentic accuracy class of the tested wheel, an intermediate solution was chosen and the values of permissible deviations for class 5 were calculated according to the standard [23]. For the remaining classes, conversion factors consistent with the standard were used. The values of the permissible deviations expressed in [mm] have been rounded to the third decimal place. Table 4 shows the formulas for calculating the permissible error values according to DIN 3961-1953 and the conversion factors.
Table 5 shows the exact values of the conversion factors for the deviations of various errors according to DIN (TGL).

3.3. Measurement Results in Statistical Analysis

The manufacturing drawings contain various types of information, both with respect to geometric deviations and, for example, the conditions for accepting the wheel from quality control. The authors verified the IT grade of the tested gears, bringing all the information on allowed errors to the DIN 3961 standard. Table 6 presents the drawing information with analysis in relation to the mentioned standard and the qualification of the wheel to the appropriate IT grade. It was needed, among others, to determine the deviation of the base pitch, the value of which is reproduced as part of the research presented in this work.
In the case of Gear No. 1 Mark 16 a calculator [24] was used to check the correctness of the drawing and, indeed, for the seventh grade g = 0.03   m m . Gear No. 4 Mark 2 is the same wheel as Gear No. 3 Mark 2, however, it has significant signs of use, including large oil stains. In the drawing of Gear No. 5 Mark 3, the deviation of the base pitch is given directly. In the drawing of Gear No. 7 Mark 11 does not provide the standard number but only the DIN letter. In the case of Gear No. 8 Mark 6, the measurement value deviation of 0.054 mm is close to the value of 0.047 mm, which corresponds to grade 9. The deviation of tooth thickness of 0.06 mm is as close to the value of 0.050 mm (grade 9) as to the value of 0.070 mm (grade 10). In both cases, the difference is 0.01 mm. However, due to the fact that the deviation of the measured value at the level of 0.054 mm is close of the right side to the normalized value of 0.047 for grade 9 and is far from the value of 0.089 mm for grade 10, the tested wheel should be classified to grade 9 according to DIN 3961.
Figures 5 to 12 show blocks of double graphs created on the basis of measurement data and statistical calculations.
Figure 5. Charts for Gear No. 1 Mark 16: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Figure 5. Charts for Gear No. 1 Mark 16: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Preprints 104939 g005
Figure 6. Charts for Gear No. 2 Mark 13: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Figure 6. Charts for Gear No. 2 Mark 13: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Preprints 104939 g006
Figure 7. Charts for Gear No. 3 Mark 2: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Figure 7. Charts for Gear No. 3 Mark 2: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Preprints 104939 g007
Figure 8. Charts for Gear No. 4 Mark 2: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Figure 8. Charts for Gear No. 4 Mark 2: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Preprints 104939 g008
Figure 9. Charts for Gear No. 5 Mark 3: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Figure 9. Charts for Gear No. 5 Mark 3: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Preprints 104939 g009
Figure 10. Charts for Gear No. 6 Mark 5: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Figure 10. Charts for Gear No. 6 Mark 5: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Preprints 104939 g010
Figure 11. Charts for Gear No. 7 Mark 11: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Figure 11. Charts for Gear No. 7 Mark 11: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Preprints 104939 g011
Figure 12. Charts for Gear No. 8 Mark 6: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Figure 12. Charts for Gear No. 8 Mark 6: a) Estimated Tolerance Areas, b) Estimator Ratio; for Models I to XI.
Preprints 104939 g012

3.4. Synthesis of Results

To begin with, all wheels should be compared in the context of the recreated manufacturing IT grade according to DIN 3961 (Table 13).
So there is 1 gear in grade 6, 2 gears in grade 7, 3 gears in grade 8 and 2 gears in grade 9.
The direction and degree of shift of the graphs relative to the nominal zone were determined. The degree of shift was determined relative to the Average and Median (two parameters that were used as the centre of the measurement data distribution). The results are presented in Table 14. Each parameter that better approximates the nominal dimension is marked in green. The following formulas were used:
δ 1 = N S A v e r a g e N S · 1000  
δ 1 = N S M e d i a n N S · 1000  
In the case of Gear No. 7 Mark 11 the Average and Median had the same value. Overall, the Average was 3 times more effective, while the Median was 4 times more effective. Due to similar results and a relatively small number of tested wheels, it is not possible to select a better parameter that estimates the nominal value of the basic pitch. The direction of shift to the right occurred 6 times. The direction of shift to the left occurred only twice. In subsequent tests, the nature and potential source of wear of the tested gears should be determined. Except in the case of Gear No. 4 Mark 2 (which, by the way, is the same wheel as Gear No. 3 Mark 2, only worn to a different extent and covered with oil stains) and, coincidentally, Gear No. 8 Mark 6 (which performed relatively well compared to the other wheels), all other cases for Average range from 0.261 to 1.141. In the future, we should try to generate average values of the so-called "straightening coefficients". The idea is to reduce the shifted graph to a nominal field. This will be done analytically (by multiplying the centre of the distribution by appropriate coefficients) or graphically (by using asymmetric distributions).
The authors assessed the degree of filling of the estimated tolerance field inside the nominal field. That is, the Estimator Ratio expressed in [%] was assessed. A point allocation matrix was proposed (Table 15).
However, the above matrix does not directly explain important issues such as: how a specific estimator relates to a given IT grade and how narrow the designed tolerance zone would have to be in relation to the nominal zone and, consequently, how expensive it would be to produce a gear with a better IT grade.
One important issue still remains, namely, what IT grade to choose for the reproduced gear. Although the research results generate a general view as to which estimator best approximates the nominal tolerance zone of the base pitch, it is still unknown within which IT grade this is done. The most effective estimator (Model X with S-F distribution and confidence level α = 60 % ) was taken and this issue was subjected to a preliminary analysis. Table 16 presents the analysis report. These results do not take into account the analysis of the IT grade of the tested wheel. For now, the authors are only preliminary checking what the applied Estimator would cause in the context of the nominal IT grade.
Some estimated IT grades are defined by a range (e.g., 7 to 8 for Gear No. 1 Mark 16), because the estimated deviation is as close to the values representing the given IT grades in the range. In general, the results are quite satisfactory from the point of view of the designer recreating the gear because the increase in the cost of producing a gear wheel using two higher IT grades occurred only in the case of one of the eight tested wheels.

4. Conclusions

The most important conclusion of this study is that the nominal module does not have to be treated as a random variable in the population. Each gear wheel, within a certain range limited by the limits of x and y values, has an infinite number of alternative sets of geometric parameters within which the designed gear remains completely identical to the gear wheel with original parameters. The authors determined EIG (Equation of Identical Gears), which allows one to obtain an alternative value of the shift factor for a new and different nominal module for the tested gear. This equation will play a significant role in future research in the context of the strength of the tooth at the base to bending loads.
The tested gears were made in 6, 7, 8 and 9 IT grades according to DIN 3961. It is important that the individual standards that determine the size of permissible gear toothing errors differ from each other. They have formulas and expressions for calculating these errors, giving divergent results. You should also remember about industry standards that are determined individually within a given company. Therefore, it is important to indicate one specific standard for each wheel tested. The IT grade of tested gears should be reduced to one common set, based on the selected standard.
The estimator with the least effectiveness in this study was the one used in Model II, which estimated the tolerance field using PCR (with a C p coefficient of 1.33). This model did not show satisfactory results in any of the gears tested. However, we expect that it could be useful for testing new and unworn gears.
The most effective estimator turned out to be the one presented in Model X, using the symmetric Student-Fisher distribution, which turned out to estimate the nominal tolerance zone by 100 or close to 100 (at least in one case of gear 53) percent. This allowed for: (i) in two cases, maintain the original IT grade of the tested wheel, (ii) leave the IT grade the same or increase it by one step in the case of two gears, (iii) increase the IT grade by one step for only three gears, and (iv) increase the IT grade of just one gear by two steps. Of course, we are only talking about the Estimator Ratio expressed in [%]. However, in the future, the possibility of restoring the nominal IT grade should be taken into account, as well as the method of "straightening" the graph of the estimated tolerance zone in relation to the nominal tolerance zone, as significant shifts in the graphs could be noticed in the tests.
The next stage will be to increase the number of tested wheels in such a way that each wheel has at least 2, optimally with a different degree of wear. These wheels will be tested on two CMMs [25]: Wenzel WGT 600 and Wenzel LH 65.
The second, completely separate stage of the research will be to determine the shift factor.

Author Contributions

Conceptualization, K.K. (Karol Konecki); methodology, K.K. (Karol Konecki) and D.W. (Dominik Wojtkowiak); software, K.K. (Karol Konecki); validation, K.K. (Karol Konecki) and D.W. (Dominik Wojtkowiak); formal analysis, K.K. (Karol Konecki); investigation, K.K. (Karol Konecki) and D.W. (Dominik Wojtkowiak); resources, K.T. (Krzysztof Talaśka); data curation, K.K. (Karol Konecki) and K.T. (Krzysztof Talaśka); writing—original draft preparation, K.K. (Karol Konecki); writing—review and editing, D.W. (Dominik Wojtkowiak) and K.T. (Krzysztof Talaśka); visualization, K.K. (Karol Konecki); supervision, D.W. (Dominik Wojtkowiak) and K.T. (Krzysztof Talaśka); project administration, K.K. (Karol Konecki); funding acquisition, K.T. (Krzysztof Talaśka). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

In this section, you can acknowledge any support given which is not covered by the author contribution or funding sections. This may include administrative and technical support, or donations in kind (e.g., materials used for experiments).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Krogstle, L.; Walter, M.S.J.; Wartzack, S.; Martinsen, K. Towards a more Comprehensive Understanding of Tolerance Engineering Research Importance. Procedia CIRP; Volume 27; (2015); ISSN 2212-8271; pp. 29-34. [CrossRef]
  2. Konecki, K.; Kołodziej, A.; Wojtkowiak, D.; Talaśka, K. Issues related to the reconstruction of the geometry of the DP6 inch cylindrical gear in the context of reverse engineering using conventional techniques and measuring instruments. AIP Conf. Proc. (2023), 2976, 020008 . [CrossRef]
  3. Sáenz-Nuño, M.A.; Lorente-Pedreille, R. ISO Tolerance specification in reverse engineering. Procedia Manufacturing; Volume 13; (2017); ISSN 2351-9789; pp. 472-479. [CrossRef]
  4. Kaisarlis, G.J. A Systematic Approach for Geometrical and Dimensional Tolerancing in Reverse Engineering. In: Reverse Engineering - Recent Advances and Applications; Telea, A.C. (Ed.); IntechOpen: 51000 Rijeka, Croatia, (2012); ISBN: 978-953-51-0158-1, EBOOK (PDF) ISBN: 978-953-51-5602-4; pp. 133-160. [CrossRef]
  5. Kaisarlis, G.J.; Diplaris, S.C.; Sfantsikopoulos, M.M. Position tolerancing in reverse engineering: The fixed fastener case. Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture. (2007); 221(3):457-464. [CrossRef]
  6. Kaisarlis, G. J.; Diplaris, S. C.; Sfantsikopoulos, M. M. Geometrical position tolerance assignment in reverse engineering, International Journal of Computer Integrated Manufacturing. (2008); 21:1, 89-96. [CrossRef]
  7. Kaisarlis, G.J.; Diplaris, S.C.; Sfantsikopoulos, M.M. A Knowledge-Based System for Tolerance Allocation in Reverse Engineering. In: Hayhurst, D.R., et al. Proceedings of the 33rd International MATADOR Conference., (2000); Springer, London. [CrossRef]
  8. Geng, Z.; Bidanda, B. Tolerance estimation and metrology for reverse engineering based remanufacturing systems, International Journal of Production Research. (2022); 60:9, 2802-2815. [CrossRef]
  9. Geng, Z.; Sabbaghi, A.; Bidanda, B. Reconstructing original design: Process planning for reverse engineering, IISE Transactions. (2023); 55:5, 509-522. [CrossRef]
  10. Standard ASME. 2009. ASME Y14.5 Dimensioning and Tolerancing.
  11. Standard ISO. 2019. ISO/TC 213. Dimensional and geometrical product specifications and verification.
  12. Konecki, K.; Wojtkowiak, D.; Talaśka, K.; Kołodziej, A.; Domek, G. Issues related to an attempt to recreate the geometry of a non-standard spur gear. IOP Conf. Ser. Mater. Sci. Eng. (2021), 1199, 012105. [CrossRef]
  13. Standard PN-82/M-53200/A1. Narzędzia pomiarowe. Przyrządy mikrometryczne. Wymagania (informal translation: Measuring tools. Micrometer instruments. Requirements).
  14. Mahr Group. Brochure: Katalog zbiorczy | Technika pomiarów produkcyjnych (informal translation: Collective catalog | Production measurement techniques), POL05 V1, 3758450-01.10.2005, Printed in Germany, pp. PG-17.
  15. Mahr Group Website. Available online: https://metrology.mahr.com/fileadmin/catalogs/OnlineCatalogs/Mahr-Metrology/index.php?catalog=Mahr-Metrology-PL&lang=pl (accessed on 19 January 2024).
  16. Jezierski, J. Analiza tolerancji i niedokładności pomiarów w budowie maszyn (informal translation: Analysis of tolerances and measurement inaccuracies in machine construction), 3rd ed.; Publishing house: Wydawnictwa Naukowo-Techniczne WNT, Warszawa, Poland, (1994); ISBN: 83-204-2906-4.
  17. Konecki, K.; Kołodziej, A.; Pytliński, R. Equation of Identity for Two Spur Gears with Straight Teeth Having Two Different Alternative Sets of Geometric Parameters. AIP Conf. Proc. (2023), 2976, 030003 . [CrossRef]
  18. Konecki, K.; Kołodziej, A.; Wojtkowiak, D.; Talaśka, K. Simulation of different geometric variants of spur gears in the context of their complete identity. AIP Conf. Proc. (2023), 2976, 030018 . [CrossRef]
  19. Spitas, C.; Spitas, V. Can non-standard involute gears of different modules mesh? Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science. (2006); Volume 220; Issue 8; pp. 1305-1313. [CrossRef]
  20. Spitas, C.; Spitas, V. Generating standard 20° involute pinions with increased fillet strength by using 25° rack cutters with non-standard module. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science. (2006); Volume 220; Issue 8; pp. 1297-1304. [CrossRef]
  21. Wojtkowiak, D.; Konecki, K. Stress analysis of the zero-backlash gear transmission. AIP Conf. Proc. (2023), 2976, 030011 . [CrossRef]
  22. Ochęduszko, K. Koła zębate t. 3 Sprawdzanie (informal translation: Gears, Volume 3 - Checking), 8th ed.; Publishing house: Wydawnictwa Naukowo-Techniczne WNT, Warszawa, Poland, (2012); ISBN: 9788363623067, EAN: 9788363623067; pp. 238.
  23. Standard DIN 3961:1953-09. Tolerances for Cylindrical Gear Teeth; Bases.
  24. DIN 3967 | Calculator (mechanicalcheck.com). Available online: https://www.mechanicalcheck.com/en/din3967/calc/ (accessed on 19 January 2024).
  25. Broszura, Instytut Inżynierii Mechanicznej Wydziału Politechnicznego „-Pracownia Inżynierii Mechanicznej-„ Uniwersytet Kaliski im. Prezydenta Stanisława Wojciechowskiego (University of Kalisz, Poland), Available online: https://uniwersytetkaliski.edu.pl/wp-content/uploads/2024/02/iim-12.02.2024.pdf (accessed on 18 March 2024).
Figure 1. Illustrating the idea of RE for a spur gear [2].
Figure 1. Illustrating the idea of RE for a spur gear [2].
Preprints 104939 g001
Figure 2. Tested gear wheels: a) Gear No. 1 Mark 16, b) Gear No. 2 Mark 13, c) Gear No. 3 Mark 2, d) Gear No. 4 Mark 2, e) Gear No. 5 Mark 3, f) Gear No. 6 Mark 5, g) Gear No. 7 Mark 11, h) Gear No. 8 Mark 6.
Figure 2. Tested gear wheels: a) Gear No. 1 Mark 16, b) Gear No. 2 Mark 13, c) Gear No. 3 Mark 2, d) Gear No. 4 Mark 2, e) Gear No. 5 Mark 3, f) Gear No. 6 Mark 5, g) Gear No. 7 Mark 11, h) Gear No. 8 Mark 6.
Preprints 104939 g002aPreprints 104939 g002b
Figure 4. Selection of the accuracy class according to the purpose. The basic class is class 5. Classes used in the aviation industry: from 5 to 10 [22].
Figure 4. Selection of the accuracy class according to the purpose. The basic class is class 5. Classes used in the aviation industry: from 5 to 10 [22].
Preprints 104939 g004
Table 1. Basic geometric parameters of the gears tested.
Table 1. Basic geometric parameters of the gears tested.
No. Mark z m α y x c*
1 16 13 3.020 21.018° 0.9924 -0.0421 0.1675
2 13 18 4.980 26.325° 0.7950 +0.0695 0.3830
3 2 42 4.029 21.104° 1.0189 +0.4184 0.1463
4 2 42 4.029 21.104° 1.0189 +0.4184 0.1463
5 3 30 4.410 25.572° 0.9012 -0.5423 0.5196
6 5 29 1.020 22.888° 0.9283 -0.5263 0.3002
7 11 17 4.980 19.358° 1.0065 -0.0641 0.1958
8 6 18 4.980 26.325° 0.8632 +0.0695 0.2525
Table 3. Comparison of the recreated wheel with the original wheel (example data).
Table 3. Comparison of the recreated wheel with the original wheel (example data).
m α y x c*
The recreated gear 1.97 17.538° 1.0339 +0.0509 0.1658
The identical gear (EIG) 1.98 18.432° 1.0278 –0.0874 0.1666
The nominal gear 1.98 18.344° 1.0280 –0.0876 0.2
Table 4. A list of formulas [22] to calculate the permissible error deviations of individual tooth elements, according to DIN 3961-1953 for IT grade 5 [23]. Values in [μm].
Table 4. A list of formulas [22] to calculate the permissible error deviations of individual tooth elements, according to DIN 3961-1953 for IT grade 5 [23]. Values in [μm].
Type of error φ Module   range   m   [ m m ]
do 0.6 od 0.6 do 1.6 od 1.6 do 4 od 4 do 10
e   t p   o   t 0 1.4 / 1.6 3 + 0.3 m + 0.2 d (19)
r z 1.4 / 1.6 2.5 + 0.25 m + 0.15 d (20)
Σ t 0 1.4 / 1.6 10 + m + 0.8 d (21)
g 1.4 1.4 + 1.85 m + 2.06 d 4 (22) 1.4 + 1.85 m + 2.18 d 4 (23) 1.4 + 1.85 m + 2.30 d 4 (24) 1.4 + 1.85 m + 2.43 d 4 (25)
r 1.4 2 + 2.57 m + 2.89 d 4 (26) 2 + 2.57 m + 3.07 d 4 (27) 2 + 2.57 m + 3.25 d 4 (28) 2 + 2.57 m + 3.45 d 4 (29)
M 1.4 1.32 + 1.7 m + 1.95 d 4 (30) 1.32 + 1.7 m + 2.06 d 4 (31) 1.32 + 1.7 m + 2.18 d 4 (32) 1.32 + 1.7 m + 2.3 d 4 (33)
Table 5. Conversion factors for the deviations of various toothing and engagement errors, according to DIN (TGL), depending on the IT grade. Only grade 5 to 10 is shown.
Table 5. Conversion factors for the deviations of various toothing and engagement errors, according to DIN (TGL), depending on the IT grade. Only grade 5 to 10 is shown.
Type of error IT grade
5 6 7 8 9 10
e   t p   o   t 0   r z   Σ t 0 1 1.4 1.96 2.74 3.84 6.15
g   r   M 1 1.4 1.96 2.74 3.84 5.36
Table 6. Qualification of the tested gears to the appropriate IT grade according to DIN 3961.
Table 6. Qualification of the tested gears to the appropriate IT grade according to DIN 3961.
Gear Drawing information Calculated errors in 5 to 10 IT grades according to DIN 3961 Qualified for IT grade
No. Mark IT grade information According to the standard Permissible errors [mm] 5 6 7 8 9 10
1 16 7 DIN 3967 ΔM = 0.028
Δg = 0.030
Δtp = ?
0.010
0.010
0.005
0.014
0.014
0.007
0.020
0.020
0.010
0.027
0.028
0.014
0.038
0.040
0.020
0.054
0.055
0.032
8
2 13 ΔM = 0.022
Δg = 0.025
Δtp = ?
0.012
0.013
0.006
0.017
0.018
0.009
0.024
0.025
0.013
0.033
0.036
0.018
0.047
0.050
0.025
0.065
0.070
0.039
7
3 2 ΔM = 0.035
Δg = 0.037
Δtp = ?
0.013
0.013
0.007
0.018
0.019
0.010
0.025
0.026
0.013
0.034
0.037
0.019
0.048
0.051
0.026
0.067
0.072
0.042
8
4 2
5 3 Δtp = 0.010 0.007 0.009 0.013 0.018 0.025 0.040 6
6 5 7 DIN 3962 Δg = 0.030
Δtp = ?
0.008
0.004
0.012
0.006
0.016
0.009
0.023
0.012
0.032
0.017
0.045
0.027
9
7 11 9 DIN ΔM = 0.024
Δtp = ?
0.012
0.006
0.017
0.009
0.024
0.012
0.033
0.017
0.046
0.024
0.065
0.039
7
8 6 I industry standard ΔM = 0.054
Δg = 0.060
Δtp = ?
0.012
0.013
0.006
0.017
0.018
0.009
0.024
0.025
0.013
0.033
0.036
0.018
0.047
0.050
0.025
0.089
0.070
0.039
9
Table 13. Recreated IT grades according to DIN 3961 for the tested gears.
Table 13. Recreated IT grades according to DIN 3961 for the tested gears.
Gear designation IT grade
5 6 7 8 9 10
No. 1 Mark 16
No. 2 Mark 13
No. 3 Mark 2
No. 4 Mark 2
No. 5 Mark 3
No. 6 Mark 5
No. 7 Mark 11
No. 8 Mark 6
Table 14. Direction and degree of shift of the graphs of the estimated tolerance zones.
Table 14. Direction and degree of shift of the graphs of the estimated tolerance zones.
Gear designation Direction Degree of shift δ1 and δ2 [‰]
Relative to
Average
Relative to
Median
No. 1 Mark 16 1.028 1.254
No. 2 Mark 13 0.918 0.704
No. 3 Mark 2 0.595 0.722
No. 4 Mark 2 5.676 4.109
No. 5 Mark 3 0.261 0.181
No. 6 Mark 5 0.464 0.972
No. 7 Mark 11 1.141 1.141
No. 8 Mark 6 0.171 0.009
Table 15. Estimator Ratio evaluation matrix, expressed in [%].
Table 15. Estimator Ratio evaluation matrix, expressed in [%].
IT grade wg DIN 3961
Model
8 7 8 8 6 9 7 9 Suma
pkt.
Occurrence in all measured gears? Verdict
XI 75 YES Second best
X 110 YES Best
IX 75
VIII 70
VII 60
VI 45
V 25
IV 30
III 50
II 0
I 30
Gear No.→ 1 2 3 4 5 6 7 8
Mark → 16 13 2 2 3 5 11 6
Explanations: Estimator Ratio > 100 %
①→ 30 pkt. Estimator Ratio < 100 %, 1st place
②→ 20 pkt. Estimator Ratio < 100 %, 2nd place
③→ 10 pkt. Estimator Ratio < 100 %, 3rd place
5 pkt. Estimator Ratio < 100 %
Table 16. Report on the preliminary analysis of the applicability of Model X.
Table 16. Report on the preliminary analysis of the applicability of Model X.
Gear True
IT Grade
Estimator Ratio
in Model X [%]
t p   according to DIN 3961 [mm] T t r u e
[mm]
T d e s i g
[mm]
t p
designed [mm]
Designed
IT Grade
No 1. Mark 16 8 81 ±0.014 0.028 0.023 ±0.012 7–8
No 2. Mark 13 7 80 ±0.013 0.026 0.021 ±0.011 6–7
No 3. Mark 2 8 65 ±0.019 0.038 0.025 ±0.013 7
No 4. Mark 2 8 93 ±0.019 0.038 0.035 ±0.018 8
No 5. Mark 3 6 100 ±0.009 0.018 0.018 ±0.009 6
No 6. Mark 5 9 53 ±0.017 0.034 0.018 ±0.009 7
No 7. Mark 11 7 83 ±0.012 0.024 0.020 ±0.010 6
No 8. Mark 6 9 81 ±0.025 0.050 0.041 ±0.021 8
Explanations: An increase of two IT grades (gear would have to be made more precisely)
An increase of one IT grade (gear would have to be made more precisely)
The IT grade remains the same
Remains the same IT grades or increases by one IT grade (more precisely)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated