The existing literature has an abundance of research works on modeling of circular reinforced concrete columns confined with FRP wraps. The compressive behavior of concrete cylinders externally confined with FRP under axial loads has been extensively studied nowadays. Numerous stress-strain models have been established during the last three decades. Some of these studies (Vintzileou and Panagiotidou, 2008; Jiang and Wu, 2012; Rousakis et al., 2012; Ozbakkaloglu et al., 2013, 2016; Rousakis and Tourtouras, 2015; Hany et al., 2016; Farahmandpour et al., 2017) concentrated on modeling the conditions during failure of FRP-confined concrete, including the ultimate compressive strength and the corresponding ultimate axial strain. Other studies simulate and describe the overall behavior of the stress-strain curve (Hosotani et al., 1997; Samaan et al., 1998; Spoelstra and Monti, 1999; Lam and Teng, 2003b; Yu et al., 2010a,b; Wang et al., 2011; Ozbakkaloglu and Lim, 2013; Eid and Paultre, 2017; Fahmy et al., 2017; Ismail et al., 2017). The available stress-strain models in literature have been grouped into two categories by Lam and Teng (2003a): (a) design-oriented models and (b) analysis-oriented models. Analysis-oriented models provide the stress-strain behavior using closed-form equations, while design-oriented models achieve that by an incremental numerical procedure. Such an analysis-oriented model for FRP-confined circular RC columns can be found in Megalooikonomou et al. (2012) and Papavasileiou and Megalooikonomou (2015). The advantage of design-oriented models over analysis-oriented models is their simplicity in application. Hence, they can be easily integrated to a structural analysis software, where they will yield results fast. To the authors' knowledge, such an analysis-oriented model for rectangular FRP-confined columns is not yet available. The proposed model was developed with the intention to be integrated in existing structural software. Furthermore, the approach of this model is unique. While other available models intend to capture the stress-strain behavior of FRP-confined concrete, the proposed model intends to simulate effectively the axial and lateral strain and, consequently, simulate the stress-strain behavior. The confining stress occurs as a reaction to the developed strain. This way, in addition to simulating the stress-strain behavior, the proposed model also yields the lateral dilation, which is not explicitly calculated in available models.

Literature on the modeling of rectangular FRP-confined RC columns is more limited. While existing studies (Karabinis and Kiousis, 1996; Teng et al., 2002; Lam and Teng, 2003a; Teng and Lam, 2004; Roussakis et al., 2008; Piscesa et al., 2018) have verified that FRP confinement can substantially enhance both the compressive strength, ductility, and energy dissipation of confined concrete in circular RC columns, the same method has been found to be much less effective for rectangular RC columns (Mirmiran et al., 1998; Rochette and Labossiere, 2000; Lam and Teng, 2003a; Megalooikonomou, 2007; Karabinis et al., 2008; De Luca et al., 2011). Corner rounding is overall recommended to enhance the confinement effectiveness in a rectangular RC column and to reduce the detrimental effect of sharp corners on the performance of FRP jackets which causes their rupture in relatively low strain (Lam and Teng, 2003a).

The main difference between the two RC column section cases is that in circular sections the confinement is uniform, while the same does not apply on rectangular sections (Mander et al., 1988). In the later, concrete is non-uniformly confined, so the effectiveness of the confinement is substantially reduced.

As illustrated in Figure 1, the imposed axial load on a circular section causes radial dilation (i.e., uniform strain) on the confining elements so the pressure applied by the jacket to the concrete as a reaction to its dilation is uniform as well. Hence, the stress-state developed in a circular section is uniform. In rectangular sections the confining element develops a concentration of stresses at the corners, where it develops its full capacity. Along the sides, its confining effectiveness is reduced as it does not have enough stiffness to fully restrain the dilation of the concrete column. Thus, a non-uniform confining stress-state is developed in the concrete section.

Previous modeling approaches concentrated mainly on dividing the rectangular section in a confined and an unconfined area (Figure 1), based on the concept of possible arching effect (e.g. Figure A1 in EN1998-3{A.4.4.2}). Then, the confined area is considered to be in a state of uniform biaxial confinement, as in the circular cross sections, thus allowing the use of formulas defined for circular FRP-confined elements. The unconfined part is considered to be unaffected.

