Abstract

Composite steel beams are commonly used element in multistorey steel buildings to enhance floor economy and serviceability and provide more clear height. Due to the low level of stress in the webs of such beams, hybrid sections are used where the flanges have higher strength than the webs. A lot of earlier research was carried out to optimize the design of the hybrid and nonhybrid composite steel beams under both static loading and dynamic behavior. However, there is still a need to develop a more practical optimization method. The aim of this research is to develop simple and practical equations to determine the optimum cross section dimensions for both shored and unshored, simply supported, hybrid and nonhybrid, composite steel beam under static loads. To achieve that goal, a research program of two phases was carried out. The first phase was generating a database of 504 composite beams with different steel grades for flanges and webs, subjected to different values of bending moment. The cross section of each beam in the database was optimized using GRG technique to minimize the cost considering the unit price of each steel grade. In the second phase, the generated database was divided into training and validation subsets and used to develop two predictive models using Nonlinear Regression (NLR) technique and Artificial Neural Network (ANN) technique to predict the optimum cross section dimensions and hence the optimum weight and cost. The accuracies of the developed models were measured in terms of average error percent. NLR and ANN models showed average error percent of 16% and 11%, respectively.

1. Introduction

Connecting steel beams to concrete decks to form composite floors is commonly used in multistorey steel buildings. The use of this system reduces both weight and cost of floors and increases their stiffness and severability. From construction point of view, composite floor could be shored or unshored. Shored floor is supported during concrete casting and till concrete hardening; accordingly, composite section supports all loads. On the other hand, unshored floor is casted without supporting, and hence weights of steel and concrete are supported by steel section only while the rest of loads are supported by the composite section. Steel part of the composite could be beam (hot rolled section, built-up section, or castellated beam), truss, or another girder type such as Vierendeel.

To achieve fully composite action, shear connectors between steel section and concrete deck should be strong enough to prevent any slippage. If slight slippage is permitted, the floor is called partially composite.

Optimizing the design of composite floors was intensively addressed in the last twenty years. Many researchers tried different optimization techniques such as Nonlinear Programming (NLP), Artificial Neural Network (ANN), Genetic Algorithm (GA), Ant Colony (AC), and Harmony Search (HS).

