Next Article in Journal
Cold Welding in Hold Down Points of Space Mechanisms Due to Fretting When Omitting Grease
Next Article in Special Issue
A Study on the Flow Resistance of Fluids Flowing in the Engine Oil-Cooler Chosen
Previous Article in Journal
Effect of Graphene Addition in Cutting Fluids Applied by MQL in End Milling of AISI 1045 Steel
Previous Article in Special Issue
Artificial Aging Experiments of Neat and Contaminated Engine Oil Samples
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impacts of a Profile Failure of the Cycloidal Drive of a Planetary Gear on Transmission Gear

Department of Machine and Product Design, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3, 1111 Budapest, Hungary
Submission received: 20 April 2021 / Revised: 10 July 2021 / Accepted: 14 July 2021 / Published: 20 July 2021
(This article belongs to the Special Issue Automotive Tribology II)

Abstract

:
Recently, cycloidal drives have been increasingly used due to their beneficial features, such as the implementation of a large transmission, efficiency, and high performance density. Production accuracy is inevitable in order to ensure the dynamically proper and smooth operation of the drive. Smaller backlash and minimal transmission fluctuations can only be achieved by improving the production accuracy, reducing the number of production failures, and shrinking the tolerance zone. This research primarily focuses on the investigation of the tolerable production accuracy of small-series and individually manufactured drives. The analysis of load distribution was calculated on the cogs of a planetary gear made with a wire EDM machine. On the other hand, we investigate how production failures affect transmission fluctuations and the backlash of a drive. The novelty of this research is based on the determined analytical equations, which can help engineers to find the right tolerances to a given gear ratio fluctuation.

1. Introduction

The investigation of the impact of production failures on kinematic features was performed by supposing different levels of production accuracy and tolerance zones created on different cog geometries [1,2,3,4]. Based on the findings of this study, we can determine the required minimum accuracy level for a particular cog geometry, which ensures achieving the predefined maximum level of transmission fluctuation [5,6,7,8].
When performing the calculations, geometrical equations describing the connection of cogs were solved by numerical methods using a special software developed for this purpose. The grade of production failures was determined according to particular tolerance zones both in the load distribution and the definition of kinematic features.
Transmission fluctuation caused by production failures was established numerically for different gear geometries. The findings suggest that the transmission changes depending on the rotation angle of the output shafts [5,6,7,8]. The required production accuracy can be estimated for the particular maximum transmission failure in advance by relying on a derived analytical formula based on the maximum transmission fluctuation calculated for a particular module, profile failure, and cog number.
The manuscript’s main sections are as follows:
  • Production failures: The profile failures of the cycloidal planetary gears are discussed.
  • Definition of the impact of the profile failure on the transmission: The effect of the cycloidal planetary gear profile failure on the transmission is discussed.
  • Examined gear geometries: The analyzed cycloidal gear geometries are presented.
  • Results: Our obtained results are presented.
  • Discussion and Conclusions: Summary of the research, with the major results highlighted.

2. Production Failures

In the main components of the drives with cycloidal cogs, diameter failures of the pitch-circle may occur [9,10]. There can be distribution failures, profile failures, diameter failures of driving bores, or hit failures of the rolling circle; however, hit failures of the planetary gear are also possible due to an eccentricity failure.
When cycloidal planetary gears are produced, the following accuracy levels are applicable: IT5 and IT6 and, in the case of individually manufactured small-series drives IT7, IT8, and, in rare cases, IT9. In this research, we examine the connection between profile failure and transmission fluctuation.
During the calculations, profile failure can be generated by shifting the theoretical cog curve without any failures by value “ d ” (failure) into the normal direction of the profile. The interpretation of the profile failure can be seen in the following figure (Figure 1).
The grade of the profile failure has been defined in due consideration of the tolerance zone determined by the diameter of the pitch-circle of the planetary gear and the particular accuracy level in advance. In Table 1, profile failures are shown.

3. Definition of the Impact of the Profile Failure on the Transmission

