Elasto-Plastic Analysis of Two-Way Reinforced Concrete Slabs Strengthened with Carbon Fiber Reinforced Polymer Laminates

: This study explores a technique for enhancing the punching strength of reinforced concrete (RC) flat slabs, namely carbon fiber reinforced polymer (CFRP). Four large-scale RC flat slabs were fabricated, to assess the efficacy of this strengthening method. One slab served as a reference and the three other specimens were strengthened with CFRP, as a method of external strengthening. These slabs, featuring identical overall dimensions and flexural steel reinforcement, underwent testing until failure, under the influence of concentrated patch loads. A concrete plastic damage constitutive model (CDP) was developed and employed to examine the strength of two-way RC slabs. Additionally, to enhance the strength of existing RC slabs, carbon fiber reinforced polymer (CFRP) strips are affixed to the tension surface of the sections. The research begins with the calibration of a numerical model, based on data from laboratory tests. The objective of this study is to constrain the plastic behavior of two-way RC slabs reinforced with CFRP, with a focus on establishing an optimal elasto-plastic analysis, aimed at controlling concrete damage plasticity using CFRP, and employing a plastic limit load multiplier. Subsequently, a series of numerical simulations, incorporating different variables, are conducted to investigate shear behavior. The numerical results indicate that an increase in the strengthening ratio has a significant impact on shear strength. Finite element simulations are carried out using Abaqus CAE ® /2018.