In this regard, Balaur S. Dhillon and Chen-Hsing Kuo [1] carried out the first trial to optimize the design of composite, hybrid plate girders using generalized geometric programming (GGP) technique. Samer Hendawi and Dan M. Frangopol [2] implement the corrosion effects in reliability-based optimization of composite hybrid plate girders using time-independent reliability-based structural optimization (RBSO) approach. Hojjat Adeli and Hongjin Kim [3] optimized the design of composite beams according to the LRFD of the AISC using the integer-discrete Nonlinear Programming and then mixed the results using neural dynamics model. Hongjin Kim and Hojjat Adeli [4] used floating-point Genetic Algorithm to optimize the LRFD design of composite floors according to AISC. S. Kravanja and S. Silih [5] used the well-known Nonlinear Programming (NLP) approach to optimize the design of both trusses and beams composite floors according to Eurocode. Christopher M. Foley and Warren K. Lucas [6] used Genetic Algorithm (GA) to optimize the design of both composite and noncomposite floors according to AISC considering vibration criteria. M. Veljkovic and B. Johansson [7] suggested some modifications in the design equations of Eurocode 3 to optimize the design of cylinder hybrid beams considering fatigue effect based on experimental and theoretical studies. Uroš Klanšek and Stojan Kravanja [8] considered the direct production costs including material costs, power consumption, and labor costs. Uroš Klanšek and Stojan Kravanja [9], in the 2nd part of their research, optimized the design of composite floors using Nonlinear Programming (NLP) approach. S. Kravanja and U. Klanšek [10] suggested an optimization technique for both truss and beam composite floors according to Eurocode using Nonlinear Programming (NLP). Ahmed B. Senouci and Mohammed S. Al-Ansari [11] presented a Genetic Algorithm (GA) model to optimize the design of composite beams according to the LRFD of the AISC considering the effect of span and loads. Kaveh and A. Shakouri Mahmud Abadi [12] optimized the design of composite floors according to the LRFD of the AISC using both traditional and improved Harmony Search algorithm. Kaveh and M. Ahangaran [13] continued their previous research using social Harmony Search technique. Victoria E. Roşca et al. [14] presented an optimization technique for composite beams design according to EN-1994-1-1/2006 using Nonlinear Programming (NLP) approach. Blachowski and W. Gutkowski [15] estimated the minimum weight of steel in composite floors to satisfy the design requirements under human induced vibrations. Roland Abspoel and Frans Bijlaard [16] suggested a new optimized design strategy for both hybrid and nonhybrid plate girder based on experimental and theoretical studies. Andrew J. Unander [17] studied the design of composite floors statistically based on a database containing 640 existing buildings to estimate the range of the optimum steel weight per unit area of floor. Hamid Eskandari and Tahereh Korouzhdeh [18] used Matlab software to generate a database of 20,000 composite beams to produce a contour map of cost with respect to affecting parameters. Rong He et al. [19] investigated the effect of optimizing the composite floors of super tall buildings on their lateral stiffness and serviceability. Pavol Juhás [20] studied the fatigue failure mechanisms of thin-walled hybrid plate girders. Tahereh Korouzhdeh et al. [21] used Improved Ant Colony (IAC) technique to optimize the design of composite beams according to the LRFD of the AISC, and the research results indicated that the efficiency of the improved technique surpassed the original one. R. Silva and T. A. Rodrigues [22] modeled the composite beam using FEM to generate a database and then applied sequential linear programming technique to optimize the design. N. M. Yossef and S. Taher [23] used Genetic Algorithm (GA) technique to optimize the design of castellated beams in composite floors with respect to section dimensions and opening size. Tahereh Korouzhdeh and Hamid Eskandari-Naddaf [24] suggested a cost-safety optimization technique based on standardized formulation using LRFD design method, and the results were verified and compared with those from AI techniques. Yong-Bo Shaoa et al. [25] presented a theoretical study for strength and behavior of laterally unrestrained S690 high-strength steel hybrid girders with corrugated webs. Ebid [26] developed an equation to estimate the optimum weight of nonhybrid unstiffened composite steel beams using the GRG optimization technique. On the other hand, Saghir et al. [27] carried out a comprehensive experimental study for the mechanical characterization of FRP composite pipes. Katsamakas et al. [28] presented a FEM-based model to study the behavior of SRG-strengthened R/C beams. Gohari et al. [29] suggested a new analytical solution for elastic flexure of thick multilayered composite hybrid plates resting on Winkler elastic foundation in air and water.

2. Objective

Despite the previous efforts in optimizing the design of hybrid composite floors beams, there are still some points that need further investigation such as the following:Consideration of unit weight cost of each type of steel partConsideration of construction sequence effect (shored or unshored)Presentation of optimization results in simple and practical equations

The objective of this research is to cover the above-mentioned points to end up with a set of equations to determine the optimum section dimensions at any location along the hybrid composite beam. This research is concerned with the weight optimization of the steel part excluding the shear connectors design.

3. Methodology

This research program is divided into two main phases. The first phase is concerned with generating a database of optimized composite beams, and the second phase is concerned with developing predictive models to predict the optimum cross section dimensions of the composite beam based on its configurations using the generated optimized dataset from the first phase. The detailed descriptions of each phase are provided in this section.

3.1. Phase 1: Generating the Optimized Database

The main difficulty of optimizing the hybrid composite beam is the large number of parameters that affect the design. These parameters could be classified as follows:Geometrical parameters: beams section dimensions, span, spacing, and type of supportsMaterial parameters: steel grads of web and flanges, slab concrete strength, and creep factorLoading parameters: values of dead, superimposed, and live loadsConstruction parameters: level of shoring (fully shored, partially shored, and unshored)Cost parameters: the unit cost of steel web and flanges besides the unit cost of the concrete