During the research, the maximum deviation of the transmission from the theoretical figure was determined with different profile failures for several modules and cog numbers. First, the geometrical connections required to define the transmission figure were derived, later we examined how profile failures affect transmissions. Equations describing the connection of gears were solved using numerical methods in a Scilab environment.
Cycloidal drives are high power-density gearboxes consisting of one, but preferably two (3) planetary gears, one (4) ring gear, and one (2) eccentric shaft driving the planetary gears. The numbers between the parentheses represents the components in Figure 2. During the calculations, the initial state were determined by the theoretical curve without any failures. An interpretation of the initial state and a schematic concept of the cycloidal drive are illustrated in Figure 2.
The first step is to define the rotation of the ( 2 ) input shaft, which is necessary for the hit ( α ) . The required rotation for the collision can be seen in Figure 3.
When the shaft rotates, the planetary gear makes an orbital motion, i.e., it does not rotate, but its center point moves along the circle determined by the ( 3 ) eccentric cam (Figure 3). The rotation angle of the ( 2 ) input shaft required for the hit always has to be determined for the first pin located left from the pole point (where the absolute velocity equals zero) applicable to the theoretical curve (without any failures).
The contact of two curves, i.e., the profile curve and the circle (pin curve) has to be examined. In the case of a hit, the center point of the pin will be at a radius distance from the profile curve and thus the profile curve can be shifted by the pin radius toward the normal direction of the profile. Therefore, only the contact between one curve and one point has to be examined. When the pin in question hits, the following connection can be established according to the figure:
r 1 r 2 ( α ) r 3 ( φ ) = 0
where r 1 is the vector pointing at the center point of the pin in question, r 2 ( α ) is the vector pointing at the center point of the planetary gear, α is the rotation angle that is necessary for the eccentric cam collision, r 3 ( φ ) is the vector describing the profile of the planetary gear shifted by the pin radius ( R c s ) into the normal direction, and φ is the parameter of the equation depicting the curve.
Therefore, we can find the figures of α and φ to which the above-mentioned equation applies, i.e., where the profile shifted by the radius R c s crosses the center point of the pin, the next step is to determine the gap between the cog pairs. To determine the gap between the cog pairs, see the following figure (Figure 4):
In the case of a particular cogpair, the following equation can be established to define the gap:
r 1 r 2 ( φ 1 ) r 3 ( R , φ ) = 0
where r 1 is the vector pointing at the center point of the pin in question, r 2 ( φ 1 ) is the vector pointing at the center point of the planetary gear, φ 1 is the angle determining the location of the ( 3 ) eccentric cam, r 3 ( R , φ ) is the vector describing the shifted cycloidal curve, and R is the value by which the interrupted curve is shifted into the normal direction of the profile.
The calculation is performed in the above-mentioned way with the difference that now R is one of the variables (unknown parameters); therefore, the gap is:
s = R R c s = 0
where R c s -is the pin radius.
The next step is to determine the transmission. To do that, it is necessary to determine the rotation angle of the ( 5 ) output shaft ( φ 2 ) depending on the rotation angle of the ( 2 ) input shaft ( φ 1 ) . The definition of the rotation angle can be seen in Figure 5.
During the motion of the planetary gear, the curve shifted by the radius from its profile always crosses the center point of the connected pin.
The equation describing the shifted curve is given in the coordinate system X Y , which is attached to the planetary gear. The vector indicated in this system ( r 3 ) will be converted into the coordinate system ( X , Y ) created for the (4) ring gear. After that, we make it equal with the vector pointing at the center point of the pin ( r 1 ) .
When two pins are connected, the following equation shall be applicable using the signs of the figure:
r 1 = T · r 3 ( φ ) + r 2 ( φ 1 )
T = cos ( φ 2 ) sin ( φ 2 ) - sin ( φ 2 ) cos ( φ 2 )
where r 1 is the vector pointing at the centerpoint of the connected pin, T is the transformation matrix, r 3 ( φ ) ) is the vector containing the parameter equations describing the shifted curve in the coordinate system X Y where φ is the parameter of the equations describing the curve, r 2 ( φ 1 ) is the vector pointing at the centerpoint of the planetary gear, φ 1 is the rotation angle of the ( 2 ) input shaft, and φ 2 ( φ 1 ) is the rotation angle of the planetary gear.
The two variables φ and φ 2 can be determined for the figure φ 1 as indicated in the above-mentioned equation. Knowing the figures φ 2 and φ 1 , the transmission can be calculated as follows (Figure 6):
If the rotation of the planetary gear ( φ 2 ) aligned to the particular rotation value of the ( 3 ) eccentric cam ( φ 1 ) is known, the location of the pole point can be determined based on these values (Figure 6).
The speed of the center point of the planetary gear can be described in two ways:
v 0 = ω 1 · e
v 0 = ω 2 · r g ( φ 1 , φ 2 )
Equations ( 6 ) and ( 7 ) are made equal and arranged:
i = ω 1 ω 2 = r g ( φ 1 , φ 2 ) e
where e is the eccentricity (axial distance) and r g ( φ 1 , φ 2 ) is the distance of the center point of the planetary gear from the pole point.
The above-mentioned equations are solved numerically. The process diagram of the algorithm is shown in the following figure (Figure 7).
Here, α is the rotation angle of the ( 3 ) eccentric cam required for the collosion, N is the number of the connected pins, s i is the gap calculated at the pin i, N ( m i n ( s i ) ) is the number of pins defined for the minimum gap, φ 1 is the rotation angle of the ( 3 ) eccentric cam ( φ 1 = φ 1 φ n ) , φ 2 is the rotation angle of the planetary gear, and Δ φ 1 is the grade of the rotation angle of the ( 3 ) eccentric cam.
We use the data of the examined gear pairs of a particular drive as an example. You can see the parameters of the examined gear pairs in the following Table 2:
As it is shown by the algorithm, the rotation angle of the ( 2 ) input shaft needed for the hit was determined first. It can be seen that pin number 2 hits first (Figure 8). On the left-hand side of the following image, the gap belonging to certain pins are marked with a circle. The gap size is proportional to the circle diameter. On the right-hand side of the picture, the column diagram also shows the gap size.
The transmission was determined supposing the rotation of the ( 3 ) eccentric cam in 360 [ ° ] ( φ 1 ) , in grades of 2 [ ° ] ( Δ φ 1 ) . Rotation angle figures of the planetary gear aligned to particular rotation angles of the ( 3 ) eccentric cam are shown in the following figure (Figure 9).
Calculation results are marked with dotted lines, and cases without any failures are marked with continuous lines in the picture. The latter one is a straight line because the transmission is constant. It can be seen in the enlarged segment that they are not identical, which is due to the profile failure.
Knowing the rotation angle of the planetary gear as-signed to the particular rotation angle of the ( 3 ) eccentric cam, the transmission change can be established using the equation shown above (8) or through derivation since the following equation applies to both rotation angles.
As shown by the algorithm, the rotation angle of the ( 2 ) input shaft needed for the hit was determined first. It can be seen that pin number 2 hits first (Figure 9). On the left-hand side of the following image, the gap belonging to certain pins are marked with a circle. The gap size is proportional to the circle diameter. On the right-hand side of the picture, the column diagram also shows the gap size.
i = d φ 1 d φ 2
This indicates that the transmission is the rise of the curve marked with dotted lines compared to the vertical axis. In this case, this means a numerical derivation according to the following equation:
i n d φ 1 ( n + 1 ) d φ 1 ( n 1 ) d φ 2 ( n + 1 ) d φ 2 ( n 1 )
The following diagram shows the transmission change, which was determined in two different ways. The transmission change calculated using Equation ( 9 ) is marked with i 1 , and the change in transmission determined using Equation ( 10 ) is marked with i 2 . The transmission was determined supposing the rotation of the ( 3 ) eccentric cam in 360 [ ° ] ( φ 1 ) , in grades of 2 [ ° ] ( Δ φ 1 ) . Rotation angle figures of the planetary gear aligned to particular rotation angles of the ( 3 ) eccentric cam are shown in the following figure (Figure 10).
It is clearly seen that transmission is not constant, it fluctuates around the value without any failures ( i = 9 ) . Furthermore, the curve describing the transmission change is not continuous. The reason for that is that there is only one cog-pair connection because of the failure, and the connected cog-pairs vary during the movement of the planetary gear. Where the curve is interrupted, cog-pairs leave the connection, and the next cog-pair will determine the motion. Therefore, each segment of the curve belongs to different cog-pairs.
Transmission fluctuations determined in two different ways have brought approximately the same results. The only difference is that the transmission ( i 2 ) determined by Equation (10) does not bring the real result at the points of the function φ 2 = f ( φ 1 ) where the derived figure can not be construed. These points are the maximum points of the above diagram; therefore, it is advisable to use Equation (11) in the future.
The percentile deviation of the real transmission from the original one can be determined by the following equation:
Δ i = i i d e a l i r e a l i i d e a l · 100
The deviation calculated by using Equation (10) is shown in the next figure (Figure 11).
It is important to investigate how further cog parameters affect transmission fluctuation. The percentile transmission fluctuation calculated by using Equation (10) is shown in the following examples. The impact of the increase of the profile shift coefficient is shown in the next figure (Figure 12).
By increasing the profile shift coefficient, the maximum value of transmission fluctuation can be raised. The impact of profile failures on transmission fluctuation is shown in the next diagram (Figure 13).
Reduction of the number of profile failures brings the expected result, meaning the transmission fluctuation decreases. The fluctuation decreases at the rate of failure reduction (approximately), i.e., if the failure is reduced by 90 % (it decreases from 0.1 [ m m ] to 0.01 [ m m ] ), the maximum value of fluctuation also decreases to its tenth.
Curves describing the fluctuation show that the curve will develop as many sections ( z 2 ) during the complete rotation that suit the cog number of the (4) ring gear. These sections coincide with each other. For this reason, the maximum fluctuation does not need to be determined in consideration of the complete rotation, it is enough to take only one cog-pair into account.
In the case of the same profile failures, the transmission of any change in the cog number modifies the transmission as follows (Figure 14).