Introduction
Recently, the use of lightweight carbon fiber reinforced composite materials that are nonmetallic has increased for the purpose of repairing and reinforcing concrete structures."Strengthening a structural member" is a process undertaken when the existing member's strength is deemed inadequate.Such deficiencies can manifest due to construction errors, such as excessive downward pressure on the steel bars into the cross section during the building process, alterations in the intended capacity, or reductions in capacity caused by human error or environmental factors [1].Civil engineering has witnessed the economical and effective implementation of CFRPs for the development and rehabilitation of deteriorating structures.Numerous variables affect the application of CFRPs for the reinforcement of RC members, including fiber type, fiber distribution, fiber orientation, and bond schemes [2].
Many researchers investigated CFRP reinforcement methods including laminates, bars, or sheets.There are numerous benefits of CFRP composites over conventional materials.The construction industry is progressively adopting these benefits for infrastructure applications-durability, transportability in any length, exceptionally high strength, lightweight, requiring less mechanical fastening, requiring less expensive maintenance, requiring shorter construction periods, minimizing interruptions, having a long service life, and facilitating rapid installation [3,4].Limam et al. [5], in 2003, provided an experimental investigation and a limit analysis modeled on RC slabs reinforced with CFRP laminate.The findings from the experimental study indicated that externally applying the CFRP plates effectively enhanced the strength of RC slabs.Mosallam and Mosalam [6], in 2003, examined the ultimate response of reinforced and unreinforced concrete slabs that had been retrofitted and repaired with CFRP laminates, through an experimental and analytical investigation.The carbon fiber/epoxy and E-glass/epoxy composite systems significantly strengthened the repaired areas.Additionally, retrofitting slabs with CFRP systems improved structural integrity.Limam et al. [7], in 2005, used 2D finite element analysis to quantify the elastic response of simply supported two-way reinforced concrete slabs with CFRP laminates, concentrating on orthotropic composite plates with a central load.Hosseini [8], in 2016, conducted an experimental program and used computational techniques to evaluate the load-deflection properties of reinforced concrete slabs, highlighting the effects of CFRP percentage, tensile longitudinal bar percentage, and CFRP elasticity modulus.The findings showed that prestressing CFRP laminates with NSM can improve RC slab fracture resistance, serviceability, yield capacity, and ultimate loads.Yılmaz et al. [9], in 2018, utilized a combination of experimental and numerical methods to examine the impact behavior of two-way RC slabs reinforced with CFRP strips.The findings indicated that reinforcing RC slabs with CFRP strips increased their rigidity and toughness and the maximal acceleration values measured from test specimens correspondingly rose.Akhundzada et al. [10], in 2019, used unidirectional CFRP laminates to test the efficiency of various strengthening techniques for internal slab-column connections and examined how laminate thickness and anchoring affect strengthening.The findings revealed that strengthening decreased maximum deflection by 50%, increased ultimate load by 25%, increased cracking load, and significantly reduced radial crack length.Mercimek et al. [11], in 2022, examined RC flat slab efficacy-enhancing methods.The strengthening methods were CFRPs, shear stud, and textile-reinforced mortar (TRM).This study used Abaqus CAE ® to simulate the experiment and compare the results to the actual results.The results showed that TRM was the best out of the three strengthening methods used to improve RC slab performance under punching loads.Türer et al. [12], in 2023, devised a CFRP strip reinforcing technique to mitigate punching strength loss from two square holes close to the column in two-way RC flat slabs.The specimen's initial stiffness, load-displacement characteristics, energy dissipation capabilities, maximum load-bearing capacities, and CFRP-enhanced failure patterns were studied at varied opening sizes and positions.For numerical analysis, Abaqus CAE ® was used to model the specimens.
On the other hand, a lot of researchers used the other types of CFRPs (sheet or bar) in strengthening the two-way reinforced concrete slabs.Enochsson et al. [13], in 2007, examined the structural response of two-way RC slabs reinforced with CFRP sheets, with sawn-up holes, under uniform loads.Both numerical and experimental methods are used to evaluate CFRPs' ability to absorb excess sectional forces.The results show that externally bonded CFRP sheets may support slabs with openings better than steel-reinforced slabs.Esfahani et al. [14], in 2009, investigated methods for using CFRP sheets to increase the flat slabs' resistance against punching shear.The test results underwent comparison with the formulas outlined in both the BS 8110 and ACI 318 Codes.The findings from this investigation align well with outcomes reported by fellow researchers, demonstrating that the addition of CFRP sheets to the slabs' tension side effectively enhances the slabs' punching.
Meisami et al. [15], in 2013, evaluated the potential of FRP rods to enhance shear strength, using an analysis of the punched shear performance of concrete slabs reinforced with epoxy resin and CFRP rods.The experimental findings showed that the selected reinforcement method improved the slabs' maximum load-bearing capacity and ability to withstand deformation, while also preventing possible brittle failures from vertical point loads.Elsanadedy et al. [16], in 2015, aimed to illustrate the efficacy of utilizing externally bonded CFRP composites for improving the flexural capacity of RC one-way slabs.The study encompassed both experimental testing and computational analysis, utilizing FE analysis.A comparison was conducted between the experimental and numerical findings, resulting in a satisfactory level of agreement.Behzard, Sharbatdar, and Kheyroddin [17], in 2015, showcased the efficacy of a novel technology that combines CFRP materials to enhance the flexural strength of existing RC two-way slabs with a thin layer of concrete covering.The test results have shown the feasibility and effectiveness of this technology in enhancing the flexural behavior of reinforced concrete two-way slabs.A 3D nonlinear numerical model was created using the FE approach, to accurately forecast the flexural behavior of the tested slabs.Mahmoud et al. [18], in 2021, studied the performance of high-strength RC slabs reinforced with CFRP sheets, to enhance their bending strength.The findings from the test were compared together to draw conclusions regarding the extent of CFRP enhancement; the load-carrying capability of the specimens improved after reinforcement using high-strength materials-thin concrete slabs.
The concrete damage plasticity (CDP) model is widely used in structural engineering to assess the behavior of concrete structures under various loads.The CDP model, a continuum-based model, takes into consideration the intricate behavior of concrete, such as cracking, crushing, and plastic deformation.This non-linear model is essential for designing and analyzing buildings that experience significant pressures, as it accurately represents the behavior of concrete after it reaches its peak strength.The CDP model was successfully applied in several situations, including earthquake engineering, blast-resistant building design, and reinforced concrete construction.This discovery might have important implications for the design and construction of reinforced concrete structures, especially in areas prone to earthquakes where ductility is a critical factor [19].Various combinations of damage and plasticity were explored by Carol et al. [20] in 2001, Ibrahim et al. [21] in 2023, and Rad et al. [22] in 2023.Gatuingt and Pijaudier-Cabot [23], in 2002, and Kratzig and Polling [24], in 2004, utilized isotropic damage in their research, elucidating various combinations of plasticity types.Lubliner et al. [25], in 1989; Imran and Pantaopoulou [26], in 2001; and Ananiev and Ozbolt [27], in 2004, explored theories that describe plasticity in nominal stress space.
This study aims to investigate the impact of varying strengthening ratios on the ultimate load, cracking load, crack width, and deflection.The recent studies did not investigate the behavior of RC two-way slabs with these specific ratios of CFRP.Furthermore, this study demonstrated the need to implement CFRP reinforcement to mitigate the damage to concrete and steel structures using Abaqus CAE ® .This research shed light on the performance of CFRP in RC two-way slabs and how its strength might be improved.Engineers and researchers can utilize the findings as a reference for determining the appropriate proportion of reinforcement ratios to employ, to attain certain structural objectives.As well as the fact that this study possesses various distinctive elements that distinguish it as a novel study, it begins by investigating two-way RC slabs experimentally, as well as carrying out an ideal elasto-plastic analysis that focuses on managing concrete damage plasticity using CFRPs and applying a plastic limit load multiplier, with the intention of controlling the plastic deformation of two-way reinforced concrete slabs that have been strengthened with CFRPs.To further our comprehension of how various strengthening ratios affect the behavior of two-way RC slabs, these percentages were derived by testing three specimens; then, their tension and compression behaviors were investigated and compared to the control slab, which did not contain CFRPs.Finally, experimental testing and numerical analysis were utilized in this study for a more thorough evaluation and validation of the findings, which increases the study's credibility and practicality.It was of the utmost importance to validate numerical models that were derived from the results of laboratory tests, with a particular emphasis on confirming the strengthening ratio.A sequence of simulations that are based on the limitations have been set.As a means of accurately representing the caused damage in the concrete, the numerical model was calibrated by making use of the CDP model.There was an evaluation of the effect of load multipliers on concrete deterioration, using both elastic and plastic analysis.All these parameters were considered and their effects were studied.
Following this introduction, Section 2 provides a comprehensive explanation of the concrete constitutive model.Then, the experimental work that includes the details of the model and the materials' properties are present in Section 3. The test set-up and discussion of the experimental data are presented in Section 4. Additionally, Section 5 presents specific information on the numerical model, including the specifics, calibration procedure, and discussion of the numerical findings.Finally, the most significant findings and conclusions are presented in Section 6.