In order to organize and facilitate the research, the following assumptions are considered:To reduce the number of the considered parameters, beams span, spacing, type of support, and load value will be presented using one parameter which is the maximum total bending moment on the beam (M)The characteristic strength of concrete cube (Fcu) is limited between 25 and 45 MPa, and the allowable compressive stress in the concrete deck is 0.25 FcuThe ratio between elastic modulus of steel and concrete (modular ratio) (n) is calculated based on the concrete strength considering the creep factor equal to 2.0The considered steel grades have ultimate strength of 240, 360, 460, and 690 MPa

Factor a presents the ratio between the bending moment Mo of the initial stage (steel section only) and the total bending moment M (i.e., α = Mo/M). If α = 0, the initial bending moment is zero, and the beam is fully shored. On the other hand, if α = 1, the initial bending moment equals the total bending moment, and the beam is noncomposite:Factor η presents the ratio between the unit cost of flanges steel and the web steel (η = Cf/).The optimization will consider the unshored construction sequence as reference to optimize the shored condition.Effective width of concrete deck is considered as 12 times its thickness. Accordingly, beam span and spacing between beams shall not be less than 48 and 12 times deck thickness, respectively.Due to construction limitations, the minimum thickness of concrete deck is considered to be 10 cm.For steel section, the minimum compression flange area is considered as 20% of the tension flange area.The considered range of total bending moment M is between 1500 and 50000 cm.ton.The units considered in all calculations and equations are ton and cm.

Generally, the design of composite beam is divided into two phases: the first before concrete hardening and the second after concrete hardening. To calculate the stresses of each phase, properties of both steel and composite sections are calculated using an MS Excel sheet. Figure 1 illustrates the terms used in the calculations.

Calculation of section properties is based on the following assumptions:Compression and tension steel flanges are modeled as steel areas Asc and Ast, respectively, and act at their centers of gravityConcrete deck is modeled as equivalent steel area (Ac = 12 ts. Ts/2n) and acts in its center of gravity, where 2n is the modular ratio considering creep factor

The web plate is modeled as rectangular area with height equal to distance between upper and lower plates centers and thickness equal to the minimum value that keeps the section in noncompacted category (Hs/ = 190/√Fyw), where Fyw is the yield stress of web steel in (ton/cm2):Thicknesses of steel plates are neglected, and hence H ≈ Hs + tsSteel section cost = Cf. γs. Asc + Ast + Asw/η, where γs is the steel density

Since both Cf and γs are constants for considered case, the optimization issue is reduced to simple question, what is the combination of Ac, Asc, Ast, and Hs for a certain given bending moment that gives minimum value of Asc + Ast + Asw/η and maintains the stresses in the section below the allowable limits? To answer this question, both actual and allowable stresses must be calculated.

Since the section is classified as noncompacted section, there is no reduction in section properties due to local buckling, and hence the actual stresses for a certain phase equal the bending moment of this phase divided by the corresponding section elastic modulus. On the other hand, the allowable stress in concrete deck shall not exceed 0.25 Fcu. Regarding the steel section, considering that the upper compression flange is kept braced out of plan till concrete deck hardening, where there will be no reduction in the allowable normal stress due to lateral torsional buckling, the allowable normal stress for noncompacted, nonhybrid section is 0.58 fy. For noncompacted, hybrid sections, AASHTO Section 10.53.2 [1] suggested a reduction factor (R) to reduce the allowable stresses in the flanges. This reduction factor could be calculated as follows:where ψ is the distance from the neutral axis to the outer fiber of tension flange/steel section depth, β is the web area/tension flange area, and δ is the yield stress of web/yield stress of flanges.

In order to find the optimum combination of Ac, Asc, Ast, and Hs for a certain bending moment, an add-in tool to MS Excel software is used. The tool is called “Solver”, and it can figure out the combination of certain variables’ values that minimize or maximize the target function under certain governing conditions. It is based on a well-known mathematical technique called “Generalized Reduced Gradient (GRG) Nonlinear Solving.” This technique depends on changing the values of the affecting variables gradually while monitoring the governing conditions until the partial derivatives of the target function equal zero. In addition to the basic GRG technique, “Solver” tool enhances the searching performance by using automatic scaling for research intervals which adjusts the magnitude of change in each variable value for each iteration step.

