Abstract

As cities grow, so does the need for housing and transportation infrastructure, necessitating deep excavations. Near the deep excavation, especially in the urban centers, there are existing structures whose loads are often different, resulting in asymmetric loading of the excavation system and different behavior on each side of the excavation. In this study the strut forces, lateral deformations, bending moments on each wall, as well as the settlements and heaves of the soil behind both walls in the asymmetrically loaded strut-supported excavation system are determined using the finite element method (FEM). The analyzes are performed using Plaxis 2D v20, a computer software based on the FEM. In addition, regression analyses are performed on the results, and correlations that predict the relevant results depending on the surcharge loads are presented.

1. Introduction

With the increasing population in recent years, high-rise buildings have been needed, especially in the city centers. With the increasing need for public transportation, subways have also become more prevalent. Deep excavations are a prerequisite for structures such as high-rise buildings and subway stations. Deep excavations must be supported so that they do not damage adjacent structures, ensure their safety and reduce the deformations. Many studies in the literature examine the deformations of the walls and surface soil and the loads on the struts and anchors in the deep excavations. Bransby and Milligan [1] established a laboratory model to determine the deformations in the soil near the cantilever sheet pile wall on dry sandy soil and presented an analytical method showing a relationship between the deformations at the soil surface and the displacements at the wall. Ou et al. [2] analyzed the case studies and showed the location and relationships between the maximum horizontal wall deformation and maximum soil surface settlement. They also proposed an empirical formula for determining the profile of soil surface settlement. In addition, Clough and O’Rouke [3] evaluated wall movements on site-built walls due to excavation, shoring, and auxiliary work separately from field observations and introduced a new understanding of modeling excavation problems using the finite element method (FEM). Cheng et al. [4] developed a simplified analytical method using the theory of elasticity to determine the ground lateral displacement resulting from foundation pit excavation and validate the outcomes of the suggested analytical solution by contrasting them with well-documented field data. Additionally, they perform a parametric analysis to identify the variables that most significantly influence ground lateral displacement. Wong and Goh [5] studied the effect of cracks in the wall on the behavior of the wall during top–down excavation and they found that the crack in the wall significantly reduced the bending moment but slightly increased the lateral displacement.

Subsequently, some researchers conducted finite element analyses to model the experimental models [6, 7] and case studies and field observations [5, 810] to calculate the soil surface settlement, horizontal displacements of the wall, and soil stresses. Ghobrial et al. [11] determined the lateral earth pressures by modeling the excavation with the FEM and found that the results were compatible with the apparent pressure diagrams developed by Peck [12] and Tschebotarioff [13]. Hsu et al. [14] calculated soil stresses occurring during the excavation in gravel soils using analytical methods [12, 13, 15] and compared them with data obtained from the case studies. With the development of computer technology, many researchers have carried out parametric studies to investigate the vertical deformations that will occur in sandy soils and the change in wall deformations and soil pressure using the 2D FEM and finite difference method. In these studies, parameters such as sand density, excavation depth and width, wall length and stiffness, and pier stiffness and spacing were considered by performing the analysis using computer-based software such as Plaxis 2D v20, Midas Gtx, Abaqus [1626].

Hsiung et al. [27] analyzed wall movements, and Chheng and Likitlersuang [28] additionally analyzed soil surface settlements using the 3D FEM, compared them with field measurements, and found that the results were reasonable and reliable. Goh et al. [29] performed both 2D and 3D finite element analyses and showed that the 3D analyses yielded lower wall deformations and proposed an equation for estimating the maximum wall deformation. Zhang et al. [30] performed both 2D and 3D finite element analyses for dense sand and gravel soils and a medium-rigid wall, and claimed that current studies in the literature did not address these soil types much. The results obtained were compatible with the field measurements, and a correlation was proposed to determine the strut force based on the parametric study results. Goh and Wong [31] performed a 3D finite element analysis to study the load change in the adjacent struts when one or two struts fail and found that transfers its load to the adjacent struts, and the bending moment on the wall does not change significantly. In addition, Pong et al. [32] proposed a procedure to optimize the strut forces determined by the 2D analyses according to the results of the 3D analyses in the case of failure of a single strut because the 2D analyses require more time than the 3D analyses. Some researchers study soft soils using the numerical methods and computer-based software.

