Flexural Strength Prediction of Welded Flange Plate Connections Based on Slenderness Ratios of Beam Elements Using ANN

,e present study was aimed to investigate the strength and ductility of welded flange plate (WFP) connections. Using the FE method, 52WFP connections with different beam overall depths and beam flanges/web slenderness ratios were analyzed. Fragility curves indicated that, for a WFP connection which is designed based on the seismic codes, its strength is of more concern than its ductility. In addition, limiting the width-to-thickness ratios of the beam flanges and web plates to 0.3 ���� E/fy 􏽱 and 2.45 ���� E/fy 􏽱 , respectively, may not always lead to the achievement of adequate connection’s strength and ductility. Proposed theoretical formulas and artificial neural network(ANN-) based models developed in this study were able to adequately predict the connection strength.


Introduction
Steel moment-resisting frames had extensive brittle fractures in their welded connections during the Northridge (1994) and Kobe (1995) earthquakes [1].Since then, modifications for design procedure of pre-Northridge connections have been introduced.E70T-4 type welding was changed to E70-TGK2 with smooth welding access holes and backing bars removed from the bottom beam flange [2][3][4].is type of connection is now known as post-Northridge connection.
ree methods are generally used to improve the strength and ductility of the post-Northridge connections: (a) strengthening the connection by adding additional elements; (b) reducing the beam section by cutting a portion of the beam flange [5] or the beam web [6][7][8]; and (c) adding slit dampers at the connection area as additional energy dissipation source [9].
Among strengthening methods, one of the most common methods to enhance the ductility of rigid connections was using plates (Figure 1). is connection is named by FEMA350 [10] as welded flange plate (WFP) connection.In this connection, the applied moment and the shear force are transferred to the column flange through the top and the bottom flange plates and a shear plate (also known as shear tab), respectively.e rectangular-shaped top and bottom flange plates have the same width and length.
Whittaker and Gilani [11] tested three WFP connections.Specimens were composed of W21 and W24 beams and I-shaped columns of Grade 50 steel with two same size rectangular-shaped flange plates.ey investigated the effect of contribution of the column panel zone (PZ) to the plastic rotation capacity of the joint.In these tests, the joint plastic rotation capacity varied between 0.028 and 0.045 radian.All specimens (except the specimen of weak PZ) experienced excessive strength degradation.SAC Joint Venture [12] tested five WFP connections.All specimens were fabricated from W30 × 99 Grade 50 beams and Grade 50 rectangularor trapezoidal-shaped reinforcing plates.
e key design variables were loading history, PZ strength and stiffness, and flange plate thickness and geometry.None of the five connections failed in a brittle manner, and the maximum beam plastic hinge rotations ranged between 2.3% and 3.9% radians.
ese specimens experienced excessive strength degradation at large rotations.In recent years, Ghobadi et al. [13] and Gholami et al. [14] tested few WFP connections.In these tests, the section used for beams was either IPE270 (beam depth 270 mm) or a built-up section with the overall depth equal to 380 mm.ese experimental tests indicated that WFP connections can reach to an acceptable level of strength and ductility.
With reference to the results summarized in the previous paragraph, it can be seen that using ange plates was effective to enhance the ductility of the pre-Northridge connections.However, it always was not the case for the connection strength.Experimental tests of the WFP connections indicated that some of these specimens experienced excessive strength degradation such that, at 4% total rotation, the measured exural resistance of the connection determined at the column face (M cf ) was below the 80 percent of the plastic moment of the connected beam (M pb ).For instance, the M cf /M pb ratios at four percent total rotation for specimens RC06 (beam section W30 × 99) and FUSD1 (beam section W21 × 50) of [12] and [11] were 0.65 and 0.74, respectively.While this parameter for specimens LF30 (beam depth 380 mm) and RC2 (beam section IPE270; beam depth 270 mm) of [14] and [13] was 1.1 and 1.25, respectively.Based on the seismic provision [15], from the strength point of view, a connection is acceptable if the ratio of M cf and M pb at 4% total rotation is greater than 0.8.Hence, it might be concluded that, for a WFP connection which is even designed based on the seismic provisions, its strength is of concern rather than its ductility.e main di erences between the sections used by these researchers were the beam overall depth and the slenderness ratios of the beam elements.In comparison to the sections used by Kim et al. [12], Gholami et al. [14], and Ghobadi et al. [13] used shallow beams of very compact ange and web.In other words, slenderness ratio of the beam elements can play a signi cant role to enhance the strength and ductility of the WFP connections.A detailed discussion about other e ective parameters can be found in [16].
Arti cial neural network (ANN) is a well-known branch of arti cial intelligence.ANNs formulate a mathematical model for a system in which no clear relationship is available between inputs and outputs.In recent years, researchers have successfully employed the mathematical models developed on the basis of ANN to predict the behavior of structural elements in di erent applications [17][18][19][20].
At the moment, many of steel moment frames, especially in Iran, have been constructed using WFP connection.As a result, it seems that having a reliable mathematical model which is able to adequately predict the behavior of WFP connections can be useful for either detailing the connected beam for new buildings or detecting the weak connections in existing buildings.Hence, the main aim of the present study was to estimate the strength of the WFP connections based on the e ective parameters.For this purpose, using the nite element method, numerous WFP connections of di erent beam overall depth and various beam ange and web slenderness ratios were created and analyzed.en, using the results of FE analyses, new theoretical and ANN-based empirical formulas were proposed to predict the strength of the WFP connections.

