Application of Component-Based Mechanical Models and Artiﬁcial Intelligence to Bolted Beam-to-Column Connections

: Top and seat beam-to-column connections are commonly designed to transfer gravitational loads of simply supported steel beams. Nevertheless, the ﬂexural resistance characteristics of these type of connections should be properly taken into account for design, when a reliable analysis of semi-rigid steel structures is desired. In this research paper, different component-based mechanical models from Eurocode 3 (EC3) and a literature proposal (by Kong and Kim, 2017) are considered to evaluate the initial stiffness ( S j,ini ) and ultimate moment capacity ( M n ) of top-seat angle connections with double web angles (TSACWs). An optimized artiﬁcial neural network (ANN) model based on the artiﬁcial bee colony (ABC) algorithm is proposed in this paper to acquire an informational model from the available literature database of experimental test measurements on TSACWs. In order to evaluate the expected effect of each input parameter (such as the thickness of top ﬂange cleat, the bolt size, etc.) on the mechanical performance and overall moment–rotation ( M– θ ) response of the selected connections, a sensitivity analysis is presented. The collected comparative results prove the potential of the optimized ANN approach for TSACWs, as well as its accuracy and reliability for the prediction of the characteristic ( M– θ ) features of similar joints. For most of the examined conﬁgurations, higher accuracy is found from the ANN estimates, compared to Eurocode 3- or Kong et al.-based formulations.


Introduction
Bolted top-seat angle connections without (TSAC) or with web angles (TSACW) have been extensively used in steel and composite structures, because of their relatively high moment capacity and easy construction. These types of connections are mainly designed to resist gravity loads of determinate steel beams. A schematic layout is proposed in Figure 1. Based on the AISC standard [1], a given TSACW is classified as a semi-grid beam-to-column connection, in which the moment-rotation (M-θ) behavior needs to be considered in the overall force distribution and structural analysis. Accordingly, there are no doubts on the necessity to accurately predict the (M-θ) behavior of TSACW joints.
The (M-θ) behavioral trends of steel connections represent essential features of their response, and in particular with regard to the initial stiffness S j,ini , maximum moment capacity M n , and maximum rotation θ u . Once these features are correctly estimated, a reliable simulation of the actual behavior toward the optimum design of structures would be possible. The behavior of TSACs and TSACWs has been investigated and equivalent equations have been proposed to simulate their mechanical behavior for applications buildings", namely, off-site prefabricated building modules in which the primary struc-ture consists of steel frames, whose construction has become increasingly popular due to the advantages, including reduced resource wastage and improved quality [7][8][9][10]. The typical connection takes the form of a bolted connection with a plug-in device or a bolted connection with a rocket-shaped tenon. Being classified as "bolted", these connections have the advantages of reduced site work and demountability.
In both the cases, the mechanical properties of connections have a great influence on the strength and stability of the structural system to which they belong. An unrealistic design of connections may negatively affect the serviceability of the structure as a whole, due to the occurrence of large deflections. Deformations in connection components as a result of local buckling phenomena is also a common reason for the failure of these connections. The key role of connections for modular buildings has been recognized to affect the overall structural performance of modular systems. Despite this importance, the design of reliable connections is still identified as a major challenge, and prevalent knowledge of their structural performance is still limited [11,12].  In this paper, an original study is presented for TSACWs. An artificial neural network (ANN) model able to capture the underlying mechanism and to accurately estimate both the initial stiffness Sj,ini and the ultimate moment capacity Mn for TSACWs is proposed. To In both the cases, the mechanical properties of connections have a great influence on the strength and stability of the structural system to which they belong. An unrealistic design of connections may negatively affect the serviceability of the structure as a whole, due to the occurrence of large deflections. Deformations in connection components as a result of local buckling phenomena is also a common reason for the failure of these connections. The key role of connections for modular buildings has been recognized to affect the overall structural performance of modular systems. Despite this importance, the design of reliable connections is still identified as a major challenge, and prevalent knowledge of their structural performance is still limited [11,12].
In this paper, an original study is presented for TSACWs. An artificial neural network (ANN) model able to capture the underlying mechanism and to accurately estimate both the initial stiffness S j,ini and the ultimate moment capacity M n for TSACWs is proposed. To this aim, two reference calculation approaches are taken into account from the literature, and investigated in detail, namely, (i) the mechanical component-based model approach proposed by Eurocode 3 ("EC3" [13]) and (ii) the analytical formulation proposed by Kong and Kim in [14] ("KK"). Since the S j,ini and M n features for a given TSACW are the most influential parameters that affect the overall mechanical characterization of a given bolted beam-to-column connection, these two selected properties are preliminarily calculated based on the selected "EC3" and "KK" formulations, and compared against each other, with the support of experimental data from the literature (77 specimens in total). Successively, a comprehensive ANN approach, combined with a metaheuristic artificial bee colony (ABC) algorithm, is developed in this paper to extract an informational model for TSACWs. A major advantage for the definition of the proposed model is given by the support of experimental data from the literature, and sensitivity analysis to assess the collected comparative results. The potential and limits of each scheme, as well as the original ABC-ANN for TSACW joints, are hence discussed in detail.