Although there are buildings and roadways close to the deep excavations, particularly in the city center and it is more likely that the surcharge loads are asymmetrical, the excavations are typically modeled using a half-model assuming the loads are symmetrical. However, the strut forces, the walls’ behavior and the soil surface’s settlements changes because each wall responds to asymmetic loading in different way hazardous and unprofiteable results may be caused. By carefully analyzing and designing the walls, it is feasible to get rid of the hazardous conditions and find a cost-effective solution. Wu [33] performed a simplified and limited analysis using Winkler’s model for the case of asymmetric surcharge load. In addition, Xu et al. [34] and Guo et al. [35] performed a finite element analysis and finite difference analysis, respectively, and found that the presence of surcharge loads on both sides of the excavation has different magnitudes, resulting in different lateral deformations and bending moments for both walls and that the asymmetry negatively affects the support system of the excavation. Zhang et al. [36] used a numerical simulation of stage excavation to study the excavation of a subway station next to a river and a main road. The asymmetrical pressure of the excavation induced by the distance between the external load and excavation is studied using parameter analysis. The findings show that asymmetrical pressure on the excavation is significant and cannot be ignored if there is a river channel and an external load within the excavation depth range. The influence of asymmetrical pressure can be reduced, and the stability of adjacent riverside layers may be enhanced by increasing the stiffness of the diaphragm wall. There are very few studies on this topic in the literature so far, so further studies on this topic will contribute to the literature.

Diaphragm walls are strong concrete structures that are framed panel by panel. Diaphragms are typically made of concrete slabs for concrete structures or metal or composite metal decks for steel structures. Diaphragm walls are mostly employed near existing structures when the earth and depth would make piled walls difficult. Diaphragm walls are simple and inexpensive to construct and available worldwide. Tunnels, deep basements, underpasses, underground car parks, and railway stations are some of the most popular applications. It has three main types flexible, rigid, and semirigid diaphragms. Flexible diaphragms withstand lateral forces based on the tributary area, regardless of the flexibility of the parts to which they transfer force. Rigid diaphragms transfer the load to frames or shear walls depending on their flexibility and position in the structure. The flexibility of a diaphragm influences the distribution of lateral forces to the vertical components of a structure’s lateral force-resisting parts.

The present study aims to investigate the change in strut loads, deformations, and bending moments on the diaphragm walls, ground surface settlements, and heave in the strut-supported deep excavations under asymmetric surcharge loads. In this context, strut-supported deep excavation systems have diaphragm walls in loose (LS) and dense (DS) sandy soils are modeled for various loading scenarios using PLAXIS 2D v20 [37], a finite element-based computer program. Parametric analyses are first performed to determine the effects of model boundaries (X, Y), mesh size, wall-soil interface strength reduction (Rinter), the width of the pit (B), embedment depth (Df), and thickness (d) of the walls to determine ascertain the main model for the study. Then, a comprehensive parametric study covering the symmetric and asymmetric surcharge load scenarios for LS and DS is carried out and the results are presented in tables and illustrated in figures. In addition, linear and polynomial second-order relationships are obtained using the SPSS Statistics v17 computer program to estimate the maximum values of the strut forces (Tst), the deformations (UxL and UxR) and bending moments (ML and MR) on the diaphragm walls, the ground surface settlements (UyL and UyR) for both the left and right sides of the excavation and the heaves (Hv) in the excavation as a function of the surcharge situations on both sides of the excavation.

2. Materials and Methods

2.1. Finite Element Method

It is difficult to solve the geotechnical correlations because there are so many factors involved, including strain, stress, shear pressures, deflections, and bending moments, makes it complicated. Because of this, numerical modeling is a crucial tool for streamlining and improving the accuracy of complex calculations. The FEM is a technique for finding solutions to engineering problems by assembling the variables into vectors and matrices in computer programs and solving the system numerically. These challenging solution processes can be easily solved using a variety of technologies, including the Plaxis program used in the study. The current study employs a two-dimensional FE model considers plane strain model, and the excavation models are modeled as a whole for symmetric and a symmetric loading scenarios. The effects of mesh size on the behavior of the excavation support system are minimized by using a “fine mesh.” The “fine” option from the program’s settings was used to generate the mesh networks depicted in Figure 1. There are 1,942 soil elements and 16,055 nodes in the linked mesh network. There are two diaphragm walls, one on the left and one on the right side of the excavation, simulated with a five-noded elastic plate and connected to the surrounding soil by an interface element. The parameters of the walls and struts are taken from Xu et al. [34]. In addition, the wall’s area and moments of inertia are changed so that the elastic modulus remains constant to study the behavior of the walls at different thicknesses. The boundary conditions allow vertical movement of the sides and the vertical and horizontal movement of the upper boundary, while no movement is permitted at the lower boundary. The horizontal boundaries are 350 m from the center of the excavation on both sides, and the lower boundary is 350 m below the ground surface to avoid the negative effects of the model boundaries. The struts are modeled with a “node-to-node anchor” at 3-m horizontal and vertical intervals. The parameters of the excavation shoring system are given in Table 1.

Excavation and wall’s embedment depths are both assumed to be 16 m, the support system is designed with five struts, and groundwater is not considered. In the analysis, the loading type “staged construction” is chosen. The excavation is carried out in eight stages, and the applied construction phases are shown in Table 2.

Each excavation is made to a depth of 1 m below the next bracing level, and the corresponding strut is activated in the next stage. The schematic model of the five-stage bracing system, diaphragm walls and surcharge loads is shown in Figure 2.

2.2. Constitutive Model

In the present work, LS and DS are used as soil materials and the hardening soil model is adopted, which takes into account the friction-induced swelling behavior in soil and allows modeling under the triaxial volume deformation conditions, where the Mohr–Coulomb failure criterion defines failure. The hardening soil parameters used in the study are taken from the studies of Brinkgreve et al. [38] and Xuan [39] and are listed in Table 3.

3. Results and Discussion

3.1. The Validation of the Model

Problems of the three research provided in the literature are modeled using Plaxis 2D v20 and the results are compared to the results presented in the relevant paper for validation. The strut forces from Zhang et al. [30] are compared in Figure 3, the ground surface settlements from Mohamed et al. [24] are compared in Figure 4, and the wall deflections for flexible and stiff walls from Goh et al. [40] are compared to the results of current models in Figure 5.

The results are found to be in fair agreement with slight fluctuation. Many parameters considered when developing the model have an impact on the results and some parameters are not explicitly stated in previous studies that are compared. The modest discrepancies in the results are assumed to be related to this uncertainty.

3.2. The Effects of the Model Parameters

The behavior of the excavation support system changes when the models are built with very fine or very coarse mesh size as well as width or narrow. The model width and mesh size cause time and labor losses and unsafety problem when they are not chosen appropriately so it is necessary and beneficial to determine the optimal model size and mesh size. Parametric studies are carried out to determine the width, depth and mesh size of the model to be used in the analyses where the effects of asymmetric surcharge load are studied, which is the main purpose of the paper. At this stage, it is assumed that a surcharge load of 60 kPa exists on both sides of the excavation, which is the largest surcharge load of this study. Here, UyL, UxL, and ML represent the results for the wall on the left side (WL) and the soil behind WL, while UyR, UxL, and MR represent the results for the wall on the right side (WR) and the soil behind WR.

3.3. The Effect of the Boundary Conditions and Mesh Size

Table 4 shows the analysis results performed at different model dimensions for LS and DS considering that the mesh size is fine, Df = 16 m, d = 40 cm, Rinter = 0.6, and B = 45 m. Model dimensions are changed to approximately 20 times the excavation width, keeping the ratio of horizontal width (X) to vertical depth (Y) constant. The results include the maximum values for Tst, UyL, UxL, ML, and Hv.

While the values of Tst, UxL, and UyL hardly change with the change in model dimensions, the values of UyL decrease with the increase in model dimensions. Because the model is supported horizontally on both sides against deformation, this condition causes bulging within the impact limits and affects vertical deformation. This change in UyL loses its effect after Y = 500 m and X = 250 m for DS and Y = 700 m, and X = 350 m for LS. Therefore, the model dimensions are considered as Y = 700 m and X = 350 m in the rest of the study.

Table 5 shows the analysis results performed at different mesh sizes, from very coarse to very fine for LS and DS considering that model dimensions are 700 m horizontally and 350 m vertically, Rinter = 0.6, Df = 16 m, d = 40 cm, and B = 45 m. The results include the maximum values for Tst, UyL, UxL, ML, and Hv.

The Tst and Hv values almost do not change with the change in mesh size. At the same time, UyL decreases slightly, and ML increases significantly from very coarse to medium and gradually stabilizes after the medium mesh size, especially in LS.

3.4. The Effect of the Thickness and Embedment Depth of the Wall

The effects of wall thickness on the deformation and bending moment of the WR as well as the deformations behind the WR and at the bottom of the excavation considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, Rinter = 0.6, Df = 16 m, and B = 45 m is shown in Figure 6.

As the wall thickness decreases, the maximum UxL and UyL increase while the maximum ML decreases, and Hv is unaffected by the wall thickness. When d decrease from 120 to 40 cm, UxL and UyL increase from 17 to 26 cm and −1 to −8 cm, respectively, while ML decreases from 1,020 to 182 kNm. When existing structures exist, thicker walls should be built, and this thickness should be increased by the sensitivity of the existing structure.

Table 6 shows the analysis results performed at various Df to He ratios for LS and DS considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, Rinter = 0.6, d = 40 cm, and B = 45 m. The results include the maximum Tst, UyL, UxL, and ML values.

ML and UxL decrease significantly until the Df to He ratio is 1 while Tst remains almost unchanged, and no significant change is observed for larger ratios. Therefore, in the present study, the ratio of Df to He is assumed to be 1. As the Df decreases the ML increases and increasing the cross-section of the wall and consequently the cost. Simultaneously, due to the increase in the settlement of the soil behind the wall, the safety and usability of existing structures may be threatened. Therefore, determining the embedment depth considering both economy and safety is necessary for the Df values smaller than the He.

3.5. The Effect of the Coefficient of the Strength Reduction in the Interface

Table 7 shows the analysis results performed for different Rinter values for LS and DS considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, d = 40 cm, Df = 16 m, and B = 45 m. The results include the maximum values for Tst, UyL, UxL, ML, and Hv.

As the Rinter increases, the maximum UxL, UyL, and ML in LS and DS decrease, whereas the maximum Hv for DS is hardly affected. Since the reduction is slight for values larger than 0.6, and the friction angle between the wall and the soil is usually assumed to be two-thirds of the internal friction angle in practice, Rinter is supposed to be 0.6 in the present study. The friction between the wall and the soil must be accurately measured and taken into account in the calculations because it can cause especially the deformation values to increase by up to 30%–40% in UxL and 300% in UyL.

3.6. The Effect of the Width of the Excavation Pit

The impact of the width of the excavation on the deformation and bending moment of the WL and the deformations behind the WL and at the bottom of the excavation considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, d = 40 cm, Df = 16 m, and Rinter = 0.6 is shown in Figure 7.

UxL, UyL, ML, and Hv increase as the B increases and the point of maximum of UxL is moving away from the bottom of the excavation. When B increases from 5 to 45 m, maximum UxL, UyL, and ML increase from 10 to 25 cm and from −3 to −8 cm, and from 160 to 185 kNm, respectively. Because a rise in pit width causes significant increases in the horizontal displacement of the wall and soil settlement, it should be operated in as small apertures as feasible. When necessary, alternatives such as progressive excavation should be proposed while using extreme caution.

3.7. The Effects of the Surcharge Loads

The analyses are performed for the no surcharge case, the symmetric surcharge case, and the asymmetric surcharge case. In the loading cases of the symmetrical situation, there are 60, 45, 30, 20, and 10 kPa surcharge loads on both sides of the excavation. When determining these loads, it is assumed that a single-story building exerts pressure on the soil of around 10–15 kPa. Furthermore, because there will be existing buildings around the excavation, particularly in towns, and because they may have varying levels, analyses were performed for equal and different surcharge load scenarios on both sides of the excavation. The asymmetric loading cases include cases where the surcharge load acting behind the WL, called qL, is kept constant at 60 kPa and the surcharge loads acting behind the wall on the excavation’s right side (qR), are 45, 30, 20, 10, and 0 kPa. The analyses use Plaxis 2D v20, a finite element-based computer software. The Tst, deformations and bending moments on the walls, deformations on the soil surface, and the excavation pit are presented in the tables and figures. Later, using SPSS Statistics v17 statistical software, linear and polynomial relationships are obtained, showing their variation as a function of surcharge loads.