4. Examined Gear Geometries

The diameter of the pitch-circle is determined by a module and the cog number. Based on the indicated cog number and module, this is 100 different gear sizes; however, in the case of m = 5 [ m m ] and z 1 = 117 [ ] , the diameter of the pitch-circle is larger than 500 [ m m ] (this is a rare application in practice), and thus calculations are only performed for 95 different gear geometries.
The accuracy levels used to define profile failures are as follows: IT5, IT6, IT7, IT8, and IT9. You can see the examined gear parameter in Table 3.

5. Results

All the given analytical results can be used in the design process. The resulting findings of the research are shown in the following Figure 15, Figure 16, Figure 17 and Figure 18.
Results belonging to the value m = 1 [ m m ] :
Equation for the fitted surface:
Δ i = 6.842 · d + 3351.9 · d · z 1 1.5
where d is the profile error, and z 1 is pinion teeth.
Results belonging to the value m = 2 [ m m ] :
Equation for the fitted surface:
Δ i = 3.464 · d + 1720.82 · d · z 1 1.5
Results belonging to the value m = 3 [ m m ] :
Equation for the fitted surface:
Δ i = 2.293 · d + 1156.47 · d · z 1 1.5
Results belonging to the value m = 5 [ m m ] :
Equation for the fitted surface:
Δ i = 1.838 · d + 680.33 · d · z 1 1.5
The results and the equation of the fitted surface also suggest that the maximum value of the transmission failure changes linearly depending on profile failures and along a hyperbole in terms of the cog number. Knowing the approximate connections (equations), the minimum production accuracy can be estimated if we know the maximum transmission fluctuation as well.

