Open Access

ARTICLE

# The Effects of the Particle Size Ratio on the Behaviors of Binary Granular Materials

School of Civil Engineering, Wuhan University, Wuhan, 430072, China

* Corresponding Author: Xihua Chu. Email:

*Computer Modeling in Engineering & Sciences* **2023**, *136*(1), 63-85. https://doi.org/10.32604/cmes.2023.025062

**Received** 20 June 2022; **Accepted** 14 September 2022; **Issue published** 05 January 2023

## Abstract

The particle size ratio (PSR) is an important parameter for binary granular materials, which may affect the microstructure and macro behaviors of granular materials. However, the effect of particle ratio on granular assemblies with different arrangements is still unclear. To explore and further clarify the effect of PSR in different packing structures, three types of numerical samples with regular, layered, and random packing are designed. Numerical results show that PSR has significant effects on binary granular samples with regular packing. The larger the PSR, the stronger the strength, the larger the modulus, and the smaller the angle between the shear band and the load direction. And a theoretical solution of the peak stress ratio*vs*. PSR is obtained for regular packing, and the results by DEM are in good agreement with the theoretical solution. Under layered packing, PSR has little effect on peak stress ratio due to similar microstructure obtained with the changing of PSR. The modulus slightly increased with the increase of PSR. Under random packing with small grain content of 50%, PSR has little effect in the range of 0.5–0.9, but in a larger range, larger PSR leads to greater modulus.

## Keywords

A granular material generally refers to a collection of a large number of discrete solid particles larger than 1 μ m in size, such as sand and soil. The unique mechanical properties of granular materials have attracted the attention of a large number of researchers. In the past two decades, a considerable number of studies have been carried out to investigate and predict the mechanical behaviors of granular materials [1–6]. It has been recognized that the macroscopic behavior of granular materials is controlled by the microstructure characteristics, such as particle packing [1], particle shape [2], and particle size [3,4].

Binary granular materials are a kind of common granular materials composed of particles of two different sizes. Generally, the size ratio of small particles to large particles (PSR) and the volume fraction of small particles (VS) are key factors affecting the properties of binary granular materials [7–13]. The application of binary granular materials with different packing methods is different. For random packing, binary granular materials can be a simplified model of a well-graded mixture [7]. Many studies have reported the mechanical properties of binary particle mixture with random packing based on Discrete Element Method (DEM) [14–26]. Ogarko et al. [8] found that in the case of isotropic compression, any polydisperse mixtures can be replaced by equivalent bidisperse mixtures when the size distribution moments are matched. Wiącek [23] experimentally and numerically studied the effects of the PSR and VS on the geometrical properties of binary granular mixtures, and found a strong relationship between the coordination number and the VS of binary mixtures. Kumar et al. [9], Zhu et al. [24] and Ng et al. [25] found that the modulus is related to PSR and VS. Ueda et al. [7] explored the PSR and VS of the binary mixture to the shear resistance, and identified the threshold VS, which determines whether the mechanical contribution of large or small particles can be negligible. In addition, they found that the contribution of the large particles disappears when their average spacing with respect to the small particle size is around 2. Nie et al. [26] explored the effects of the PSR and VS on the internal stability of binary mixtures, and demonstrated that the internal stability of binary mixtures generally decreases with decreasing PSR and increases with increasing VS.

For regular packing, binary granular materials are called binary granular crystals or metamaterials. At present, the properties of wave or acoustics in binary granular crystals are mainly focused [12,13], which is widely used in building shock absorption and sound insulation. However, in some engineering applications, such as for building exterior walls, a certain bearing capacity of binary granular metamaterials is required. However, there are few reports on the static properties of binary granular metamaterials. Due to the need for contact between adjacent particles and the limitation of geometric space, the PSR of binary granular metamaterials varies within a small range. In this case, whether PSR is a key factor affecting the behavior of binary granular materials remains to be investigated.

In this study, the strength and deformation characteristics of binary granular materials with different packing structures in a quasi-static state are concerned. Which is of great significance for the design of metamaterials and the study of well-graded mixtures. The PSR ranges from 0.5 to 0.9 in binary granular materials according to the actual situation in binary granular crystals. Binary granular materials with different packing structures are focused on, which are regular packing, layered packing, and random packing. Their quasi-static mechanical properties are investigated based on DEM under biaxial loading in 2-dimensional. The effect of PSR on their strength and deformation characteristics at different confining stress are discussed.

