Cycle Time Analysis of Open Pit Mining Dump Trucks

Abstract

This study demonstrates a practical cycle time analysis of dump truck haulage system of “Ukhaa Khudag” open-pit coal mine located in Umnugobi Province, Mongolia. It examines the possibility of minimizing the cycle time of the haulage system as well as factors impacting the speed of the dump truck. The current study divides the open pit mine road for the dump trucks into five sections which are bench road, ramp, surface road, dump road uphill, and dump road. Meanwhile, it investigates the influence of the length, the grade, and the rolling resistance of the road section on the cycle time. The data is analyzed using mathematical regression methods via Microsoft Excel program. For each of the five road sections, we compare the statistical calculations of three regression models: linear, quadratic and exponential; thus, a total of thirty regression models are obtained in this research. Accordingly, the cycle time for each road section is predicted by the most accountable model. The loaded and empty direction of the movement is measured and calculated for each road section, and it appears that the difference between the calculated mean value and the actual cycle time of the models is 0.82 seconds with a relative error of 2.51 percent.

Share and Cite:

Enkhchuluun, B. , Batgerel, B. and Ping, C. (2023) Cycle Time Analysis of Open Pit Mining Dump Trucks. International Journal of Geosciences, 14, 689-709. doi: 10.4236/ijg.2023.148037.

1. Introduction

Mining is one of the costly and complicated industries. Widespread studies have been conducted in relevant sections such as geology, canalization, and excavation planning and operation process [1] . Particularly, transportation costs amount to 50 to 60 percent of the total investment costs and 70 percent of the operation costs [2] [3] [4] . In open pit mining, the movement of raw materials is considered one of the most challenging tasks with the truck transportation representing the most influential factor of mining costs [5] . The deeper the open pit, the longer the open pit road is needed; thus leading to higher transportation cost. Accordingly, in order to reduce the transportation cost, it is necessary to investigate the potential strategies to improve the open-pit mine cycle time of the dump trucks.

The current study analyzes the dump truck haulage system of “Ukhaa Khudag” Open Pit Mine in Umnugobi Province, Mongolia which started its operations in 2009 and produced ten million tons of coals in 2018. Due to its potential growth and expected long life cycle, we aim to conduct a detailed survey on cycle time of the dump truck. The main contribution is to identify the potential opportunities to reduce the costs of dump truck systems by evaluating the mining road conditions in Mongolia which had not been rigorously studied. Specifically, this study surveys 13 trucks of CAT 785 model. Two main channels to improve the efficiency of the dump trucks are to increase truck capacity and to reduce cycle time [6] . As the movement time is the most important factor of overall cycle time, it is crucial to understand which factors significantly affect the speed of the movement so that we can optimize the speed and increase the efficiency of the dump truck [6] . To improve the efficiency of the dump truck, an optimal use regime for open-pit mining trucks should be established [7] [8] . The determinant of the optimal use regime includes the movement regime which is directly related to the speed of the dump truck.

The type of the road pavement and the grade of the road have a great influence on the speed of the dump truck [9] [10] . If the speed of the dump truck is accurately determined, the production amount as well as the number of trucks required will be calculated correctly [6] [2] . The speed of the dump truck can be controlled by optimizing the condition of each road section. Examining the effects of the open pit mine dump truck in the aspect of the movement regime can help us determine which parts bear the most cost and which strategies should be implemented to reduce the costs [7] .

The road pavement expresses the direct rolling resistance of the road. Rolling resistance is a measure of the extra resistance to the motion that a haul truck experiences, and it is influenced by tyre flexibility, internal friction and most importantly, wheel load and road conditions [9] .

Roger J. Thompson (2013) has thoroughly researched various mining road conditions and reported that the optimal overall road resistance of the dump truck is at 8 to 11 percent.

The current study employs regression technique. Our goal is that the equations resulting from our calculations are only applicable to this open pit mine.

2. Methods

2.1. Cycle Time of Dump Truck

Cycle time is defined as the time required for any equipment to complete one cycle of operation. For a truck, cycle time comprises the time to spot and load, travel to the dump site, maneuver, spot and dump, and drive back to the loading point, also inclusive of predictable delays and unpredictable waiting times (Lineberry, 1985) [11] [12] .