Finite Element Modeling
In order to verify the accuracy of nite element modeling, specimens LF30 of [14] and RC06 of [12] were modeled using ANSYS software [21].Specimen LF30 was composed of a shallow beam, a built-up I-380 × 200 × 8 × 12 mm beam, and a built-up B-400 × 400 × 20 × 20 (mm) box column.e bottom ange was a rectangular-shaped plate 250 × 300 × 15 (length 300 mm; width 250 mm).e width of the top ange plate at the column face level was 200 mm and at the end of the ange plate was 170 mm. e length and thickness of the top ange plate were 300 mm and 20 mm, respectively.In this specimen, the steel grade used for both beam and column was ASTM A36.Specimen RC06 was composed of a W30 × 99 beam and a W14 × 176 column.For this specimen, the material used for both beam and column were made of A572 Grade 50 steel [12], and the length of the column and the beam were 3450 mm and 3400 mm, respectively.e thickness of the continuity plates and the web doubler plate are 25 mm and 10 mm, respectively.e high fracture toughness weld metal used was E70.
In nite element models, all connection's parts were modeled using shell elements (SHELL181 in ANSYS program).Interaction between plates (e.g., interaction between the anges of beam and the top or the bottom ange plates) was modeled using contact elements.SHELL181 is a multilayer four-nodes shell element which has the ability to model plasticity, large de ection, and large strain phenomena.In this case, each element was separated into ve layers across the thickness.e number of layers was selected based on the nite element study carried out by Gilton and Uang [22].In order to determine the appropriate mesh density, a study on mesh sensitivity was carried out based on the recommendations given by ANSYS software.Results were then compared with the related experimental results of [12,14].
To perform material nonlinearity analyses, plasticity behavior was based on the von Mises yielding criteria and the associated ow rule.Isotropic hardening was assumed for the monotonic analysis, whereas kinematic hardening was assumed for the cyclic analysis as used by Mao et al. [23] and Ricles et al. [24].For base metals, both bilinear and 2 Advances in Civil Engineering multilinear material responses were investigated.However, the analytical results showed that the use of multilinear response agreed well with the experimental results.For weld metals, a multilinear material response based on the material property given by Mao et al. [23] and Ricles et al. [24] was used.Nonlinear geometric analyses were performed through a small strain, large displacement formulation.e monotonic analyses were conducted by applying a monotonic vertical displacement load to the beam tip until more than 4 percent total rotation at column web center was achieved and the load history recommended by the AISC [15] was utilized for cyclic analyses.When loads are applied only in the vertical direction, the out-of-plane deformations (normal to the web) may not occur.erefore, in order to ensure that buckling would occur due to instability of model, the imperfect model was used.In order to determine the imperfect model, rst the buckling mode shapes were computed in a separate buckling analysis and then were implemented to perturb the original perfect geometry of the model.
In order to verify the accuracy of the models, the experimental results of specimens RC06 and LF30 were compared to the analytical results obtained from the nite element models in terms of load-rotation and moment-rotation curves.Figures 2 and 3 show these comparisons for specimens RC06 and LF30, respectively.As Figures 2 and 3 show, the applied load-rotation and moment-rotation curves obtained from the nite element analyses are in good agreement with the experimental ones.Figure 4 also shows the plastic equivalent strain distribution for specimen RC06 at 5.5 percent total rotation.In Figure 4, high concentration of the plastic equivalent strain at the end of the ange plates of the beam is quite evident, and similar results were achieved during the experimental tests of Kim et al. [12].Furthermore, as Figure 4(b) shows, for the nite element model, the beam anges experienced a signi cant local buckling which is completely consistent with those that occurred during the experimental tests.