3.8. The Effect of Surcharge Loads on the Strut Forces

Table 8 shows the maximum Tst results from the analysis performed for the different symmetrical and asymmetrical surcharge loading conditions for LS and DS considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, d = 40 cm, Rinter = 0.6, Df = 16 m, and B = 45 m. The maximum Tst increases with increasing surcharge load for both asymmetric and symmetric loading for LS and DS. The maximum Tst is 385, 453, and 506 kN respectively, when both sides are empty, one side is loaded with 60 kPa and the other is empty, and both sides are loaded at 60 kPa for DS.

Tstrut increases as the surcharge loads increase, and the asymmetry of the surcharge loads has no additional negative effects on the Tstrut for LS and DS. The maximum Tst values are almost twice as high for LS as for DS. In addition, the correlation between the maximum Tst and the surcharges can be expressed by Equations (1)–(3) for LS and Equations (4)–(6) for DS. It should be emphasized that these equations will only produce valid results for loose and dense sandy soils with the same or close values of wall and sturt properties employed in the present paper, and for surcharge loads of less than 60 kPa.

The Tst increase with depth but decrease especially at the last strut, and the difference between the last two struts decreases as the distance between the previous strut and the bottom of the excavation changes from 3 to 4 m by increasing He to 17 m (Figure 8). This result is consistent with the approach that assumes the lateral earth pressure in braced excavations is trapezoidal. Another possible explanation is console-like behavior and a rapid decrease in wall lateral displacement toward the bottom.

3.9. The Effect of Surcharge Loads on the Lateral Deformations of the Walls

Table 9 shows the maximum UxL and UxR results from the analysis performed for the different symmetric and asymmetric surcharge conditions for LS and DS considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, d = 40 cm, Rinter = 0.6, Df = 16 m, and B = 45 m. The UyL values are almost twice as large in LS soils compared to DS soils for symmetrical load cases, but this difference gradually decreases as the surcharge load acting on WR decreases against a constant surcharge load on LR, and the asymmetry ratio increases.

The effects of surcharge loading with asymmetric conditions on lateral deformations at WL and WR along the walls for DS considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, d = 40 cm, Rinter = 0.6, Df = 16 m, and B = 45 m is shown in Figure 9. In the case of symmetrical loading, both UxL and UxR increase with the increase in surcharge load. However, in the case of asymmetric loading with a constant value of qR = 0 kPa, UxL increases while UxR decreases with the rise of qL. For an asymmetric load with a constant value of qL = 60 kPa, UxL decreases while UxR increases with the increase of qR. For symmetric conditions, the load increases the horizontal deformation of the wall, but while the asymmetry ratio increases because the load behind only one wall decreases, the deformation of the wall with decreasing load behind it decreases, whereas the deformation of the other wall increases even though the load behind it is constant. Therefore, when evaluating wall displacement, the influence of asymmetry, as well as the load, should be considered.

Except for a few cases for WR, less deformation occurred in asymmetric loading cases than in the unloaded case. It can be seen that the walls can affect each other positively and negatively through the struts. In addition, the correlation between the maximum UxL and UxR and surcharge loads can be expressed below in Equations (7)–(12) for LS and Equations (13)–(18) for DS. It should be emphasized that these equations will only produce valid results for loose and dense sandy soils with the same or close values of wall and sturt properties employed in the present paper, and for surcharge loads of less than 60 kPa.

3.10. The Effect of Surcharge Loads on the Bending Moments of the Walls

Table 10 shows the maximum ML and MR results from the analysis performed for the different symmetric and asymmetric surcharge conditions for LS and DS considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, d = 40 cm, Rinter = 0.6, Df = 16 m, and B = 45 m. Figure 10 shows the effects of surcharge loading with asymmetric conditions on the bending moments at WL and WR along the walls for DS.

Both ML and MR increase with increasing surcharge load for both symmetric and asymmetric loading and it is more significant in LS compared to DS. In addition, the maximum values of ML and MR are more than twice as high at LS as at DS. Asymmetry in the loading causes an increase in the bending moment value for both walls as the asymmetry ratio increases. Even on the wall where the load does not change, the change in bending moment shows the interaction of both walls in braced systems and the importance of evaluating them together.