Using the previously described “Solver” tool, the optimum combination of Ac, Asc, Ast, and Hs values is investigated for certain bending moment under the following governing conditions:Stress in upper plate shall not exceed 0.58 FyfR during the 1st phaseSummation of stress in upper plate shall not exceed 0.58 FyfR for both 1st and 2nd phasesSummation of stress in lower plate shall not exceed 0.58 FyfR for both 1st and 2nd phasesStress in concrete deck shall not exceed 0.25 Fcu during the 2nd phase

3.2. Phase 2: Developing the Predictive Models

A complete database of 504 records is generated using the previous discussed technique. The database was designed to cover the following variables’ ranges:Total bending moment (cm.ton) [M = 1500, 3000, 6000, 12500, 25000, 50000]Initial to total bending moment [α = Mo/M = 0.00, 0.25, 0.50, 1.00]Web and flange yield stress (t/cm2) [2.4/2.4, 2.4/3.6, 3.6/3.6, 3.6/4.6, 4.6/4.6, 4.6/6.9, 6.9/6.9]Flange steel to web steel unit cost [η = Cf/Cw = 1.00, 1.25, 1.50]

Each record contains the values of M, a, η, Fyw, Fyf, and Fcu and the corresponding optimized values of Hs, Asc, Asw, Ast, and ts. The complete database is presented in Table 1.

Finally, the generated database was used to find the correlations between inputs and outputs using both Nonlinear Regression (NLR) and Artificial Neural Network (ANN).

The generated database was divided into a training set (354 records) and validation set (150 records) as presented in Table 1. The training set is used to develop the predictive models while the validation set is used to test the efficiency (validity) of the developed models. Tables 2 and 3 summarize their statistical characteristics and the Pearson correlation matrix. Finally, Figure 2 shows the histograms for both inputs and outputs. The results of all developed models are summarized in Table 4.

4. Predicting the Optimized Cross Section

4.1. Model (1): Using NLR Technique

Studying the relation shapes between each input and output graphically indicated that power function is the suitable choice to describe the correlations with M, Fyw, and Fyf, while exponential function is the suitable choice to describe the correlations with a and η. Hence, the proposed general formula for these correlations is

The same add-in tool, “Solver,” was used to figure out the values of x2 to x6 for each output by minimizing the “Sum of Squared Errors” (SSE), while the value of x1 was used to adjust the function bias (the ratio between actual and calculated values). Applying the previously described methodology led to the following formulas:

This technique showed fair level of accuracy where error % and R2 values were 8.8%, 0.966; 27.4%, 0.941; 19.1%, 0.950; 18.1%, 0.950; and 5.6%, 0.991 for Hs, Asc, Asw, Ast, and ts, respectively, as shown in Figures 37. Figures 3(a) to 7(a) presents the correlations between optimized and calculated values.

4.2. Model (2): Using ANN Technique

Artificial Neural Network (ANN) is a type of soft computing that is fashioned to simulate the way the human brain analyzes and processes information. ANN consists of a number of nodes (neurons) arranged in layers. The first layer is called “input layer,” the last one is called “output layer,” and all intermediate layers are “hidden layers.” All neurons of a certain layer are connected to all neurons of the nest and previous layers by “links.” Each link has a weight representing the relation strength between the connected neurons, and each neuron has an “activation function” to convert the summation of inputs into outputs. The prediction process, “forward propagation,” starts by normalizing the input values (scaling them to lie within (0, 1) or (−1, 1) according to type of activation function), and then the normalized values are introduced to the neurons of the input layer. These values are transferred to the neurons of the first hidden layer through the links. The input of each neuron in the hidden layer is the summation of each connected neuron from the previous layer multiplied by the corresponding link weight, and the output of this neuron is the results of applying the activation function to the summation of inputs. This process propagates from hidden layer to the next till the output layer. The output of the output layer is the normalized values of the prediction; hence, these values are denormalized (scaled back to their original values) to get the ANN predicted values. Prediction accuracy depends on the learning ability of the ANN (adjusting the weights of the links). This process is carried out using a technique called “backward propagation” where the error (difference between predicted value and correct value) is distributed on the links of the previous layer based on the weight of each link, and then the weight of each link is adjusted according to its error share. This process propagates from the last hidden layer to the input layer.