Failure Criterion
It should be noted that the fracture prediction is the most questionable part of a nite element study because it is inherently a complicated phenomenon and is dependent on many parameters such as weld and base metal properties, weld defects, notch e ects, weld quality, and weld toughness.In this study, it was assumed that the quali ed welders and fabricators are employed, and high fracture toughness weld metals are used (as these should be based on the AISC 360-10 [25], Speci cation for Structural Steel Buildings).Considering the fact that the locations of high level of strains have signi cant probability of premature fracture and by comparing the nite element results of specimens RC06 and LF30 with the experimental results presented by Kim et al. [12] and Gholami et al. [14], the following failure criterion was de ned.
e connection fracture may occur when the von Mises strains at the whole width of either the beam ange (at the end of ange plate length) or the ange plates (at the column  Experimental [14] Finite element monotonic loading

Advances in Civil Engineering
face level) exceed the strain associated with the ultimate strength of the materials.A similar failure criterion was also used by other researchers [26].Note that the connection failure can also be predicted using rupture index (RI) (e.g., [27,28]).is index is defined as the ratio between the plastic equivalent strain index and the ductile fracture strain.However, this method may be better when the fracture initiates from the beam web rather than the beam flange.Note that none of the mentioned methods are intended to predict exact rotation capacities for connections; rather, they provide a tool for comparing various models.

Sections Used to Design WFP Connections
Proposed design procedures to modify pre-and post-Northridge connections are usually based on three beam sections: W36 × 150, W30 × 99, and W24 × 68.erefore, all parametric studies were done on specimens SPE07 (beam: W36 × 150; column: W14 × 257), SPE05 (beam: W30 × 99; column: W14 × 176), and SPE03 (beam: W24 × 68; column: W14 × 120).ese specimen sizes were chosen since they may be good representatives of the conventional pre/post Northridge specimen sizes, large, medium, and small [3] and were also tested in phase 1 of SAC steel projects [29].e beam and column lengths of these specimens were similar to those that were used for specimen RC06.In addition to these specimens, specimen SPE01 which is similar to the specimen LF30 of [14] was also used.
is specimen was composed of a built-up I-380 × 200 × 8 × 12 (mm) beam and a built-up box column (B-400 × 400 × 20 × 20 mm). is specimen was used so that we could have a wider range of the beam depths, from a very deep beam (W36 × 150) to a very shallow beam (I-380 × 200 × 8 × 12).
All WFP connections were designed based on the AISC [15,25].To design the flange plates, the maximum applied moment at the column face (M cf ) was determined using (1).In this equation, Z b and f ye are the beam plastic section modulus and the expected yield stress of the material, respectively.S c is equal to the length of the flange plate plus ¼ of the beam overall depth, and V exp is the expected shear force at the plastic hinge location.C pr is a factor to account for the peak connection strength, including strain hardening, local restraint, and additional reinforcement.
In FEMA350 [10], the C pr factor is given by equation ( f y + f u )/2f y , where f u is the specified minimum tensile stress of the material.FEMA350 [10] proposes the use of value 1.2 for any case of modified connections except where otherwise noted in the individual connection design procedures.In the present study, the C pr factor used to design the WFP connections was 1.1 which is as same as the one that was used by Kim et al. [12]. (1)

Finite Element Models and Results
In general, in WFP connections, the thickness of the flange plates is larger than those of the beam flanges. is subject is much more pronounced when the flange plates are designed based on the seismic provision (i.e., ( 1) is used to detail the flange plates).In this case, an elastic region with a length equal to the length of the flange plates forms at the beginning of the beam length.It therefore causes the majority of nonlinear deformations to be concentrated at the plastic hinge area which is the beam flanges and web plates located after the tip of the flange plates (Figure 4).In this case, the behavior of a WFP connection is mainly controlled by the slenderness of the beam elements and further strengthening of the connection area by adding additional elements such as rib plates may not be effective.Here, when the beam elements are slender, an excessive local buckling will occur in the beam flange and web plates, and consequently significant strength degradation can be expected.On the other hand, when the beam elements are too compact, the beam elements will not experience any controlled local buckling.In this case, the main energy dissipation source (i.e., beam end) will be lost, and consequently the majority of stresses and strains concentrate at the column face level.is finally can lead to the premature fracture of the connection and to the excessive reduction in the connection's strength and ductility.Hence, by controlling the local buckling of the beam flange and web plates, higher connection strength and ductility can be expected.e beam flange and web local buckling can be controlled by limiting the width-to-thickness ratios of these  4 Advances in Civil Engineering elements.ese ratios for the beam flange and web can be expressed by equations b/t � b fb /(2•t fb ) and h/t � (d b -2k)/t wb , respectively.In these equations, b fb , t fb , t wb , d b , and k are the beam flange width, the beam flange thickness, the beam web thickness, the beam overall depth, and the distance from the outer face of the beam flange to the web toe of the fillet, respectively.Based on the seismic provision for structural steel buildings [15], the width-to-thickness ratios of the flange and web of a beam used in a special moment frame, must be limited to 0.3

􏽱
, respectively, where E s and f y are the modulus of elasticity and the yield stress of the material, respectively.In order to realize the effect of the parameters b/t and h/t on the strength and ductility of the WFP specimens, these ratios were defined as a function of . Hence, based on the AISC [15] and using these definitions, the values of parameters α f and α w for a beam used in a special moment frame must be at most equal to 0.3 and 2.45, respectively. ( Table 1 summarizes the effects of different values used for parameters α f , α w , and L b /d b on the strength and ductility of WFP specimens SPE07, SPE05, SPE03, and SPE01.In Table 1, L b is the beam length measured between the faces of the two adjacent columns.M cf /M pb is the normalized connection strength at four percent total rotation under cyclic loading and θ is the total rotation of the connection at the failure time.Figure 5 graphically compares the effect of parameters α f , α w , and L b /d b on the strength of some of WFP specimens SPE07, SPE05, SPE03, and SPE01.Figure 5 was drawn as such that, in each case, one parameter which was considered as the main parameter varied in its total range, while the other two parameters were kept constant or varied in a limited range.As Figure 5 shows, by reducing the values of parameters α f and α w and increasing the value of parameter L b /d b the connection strength increased.However, as the results of Table 1 indicate, the high value of strength is not always accompanied by the high value of ductility such that excessive decrease in the value of parameters α f and α w was detrimental and caused a significant reduction in the connection ductility.ese connections generally failed due to the fracture of flange plates at the column face level.Based on the FE results of this study and available experimental results of other researchers, following questions might be answered.However, more accurate answers need more numerical and experimental analyses.

First Question.
How reliable is a WFP connection? is question should be answered from both strength and ductility points of view.From ductility point of view, a connection is reliable if it is able to accommodate a story drift angle of at least 0.04 rad [15].Using the available test results, FEMA355D [30] has presented related relationships to predict the connection plastic rotational capacity in terms of the connection type and the overall depth of the connected beam.In this reference, using the experimental results of only 5 specimens tested by Kim et al. [12], the rotational capacity of the WFP connections has been estimated using the regression analysis.
e recommended rotation limit presented by this reference is 3% plastic rotation (or 4% total rotation, assuming the elastic rotation is 1%). is means that, from ductility point of view, this connection is one hundred percent reliable.However, this question might be more accurately answered using the fragility curves.Here, finite element and experimental results presented in Table 1 were used to find the cumulative distribution function (CDF) also known as fragility functions for the normalized strength (M cf /M pb ) and ductility (θ) of WFP connections.Fragility curves were obtained by assuming that the data has lognormal distribution.It has been reported by many of researchers [31][32][33] that the lognormal distribution provides relatively good fit to empirical cumulative distributions computed from experimental or numerical data.As the fragility curves shown in Figures 6(a) and 6(b) indicate, the probabilities of achieving adequate strength (M cf /M pb ≥ 0.8) and ductility (θ ≥ 4%) for WFP connections are 49% and 92%, respectively.is might be the reason of why this connection type has not been prequalified by the AISC358-10 [34].ese results also indicate that, for a WFP connection, its strength is of more concern rather than its ductility.

Second
Question.Are 0.3 and 2.45 the appropriate values for parameters α f and α w , respectively, to achieve both adequate connection's strength and ductility simultaneously when WFP connection is used in a special moment frame? is question can be answered by comparing the results of specimen numbers 9, 30, 43, and 52 presented in Table 1.As the results show, only specimen number 52 in which the beam depth is equal to 380 mm could achieve adequate connection strength (i.e., M cf /M pb ≥ 0.8).In other words, the limiting values for parameters α f and α w (i.e., 0.3 and 2.45, resp.) must be modified for beams deeper than 400 mm.Similar conclusion might be drawn by investigating the experimental results of specimens RC06 and FUSD-1 tested by Kim et al. [12] and Whittaker and Gilani [11], respectively.Specimen FUSD-1 was composed of W21 ×50 beam (α f � 0.27 and α w � 2.10).
e value of M cf /M pb at four percent total rotation for this connection was 0.74.Based on the result obtained from specimen RC06 (M cf /M pb � 0.65, α f � 0.33 and α w � 2.21), Kim et al. [12] recommended a reduction in the maximum width-to-thickness ratio of the beam web from 520/ , where the latter is corresponded to α w equal to 2.05.

ird Question.
What can be the maximum values for parameters α f and α w to achieve the minimum required connection's strength and ductility requested by AISC [15]?Based on the results presented in Table 1, it can be concluded that 0.3 and 2.2 are the maximum permissible values for parameters α f and α w , respectively, when the beam depth is less than 800 mm.  1, it might be concluded that 0.3 and 1.9 are the best values for parameters α f and α w , respectively, if the beam depth is less than 800 mm.While for deeper beams, 0.275 and 1.9 might be the most suitable values for parameters α f and α w , respectively.In this study, in order to predict the strength of a WFP connection two types of mathematical formulas are proposed.First type is based on the theoretical formulations relevant to the exural strength of beams with noncompact elements which is discussed here.e second type is on the basis of the arti cial neural network-based models which is discussed in Section 6.3.FE results indicated that the behavior of a WFP connection which is designed based on the seismic codes (i.e., (1) is used to detail the ange plates) is mainly controlled by the slenderness of the beam elements.Hence, it seems that equation F4-13 of AISC 360-10 [25] which is relevant to the exural strength of beams with noncompact elements can be used to estimate the amount of strength reduction of a WFP connection at four percent total rotation when its beam undergoes ange and web local buckling.By using this equation and with respect to the answer of question 4 and assuming that the ratio of elastic to plastic section modulus of a beam section (S b /Z b ) is 0.87, the normalized strength of a WFP connection at 4% total rotation (M/M P ) can be estimated using ( 4)- (6).Note that these equations were derived by changing (1) the value of parameter λ pf (limiting slenderness parameter for compact anges) from 0.38 E/f y to 0.3 E/f y for beam depth less than or equal to 800 mm; (2) the value of parameter λ pf from 0.38 E/f y to 0.275 E/f y for beam depth greater than 800 mm; (3) the value of parameter λ pw (limiting slenderness parameter for compact web) from 3.76 E/f y to 1.9 E/f y ; (4) λ pf from 0.38 E/f y to 0.3 E/f y for beam depth less than or equal to 800 mm; (5)   Advances in Civil Engineering (limiting slenderness parameter for noncompact ange) from 1 E/f y to 0.38 E/f y ; and (6) the value of parameter λ rw (limiting slenderness parameter for noncompact web) from 5.7 E/f y to 3.76 E/f y .ese changes were made mainly based on the answer of question 4. In order to investigate the accuracy of ( 4), for all specimens presented in Table 1, the value of M/M P was calculated using this equation and then results were compared to the related experimental and numerical results.Figure 7 compares the related absolute errors for these specimens.As Figure 7 indicates, the proposed equation was able to adequately predict the connection strength with an acceptable level of error.e average of absolute error was 12.5%.a ect the result more? is question and also the previous question might be adequately answered using the ANNs which are discussed in the following section.

Artificial Neural Network (ANN)
ANN is inspired in natural nervous systems.A natural nervous system is composed of several natural neurons (or nodes) where they receive signals through synapses (or connections) and process them to obtain an output.As Figure 8 shows, an ANN is composed of three layers: one input layer, one output layer, and several hidden layers.Layers are made up of a number of interconnected "nodes" which contain an "activation function."Hence, nodes can be seen as computational units, and connections determine the information ow between nodes.Patterns are presented to the network via the "input layer," which communicates to one or more "hidden layers"   where the actual processing is done via a system of weighted "connections."e higher a weight of a node is, the stronger the input.e hidden layers then link to an "output layer" where the answer is output as shown in Figure 8. e number of nodes in the input layer is equal to the number of independent variables which are believed to have an e ect on the response of the network.e number of nodes in the output layer is associated to the number of dependent variables (i.e., outputs) which the network is going to predict them.ere is no speci c role to determine the number of nodes in the hidden layers.In fact, they can be determined by a try and error method to achieve the highest performance of the network which corresponds to the achievement of the lowest value of error.

Hidden layers Connections
Most ANNs contain some form of "learning rule" which modi es the weights of the connections according to the input patterns that it is presented with.Although there are many di erent kinds of learning rules used by neural networks, this demonstration is concerned only with one: the delta rule.e delta rule is often utilized by the most common class of ANNs called "backpropagational neural networks" (BPNNs).Backpropagation is an abbreviation for the backwards propagation of error [35].With the delta rule, as with other types of backpropagation, "learning" is a supervised process that occurs with each cycle or "epoch" (i.e., each time the network is presented with a new input pattern) through a forward activation ow of outputs and the backwards error propagation of weight adjustments [36].More simply, when a neural network is initially presented with a pattern it makes a random "guess" as to what it might be.It then sees how far its answer was from the actual one and makes an appropriate adjustment to its connection weights. is process can be graphically shown in Figure 9.In this study, Levenberg-Marquardt (LM) algorithm [37] is used as the learning rule.Generally, the LM algorithm needs less time and epochs for convergence, and it performs more e ciently compared to other learning rules, which in turn makes it ideal for training of di erent networks [37].e LM algorithm is an approximation of Newton's method, and it is very e cient for training networks which have up to a few hundred weights.In many cases, the Marquardt algorithm was found to converge when other backpropagation techniques diverged [20].Once a neural network is "trained" to a satisfactory level, it can be used as an analytical tool on other data.In this case, the user no longer speci es any training runs and instead allows the network to work in forward propagation mode only to predict the output as follows.
At rst each input value is multiplied by the connected weight value and summed.
en, the result is added to a constant value.is constant value is named bias and is

Advances in Civil Engineering
only connected to nodes in the hidden and output layers.e input into a node can be expressed by (6).
In this equation, net j is the collection of information from each node, x i is the input value, and w ij is the connection weight value.en, the net j is put through an activation function to produce the output of the neuron.Activation functions which serve to introduce nonlinearity into neural networks include a linear manner as step and rampage functions or nonlinear like Gaussian, sigmoid, hyperbolic tangent functions, and others.e outputs using a logistic activation function can be expressed by (7).out j f net j .(7) In this equation, out j is the output of the jth node and f is the activation function.Figure 10 graphically shows this procedure for a simple ANN network consisting of the input layer consisted of several inputs (X 1 through X n ), one hidden layer made up of only one node and one output in the output layer.Equation ( 8) expresses the output of such ANN.
In this expression, W 1 denotes the input layer weight matrix, W 2 represents the hidden layer weight value, and b 1 and b 2 are the hidden and the output layer bias values, respectively.

ANN Structure for Predicting the Strength of WFP
Connections.In this study, using the nite element results presented in Section 5, one ANN-based model was developed for predicting the strength of the WFP connections.is model consisted of three layers, one input layer, one hidden layer, and one output layer with either LogSigmoid or TanSigmoid transfer functions in the hidden layer and identity (or linear) transfer functions in the output layer.e number of nodes in the input layer was equal to the number of input variables to the ANN model which was equal to the number of e ective parameters.
ese e ective parameters are those used in the parametric study presented in Section 5 and are α f , α w , and L b /d b .As a result, the input layer of the ANN-based model was made up of 3 nodes.
e out layer of the ANN-based model consisted of one node.e number of nodes in the hidden layer was determined using a sensitivity analysis.For this purpose, many predicted models were developed by changing the number of neurons or nodes in the hidden layer from one to 10.
en, the performance of each predicted model was measured using the correlation coe cient, R, and compared together to obtain the optimum number of nodes in the hidden layer (all performance measures used in this study are explained in the following section).is comparison for both cases where LogSigmoid or TanSigmoid transfer functions are used in the hidden layer is shown in Figure 11.As Figure 11 indicates, generally by increasing the number of nodes in the hidden layer the accuracy of the predicted models increased.For both Log-Sigmoid and TanSigmoid transfer functions, the highest performance was achieved for having 8 nodes in the hidden layer.However, as Figure 11 indicates, the di erence between the R values corresponding to the ANN-based models of having three and eight nodes in the hidden layer was negligible such that the percentage of di erence corresponding to the Log-Sigmoid and the TanSigmoid transfer functions was 0.3% and 0.4%, respectively.It should be noted that by increasing the number of nodes in the hidden layer, the number of terms of the mathematical formula of the predicted model increases.It therefore makes the application of such formula di cult for the users.Hence, in this study in order to propose more practical mathematical formulas, one, two, and three nodes were utilized to make up the hidden layer.It is clear that among these three models, the one of three nodes in the hidden layer has the highest accuracy but the lowest applicability.
In this study, performance of the ANN-based models was evaluated using four parameters, including the average absolute percentage error (Err), mean squared error (MSE), mean absolute error (MAE), and correlation coe cient (R).
ese parameters are given by (9).In these equations, Y i , T i ,  Advances in Civil Engineering Y, T, and N are the predicted output, the actual output obtained from nite element models, the average of the predicted outputs, the average of the actual outputs, and the number of samples, respectively.
Over tting is one of the most common problems during neural network training which potentially can lead to the achievement of a not optimally tted network.In this case, despite that the error associated to the training data set is a small value, the trained network is not able to adequately predict the target values when new inputs are given to the network.In this study, in order to overcome this problem, the early stopping technique was utilized [38].In this technique, the total data is divided into three sets: train, validation, and test.In this case, training of the network will be stopped when for the validation data set the network gives the lowest value of error even if further training of the network leads to the lower error with respect to the training data set.In such case, the network is optimally tted, and the weights and biases at the minimum of the validation error are returned and used to develop the mathematical formulas.is technique is graphically shown in Figure 12.In the present study, 60% of the total data set which were selected randomly was used for training of the network, 20% of the total data set was used for validating the network, and the rest of data were utilized for testing.In addition, in order to speed up the learning process and to achieve more reliable ANN-based model, input and output parameters were standardized using (10).In this equation, X std , X i , X, and SD are standardized parameter, initial value of the parameter, mean value, and standard deviation of the parameter, respectively.
In this study, a 3-layer network based on the LM/BP learning algorithm with 3 input parameters and 1 hidden layer made up of 1, 2, and 3 neurons were constructed using the neural network toolbox in MATLAB.For instance, the results of the predicted strength values using the ANN method for the training, the validation, and the testing data in the case of having two nodes in the hidden layer with LogSigmoid transfer functions are shown in Figure 13.As Figure 13 indicate for all three sets of data (training, validation, and testing data), the R values are close to unity, indicating the high performance of the proposed model.Other performance measures for all other predicted models are summarized in Table 2.As the results presented in Table 2 indicate, all ANN-based models are able to predict the connection strength with an average absolute percentage error less than 6%.

Investigation of the Adequacy of Sample Size.
In addition to the number of nodes in the hidden layer, one factor that can also in uence network performance is the quantity of the samples.In this study, the results of 52 nite element models which were analyzed under cyclic loading and 8 experimental models were utilized to make up the sample space. is number of nite element models was obtained by changing the input parameters (α f , α w , and L b /d b ) in their reasonable ranges.By considering this fact that the reasonable range of the input parameters is very limited and also running such nite element analysis is very timeconsuming and also due to the limited number of credible experimental tests, the sample size was limited to 60.It should be noted that, in some cases, this number was within the acceptable range proposed in [39] for sample space used in ANN analyses.In [39], a value of 10 times of the number of degree of freedom is proposed as the minimum size of the sample space.In the present study, by having 3 input variables in the input layer and 1, 2, or 3 nodes in the hidden layer, this requirement leads to a sample space of 30, 60, and 90 data, respectively.Hence, it can be concluded that this requirement has been met in the case of using one and two nodes in the hidden layer, but it was not the case when 3 nodes were used in the hidden layer.However, it is worth to mention that the number of data in the sample space can be quali ed by investigating the performance of the network.In this study, the performance was evaluated using the coe cient of correlation, R. Since for both of the networks having two or three nodes in the hidden layer the R values were greater than 0.95 (R 0.98 (0.97)) for three (two) nodes in the hidden layer, it was concluded that the number of data in sample space (60) was adequate.
is procedure was also adopted by other researchers [40][41][42] where, due to limitations of the problem under study, the size of sample space was remarkably less than the minimum sample size proposed in [39].Early stop technique to overcome over tting problem [38].
Advances in Civil Engineering Hidden layer made up of two neurons: Hidden layer made up of three neurons: In these equations, M P is the beam nominal plastic moment capacity.Constants C 1 through C 16 are summarized in Table 3.

Sensitivity Analysis.
At the end of Section 5, several questions were asked.e sixth one was about the most effective parameters that impact on the connection's strength.Here, in order to answer this question, the importance of each input parameter in the ANN-based models is evaluated through a sensitivity analysis using Garson's method [43]. is method was also used by other researchers (e.g., [44,45]) in different applications.For instance, Figure 14 graphically summarizes the steps of Garson's method for an ANN Advances in Civil Engineering composed of three inputs with two nodes in its hidden layer.Following this method, the input-hidden and hidden-output weights of the trained ANN models are partitioned, and the absolute values of the weights are taken to calculate the relative importance values [45].For connection's strength, the importance of each independent parameter (i.e., parameters α f , α w , and L b /d b ) was measured and is summarized in Table 4.As the results of Table 4 indicate, values of the relative importance are remarkably depended on the number of nodes in the hidden layer, while sensitivity to the type of the transfer function is negligible.For any type of the ANN-based models, parameter L b /d b was the most e ective parameter, while the less e ective parameter, depending on the number of nodes in the hidden layer, was either α f or α w .

Veri cation of the Proposed ANN-Based Models.
In this section, the accuracy of the ANN-based models for predicting the strength of the WFP connections was evaluated using the numerical and available test results presented in Table 1.For instance, Table 5 summarizes the key parameters of the experimentally tested specimens (specimens 53 through 60) along with the average absolute percentage error obtained using three of the proposed ANN-based models.
As the results presented in Table 5 indicate, all of the proposed models were able to adequately predict the connection's strength with the average absolute percentage error (Err) of 15.8%, 9.6%, and 9% for ANN-based models of having one, two, and three nodes in the hidden layer, respectively.However, as it is evident from these results, (22) which is corresponding to the ANN-based models of three nodes in the hidden layer is more reliable than ( 15) and (19).Furthermore, in order to have a comparison between the accuracy of ( 22) and the one on basis of the theoretical

Hidden layer
Steps of Garson's method [43] for determining the relative importance of input variables.4)), the value of connection strength at 4% total rotation (M/M P ) for all specimens presented in Table 1 was calculated using both equations and results are compared to the related experimental and numerical ones and associated absolute errors are calculated.Figure 7 shows this comparison.As Figure 7 indicates, ( 22) is more reliable than (4) such that the average of absolute error corresponding to this equation was only 3.74%.It is because of this fact that, as opposed to ( 22), ( 4) is not able to take into account the effect of parameter L b /d b .Note that, as the results of Section 6.4 indicated, this parameter has a significant effect on the value of connection strength.

Conclusions
is study was aimed to investigate the effects of the beam flanges and web slenderness and the beam length to the beam overall depth ratio (L b /d b ) on the strength and ductility of the welded flange plate (WFP) connections.
e conclusions drawn from this finite element study can be summarized as follows: (1) Fragility curves obtained on basis of the numerical and the experimental results indicated that the probabilities of achieving adequate strength (M cf /M pb ≥ 0.8) and ductility (θ ≥ 4%) for a WFP connection are 49% and 92%, respectively.In other words, for a WFP connection, its strength is of more concern than its ductility.(2) For a WFP connection, which is used in a special moment frame, limiting the width-to-thickness ratios of the beam flange and web plates to 0.3 ���� � E/f y  and 2.45 ���� � E/f y  , respectively, may not lead to the achievement of the adequate connection's strength if the beam overall depth is greater than 400 mm.
(3) In WFP connections, when the beam depth is less than 800 mm, 0.3 and 2.2 might be the maximum appropriate values for parameters α f and α w , respectively, to achieve the minimum required connection strength requested by AISC [15].While for deeper beams, 0.275 and 2.2 can be proposed as the maximum allowable values for parameters α f and α w , respectively.(4) Based on the finite element results, when the beam depth is less than 800 mm, 0.3 and 1.9 were the best values for parameters α f and α w , respectively, to achieve the highest connection performance.While for deeper beams, 0.275 and 1.9 might be the most suitable values for parameters α f and α w , respectively.(5) e proposed theoretical equation (i.e., (4)) which was based on the theoretical formulations relevant to the flexural strength of beams with noncompact elements was able to adequately predict the connection strength with an average of absolute error of 12.5%.is equation was able only to take into account parameters α f and α w .( 6) All six ANN-based models which were developed in terms of the beam length to the beam overall depth ratio and the slenderness ratios of the beam flanges and web plates showed adequate accuracy to predict the strength of the WFP connections at 4% total rotation.ese new formulates can be employed for assessment and design of WFP connections by practical engineers.(7) Among seven predicting models (six ANN-based models and one theoretical formula), the one on basis of the ANN with three nodes in the hidden layer and with the LogSigmoid transfer function was the most accurate predicting model.(8) For any type of the ANN-based models, parameter L b /d b was the most effective parameter.Its relative importance, depending on the number of nodes in the hidden layer, was 74%, 50%, and 61% for one, two, and three nodes in the hidden layer, respectively.(9) e less effective parameter, depending on the number of nodes in the hidden layer, was either α f or α w .(10) e relative importance of parameter α f was 16%, 23%, and 18% for one, two, and three nodes in the hidden layer, respectively.(11) e relative importance of parameter α w was 10%, 27%, and 21% for one, two, and three nodes in the hidden layer, respectively.

Figure 3 :
Figure 3: Comparison between experimental and analytical results of specimen LF30.
(a) Experimental test done by Kim et al.[12](b) Finite element model

5. 5 .
Fifth Question.How can we estimate the strength of a WFP connection based on parameters α f , α w , and L b /d b ?

45 MFigure 5 :
Figure 5: E ect of parameters α f , α w , and L b /d b on the connection strength.

5 . 6 .
Sixth Question.In order to reach adequate connection's strength, which one of the parameters α f , α w , and L b /d b can prediction-Eq (4) Average absolute error-Eq (4) ANN prediction -3 nodes in the hidden layer-Logsig Average absolute error-ANN

Figure 7 :
Figure 7: Comparison between proposed theoretical and ANN-based equations.

Figure 8 :
Figure 8: Structure of a typical ANN.

Figure 11 :
Figure 11: Relation between R and the number of nodes in the hidden layer for di erent transfer functions.

Table 1 :
Effect of parameters α f and α w on the connection's strength and ductility.
6A � not available.6AdvancesinCivilEngineeringWhilefor deeper beams, 0.275 and 2.2 might be the maximum appropriate values for parameters α f and α w , respectively.5.4.Forth Question.What can be the best values for parameters α f and α w to achieve the highest connection performance?By comparing the results presented in Table

Table 2 :
Performance statistics for strength's predicted model.Figure 13: Results of predicted connection strength using ANN model of two nodes in the hidden layer with Logsig transfer function.

Table 4 :
Relative importance values for di erent ANN-based models.

Table 5 :
Comparison between test and ANN-based model results in terms of average absolute percentage error.