Principally, cycle time is a direct measure of process and equipment performance in waste material transport for open-pit mining and other processes [13] [14] .

Cycle time analysis was conducted starting with identification of the activities to be analyzed, for which the cycle times of seven activities were determined and recorded. In accordance with Samwel Victotmanyele (2017) [13] , the total cycle time for the truck was appointed as per Equation

T c t = t q + t s p + t l d + t f h + t d p + t h g + t e h (1)

where t q : queuing time at the excavator, t s p : spotting time, t l d : loading time, t f h : full haul time, t d p : dumping time, t h g : hanging time and t e h : empty haul time.

2.2. Determine the Road Parameters and the Total Resistance

Road conditions affect the technical and economic performance of open pit mine dump truck. The road of the open-pit mine is classified by its character. Relevant characters include:

Structure: paved and unpaved;

Movement direction: one lanes, two or more lanes in traffic;

Time using: permanent and temporary;

Location: bench road, ramp, surface, dump road.

Rolling resistance is defined as the force required to maintain a vehicle at a steady speed on the ground level, and is a function of not only the gross vehicle mass and driving characteristics, but also the type and conditions of the tires and the road surface on which the vehicle is operated [9] . The characteristics of the road pavement are expressed as a coefficient of the rolling resistance.

Empirical estimations of the rolling resistance are based on tire penetration, and it turns out that 0.6 percent increase in rolling resistance per centimeter of tire penetration into the road typically result in 1.5 to 2 percent minimum resistance [9] .

Figure 1 shows that rolling resistance for mine haul road depends on road conditions. Calculations of the rolling resistance are categorized as shown in Table 1.

According to B. Purevtogtokh (2018), the average value of the road rolling resistance in the direction of the load section j is defined as follow.

w l j = l 1 w 1 + l 2 w 2 + + l n w n l 1 + l 2 + + l n = i = 1 n l i w i i = 1 n l i , N/kN (2)

The average value of the road rolling resistance in the direction of the empty section j is set as:

Figure 1. Rolling resistance.

Table 1. Calculation of rolling resistance for mine haul road.

w e j = w l j 1.25 , N/kN (3)

At the same time, the average slope of the road in the load and empty section j is specified as:

i j = l 1 i 1 + l 2 i 2 + + l n i n l 1 + l 2 + + l n = i = 1 n l i i i i = 1 n l i , ‰ (4)

where

w i : the rolling resistance of the road in the load of direction;

l i : the length of the road section i for traffic j, km;

i i : the slope of the road section i for traffic j, ‰.

Meanwhile, the total resistance of the movement is described as:

W j = P l 1 ( w l 1 + i l 1 ) l l 1 + P e 1 ( w e 1 + i e 1 ) l e 1 L 1 , N (5)

where: P l 1 , P e 1 : the weight of the dump truck in the load and empty direction, tons;

w l 1 , w e 1 : the rolling resistance of the road in the load and empty direction, N/kN;

l l 1 , l e 1 : the length of the road in the load and empty direction, km;

i l 1 , i e 1 : the slope of the road in the load and empty direction, ‰;

L 1 : the total distance of road.

It is essential to develop a methodology to optimize the utilization regime in connection with the increase in the volume of transport. Consequently, a scientific assessment of the organizational level, operating conditions and quality characteristics should be conducted; otherwise it is impossible to maximize the efficiency of the dump truck.

2.3. Multiple Linear Regression

Multiple linear regression analysis is a statistical method used for determining a formula that explains the prediction of a dependent (unknown) variable by its relationship to a set of independent (known) variables. A weighting is assigned to each independent variable to reflect the portion of its impacts on the value of the dependent variable; thus, the weighting is referred to as the coefficient of the variable in the equations [15] [16] .

Multiple R. Multiple R is the Correlation Coefficient, illustrating the strength of the linear relationship. For instance, a value of 1 signifies a perfectly positive relationship while a value of 0 expresses no relationship between a certain independent variable and the dependent variable. Multiple R is the square root of R squared.

R squared. R squared (r2) is the Coefficient of Determination, representing the proportion of the variance for a dependent variable that is explained by the independent variables included in a regression model. For instance, R squared value of 80 %describes that 80 percent of the variations of the value of y fall on or closed to the regression line (the mean) are explained by the value of x; thus, implying that 80 percent of the values calculated are fit to the model.