6. Discussion and Conclusions

The results and the equation of the fitted surface also suggest that the maximum value of the transmission failure changes linearly depending on profile failures and along a hyperbole in terms of the cog number. The minimum production accuracy can be estimated if we know the maximum transmission fluctuation. Using the method of the smallest squares, the analytical equation of the transmission fluctuation can be created when fitting a surface to the points representing the results. The percentile deviation of the real transmission from the original one can be determined with the derived equation.

Funding

The research reported in this paper and carried out at BME was supported by the NRDI Fund (TKP2020 NC, Grant No. BME-NCS) based on the charter of bolster issued by the NRDI Office under the auspices of the Ministry for Innovation and Technology.

Acknowledgments

Hereby, I would like to express my thanks to the engineering student Dávid Roboz who made a considerable contribution to my research by doing a great job.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cyclo Drive 6000 Reducer. Available online: https://emeia.sumitomodrive.com/en-de/product/cyclo-drive-6000-reducer (accessed on 19 April 2021).
  2. Fine Cyclo Zero Backlash Speed Reducers. Available online: https://emeia.sumitomodrive.com/en-de/fine-cyclo-brand (accessed on 19 April 2021).
  3. Precision Reduction Gear RV. Available online: http://www.nabtescomotioncontrol.com/pdfs/af-series-brochure-2016.pdf (accessed on 19 April 2021).
  4. Roladrive Speed Reducer. Available online: http://www.kotec.com.tw/PDF/KOTECCatalogueCHT.pdf (accessed on 19 April 2021).
  5. Chmurawa, M.; Lokiec, A. Distribution of loads in cycloidal planetary gear (CYCLO) including modification of equidistant. In Proceedings of the 16th European Mechanical Dynamics User Conference, Berchtesgden, Germany, 14–15 November 2001. [Google Scholar]
  6. Thube, S.V.; Bobak, T.R. Dynamic Analysis of a Cycloidal Gearbox Using Finite Element Method; AGMA Technical Paper; AGMA: Alexandria, VA, USA,, 2012. [Google Scholar]
  7. Kumar, N. Investigation of Drive-Train Dynamics of Mechanical Transmission Incorporating Cycloidal Drives. Ph.D. Thesis, Queensland University of Technology, Brisbane, Australia, 2015; pp. 261–268. [Google Scholar]
  8. Borislavov, B.; Borisov, I.; Panchev, V. Design of a Planetary-Cyclo-Drive Speed Reducer Cycloid Stage, Geometry, Element Analyses. Ph.D. Thesis, Linnaeus University, Växjü, Sweden, 2012; pp. 37–44. [Google Scholar]
  9. Tsai, S.J.; Huang, C.H.; Yeh, H.-Y.; Huang, W.J. Loaded Tooth Contact Analysis of Cycloid Planetary Gear Drives. In Proceedings of the The 14th IFToMM World Congress, Taipei, Taiwan, 25–30 October 2015. [Google Scholar]
  10. Yang, D.C.H.; Blanche, J.G. Design and application guidelines for cycloid drives with machining tolerances. Mech. Mach. Theory 1990, 25, 487–501. [Google Scholar] [CrossRef]
