Predicting Crack Width in CFRP-Strengthened RC One-Way Slabs Using Hybrid Grey Wolf Optimizer Neural Network Model

: This study deploys a hybrid Grey Wolf Optimizer Neural Network Model for predicting the crack width in reinforced concrete slabs strengthened with carbon ﬁber-reinforced polymers (CFRP). Reinforced concrete (RC) one-way slabs (1800 × 400 × 120 mm in size) were strengthened with CFRP with various lengths (1800, 1100, and 700 mm) and subjected to four-point bending. The experimental results were compared to corresponding values for conventional RC slabs. The observed crack width results were recorded, and subsequently examined against the expression recommended by Eurocode 2. To estimate the crack width of CFRP-reinforced slabs, ANN combined with the Grey Wolf Optimizer algorithm was employed whereby the applied load, CFRP width/length, X/Y crack positions, and stress in steel reinforcement and concrete were deﬁned as the input parameters. Experimental results showed that the larger the length and width of the carbon ﬁber, the smaller the maximum crack width in the tensile area of the slab at the ﬁnal load step. On average, the crack width in slabs retroﬁtted with CFRP laminates increased by around 80% compared to a slab without CFRP. The results conﬁrm that the equation provided by Eurocode 2 provides an unconservative estimation of crack widths for RC slabs strengthened with CFRP laminates. On the other hand, the results also conﬁrm that the proposed informational model could be used as a reliable tool for estimating the crack width in RC slabs. The ﬁndings provide valuable insight into the design approaches for RC slabs and rehabilitation strategies for existing deﬁcient RC slabs using CFRP.


Introduction
Reinforced concrete (RC) flat slabs are two-way reinforced horizontal structural elements that carry relatively light floor loads and transfer them to the structure's columns without the use of beams or girders. It is a favorable construction system offering several benefits, including (i) reducing the building height; (ii) flexibility in the design layout; (iii) decreasing the construction cost and time; and (iv) ease of installation of the electrical and mechanical networks.
Cracking of concrete generally occurs when the tensile strength of concrete is exceeded. This is inevitable in conventional RC flat slabs; once formed, the cracks will be pertinent to the slab's service lifetime. Since the cracks influence the serviceability of the slab, the limit state of excessive crack width needs to be considered in the design. Concrete may crack early in the load history. Most of the cracks in concrete initiate as a result of the following conditions: structure of the biological brain. They can establish a functional relationship between two datasets over the learning process and replicate that relationship over a recall process. An ANN comprises an adaptable methodology that can accurately estimate highly non-linear functions over the entire domain with adequate precision. Mangalathu et al. [17] studied several machine learning methodologies (support vector regression, ridge regression, random forest, decision tree, K-nearest neighbors, adaptive boosting, and extreme gradient boosting) for estimating the punching shear strength of flat slabs and identifying the most reliable tool. The study concluded that the extreme gradient boosting model provided the lowest mean square error and the highest coefficient of determination among the examined machine learning methodologies.
ANN predictive models have been applied to many problems in economics, engineering and other scientific fields. Particularly in structural engineering [18,19] and concrete structures, they have been successfully used for predicting the compressive strength of concrete containing recycled aggregate [20,21], predicting the compressive strength of cement replacement material (CRM) samples [22], determining the nominal shear capacity of steel fiber-reinforced concrete beams [23], predicting the capacity of concrete walls [24], predicting the compressive strength of lightweight foamed concrete [25], predicting the properties of FRP-confined concrete cylinders [26], designing reinforced concrete footings [27], analyzing the vulnerability of large concrete dams [28], measuring the mechanical properties of hybrid concrete through image processing [29], and predicting the creep and shrinkage deflection of reinforced concrete beams containing GGBFS [30], among other interesting and innovative applications [31,32].
Even though several design models have been proposed to estimate the crack width of RC flat slabs, such formulas are still perceived as a work-in-progress, especially for slabs reinforced with FRP laminates. While FRP laminates can benefit a RC flat slab in terms of enhanced sectional strength and long-term durability, the uncertainty regarding their mechanical properties remains a hurdle in guaranteeing the reliability of pertinent crack width equations. For instance, most available equations are either empirical, based on fitting the available test data, or are based on the modification of existing design code estimations for steel-reinforced RC slabs by considering the lower elastic modulus of FRP laminates. Nevertheless, the accuracy of the above-mentioned approaches is still questionable due to the differences in the mechanical properties and structural behavior between FRP and conventional steel reinforcement. As documented by several researchers, FRP laminates exhibit a fragile linear elastic response and different bond characteristics, resulting in different cracking behavior than conventional steel reinforcement [33][34][35].
The present study aims at developing a generalized informational model for determining the crack width of RC flat slabs reinforced with FRP laminates. Unlike existing empirical equations developed based on well-established fundamental mechanical theories using geometric and material properties, the proposed comprehensive informational model can inherently and effectively capture the underlying mechanisms of cracks in RC flat slabs. Using an informational model, the information about the underlying mechanisms of crack behavior is directly extracted from existing experimental datasets and is embedded in neural networks. This implies that the informational model does not need a pre-defined mathematical expression, in contrast with existing empirical equations. Considering the complexity of crack behavior, the proposed model can offer a reliable alternative approach. Its primary advantage lies in its ability to infer a general rule from the experimental dataset with greater efficiency than developing closed mathematical equations that, in some cases, may be impractical. To this end, by retrieving pertinent experimental results reported in the open literature, an experimental database was developed, which is representative of a large proportion of RC slabs (either with FRP or conventional steel reinforcement) under pure bending test conditions. Direct comparisons with existing theoretical models and international standards were performed using this extensive, custom-made dataset. Subsequently, an artificial neural network model was developed, coupled with a metaheuristic Grey Wolf Optimization algorithm (GWO-ANN) to form a generalized hybrid Buildings 2022, 12, 1870 4 of 26 model for estimating the crack width of RC flat slabs reinforced with FRP laminates. This model considers the key slab parameters and is designed specifically for RC flat slabs without shear reinforcement. The accuracy of the proposed informational model versus that of existing empirical equations was further examined using several statistical metrics to ensure the model's robustness and reliability.