2 Discrete Element Method and Numerical Samples

Discrete element method (DEM) was initiated by Cundall et al. in 1979 [27]. The cores of discrete element method are particle motion equations and contact models between particles.

The motion equation of a particle A can be written according to Newton’s second law as:

where

In this study, Hertz-Mindlin contact law is used to calculate the contact force between particles. As shown in Fig. 1, the contact force is resolved into normal force

The normal force is related to the overlap of two particles:

where

with

where

The calculation of

where

where

In this section, Numerical samples of binary granular materials with three packing structures were established. The size of the samples is about 1.5 m × 2 m. In all samples, the radius of large particles is R = 0.01 m, and the radius of small particles increases from r = 0.5 R to r = 0.9 R. With the decrease in the regularity of the packing structure, three different packing arrangements are generated, which are regular packing, layered packing, and random packing, as shown in Fig. 2. In the regular packing, we adopted the same packing method as Wang et al. [12], filled the small particles in the square array of large particles to generate regular packing structure, as shown in Fig. 2a. When PSR is 1, the packing structure is hexagonal packing with single particle size. The VS for regular binary mixtures are 18.32%, 25.97%, 32.28%, 38.35% and 44.01% for PSR of 0.5, 0.6, 0.7, 0.8 and 0.9, respectively. In the layered packing, the sample is generated by stacking three large particle layers and two small particle layers, and each particle layer results from a hexagonal packing of particles, as shown in Fig. 2b. And the layered samples with different PSR have the same VS of 40%. In the random packing, the samples are generated by random packing of large and small particles with VS of 50%, and samples with different PSR have the same initial porosity of 0.15, as shown in Fig. 2c.

The biaxial test is widely used in the study of granular materials, especially in DEM numerical simulation, as it is easier to visualize particle interactions in two dimensions, and the results are valid for 3-dimensional analysis. The biaxial test is regarded as a simplification of the triaxial test, so a series of biaxial tests for granular samples is performed, as shown in Fig. 3. During biaxial tests based on DEM, firstly rigid plates were manipulated by a servo-controlled mechanism to reach the initial confining pressure and complete the consolidation process. And then the loading was applied by the displacement of the rigid plates at top and bottom, and the rigid plates on the sides keep the confining pressure constant. Loading ends when the target axial strain is reached. The axial strain (

3 Numerical Results and Discussion

The numerical results for each granular sample, including the axial stress-strain curve, axial-volumetric strain curve, and the distribution of the effective strain [28,29] were presented. For the definition of the effective strain, see Appendix A.

A series of DEM simulations were carried out for regular granular samples with different PSR at different confining stress. The deviatoric stress-axial strain curves of regular granular samples with different PSR under different confining stress levels are shown in Figs. 4a–4d, and the relationship between peak stress and PSR is shown in Fig. 5a. It can be seen that with the increase of PSR, the peak deviatoric stress and corresponding axial strain increase. After reaching the peak deviatoric stress, with the increase of axial strain, the deviatoric stress decreases rapidly. The initial modulus increases with increasing PSR, and the peak deviatoric stress and corresponding axial strain increase with the increase of confining pressure stress level. Notice that regular granular samples exhibit strong brittle behavior, a similar phenomenon can be found in the study of O’Sullivan et al. [30], which is related to the rapid destruction of regular granular samples.

For the biaxial tests, there is an analytical solution for the strength of the regular packing structure. Utilizing the condition for sliding at contact points and the symmetry of the packing structure, the analytical solution of the maximum stress ratio was given by Rowe [31]:

where

where

It can be seen from Fig. 5b, the DEM simulation results agree well with the theoretical solution.

Figs. 7a–7d show the relationship between the volumetric strain and the axial strain of the regular granular samples with different PSR under different confining stress. It can be seen that the regular granular samples with different PSR show the phenomenon of first volumetric contraction to volumetric dilation. The peak volumetric contraction increases with the increase of PSR and confining stress level. At the same time, the higher confining stress level and PSR, the higher the axial strain of the regular granular samples at peak volumetric contraction.

Fig. 8 shows the effective strain distribution of regular granular samples at confining stress level of 100 kPa after peak deviatoric stress. Both of them show that shear bands appear after peak deviatoric stress. In the regular packing structure shown in Fig. 6, the angle between the sliding particles and the loading direction is β when the particles slide, which decreases with the increase of PSR. Therefore, with the increase of PSR, the angle between the main shear band and the loading direction decreases. The formation and shape of the shear band have certain randomness, and there is a bifurcation phenomenon.

Layered granular samples whose PSR increases from 0.5 to 0.9 are considered, as shown in Fig. 2b. A series of DEM simulations were carried out at different confining stress. The deviatoric stress-axial strain curves of layered granular samples with different PSR under different confining stress levels are shown in Figs. 9a–9d, and the relationship between peak deviatoric stress and PSR is shown in Fig. 10. It can be seen that the layered granular samples with different PSR show similar deviatoric stress-axial strain curves and peak deviatoric stress. It is because although the granular samples are layered packing with different PSR, each particle layer is regularly hexagonal packing. Similar microstructures result in similar strengths for particle slip failure. Note that the modulus of the layered granular samples increases slightly with increasing PSR, a plausible explanation is the effect of Hertz-Mindlin contact law. According to Eq. (4), the normal contact stiffness is

It can be seen that the normal contact stiffness is related to particle radius, the larger the particle radius, the greater the normal contact stiffness. With the increase of PSR, the particle radius of small particle increases, and the normal contact stiffness between small particles increase, resulting in a slight increase in modulus. To demonstrate this explanation, we use the linear contact model for layered granular samples with the same contact stiffness kn = 3e7 N/m and ks = 1e7 N/m. The layered granular samples are loaded in the same way at the confining stress of 100 kPa. The result is shown in Fig. 11, it can be seen that the layered granular samples with different PSR have the same modulus, which demonstrated the increase in the modulus of layered granular samples is related to the Hertz-Mindlin contact model.

Fig. 12 shows the relationship between the volumetric strain and the axial strain of layered granular samples at different confining stress. It can be observed that similar to the regular granular samples, the volumetric contraction occurs first, and with the increase of axial strain, the volumetric dilatation occurs. Similar volumetric strain-axial strain curves are observed for layered granular samples with different PSR.

Fig. 13 shows the effective strain distribution of layered granular samples at confining stress level of 100 kPa at the axial strain of 2%. Both of them show that shear bands appear. The main shear bands run through several grain layers, and several shear bands occur at the layer boundaries. Although the PSR of layered granular samples is different, the angles between the shear band and the loading direction are similar.

In summary, the PSR has little influence on the strength and deformation characteristics of layered granular samples. It may be because although the PSR of the layered granular samples is different, they have very similar microstructures, resulting in similar strength and deformation properties.

Random granular samples with the same porosity whose PSR increases from 0.5 to 0.9 are considered, as shown in Fig. 2c. And VS of random granular samples is 50%. A series of DEM simulations were carried out at different confining stress. The deviatoric stress-axial strain curves of random granular samples with different PSR under different confining stress levels are shown in Figs. 14a–14d. It can be seen that the particle assemblies show strain softening. The deviatoric stress tends to stabilize with the increase of axial strain.

Fig. 15 shows the volumetric strain-axial strain curves of random granular samples at different confining stress. Similar to the particle assemblies with the regular and layered arrangement, the volumetric contraction occurs first, and with the increase of axial strain, the volumetric dilatation occurs. Fig. 16 shows the effective strain distribution of randomly arranged particle assemblies at confining stress level of 100 kPa at the axial strain of 15%. Compared with particle assemblies with the regular and layered arrangement, particle assemblies with the random arrangement also show that shear bands appear, but the locations of shear bands are random and not very concentrated.

It can be seen that within the range of 0.5–0.9, PSR has little influence on the mechanical and deformation properties of random granular samples. It may be related to the similar microstructure. As indicated by Zhu et al. [24] and Gong et al. [32], Young’s modulus and the shear modulus are strongly linked to the coordination number. In this study, the coordination number

where

Although the mechanical and deformation properties of random granular samples are similar within such a small range of PSR, the contact behavior is changed to some extent. For binary mixture with random packing, the contacts in a particle assembly can be classified as big particle-big particle (

where

Within a limited range of 0.5–0.9, the microstructure changes caused by PSR are very limited. We generated random granular samples with a larger range of PSR with VS of 50% in the same way, in which PSR is 1/2, 1/4, 1/8, and 1/16, and performed the same loadings at confining stress of100 kPa. If the radius of large particles is fixed, the PSR will cause a large computational cost when the PSR is relatively small, so we fix the radius of small particles to 0.008 m. The deviatoric stress-axial strain curves of random granular samples with different PSR at confining stress of 100 kPa are shown in Fig. 18. It can be seen that the modulus of random granular samples increases with the decrease of PSR, which is consistent with the result obtained by Zhu et al. [24]. They also discussed the contribution of contact type to modulus and found that the

The effect of PSR on the strength and deformation of binary granular materials with different packing methods was studied. PFC2D was used to simulate the influence of PSR on basic mechanical properties such as shear strength, volumetric strain, and so on. Besides, the effective strain distribution and the fabric tensors of the granular materials are analyzed. The main conclusions are as follows:

1) The macroscopic mechanical properties of binary granular materials with regular packing are greatly affected by the PSR. It is mainly because the effect of PSR mainly comes from the changes in microstructure induced by PSR. The theoretical formula for the peak stress ratio and PSR is obtained, which is in good agreement with the DEM simulation. The strength and modulus of binary granular materials with regular packing increase with increasing PSR, which is mainly due to the change in the structural slip angle. Clear shear bands appear in binary granular materials with regular packing after peak deviatoric stress. With the increase of PSR, the angle between the main shear band and the loading direction decreases.

2) The PSR has little influence on the macroscopic mechanical properties of binary granular materials with layered packing, which is caused by similar microstructure. The modulus of binary granular materials with layered packing increases with increasing PSR, and we suspect that it is due to the Hertz-Mindlin contact law used in the DEM simulations, which results in greater contact stiffness with an increase of PSR.