State of the Art
Different parameters contribute to the (M-θ) behavioral characterization of bolted connections, as shown in Figure 2 [15][16][17][18][19][20], and it is thus first necessary to use reliable modeling techniques. In the case of modular steel buildings, non-linear link elements are generally incorporated for computational efficiency, where the bending moment, shear, axial force, and corresponding moment-rotation (M-θ) curves are simplified and included as spring elements. Nevertheless, since the knowledge of the structural behavior of modular steel buildings is limited, the simplified models for the force-displacement and momentrotation curves are not well established. Therefore, the intrinsic uncertainty can lead to excessively conservative design practices, requiring additional investigations on the topic. this aim, two reference calculation approaches are taken into account from the literature and investigated in detail, namely, (i) the mechanical component-based model approach proposed by Eurocode 3 ("EC3" [13]) and (ii) the analytical formulation proposed by Kong and Kim in [14] ("KK"). Since the Sj,ini and Mn features for a given TSACW are the mos influential parameters that affect the overall mechanical characterization of a given bolted beam-to-column connection, these two selected properties are preliminarily calculated based on the selected "EC3" and "KK" formulations, and compared against each other with the support of experimental data from the literature (77 specimens in total). Succes sively, a comprehensive ANN approach, combined with a metaheuristic artificial bee col ony (ABC) algorithm, is developed in this paper to extract an informational model fo TSACWs. A major advantage for the definition of the proposed model is given by th support of experimental data from the literature, and sensitivity analysis to assess the col lected comparative results. The potential and limits of each scheme, as well as the origina ABC-ANN for TSACW joints, are hence discussed in detail.

State of the Art
Different parameters contribute to the (M-θ) behavioral characterization of bolted connections, as shown in Figure 2 [15][16][17][18][19][20], and it is thus first necessary to use reliabl modeling techniques. In the case of modular steel buildings, non-linear link elements ar generally incorporated for computational efficiency, where the bending moment, shear axial force, and corresponding moment-rotation (M-θ) curves are simplified and in cluded as spring elements. Nevertheless, since the knowledge of the structural behavio of modular steel buildings is limited, the simplified models for the force-displacemen and moment-rotation curves are not well established. Therefore, the intrinsic uncertainty can lead to excessively conservative design practices, requiring additional investigation on the topic.  Extended literature studies are indeed available for TSACs and TSACWs. To simulate the beam-to-column connection behavior, apart from experimental tests, three different modeling options are introduced in the form of analytical or empirical models [21], advanced finite element (FE) models [22][23][24], and component mechanical models [25,26].
Analytical and empirical models naturally extract simple mathematical expressions, and thus typically facilitate their application in practice, with a reasonable level of accuracy. On the other hand, by applying advanced FE methods, it is possible to simulate more precisely the complex non-linear beam-to-column connection response. Some limitations are unavoidably introduced in terms of high computational cost and sophisticated preparation process (assembly, calibration, etc.). Component-based mechanical models, finally, consider an assembly of rigid elements and equivalent springs relay on the two previous modeling approaches, considering computational complexity and reliability.
Several research studies have documented the (M-θ) performance of TSACWs, especially through extended experimental investigations. In 1985 and 2000, respectively, Azizinamini and Radziminski [25] and Calado et al. [15] experimentally investigated the (M-θ) behavior of TSACW specimens under monotonic and cyclic loading. The hysteresis behavior of TSACW specimens was investigated in several cited projects. The presented results indicated that the main sources of plastic deformation are typically located in the top and bottom cleats, which represent the governing components. In addition, the same studies proved a limited influence of the column size on the overall hysteretic performance of the examined connections. A considerable number of literature studies has also been published in support of the prediction of the (M-θ) behavior, initial stiffness, and ultimate strength of TSACW systems through various mechanical and analytical models. In 1987, Azizinamini et al. [27] proposed a set of relevant equations to predict initial stiffness. Kishi et al. [28,29] presented an analytical model to estimate both the initial stiffness and ultimate moment capacity of TSACWs through cantilever beam theory. In more detail, a three-parameter power model was discussed in [29], in which the elastic stiffness and ultimate moment capacity of connection of TSACWs were estimated by a simple analytical method. Pucinotti [30] proposed a simplified mechanical model along with relevant initial stiffness and ultimate moment equations. A three-linear (M-θ) constitutive law based on the fiber element formulation was proposed by Shen and Astaneh-Asl [31]. In that study, two possible mechanisms were considered, according to the strength of the angles relative to the bolt. They concluded that plastic hinges were crated at the leg of the angles in thinner top cleats, while the plastic hinge was generally found to appear at the central line of the column bolts, along with plastic deformations of the bolts. A mechanical model for estimating the inelastic cyclic (M-θ) relationship of TSACWs was proposed by De Stefano et al. [32]. Currently, by the advances in computer technologies, there is a large volume of published studies on (M-θ) behavior of beam-to-column connections through the FE method. For example, Danesh et al. [33] investigated the effect of shear force on the expected initial stiffness. Pirmoz et al. [34] focused on the effect of web angle on the (M-θ) behavior. Salem et al. [35] conducted parametric research on the prying action of bolts. Their study indicated that the prying force of top cleat and column flange bolts is found to increase by decreasing the vertical leg of the top cleat. Collectively, all the mentioned studies (and others) outline a set of different methodologies in support of a reliable (M-θ) modeling and mechanical characterization of beam-to-column connections.
Among the multitude of research studies that have been done on TSACs or TSACWs, traditional experimental or numerical methods are used in most of the cases. The current investigation, in this regard, applies ANN techniques for the analysis of this type of connection. ANNs are largely used for solving complex civil engineering (and other) problems, and their application has increased significantly in the last few years. ANNs can be effectively applied for predictive modeling in different engineering fields, especially in those cases where some prior (experimental or numerical) analyses are already available. The origins of ANNs can be found in the field of biology, where the biological brain consists of billions of highly interconnected neurons, forming a neural network.
A validated informational method is developed in this paper, to provide an alternative tool that could be used to model the complex structural and material behavior in TSACWs that cannot be easily approximated and generalized by conventional approaches. The information about the underlying mechanics for TSACWs is extracted from the observed ex-perimental data of the literature, which are processed in neural networks and subsequently trained with optimization techniques.