In addition, the correlation between the maximum ML and MR and surcharge loads can be expressed in Equations (19)–(24) for LS and Equations (25)–(30) for DS. It should be emphasized that these equations will only produce valid results for loose and dense sandy soils with the same or close values of wall and sturt properties employed in the present paper, and for surcharge loads of less than 60 kPa.

3.11. The Effect of Surcharge Loads on the Surface Soil Settlements and Heaves

Table 11 shows the maximum UyL and UyR results for the different symmetric and asymmetric loading conditions for LS and DS considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, d = 40 cm, Rinter = 0.6, Df = 16 m, and B = 45 m. Also, Figure 11 shows the effects of surcharge loading with asymmetric conditions on the surface settlements behind the WL and WR along the walls for DS.

While the asymmetry ratio increases because the surcharge load behind only one wall decreases, the deformation of the wall with decreasing load behind it decreases, whereas the deformation of the other wall increases even though the load behind it is constant. Therefore, when evaluating surface settlements, the influence of asymmetry, as well as the load, should be considered. Even on the wall where the load does not change, the change in surface settlements shows the interaction of both walls in braced systems and the importance of evaluating them together.

Furthermore, the correlation between maximum UyL and UyR and surcharge loads can be expressed in Equations (31)–(36) for LS and Equations (37)–(42) for DS. It should be emphasized that these equations will only produce valid results for loose and dense sandy soils with the same or close values of wall and sturt properties employed in the present paper, and for surcharge loads of less than 60 kPa.

Table 12 shows the maximum Hv for the different symmetrical and asymmetrical surcharge conditions for LS and DS considering that the model size is 700 m horizontally and 350 m vertically and the mesh size is fine, d = 40 cm, Rinter = 0.6, Df = 16 m, and B = 45 m. The maximum Hv values are not affected by the loading conditions and are almost twice as large for LS as for DS.

4. Conclusions

The current study intends to evaluate the link between surcharge loads and changes in strut loads, deformations, and bending moments on walls, soil surface settlements, and excavation heave in the deep excavations. Using PLAXIS 2D v20, a finite element-based computer program, deep excavation pits with diaphragm walls on loose (LS), and dense (DS) sandy soils are modeled for various symmetric and asymmetric loading scenarios. The maximum values of Tst, UxL, UxR, ML, MR, UyL, and UyR for both the left and right sides of the excavation, and the heaves (Hv) in the excavation are estimated using linear and polynomial second-order relationships using the SPSS Statistics v17 computer program.

The following conclusions were drawn from the study:(i)With increasing surcharge load, both for LS and DS, under symmetric and asymmetric loading, the maximum of Tst rises.(ii)Both UxL and UxR rise when the surcharge load increases in the case of symmetrical loading. In asymmetric loading, while the load on one side is constant, the increase in the load on the other side increases the deformation of the wall on the side where the load increases, while it reduces the deformation of the wall on the side where the load remains constant.(iii)ML and MR increase with increasing surcharge load for both symmetrical and asymmetrical loading.(iv)The maximum Hv values are unaffected by the loading conditions.(v)The ground surface settlements increase as the surcharge load increase in symmetrical load scenarios. In asymmetrical loading scenarios, the swelling increase on the side on which the smaller load acts, whereas the settlement increases on the side which on the higher load as the asymmetry ratio increases.(vi)The behavior of the retaining walls on the left and right sides are influenced by each other, therefore the strut-supported excavation systems should be designed as a whole, but not individually.(vii)It should be emphasized that these equations only produce valid results for loose and dense sandy soils with the same or close values of wall and sturt properties employed in the present paper, and for surcharge loads of less than 60 kPa.

Consequently, it has been demonstrated that the estimating models produced by regression analyses and the outcomes from the finite element approach agree with one another. This gives promise that by broadening the study’s scope in the future, viable models can be created using multiple regression analyses. All the results obtained are limited to the scope of the study, and it is necessary to continue and expand the studies to give more inclusive outcomes.

Data Availability

The data supporting the findings of this study are available within the article.

Conflicts of Interest

The author declares that there is no conflicts of interest.