Adjusted R square. Adjusted R-square (adjusted R2) shows how well the data points fit the regression line, but adjusts for the number of terms in a regression model.

Standard Error of the Regression. Standard Error of the Regression is an estimate of the standard deviation of the error μ, which is not similar to the standard error in descriptive statistics. The standard error of the regression is the precision that the regression coefficient is measured; if the coefficient is large compared to the standard error, then the coefficient is probably different from 0.

Observations. Observations report the number of observations in the sample.

2.4. Simplification of the Road Scheme

In order to minimize the amount of road sections, a summary of the road parameter must be checked to ensure the road conditions for possible tractive forces, define the speed of movement and the time as well as the efficient usage frequency. Afterwards, the principle road scheme is simplified in accordance with associated parameters.

Before the road is simplified, following road conditions must be taken into account.

· Summarized roads must be a same type of road section (bench road, ramp, surface road, dump road).

· Pavement and structure of road must be the same.

· Road gradient and rolling resistance are nearby for dump trucks operating in the same regime.

· The speed of the dump truck is nearby.

· Safety conditions must be the same.

The R radius curved sections is converted to the slope of straight section. The curved sections of radius R are defined as the curve resistance according to the equation:

Δ i R = 30 200 R 200 , ‰ (per mille) (6)

where: R: turn radius, m.

The radius R curve in section j is defined according to the equations:

i j = i j + Δ i R , ‰ (7)

where: i j : slope of the j section curve converted into slope, ‰;

i j : Slope of the j section curve, ‰.

The road has to be subdivided into sections of the road, namely bench road (BR), ramp (R), surface road (SR), dump road uphill (DRU), and dump road (DR).

For each section of road, the average length, the average slope and the average rolling resistance are defined. In addition, the amount of charge is considered.

For example; weighted average length of the bench road:

L B R = L B R j Q j Q j , km (8)

Weighted average slope of the bench road:

i B R = L B R j i B R j L B R j , km (9)

Weighted average rolling resistance of the bench road:

w B R = L B R j w B R j L B R j , N/kN (10)

where; LBRj: length of bench road in j section, km;

IBRj: slope of bench road in j section, ‰;

Qj: quantity of transport load by section j, tn/m3;

WBR: rolling resistance bench road in j section, N/kN (‰).

3. Practical Experiments of Transport Process

The test of dump truck regime determines the dimensions of the road and the speed (time) of each road section. Testing process consists of measuring chips in each dump truck and recording the cycle time. The determination of rolling resistance is based on the knowledge on the road surface and the current road defect. Meanwhile, road gradient and road length are determined by the measurement of the open-pit mine surface.

The road description from the excavation point to the discharge point is shown in Table 2.

Figure 2 exhibits the scheme of an open-pit mine road for a dump truck. Points 1, 2, 3, 4, 5 are loading positions and Point 12, 15, 18 are discharge positions.

Simplifying the Calculation of the Road Section

For the sake of simplicity, the same type of road is considered for this open pit mine. The road with five bench roads, three ramps, three surface roads, three dump uphill roads and three dump roads are combined and calculated for this open pit. The average length, average slope, and average rolling resistance of each road are calculated and the following equations are derived. The results of the simplified transport system calculations are shown in Table 3.

Weighted average length of the bench road:

L B R = Q 1 L 1 6 + Q 2 L 2 7 + Q 3 L 3 8 + Q 4 L 5 8 + Q 5 L 4 7 Q = 0.26 km (11)

Weighted average length of the ramp:

L R = Q 1 L 6 7 + Q 2 L 7 8 + Q 4 L 8 9 Q = 0.358 km (12)

Weighted average length of the surface road:

Table 2. Haulage systems description.

Figure 2. Scheme of the road section.

L S R = Q ( L 9 13 + L 9 10 + L 9 16 ) Q = 0.600 km (13)

Weighted average length of the dump road uphill:

L D R U = Q ( L 10 11 + L 13 14 + L 16 17 ) Q = 0.334 km (14)

Weighted average length of the dump road:

L D R = Q ( L 11 12 + L 14 15 + L 17 18 ) Q = 0.387 km (15)

Weighted average slope of the bench road:

i B R = i 1 6 L 1 6 + i 2 7 L 2 7 + i 3 8 L 3 8 + i 5 8 L 5 8 + i 4 7 L 4 7 L = 8.39 (16)