Figure 1. Interpretation of profile failure and the planetary gear profile.
Figure 1. Interpretation of profile failure and the planetary gear profile.
Lubricants 09 00071 g001
Figure 2. Interpretation of the initial state and on the right schematic concept of the cycloidal drive.
Figure 2. Interpretation of the initial state and on the right schematic concept of the cycloidal drive.
Lubricants 09 00071 g002
Figure 3. The rotation of the input shaft required for the collision ( α ) .
Figure 3. The rotation of the input shaft required for the collision ( α ) .
Lubricants 09 00071 g003
Figure 4. Determination of the gap.
Figure 4. Determination of the gap.
Lubricants 09 00071 g004
Figure 5. Definition of the rotation angle of the planetary gear.
Figure 5. Definition of the rotation angle of the planetary gear.
Lubricants 09 00071 g005
Figure 6. Determination of the gear ratio.
Figure 6. Determination of the gear ratio.
Lubricants 09 00071 g006
Figure 7. Flow chart of the calculation process.
Figure 7. Flow chart of the calculation process.
Lubricants 09 00071 g007
Figure 8. Gaps in the case of certain cog-pairs.
Figure 8. Gaps in the case of certain cog-pairs.
Lubricants 09 00071 g008
Figure 9. Rotation angle of the planetary gear ( φ 2 ) depending on the rotation angle of the ( 3 ) eccentric cam ( φ 1 ) .
Figure 9. Rotation angle of the planetary gear ( φ 2 ) depending on the rotation angle of the ( 3 ) eccentric cam ( φ 1 ) .
Lubricants 09 00071 g009
Figure 10. Transmission change and fluctuation.
Figure 10. Transmission change and fluctuation.
Lubricants 09 00071 g010
Figure 11. Percentile deviation of the transmission.
Figure 11. Percentile deviation of the transmission.
Lubricants 09 00071 g011
Figure 12. Impact of the profile shift coefficient ( x ) on transmission fluctuation ( m = 1 ; z 1 = 9 ; d = 0.1 ) .
Figure 12. Impact of the profile shift coefficient ( x ) on transmission fluctuation ( m = 1 ; z 1 = 9 ; d = 0.1 ) .
Lubricants 09 00071 g012
Figure 13. Impact of profile failures ( d ) on the transmission fluctuation ( m = 1 ; z 1 = 9 ; x = 0.3 ) .
Figure 13. Impact of profile failures ( d ) on the transmission fluctuation ( m = 1 ; z 1 = 9 ; x = 0.3 ) .
Lubricants 09 00071 g013
Figure 14. The impact of different gear teeth with a constant profile error on the transmission fluctuation ( d = 0.01 ; x = 0.3 ) .
Figure 14. The impact of different gear teeth with a constant profile error on the transmission fluctuation ( d = 0.01 ; x = 0.3 ) .
Lubricants 09 00071 g014
Figure 15. Maximum transmission fluctuation in the case of m = 1 [ m m ] .
Figure 15. Maximum transmission fluctuation in the case of m = 1 [ m m ] .
Lubricants 09 00071 g015
Figure 16. Maximum transmission fluctuation in the case of m = 2 [ m m ] .
Figure 16. Maximum transmission fluctuation in the case of m = 2 [ m m ] .
Lubricants 09 00071 g016
Figure 17. Maximum transmission fluctuation in the case of m = 3 [ m m ] .
Figure 17. Maximum transmission fluctuation in the case of m = 3 [ m m ] .
Lubricants 09 00071 g017
Figure 18. Maximum transmission fluctuation in the case of m = 5 [ m m ] .
Figure 18. Maximum transmission fluctuation in the case of m = 5 [ m m ] .
Lubricants 09 00071 g018
Table 1. Profile failures in the case of m = 3 [ m m ] .
Table 1. Profile failures in the case of m = 3 [ m m ] .
z 1 9254353117
5 0.0045 0.0065 0.009 0.009 0.0125
6 0.0065 0.0095 0.0125 0.0125 0.018
I T 7 0.0105 0.015 0.02 0.02 0.0285
8 0.0165 0.023 0.0315 0.0315 0.0445
9 0.026 0.037 0.05 0.05 0.07
Table 2. Gear parameters.
Table 2. Gear parameters.
Module ( m ) : 1 [ m m ]
Number of teeth planet-gear ( z 1 ) : 9 [ ]
Number of teeth ring-gear ( z 2 = z 1 + 1 ) : 10 [ ]
Addendum modification ( x ) : 0.3 [ ]
Generating circle radius factor ( r c * ) : 1 [ ]
Profile failure ( d ) : 0.1 [ m m ]
Table 3. Gear parameters.
Table 3. Gear parameters.
Module range ( m ) : 1 , 2 , 3 , 5 [ m m ]
Pinion teeth range ( z 1 ) : 9 , 25 , 43 , 53 , 117 [ ]
Addendum modification ( x ) : 0.3 [ ]
Generating circle radius factor ( r c * ) : 1 [ ]

Short Biography of Author

Lubricants 09 00071 i001Attila Csobán Assistant professor at Budapest University of Technology and Economics. Member of the Association of Hungarian Inventors since 2000. Member of the Entrepreneurship Council of the Hungarian Research Student Association since 2006. Member of public body of Hungarian Academy of Sciences (MTA) since 2012. Gold level member of the European Who is Who Association since 2013. Research field: gear drives, gearboxes, planetary gear drives, cycloidal drives.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Csobán, A. Impacts of a Profile Failure of the Cycloidal Drive of a Planetary Gear on Transmission Gear. Lubricants 2021, 9, 71. https://0-doi-org.brum.beds.ac.uk/10.3390/lubricants9070071

AMA Style

Csobán A. Impacts of a Profile Failure of the Cycloidal Drive of a Planetary Gear on Transmission Gear. Lubricants. 2021; 9(7):71. https://0-doi-org.brum.beds.ac.uk/10.3390/lubricants9070071

Chicago/Turabian Style

Csobán, Attila. 2021. "Impacts of a Profile Failure of the Cycloidal Drive of a Planetary Gear on Transmission Gear" Lubricants 9, no. 7: 71. https://0-doi-org.brum.beds.ac.uk/10.3390/lubricants9070071

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop