Design, Fabrication and Performance Evaluation of a Multipurpose Cassava Grater

Abstract

Food security has been an issue of global concern and this has attracted a lot of research interest. Cassava is an extremely popular crop and is becoming the cornerstone for addressing food security in many parts of the world. The competing needs for cassava cuts across both human and animal consumption. It serves as a raw material in textile industry and is now one of the preferred materials for making biofuels. As the world’s population continues to grow, the demand for drought resistant crops such as cassava is increasing. The high demand for various forms of processed cassava will continue to increase cassava prices making it an attractive business venture. Several small-scale cassava farmers are making a fortune and additional income through this business. Preliminary investigations show that the profitable way for a farmer to market his cassava is to add value to it. Unfortunately, cassava undergoes post-harvest physiological deterioration (PPD) after three days of harvest. In order to make cassava farming even more profitable, there is a need to process it within the shelf life of 2 - 3 days after harvesting. One way to preserve cassava is by grating it into pulp and drying it into pellets or chips. In this study, an electrically powered multi-purpose cassava grating machine with grater blades inclined at two different tooth angles, 25˚ and 30˚ was designed, fabricated and its performance characteristics investigated. The results showed that the plate with a tooth angle of 30˚ resulted in higher grating efficiency. This was attributed to better grip on the cassava when perforations on the plate inclined at 30˚. The grating capacity was also significantly improved as very small amounts of cassava slipped out un-grated.

Share and Cite:

Mose, B. , Odhiambo, V. , Maryvelma, N. and Alunda, B. (2025) Design, Fabrication and Performance Evaluation of a Multipurpose Cassava Grater. World Journal of Engineering and Technology, 13, 49-63. doi: 10.4236/wjet.2025.131004.

1. Introduction

Cassava was introduced to Africa by the Portuguese more than 300 years ago and is a major source of carbohydrate in sub-Saharan Africa. Compared with other carbohydrate sources such as rice and maize, cassava has the advantage of high productivity under marginal climatic conditions. Cassava has the capacity to yield under marginal soil conditions and is tolerant to drought [1]. The crop, therefore, plays a vital role in the food security of the world [2]-[11].

Over the past 25 years, significant market opportunities for cassava have opened up in the animal feeding industry. Moreover, cassava starch is an important ingredient in the manufacture of dyes in textile industry, drugs, chemicals, carpets and in coagulation of rubber latex. This is in addition to the capacity of cassava to feed the rapidly increasing population.

Despite the numerous advantages as a source of food for humans, livestock feed and for industrial applications, cassava undergoes post-harvest physiological deterioration (PPD) after three days of harvest due to its high moisture content of about 70%. Therefore, there is a need to process cassava within three days after harvesting in order to obtain the multiple economic benefits from the crop. The best form of preservation and reduction of post-harvest losses is the immediate processing of cassava into a stable product such as pulp which can be easily and properly dried.

Processing cassava by traditional methods suffers low productivities and a susceptibility to food contamination. Various attempts to design and fabricate cassava-grating machines have been made with a view to overcome the negative attributes of the traditional processing techniques and promote timely, large-scale processing of the tubers in a hygienic environment. These machines, however, are usually found to suffer significant corrosion due to the acidic nature of the cassava fluid [12]. Non-food grade materials used for the fabrication of cassava processing machines have been found to react with the fluid to cause contamination.

It is, therefore, essential to come up with alternative machine designs with materials that do not easily corrode. This is because cassava is becoming a vital means of livelihood to feed the global population with multiple other economic benefits, thus, having the ability to generate revenue and income for developing countries.

In this study, an electrically powered multi-purpose cassava grating machine with grater blades inclined at two angles, 25˚ and 30˚ was designed, fabricated and its performance characteristics investigated. The developed multipurpose grating machine can be used for grating not only cassava tubers but also sweet potatoes, irish potatoes, yams, arrow roots, carrots, beetroots, bananas as well as pumpkins. Moreover, the grater can be used to process medicinal herbs such as aloe Vera, and other fruits such as avocados, raw mangoes. As a result, the proposed machine can trigger and promote processing of cassava tubers as well as other farm products into commercially acceptable products including medicinal herbs, thereby increasing income obtainable by local farmers.

2. Materials and Method

The schematic of the main components of the cassava grating machine is shown in Figure 1.

Figure 1. Main components of the grating machine.

2.1. Design Considerations

A number of considerations were put in place while designing the various components of the grater. These considerations include the following;

1) Attain high throughput of cassava grating by reducing grating time.

2) Reduce high labor requirements of traditional manually operated graters.

3) Reduce corrosion and enhance hygienic conditions by use of food grade material.

4) Enable easy disassembly, cleaning and assembly of the machine.

5) Achieve economic viability and reduce cost.

2.2. Design of Various Components

2.2.1. Hopper Design

The hopper was designed based on the quantity of cassava to be processed to meet the farmer/market demands. It was found that an average cassava farmer with acreage of 1 - 2 acres [13] required 336 kg of cassava to be grated per hour, which translates to 5.6 kg of cassava per minute. Using Equation (1), for a mass of 5.6 kg/min and an average cassava density of 670 kg/m3, the volume of the hopper was estimated as follows;

Volume= Mass Density (1)

V= 5.6 kg/ min 670 kg/ m 3 =8.3582× 10 3 m 3

The volume of hopper to grate 5.6 kg of cassava per minute was found to be at least 8,358,000 mm3. The next critical item concerning the hopper was the shape. The hopper was designed to have a pyramidal shape with an angular inclination of 50 degrees to allow cassava to be fed into the grating chamber by the force of gravity so as to reduce human effort [14].

2.2.2. Belt Design

Determination of Center Distance (C)

The Rubber Manufacturing Association recommends that the distance between the center of the driving pulley and that of the driven pulley be obtained from consideration of the speed ratio,

Speed Ratio= D d = 0.2 0.1 =2

For speed ratios less than three:

c= D+d 2 +d (2)

where:

D = Diameter of the driving pulley= 0.27m;

d = Diameter of the driven pulley = 0.135m;

Center distance was computed and found as;

C=0.34m

Length of Belt

A V-belt was chosen to transmit power generated by the motor to the shaft since it requires less installation work and time, absorbs shock and operates at low bearing pressure. The formula in Equation (3) was used to compute the length of the belt;

L= π( D+d ) 2 +2C+ ( Dd ) 2 4C (3)

where:

C = Distance between the center of the two pulleys = 0.34;

L = Length of Belt Required was upon substitution found to be 1.3177 m.

A50 V belt of length 1.3 m was selected due to its low vibration, less noise, high efficiency which facilitates longer motor and V-belt life.

2.2.3. Motor Design

Capacity of the electric motor

Power transmitted by motor; P is given by Equation (4);

P= πF N 1 D 1 60 (4)

where:

F is the force to crush fresh cassava 105N;

N1 is the rotational speed of motor 1140 rpm;

D1 is the diameter of motor pulley 0.135 m.

Applying Equation (4), the power rating is obtained as 1068 Watt. A motor with a power rating of 1.5 HP = 1119-Watt, 13 Ampere, 240 volts, 1440 rpm, 50 Hz was thus selected.

2.2.4. Shaft Design

To transmit power from the shaft, various members/elements such as the pulleys, bearings and the grating unit were mounted on it. These members along with the force exerted upon them cause the shaft to bend. Therefore, the shaft is subjected to bending and torsional forces since it is utilized for transmission of torque and bending moment.

Belt Tension in Open Belt Drive

Belt tension exists when power is transmitted from one shaft to the other as indicated in Figure 2.

Figure 2. Open Belt Drive for power transmission.

The maximum tension on the belt was determined using Equation (5) [15]

T m = τ all ×cross sectional area of belt (5)

where; τ all is the allowable shear stress on the belt.

For belt materials, allowable shear stress is 2.8 MPa. For A50 V-belt, cross sectional area is 20 mm2. The maximum tension was found to be 56N.

Elements Mounted on Shaft

The elements mounted on the shaft along with the force exerted upon them cause bending and torsional forces.

Determination of the Bending Moment at each Point of Loading

In order to determine the bending moment at each point of loading of the shaft, it was necessary to prepare the bending moment diagram.

Figure 3. The main machine elements mounted on the shaft.

Force Exerted on Shaft (Vertical Force)

The main machine elements that exert forces on the shaft are the driven belt pulley and grating unit as shown in Figure 3.

The weights of pulley and grating unit were obtained by getting the product of material density and volume as follows.

Weight of the pulley in kg, W p = ρ p V p =2700×3.5× 10 4 =1.0kg

Weight of the grating unit in kg, W g = ρ g V g =2700×1.1× 10 3 =3kg

where ρ p and ρ g are densities of pulley and grater materials, V p and V g are pulley and grater volumes. Point loading, W d on shaft due to hopper and cassava was determined,

W d =g [ ( V c × ρ c )+( ρ H × V H ) ]=404.94N

W d =9.81[ ( 4.5× 10 3 )×675+7280×( 4.5× 10 3 ) ]=351.2N

where; ρ H , ρ c = Density of hopper material and density of cassava respectively.

For cassava, ρ c =675 kg/m3 while for hopper material, ρ H =7280 kg/m3.

Reactions at the Bearings due to Vertical Loading

In Figure 4, the expected free body diagram of the vertical forces acting on the shaft is shown. The forces acting on the shaft due to the pulley, Fp and due to the grater Fg were obtained by multiplying the weights, Wp and Wg by the acceleration due to gravity, 9.81 m/s2. Note that RBV and RDV denote the reactions in the vertical direction at point B and D respectively.

Figure 4. Free body diagram of the main forces and reactions acting on the shaft.

To determine the reactions at each bearing, moments about the two expected bearing points (B and D) were determined. Taking moments about point B;

RDV×0.3+9.81×0.05=351.2×0.15+29.43×0.35

RCV=208.3N

Taking Moments about D we get

RBV×0.3+29.43×0.05=351.2×0.15+9.81×0.35

RBV=182.14N

Shear Force Diagram

The shear force diagram was drawn as shown in Figure 5.

Figure 5. Shear force diagram.

Bending Moment Diagram

From the evaluation of the forces and determination of the bearing reactions, the maximum bending moment ( M max ) for the shaft was evaluated. The following equations were used for bending moment diagram.

Section A-B: M=9.81x ;

Section B-C: M=172.33x9.107 ;

Section C-D: M=178.87x+61.133 ;

Section D-E: M=29.43x11.772 ;

where x is the distance from point A in meters.

The bending moment diagram (BMD) was drawn as shown in Figure 6. As noted, the maximum bending moment was found to be 25.36 Nm.

Figure 6. Bending moment diagram.

For a solid shaft having little or no axial load, the diameter of the shaft was determined from Equation (6) [16].

D s 3 = 16 π S s [ ( K b M b ) 2 + ( K t M t ) 2 ] 1 2 (6)

where; Ds is diameter of the shaft, Mt is torsional moment, Mb is bending moment, Kb is the combined shock and fatigue applied to bending moment, Kt is the combined shock & fatigue applied to torsional moment and Ss = allowable stress.

For a shaft transmitting power (W) at a given rotational speed (rpm), the torque is given by Equation (7a).

T= P ω (7a)

where; ω is the angular velocity and is given by;

ω= 2πN 60 rad/s

From motor design, it was found that power rating of motor is 1.5 HP = 1119-Watt, 13 amperes, 240 volts, 1440 rpm, 50 Hz and with a speed reduction of 1:2 to give the driven speed of 720 rpm, torsional moment, Mt is given by Equation (7b);

M t = Power/ Speed = 1119×60 2π×720 =14.84Nm (7b)

Given that Ss = 20 × 106 N/m2, Kb = 1.5, Kt = 1, Mb = 25.36 Nm, Mt = 14.84 Nm and π = 3.142; and upon substituting these values in Equation (6) we get;

D s 3 =10388.77 mm 3 fromwhich D s =21.82mm

From table of standard shaft sizes, the 25 mm diameter was selected [15]. The shaft was joined to the driving motor by means of keys, and aligned using ball bearings. The selection of this diameter also takes into account of any frictional forces at the points where bearings are located.

2.2.5. Selection of Ball Bearings for Pillow Block

Deep Groove Ball bearings were used since they can withstand the thrust action in the machine and are easy to repair and assemble. To determine the equivalent load when only a radial load, R, is applied Equation (8) was used.

F eq =VR (8)

where V is rotation factor and V = 1.0 if the inner race of the bearing rotates, or V = 1.2 if the outer race rotates. Noting that this is a pure radial load, the inner race is to be pressed onto the shaft and rotate with it. Therefore, the factor V= 1.0 in Equation (8) was applied.

The deep-groove ball bearing is to carry 100 kg of pure radial load from a shaft that rotates at 720 rpm. The design load is obtained as 100 kg. The bearing is to be mounted on a shaft with a diameter of 25 mm. Based on these parameters, bearing number 6205 was selected. The summary of data for the selected bearing: bearing number: 6205, single-row, deep-groove ball bearing, bore: d = 25 mm, outside diameter: D = 52 mm, Width: B = 15 mm.

2.2.6. Perforated Plate

The perforated plate, which is a critical component in cassava grating machines was made from food grade stainless steel. The perforated plate was fixed on the grating unit which was in turn mounted on the shaft. The force acting on the plate is given by Equation (9).

F= V p ρg (9)

where g is acceleration due to gravity (9.81 m/s2), Vp is volume of the plate and ρ is the material density, F is the force required to crush cassava. Rearranging Equation (9), we can get the diameter of the perforated plate using Equation (10).

π D 2 4 = F ρgt (10)

Substituting the values, diameter of perforated plate was found to be 306 mm.

Tooth Angle of the Perforations

The angle of inclination of the perforations was calculated based on the radius of the perforations and the height of the inclination to be made. For this case a radius

of 4 mm and a height of 8.5 mm were used. Applying the formula tanθ= r h an

angle of 25˚ was obtained. 60˚ staggered round holes selected for the perforation pattern because they are the most economical, strongest and most versatile as compared to rectangular, square and slotted holes.

3. Results and Discussions

After designing and fabricating the motorized cassava grater, experimental work was carried out to investigate the performance characteristics to determine the Grating efficiency, Grating Capacity and Grating Losses of an electrically powered multi-purpose cassava grating machine with grater blades inclined at tooth angles of 25˚ and 30˚.

In Figure 7(a) the perforated plate with inclined blades used for experimental work is shown. On the other hand, Figure 7(b) shows grating operation of the cassava grater during performance testing. In Figure 7(c) weighing scale measuring the mass of grated cassava is illustrated.

(a)

(b)

(c)

Figure 7. (a) Perforated plate, (b) performance testing of cassava grater, (c) mass measurement.

For comparison, performance characteristics of a grating plate whose blades were not inclined at any special angle was investigated. The results are shown in Table 1.

Table 1. Grating time for various masses of cassava and losses for grater blades inclined haphazardly.

Loading Number

Initial mass of un-grated Cassava (kg)

Time taken to grate (sec)

Mass of grated cassava (kg)

Losses (kg)

1

1.25

24.59

1.064

0.186

2

1.20

30.86

1.041

0.209

3

1.05

27.03

0.960

0.090

4

1.17

23.41

0.854

0.316

Total

4.67

105.89

3.919

0.801

Average mass of initial cassava

= Total Mass of initial Cassava Number of Loadings = 4.67 4 =1.1675kg

Average mass of grated cassava

= Total Mass of Grated Cassava Number of Loadings = 3.919 4 =0.97975kg

Average Time Taken to Grate 1 kg

= Total Time Taken to Grate( sec ) Number of Loadings = 105.89 3.919 =27.020seconds

Grating efficiency, ηg

= Mass of well grated cassava Initial mass of cassava tuber ×100%= 0.97975 1.1675 ×100%=83.92%

Grating Capacity

= Mass of total grated cassava Grating time kg hr = 0.97975 0.007353 kg hr =133.2 kg/ hr

Grating Losses

= Initial mass of cassava tuberMass of total grated cassava Initial mass of cassava tuber ×100% = 1.16750.97975 1.1675 ×100%=16.05%

For the case when the blades were not inclined at a specific angle, the results show that it requires an average of 27.02 seconds to grate 1.0 kg of cassava. The grating efficiency when blades are haphazardly inclined was found to be 83.92% and the grating capacity was 133.2 kg/hr.

When the blades were inclined at an angle of 30˚, the measurements of grating time for different masses of cassava as well as losses were recorded in Table 2. It was shown that average grating time was 14.183 seconds for 1.0 kg of cassava when blades are inclined at 30˚ angle. The grating efficiency increased from 83.92% when blades are haphazardly inclined to 93.54%. In addition, it was observed that inclining the grating blades at 30˚ increased the grating capacity from was 133.2 kg/hr to 398.33 kg/hr.

Table 2. Grating time for various masses of cassava and losses for grater blades inclined at 30˚.

Loading Number

Initial mass of un-grated Cassava (kg)

Time taken to grate (sec)

Mass of grated cassava (kg)

Losses (kg)

1

2.00

27.07

1.980

0.020

2

2.09

28.14

1.812

0.278

3

2.05

26.87

1.898

0.152

4

2.03

27.22

1.920

0.110

5

2.05

26.29

1.950

0.100

Total

10.22

135.59

9.560

0.660

Average mass of initial cassava

= Total Mass of Ungrated Cassava Number of Loadings = 10.22 5 = 2.044kg

Average mass of grated cassava

= Total Mass of Grated Cassava Number of Loadings = 9.560 5 =1.912kg

Average Time Taken to Grate 1 kg

= Total Time Taken to Grate( sec ) Number of Loadings = 135.59 9.56 =14.183seconds

Grating efficiency, η g

= Mass of well grated cassava Initial mass of cassava tuber ×100%= 1.912 2.044 ×100%=93.54%

Grating Capacity

= Mass of total grated cassava Grating time kg hr = 1.912 0.0048 kg hr =398.33 kg/ hr

Grating Losses

= Initial mass of cassava tuberMass of total grated cassava Initial mass of cassava tuber ×100% = 2.0441.912 2.044 ×100%=6.458%

When the blades were inclined at an angle of 25˚, the grating time for various masses of cassava and losses were recorded in Table 3.

Table 3. Grating time for various masses of cassava and losses for grater blades inclined at 25˚.

Loading Number

Initial mass of un-grated Cassava (kg)

Time taken to grate (sec)

Mass of grated cassava (kg)

Losses (kg)

1

2.02

27.25

1.947

0.073

2

1.98

26.98

1.775

0.205

3

2.13

28.15

1.869

0.261

4

2.05

27.07

1.863

0.187

5

2.01

26.58

1.924

0.086

Total

10.19

136.03

9.378

0.812

Average mass of un-grated cassava

= Total Mass of initial Ungrated Cassava Number of Loadings = 10.19 5 =2.038kg

Average mass of grated cassava

= Total Mass of Grated Cassava Number of Loadings = 9.378 5 =1.8756kg

Average Time Taken to Grate 1 kg

= Total Time Taken to Grate( sec ) Number of Loadings = 136.03 9.378 =14.052seconds

Grating efficiency, η g

η g = Mass of well grated cassava Initial mass of cassava tuber ×100%= 1.8756 2.038 ×100%=92.03%

Grating Capacity

= Mass of total grated cassava Grating time kg hr = 1.8756 0.004735 kg hr =396.12 kg/ hr

Grating Losses

= Initial mass of cassava tuberMass of total grated cassava Initial mass of cassava tuber ×100% = 2.0381.8756 2.038 ×100%=7.949%

At 25˚ inclination of the grating blades, it was possible to grate 1.0 kg of cassava within an average period of 14.052 seconds. The efficiency increased to 92.03% from 83.92% when no inclination was included on the grating plate. The losses also decreased from 17.15% to 7.949% as a result of inclining the perforations at 25˚. The results obtained from this study are summarized in Table 4.

Table 4. Summary of results.

TOOTH ANGLE

25˚

30˚

Haphazard Angle

Average Mass of Ungrated Cassava (kg)

2.038

2.044

1.1675

Average Mass of Grated Cassava (kg)

1.8756

1.912

0.97975

Average Grating Time (s) for 1 kg

14.052

14.183

27.020

Grating Capacity (kg/hr.)

248.19

253.82

133.2

Grating Efficiency (%)

92.03

93.54

83.92

Grating Losses (%)

7.949

6.458

17.15

From the results in Table 4, it was noted that inclining the perforations on the grating plate at angles of 25˚ and 30˚ has a significant effect on reducing the grating time required to process 1 kg of cassava. This leads to increased grating capacity and efficiency. These findings show that design of cassava grating machine with higher performance characteristics such as grating efficiency, grating capacity and grating Losses, has the potential to accelerate cassava production as a vital means of livelihood for feeding the global population and spur economic development.

4. Conclusions

In this study, an electrically powered multi-purpose cassava grating machine with grater blades inclined at two angles, 25˚ and 30˚ was designed, fabricated and its performance characteristics investigated. From the results of the experiment, the following conclusions were made.