Data Bank Development
Having recognized that different parameters contribute to the (M-θ) behavioral characterization of bolted beam-to-column connections, it is first fundamental to develop a consistent data bank of test results for the examined connection typology. Experimental investigations and their results-when correctly extracted-notoriously allow a more robust and accurate classification of different behavioral features for beam-to-column connections, including S j,ini and M n , but also hardening, non-linearities, progressive damage and degradation of mechanical parameters, rotation capacity, failure mechanism, and sequence. The same database is also strictly necessary for the development of the ANN model herein proposed.
In this paper, the preliminary verification of the component-based mechanical model for bolted TSACWs, as well as the training data to develop the ANN model, is carried out with the support of experimental data from the literature. Table 1 presents a short summary of the major geometrical and material properties for the examined TSACWs. The full data of TSACW specimens are available in [27,36,37].

Initial Stiffness
The initial stiffness S j,ini of bolted connections is one of the most important parameters for the characterization of their (M-θ) response. There is a large volume of published studies investigating S j,ini of bolted beam-to-column connections. The first deep discussion and analysis about S j,ini emerged in the 1980s with the research studies by Azizinamini et al. [27]. The authors proposed a beam model dividing the leg of the top cleat flange into two types of beam segments, flexible and rigid sections. Relevant equations of practical use have been presented to estimate S j,ini for TSACWs. Later, Kishi and Chen [29] presented a new equation that was developed according to the cantilever beam model.
The EC3 standard document [13]-Annex J-implemented a component method to estimate the S j,ini of TSACs. In this method, the connection behavior is simulated by a series of different components, each one represented in the form of an elastic spring with specific stiffness and strength. The overall initial stiffness of the connection is hence notoriously estimated by assembling all the springs in a parallel series configuration, as shown in Figure 3. compression (K2), the flexural stiffness of column flange (K4), the top cleat flexural stiffness (K6), the tensile stiffness of bolts (K10), and, for non-preloaded bolts, their shear stiffness (K11) and their bearing stiffness (K12). In conclusion, Sj,ini of TSACs is given by: where E is Young's modulus, z is the lever arm, Ki is the i-th component stiffness coefficient, and n is the number of joint components, as in Figure 3. Finally, z should be taken as the distance from the bolt-row in tension and mid-thickness of the leg of the seat cleat on the compression flange. While the above formulation is widely used in practice, the EC3 formulation does not include a mechanical model for TSACWs. Accordingly, an extension of EC3 for TSACWs is proposed in this paper. For a bolt-row in a web cleat, the stiffness illustrated in Figure  4 should be considered. In the EC3, the input stiffness coefficients that must be considered are associated with the column shear panel zone (K 1 ), the column flange in tension (K 3 ), the column flange in compression (K 2 ), the flexural stiffness of column flange (K 4 ), the top cleat flexural stiffness (K 6 ), the tensile stiffness of bolts (K 10 ), and, for non-preloaded bolts, their shear stiffness (K 11 ) and their bearing stiffness (K 12 ). In conclusion, S j,ini of TSACs is given by: where E is Young's modulus, z is the lever arm, K i is the i-th component stiffness coefficient, and n is the number of joint components, as in Figure 3. Finally, z should be taken as the distance from the bolt-row in tension and mid-thickness of the leg of the seat cleat on the compression flange. While the above formulation is widely used in practice, the EC3 formulation does not include a mechanical model for TSACWs. Accordingly, an extension of EC3 for TSACWs is proposed in this paper. For a bolt-row in a web cleat, the stiffness illustrated in Figure 4 should be considered.
The overall stiffness of basic components illustrated in Figure 4 is represented by a single equivalent stiffness coefficient k eq , that can be calculated as: where S j,ini is again the initial stiffness for TSACs, while z r represents the distance between the center of the compression cleat and the bolt-row r of the web cleat. Moreover : where k i,r is the stiffness coefficient representing component i relative to bolt-row r, and: In conclusion, the initial stiffness of TSACWs can be calculated as: The overall stiffness of basic components illustrated in Figure 4 is represented by a single equivalent stiffness coefficient keq, that can be calculated as: where Sj,ini is again the initial stiffness for TSACs, while zr represents the distance between the center of the compression cleat and the bolt-row r of the web cleat. Moreover: where ki,r is the stiffness coefficient representing component i relative to bolt-row r, and: In conclusion, the initial stiffness of TSACWs can be calculated as: Kong and Kim [14] have used curve fitting software to obtain the effects of top and seat angles on the Sj,ini value of TSACWs. Past literature studies also recognized that Sj,ini Kong and Kim [14] have used curve fitting software to obtain the effects of top and seat angles on the S j,ini value of TSACWs. Past literature studies also recognized that S j,ini correlates well with correlates the angle thickness, the length of the top angle, the height of the beam, the thickness of the column flange, the fillet size of the top angle, the thickness of the beam web, the gauge distance, and the diameter of bolts. In [14], therefore, a semi-empirical equation has been proposed as: In Equation (6), E is the previously defined Young's modulus; l t is the length of top angle; t t is the top angle thickness; t s is the thickness of seat angle; d is the height of beam; k t is the fillet size of top angle; t cf is the thickness of column flange; t bw is the thickness of beam web; g t is the distance from top angle heel to the center of bolts; d b is the diameter of bolts; n is the number of bolts; α = 1.0 mm; l p is the angle length of the web; ta is the angle thickness of web; υ is Poisson's ratio; and g c = g 1 − t a , as shown in Figure 5.