3) For random binary granular materials with VS of 50%, PSR has little effect on mechanical and deformation properties within the range of 0.5–0.9. In the larger range of PSR, with the increase of PSR, the modulus of the random binary granular materials increases. Strong force chains are mainly distributed in

Funding Statement: This work was supported by the National Natural Science Foundation of China (Nos. 12172263, 11772237).

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

## References

1. Chang, C. S., Misra, A. (1990). Packing structure and mechanical-properties of granulates. Journal of Engineering Mechanics, 116(5), 1077–1093. DOI 10.1061/(ASCE)0733-9399(1990)116:5(1077). [Google Scholar] [CrossRef]

2. Al-Raoush, R. (2007). Microstructure characterization of granular materials. Physica A: Statistical Mechanics and Its Applications, 377(2), 545–558. DOI 10.1016/j.physa.2006.11.090. [Google Scholar] [CrossRef]

3. Hare, C. L., Ghadiri, M. (2013). Influence of measurement cell size on predicted attrition by the distinct element method. Powder Technology, 236(1), 100–106. DOI 10.1016/j.powtec.2012.04.061. [Google Scholar] [CrossRef]

4. Obregon, L., Realpe, A., Velazquez, C. (2010). Mixing of granular materials, part II: Effect of particle size under periodic shear. Powder Technology, 201(3), 193–200. DOI 10.1016/j.powtec.2010.03.006. [Google Scholar] [CrossRef]

5. Wang, J., Chu, X., Wang, J. (2022). The material deformation and internal structure development of granular materials under different cyclic loadings. Computer Modeling in Engineering & Sciences, 130(2), 653–670. DOI 10.32604/cmes.2022.018207. [Google Scholar] [CrossRef]

6. Qu, T., Di, S., Feng, Y. T., Wang, M., Zhao, T. et al. (2021). Deep learning predicts stress-strain relations of granular materials based on triaxial testing data. Computer Modeling in Engineering & Sciences, 128(1), 129–144. DOI 10.32604/cmes.2021.016172. [Google Scholar] [CrossRef]

7. Ueda, T., Matsushima, T., Yamada, Y. (2011). Effect of particle size ratio and volume fraction on shear strength of binary granular mixture. Granular Matter, 13(6), 731–742. DOI 10.1007/s10035-011-0292-1. [Google Scholar] [CrossRef]

8. Ogarko, V., Luding, S. (2012). Equation of state and jamming density for equivalent bi- and polydisperse, smooth, hard sphere systems. Journal of Chemical Physics, 136(12), 124508. DOI 10.1063/1.3694030. [Google Scholar] [CrossRef]