Crack Width Calculation for Concrete Flat Slabs
The control of surface cracking in concrete is particularly important in certain situations. The most common of these is in cases where the surface will be visible, as excessive crack widths can give an overall impression of poor quality and limit the types of floor coverings that can be successfully applied on the slab. Crack control is also essential for durability, as the cracks will provide pathways for ingressive corrosive substances into the concrete, such as water or other chemicals.
The Eurocode 2 provides engineers with two methods for controlling cracking in reinforced concrete slabs. The overall thickness of the slab can influence the design procedure. The two methods are described as follows. Beeby and Narayanan [36] can be referred to for a more detailed account of the derivation of the crack width formulae in Eurocode 2.
where W k is the design crack width; β is a factor that relates the mean crack width in tests to the design value, e.g., it equals 1.7 for cracking due to direct loading; S crm is the average final crack spacing; and ε sm is the average difference in strain between the steel and the concrete, including the effects of bond stress, tension stiffening, concrete shrinkage, etc. The average final crack spacing S crm can be estimated using the following equation.
where K 1 is a factor that takes account of the bar bond properties (0.8 is specified in Eurocode 2); K 2 is a factor that takes account of the form of the stress distribution (0.5 for pure bending); and ρ r is the effective reinforcement ratio of the slab. The average strain, ε sm , is calculated at the section being considered, as where f s is the stress in the tension steel under the serviceability condition (calculated on the basis of a cracked section); f sr is the stress in the tension steel under the relevant condition that just causes the tensile strength of the concrete to be reached (calculated on the basis of a cracked section); β 1 is a factor that accounts for the bond properties of the reinforcement (1.0 for deformed bars); and β 2 is a factor that accounts for repeated stressing of the bars (0.5 for repeated stressing as the normal design situation). The minimum area of reinforcement placed at the cross-section being designed is given by Equation (4), as a basic requirement of the design rules of Eurocode 2 for controlling cracking without requiring a direct calculation.
where K 4 = z t A ct z ; f t is the mean value of the tensile strength of concrete at the critical time when the cracks might occur (a value of 3.0 MPa is recommended in Eurocode 2 for normal use); and A ct is the concrete area in the tensile zone at the section of concern before cracking.
Considering a simple rectangular slab, ignoring the presence of reinforcement when calculating z t and A ct , i.e., z t = bD 2 /6 and A ct = bD/2, and assuming z = 0.8D gives K 4 = 0.42. This explains the value of 0.4 given in Eurocode 2.
Meanwhile, the simplified design rules in Eurocode 2 are based on choosing an appropriate bar diameter or bar spacing. Design engineers usually prefer this approach rather than having to calculate the crack width directly. Limits are placed on bar diameters and bar spacing to ensure that the crack widths will not generally exceed 0.3 mm for reinforced concrete elements, as shown in Table 1.

Material and Methods
This section describes the characteristics of the materials utilized in this research work as well as the experimental program employed to collect data for developing the informational model.