Considering all the above remarks, an iterative procedure is being proposed in this research study based on the results of a 3D FEM analysis performed by the authors which indicates that the arching effect does not apply. The unconfined areas shown in Figure 1 are indeed partially confined and they contribute to the overall capacity of the section until they reach their maximum strength which is significantly reduced compared to that of the confined areas. In this work, a simplified analytical procedure to model the stress-strain behavior of a rectangular concrete section under concentric load is proposed which employs a system of “generalized” springs, well-known stress-strain laws and a failure criterion.

Figure 2 (Campione and Miraglia, 2003) shows the cross-section of an experimentally tested FRP-confined square concrete specimen under concentric load and the final concrete state along with the effective concrete core after FRP failure at the corner at the end of the experiment. Uneven damage can be observed throughout the section. Two different regions can be identified having different confining stress states. To determine the confining stresses and define the confined and unconfined regions in a rectangular section, a 3D Finite Element Model was developed in SAP 2000 (Computers Structures Inc., 2016). The model consists of a square concrete section 200 × 200 mm (Figure 2) with rounded corners. Solid elements are used to model concrete in the section. The FRP wrap is modeled using shell elements applied on the solid elements in the perimeter. A typical slice with thickness 10 mm was simulated. The CFRP wrap was modeled using an orthotropic material model with Young's modulus 3.61 10⁵ MPa. The confined concrete was modeled with a “concrete” material with f_ck = 17 MPa. The applied axial displacement on the FEM model was increased until the lateral strain at the slice was equal to the maximum lateral strain at a cylindrical specimen, when it reaches the maximum strength. The maximum axial displacement was defined so that the lateral strain developed in the slice is equal to that an unconfined cylindrical specimen at the ultimate stress. The FEM results were plotted in three-dimensional graphs (Figure 3) to allow further study of the stress field. Verification of the results yielded by the FEM model against experimental results is available in Teng et al. (2015).

Graphs in Figure 3 depict various stress fields in a quadrant of this square section. Stresses near the rounded corner are not shown since they represent local stress state. In Figures 3A,B, the normal stresses parallel to the diagonals (SD1 and SD2) are presented, while in Figures 3C,D normal stresses parallel to the section's sides (SO1 and SO2) are shown. Stresses parallel to the diagonal of the rounded corner, are increased near the corner, but moving inwards to the center of the section they reduce significantly (Megalooikonomou, 2007; Nisticò and Monti, 2014). Figure 3 also illustrates that the stresses perpendicular to the rounder corner diagonal direction and close to the sides have much lower values compared to those parallel to rounded corner diagonal, while toward the center they both become almost equal. This confirms that close to the center the confinement stress state is similar to circular sections. It is apparent from the plots of the stress field parallel to the orthogonal directions (Figures 3C,D) that some confining stresses are present along the side, while in the central part of the perpendicular side they are close to zero, as the confining device has minimal flexural stiffness.

Based on the aforementioned, the following main remarks can be made:

Accepting some tolerance, the regions where a biaxial and a uniaxial confinement exist can thoroughly be defined based on the ratio of the principal stresses of the two perpendicular directions in the joints of the FEM. Regarding the stress output in SAP 2000 (Computers Structures Inc., 2016), it should be noted that the direction of the middle principal stress (S_mid) is perpendicular to the maximum (S_max) and minimum (S_min) principal directions.

For a specific tolerance (i.e., the ratio of the confining stresses is < 15%), the width of the biaxial stress state region is independent of the stiffness of the FRP (Figure 4). The radius of the rounded corner affects more the diagonal dimension of this region, while parallel to the sides the width remains the same (Figure 4). Unlike reinforced concrete, where the spacing of stirrups allows for such regions to develop, not only on the concrete cover but also between two consecutive stirrups, FRP-confined members are totally inside the FRP wrap. Due to the FRP's non-infinite stiffness, the FRP wrap can deform laterally to the cross-section. A finite slice of the FRP wrap performs as a cantilever beam with its ends at the rounded corners of the concrete section. The larger the deformation might get (i.e., a weaker FRP slice), the smaller the stress applied by the FRP on the confined concrete core is. No matter how small it might get, there is always stress applied on the concrete, so it does not perform as unconfined-concrete. The model considers the confinement pressure to be the same around the column. Hence, the smaller the side which this pressure distributes on, the more effective the confinement provided. This leads to the conclusion that, unlike unconfined concrete, in FRP-confined concrete the part of the section that is under biaxial confinement is not proportional to the relevant side, but proportional to the opposite side. This assumption is confirmed in previous works (Karam and Tabbara, 2005; Megalooikonomou, 2007). Based on these remarks, the width of the uniaxially confined region can be calculated using Equation (1), where a and b are the width and the depth of a rectangular section accordingly (Figure 8A). The dimensions of the regions are then directly related and determined according to this width (Figure 10).