Concrete Constitutive Model
Concrete's nonlinear reaction to stress, compression, and bending is modeled using the CDP model.The concepts of continuum damage mechanics state that loading damages and degrades material.The model considers compressive and tensile damage.Micro-cracking in the material causes strength loss in tensile damage and stiffness, according to the model.The model describes compressive damage as crushing, which reduces strength and stiffness.The model also accounts for plasticity, which causes the material to deform permanently once the stress is removed.With isotropic hardening, the model properly depicts the material's plastic reaction.The modulus of elasticity, which describes compressive and tensile concrete behavior inside the calibrated model, must be obtained experimentally to implement the CDP model.Monographs and the relevant literature provide detailed information on the procedure.This article briefly describes the concrete constitutive model and then elaborates on essential points.Elastic and plastic components may be examined in the strain tensor value (ϵij) using Prandtl-Reuss theory and elasto-plastic deformations, as explained below [19,28].
The equation for scalar damage elasticity provides us with understanding regarding the interplay between strains and stresses within a material.
Hence, D el ijkl serves as an abbreviation indicating a reduction in elastic stiffness.
where D el 0 represents the initial stiffness in elasticity of the material before any damage occurs, while "d" serves as a single degradation parameter ranging from a state of 0, indicating no damage, to a state of 1, representing complete damage.Following the hypothesis of scalar damage, a single degradation variable, denoted as d, is employed to characterize the reduction in stiffness in all directions.The effective internal stress is described utilizing fundamental concepts from continuum damage mechanics, as outlined below: The relation of scalar reduction correlates with the internal stress, σij , to the internal stress's effectiveness, When d equals 0, the effective internal stress, σ ij , coincides with the stress within the structure, σij .As damage occurs, the zone characterized by the effective internal stress, which opposes external loads, gains greater significance than the internal stress itself.Expressing the diminished elastic tensor and the integration of nominal stress, as defined in Equation ( 4), into a modified version of Equation (2) results in the following equation: The interaction between internal stress and strain constitutes the damage plasticity constitutive model.
The cyclical performance of concrete under uniaxial stress, characterized by the formation and sealing of micro-cracks, was explored to understand the intricate degradation mechanism.Damage models for concrete traditionally considered compressive crushing and tensile cracking, with d c and d t representing compression and tension damage ranging from a state of 0, indicating no damage, to a state of 1, representing complete damage, and σ t and σ c denoting internal stress that is efficacious under tension and compression, respectively.Figure 1 [19] illustrates the anticipated impact of plastic deformation on the concrete's reactions to uniaxial tension and compression.The interaction between internal stress and strain constitutes the damage plasticity constitutive model.

𝜎 ̂𝑖𝑗 = ( 1 − 𝑑 ). 𝜎 ̅ 𝑖𝑗 , 𝜎 ̂𝑖𝑗 = ( 1 − ⅆ 𝑡 ). 𝜎 ̅ 𝑡 𝑖𝑗 + ( 1 − ⅆ 𝑐 ). 𝜎 ̅ 𝐶 𝑖𝑗 (7)
The cyclical performance of concrete under uniaxial stress, characterized by the formation and sealing of micro-cracks, was explored to understand the intricate degradation mechanism.Damage models for concrete traditionally considered compressive crushing and tensile cracking, with ⅆ  and ⅆ  representing compression and tension damage ranging from a state of 0, indicating no damage, to a state of 1, representing complete damage, and  ̅  and  ̅  denoting internal stress that is efficacious under tension and compression, respectively.Figure 1 [19] illustrates the anticipated impact of plastic deformation on the concrete's reactions to uniaxial tension and compression.The manifestation of concrete's response to compression and tension stress in its uniaxial compressive and tensile behavior, as elucidated using the CDP model, is depicted by the follow equations: Here,  0 signifies the modulus of elasticity of the pristine material, while   ,ℎ and   ,ℎ denote the corresponding plastic deformations in compression and tension.This enables the derivation of the effective uniaxial compressive stress,  ̅  , and tensile stress,  ̅  .
The tensile strain,   , is equivalent to   ,ℎ +    ; the compressive strain,   , is equivalent to   ,ℎ +    .Consequently,    and    represent the equivalent elastic strains under compression and tension.Through the utilization of the expressions   =   ,ℎ +    and   =   ,ℎ +    for tensile strain and compressive strain, respectively, during tension and compression, the elastic strains are represented by    and    , respectively.Ultimately, the model that has been calibrated employs the CDP constitutive model to depict the concrete's behavior under both compression and tension.This depiction aligns with the properties of concrete in compression and tension, as illustrated in Figure The manifestation of concrete's response to compression and tension stress in its uniaxial compressive and tensile behavior, as elucidated using the CDP model, is depicted by the follow equations: Here, E 0 signifies the modulus of elasticity of the pristine material, while ϵ pl,h t and ϵ pl,h c denote the corresponding plastic deformations in compression and tension.This enables the derivation of the effective uniaxial compressive stress, σ t , and tensile stress, σ c .
The tensile strain, ϵ t , is equivalent to ϵ pl,h t Consequently, ϵ el t and ϵ el c represent the equivalent elastic strains under compression and tension.Through the utilization of the expressions ϵ t = ϵ pl,h t + ϵ el t and ϵ c = ϵ pl,h c + ϵ el c for tensile strain and compressive strain, respectively, during tension and compression, the elastic strains are represented by ϵ el t and ϵ el c , respectively.Ultimately, the model that has been calibrated employs the CDP constitutive model to depict the concrete's behavior under both compression and tension.This depiction aligns with the properties of concrete in compression and tension, as illustrated in Figure 1.
Subsequently, the validation of these properties was conducted using Abaqus CAE ® and, after running Abaqus CAE ® , the parameters for the CDP model were acquired, preparing the model for testing with various concrete mixes.