Cement
ASTM C150 Type I Ordinary Portland Cement (OPC), produced from a single source with a specific surface of 3310 cm 2 /g and a relative density of 1.44 g/cm 3 , was utilized for all mixed materials. Table 2 shows the chemical properties of used OPC. Potable water, free from chemical contaminants, was used for both mixing and curing. It satisfies the requirements of BS 8110 (BS, 1997) [38].

Fine and Coarse Aggregates
The mining sand utilized in this study was sieved to a particle size range of 0.15 mm to 2.36 mm. The fine aggregate specific gravity in the saturated surface dry state was 2.61. In addition, crushed granite with a maximum size of 19 mm and a relative density of 2.65 was employed as the coarse material in this study. Figure 1 depicts the sand and coarse aggregate sieve analyses.

Water
Potable water, free from chemical contaminants, was used for both mixing and curing. It satisfies the requirements of BS 8110 (BS, 1997) [38].

Fine and Coarse Aggregates
The mining sand utilized in this study was sieved to a particle size range of 0.15 mm to 2.36 mm. The fine aggregate specific gravity in the saturated surface dry state was 2.61. In addition, crushed granite with a maximum size of 19 mm and a relative density of 2.65 was employed as the coarse material in this study. Figure 1 depicts the sand and coarse aggregate sieve analyses.

Carbon Fiber-Reinforced Polymer
CFRP is composed of carbon atoms and is encased in an extremely thin fiber with a diameter of 0.005-0.010 mm. CFRP has a high strength/weight ratio, good fatigue performance, and good electrochemical corrosion resistance, making it ideal for concrete structural applications [39]. It is commonly utilized in applications that require both strong mechanical characteristics and minimal weight. Sika's pultruded carbon fiber-reinforced polymer laminates were utilized in this study to enhance the flexural performance of flat concrete slabs. Table 3 presents the mechanical characteristics of the CFRP laminate provided by the manufacturer.

Carbon Fiber-Reinforced Polymer
CFRP is composed of carbon atoms and is encased in an extremely thin fiber with a diameter of 0.005-0.010 mm. CFRP has a high strength/weight ratio, good fatigue performance, and good electrochemical corrosion resistance, making it ideal for concrete structural applications [39]. It is commonly utilized in applications that require both strong mechanical characteristics and minimal weight. Sika's pultruded carbon fiber-reinforced polymer laminates were utilized in this study to enhance the flexural performance of flat concrete slabs. Table 3 presents the mechanical characteristics of the CFRP laminate provided by the manufacturer.

Adhesive
The Sikadur-30 adhesive from Sika was used to attach the CFRP to the RC slab. It is made of epoxy resins and a specific filler and is intended for usage at typical temperatures ranging from 8 • C to 35 • C. Table 4 shows the mechanical characteristics of the adhesive Sikadur-30 utilized in this study.

Concrete Mixture Design and Preparation
Concrete mix design calculates the proportions of the concrete components to obtain the required characteristics for performance and cost-efficiency. The design of experiments (DOE technique) provided by the British Department of the Environment [40] was utilized in the experimental mix design. The calculated concrete materials were combined, and the concrete specimens were cast and put under appropriate curing conditions. The tensile strength, compressive strength, and modulus of elasticity of all developed mixes were calculated based on the standards listed in Table 5. Table 6 summarizes the fresh and mechanical property test results for all tested slabs. Figure 2 shows the stress-strain curve of all tested specimens.   Figure 3 shows the wooden formwork constructed for the specimen casting. The internal dimensions of the formwork were 400 mm wide and 120 mm thick, while the length varied at 860, 1350, 1800, and 2400 mm. All the formwork was coated with form-release oil to achieve easier specimen removal after the concrete casting and curing.

Reinforcing Bar and Concrete Casting
The one-way RC slabs (either strengthened with CFRP or conventional non-strengthened ones) were reinforced longitudinally with 10 mm-diameter hot-rolled high-tensile strength deformed steel bars. The steel bars' modulus of elasticity and yield strength were 215 GPa and 610 MPa, respectively. Figure 4 depicts the longitudinal reinforcement of an RC slab.  Figure 3 shows the wooden formwork constructed for the specimen casting. The internal dimensions of the formwork were 400 mm wide and 120 mm thick, while the length varied at 860, 1350, 1800, and 2400 mm. All the formwork was coated with formrelease oil to achieve easier specimen removal after the concrete casting and curing.  Figure 3 shows the wooden formwork constructed for the specimen casting. The internal dimensions of the formwork were 400 mm wide and 120 mm thick, while the length varied at 860, 1350, 1800, and 2400 mm. All the formwork was coated with form-release oil to achieve easier specimen removal after the concrete casting and curing.