The equations used in the proposed model was initially developed based on results yielded by modeling square FRP-confined concrete columns. Having assessed its effectiveness in modeling square columns, its applicability on rectangular columns was also assessed against experimental results and found to be adequate for a simplified analytical model.

A simplified mechanical model which describes the stress-strain behavior of a rectangular concrete cross-section under concentric load is proposed in this Section. A series of “generalized” springs (Figure 5) is used to describe the confinement mechanism. Compressed concrete expands laterally according to its confinement state. Such expansion activates the confining elements. The confining forces are applied to the section corners and transferred along the diagonal. In the proposed model, the contribution of the stirrups, as it is considered to be minimal and is not taken into account. This refers to cases that the stirrup spacing is adequately large to minimize the confinement effect due to stirrups.

Figure 5 shows the springs used to model the axial and the lateral behavior of the section. Axially, the springs receive the same displacements and their cumulative strength is the summation of the strength of the individual springs (parallel system). Laterally, the springs receive the same force and their displacements are added (series system).

The “generalized” springs illustrated in Figure 5 use the constitutive law for concrete proposed by Pantazopoulou and Mills (1995). The model which relates the volumetric strain to the axial strain is shown in Figure 6 (right). In the same figure (left), the corresponding axial stress vs. axial strain is also demonstrated. It is evident that adequate confinement increases the ductility of the RC member which develops considerably larger strain until its stiffness and strength are particularly reduced leading to loss of stability. This is shown in the relationship between volumetric strain ε_v and axial strain ε_c illustrated in Figure 6.

The initial slope of the curve (Equation 2, for εc≤εllim) is characteristic of a perfectly elastic condition. Both curves (for confined and unconfined concrete) deviate from this idealized situation. However, confined concrete develops large volumetric strain in much larger axial strain than unconfined concrete. Experimentally, in conditions of high confining stress, it has been observed that the ε_V vs. ε_c curve might even remain at negative values of ε_V throughout the test.

Prior to evident surface cracking, the ε_V–ε_c relationship is practically linear, with the lateral strain ε_l being equal to v·ε_c (v is typically within the range of 0.15–0.25). Beyond the limit axial strain (ε_c = ε_c,lim) that corresponds to lateral strain ε_l in excess of the tensile cracking of concrete ε_cr, the relationship between ε_V and ε_c shows substantial deviation from the idealized linear response and appears to be well approximated by a parabolic expression. For confined concrete under uniaxial load, the model is:

The product of α_ε · ε_co is the compressive axial strain at zero volumetric strain. For standard concrete strength, it is typically observed at axial strain from 2 to 3.5%0, i.e., at 80–100% the strain at peak stress ε_co. Coefficient b_ε denotes the degree of passive confinement of the concrete. For unconfined concrete loaded uniaxially, b_ε receives a value of 1, while smaller values are used as passive confining pressure increases. Coefficients α_ε and c_ε increase for higher strength concrete. Higher values of α_ε are used for higher nominal strength, up to 1, whereas the post peak response becomes more brittle (c_ε tends toward or exceeds 3). The same behavior (linear and parabolic part) has also been observed in the relationship between volumetric strain and axial stress (Figure 7).

To comply with the mechanical model of the “generalized” springs and the regions with varying confinement, this model has been modified to correlate volumetric strain ε_V to axial stress σ_c according to Equation (3) (axial stress in MPa). Initially, change of volume occurs due to compaction and is practically linear up to the point of critical stress α_σ·f_co (unconfined concrete strength, usually α_σ is taken as 0.7). For this axial stress level the Poisson's ratio ν remains within the range of 0.15–0.25 (here, the initial Young's modulus of concrete is determined as well). At this point, volume change is reversed resulting in volumetric expansion called (near-strength or at-peak-strength) dilatancy. A point can be defined where the compression rate of the specimen equates the expansion rate, thus resulting in zero volumetric strain. This point is considered to appear when the ultimate strength of the uniaxially confined region is reached (biaxial stress state, b_σ = 1.2) (Kupfer et al., 1969). After the deterioration of this region, the expansion rate increases faster than the compression rate (second order parabola, c_σ = 2) due to reduced effective confinement. The expansion becomes unstable during the crushing phase beyond the ultimate strength.

