SciELO - Scientific Electronic Library Online

 
vol.90 número226Sostenibilidad y futuro de la industria del petróleo y gas: una mini-revisiónEstudio sobre la ley de variación del estrés del carbón de concreción dura bajo la acción del pico basado en DYNA índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • En proceso de indezaciónCitado por Google
  • No hay articulos similaresSimilares en SciELO
  • En proceso de indezaciónSimilares en Google

Compartir


DYNA

versión impresa ISSN 0012-7353versión On-line ISSN 2346-2183

Dyna rev.fac.nac.minas vol.90 no.226 Medellín abr./jun. 2023  Epub 13-Feb-2024

https://doi.org/10.15446/dyna.v90n226.106462 

Article

Fatigue evaluation of asphalts highly modified with recycled tire rubber

Evaluación a fatiga de asfaltos altamente modificados con polvo de neumático fuera de uso

aLEMaC Centro de Investigaciones Viales UTN - FRLP CIC PBA, Universidad Tecnológica Nacional, La Plata, Argentina. izapata@frlp.utn.edu.ar, orebollo@frlp.utn.edu.ar, jrivera@frlp.utn.edu.ar, gbotasso@frlp.utn.edu.ar


Abstract

Fatigue resistance is one of the most frequent and most complex problems to predict in asphalt mixtures. Due to its viscoelastic nature, asphalt is the material that presents the greatest contribution to fatigue within the asphalt mixture. The G*sinδ parameter proposed by SUPERPAVE is the most widely used for fatigue evaluation; however, it has its limitations since it does not allow to observe if the sample is damaged, for example.

The results of the present work allow to demonstrate the benefit of the incorporation of the recycled tire rubber as a modifier residue of the asphalt, through the analysis with the rheometer through the proposed test. Finally, the cycle curves are made based on the imposed deformations, which show a greater resistance to fatigue in the asphalt, with increasing rates of tire rubber out of use.

Keywords: asphalt; fatigue; recycled tire rubber; rheology

Resumen

La resistencia a fatiga es una de las problemáticas más frecuentes y más compleja de predecir en las mezclas asfálticas. Por su carácter viscoelástico, el asfalto es el material que presenta mayor aporte a la fatiga dentro de la mezcla asfáltica. El parámetro G*senδ propuesto por SUPERPAVE es el más ampliamente utilizado para la evaluación de fatiga; sin embargo, tiene sus limitaciones ya que no permite observar si la muestra está dañada, por ejemplo.

Los resultados del presente trabajo permiten evidenciar el beneficio de la incorporación del neumático fuera de uso como un residuo modificador del asfalto, mediante el análisis con el reómetro por medio del ensayo planteado. Por último, se confeccionan las curvas de ciclos en función de las deformaciones impuestas, las cuales permiten evidenciar una mayor resistencia a fatiga en el asfalto, con crecientes tasas de caucho de neumático fuera de uso.

Palabras clave: asfalto; fatiga; neumático fuera de uso; reología

1. Introduction

Fatigue in flexible pavements is one of the main causes of reduced service life. The cracking initiated in the lower part of the asphalt mixture is due to reaching a limit of the value of the specific deformation of the asphalt layer [1,2]. The damage caused by fatigue is caused by a combination of the structural design of the package used and the selection of intervening materials. In turn, it has a great impact, the thickness of the pavement that comes into play, as well as the volumetric properties of the mixture used. Within the asphalt mix, asphalt has a significant role in fatigue failure [3,4]. Studies carried out by means of fracture mechanics and microsimulation indicate that fatigue cracking initially propagates through the asphalt binder, then passing to the mix matrix.

The methodologies proposed by SHRP (Strategic Highway Research Program), such as the evaluation of the PG grade for the classification of asphalts, were useful in the 1990, but currently the advances in the studies show that several criteria do not reflect the reality of the materials. The rheological characterization of asphalts in the dynamic shear rheometer (DSR) evaluates the material in the linear viscoelastic region, both for rutting and fatigue parameters in the guidelines set by SHRP. In particular, the fatigue damage suffered by asphalt pavements is strongly influenced by the nature of the asphalt bitumen that composes it, said failure occurs outside the linear viscoelastic zone. Due to this, the traditional methods proposed by SUPERPAVE, such as G*sinδ, do not fully evaluate this problem, since the analysis carried out is absent from strain stress variations, which translate to only one load state evaluated in the rheometer [5-7].

Based on these limitations, the time sweep test in the dynamic shear rheometer is a more accurate way to quantify the fatigue in the asphalt binder. This methodology contemplates the evaluation of the complex modulus of the sample at a defined frequency, amplitude and temperature. The failure criterion of the material according to this test is proposed as a general criterion when the complex modulus drops to 50 % of the initial modulus. By repeating the procedure at different amplitudes, the relationship between the imposed deformation and the failure cycles of the material is established, which is represented by eq. (1).

Where Nf is the number of cycles to failure, γ the specific shear strain, A and B the model fit parameters that depend on the material. This test scheme and representation model is the one that is frequently used to evaluate fatigue in asphalt mixtures by means of the four-point beam.

However, despite being a proposal that surpasses the fatigue failure methodology in asphalt proposed by SUPERPAVE, the time sweep test in the DSR implies long test times for the material to be evaluated, depending on the deformation applied to the surface. sample [2]. On the other hand, it implies a test for each point of specific strain applied, resulting in many samples to be able to conform the failure curve presented above. In turn, the reduction to 50 % of the initial modulus as a material failure criterion is arbitrary since there is still less than half of the useful life of the material to be used [3].

1.1. Linear Amplitude Sweep (LAS)

In recent years, in search of an alternative to the methodologies described above, the linear amplitude sweep test has been proposed to evaluate fatigue in asphalt cements. This method considers evaluating the sample with an increasing deformation stress to simulate the state of solicitations that are generated in a pavement.

The results of the linear amplitude scan are applicable to the Simplified Theory of Viscoelastic Continuous Damage (S-VECD), [5,6,8] which models the internal work of the material and quantifies the damage by means of microstructural changes in the sample [9, 10]. Based on the study carried out by Schapery and Park [11], the damage accumulation law is reached, represented by eq. (2).

Where S is the damage growth in the sample, W the work done, t the time, and α the sample integrity parameter. This theory is recently applied to studies to predict fatigue in both asphalts and asphalt mixtures [4,9].

The sample integrity parameter α was previously estimated from the curve of relaxation modulus (G*(t)) as a function of time. To obtain it, a conversion of the G* module from the frequency domain to the time domain has to be carried out. This conversion is not used regularly due to the complexity involved in its application. To simplify obtaining α, it is possible to quantify it by means of the slope of the frequency sweep test, in which the sample is considered without damage [5]. Finally, it is expressed as α=1/m, where m is the slope of the frequency sweep test, with which the sample without damage is characterized. On the other hand, the work done by the sample can be modeled according to eq. (3).

Where W is the work done expressed as energy dissipated, γ0 the specific applied shear strain, G* the complex shear modulus and δ the phase angle. In this way, eq. (2) can be integrated numerically using eq. (3) in order to arrive at the expression for damage growth, shown in eq. (4).

Where Cj is the ratio between the modulus at each point and the initial modulus of the material without damage, γ0 is the applied shear strain, and t is time.

To quantify the damage presented by the samples during the test, the C-S damage curves are plotted, which represent the integrity of the sample as a function of the accumulated damage. Said curve can be adjusted by means of eq. (5).

Where C0 =1 and C1 and C2 are the fit coefficients of the model. Differentiating eq. (5) and substituting it into eq. (2), we can finally arrive at eq. (1), which in its extended form appears as eq. (6).

Where f is the test frequency (defined at 10 Hz), Df is the failure criterion adopted, k = 1+ (1- C2) α, and B = -2α. Based on the study carried out by Masad [12], it is considered that the specific strain provided by the asphalt represents between 12 and 90 times the specific strain of the asphalt mixture. Finally, Bahía sets this parameter at 50 times and proposes a specific shear deformation of 2,5 % for pavements greater than 10 cm and 5,0% for pavements less than 10 cm [13-15].

Several authors propose different failure criteria Df [1,2,7,9,16] to define the fatigue limit of asphalts contemplating this test. The present work proposes to consider the two most used criteria to obtain the resulting equivalent axis curves. In the first place, the drop in the value of the C-S curve by 35 % from the initial value, data that unifies the position for the evaluated binders. The second corresponds to the peak value of the shear stress exerted on the sample evaluated in the C-S curve [15].

2. Materials and methods

For the study, an asphalt cement classified by viscosity as CA-30 and recycled tire rubber powder (RTR, Fig. 1) is used as asphalt modifier. The RTR is added to the asphalt in percentages of 15 %, 20 % and 25 %, according to the particle size of the sample as can be seen in Fig. 2. The modification of the asphalts is carried out in a high-speed disperser for one hour to obtain a homogeneous sample. Base asphalt and RTR-modified asphalts are subjected to a thin-film oven (RTFO) to simulate the short-term aging conditions that occur in asphalt. Asphalts are characterized by rotational viscosity, torsional elastic recovery and softening point tests before and after aging the samples.

Source: Prepared by the author.

Figure 1 Recycled tire rubber. 

Source: Prepared by the author.

Figure 2 Particle Size distribution of RTR 

Source: Prepared by the author.

Figure 3 Dynamic shear rheometer (left), 8 mm geometry with the sample (right) 

For the fatigue characterization of asphalts by means of the LAS test, a Discovery HR-2 Hybrid Rheometer brand DSR from TA Instruments is used. It has an ambient test chamber (ETC, for its acronym for Environmental Test Chamber) as shown in Fig. 3. To carry out the test, the geometry of parallel plates with 8 mm diameter is used, with a gap of 2 mm between them.

The LAS test is regulated by the AASHTO 101-14 standard [17,18] to characterize fatigue damage in asphalt. However, to date there is no asphalt classification specification that has incorporated this evaluation [19], as if the Multiple Stress Creep and Recovery (MSCR) test is present for the evaluation of rutting in asphalt bitumens.

The methodology includes two stages within the trial itself. The first consists of characterizing the sample by means of a frequency sweep in the range of 0,2 to 30 Hz at a specific shear strain of 0,1 %. This stage is used to characterize the material without damage. The second part lies in subjecting the sample to a state of increasing deformations. This stage is carried out at a frequency of 10 Hz and at a constant increase in specific strains from 0,1 % to 30 %. At each load cycle, the peak stress associated with the imposed strain is recorded, along with the resulting phase angle. The standardized procedure proposes the increase of deformations in a time of 310 seconds, finding in the bibliography alternatives for this load time. The test temperature at which the samples are evaluated is 20 °C, 25 °C and 30 °C.

With the data collected from the tests, the C-S curves are drawn for each of the asphalts studied and for each test temperature. Subsequently, the fatigue life curves are evaluated at the temperatures considered and, finally, a comparison is made between the two failure criteria stated for the asphalts evaluated at 20 °C.

3. Results and discusion

3.1. Damage C-S Curves

The asphalt binders used are a conventional asphalt classified by viscosity (CA-30) and three binders modified from the base asphalt with increasing rates of rubber from RTR (15%RTR, 20%RTR and 25%RTR). Table 1 shows the properties of the asphalts used.

Table 1 Characterization of the asphalts used 

S.P.: Softening Point

TER: torsional elastic recovery

Source: Prepared by the author

Table 2 Parameters of the model for the asphalts evaluated. 

Source: Prepared by the author.

After aging the asphalts in the RTFO, the LAS test is carried out in the DSR. Through the previously performed analysis, the parameters of the model are calculated according to the S-VECD to evaluate the coefficient α of integrity of the material and the coefficients C1 and C2 of the damage curve. Table 2 shows the tests at 20 °C, 25 °C and 30 °C below.

In the first instance, an analysis of the alpha parameter is carried out, which indicates the integrity of the sample without damaging it. Said parameter strongly depends on the stiffness of the sample, increasing the value of α when the stiffness increases. Therefore, the value of α shows a first indication of fatigue resistance of the samples, increasing with the increase of RTR as can be seen in Fig. 4.

Source: Prepared by the author.

Figure 4 Parameter α as function percentage RTR 

With the adjustment parameters found, the C-S curves of material integrity are plotted as a function of the accumulated damage. Figs. 5-7 show the curves for 20 °C, 25 °C and 30 °C, respectively. For the three temperatures, it can be seen in the graphs that the asphalts modified with RTR maintain a greater integrity of the material for damage intensities greater than 50. In the same way, a superior performance is evidenced when the percentage of incorporation of RTR increases, although in a less pronounced. It can also be shown that the asphalts with RTR have a 60 % material life based on the damage intensity limit value for the unmodified asphalt, for the evaluations at 20 °C.

Source: Prepared by the author.

Figure 5 C-S curve at 20 °C for the evaluated asphalts. 

Source: Prepared by the author.

Figure 6 C-S curve at 25 °C for the evaluated asphalts. 

Source: Prepared by the author.

Figure 7 C-S curve at 30 °C for the evaluated asphalts. 

3.2. Fatigue laws

From the developed C-S graphs, the models are assembled to estimate the fatigue life based on the maximum specific deformation of the material and the failure criterion adopted for its evaluation. In the first instance, these curves are evaluated with the 35 % reduction of the initial integrity to develop the model seen in eq. (1). As can be seen, in principle, for 20 °C there is a considerable benefit from the incorporation of RTR as shown. see in Figs. 8-10.

3.3. Comparison between failure criteria contemplated

As previously mentioned, based on the tests carried out, the comparison between failure criteria is proposed to obtain the fatigue life curves. With the curves already exposed with a failure criterion of reduction of 35 % of the initial integrity, the fatigue life curves are formed with the criterion of the maximum shear stress peak, together with its corresponding value in the C-S curve. Fig. 11 shows comparatively the asphalts evaluated at 20 °C and with an applied deformation of 2,5 %.

The graph at 25 °C follows the same trend, however, a difference is evident between the different proportions of rubber incorporated. The fatigue phenomenon manifests itself at intermediate temperatures (approximately 10 °C - 30 °C).

Source: Prepared by the author.

Figure 8 Fatigue life at 20 °C for the asphalts evaluated with Df = 35% damage. 

Source: Prepared by the author.

Figure 9 Fatigue life at 25 °C for the asphalts evaluated with Df = 35% damage. 

Source: Prepared by the author.

Figure 10 Fatigue life at 30 °C for the asphalts evaluated with Df = 35% damage. 

Source: Prepared by the author.

Figure 11 Fatigue life at 20 °C for different failure criteria (2,5 % shear strain) 

Therefore, as it is tested at a higher temperature, improvements in performance are evidenced and, in turn, materials with a higher degree of modification present better performance. In this last graph you can see an improvement in the incorporation of 25 % compared to the other alternatives.

As can be seen, there is a significant difference between the proposed criteria, a gap that widens with the growth of RTR incorporation. On the other hand, the 35 % criterion, although it is conservative, and is not something minor when evaluating fatigue in pavements, is a fixed criterion for all types of asphalt. In this way, for asphalts modified with RTR, the life of the material is being underestimated, since it is considered to be in a state of failure, when its remaining life is considerable. The criterion of reduction of an initial value is equivalent to the criterion of the reduction of 50 % of the initial modulus to evaluate asphalt mixtures by means of the fatigue beam at four points. As previously mentioned, these criteria are being questioned due to the randomness when it comes to pinpointing a material failure point.

On the other hand, although the life data considering the peak shear stress are much higher than those evaluated by the 35% reduction, it is a criterion that reveals the variation of the failure point of the material and that overestimates the useful life of the material.

Source: Prepared by the author.

Figure 12 Shear stress as a function of the specific deformation at 20°C. 

As can be seen in Fig. 12, despite obtaining a higher peak stress, the CA-30 sample reaches this value prematurely and falls sharply compared to asphalts with RTR. The modified asphalts show a more stable curve as the percentage of RTR increases. Said effect is contemplated in the C-S curves that evaluate the internal work carried out by the material. On the other hand, the displacement of the stress peaks in the samples translates into large displacements in the C-S curves. Finally, this has a considerable effect when evaluating the number of axes considered to be fatigued.

4. Conclusions

In the present work, the characterization of four asphalts (one conventional and three modified with increasing contents of recycled tire rubber) aged in the rotary thin-film oven is carried out, by means of the fatigue test in a dynamic shear rheometer, known as sweep of linear amplitude, at 20 °C, 25 °C and 30 °C. This test allows to print a considerable damage to the samples, applying deformations of up to 30 %, effectively leading to fatigue. The test methodology is consistent with the simplified theory of continuous viscoelastic damage and allows modeling the C-S curves (Cumulative Damage versus Material Integrity) and the fatigue life curves to estimate the number of axes that the material is capable of resisting (Nf).

For the preparation of these curves, two commonly used failure criteria were studied. On the one hand, the damage corresponding to a 35 % drop in the damage curve, and on the other hand, the value of the damage curve resulting from the shear stress peak during the test. This last criterion converges to fatigue life results that allow evaluating the material with a higher performance than that considered by the 35 % reduction of the C-S curve.

References

[1] Kröger, I., and Pfeiff, C., Ensayos ágiles en DSR para estimar desempeño de mezclas asfálticas. Conference Paper, December 2020, 13 P. [ Links ]

[2] Saboo, N., New damage parameter for fatigue analysis of asphalt binders in linear amplitude sweep test. Journal of Materials in Civil Engineering, 32(6), art. 04020126, 2020. DOI: https://doi.org/10.1061/(ASCE)MT.1943-5533.0003225Links ]

[3] Wang, C., Castorena, C., Zhang, J., and Richard-Kim, Y., (). Unified failure criterion for asphalt binder under cyclic fatigue loading. Road Materials and Pavement Design, 16(sup2), pp. 125-148, 2015. DOI: https://doi.org/10.1080/14680629.2015.1077010Links ]

[4] Zhang, F., Wang, L., Li, C., and Xing, Y., Predict the phase angle master curve and study the viscoelastic properties of warm mix crumb rubber-modified asphalt mixture. Materials, 13(21), art. 5051, 2020. DOI: https://doi.org/10.3390/ma13215051Links ]

[5] Hintz, C., Velasquez, R., Johnson, C., and Bahia, H., (). Modification and validation of linear amplitude sweep test for binder fatigue specification. Transportation Research Record, 2207(1), pp. 99-106, 2011. DOI: https://doi.org/10.3141/2207-13Links ]

[6] Hintz, C., and Bahia, H., Simplification of linear amplitude sweep test and specification parameter. Transportation Research Record, 2370(1), pp. 10-16, 2013. DOI: https://doi.org/10.3141/2370-02Links ]

[7] Yue, M., Yue, J., Wang, R., and Xiong, Y., Evaluating the fatigue characteristics and healing potential of asphalt binder modified with Sasobit® and polymers using linear amplitude sweep test. Construction and Building Materials, 289, art. 123054, 2021. DOI: https://doi.org/10.1016/j.conbuildmat.2021.123054Links ]

[8] Klug, A., Ng, A., and Faxina, A., Application of the Viscoelastic continuum damage theory to study the fatigue performance of asphalt mixtures-A literature review. Sustainability, 14(9), art. 4973, 2022. DOI: https://doi.org/10.3390/su14094973Links ]

[9] Safaei, F., Lee, J.S., Nascimento, L.A.H.D., Hintz, C., and Kim, Y.R., Implications of warm-mix asphalt on long-term oxidative ageing and fatigue performance of asphalt binders and mixtures. Road Materials and Pavement Design, 15(sup1), pp. 45-61, 2014. DOI: https://doi.org/10.1080/14680629.2014.927050Links ]

[10] Chen, H., Zhang, Y., and Bahia, H.U., Estimating asphalt binder fatigue at multiple temperatures using a simplified pseudo-strain energy analysis approach in the LAS test. Construction and Building Materials, 266, art. 120911, 2021. DOI: https://doi.org/10.1016/j.conbuildmat.2020.120911Links ]

[11] Schapery, R.A., and Park, S.W., Methods of interconversion between linear viscoelastic material functions. Part II-An approximate analytical method. International Journal of Solids and Structures, 36(11), pp. 1677-1699, 1999. DOI: https://doi.org/10.1016/S0020-7683(98)00060-2Links ]

[12] Masad, E., Somadevan, N., Bahia, H.U., and Kose, S., Modeling and experimental measurements of strain distribution in asphalt mixes. Journal of Transportation Engineering, 127(6), pp. 477-485, 2001. DOI: https://doi.org/10.1061/(ASCE)0733-947X(2001)127:6(477)Links ]

[13] Chen, H., Zhang, Y., and Bahia, H.U., Modelling asphalt binder fatigue at multiple temperatures using complex modulus and the LAS test. International Journal of Pavement Engineering, 23(13), pp. 4600-4609. 2022. DOI: https://doi.org/10.1080/10298436.2021.1968394Links ]