Reinforcing Bar and Concrete Casting
The one-way RC slabs (either strengthened with CFRP or conventional non-strengthened ones) were reinforced longitudinally with 10 mm-diameter hot-rolled high-tensile strength deformed steel bars. The steel bars' modulus of elasticity and yield strength were 215 GPa and 610 MPa, respectively. Figure 4 depicts the longitudinal reinforcement of an RC slab.

Reinforcing Bar and Concrete Casting
The one-way RC slabs (either strengthened with CFRP or conventional non-strengthened ones) were reinforced longitudinally with 10 mm-diameter hot-rolled high-tensile strength deformed steel bars. The steel bars' modulus of elasticity and yield strength were 215 GPa and 610 MPa, respectively. Figure 4 depicts the longitudinal reinforcement of an RC slab.  Before casting the concrete, the formwork's transparent cover and internal size were ensured by utilizing a 25 mm mortar block and verifying with a measuring tape, respectively. The slabs were cured in the formwork for three days after the concrete was poured. Subsequently, as shown in Figure 5, the slabs were wet gunny cured for seven days before being stored in an unregulated concrete laboratory until the testing day.

Instrumentation
The following testing equipment and instruments were deployed in the experimental program: • Linear variable differential transducers (LVDT) The deflection at the center of the slab was measured using an LVDT with a 50 mm maximum deflection. The LVDT was linked to a data logger, which recorded the midspan deflection as the load increased. •

Data Logger
The measurement mid-span deflection was recorded with a data logger (TS-TDS-302). •

Handheld microscope
A portable microscope was used to measure the flexural fractures at the level of the primary steel bar. The crack-measuring microscope had a precision of 0.02 mm and offers a 40× magnification.
Six RC one-way slabs with dimensions of 1800 × 400 × 120 mm and strengthened with various lengths and widths of CFRP were evaluated and compared to equivalent samples Before casting the concrete, the formwork's transparent cover and internal size were ensured by utilizing a 25 mm mortar block and verifying with a measuring tape, respectively. The slabs were cured in the formwork for three days after the concrete was poured. Subsequently, as shown in Figure 5, the slabs were wet gunny cured for seven days before being stored in an unregulated concrete laboratory until the testing day.  Before casting the concrete, the formwork's transparent cover and internal size were ensured by utilizing a 25 mm mortar block and verifying with a measuring tape, respectively. The slabs were cured in the formwork for three days after the concrete was poured. Subsequently, as shown in Figure 5, the slabs were wet gunny cured for seven days before being stored in an unregulated concrete laboratory until the testing day.

Instrumentation
The following testing equipment and instruments were deployed in the experimental program: • Linear variable differential transducers (LVDT) The deflection at the center of the slab was measured using an LVDT with a 50 mm maximum deflection. The LVDT was linked to a data logger, which recorded the midspan deflection as the load increased. •

Data Logger
The measurement mid-span deflection was recorded with a data logger (TS-TDS-302). •

Handheld microscope
A portable microscope was used to measure the flexural fractures at the level of the primary steel bar. The crack-measuring microscope had a precision of 0.02 mm and offers a 40× magnification.
Six RC one-way slabs with dimensions of 1800 × 400 × 120 mm and strengthened with various lengths and widths of CFRP were evaluated and compared to equivalent samples

Instrumentation
The following testing equipment and instruments were deployed in the experimental program: The deflection at the center of the slab was measured using an LVDT with a 50 mm maximum deflection. The LVDT was linked to a data logger, which recorded the mid-span deflection as the load increased.

• Data Logger
The measurement mid-span deflection was recorded with a data logger (TS-TDS-302).