Utilizing the Static Principles of Plastic Limit Analysis
Incorporating the static framework of plastic limit analysis into the method of optimization enables the consideration of concrete slabs' plastic behavior under different loading conditions, ensuring that the designs that have been optimized meet the state of the plastic limit requirements.An in-depth study of the original development of the static limit analysis principle is essential to understand the fundamental notion of the newly created technique.Contemplating an elasto-plastic body subjected to a steadily increasing force, F i , to comprehend the formula for proportional loading, also known as one-parameter loading [29,30], the following is obtained: Here, F 0 denotes the constant initial applied force externally and m i is a scalar factor for load amplification that steadily rises.As m i is further elevated, the zones of plastic deformation of the body gradually expand, reaching a well-defined intensity, m p .At this point, unconstrained plastic deformation occurs, allowing for constant expansion in displacements and plastic deformations to be accomplished through external forces throughout the loading phase.Under a sustained external stress, an elasto-plastic body reaches its plastic limit state, wherein the entire body or specific sections undergo unrestricted deformation into a plastic form.The ultimate load at plastic deformation, F p , equals m p times the load at the elastic limit, F 0 .Equilibrium equations play a crucial role during the plastic limit phase, as forces and stresses can persist in static equilibrium.For instance, we demonstrate by contrasting the stress, σ ij , in a body experiencing quasi-static equilibrium at the plastic limit load F p = m p F 0 to the stress, σ s ij , within arbitrary static permissibility force and stress F is = m s F 0 , meeting the yield condition.
where the symbol k represents the plastic characteristics of the material.The concept of imaginary velocities can be extended to the force and stress fields by considering a body capable of deformation with volume V and surface of loading S q and also by incorporating kinematically permissible strain rate .
These two equations are subtracted to give the following: Each surface point needs to be normalized to ensure compliance with both the convexity of the yield surface and the normality criterion, resulting in all points falling within the normal distribution.
hence, Equation ( 16) yields: In this depiction, the integral symbolizes that external forces perform work on the actual body's velocities.The plastic limit condition avoids the occurrence of negative work, the condition m s − m p ≤ 0 must remain positive.The enhanced aspect of uniqueness and convexity in the procedure is demonstrated by the outlined steps.The incorporation of this element into topology design does not escalate the mathematical complexity of the method.

Experimental Work
An experimental program was undertaken to investigate the performance of two-way RC slabs strengthened with CFRP laminates.The investigation entails the construction of four two-way slab specimens supported at each end and subjected to concentrated load.All slab specimens had identical dimensions, 1200 mm × 1200 mm × 130 mm, as illustrated in Figure 2a.The specimens were categorized as GA, GB, or GC, with slab G0 serving as the control slab specimen.Furthermore, experimental measurements were conducted in the laboratory to assess concrete properties.The test specimens were produced using concrete that had a compressive strength of 30 MPa, as measured by its typical cylindrical strength.For the purpose of evaluating the compressive strength of concrete, a total of three cylinders measuring a diameter of 150 mm and a height of 300 mm and three cubes measuring side length of 150 mm were employed.These specimens were subjected to a 28-day curing period, after which their qualities in both tension and compression were examined, resulting in average concrete values.Table 1 displays the compressive strength results of concrete for both cylinders and cubes.The experimental measurements conducted in the laboratory encompassed evaluations of both concrete compressive behavior and tensile behavior, as illustrated in Figure 3.
In this depiction, the integral symbolizes that external forces perform work on the actual body's velocities.The plastic limit condition avoids the occurrence of negative work, the condition   −   ≤ 0 must remain positive.The enhanced aspect of uniqueness and convexity in the procedure is demonstrated by the outlined steps.The incorporation of this element into topology design does not escalate the mathematical complexity of the method.

Experimental Work
An experimental program was undertaken to investigate the performance of twoway RC slabs strengthened with CFRP laminates.The investigation entails the construction of four two-way slab specimens supported at each end and subjected to concentrated load.All slab specimens had identical dimensions, 1200 × 1200 × 130 mm, as illustrated in Figure 2a.The specimens were categorized as GA, GB, or GC, with slab G0 serving as the control slab specimen.Furthermore, experimental measurements were conducted in the laboratory to assess concrete properties.The test specimens were produced using concrete that had a compressive strength of 30 MPa, as measured by its typical cylindrical strength.For the purpose of evaluating the compressive strength of concrete, a total of three cylinders measuring 150 × 300 mm and three cubes measuring 150 mm were employed.These specimens were subjected to a 28-day curing period, after which their qualities in both tension and compression were examined, resulting in average concrete values.Table 1 displays the compressive strength results of concrete for both cylinders and cubes.The experimental measurements conducted in the laboratory encompassed evaluations of both concrete compressive behavior and tensile behavior, as illustrated in Figure 3.In all specimens, the primary tensile reinforcement is composed of fifteen steel bars, each with 10 mm diameter, positioned with 61 mm spacing center-to-center and the compression reinforcement consisted of twelve steel bars, each with 6 mm diameter, positioned with 85 mm spacing center-to-center.The characteristics of the reinforcements are detailed in Table 2.The parametric study was designed in relation to NSM CFRP laminates strengthening.This procedure involves affixing CFRP laminate strips to the underside of the slab using an adhesive substance known as epoxy resin.The mechanical properties of CFRP and epoxy are displayed in Table 3.The designation of this group is as follows: Slab GA was reinforced with five strips of CFRP laminates in both directions with 1.2 mm thickness and 15 mm width, as shown in Figure 2b.Slab GB was reinforced with five strips of CFRP laminates in both directions with 1.2 mm thickness and 30 mm width.Slab GC was reinforced with five strips of CFRP laminates in both directions with 1.2 mm thickness and 45 mm width.