Based on Figure 8A, from the volumetric strain ε_V, both area strain (ε_A) and side strains (ε_a and ε_b) can be calculated as shown in Equations (3–7) (compressive axial strains are taken as negative).

Based on the commonly applied assumption that no friction is developing between the FRP jacket and the concrete surface, the elongation strain of the sides can be assumed equal to the jacket strain. Therefore, the diagonal force of the jacket applied laterally from the corners to the springs in series can be determined by projection (Figure 8B):

The coefficient k_e is a confinement efficiency factor by Karam and Tabbara (2005) which takes into account the increasing confinement effectiveness as the corner radius increases and is reduced when the ratio of the larger side to the smaller side of the cross section is increased. Figure 9 shows the considerations made for the definition of this factor.

In Figure 9 (left), a confining FRP wrap is shown acting on a generic rectangular section with long side a, a short side b and radius of the rounded corner R. The concrete is assumed to be subjected to average (uniform) confining stresses at its middle sections with f_a acting along the long side and f_b acting along the short side as illustrated in Figure 9. The relationship between the f_a, f_b and the jacket tensile stress f_j is determined statically (t is the thickness of the FRP wrap).

In Figure 9 (right), the FRP wrap is considered to act as a cable around the corner similarly to a pulley. Assuming no friction between the FRP jacket and the concrete surface, the relationship between f_j and the confining stress at the corner f_r is also defined statically:

Combining Equations (9) and (10), the following Equation can be defined:

Hence, the sharper the corner radius is, the higher the confining stress at the corner with respect to the average confining stress inside the cross section. This is supported by finite element analysis results presented by Parvin and Wang (2001) and experimental findings by Parvin and Wang (2002), and Chaallal et al. (2003).

Based to the aforementioned, the ratio of the average confinement stress over the maximum confinement stress attained in the cross section can be defined as the geometric confinement effectiveness factor (k_e). In a rectangular section, the maximum confinement stress (f_r) develops at the corners where stress concentration occurs on the jacket due to dilation of concrete. The average confinement stress is the average of the individual stresses f_a and f_b which act at the center of the section. The confinement effectiveness factor (k_e) can be calculated as:

In a circular cross section, this factor is equal to 1, while for a square cross section, it is: k_e = 2R/a. For an elongated rectangular cross section with semi-circular ends b ≫ a and 2R = a, the factor k_e tends asymptotically to 0.5, which corresponds to confinement in a single direction in the cross-sectional plane.

After the determination of the diagonal force applied to the lateral springs in series from the corners, the lateral pressures for each region can also be calculated (Figure 10).

The springs are in series so they develop the same force. For the biaxial stress state region, the force starts from the corners and moves parallel to the sides. Therefore, the confining pressure (assumed uniform) can be determined by the following equation:

The stress-strain model proposed by Popovics (1973) modified by Mander et al. (1988) is used to describe the behavior of the triaxial stress state region. Based on the lateral pressures calculated according to the Figure 10 the corresponding axial stress can be determined by the use of a stress-strain model corresponding to the confinement stress state of this region. The equation of the model by Mander et al. (1988) for the maximum axial stress σ_cc is not used in this case due to the fact that it describes the performance of uniform biaxial confining pressure. For the region in triaxial stress state, the uniform confining pressures on the sides of the section can be determined based on the geometry of the region (Figure 11) as follows:

To comply with the above modeling of the triaxial stress-state region where the confining pressures are different in the two lateral directions (only in the case of a square section they are the same), a failure criterion where all the lateral confining stress state cases are considered is applied. Specifically, the failure surface by Ottosen (1977) is used (Figure 11). The latter failure criterion corresponds to a smooth convex failure surface with curved peaks. This surface expands in the negative direction of the hydrostatic axis, while its projection to the deviatoric plane (perpendicular to the hydrostatic axis) changes from nearly triangular to a more circular as the hydrostatic pressure increases. The behavior of concrete can be modeled using octahedral normal (hydrostatic) and shear (deviatoric) stresses σ_o and τ_o. Any point in the stress space is described by the coordinates (ξ, ρ, θ), in which ξ is the projection in the hydrostatic axis (σ₁ = σ₂ = σ₃) and (ρ, θ) are the polar coordinates in the deviatory plane.