• Handheld microscope
A portable microscope was used to measure the flexural fractures at the level of the primary steel bar. The crack-measuring microscope had a precision of 0.02 mm and offers a 40× magnification.
Six RC one-way slabs with dimensions of 1800 × 400 × 120 mm and strengthened with various lengths and widths of CFRP were evaluated and compared to equivalent samples without CFRP. All the slabs had two steel bars with a diameter of 10 mm as longitudinal reinforcement, and the concrete cover was 25 mm. Table 7 presents the characteristics of the specimens, i.e., unique slab code, CFRP width, and CFRP length. All of the prepared slabs were tested under four-point bending loading conditions. Figure 6 depicts the instrument used and the loading configuration. As shown in the figure, mechanical discs (pre-drilled stainless-steel discs) known as Demountable Mechanical (DEMEC) were glued to the concrete surface. A DEMEC strain gauge was used to measure the changes in length between the two DEMEC points.
Buildings 2022, 12, x FOR PEER REVIEW 10 of 26 without CFRP. All the slabs had two steel bars with a diameter of 10 mm as longitudinal reinforcement, and the concrete cover was 25 mm. Table 7 presents the characteristics of the specimens, i.e., unique slab code, CFRP width, and CFRP length. All of the prepared slabs were tested under four-point bending loading conditions. Figure 6 depicts the instrument used and the loading configuration. As shown in the figure, mechanical discs (pre-drilled stainless-steel discs) known as Demountable Mechanical (DEMEC) were glued to the concrete surface. A DEMEC strain gauge was used to measure the changes in length between the two DEMEC points.

Experimental Results
Six CFRP-reinforced one-way RC slabs were designed, manufactured, and tested in flexure, as shown in Figure 7. All of the prepared slabs had the same shape, dimensions (width, length, depth), and internal steel reinforcement, with the width and length of the

Experimental Results
Six CFRP-reinforced one-way RC slabs were designed, manufactured, and tested in flexure, as shown in Figure 7. All of the prepared slabs had the same shape, dimensions (width, length, depth), and internal steel reinforcement, with the width and length of the CFRP reinforcement serving as the distinguishing characteristic for each one of them. Figure 8 depicts the fracture patterns of the tested slabs. The results confirm that the steel yielding occurred prior to the collapse of the concrete in the compression zone in the slab failure mode.   As shown in Figure 9, before yielding of the steel reinforcement, the CFRP plate was de-bonded at the CFRP-concrete interface. Figure 10 shows the meticulous work of an   As shown in Figure 9, before yielding of the steel reinforcement, the CFRP plate was de-bonded at the CFRP-concrete interface. Figure 10 shows the meticulous work of an As shown in Figure 9, before yielding of the steel reinforcement, the CFRP plate was de-bonded at the CFRP-concrete interface. Figure 10 shows the meticulous work of an engineer who used a portable microscope to study and quantify the flexural fractures in the CFRP-enhanced RC slabs. This was undertaken carefully at each stage of the applied stress. similar benchmark slabs (in terms of materials used and dimensions), but without CFRP. For this purpose, Table 8 summarizes the bending moment at the first crack and the ultimate stage and crack width at the ultimate stage. The EC2 predicted the crack width at the service load. The results confirm the test crack width characteristics presented were reached at the post yielding load level. This behavior can be explained by the fact that attaching the CFRP laminates to the RC slabs leads to increased stiffness and bending moment capacity for the section.

Slab
First Crack Ultimate Stage The structural behavior of the RC slabs strengthened with CFRP was compared to similar benchmark slabs (in terms of materials used and dimensions), but without CFRP. For this purpose, Table 8 summarizes the bending moment at the first crack and the ultimate stage and crack width at the ultimate stage. The EC2 predicted the crack width at the service load. The results confirm the test crack width characteristics presented were reached at the post yielding load level. This behavior can be explained by the fact that attaching the CFRP laminates to the RC slabs leads to increased stiffness and bending moment capacity for the section.

Slab
First Crack Ultimate Stage The structural behavior of the RC slabs strengthened with CFRP was compared to similar benchmark slabs (in terms of materials used and dimensions), but without CFRP. For this purpose, Table 8 summarizes the bending moment at the first crack and the ultimate stage and crack width at the ultimate stage. The EC2 predicted the crack width at the service load. The results confirm the test crack width characteristics presented were reached at the post yielding load level. This behavior can be explained by the fact that attaching the CFRP laminates to the RC slabs leads to increased stiffness and bending moment capacity for the section. The analytical first crack moment is calculated based on the following equation: where f cr = 0.4 f c (in MPa units); I t is the moment of inertia of the transformed reinforced section (un-cracked); and y t is the distance from the center of the un-cracked transformed section to the extreme tension fiber. Furthermore, the position and width of the fractures formed by constant moment were examined and measured over the loading protocol. Figure 11 shows the measurement reference for recording the cracks, determined using the left-hand coordinate. Tables 9-11 list the recorded crack widths and locations for RC slabs retrofitted with CFRP laminates (all six specimens) during the pure bending testing. forced section (un-cracked); and yt is the distance from the center of the un-cracked transformed section to the extreme tension fiber. Furthermore, the position and width of the fractures formed by constant moment were examined and measured over the loading protocol. Figure 11 shows the measurement reference for recording the cracks, determined using the left-hand coordinate. Tables 9-11 list the recorded crack widths and locations for RC slabs retrofitted with CFRP laminates (all six specimens) during the pure bending testing.