Ultimate Moment Capacity
The ultimate moment capacity M n is one of the well-known fundamental parameters that affects the overall response of bolted beam-to-column connections. Several literature studies on TSACW specimens suggest that a rotation of 0.03 rad is sufficient to experience the full plastic moment capacity M p of the connected beam [38]. The American Institute of Steel Construction [39], in this regard, recommends considering a relatively low rotation during the analysis and design process for this connection typology. Accordingly, the amplitude of 0.03 rad is defined as the maximum rotation capacity of TSACWs.

Ultimate Moment Capacity
The ultimate moment capacity Mn is one of the well-known fundamental parameters that affects the overall response of bolted beam-to-column connections. Several literature studies on TSACW specimens suggest that a rotation of 0.03 rad is sufficient to experience the full plastic moment capacity Mp of the connected beam [38]. The American Institute of Steel Construction [39], in this regard, recommends considering a relatively low rotation during the analysis and design process for this connection typology. Accordingly, the amplitude of 0.03 rad is defined as the maximum rotation capacity of TSACWs.
The component-based method can be extended to estimate the ultimate moment capacity Mn of TSACWs, as far as the contribution of influencing components is properly considered. These all affect the overall TSACW flexural resistance and are represented by: the column web panel in shear, the column web in compression, the column web in tension, the column flange in bending, the top cleat in bending, the web cleat in bending, the bolts in tension or in shear, the beam web in tension, and the beam flange and web in compression. Past research studies proposed different mechanical models to calculate Mn [30,40,41]. In most of those proposals, the failure of the top, seat, and web cleat's leg under maximum moment was only recognized, while other failure mechanisms (such as the failure of one or more bolts) were fully disregarded. In 2005, the EC3 standard [13] defined the required ultimate moment capacity Mn of TSACs and TSACWs as: where Z is the lever arm and FRd is the design resistance of weak joint components. The second term can be governed by one of the following mechanisms: top cleat in bending (Ftc,Rd), bolts in tension, beam flange in tension and compression, and beam web in tension. For a TSAC with a single-web cleat, Kishi and Chen [29], in 1990, elaborated a model proposal for the detection of collapse mechanisms. Kong and Kim [14] further extended the method from [29], and presented a formulation for TSACWs: In Equation (8), Mos is the plastic moment capacity of the seat angle; Mp is the plastic moment capacity of the top angle; Vpt is the ultimate shear force acting on the top angle; The component-based method can be extended to estimate the ultimate moment capacity M n of TSACWs, as far as the contribution of influencing components is properly considered. These all affect the overall TSACW flexural resistance and are represented by: the column web panel in shear, the column web in compression, the column web in tension, the column flange in bending, the top cleat in bending, the web cleat in bending, the bolts in tension or in shear, the beam web in tension, and the beam flange and web in compression. Past research studies proposed different mechanical models to calculate M n [30,40,41]. In most of those proposals, the failure of the top, seat, and web cleat's leg under maximum moment was only recognized, while other failure mechanisms (such as the failure of one or more bolts) were fully disregarded. In 2005, the EC3 standard [13] defined the required ultimate moment capacity M n of TSACs and TSACWs as: where Z is the lever arm and F Rd is the design resistance of weak joint components. The second term can be governed by one of the following mechanisms: top cleat in bending (F tc,Rd ), bolts in tension, beam flange in tension and compression, and beam web in tension. For a TSAC with a single-web cleat, Kishi and Chen [29], in 1990, elaborated a model proposal for the detection of collapse mechanisms. Kong and Kim [14] further extended the method from [29], and presented a formulation for TSACWs: In Equation (8), M os is the plastic moment capacity of the seat angle; M p is the plastic moment capacity of the top angle; V pt is the ultimate shear force acting on the top angle; d 2 is a parameter related to the depth of the beam and thickness of top and seat cleats; V pa is a parameter that depends on the ultimate shear force at the upper and lower edges of the web cleat; and d 4 , finally, is the distance between plastic shear at the lower edge of the web cleat to the center of compression. These parameters are schematized in Figure 6.
d2 is a parameter related to the depth of the beam and thickness of is a parameter that depends on the ultimate shear force at the upp the web cleat; and d4, finally, is the distance between plastic shear a web cleat to the center of compression. These parameters are schem Figure 6. Collapse mechanism for TSACW. Reproduced from [14] with pe copyright agreement with license number 5011231057466, February 2021.