Weighted average slope of the ramp:

i R = i 6 7 L 6 7 + i 7 8 L 7 8 + i 8 9 L 8 9 L = 78.43 (17)

Weighted average slope of the surface road:

i S R = i 9 13 L 9 13 + i 9 10 L 9 10 + i 9 16 L 9 16 L = 1.33 (18)

Weighted average slope of the dump road uphill:

i D R U = i 10 11 L 10 11 + i 13 14 L 13 14 + i 16 17 L 16 17 L = 59.62 (19)

Weighted average rolling resistance of the bench road:

w B R = L B R j w B R j L B R j = 75.46 , N/kN (20)

Weighted average rolling resistance of the ramp:

w R = L R j w R j L R j = 40.99 , N/kN (21)

Weighted average rolling resistance of the surface road:

w S R = L S R j w S R j L S R j = 36.41 , N/kN (22)

Weighted average rolling resistance of the dump road uphill:

w D R U = L D R U j w D R U j L D R U j = 152 , N/kN (23)

Weighted average rolling resistance of the dump road:

w D R U = L D R j w D R j L D R j = 152 , N/kN (24)

In this case, three roads have a curve radius and have been converted to a straight road gradient as illustrated in Table 4.

Figure 3 shows a simplified representation of road sections. Measured values of road are grouped into bench road, ramp, surface road and dump road.

4. Result and Discussion

4.1. Factors That Influence the Movement Process of Dump Trucks and the Mathematical Model

The current study observes every road section and dump truck movement. Subsequently, the data is analyzed via multiple mathematical regression models using mathematical statistical methods. The data regarding five bench roads, three ramps, three surface roads, three dump uphill roads and three dump roads are combined and calculated and accordingly the measurement of “UkhaaKhudag” open-pit mine has been summarized. The results show the estimated ranges obtained via Microsoft Excel application as outlined in Table 5.

The Microsoft Excel program compares statistical calculations with of regression models: linear, quadratic and exponential. Factors are validated by the F-test. Table 6 shows comparisons of road moving time regression models in each section in the loaded direction.

Table 3. Simplified Haulage systems description.

Table 4. Converted slope for the curved sections of radius R.

Figure 3. Simplified scheme for road sections.

Table 5. Descriptive statistics of road sections with duration of movement.

Table 6. Comparisons of the moving times regression model of the road in each section (loaded).

Table 7. Comparisons of the moving times regression model of the road in each section (empty).

Table 7 shows comparisons of road moving time regression models in each section in the empty direction.

By comparing the regression models described in each road section, the models with greater correlation coefficient and smaller standard errors are selected; therefore the prediction models are acquired.

For the modeling variables: L: the length of the section on the road, km;

i: the grade of the section on the road (‰);

w: rolling resistance of the section on the road, N/mN (m-mega).

4.2. Demonstrating the Model of Each Road Section

Bench road in loaded direction

T = 21.12 + 140.50 L + 0.043 i 279.68 w (25)

Road length and grade have a direct effect on moving time of cycle whereas rolling resistance has an indirect effect.

Figure 4 and Figure 5 prove that the regression model described in the direction of the bench road is similar to the value of the actual measurement.

Bench road in empty direction

T = 101.37 + 66.31 e L + 0.0004 e i 147.25 e w (26)

The model indicates that an increase of 179.7 seconds movement time results

Figure 4. Comparison of measured time and calculated time of bench road loaded direction.

Figure 5. Comparison of measured time and calculated time of bench road empty direction.

in one unit increase of the road length. Furthermore, one unit increase of pavement road leads to a decrease of 399 seconds movement time.

Ramp in loaded direction

T = 36.99 + 134.67 L + 0.79 i + 121.69 w (27)

The regression model reveals that one unit increase of the rolling resistance causes an increase of 121.69 seconds movement time. One kilometer increase of the road length contributes to an increase of 134.67 seconds movement time. In the meantime, one unit increase of the road slope per mille conduces to an increase of 0.79 seconds movement time.

In summary, three factors included in this model significantly impact the cycle time.

Figure 6 reflects the similarities between the regression model described in the direction of ramp and the value of the actual measurement.

Ramp in empty direction