A backward propagation ANN with two hidden layers and (Hypertan) activation function was used to predict the same output values. The used network layout and its connection weights are illustrated in Figure 8 and Table 5. Error % and R2 values were 6.7% and 0.979 for Hs in Figure 3, 16.0% and 0.982 for Asc in Figure 4, 12.2% and 0.980 for Asw in Figure 5, 11.9% and 0.977 for Ast in Figure 6, and 7.0% and 0.974 for ts in Figure 7, and the predicted values are presented in Figures 3(b) to 7(b).

5. Discussion

Studying (3)–(7) leads to the following notes:For noncomposite, nonhybrid steel section (Fyf = Fyw = Fy, α = 1, η = 1),This means that the optimum design for nonhybrid steel section neglecting the lateral torsional bulking is achieved when Asc = Ast = Asw.For fully shored composite, nonhybrid steel section (Fyf = Fyw = Fy, α = 0, η = 1),This means that the optimum design for fully shored composite, nonhybrid steel section is achieved when Ast = Asw = 3Asc. Moreover, it could be noted that the optimum depth and total steel weight are about 90% and 62% of those of the equivalent pure steel section, respectively.For unshored composite, nonhybrid steel section (Fyf = Fyw = Fy, α = 0.33, η = 1),

For unshored composite, nonhybrid steel section, where the initial load is about 33% of the total load (which is the common case in residential and office buildings), the optimum design is achieved when Ast = Asw = 2Asc. Furthermore, it could be noted that the optimum depth and total steel weight are about 95% and 73% of those of the equivalent pure steel section, respectively.

6. Conclusions

This research presents two predictive models to determine the optimum cross section dimensions of fully composite, unstiffened, built-up, hybrid steel girders. The research program was divided into two phases. The first is to generate a dataset of composite beams with optimized cross sections based on their steel grades, construction sequences, and material prices besides the bending moments. This phase was carried out using GRG technique and produced a dataset of 504 records. In the second phase, the generated dataset is divided into a training set used to develop the predictive models and a validation set used to test the validity of these models. Two techniques were used to develop the models, one mathematical technique (NLR) and one soft computing technique (ANN). The results of the present research revealed the following:(a)The average prediction accuracy of NLR and ANN models is about 92.8% and 93.1% for section dimensions (Hs, ts) and 78.5% and 86.5% for section weights (Asc, Asw, Ast), respectively.(b)Although, the ANN model is more accurate than NLR model, it is more complicated and cannot be manually used, which gives an advantage to NLR model which introduces a set of simple closed form equations.(c)The summation of the absolute weights of each neuron in the input layer of the ANN models indicates the importance of each input parameter. Accordingly, bending moment (M) is the most important parameter, and both Fyf and a shared the second ranking. Finally, Fyw, η, and Fcu are less important parameters.(d)The NLR developed equations are versatile and could be used to determine the optimum cross section dimensions for any design and construction parameter combinations as follows:Shoring condition, presented by the ratio of initial to total bending moment (α)Even pure steel girders (noncomposite), where α = 1Hybrid condition, presented by yield stresses of both flanges (Fyf) and web (Fyw)Economic condition, presented by the ratio of flange to web unit cost (η)(e)Testing the NLR developed equations for well-known common cases showed perfect matching as follows:For noncomposite, nonhybrid steel section (Fyf = Fyw = Fy, α = 1, η = 1), the optimum steel section cost occurs when each flange has the same area as webFor fully shored composite, nonhybrid steel section (Fyf = Fyw = Fy, α = 0, η = 1), the optimum steel section cost occurs when both tension flange and web areas equal about 45% of the total steel section area, and the area of compression flange area is about one-third of the area of tension flangeFor unshored composite, nonhybrid steel section (Fyf = Fyw = Fy, α = 0.33, η = 1), the optimum steel section cost occurs when both tension flange and web areas equal about 40% of the total steel section area, and the area of compression flange area is about half the area of tension flange(f)Like any other regression technique, the generated formulas are valid within the considered range of parameter values; beyond this range, the prediction accuracy should be verified.

Data Availability

The data supporting the outcome of this research work are included within the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.