Informational Based Modeling
An informational-based method recognizes, as an alternative m complex structural and material behavior, that conventional metho mated. This alternative method simulates the behavior by the inf selected tested specimens. Accordingly, this is an essential evoluti equations to preserve data that contain the necessary information teristics. In this approach, the underlying mechanics information is perimental test data and processed in the ANN program.
Among numerous structures of ANNs that have been studied one is represented by the multi-layered feed-forward (FF) network implemented in the current investigation. Among others, FF is the existing ANNs, representing the most popular class of ANNs for ciency [42]. The main structure of a neural network is usually ma layers, and never cycles. According to the literature, three-layer FF be sufficient in civil engineering practices [43]. The input layer is w Figure 6. Collapse mechanism for TSACW. Reproduced from [14] with permission from Elsevier, copyright agreement with license number 5011231057466, February 2021.

Informational Based Modeling
An informational-based method recognizes, as an alternative method to simulate the complex structural and material behavior, that conventional methods are not simply estimated. This alternative method simulates the behavior by the information provided by selected tested specimens. Accordingly, this is an essential evolution from mathematical equations to preserve data that contain the necessary information on mechanical characteristics. In this approach, the underlying mechanics information is extracted from the experimental test data and processed in the ANN program.
Among numerous structures of ANNs that have been studied, the most widely used one is represented by the multi-layered feed-forward (FF) network, which is the structure implemented in the current investigation. Among others, FF is the oldest and simplest of existing ANNs, representing the most popular class of ANNs for its computational efficiency [42]. The main structure of a neural network is usually made up of three distinct layers, and never cycles. According to the literature, three-layer FF networks are found to be sufficient in civil engineering practices [43]. The input layer is where the data are introduced to the model, the hidden layer is where the data are processed, and the output layer is where the model results are generated. Each layer is made up of nodes called neurons [44]. Apart from the neurons in the input layer, which only receive and transmit incoming signals to other neurons in the hidden layer, each neuron in the other layers consists of three main components, weights, bias, and an activation function, which can be linear or non-linear. After determining the network topology, including the number of layers, the number of neurons, the type of transmission function, and the network learning algorithm, modeling for data can be started [45]. During the analysis of the model using learning algorithms, the number of errors can be reduced by adjusting weights and bias in each neuron. For this purpose, the Levenberg-Marquardt algorithm has been used in modeling the neural network due to its high speed and accuracy. The training algorithm distributes the network error to achieve optimal or minimum error [46,47].