T = 5 + 64.74 L 0.29 i + 52.44 w (28)

The regression model demonstrates that one unit increase of the rolling resistance and the length of the road lead to the increase of the movement time for 52.44 and 64.74 seconds respectively.

As evidenced in Figure 7, the differences in the values of the errors have a higher value compared with the other figures.

Surface road in loaded direction

T = 217.12 + 23.23 e L 6.67 e i + 239.17 e w (29)

Surface road in empty direction

T = 34.25 + 17.25 e L 6.12 e i + 62.36 e w (30)

The regression models shown in Figure 8 and Figure 9 are drawn from

Figure 6. Comparison of measured time and calculated time of ramp loaded direction.

Figure 7. Comparison of measured time and calculated time of ramp empty direction.

Figure 8. Comparison of measured time and calculated time of surface road in loaded direction.

Figure 9. Comparison of measured time and calculated time of surface road in empty direction.

related exponential functions. It is obvious that the road pavement has larger influence on the movement time. Hence, the pavement should be highly considered as the surface road length increases.

Dump road uphill in loaded direction

T = 25287.04 256924.76 L 2 1.03 D 2 + 570.85 L D + 160942.16 L 156.76 D (31)

Dump road uphill in empty direction

T = + 578.552 13619.027 L 2 0.1684488 D 2 288.4464 L D + 2913.8509 L 95.6335 D (32)

The regression models portray the values of the slope while main resistance factors on the road are expressed in the term of D;

D = ( i + w ) / 10 (33)

where: i: the grade of the section on the road (‰);

w: the rolling resistance of the section on the road, N/kN.

Figure 10 and Figure 11 manifest that the difference between the calculations and actual measured values varies from 0.32 to 0.55 in the loaded direction and 0.34 to 0.52 in the empty direction; however, the proportion is comparatively small.

Dump road in loaded direction

T = 254.87 + 72.88 L 1457.53 w (34)

Dump road in empty direction

T = 1040.43 + 448.69 L 2 16221.3 w 2 + 1304.72 L 8663.95 L w + 9209.55 w (35)

The turning radius and the slope of the dump road are relatively small, for this

Figure 10. Comparison of measured time and calculated time of dump uphill road in loaded direction.

Figure 11. Comparison of measured time and calculated time of dump uphill road in empty direction.

Figure 12. Comparison of measured time and calculated time of dump road in loaded direction.

reason, the regression model should be determined by including only two variables: rolling resistance and the road length.

The regression models emphasize that one unit increase of rolling resistance is attributed to a reduction of 1457.53 seconds movement time in loaded direction. Simultaneously, one kilometer increase of the road length is the cause of an increase in the movement time for 72.88 seconds and 36.25 seconds in the loaded direction and the empty direction respectively. As a result, it is necessary to define the road exactly as the length increases. Figure 12 and Figure 13 show the comparison of the calculated and measured values of the defined models in the loaded and empty direction.

4.3. Research of the Dump Truck Stop Times

This section measures the waiting time of the excavator, spotting time, loading

Figure 13. Comparison of measured time and calculated time of dump road in empty direction.

Table 8. Load and dump time of dump trucks (seconds).

time, dumping time, and dumping spot time. The experimental measurement exemplifies the combination of Excavator R9250 and Truck CAT785, and the results are displayed in Table 8.

It is affirmed that waiting usually takes longer time than most of the other processes. Moreover, the difference between the minimum and the maximum records of the waiting time appears extremely large. Thereby, it is necessary to optimize the combination of trucks and shovels.

Corresponding to Table 9, the difference between the actual and mathematical model of the mean time is 0.82 seconds. Regarding the errors, the absolute error of 0.82 seconds with the relative error of 2.51 percent confirms the probability value of the mathematical models. Table 10 shows the total cycle times of the dump trucks.

5. Conclusion

· The results of the study on movement duration reveal that it is possible to increase the productivity by 34.19 percent according to the comparison between the maximum duration and average duration. Regarding the practical experiments, the regression equations are derived and can be used for defining the movement time of both loaded and empty direction for each road section.

Table 9. Comparison of actual and calculated time for road section (seconds).

Table 10. Cycle time of dump trucks (seconds).

· Experimental results disclose that it is feasible to reduce the waiting time while emphasizing the advancement of spotting method and the combination of trucks and shovels.