ANN and Grey Wolf Optimization Algorithm
ANN is a computational model that simulates the practical features of biological neural networks [41]. As a parallel structure, an ANN consists of simple processing units similar to the human brain's structure, referred to as artificial neurons. The artificial neuron comprises bias, weights, and an activation function as per Equation (6), where f is the activation function, W m is the weight matrix, X m is the input vector, b is the bias vector, and Y is the output.
Grey wolves (which belong to the Canidae family) are recognized as apex predators located at the top of the food chain. In general, they live in packs with a group size of 5-12 with a strict community dominant hierarchy, as shown in Figure 12. The leaders are either female or male, and these so-called alphas (a, or leading wolf) principally make decisions regarding where/when to sleep, when to wake, and other activities, such as hunting. Every decision made by the a is fully dictated to the whole pack, though some cases of social behavior have also been observed, in which a may follow the other pack members. When they are near to each other, all of the group members admit the a as a leader by keeping their tails down. Only the a is allowed to mate in the pack. However, the a is not necessarily the strongest member, but rather the most skilled one in managing the group. In other words, the organization's discipline of the pack is of higher significance than the power of a group member. comprises bias, weights, and an activation function as per Equation (6), where is the activation function, is the weight matrix, is the input vector, is the bias vector, and is the output.
Grey wolves (which belong to the Canidae family) are recognized as apex predators located at the top of the food chain. In general, they live in packs with a group size of 5-12 with a strict community dominant hierarchy, as shown in Figure 12. The leaders are either female or male, and these so-called alphas (a, or leading wolf) principally make decisions regarding where/when to sleep, when to wake, and other activities, such as hunting. Every decision made by the a is fully dictated to the whole pack, though some cases of social behavior have also been observed, in which a may follow the other pack members. When they are near to each other, all of the group members admit the a as a leader by keeping their tails down. Only the a is allowed to mate in the pack. However, the a is not necessarily the strongest member, but rather the most skilled one in managing the group. In other words, the organization's discipline of the pack is of higher significance than the power of a group member.
In this hierarchy, the second level belongs to the beta, β, which is a subsidiary member that aids a in performing its activities and making decisions. The β can either be female or male and is possibly the most suitable candidate to be the a if one of the a wolves becomes very old or passes away. The β should obey the a, but commands the other wolves in the pack. β wolves play the role of a counselor to the a and a disciplinarian for the pack. The β reinforces the a's orders throughout the pack and advises the a.
The lowest position in the grey wolves' pack is the omega (ω). The ω is representative of a scapegoat, should always obey all of the other leading wolves, and comprises the wolves in the pack least allowed to eat. It may appear that the ω is not a group member; nevertheless, several observations have shown that the entire pack faces problems and internal fighting once losing the ω. It therefore helps to satisfy the whole pack and preserves the dominant structure. In the case when a wolf is neither α nor β or ω, it is a subordinate (or delta, δ). δs obey and , but are still dominant over . Sentinels, elders, hunters, caretakers, and scouts are categorized in this group. In general, scouts warn the pack when danger lies ahead and inspect the territory's boundaries, while sentinels defend the pack's safety. Elders are α β δ ω In this hierarchy, the second level belongs to the beta, β, which is a subsidiary member that aids a in performing its activities and making decisions. The β can either be female or male and is possibly the most suitable candidate to be the a if one of the a wolves becomes very old or passes away. The β should obey the a, but commands the other wolves in the pack. β wolves play the role of a counselor to the a and a disciplinarian for the pack. The β reinforces the a's orders throughout the pack and advises the a.
The lowest position in the grey wolves' pack is the omega (ω). The ω is representative of a scapegoat, should always obey all of the other leading wolves, and comprises the wolves in the pack least allowed to eat. It may appear that the ω is not a group member; nevertheless, several observations have shown that the entire pack faces problems and internal fighting once losing the ω. It therefore helps to satisfy the whole pack and preserves the dominant structure.
In the case when a wolf is neither α nor β or ω, it is a subordinate (or delta, δ). δs obey α and β, but are still dominant over ωs. Sentinels, elders, hunters, caretakers, and scouts are categorized in this group. In general, scouts warn the pack when danger lies ahead and inspect the territory's boundaries, while sentinels defend the pack's safety. Elders are skilled wolves that have formerly been an α or β. Hunters assist the α and β in hunting and preparing food. Finally, caretakers care for the sick or weak members of the group. Furthermore, group hunting is another interesting aspect of the social hierarchy of a wolf' pack. The fundamental steps involved in grey wolves' hunting process have been summarized by Muro et al. [42] as follows: (i) chasing, tracking, and pursuing the prey; (ii) surrounding and harassing the prey until it stops running; and (iii) attacking the prey. Figure 13 illustrates all of these steps. Further information about the meta-heuristic Grey Wolf Optimizer can be found in [43], while the algorithm's performance in structural optimization problems compared to other metaheuristic algorithms has been studied in [44]. surrounding and harassing the prey until it stops running; and (iii) attacking the prey. Figure 13 illustrates all of these steps. Further information about the meta-heuristic Grey Wolf Optimizer can be found in [43], while the algorithm's performance in structural optimization problems compared to other metaheuristic algorithms has been studied in [44]. The steps of the optimization methodology using the Grey Wolf algorithm are shown in Figure 14. The iterations continue to find constant value as the stopping criteria, which are defined by the user. The objective function is the crack width of the one-way slab. The steps of the optimization methodology using the Grey Wolf algorithm are shown in Figure 14. The iterations continue to find constant value as the stopping criteria, which are defined by the user. The objective function is the crack width of the one-way slab. Figure 13. Hunting strategy of grey wolves: (A) chasing, approaching, and pursuing the prey, (B-D) pursuing, harassing, and surrounding, (E) inactive situation and attack [42] (with permission).
The steps of the optimization methodology using the Grey Wolf algorithm are shown in Figure 14. The iterations continue to find constant value as the stopping criteria, which are defined by the user. The objective function is the crack width of the one-way slab.  Table 12 illustrates the properties of the dataset used in this study, which is given in more detail in Appendix A. As mentioned earlier, the independent input parameters in each datum include the load magnitude, the ratio of the CFRP length to the slab bay, the ratio of the CFRP width to the slab bay, the crack location in the x-direction, the crack location in the y-direction, the stress in the steel bar, and the stress in the concrete, which form a 7 × 1 matrix, while the dependent output parameter is the crack width, which forms a 1 × 1 matrix. The data, including crack width, the stress in concrete and steel reinforcement, were generated at a given load for each specimen. Recording of the data began with the appearance of the first crack and continued until the ultimate capacity of the slab was reached. Any algebraic relationship, either causal or not, between two causal variables is known as the correlation or dependency. On a large scale, correlation denotes the degree to which a two or more variables are linearly correlated. A correlation matrix is a table presenting the relationship coefficients amongst the input variables where there is a correlation among the two parameters in each cell of the table. To analyze the data at an advanced level, a correlation matrix aids in summarizing the relationships between the data. The present paper's correlation matrix for the input/output parameters is displayed in Figure 15. Considering each input parameter's domain and evading any divergence in the results, Equation 7 was used to normalize each parameter in the range of −1 to 1, where X n is the normalized value of the parameter, X max is its maximum value, and X min is its minimum value. X is the original (non-transformed) value of the variable. According to the number of input parameters, which is seven, Equation (8) resulted in 15 neurons within the hidden layers. Therefore, different architectures with two hidden layers and a maximum of 15 neurons in total were investigated. Overall, 30 different network topologies were evaluated. It was revealed that the network with a topology with 7-7-4-1 layers (two hidden layers) reached the lowest error values for RMSE, AAE, VAF, and the highest value of R 2 to estimate the crack width. Table 13 depicts the statistical metrics of selected topology which provided the most accurate results compared to the other two topologies. The ANN developed in this research was the Newff feedforward, where 70% of the dataset was considered for training, and the remaining 30% was used for testing the network. The Grey Wolf Optimizer (GWO) algorithm provided the least prediction error for the trained structure and optimized both the biases and weights of the ANN. Table 14 shows the parameters of the employed GWO algorithm.  Figure 15 confirms that the output parameter, the crack width, mainly correlated with the depth of the beam (y location, stress in concrete and steel bars, and loading magnitude. Moreover, Figure 16 shows the histogram and half violin diagram of the output parameter.   10

Multiple Linear Regression and Imperialist Competitive Algorithm Models
To examine the reliability of the proposed hybrid GWO-ANN model, a multiple linear regression (MLR) model was also developed in this research. In the MLR model, some independent parameters mainly influence the dependent variable, as per Equation (10), where y is a dependent or output parameter and x1, x2, … are independent input parameters. a1, a2, … are coefficients of the equation.
Equation (11) shows the most appropriate coefficients for the MLR model for estimating the studied specimens' crack width, and

Topology
Train Test The trial-and-error approach was adopted to develop the most effective architecture of the ANN model, which can best reproduce the features of the experimental dataset. In the present paper, an advanced technique was implemented for determining the number of neurons within the hidden layers based on Equation (8).
where N H signifies the number of neurons within the hidden layers, and N I stands for the number of input variables. The Levenberg-Marquardt training algorithm and hyperbolic tangent stimulation function were used in all developed networks. Moreover, the statistical indices, including the root mean squared error (RMSE), average absolute error (AAE), variance account factor (VAF), and correlation coefficient R, as expressed in Equations (9)- (12), were used to evaluate the performance of developed topologies.
According to the number of input parameters, which is seven, Equation (8) resulted in 15 neurons within the hidden layers. Therefore, different architectures with two hidden layers and a maximum of 15 neurons in total were investigated. Overall, 30 different network topologies were evaluated. It was revealed that the network with a topology with 7-7-4-1 layers (two hidden layers) reached the lowest error values for RMSE, AAE, VAF, and the highest value of R 2 to estimate the crack width. Table 13 depicts the statistical metrics of selected topology which provided the most accurate results compared to the other two topologies. The ANN developed in this research was the Newff feedforward, where 70% of the dataset was considered for training, and the remaining 30% was used for testing the network. The Grey Wolf Optimizer (GWO) algorithm provided the least prediction error for the trained structure and optimized both the biases and weights of the ANN. Table 14 shows the parameters of the employed GWO algorithm.

Multiple Linear Regression and Imperialist Competitive Algorithm Models
To examine the reliability of the proposed hybrid GWO-ANN model, a multiple linear regression (MLR) model was also developed in this research. In the MLR model, some independent parameters mainly influence the dependent variable, as per Equation (10), where y is a dependent or output parameter and x 1 , x 2 , . . . are independent input parameters. a 1 , a 2 , . . . are coefficients of the equation.
Equation (11) shows the most appropriate coefficients for the MLR model for estimating the studied specimens' crack width, and Table 15

Conclusions
This research work examines the effects of carbon fiber-reinforced polymer (CFRP) strengthening on the crack development in reinforced concrete slabs and develops an informational model using a hybrid Grey Wolf Optimizer-Artificial Neural Network algorithm (GWO-ANN) to estimate the crack width. Six reinforced concrete (RC) slabs strengthened with CFRP of various lengths were tested, and the results were compared to those of conventional RC slabs. Applied loading, CFRP width/length, X/Y crack positions, and steel reinforcement and concrete stress were defined as the input parameters for the development of the informational model. The main findings of the research work are summarized below: • Before steel reinforcement's yielding, the CFRP plate was de-bonded at the CFRP/concrete contact. • EC2 provides an unconservative estimation for the RC slabs' crack widths when CFRP laminates are attached to the slab for strengthening purposes. This behavior can be explained by the fact that attaching the CFRP laminates to the RC slabs leads to increased stiffness and bending moment capacity of the section, which is not accounted for by the EC2 formulas and may be associated with the unconservative estimation for the crack widths.

Conclusions
This research work examines the effects of carbon fiber-reinforced polymer (CFRP) strengthening on the crack development in reinforced concrete slabs and develops an informational model using a hybrid Grey Wolf Optimizer-Artificial Neural Network algorithm (GWO-ANN) to estimate the crack width. Six reinforced concrete (RC) slabs strengthened with CFRP of various lengths were tested, and the results were compared to those of conventional RC slabs. Applied loading, CFRP width/length, X/Y crack positions, and steel reinforcement and concrete stress were defined as the input parameters for the development of the informational model. The main findings of the research work are summarized below: • Before steel reinforcement's yielding, the CFRP plate was de-bonded at the CFRP/concrete contact. • EC2 provides an unconservative estimation for the RC slabs' crack widths when CFRP laminates are attached to the slab for strengthening purposes. This behavior can be explained by the fact that attaching the CFRP laminates to the RC slabs leads to increased stiffness and bending moment capacity of the section, which is not accounted for by the EC2 formulas and may be associated with the unconservative estimation for the crack widths.
• On average, the crack width in slabs retrofitted with CFRP laminates increased by around 80% compared to a slab without CFRP. Nevertheless, increasing the length and width of CFRP laminates had a minor effect on strength and crack development. The corresponding author can provide any other details of this study upon request.

Conflicts of Interest:
The authors declare no conflict of interest.