Artificial Bee Colony (ABC) Algorithm
The artificial bee colony (ABC) algorithm is an optimization algorithm based on the bee population's collective intelligence and intelligent behavior in finding food [44]. In its primary model, the algorithm performs a neighborhood search combined with a random search, and can be used for either combined or functional optimization [48].
In the ABC algorithm, the colony consists of three groups of bees: employed, onlookers, and scout bees. The first and the second half of the colony consist of employed artificial bees and onlookers, respectively. There is only one employed bee for every food source. In this algorithm, moreover, a scout is representative of an employed bee of an abandoned food source. Overall, the search process in ABC can be summarized as follows [49,50]: • A food source is determined by employed bees in their memory within the neighborhood.

•
The collected information of food sources by employed bees is shared with onlookers within the hive, and subsequently, the optimum food sources will be selected by onlookers.

•
A food source will be selected by onlookers themselves within the neighborhood of the food sources. • An employed bee becomes a scout once the food source has been abandoned and starts to search for a new food source randomly.

Training the ANN and Methodology
To train the artificial neural network, a total number of 77 specimens from [36] and [37] was considered in this study, in order to determine two outputs of S j,ini and (M n /M p , beam ). Up to ≈ 80% of the selected samples (62) was used for training, while 20% (15 specimens) was used for testing the network.
Accordingly, several variables were introduced as input parameters for the model, including: The statistical characteristics of these variables are summarized in Table 2. The number of hidden layers and the total number of neurons in the hidden layers in an ANN depend on the nature of the problem [51]. Generally, the trial and error method is used to obtain the ideal architecture that best reflects the characteristics of laboratory data. In this paper, an innovative method for calculating the number of neurons in hidden layers is taken into account, namely: where N H represents the number of neurons in the hidden layers and N I is the number of input variables. Since the number of influential input variables for the current study is equal to 6, the empirical Equation (9) shows that the number of neurons in hidden layers can be less than 13. Therefore, several networks with different topologies (with a maximum of two hidden layers and a maximum of 13 neurons) were trained and studied in this paper.
The hyperbolic tangent stimulation function and Levenberg-Marquardt training algorithm were used in all networks. The statistical indices used to evaluate the performance of different topologies are the root mean squared error (RMSE), the average absolute error (AAE%), the model efficiency (EF), and the variance account factor (VAF%) that are defined as follows [52]: After examining different topologies of networks, it was found that the network with a 6-7-6-1 topology is characterized by the lowest value of error in RMSE, AAE%, EF, and VAF%, and by the highest value of R 2 to estimate the two output parameters S j,ini and (M n /M p,beam ), as also shown in Table 3. It is necessary to mention that the error criteria for training and testing the selected data are calculated in the main range of variables and not in the normal range.
In this study, as it is for several structural engineering practice applications, the ABC algorithm is used as a new metaheuristic algorithm to determine the weight optimization of each ANN model. In more detail, the (M n /M p , beam ) output has a numerical range of 0-1 while the S j,ini output is characterized by a numerical range of 1600-37,000 (kNm/rad). This means that there is a big difference between the two target outputs. For modulation, two separate ANNs were hence used in this study, each one with one output. Their basic characteristics are summarized in Table 4. Figure 7 shows the optimal topology of an FF network with two hidden layers, six input variables (neurons), and one output parameter. The ABC has been also used to provide the least prediction error for the trained structure, in order to optimize the weights and biases of the ANN model. The ABC parameters are also presented in Table 5. In this study, as it is for several structural engineering practice applications, the ABC algorithm is used as a new metaheuristic algorithm to determine the weight optimization of each ANN model . In more detail, the (Mn/Mp,beam) output has a numerical range of 0-1 while the Sj,ini output is characterized by a numerical range of 1600-37,000 (kNm/rad). This means that there is a big difference between the two target outputs. For modulation, two separate ANNs were hence used in this study, each one with one output. Their basic characteristics are summarized in Table 4.  Figure 7 shows the optimal topology of an FF network with two hidden layers, six input variables (neurons), and one output parameter. The ABC has been also used to provide the least prediction error for the trained structure, in order to optimize the weights and biases of the ANN model. The ABC parameters are also presented in Table 5.     Table 6 shows the comparison between different models with experimental data for estimating S j,ini and M n . Concerning the calculated average ("Avg.") and standard deviation ("STD") values, the results in Table 6 indicate that the ABC-ANN model provides more reliable predictions for both, compared to the EC3 or KK formulations described earlier.
Using the existing empirical models, in particular, the S j,ini and M n predictions for some specimens are either underestimated or overestimated, and this suggests the limitation of mechanical models to capture the underlying mechanism that governs both the parameters. The ABC-ANN predictions, conversely, are characterized by minimum deviation.  Figures 8 and 9 show the scatter graph that provided the relationship between test results and the proposed ABC-ANN model for estimating the S j,ini and M n parameters, respectively. In this case, the comparative results also indicate that the ABC-ANN model offers a reliable value for the ratio of experimental to computational predictions (R 2 ), for both the examined mechanical parameters, and thus confirming further the high potential and accuracy of the proposed model. Appl. Sci. 2021, 11, x FOR PEER REVIEW 14 of 21     In conclusion, Tables 7 and 8 provide the final weights and biases for both hidden layers, as estimated by the ABC-ANN model. Using the values of the weights and biases between the different ANN layers, the two output parameters (S j,ini and M n /M p,beam ) can be determined and predicted. To formulate ANN results, weights and biases from Tables 7 and 8 should be normalized as: By substituting the normalized values of Table 3 for each one of the six input parameters, they are represented by a 6 × 1 vector that is herein labeled as a (1) . Then, by using the following equations, the values of S j,ini and (M n /M p,beam ) can be calculated from: The parameters IW, LW1, LW2, b 1 , b 2 , and b 3 are shown as vector matrices in Tables 7 and 8