· The mathematical models developed during the current study can be utilized for future open pit mine planning of which the key factors generally include the type of pavement, the slope of the road, and the cycle time.

· In general, improvement of the road increases the speed; however, for mining, the roads are usually for temporary use so that not much investment is dedicated to road construction. As reported by this study, it is crucial to pay attention the road built by focusing on its expected life cycle.

· Specifically, the equations obtained via the regression analysis clearly identify which road sections significantly influence the speed of the trucks.

Acknowledgements

My thanks go to Dr. Boryo Orkhontuul from the Mining Department of the Mongolian University of Science and Technology, who gave me great help in writing this project. My deepest thanks go to the mining engineers who worked at the Ukhaa Khudag (UHG) open pit mine.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Eskandari, H., Darabi, H. and Hosseinzadeh, S.A. (2013) Simulation and Optimization of Haulage System of an Open-Pit Min. Summer Computer Simulation Conference, Toronto, July 2013, Article No. 13.
[2] Purevtogtokh, B. and Orkhontuul, B. (2015) Determination of Operating Regimes of Dump Truck on “Ukhaakhudag” Open Pit Coal Mine. Research Report.
[3] Thompson, R.J. and Visser, A.T. (1999) Designing and Management Unpaved Opencast Mine Haul Roads for Optimum Performance. SME Annual Meeting, Denver, 1-3 March 2013, 90-99.
[4] Liu, G.W. and Chai, S.L. (2019) Optimizing Open-Pit Truck Route Based on Minimization of Time-Varying Transport Energy Consumption. Mathematical Problems in Engineering, 2019, Article ID: 6987108.
https://doi.org/10.1155/2019/6987108
[5] Alegre, D.A.G., Peroni, R.L. and Aquino, E.R. (2019) The Impact of Haulroad Geometric Parameters on Open Pit Mine Strip Ratio. REM, International Engineering Journal, 72, 25-31.
https://doi.org/10.1590/0370-44672018720136
[6] Purevtogtokh, B. (2018) Transport of Open Pit Mining Trucks. Study Book.
[7] Orkhontuul, B. (2006) Definition Methods of Parameters of Truck Exploitation in Open Pit Mining. Ph.D. Thesis, Mongolian University of Science and Technology, Ulaanbaatar.
[8] Orkhontuul, B. (2007) Analysis of Truck Dismounting Process. IEEE International Forum on Strategic Technology, Ulaanbaatar, 3-6 October 2007, 457-458.
https://doi.org/10.1109/IFOST.2007.4798630
[9] Thompson, R.J. (2015) Principles of Mine Haul Road Design and Construction. Study Book.
[10] (2018) Caterpillar Performance Handbook 48. Peoria, Illinois, U.S.A. Mining & Off-Highway Trucks 10-21, Mining and Earthmoving-28-6.
https://wheelercat.com/wp-content/uploads/2018/07/SEBD0351_ED48.pdf
[11] Coronado, V. and Pablo, P. (2014) Optimization of the Haulage Cycle Model for Open Pit Mining Using a Discrete-Event Simulator and Context-Based Alert System. Master Thesis, The University of Arizona, Tucson.
[12] Rahimel, M.J., Mirzaei, M., Sattervand, J. and Hoseinie, S.H. (2017) Health Risk of Whole Body Vibration in Mining Trucks during Various Operational Conditions. Journal of Central South University, 24, 1808-1816.
https://doi.org/10.1007/s11771-017-3589-3
[13] Manyele, S.V. (2017) Analysis of Waste-Rock Transportation Process Performance in an Open-Pit Mine Based on Statistical Analysis of Cycle Times Data. Engineering, 9, 649-679.
https://doi.org/10.4236/eng.2017.97040
[14] Krause, A.J. (2006) Shovel-Truck Cycle Simulation Methods in Surface Mining. Master Thesis, University of the Witwatersrand, Johannesburg.
[15] Chatterjee, S. and Hadi, A.S. (2006) Regression Analysis by Example 4th Edition. Study Book.
[16] Chanda, E.K. and Gardiner, S. (2010) A Comparative Study of Truck Cycle Time Prediction Methods in Open-Pit Mining. Engineering, Construction and Architectural Management, 17, 446-460.
https://doi.org/10.1108/09699981011074556

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.