Test Procedure and the Response of Concrete Slabs
The structural performance of the two-way RC slab was monitored at each loading stage, using instrumentation.Data on various parameters including central deflection, fracture widths, initial crack load, and eventual failure load were collected during each loading increment.The testing continued until failure was observed; the maximum burden was recorded.Using a load cell with a capacity of 700 kN, deflection measurements were taken at the center of the bottom tension face of the slab specimens.A dial gauge graduated to 0.01 mm increments was employed for the load measurements.
The specimens undergoing examination were meticulously positioned within the specialized frame apparatus, disregarding their own weight.In the initial phases of the experiment, a load of 5 kN was progressively applied.At zero loading, initial dial gauge readings were obtained.As the applied force escalated by 5 kN at each phase, cracks form on the concrete slabs' surfaces, documented using a pencil.Concurrently, the initial cracking load was recorded at each load increment, fracture breadth was assessed utilizing steel fillers, and deflection data were gathered via the dial gauge.Following the completion of this procedure, loading continued with the subsequent load step, in accordance with the same protocol.The burden continued to increase until it reached its maximum value, at which juncture it was documented.
Table 4 presents the load of cracking (Pcr), ultimate load (Pu), width of cracks (Wcr) at the ultimate load, and the central deflections (∆v) at the ultimate load.The control slab is designed to undergo flexural failure.For comparative analysis, the following was performed: The control slab specimens underwent testing without any strengthening technique.Its design aimed for flexural failure.The load was systematically increased and the slab's behavior was closely monitored.During the initial phase of testing, the slab exhibited elastic behavior until the emergence of the first crack.The first crack became evident in the constant moment region at a load of 120 kN.With continued loading, flexural cracks commenced in the specimens, extending towards the edges, and widening rapidly.At a load of 241 kN, slab G0 experienced flexural failure, marked by a sudden increase in deflection, while the load remained constant.
External CFRP laminates strengthening-in this procedure, CFRP laminate strips were affixed to the underside of the slab using epoxy resin as the adhesive material.The initiation of the initial crack occurred at a load of 190.6 kN for GA, 208.5 kN for GB, and 237.28 kN for GC.At loads of 455 kN, 498 kN, and 564 kN, GA, GB, and GC slabs experienced failure through a punching mode, respectively.Comparing to the control slab G0, there were increases in ultimate loads by 88.7%, 107%, and 134% for slabs GA, GB, and GC, emphasizing the notable influence of augmenting the reinforcement ratio within CFRP laminates.While the load remained constant, there was a rapid increase in deflection.For comparative analysis, the control slab was compared to these three specimens; Figure 4 illustrates the load deflection behavior for the control slab alongside the GA, GB, and GC, experimentally.Observably, it is apparent that, by increasing the applied load, there is a proportional rise in damage intensity, signifying that the additional load induces supplementary initial stresses in both concrete and steel components.This emphasizes the importance of complementary strain energy in representing plastic damage.Ultimately, the adept management of this factor alleviates instances of plastic damage and failure.Table 4 reveals that using CFRP laminate in strengthening the concrete slabs enhanced the ultimate load, in contrast to the control slab.Furthermore, slab GC demonstrated the maximum ultimate load and the smallest crack width.The data also indicate that slabs strengthened with CFRP laminates experienced punching failure, while slab GA failed due to debonding.These results underscore the effective capacity of CFRP in enhancing the ultimate load bearing capability of slabs, concurrently reducing the central deflection and the crack width under equivalent loads.