9. Kumar, N., Magnanimo, V., Ramaioli, M., Luding, S. (2016). Tuning the bulk properties of bidisperse granular mixtures by small amount of fines. Powder Technology, 293, 94–112. DOI 10.1016/j.powtec.2015.11.042. [Google Scholar] [CrossRef]

10. Aarons, L. R., Balachandar, S., Horie, Y. (2013). The mixing of cohesive granular materials featuring a large size range in the absence of gravity. Powder Technology, 235(2), 18–26. DOI 10.1016/j.powtec.2012.09.044. [Google Scholar] [CrossRef]

11. Shire, T., O’Sullivan, C., Hanley, K. J. (2016). The influence of fines content and size-ratio on the micro-scale properties of dense bimodal materials. Granular Matter, 18(3), 52. DOI 10.1007/s10035-016-0654-9. [Google Scholar] [CrossRef]

12. Wang, J., Chu, X. H., Zhang, J. B., Liu, H. (2019). The effects of microstructure on wave velocity and wavefront in granular assemblies with binary-sized particles. International Journal of Solids and Structures, 159(2), 156–162. DOI 10.1016/j.ijsolstr.2018.09.026. [Google Scholar] [CrossRef]

13. Wang, J., Chu, X. H., Jiang, Q. H., Xiu, C. X. (2019). Energy transfer and influence of excitation frequency in granular materials from the perspective of Fourier transform. Powder Technology, 356(2), 493–499. DOI 10.1016/j.powtec.2019.08.061. [Google Scholar] [CrossRef]

14. Wang, D. F., An, X. Z., Gou, D. Z., Zhao, H. Y., Wang, L. et al. (2018). DEM construction of binary hard sphere crystals and radical tessellation. AIP Advances, 8(10), 105203. DOI 10.1063/1.5052478. [Google Scholar] [CrossRef]

15. Martin, C. L., Bouvard, D. (2004). Isostatic compaction of bimodal powder mixtures and composites. International Journal of Mechanical Sciences, 46(6), 907–927. DOI 10.1016/j.ijmecsci.2004.05.012. [Google Scholar] [CrossRef]

16. Minh, N. H., Cheng, Y. P., Thornton, C. (2014). Strong force networks in granular mixtures. Granular Matter, 16(1), 69–78. DOI 10.1007/s10035-013-0455-3. [Google Scholar] [CrossRef]

17. Nordstrom, J., Alderborn, G., Frenning, G. (2018). Compressibility and tablet forming ability of bimodal granule mixtures: Experiments and DEM simulations. International Journal of Pharmaceutics, 540(1–2), 120–131. DOI 10.1016/j.ijpharm.2018.02.006. [Google Scholar] [CrossRef]

18. Gong, J., Liu, J. (2017). Mechanical transitional behavior of binary mixtures via DEM: Effect of differences in contact-type friction coefficients. Computers and Geotechnics, 85(4), 1–14. DOI 10.1016/j.compgeo.2016.12.009. [Google Scholar] [CrossRef]

19. Cai, R. H., Xiao, H. Y., Zheng, J. Y., Zhao, Y. Z. (2019). Diffusion of size bidisperse spheres in dense granular shear flow. Physical Review, E, 99(3), 032902. DOI 10.1103/PhysRevE.99.032902. [Google Scholar] [CrossRef]

20. Tunuguntla, D. R., Weinhart, T., Thornton, A. R. (2017). Comparing and contrasting size-based particle segregation models: Applying coarse-graining to perfectly bidisperse systems. Computational Particle Mechanics, 4(4), 387–405. DOI 10.1007/s40571-016-0136-1. [Google Scholar] [CrossRef]

21. Staron, L., Phillips, J. C. (2016). How large grains increase bulk friction in bi-disperse granular chute flows. Computational Particle Mechanics, 3(3), 367–372. DOI 10.1007/s40571-015-0068-1. [Google Scholar] [CrossRef]

22. Tunuguntla, D. R., Thornton, A. R., Weinhart, T. (2016). From discrete elements to continuum fields: Extension to bidisperse systems. Computational Particle Mechanics, 3(3), 349–365. DOI 10.1007/s40571-015-0087-y. [Google Scholar] [CrossRef]