1) The plate with the tooth angle of 30˚ gave better results than that with a tooth angle of 25˚ because the former provided a better grip on the cassava ensuring higher grating efficiency and a reduction in losses.

2) The cassava grating machine will encourage farmers who have not ventured into cassava farming to do so and those practicing in small-scale to expand their acreage due to low post-harvest losses.

3) The machine is also an avenue for job creation for people in and around cassava growing areas for post-harvest processes which include sorting, peeling, cleaning, grating, drying and packaging.

Acknowledgements

The research was supported by the Kenyan German Centre for Mining, Environmental Engineering and Resource Management (CEMEREM) based at Taita Taveta University and also Jomo Kenyatta University of Agriculture & Technology.

Conflicts of Interest

The authors declare no conflict of interest.

References

[1] Burns, A., Gleadow, R., Cliff, J., Zacarias, A. and Cavagnaro, T. (2010) Cassava: The Drought, War and Famine Crop in a Changing World. Sustainability, 2, 3572-3607.
https://doi.org/10.3390/su2113572
[2] Agbetoye, L.A.S. and Oyedele, O.A. (2013) Investigations into Some Engineering Properties of Gari Produced in South—Western Nigeria. International Journal of AgriScience, 3, 728-742.
[3] Ikechukwu, G.A. and Lorreta, E.C. (2015) Development of a Motorized Akpu Milling Machine Design to Improve Poverty Eradication in Nigeria. Proceedings of the World Congress on Engineering and Computer Science 2015, San Francisco, 21-23 October 2015, 945-952.
[4] Bello, S.K., Lamidi, S. and Oshinlaja, S. (2020) Design and Fabrication of Cassava Grating Machine. International Journal of Advances in Scientific Research and Engineering, 6, 162-167
https://doi.org/10.31695/ijasre.2020.33915
[5] Githunguri, C.M., Gatheru, M. and Ragwa, S.M. (2017) Cassava Production and Utilization in the Coastal, Eastern and Western Regions of Kenya.
https://www.researchgate.net/publication/320232890
[6] Hahn, S.K., Mahungu, N.M., Otoo, J.A., Msabaha, M.A.M., Lutaladio, N.B. and Dahniya, M.T. (1986) Cassava and the African Food Crisis, Tropical Root Crops. International Institute of Tropical Agriculture (IITA), 24-26.
[7] Thresh, J.M., Otim Nape, G.W. and Fargette, D. (1997) The Current Pandemic of Cassava Mosaic Virus Disease in East Africa and Its Control. APS Press, 670-677.
[8] Ugoamadi, C.C. (2012) Development of a Cassava Pelleting Machine. Nigerian Journal of Technology (NIJOTECH), 31, 233-240.
[9] Esteves, U., Pantuhan, G.P., Servi, M.O. and Malasador, J.S. (2019) Design, Fabrication and Performance Evaluation of Motor-Operated Cassava Grater. Mindanao Journal of Science and Technology, 17, 227-241.
[10] Okoli, I.G. and Okonkwo, U.C. (2021) Design of an Improved Double Barrel Cassava Grating Machine. International Journal of Engineering Research & Technology (IJERT), 10, 112-117.
[11] Oriaku, E.C., Ossai, E.N. and Agulanna, C.N. (2015) Design and Performance Evaluation of a Double Action Cassava Grating Machine. Journal of Emerging Trends in Engineering and Applied Sciences, 6, 196-203.
[12] Olajide, R.A. and Quadri, A. (2010) Design and Fabrication of an Improved Cassava Grater. The Pacific Journal of Science and Technology, 12, 120-129.
[13] CUTS Africa Resource Centre, Nairobi (2019) Catalysing Private Sector Investment in Kenya’s Cassava Value Chain. CUTS Nairobi.
https://www.cuts-nairobi.org
[14] Peter, K.O., Leo, A.S. and Simeon, O.A. (2010) Cassava Mash Process Handling.
https://www.researchgate.net/publication/272686477
[15] Khurmi, J. (2005) Machine Design. Eurasia Publishing House.
[16] Budynas, R. and Nisbett, K. (2006) Shigley’s Mechanical Engineering Design. 8th Edition, McGraw-Hill.

Copyright © 2025 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.