The following equations describe the surface. Having the lateral confining pressures of this region the value of maximum axial stress σ_{cc, triaxial} can be reached through iteration.

where:

The remaining parameters of the model are calibrated by the values proposed in the original work of Ottosen (1977) for different ratios of tensile concrete strength over compressive concrete strength (k = f_to /f_co) given in (Table 1).

Based on the σ_{cc, triaxial} determined above, the stress-strain law by Mander et al. (1988) is applied. The following equations describe the model's behavior:

(Richartet al., 1928)

The model for concrete under biaxial stress-state by Liu et al. (1972) is used for the biaxial stress state region:

where:

Note: For a₁ > 1, σ_c,biaxial is constant, equal to 1.2·f_c0.

Based on the areas of the different regions the total averaged axial stress of the cross section can be obtained:

The proposed model (Figure 12) is an iterative procedure where an assumed value of axial stress corresponding to an imposed axial strain is brought to convergence. After converging the assumed axial stress calculated with the above considerations, the resulted elongation strain through the iterative procedure should be compared to the ultimate rupture strain of the jacketing. It has been observed from experimental results that the average failure strains of the FRP wraps are of the order of 50–80% of the failure strain of the tensile coupons made from the same material and tested before the application of the material. This actual value of factor k (ranging between 50 and 80%) depends on the type of FRP used (Lam and Teng, 2003a).

The proposed algorithm's efficiency was assessed against small scale experiments, the results of which are presented in Megalooikonomou (2007) and Megalooikonomou et al. (2007). In this paper, further assessment against experimental results in large scale experiments takes place. Recent experimental data collected by large-scale tests on FRP-confined rectangular RC columns (Zeng et al., 2018) were used for the simulation of monotonic loading. Zeng et al. (2018) presented the test results of an experimental study consisting of nine large-scale rectangular RC columns with a cross-section of 435 mm in depth and 290 mm in width, including eight FRP-confined RC columns and one RC column without FRP jacketing as the control specimen, tested under axial compression. The experimental program examined the sectional corner radius and the FRP jacket thickness as the key test variables. The proposed algorithm was assessed against three of these specimens, i.e., the specimens with corner radii 25 or 45 mm and one or two layers of Carbon Fiber Reinforced Polymer (CFRP) wrap. The corner radius of 65 mm was not considered in the assessment, since in buildings designed using obsolete design codes the concrete cover thickness is typically small. Hence, a cover of at least 65 mm which would allow the formation of the round corners in such a column is highly unlikely to be found. Also, Zeng et al. (2018) assessed the effectiveness of material models available in literature against the defined experimental results.

The proposed model does not take into account the contribution of the stirrups, as their spacing is considered to be adequately large to minimize their confinement effect. This is in accordance with the results presented by Zeng et al. (2018), where the contribution of the reinforcement is also not taken into account due to the large stirrup spacing. Figure 13 shows the correlation of the proposed material model with these experimental results and overall the numerical response can be characterized as satisfactory. It should be noted that for this comparison the iterative procedure was terminated for FRP rupture strain equal to 50% of that of the experimentally tested tensile FRP coupons (Zeng et al., 2018).

A further evaluation of the model's performance under cyclic lateral loading which simulates earthquake loads and simultaneous constant axial compression was performed. This was achieved by comparison against the experimental tests of FRP-retrofitted square RC columns performed by Memon and Sheikh (2005). This experimental study evaluates the effectiveness of Glass Fiber Reinforced Polymer (GFRP) wraps in strengthening deficient and repairing damaged square RC columns. Each of the eight specimens tested, representing columns of buildings and bridges constructed before 1971, consisted of a 305 × 305 × 1473 mm column connected to a 508 × 762 × 813 mm stub. Specimens were tested under constant axial compression and cyclic lateral displacement excursions simulating earthquake loads.