Numerical Modeling
FE analysis was employed to replicate the nonlinear response exhibited by two-way RC slabs strengthened with CFRP laminates, as depicted in Figure 5.The damage plasticity model was adopted in this study, to represent concrete inside the Abaqus CAE ® .The modeling process involved four simply supported two-way RC slabs, derived from the experimental outcomes of one specimen-tested slab as a control slab, without using CFRP, as well as three strengthened concrete slabs with CFRP in different ratios.The concrete slab was modeled using the 8-node solid (C3D8R) element; as for the reinforcement bars, they were depicted as 2-node linear beam elements (T3D2).An embedded region was employed to simulate the bond between concrete and reinforcements.Additionally, the bond between concrete and CFRP was modeled as a surface-to-surface contact, as specified in the Abaqus CAE ® manual.Steel-bearing plates (250 × 250 × 25 mm) were strategically placed to prevent local failure caused by crushing, utilizing the (C3D8R) 8-node solid element for modeling.Boundary conditions were selected to include a pin support at the first end for rotation and a roller support at the other end for rotation and horizontal movement.Additionally, at the center of the slab a concentrated vertical load was applied to simulate the experimental conditions.To ensure accuracy, a fine mesh was applied, resulting in a total element count for the slabs approximately equal to 21,600 elements.The impact of varying mesh sizes on both numerical precision and computational efficiency was analyzed, revealing that a mesh size of 20 yielded the highest level of accuracy.
For model validation, the slab was compared to a recent study with similar properties and analysis approaches.The validation process involved 3D finite element models, using Abaqus CAE ® and the CDP model.Experimental results from four simply supported twoway slabs were used for validation.Concrete material properties were determined through standard cube tests for compressive behavior and split-cylinder tests for tensile behavior, revealing an average compressive strength of f c` = 30 MPa.

Numerical Modeling
FE analysis was employed to replicate the nonlinear response exhibited by two-way RC slabs strengthened with CFRP laminates, as depicted in Figure 5.The damage plasticity model was adopted in this study, to represent concrete inside the Abaqus CAE ® .The modeling process involved four simply supported two-way RC slabs, derived from the experimental outcomes of one specimen-tested slab as a control slab, without using CFRP, as well as three strengthened concrete slabs with CFRP in different ratios.The concrete slab was modeled using the 8-node solid (C3D8R) element; as for the reinforcement bars, they were depicted as 2-node linear beam elements (T3D2).An embedded region was employed to simulate the bond between concrete and reinforcements.Additionally, the bond between concrete and CFRP was modeled as a surface-to-surface contact, as specified in the Abaqus CAE ® manual.Steel-bearing plates (250 mm × 250 mm × 25 mm) were strategically placed to prevent local failure caused by crushing, utilizing the (C3D8R) 8-node solid element for modeling.Boundary conditions were selected to include a pin support at the first end for rotation and a roller support at the other end for rotation and horizontal movement.Additionally, at the center of the slab a concentrated vertical load was applied to simulate the experimental conditions.To ensure accuracy, a fine mesh was applied, resulting in a total element count for the slabs approximately equal to 21,600 elements.The impact of varying mesh sizes on both numerical precision and computational efficiency was analyzed, revealing that a mesh size of 20 yielded the highest level of accuracy.
For model validation, the slab was compared to a recent study with similar properties and analysis approaches.The validation process involved 3D finite element models, using Abaqus CAE ® and the CDP model.Experimental results from four simply supported two-way slabs were used for validation.Concrete material properties were determined through standard cube tests for compressive behavior and split-cylinder tests for tensile behavior, revealing an average compressive strength of f c' = 30 MPa.
the tests of the mechanical properties on specimens.The integration of these traits within Abaqus CAE ® facilitates the determination of CDP parameters, ensuring an accurate representation of concrete's desired damage behavior.Following sensitivity analyses, the CDP plasticity parameters are determined using the information provided in Table 5.However, it is important to note that these CDP parameters stay consistent throughout the following optimization iterations.In Figure 5, the numerical model is illustrated.6 highlights the congruence in tension damage patterns showcasing a significant reduction in red damaged areas in the flexural portion of the slabs, with the use of CFRP; the Abaqus CAE ® calibrated numerical models exhibited strong agreement with experimental findings, affirming the model's efficacy in predicting the performance of two-way RC slabs reinforced with CFRPs.The comparable results produced by both experimental and numerical methods are demonstrated in Figure 7.Both compressive crushing and tensile cracking are acknowledged as forms of failure mechanisms within the dataset, pertaining to concrete damage plasticity, derived from the tests of the mechanical properties on specimens.The integration of these traits within Abaqus CAE ® facilitates the determination of CDP parameters, ensuring an accurate representation of concrete's desired damage behavior.Following sensitivity analyses, the CDP plasticity parameters are determined using the information provided in Table 5.However, it is important to note that these CDP parameters stay consistent throughout the following optimization iterations.In Figure 5, the numerical model is illustrated.Figure 6 highlights the congruence in tension damage patterns showcasing a significant reduction in red damaged areas in the flexural portion of the slabs, with the use of CFRP; the Abaqus CAE ® calibrated numerical models exhibited strong agreement with experimental findings, affirming the model's efficacy in predicting the performance of two-way RC slabs reinforced with CFRPs.The comparable results produced by both experimental and numerical methods are demonstrated in Figure 7. Furthermore, Table 6 illustrates alternative findings, presenting the percentages of tension concrete damage (  %) observed in each case at different loads in specimen GC.The percentages are graphically depicted in Table 6, with red denoting areas of high damage intensity and blue indicating undamaged regions.Damage intensity rises with increasing load, suggesting that additional stress is introduced into the concrete, steel, and CFRP materials.This highlights the importance of incorporating complementary strain energy to accurately depict plastic damage.Overall, the effective control of this role mitigates plastic damage and failure occurrences.Furthermore, Table 6 illustrates alternative findings, presenting the percentages of tension concrete damage (d t %) observed in each case at different loads in specimen GC.The percentages are graphically depicted in Table 6, with red denoting areas of high damage intensity and blue indicating undamaged regions.Damage intensity rises with increasing load, suggesting that additional stress is introduced into the concrete, steel, and CFRP materials.This highlights the importance of incorporating complementary strain energy to accurately depict plastic damage.Overall, the effective control of this role mitigates plastic damage and failure occurrences.Furthermore, Table 6 illustrates alternative findings, presenting the percentages of tension concrete damage (  %) observed in each case at different loads in specimen GC.The percentages are graphically depicted in Table 6, with red denoting areas of high damage intensity and blue indicating undamaged regions.Damage intensity rises with increasing load, suggesting that additional stress is introduced into the concrete, steel, and CFRP materials.This highlights the importance of incorporating complementary strain energy to accurately depict plastic damage.Overall, the effective control of this role mitigates plastic damage and failure occurrences.Furthermore, Table 6 illustrates alternative findings, presenting the percentages of tension concrete damage (  %) observed in each case at different loads in specimen GC.The percentages are graphically depicted in Table 6, with red denoting areas of high damage intensity and blue indicating undamaged regions.Damage intensity rises with increasing load, suggesting that additional stress is introduced into the concrete, steel, and CFRP materials.This highlights the importance of incorporating complementary strain energy to accurately depict plastic damage.Overall, the effective control of this role mitigates plastic damage and failure occurrences.Furthermore, Table 6 illustrates alternative findings, presenting the percentages of tension concrete damage (  %) observed in each case at different loads in specimen GC.The percentages are graphically depicted in Table 6, with red denoting areas of high damage intensity and blue indicating undamaged regions.Damage intensity rises with increasing load, suggesting that additional stress is introduced into the concrete, steel, and CFRP materials.This highlights the importance of incorporating complementary strain energy to accurately depict plastic damage.Overall, the effective control of this role mitigates plastic damage and failure occurrences.Furthermore, Table 6 illustrates alternative findings, presenting the percentages of tension concrete damage (  %) observed in each case at different loads in specimen GC.The percentages are graphically depicted in Table 6, with red denoting areas of high damage intensity and blue indicating undamaged regions.Damage intensity rises with increasing load, suggesting that additional stress is introduced into the concrete, steel, and CFRP materials.This highlights the importance of incorporating complementary strain energy to accurately depict plastic damage.Overall, the effective control of this role mitigates plastic damage and failure occurrences.

Results and Conclusions
This research investigated the elasto-plastic numerical modeling of reinforced concrete slabs strengthened with CFRP laminates.In order to accomplish this objective, the numerical model was validated by utilizing experimental data and the concrete damage plasticity (CDP) constitutive model in Abaqus CAE ® , which was employed to depict the behavior of the concrete.Consequently, it is feasible to draw the subsequent conclusions regarding this endeavor.

•
The performance of the reinforced slabs was distinctly influenced by the strengthening ratios.It was observed that an elevation in the strengthening ratio consistently led to a higher ultimate strength.

•
The enhancement of reinforcement ratios across all the three specimens led to a rise in initial cracking, alongside a simultaneous decrease in deflection.

•
Strengthening through CFRP laminate exhibited superior ultimate strength, cracking load, and minimized crack width, when compared to the control slab.

•
The damage pattern presented in the data reflects the percentages of concrete damage due to tension.It can be deduced that, as the applied load rises, the severity of damage also increases.This suggests that the additional load imparts additional initial stresses to both the steel and concrete components.The function of complementary strain energy is emphasized by reflecting the plastic damage and, by controlling it, the plastic damage and failure may also be managed.

Results and Conclusions
This research investigated the elasto-plastic numerical modeling of reinforced concrete slabs strengthened with CFRP laminates.In order to accomplish this objective, the numerical model was validated by utilizing experimental data and the concrete damage plasticity (CDP) constitutive model in Abaqus CAE ® , which was employed to depict the behavior of the concrete.Consequently, it is feasible to draw the subsequent conclusions regarding this endeavor.

•
The performance of the reinforced slabs was distinctly influenced by the strengthening ratios.It was observed that an elevation in the strengthening ratio consistently led to a higher ultimate strength.

•
The enhancement of reinforcement ratios across all the three specimens led to a rise in initial cracking, alongside a simultaneous decrease in deflection.

•
Strengthening through CFRP laminate exhibited superior ultimate strength, cracking load, and minimized crack width, when compared to the control slab.

•
The damage pattern presented in the data reflects the percentages of concrete damage due to tension.It can be deduced that, as the applied load rises, the severity of damage also increases.This suggests that the additional load imparts additional initial stresses to both the steel and concrete components.The function of complementary strain energy is emphasized by reflecting the plastic damage and, by controlling it, the plastic damage and failure may also be managed.

Results and Conclusions
This research investigated the elasto-plastic numerical modeling of reinforced concrete slabs strengthened with CFRP laminates.In order to accomplish this objective, the numerical model was validated by utilizing experimental data and the concrete damage plasticity (CDP) constitutive model in Abaqus CAE ® , which was employed to depict the behavior of the concrete.Consequently, it is feasible to draw the subsequent conclusions regarding this endeavor.

•
The performance of the reinforced slabs was distinctly influenced by the strengthening ratios.It was observed that an elevation in the strengthening ratio consistently led to a higher ultimate strength.

•
The enhancement of reinforcement ratios across all the three specimens led to a rise in initial cracking, alongside a simultaneous decrease in deflection.

•
Strengthening through CFRP laminate exhibited superior ultimate strength, cracking load, and minimized crack width, when compared to the control slab.

•
The damage pattern presented in the data reflects the percentages of concrete damage due to tension.It can be deduced that, as the applied load rises, the severity of damage also increases.This suggests that the additional load imparts additional initial stresses to both the steel and concrete components.The function of complementary strain energy is emphasized by reflecting the plastic damage and, by controlling it, the plastic damage and failure may also be managed.

Results and Conclusions
This research investigated the elasto-plastic numerical modeling of reinforced concrete slabs strengthened with CFRP laminates.In order to accomplish this objective, the numerical model was validated by utilizing experimental data and the concrete damage plasticity (CDP) constitutive model in Abaqus CAE ® , which was employed to depict the behavior of the concrete.Consequently, it is feasible to draw the subsequent conclusions regarding this endeavor.

•
The performance of the reinforced slabs was distinctly influenced by the strengthening ratios.It was observed that an elevation in the strengthening ratio consistently led to a higher ultimate strength.

•
The enhancement of reinforcement ratios across all the three specimens led to a rise in initial cracking, alongside a simultaneous decrease in deflection.

•
Strengthening through CFRP laminate exhibited superior ultimate strength, cracking load, and minimized crack width, when compared to the control slab.

•
The damage pattern presented in the data reflects the percentages of concrete damage due to tension.It can be deduced that, as the applied load rises, the severity of damage also increases.This suggests that the additional load imparts additional initial stresses to both the steel and concrete components.The function of complementary strain energy is emphasized by reflecting the plastic damage and, by controlling it, the plastic damage and failure may also be managed.

Results and Conclusions
This research investigated the elasto-plastic numerical modeling of reinforced concrete slabs strengthened with CFRP laminates.In order to accomplish this objective, the numerical model was validated by utilizing experimental data and the concrete damage plasticity (CDP) constitutive model in Abaqus CAE ® , which was employed to depict the behavior of the concrete.Consequently, it is feasible to draw the subsequent conclusions regarding this endeavor.

•
The performance of the reinforced slabs was distinctly influenced by the strengthening ratios.It was observed that an elevation in the strengthening ratio consistently led to a higher ultimate strength.

•
The enhancement of reinforcement ratios across all the three specimens led to a rise in initial cracking, alongside a simultaneous decrease in deflection.

•
Strengthening through CFRP laminate exhibited superior ultimate strength, cracking load, and minimized crack width, when compared to the control slab.

•
The damage pattern presented in the data reflects the percentages of concrete damage due to tension.It can be deduced that, as the applied load rises, the severity of damage also increases.This suggests that the additional load imparts additional initial stresses to both the steel and concrete components.The function of complementary strain energy is emphasized by reflecting the plastic damage and, by controlling it, the plastic damage and failure may also be managed.

Results and Conclusions
This research investigated the elasto-plastic numerical modeling of reinforced concrete slabs strengthened with CFRP laminates.In order to accomplish this objective, the numerical model was validated by utilizing experimental data and the concrete damage plasticity (CDP) constitutive model in Abaqus CAE ® , which was employed to depict the behavior of the concrete.Consequently, it is feasible to draw the subsequent conclusions regarding this endeavor.

•
The performance of the reinforced slabs was distinctly influenced by the strengthening ratios.It was observed that an elevation in the strengthening ratio consistently led to a higher ultimate strength.

•
The enhancement of reinforcement ratios across all the three specimens led to a rise in initial cracking, alongside a simultaneous decrease in deflection.

•
Strengthening through CFRP laminate exhibited superior ultimate strength, cracking load, and minimized crack width, when compared to the control slab.

•
The damage pattern presented in the data reflects the percentages of concrete damage due to tension.It can be deduced that, as the applied load rises, the severity of damage also increases.This suggests that the additional load imparts additional initial stresses to both the steel and concrete components.The function of complementary strain energy is emphasized by reflecting the plastic damage and, by controlling it, the plastic damage and failure may also be managed.

Results and Conclusions
This research investigated the elasto-plastic numerical modeling of reinforced concrete slabs strengthened with CFRP laminates.In order to accomplish this objective, the numerical model was validated by utilizing experimental data and the concrete damage plasticity (CDP) constitutive model in Abaqus CAE ® , which was employed to depict the behavior of the concrete.Consequently, it is feasible to draw the subsequent conclusions regarding this endeavor.
• The performance of the reinforced slabs was distinctly influenced by the strengthening ratios.It was observed that an elevation in the strengthening ratio consistently led to a higher ultimate strength.• The enhancement of reinforcement ratios across all the three specimens led to a rise in initial cracking, alongside a simultaneous decrease in deflection.

•
Strengthening through CFRP laminate exhibited superior ultimate strength, cracking load, and minimized crack width, when compared to the control slab.• The damage pattern presented in the data reflects the percentages of concrete damage due to tension.It can be deduced that, as the applied load rises, the severity of damage also increases.This suggests that the additional load imparts additional initial stresses to both the steel and concrete components.The function of complementary strain energy is emphasized by reflecting the plastic damage and, by controlling it, the plastic damage and failure may also be managed.• The outcomes underscored the commendable efficacy of CFRPs in enhancing the ultimate load capacity of slabs, while simultaneously decreasing crack width and central deflection under equivalent loads.• Alternatively, based on the numerical findings, varying quantities of CFRP strips have been observed to influence the behavior of the slabs.A comparison between experimental and numerical data reveals that the presence of CFRPs results in an increase in the slabs' ultimate load capacity, serving as additional support to withstand shear stress and confine the concrete section, thereby inducing a more brittle and abrupt failure mode.

•
In general, a reduction in the number of CFRP strips tends to exacerbate the extent of red area propagation, suggesting potential damage to either the steel or concrete.Conversely, employing more CFRP strips tends to mitigate the expansion of these red areas, akin to the effect observed with lower reinforcement ratios.In such cases, structures endure elevated stresses under lighter loads, when lower strengthening ratios are applied, resulting in failure at correspondingly lower loads.

•
Linear behavior was observed in all specimens up to the first cracking load, with the load deflection curve at mid-span potentially approximated by a bilinear or trilinear curve.Funding: This research received no external funding.

Figure 1 .
Figure 1.Response of concrete under the uniaxial loading conditions: (a) in compression, (b) in tension.

Figure 1 .
Figure 1.Response of concrete under the uniaxial loading conditions: (a) in compression, (b) in tension.

Figure 4 .
Figure 4. Experimentally load deflection relationships for all models.

Figure 4 .
Figure 4. Experimentally load deflection relationships for all models.

Figure 5 .
Figure 5. Details of the numerical model: (a) supporting and loading conditions, (b) details of the mesh, (c) strengthening details, and (d) details of reinforcement.

Figure 5 .
Figure 5. Details of the numerical model: (a) supporting and loading conditions, (b) details of the mesh, (c) strengthening details, and (d) details of reinforcement.

Figure 7 .
Figure 7. Experimental and numerical load deflection relationships for all four models.

Figure 7 .
Figure 7. Experimental and numerical load deflection relationships for all four models.

Figure 7 .
Figure 7. Experimental and numerical load deflection relationships for all four models.

Figure 7 .
Figure 7. Experimental and numerical load deflection relationships for all four models.

Figure 7 .
Figure 7. Experimental and numerical load deflection relationships for all four models.

Figure 7 .
Figure 7. Experimental and numerical load deflection relationships for all four models.

Table 1 .
Compressive strength of concrete.

Table 4 .
The experimental findings from the slab testing.

Table 4 .
The experimental findings from the slab testing.

Table 5 .
Inserted CDP data for concrete.

Table 5 .
Inserted CDP data for concrete.

Table 6 .
Comparison of damage patterns for slab GC.

Table 6 .
Comparison of damage patterns for slab GC.
i ) CDP Model Percentage of Damage (d t %) GC 1 Computation 2024, 12, x FOR PEER REVIEW 13 of 16

Table 6 .
Comparison of damage patterns for slab GC.

Table 6 .
Comparison of damage patterns for slab GC.

Table 6 .
Comparison of damage patterns for slab GC.

Table 6 .
Comparison of damage patterns for slab GC.