Sensitivity Analysis
As previously discussed, the accuracy and potential of the proposed ABC-ANN model for the estimation of the Sj,ini and Mn parameters in TSACWs was acknowledged b comparison with the EC3 and KK component-based formulations. In order to investigat in more detail the effect of all the required input parameters, a sensitivity analysis (SA was carried out for the selected TSACW specimens. The SA reveals how significantly th model output can be affected by changes within input variables. There are two main type of SA, known as "global" and "local" sensitivity analyses, where the local sensitivity anal ysis (LSA) concentrates on the local impact of individual input parameters on the overal performance. The global sensitivity analysis (GSA), on the other hand, evaluates the in fluence of individual input parameters over their entire spatial range and measures th uncertainty of the overall performance (output) caused by input uncertainty, over inter action with other parameters, or also taken individually. Therefore, considering the natur of the complex non-linear behavior and variation of Sj,ini and Mn parameters in the curren study, the GSA was selected as much more rational for investigating the impact of inpu parameters on the overall performance.
Amongst GSA methods, a variance-based approach has been primarily considered in past literature for SA. The method provides a specific methodology for defining tota and first-order sensitivity indices for each input parameter of the ANN model. Assumin a model in the form Y = f (X1, X2,…, Xk), where Y is scalar, the variance-based techniqu takes a variance ratio to evaluate the impact of individual parameters using variance de composition as: The results shown in Figures 10 and 11 clearly indicate that the closest prediction for both the S j,ini and M n input parameters, to the point representing the experimental data of the literature, are provided by the herein developed ABC-ANN model. The EC3 componentbased model, as shown, also results in high values of root mean square-centered difference and standard deviation, thus further suggesting a good accuracy of the formulation over the selected experimental data. The same comparative parameters, conversely, are rather low regarding the application of the KK model to the selected experimental specimens.

Sensitivity Analysis
As previously discussed, the accuracy and potential of the proposed ABC-ANN model for the estimation of the S j,ini and M n parameters in TSACWs was acknowledged by comparison with the EC3 and KK component-based formulations. In order to investigate in more detail the effect of all the required input parameters, a sensitivity analysis (SA) was carried out for the selected TSACW specimens. The SA reveals how significantly the model output can be affected by changes within input variables. There are two main types of SA, known as "global" and "local" sensitivity analyses, where the local sensitivity analysis (LSA) concentrates on the local impact of individual input parameters on the overall performance. The global sensitivity analysis (GSA), on the other hand, evaluates the influence of individual input parameters over their entire spatial range and measures the uncertainty of the overall performance (output) caused by input uncertainty, over interaction with other parameters, or also taken individually. Therefore, considering the nature of the complex non-linear behavior and variation of S j,ini and M n parameters in the current study, the GSA was selected as much more rational for investigating the impact of input parameters on the overall performance.
Amongst GSA methods, a variance-based approach has been primarily considered in past literature for SA. The method provides a specific methodology for defining total and first-order sensitivity indices for each input parameter of the ANN model. Assuming a model in the form Y = f (X 1 , X 2 , . . . , X k ), where Y is scalar, the variance-based technique takes a variance ratio to evaluate the impact of individual parameters using variance decomposition as: where V is the variance of the ANN model output; V i is the first-order variance for the input X i ; V ij to V 1,2, . . . , k correspond to the variance of the interaction of the k parameters. V i and V ij , which denote the significance of the individual input to the variance of the output, are a function of the conditional anticipation variance: where the suffix x ∼i designates the set of all input variables apart from X i . The first-order sensitivity index (S i ) represents the first-order impact of an input X i on the overall output provided by: The abovementioned methodology for calculating the first-order sensitivity index was used in this research study. Major results of the SA are presented in Figure 12. where V is the variance of the ANN model output; Vi is the first-order variance for the input Xi; Vij to V1,2, …, k correspond to the variance of the interaction of the k parameters.
Vi and Vij, which denote the significance of the individual input to the variance of the output, are a function of the conditional anticipation variance: where the suffix x∼ designates the set of all input variables apart from Xi.
The first-order sensitivity index (Si) represents the first-order impact of an input Xi on the overall output provided by: The abovementioned methodology for calculating the first-order sensitivity index was used in this research study. Major results of the SA are presented in Figure 12.
Apart from the yield strength (fy) that depends on material properties, the comparative results indicate that the thickness of top flange (thtc) has the most influence, while the moment inertia ratio of column to beam (Icol/Ib) has the least effect on both the output parameters, Sj,ini and Mn. The thickness of the bottom flange cleat (thbc) can be classified as the second most influential input variable, especially in terms of maximum moment capacity Mn.

Concluding Remarks
Modeling the plastic response of different components in beam-to-column bolted connections and their interactions is a challenging issue in the structural engineering community. In this paper, an informational artificial neural network (ANN) model combined with the metaheuristic artificial bee colony (ABC) algorithm was developed to model the initial stiffness (Sj,ini) and maximum moment capacity (Mn) of top-seat angle connections with double web angles (TSACWs). Two different formulations of component-based mechanical models of the literature (as proposed by Eurocode 3 ("EC3") or by Kong and Kim in 2017 ("KK")) were also investigated in detail.
The discussed comparisons and results confirmed that the efficiency of the component-based mechanical model depends on the number and accuracy of the relationships of the constitutive components. By defining a sufficient number of components and subsequently idealizing the physical behavior in analytical equations, the reliable application of the mechanical model to different connection configurations is possible. Apart from the yield strength (f y ) that depends on material properties, the comparative results indicate that the thickness of top flange (th tc ) has the most influence, while the moment inertia ratio of column to beam (I col /I b ) has the least effect on both the output parameters, S j,ini and M n . The thickness of the bottom flange cleat (th bc ) can be classified as the second most influential input variable, especially in terms of maximum moment capacity M n .

Concluding Remarks
Modeling the plastic response of different components in beam-to-column bolted connections and their interactions is a challenging issue in the structural engineering community. In this paper, an informational artificial neural network (ANN) model combined with the metaheuristic artificial bee colony (ABC) algorithm was developed to model the initial stiffness (S j,ini ) and maximum moment capacity (M n ) of top-seat angle connections with double web angles (TSACWs). Two different formulations of component-based mechanical models of the literature (as proposed by Eurocode 3 ("EC3") or by Kong and Kim in 2017 ("KK")) were also investigated in detail.
The discussed comparisons and results confirmed that the efficiency of the componentbased mechanical model depends on the number and accuracy of the relationships of the constitutive components. By defining a sufficient number of components and subsequently idealizing the physical behavior in analytical equations, the reliable application of the mechanical model to different connection configurations is possible.
Nevertheless, idealization typically resulted in equations that excluded several important physical behavioral features of TSACWs, such as slippage. On the other hand, unlike the conventional mechanical modeling process that involves idealization from the observed behavior to the mathematical equations, in the informational base method, the information about essential behavior is extracted from available experimental test data and processed using an ANN. Nevertheless, the ANN model was limited to providing the global response only for bolted connections that include the contributions of all the constitutive components. Using this method, as shown, it is impossible to represent individual components and their actual contribution. Therefore, the model does not offer an insight into the underlying mechanics. Overall, the following conclusions can thus be derived: • both the EC3 and KK component-based models failed to capture the underlying mechanism for estimating S j,ini and M n parameters. As a result, these were either underestimated or overestimated for the reference specimens. On the other hand, the herein developed ABC-ANN model proved to offer a reliable prediction of required parameters, as also emphasized by the ratio of observational to computational values (R 2 ), and thus suggesting the high potential and accuracy of the proposal. • The ANN model combined with the ABC algorithm established an excellent agreement with the available experimental database. The results highlighted that the ANN model may be a reliable alternative to a component-based mechanical model to estimate the mechanical behavior of bolted beam-to-column connections. Using the values of weights and biases between the different ANN layers, the two output parameters (S j,ini and M n ) can be accurately predicted.

•
The sensitivity analysis confirmed that (apart from the yield strength f y that necessarily depends on material properties) the thickness of the top flange (th tc ) has a significant influence, while the moment inertia ratio of column to beam (I col /I b ) has the least effect on both the predicted output parameters, S j,ini and M n .
Author Contributions: This research paper results from a joint collaboration of all the involved authors. All authors contributed to the paper drafting and review, I.F., M.N., R.P. and C.B. All authors have read and agreed to the published version of the manuscript.
Funding: This work was supported by funding from the Federal State Autonomous Educational Institution of Higher Education South Ural State University (National Research University). Further, MDPI is also acknowledged for providing waived fees (Guest Editor C.B.) in support of standard publication APCs.