[14] Mandal, T., and Bahia, H.U., Effect of crumb rubber on rheological properties of asphalt binder and aggregate packing of asphalt mixtures. Conference Paper, Rubberized Asphalt - Asphalt Rubber, October 2015. [ Links ]

[15] Foroutan-Mirhosseini, A., Kavussi, A., Jalal-Kamali, M.H., Khabiri, M.M., and Hassani, A., Evaluating fatigue behavior of asphalt binders and mixes containing date seed ash. Journal of Civil Engineering and Management, 23(8), pp. 1164-1175, 2017. DOI: https://doi.org/10.3846/13923730.2017.1396560Links ]

[16] Coutinho, R., Castelo-Branco, V.T.F., Babadopulos, L.F.A.L., and Soares, J., (). The use of linear amplitude sweep tests to characterize fatigue damage in fine aggregate matrices. In 1st Conference on Rheology and Processing of Construction Materials. Paris, France, 2013. [ Links ]

[17] AASHTO TP101., Estimating damage tolerance of asphalt binders using the linear amplitude sweep. 2014. [ Links ]

[18] AASHTO TP101., Estimating damage tolerance of asphalt binders using the linear amplitude sweep. 2012. [ Links ]

[19] Chen, H., and Bahia, H.U., Proposed asphalt binder fatigue criteria for various traffic conditions using the LAS or the GR parameters. Materials and Structures, 55(1), pp. 1-12, 2022. DOI: https://doi.org/10.1617/s11527-022-01883-2Links ]

How to cite: Zapata-Ferrero, I., Rebollo, O., Rivera, J. and Botasso, G., Fatigue evaluation of asphalts highly modified with recycled tire rubber. DYNA, 90(226), pp. 139-145, April - June, 2023.

I. Zapata-Ferrero, Is BSc. Eng. in Civil Engineer from the National University of Technology, Argentina. PhD student, materials mention in National University of Technology, Argentina. Research in LEMaC, Center for Road Research of the National University of Technology, Argentina. ORCID: 0000-0002-8779-4645

O. Rebollo, Is BSc. Eng. in Civil Engineer from the National University of Technology, Argentina. PhD student, materials mention in National University of Technology, Argentina. Research in LEMaC, Center for Road Research of the National University of Technology, Argentina. ORCID: 0000-0002-0728-8794

J. Rivera, Is Doctor in Engineering, MSc. in Transportation and Logistics, and BSc. Eng. in Civil Engineer from the National University of Technology, Argentina. Mention to the merit of the College of Engineers of the Province of Buenos Aires, Second Prize of the International Award Juan Antonio Fernández del Campo for Innovations in Highway from the Spanish Highway Association, and Award from the Argentine Society of Traffic Engineers for Road Safety. Subdirector of the LEMaC, Center for Road Research of the National University of Technology, Argentina. Exchange researcher at the Universitat Politècnica de Catalunya, and the University of Zaragoza, Spain. Professor of postgrads in Road Engineering from National University of Technology, and University of Buenos Aires, Argentina, Pontifical Catholic University, and Private Technical University of Loja, Ecuador, University of Piura, Peru, University of the Republic, Uruguay, and Catholic University of Our Lady of the Assumption, Paraguay. ORCID: 0000-0001-7391-4469

Gerardo Botasso, Civil Engineer from the National University of Technology, Argentina. PhD student, materials mention in National University of Technology, Argentina. Research in LEMaC, Center for Road Research of the National University of Technology (Argentina). Civil Engineer, Master in Environmental Engineering and PhD in Materials from the National Technological University. Director of the LEMaC Center Highway Research Center - National Technological University La Plata Regional Faculty CIC PBA. Secretary of Science Technology and Postgraduate of the National Technological University, La Plata Regional Faculty. Teacher in master's and doctoral programs of 10 Master's and Doctoral programs in Latin America. Researcher Category A of the National Technological University. Consultant in road infrastructure program in Argentina, specialized in technology of road materials. Author of more than 80 articles published internationally in specialized magazines and 2 books on asphalt mix technology with the incorporation of recycled rubber. Editor of “Ingenio Tecnológico”. ORCID:0000-0002-8859-7256

Received: December 21, 2022; Revised: May 25, 2023; Accepted: May 30, 2023

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License