Abstract
In this study, with the objective to develop a reliability-based seismic design tool, ductility and dissipated hysteretic energy uniform annual failure rate (UAFR) spectra are obtained and compared using the spectral acceleration at first mode of vibration of the structure Sa(T1) and the well-known spectral shape-based intensity measure INp. Notice that this is the first time in the literature that UAFR spectra are obtained for the advanced spectral shape intensity measure INp. For this aim, 110 simulated ground motions recorded from the soft soil of Mexico City were selected due to their large energy amount demanded to the structures; moreover, four elastoplastic hysteretic behavior models are considered for the dynamic analyses with post-yielding stiffness of 0, 3, 5, and 10%. It is observed that the use of elasto-perfectly plastic models provided similar UAFR spectra in comparison with hysteretic models with different post-yielding stiffness. This conclusion is valid for the two selected intensity measures. In addition, the lateral resistance required to achieve similar structural reliability levels is larger when the INp intensity measure is used, especially for buildings with vibration periods equal or larger than the soil period, in such a way that the traditional use of Sa(T1) could provide structures with less structural reliability levels.
1. Introduction
One of the main features to design structures subjected to earthquakes is the use of design or response spectra. The spectra provided by the seismic regulations are determined in most of the cases through single-degree-of-freedom (SDOF) systems with elastic behavior subjected to a set of seismic records, which are reduced to account for nonlinear behavior via ductility reduction factors. The current version of Mexico City Building Code (MCBC) and most of the codes around the world use the spectral acceleration at first mode of vibration of the structure Sa(T1) to estimate the design lateral resistance and stiffness demand of buildings under earthquakes and to control the maximum lateral displacement demand. However, it is known that Sa(T1) presents some limitations when it is used as intensity measure due to its lack of efficiency to predict the nonlinear structural response [1–4]. For example, Sa(T1) does not consider the effect of the elongation of the vibration period when nonlinear behavior occurs. For this reason, Bojórquez and Iervolino [5] proposed a spectral shape parameter named Np and the INp intensity measure toward more efficient parameters. Several researchers have demonstrated the great potential of the INp intensity measure [1, 6–11]. In general, the studies have concluded that the most efficient scalar or vector-valued intensity measures are those based on the Np spectral shape parameter as in the case of INp. Thus, it is important to provide new design or response spectra based on more efficient parameters such as INp. In addition, another two important issues should be accounted for seismic-resistant design of structures. The first one is the inclusion of plastic deformation demands through energy concepts. The use of energy for this purpose was initially discussed by Housner [12]. This concept has been used by various researchers to propose energy-based methodologies aimed to provide the structure with a capacity higher or equal to that demanded by earthquakes [13–20]. In particular, the hysteretic or normalized hysteretic energy has been selected for energy-based design because it is related to the structural damage [21, 22]. On the other hand, the current seismic designs do not take into account explicitly the specific reliability levels in the designed structures. Most of the regulations around the world are mainly based on studies of SDOF systems with elasto-perfectly plastic hysteretic behavior for seismic design of structures, and they do not guarantee the same failure rate on structures as observed in previous studies [23–25]. The seismic design spectra of structures proposed by current codes are not usually associated with specific reliability levels or annual failure rates [26, 27].
This work is motivated by the need to consider the cumulative damage, the structural reliability level, and efficient ground motion intensity measures for earthquake-resistant design of buildings based on the use of uniform annual failure rate spectra. For this reason, in this study ductility and normalized hysteretic energy uniform annual failure rate spectra based on INp and Sa(T1) are obtained and compared. To achieve the objectives of this study, nonlinear systems with different post-yielding stiffness, ductility, and normalized hysteretic energy capacities and several structural reliability levels are studied. To compute the uniform annual failure rate spectra, the structures are subjected to 110 simulated narrowband ground motions of the soft soil of the Valley of Mexico, and the results of UAFR spectra are presented for the case of ductility and normalized hysteretic energy. Notice that the study can be expanded using other types of ground motion records, such as those recommended in several works or ground motion selection procedures [28–33]. This study is limited to these seismic records because in order to obtain just one spectrum for a specific intensity measure, ductility value, or hysteretic energy and hysteretic curve, the authors required more than 5 million of nonlinear dynamic analyses, which is a high computational work. In addition, the implications of using simplified models such as the well-known elasto-perfectly plastic behavior to obtain the lateral resistance requirements when the two intensities Sa(T1) vs. INp are used are estimated. Finally, it is important to say that this is the first time in the literature that UAFR spectra are obtained for the advanced spectral shape intensity measure INp toward future energy-based seismic design taken into account the structural reliability.
2. Methodology
2.1. Estimation of UAFR Spectra
The ductility and hysteretic energy UAFR spectra were estimated via the previous works developed by Esteva [34] and Cornell [35] based on the total probability theorem, where the annual failure rate can be defined as next:where corresponds to the absolute value of the derivative of the site seismic hazard curve; Q expressed the relationship between demand and capacity; and is the conditional probability of failure given a seismic intensity . Notice that the structural failure occurs when the capacity is smaller than demand, and in other words:
The procedure to compute the UAFR spectra is as follows:(1)Selection of the structural systems For this study, four elastoplastic-based hysteretic behaviors are considered for the nonlinear models with different post-yielding stiffness ratios. The first single-degree-of-freedom system selected corresponds to the elasto-perfectly plastic hysteretic behavior, and this model can exhibit a linear elastic behavior until a value of strength called yielding point and then behave in a plastic mode until a maximum displacement is obtained. In addition, three elastoplastic models with different post-yielding stiffness ratios corresponding to 3, 5, and 10% have been selected (see Figure 1).(2)Selection of the simulated earthquake ground motion records In the seismic design of structures, it is necessary to have seismic records of different magnitudes, and in a sufficient quantity so that it is possible to characterize, in a precise way, the seismic hazard; for this, it is required to identify and take into account all the seismic sources with potential influence in the study site. Unfortunately, the short time of observation and the lack of instrumentation in many regions make it difficult to obtain such records, particularly for large earthquakes. This has motivated the development of several methodologies of seismic simulation to obtain synthetic records of great magnitude earthquakes [36, 37]. In particular, the methods of seismic simulation that use a stochastic summation approach with Green’s empirical functions (FGE) are widely accepted in structural engineering, due to their easy application to generate synthetic records. Such methodologies have their origin in 1978 thanks to Hartzell [38]. In this study, it was selected to use the methodology developed by Niño et al. [39]. They propose a source model defined by two corner frequencies and a summation scheme divided into two stages [40, 41], to obtain more accurate spectra to the conditions that affect Mexico City, as well as to improve the definition of the changes in the amplitudes, according to the size of the source and the total duration of the seismic event [42]. With this method, 110 seismic records were simulated for different moment magnitudes Mw from 7.2 to 8.2, with a magnitude increase . For the purpose of this study, the record of the seismic event of April 29, 1989, obtained in the SCT station was used as seed for the simulations; in addition, a stress drop bars was considered. The corresponding parameters of moment of magnitude, Mw, and the two corner frequencies for each stage for the events employed as seed and each of the simulated magnitudes are presented in Table 1. The average response spectra of the simulated ground motion records are illustrated in Figure 2 for each selected magnitude. Notice the large amplifications of the pseudo-acceleration in a well-defined region of the spectra (1.5–2.5 s). Figure 3 shows the elastic pseudo-acceleration response spectra for all the simulated ground motion records and 5% of damping ratio; moreover, the mean, 25th percentile, and 75th percentile earthquake response spectra also are illustrated.(3)Selection of the ground motion intensity measures As it was previously indicated, the most used intensity measure by most of the seismic design codes is the spectral acceleration at first mode of vibration (especially for low-rise structures) [2, 5], because it is the perfect predictor of the response of elastic single-degree-of-freedom systems and multi-degree-of-freedom systems dominated by the first mode of vibration. Nevertheless, Sa(T1) does not provide information about spectral shape in other regions of the spectrum, which is essential for nonlinear behavior of structures or for which their behavior is dominated by higher modes (modes with periods below than T1). On the other hand, in recent studies some researchers have proposed parameters to define appropriate intensity measures, which are related to the spectral shape due to its relationship with the structural response. Recently, following the approach to develop spectral shape parameter as intensity measure, Bojórquez and Iervolino [5] proposed the well-known spectral shape parameter Np, which takes into account the nonlinear behavior of the structures by including several points of a response spectrum. The Np spectral shape parameter is defined as the ratio of the geometrical mean between the periods T1 and TN normalized with respect to Sa(T1) as it is illustrated in the following equation: The information given by the Np equation is that if we have one or n records with a mean Np value close to one, we can expect an average spectrum near flat in the period range between T1 and TN. For a mean Np lower than one, it is expected an average spectrum with a negative slope, and finally, Np values larger than one are representative of spectra with positive slope increasing the spectral ordinates. To incorporate the influence of nonlinear behavior in the prediction of structural response, Bojórquez and Iervolino proposed a new scalar ground motion intensity measure based on Sa(T1) and Np, which is described by the following equation: In (3), INp is a scalar ground motion intensity measure, Sa(T1) is the commonly known intensity measure, and Np is the spectral shape parameter, and the α value must be calibrated according to the structure and the selected seismic demand parameter. Different researchers have demonstrated the great potential of the INp intensity measure [1, 6, 7, 43]. For this reason, in this study, the two selected ground motion intensity measures are Sa(T1) and INp to propose new UAFR spectra based on ductility and dissipated hysteretic energy.(4)Selection of the performance parameters (ductility and hysteretic energy) In this work, the ductility and normalized hysteretic energy were selected as performance parameters. Ductility The ductility is defined as the ability of a system to be deformed without significant loss of resistance. The term ductility factor is defined as a measure of the amount of energy dissipation capacity as follows: where is the maximum absolute displacement achieved by the system and is the yielding displacement. For this study, ductility capacity values of 2, 3, and 4 have been selected according to the recommendations of the Mexican City Building Code. Normalized hysteretic energy Energy-based methodologies focus on providing structures with energy dissipation capacities larger than or equal to their expected energy demands [13, 14]. The energy most related to the structural damage is the hysteretic energy . The dissipated hysteretic energy can be interpreted physically by considering that it is equal to the area enclosed by each of the hysteresis loops that the structure develops during a seismic excitation. Although hysteretic energy provides an approximate idea of the accumulated plastic deformation demands, this response parameter by itself does not provide sufficient information to evaluate the structural behavior, so it is convenient to normalize it as follows: where and are the strength and displacement at first yield, respectively. is a parameter that correlates best with the structural damage [44, 45]. In this work, the normalized hysteretic energy is considered as a parameter to control the accumulated damage. Thus, UAFR spectra are computed.(5)Estimation of the structural response of the system by incremental dynamic analysis (IDAs) [46] for a selected period T1 and yield force coefficient Cy (defined as the yielding force divided by the total weight of the system) subjected to ground motion records scaled in terms of Sa(T1) or INp. The seismic responses obtained in this study are the ductility and normalized dissipated hysteretic energy as it was previously discussed.(6)Proposing specific values of ductility capacity or normalized dissipated hysteretic energy. In this work, the ductility capacity values correspond to those indicated by the Mexico City Building Code. On the other hand, different values for the normalized hysteretic energy capacity are suggested in Bojórquez et al. [18]. For this reason, in this study the values of 3, 6, 9, and 12 of the normalized hysteretic energy capacity have been selected, which are representative of reinforced concrete and steel buildings according to Terán and Jirsa [21] and Bojórquez et al. [18, 47].(7)Assessment of the failure probability by dividing the number of ground motion records in which the ductility of energy demanded is larger than the capacity (step 6) between the total numbers of records used for a specific intensity level.(8)Selection of the seismic hazard curve based on the seismic intensity measures Sa(T1) and INp. For the present word, the seismic hazard curves computed by Rodríguez et al. have been used [48, 49].(9)Numerical assessment of the annual failure rate using (1). The steps for other structural periods and Cy to obtain the structural annual failure rate curves are repeated.(10)The UAFR spectra for a specific value are calculated for each of the parameters here studied.
Finally, Figure 4 shows a flow chart of the procedure to compute the UAFR spectra.
3. Numerical Results
The numerical results obtained in this work are presented in this section. It is important to say that all the UAFR spectra illustrated correspond to the spectra obtained via the 110 simulated narrowband earthquake ground motions. Notice that the selection of the seismic hazard curves based on the intensity measure Sa(T1) and INp is a very important issue to compute the UAFR spectra. The seismic hazard curves are used frequently to represent the seismic hazard of a specific site. They indicate the annual rate of exceeding a variety of intensity levels of a ground motion parameter at a site of interest (i.e., Sa(T1) and INp). The procedure to compute a ground motion hazard curve is based on the total probability theorem [34, 35, 50–52]. As it was indicated previously, the seismic hazard curves for Sa(T1) and INp developed by Rodríguez et al. 2021 [40, 41] have been selected to compute the uniform annual failure rate spectra. For example, Figure 5 illustrates the seismic hazard curves for the SCT site corresponding to soft soil of Mexico City for periods ranging from 0.2 to 2 seconds and both selected intensity measures.
(a)
(b)
3.1. Influence of Ductility in the UAFR Spectra in terms of Sa(T1)
The ductility UAFR spectra when Sa(T1) is used as intensity measure and for all the hysteretic models under consideration are illustrated in Figure 6. Furthermore, the UAFR spectra are compared for different ductility capacity levels. While Figure 6(a) shows the results for the elasto-perfectly plastic hysteretic model, Figures 6(b)–6(d) correspond to elastoplastic with 3%, 5%, and 10% of post-yielding stiffness, respectively.
(a)
(b)
(c)
(d)
All the spectra are associated with a UAFR . It is observed that as the ductility increases, the required lateral resistance or seismic coefficient (Cy) decreases, especially for structures close to the soil period (Ts = 2s). The largest differences are observed in the intervals of the structural periods ranging from 1.3 to 2 seconds. For example, in Figure 6(a), when T = 0.8 seconds (structures with periods smaller than the soil period, which is equal to 2s) and , the Cy is equal to 0.25, and for , the Cy is equal to 0.21, indicating a reduction of 16% when the ductility increases from 2 to 4. On the other hand, if the structure is close to the soil period, for example, for a system with T = 1.6s and , the Cy is equal to 0.41 and for the Cy is equal to 0.24. Thus, a reduction of 41.4% is obtained if the ductility capacity increases. Therefore, the structures with vibration period close to the soil period require special attention; in particular, this type of building is more sensitive to the ductility capacity. Similar results are observed for the elastoplastic models with different post-yielding stiffness analyzed in this study as observed in Figures 6(b)–6(d).
3.2. Influence of Post-Yielding Stiffness in the UAFR Spectra in terms of Sa Intensity Measure
The influence of post-yielding stiffness is obtained by comparison of the elasto-perfectly plastic and elastoplastic model with 3% (BL03), 5% (BL05), and 15% (BL15) of post-yielding stiffness. Figures 7 and 8 compare the UAFR spectra with and for two different ductility capacities ( and ). The selected ductility values are representative of structures with low and high levels of ductility capacity. It is observed that in both cases, the spectra are similar for all the post-yielding stiffness; therefore, it can be concluded that the post-yielding stiffness has low influence to compute the UAFR spectra, except in the case of structures with period near or about the soil period. For these cases, it is observed that as the post-yielding stiffness increases the seismic coefficient (Cy) decreases in low ratio for structures with moderately low ductility capacity (see Figures 7(a) and 8(a)). On the other hand, it is shown that for structures with high ductility (see Figures 7(b) and 8(b)), the effect of post-yielding stiffness in all the analyzed models is more evident. However, the elasto-perfectly plastic hysteretic behavior model provides reasonable results to estimate the required lateral strength in comparison with the elastoplastic model with different post-yielding stiffness.
(a)
(b)
(a)
(b)
3.3. Effect of the Annual Failure Rate () in the Spectra in terms of Sa Intensity Measure
In this section, the influence of the annual failure rate in the UAFR spectra is studied. For this aim, three failure rates are selected: , respectively. Figure 9 illustrates for a ductility value equal to 2 and the elasto-perfectly plastic hysteretic model the effect to compute spectra for the selected annual failure rates. It is observed that as the annual failure rate decreases, the seismic coefficient tends to increase, especially for structural periods smaller or equal to the soil period. For this reason, if the designed buildings require larger structural reliability levels, it is required to increase the lateral resistance. For example, when T = 1.0 seconds with , Cy is equal to 0.20; on the other hand, if , Cy = 0.33. Therefore, Cy increases about 65%, while for a T = 3.0 seconds with , Cy is equal to 0.09 and for the Cy = 0.13, and for this case, the percentage is 44.4% larger. In conclusion, it is observed that for structural periods larger or very larger than the soil period, the influence of the uniform annual failure rate selected to compute Cy tends to be despicable.
3.4. Influence of the Hysteretic Energy in the UAFR Spectra in terms of the Sa Intensity Measure
The normalized hysteretic energy UAFR spectra when Sa is used as intensity measure for all the models with different post-yielding stiffness under consideration are illustrated in Figure 10. The UAFR spectra are compared for different levels of energy capacity for the nonlinear models. Figure 10(a) shows the results for the elasto-perfectly plastic model and Figures 10(b)–10(d) for 3%, 5%, and 10% of post-yielding stiffness. Notice that all the spectra are associated with a UAFR . It is observed that as the energy increases, the required lateral resistance decreases, particularly in the region of periods close to the soil period. The same trend is observed for the different post-yielding stiffness selected in this study.
(a)
(b)
(c)
(d)
The influence of post-yielding stiffness on the behavior models in the UAFR spectra is compared in Figures 11 and 12, with and . It is very interesting to observe that the figures suggest that the UAFR spectra are similar for the different post-yielding stiffness. This conclusion is also valid if the energy capacity selected or annual failure rate is different.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
3.5. Effect of Failure Rates in Hysteretic Energy UAFR Spectra in terms of Sa Intensity Measure
In this section, the influence of the annual failure rate in the UAFR spectra in terms of the dissipated hysteretic energy capacity is studied. For this aim, three failure rates are selected: respectively; see Figure 13. It is observed that as the failure rates decrease, the required lateral resistance increases. Similar results were found for all the periods under consideration. It is important to say that similar results were obtained when comparing the UAFR spectra in terms of INp as intensity measure with respect to the previously mentioned parameters. The following section compares the results obtained for Sa and INp intensity measures.
4. Comparison of Ductility UAFR Spectra: Sa vs. INp
In this part, the results of ductility UAFR spectra for the previously mentioned intensity measures are compared. In Figures 14 and 15, the elasto-perfectly plastic behavior for different ductility values with annual failure rates equal to 0.004 and 0.008 is used to compare the UAFR spectra. Figure 14(a) shows the UAFR spectra for a ductility value equal to 2. It is observed that the lateral resistance required for both intensity parameters is similar for periods less than 1.6 seconds, while for periods between the ranges of 1.6 and 2.5 seconds a significant difference is observed when using both intensity measures selected to compute the lateral resistance required. For example, from Figure 14(a) and a period equal to 1.3 seconds, the values of Cy = 0.31 for INp and Cy = 0.33 for Sa have a difference of 6%, and for a period equal to 2.0 seconds, a Cy = 0.23 is obtained for Sa and Cy = 0.30 for INp, and this difference represents about 30%. On the other hand, Figure 14(b) shows the UAFR spectra for ductility equal to 4, where both intensity measures are compared. The results suggest that the value of Cy is very similar for all the periods of the systems under consideration.
(a)
(b)
(a)
(b)
The effect of the influence of the ductility value is compared in Figures 14(a) and 14(b). It is found that as the ductility increases, the required lateral resistance decreases considerably for all periods of the systems analyzed in this study, varying from 30% to 60% as the ductility value increases from to , respectively. Similar results are found in Figures 15(a) and 15(b) for annual failure rate equal to 0.008. In general, the results suggest that for structures with low ductility capacity such as masonry or reinforced concrete buildings, the selection of an appropriate intensity measure is very important in the case of ductility UAFR spectra at least when several ground motion records are used to compute the average spectrum. It is important to say that the authors will study in future works the importance of efficient intensity measures and the number of records used for nonlinear dynamic analysis to compute the UAFR spectra.
Figures 16 and 17 show the UAFR spectra for ductility values 2 and 4, respectively, where the influence of failure rates ( and ) also is compared. It is observed that as the values of the annual failure rate decrease, the required lateral resistance increases up to 30% for structural periods less or equal to the soil period, and the increment is up to 20% higher for periods larger than the soil period. This observation applies for both ground motion intensity measures studied in this work (see Figures 16 and 17).
(a)
(b)
(a)
(b)
5. Comparison of Normalized Hysteretic Energy UAFR Spectra: Sa vs. INp
The normalized hysteretic energy UAFR spectra are compared in Figures 18 and 19 for the intensity measures Sa and INp and the elasto-perfectly plastic behavior model with one annual failure rate and different levels of energy capacities. Figures 18(a) and 18(b) show the UAFR spectra associated with annual failure rate equal to 0.004. A similar effect is observed for both intensity measures for the structures with periods less than 1.5 seconds. However, for structural periods between 1.6 and 2.3 seconds (particularly in the region of periods close to the soil period), an important effect is found when using Sa or INp, where the seismic coefficient Cy increases more in the case of INp compared with Sa. This conclusion is valid for the annual failure rate equal to 0.004 and 0.008 (see Figures 18 and 19).
(a)
(b)
(a)
(b)
The following figures (Figures 20 and 21) compare the influence of the annual failure rate at different normalized hysteretic energy values. Figures 20 and 21 show the UAFR spectra for EN equal to 6 and 9, respectively, where the influence of failure rates ( and ) is compared. The figures illustrate that as the values of decrease, the required lateral resistance increases up to 20% for structural periods close to the soil period. Finally, because INp is a more efficient intensity measure in comparison with the most used intensity measure Sa(T1), the results indicate that in the case of cumulative demands such as hysteretic energy, the lateral resistance required is larger when an appropriate intensity measure such as INp is selected. In fact, for example for a structural period equal to 2 (the soil period), the seismic coefficient could be up to 34% larger when using INp instead of Sa. Thus, the use of the intensity measure Sa(T1) could provide a lateral resistance that produces unsafety structural designs accounting for cumulative plastic deformation demands. In other words, for energy-based design or energy-based spectra it is concluded that INp should be used as intensity measure. Finally, Figure 22 compares the ratio of the seismic coefficient obtained for the normalized hysteretic energy UAFR spectra in terms of INp divided by that obtained via the spectral acceleration in the case of the elasto-perfectly plastic model with annual failure rate equal to 0.004 and different values of EN. This figure confirms that the seismic coefficient could be as larger as 45% when using INp in comparison with Sa(T1).
(a)
(b)
(a)
(b)
6. Conclusions
In this study, ductility and dissipated hysteretic energy UAFR spectra have been computed and compared using spectral acceleration at first mode of vibration and the well-known INp intensity measures. The results obtained for ductility and hysteretic energy UAFR spectra are as follows.
6.1. Ductility
(i)The ductility UAFR spectra when Sa and INp are used as intensity measures were computed. It was observed that as the ductility increases, the required lateral resistance decreases, where the largest differences occurred for structures close to the soil period. Therefore, the structures with vibration period close to the soil period require special attention; in particular, this type of building is more sensitive to the ductility capacity.(ii)The influence of post-yielding stiffness also is analyzed. It was observed that the ductility spectra are very similar for the selected hysteretic models with different post-yielding stiffness, in such a way that the elasto-perfectly plastic behavior model can give reasonable results for estimating the Cy with bilinear behavior and different percentages of post-yielding stiffness.(iii)In the case of the influence of the annual failure rate in the UAFR spectra, it was found that as the annual failure rates are reduced, the Cy increases in a reasonable rate. For this reason, if the designed buildings require larger structural reliability levels, it is necessary to increase the required lateral resistance.(iv)Finally, for the case of ductility UAFR spectra, the results suggest that for structures with low ductility capacity such as masonry or reinforced concrete buildings, the selection of an appropriate intensity measure is very important at least when several ground motion records are used to compute the average spectrum.6.2. Normalized Hysteretic Energy
(i)For the case of normalized hysteretic energy UAFR spectra and both intensity measures used, it is concluded that as the EN increases, the required lateral resistance decreases, particularly in the region of periods close or larger to the soil period. Notice that the same trend is observed for all the selected post-yielding stiffness under consideration.(ii)As in the case of ductility, no influence of the post-yielding stiffness was observed to compute the average normalized hysteretic energy UAFR spectra.(iii)As it was expected, as the annual failure rate decreases (structures with large reliability levels) the required lateral resistance increases.(iv)Finally, the use of the intensity measure Sa(T1) could provide a lateral resistance that produces unsafety structural designs accounting for cumulative plastic deformation demands. In other words, for energy-based design or energy-based spectra it is concluded that INp should be used as intensity measure.Data Availability
The data used to support the findings of this study are available from the corresponding authors upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this study.
Acknowledgments
The financial support given by the Universidad Autónoma de Sinaloa under Grant PROFAPI 2022 is appreciated. The authors express their gratitude to the Consejo Nacional de Ciencia y Tecnología (CONACYT) in Mexico for funding the research reported in this study under Grant Ciencia Básica 287103 and for the scholarships given to the PhD students. The support of the UNAM-DGAPA-PAPIIT under Project No. IN100320 is appreciated.