23. Wiącek, J. (2016). Geometrical parameters of binary granular mixtures with size ratio and volume fraction: experiments and DEM simulations. Granular Matter, 18(3), 1–10. DOI 10.1007/s10035-016-0642-0. [Google Scholar] [CrossRef]

24. Zhu, Y. G., Gong, L., Nie, Z. H. (2020). Numerical investigation of the elastic properties of binary mixtures as a function of the size ratio and fines content. International Journal of Geomechanics, 20(9), 04020155. DOI 10.1061/(ASCE)GM.1943-5622.0001792. [Google Scholar] [CrossRef]

25. Ng, T. T., Zhou, W., Ma, G., Chang, X. L. (2018). Macroscopic and microscopic behaviors of binary mixtures of different particle shapes and particle sizes. International Journal of Solids and Structures, 135(3), 74–84. DOI 10.1016/j.ijsolstr.2017.11.011. [Google Scholar] [CrossRef]

26. Nie, Z., Zhu, Y., Zou, J., Gong, J., Liu, S. (2019). DEM study of the microscopic characteristics and internal stability of binary mixtures. Powder Technology, 352, 314–324. DOI 10.1016/j.powtec.2019.04.077. [Google Scholar] [CrossRef]

27. Cundall, P. A., Strack, O. D. L. (1979). A discrete numerical mode for granular assemblies. Géotechnique, 29(1), 47–65. DOI 10.1680/geot.1979.29.1.47. [Google Scholar] [CrossRef]

28. Li, X. K., Chu, X. H., Feng, Y. T. (2005). A discrete particle model and numerical modeling of the failure modes of granular materials. Engineering Computations, 22(8), 894–920. DOI 10.1108/02644400510626479. [Google Scholar] [CrossRef]

29. Tang, H. X., Dong, Y. F., Chu, X. H., Zhang, X. (2016). The influence of particle rolling and imperfections on the formation of shear bands in granular material. Granular Matter, 18(1), 12. DOI 10.1007/s10035-016-0607-3. [Google Scholar] [CrossRef]

30. O’Sullivan, C., Bray, J. D., Riemer, M. (2004). Examination of the response of regularly packed specimens of spherical particles using physical tests and discrete element simulations. Journal of Engineering Mechanics, 130(10), 1140–1150. DOI 10.1061/(ASCE)0733-9399(2004)130:10(1140). [Google Scholar] [CrossRef]

31. Rowe, P. W. (1962). The stress-dilatancy relation for static equilibrium of an assembly of particles in contact. Proceedings of the Royal Society of London A, 269(1339), 500–527. DOI 10.1098/rspa.1962.0193. [Google Scholar] [CrossRef]

32. Gong, J., Wang, X., Li, L., Nie, Z. (2019). DEM study of the effect of fines content on the small-strain stiffness of gap-graded soils. Computers and Geotechnics, 112(124), 35–40. DOI 10.1016/j.compgeo.2019.04.008. [Google Scholar] [CrossRef]

33. Thornton, C. (2000). Numerical simulations of deviatoric shear deformation of granular media. Géotechnique, 50(1), 43–53. DOI 10.1680/geot.2000.50.1.43. [Google Scholar] [CrossRef]

34. Shaebani, M. R., Madadi, M., Luding, S., Wolf, D. E. (2012). Influence of polydispersity on micromechanics of granular materials. Physical Review E, 85(1), 011301. DOI 10.1103/PhysRevE.85.011301. [Google Scholar] [CrossRef]

Appendix A. The Definition of the Effective Strain

To describe the change of relative position between particles, Li et al. [28] and Tang et al. [29] defined the nominal strain based on the change of particle position. Fig. A1 shows the position change of particle A and one of its neighboring particles B.

The coordinates of particles A and B at time

In the local coordinate system (

where

Refer to local coordinate system (

where

where

where

Then the displacement derivative matrix can be defined as

where

where

## Cite This Article

**APA Style**

*Computer Modeling in Engineering & Sciences*,

*136*

*(1)*, 63-85. https://doi.org/10.32604/cmes.2023.025062

**Vancouver Style**

**IEEE Style**

*Comput. Model. Eng. Sci.*, vol. 136, no. 1, pp. 63-85. 2023. https://doi.org/10.32604/cmes.2023.025062