The modeling of these FRP-confined square RC columns has been performed using the MatLab toolbox FEDEAS lab “Finite Elements for Design Evaluation and Analysis of Structures” (Filippou and Constantinides, 2004). The experimental moment-curvature responses within the plastic hinge regions are reported along with the numerical results in Figure 14. The simulation of the cantilever columns has been applied using a unique fiber beam-column element (Spacone et al., 1996) with force formulation for the entire column, in which the FRP-confined concrete was modeled using the proposed material model with degraded linear unloading/reloading stiffness according to the work of Karsan and Jirsa (1969) and no tensile strength. The constitutive model by Menengotto (1973) is used to model the longitudinal steel behavior. The moment-curvature response of the most critical fiber section of the applied non-linear fiber element was then reported. It can be seen that the agreement is very close to the experimental one, with some deviation concentrated on the parts of reloading after reversal of the imposed displacement. This difference in response in terms of modeling can be explained based on the way the cracks on the concrete surface are described in the level of the material model. Because the crack is described as a two-event phenomenon (open or closed cracks), when the longitudinal steel reinforcement is in compression and the crack is closing, the concrete contributes to the total strength of the column, creating this deviation in the response. In reality, this is not the case due to imperfect crack closure.

An iterative approach was proposed to model both the axial and lateral stress-strain response of axially loaded FRP-confined rectangular and square reinforced concrete columns.

In FRP-confined square or rectangular sections, no unconfined concrete regions are observed, as assumed in many models. These sectors along the sides between adjacent corners are confined from forces coming from the corners and moving parallel to the sides. Therefore, the areas where arching effect is assumed in the section are actually partially confined, so they contribute to the column's total strength until their maximum strength (which is lower than the inner part of the section). Thus, two different regions with different confining stress-states are identified.

The two regions are uniaxially and biaxially confined (biaxial and triaxial stress-state, respectively). Therefore, the contribution of each region to the total section strength can be modeled as a system of parallel springs, whose axial stresses are added based on the corresponding constitutive law under biaxial or triaxial stress state.

The lateral behavior develops along the diagonals of the section and can be represented by a system of springs in series. It was shown that both sides' lateral strains in the rectangular sections are equal, regardless of their aspect ratio. The reacting force of the confining device applied from the corners can be shared among the regions based on the defined path of the confining forces and the geometry of the regions. The resulting lateral uniform pressures lead to the corresponding axial strength of the regions.

The algorithm takes into account all parameters available to designers, such as corner rounding radius, stiffness of the FRP, and concrete strength, while it can be easily understood and implemented. Its results are found to correlate adequately to recent experimental data yielded by large-scale tests on FRP-confined rectangular RC columns. Assessment against small-scale experimental results can be found in Megalooikonomou (2007). Finally, the performance of this material model was further investigated by its implementation to the simulation of a series of experimental tests of FRP-retrofitted square RC columns under cyclic lateral loading simulating earthquake loads and simultaneous constant axial compression. In particular, all specimens were simulated using non-linear fiber elements, in which the FRP-confined concrete was modeled using the developed material model.

Comparison between the numerical and experimental hysteresis of the column is indicative of the effectiveness of the implemented modeling. The model is effective in the modeling of square FRP-confined concrete sections, while it can be applied on the modeling of rectangular sections as well. Its effectiveness can be reduced for large ratio of sides. In particular, for large ratios of the larger side over the smaller side of the column, the proposed model seems to simulate effectively the initial branch of the stress-strain curve, but deviation from the experimental results was observed at the second branch, for larger deformations. Furthermore, it can be applied for various types of FRP wraps such as GFRP and CFRP. It is effective for a large range of FRP thickness that is expected to be met in practice. For sections with large rounding radius at the corners, the proposed model seems to overestimate the contribution of the confinement to the specimen's overall capacity. However, it should be pointed out that such cases are not expected to be met in practice, as the reduced concrete cover sizes defined by the previous design codes, do not allow for a large radius of rounding to be applied on existing reinforced concrete columns. Even though the confinement effectiveness factor (k_e) considers the case of a large radius, use of this model in circular sections is not recommended. For circular FRP-confined sections, the readers are referred to Megalooikonomou et al. (2012) and Papavasileiou and Megalooikonomou (2015).

KGM developed the iterative procedure and performed the comparison and interpretation. GSP contributed to the development of the cyclic constitutive law and its implementation in performing cyclic non-linear static analyses for comparison with experimental results. All the authors read and approved the final manuscript.