ANN-based optimized design of doubly reinforced rectangular concrete beams based on multi-objective functions

ABSTRACT In the structural engineering field, multi-objective optimization is a difficult task as demands of objective functions in structural designs are sometimes conflicting, such as costs and structural weight. Some previous studies which are mainly based on metaheuristics have investigated multi-objective optimizations to deal with such conflicts. However, the optimized results obtained by the populations-based methods do not present the tradeoff ratios contributed by each objective function. This study provides a five-step multi-objective optimization that can be used by structural engineers and decision-makers who are not familiar with optimization algorithms. The unified function of objectives (UFO) is optimized to derive the Pareto frontier. Design requirements are imposed by equality and equality constraints. The multiple design variables are obtained using Newton-Raphson iterations which solves huge differential equations representing Jacobi equations of UFO. The five steps to optimize the design of a doubly reinforced concrete beam proposed in this study are based on an ANN-based Lagrange optimization algorithm developed in previous studies by the authors. Optimized designs are presented and verified, showing trade-off ratios of objective functions. An example of decision-making based on the UFO is presented, which engineers and decision-makers can use as a guide throughout a preliminary design stage. GRAPHICAL ABSTRACT


Introduction and previous studies
Over the last few decades, an optimization in structural design has gained an interest of researchers and engineers. According to Mei, L. and Wang, Q. et al. (Mei and Wang 2021), structural optimization objective functions can be divided into four major categories, including cost minimization, structural performance improvement, environmental impact minimization, and multi-objective optimization. The objective functions are typically conflicting, requiring careful steering to achieve a balance. Multi-objective-based optimization challenges arise in many areas of the structural engineering field. For example, a common task based on multiple objective functions is to minimize cost while reducing the amount of carbon dioxide (CO 2 ) emissions. These objective functions are reasonably prevailing when the demands of the architect, structural engineers, construction engineers, and the government are conflicting. Choi, S.W. et al. (Choi, Oh, and Park 2017) proposed design strategies for reducing high-rise building costs and CO 2 emissions. Martinez-Martin et al. (Martinez-Martin et al. 2012) included the economic cost, embedded CO 2 emissions, and reinforcing steel congestion as three objective functions to design bridge piers with hybrid heuristic algorithms. Munk, D.J. et al. (Munk, Vio, and Steven 2015) reviewed studies of bi-objective-based optimization, including topology and shape optimization using evolutionary algorithms. Numerous researchers have investigated and developed optimization algorithms to solve multiple-objective optimization (MOO) challenges. A prevailing MOO algorithm is based on an evolutionary algorithm such as Nondominated Sorting Genetic Algorithm I and II (NSGA and NSGA II) by Deb et al. (Deb et al. 2002), Self-regulated Particle Swarm Multi-task optimization by Zheng et al. (Zheng et al. 2021) and Mohd Zain et al. (Mohd Zain et al. 2018), Multi-objective Tabu Search algorithm by D Jaeggi et al. (Jaeggi et al. 2005), MOO for diversity and performance in conceptual structural design by Brown (Brown, Tseranidis, and Mueller 2015;Brown 2016). Hu in 2018 (Zheng andHu 2018) studied multi-objective optimal design on vibration suppression of building structures with active mass damper. MOO design of truss structures has been studied by Nan (Nan, Bai, and Wu 2020) and Kaven (Kaveh and Mahdavi 2019). Afshari et al. (Afshari, Hare, and Tesfamariam 2019) reviewed and compared MOO algorithms with constraints in reinforced concrete (RC) structures. Metaheuristics-based algorithms applied to multi-objective functions optimizing a design of RC frames have been studied in Babaei and Mollayi (2021;Bekdaş and Nigdeli (2017); Kaveh and Sabzi (2012). Barraza M. and his partners (Barraza et al. 2017) investigated MOO of structural steel buildings. Some other MOO studies of RC structures have been done in Arama et al. (2020); Lee, Hong, and Choi n.d.;Park et al. (2013); Kayabekir et al. (2020). Researchers have investigated populations-based MOO algorithms in a design of RC beams (Shaqfa and Orbán 2019) in Tahmassebi et al. (2020); Zhang et al. (2021) ;Ferreira, Barros, and Barros (2003); Bekdas and Nigdeli (2013); Coello, Hernández, and Farrera (1997) ;Jahjouh, Arafa, and Alqedra (2013).
A Pareto frontier or a Pareto front is regarded as a set of the multi-objective optimized results. One of the engineers' interests is to investigate the particular tradeoff ratio estimating how much sacrifice each objective function makes. The evolutionarybased MOO algorithms such as NSGA and Particle Swarm Optimization require post-processing techniques to show the tradeoff ratio, making them known as approaches with a posterior articulation of preferences (Marler and Arora 2004). This seems unwelcome by engineers and decision-makers who are unfamiliar with optimization algorithms while having more interests in trade-off ratios among several objective functions that prevent applications of evolutionary-based MOO algorithms to real-life situations. Furthermore, the methodology of the population-based MOO algorithms is based on biological rules such as crossover and mutation in NSGA I and II that are unfamiliar to structural engineers. When a full understanding of the evolutionary-based optimization algorithms is absent, the evolutionary-based MOO studies can be difficult for engineers to apply to real-life structural design cases. Few studies have investigated the gradientbased MOO algorithms in which the first gradient of objective functions become zero at a saddle point, yielding optimized results. Xin Liu et al. (Liu and Reynolds 2016) proposed a gradient-based multiobjective-based optimization implementing a weighted sum approach.
In the present study, the authors implemented the AI-based Lagrange optimization algorithm incorporating multiple objective functions and its application in design of RC beams following the Building Code Requirements for Structural Concrete ACI 318-19 (2019). The ANN-based Lagrange algorithm based on a weighted sum approach (Yang 2014) to formulate the unified function of objective including a five-step optimization is used to solve MOO problems in this study. Lagrange multipliers are used with equality and inequality conditions based on KKT conditions to handle constrained conditions following design requirements. The unified function of objectives (UFO) is used to globalize multiple objective functions so that multiple objective functions can be optimized simultaneously to derive the Pareto frontier which cannot be found using conventional design methods. Design requirements are imposed by equality and equality constraints. The multiple design variables are obtained using Newton-Raphson iteration which solves huge differential equations representing Jacobi equations of UFO. Saddle points are, then, identified using initial trial input variables which are updated until the solutions converge. Thus, ANN-based Lagrange optimization algorithm calculates saddle points at which the first derivative of Lagrange function becomes zero. An application to a design of doubly RC beams minimizing cost, environmental impact, and beam weight is conducted based on the five-step optimization using the optimization and training toolbox provided by MATLAB (MathWorks 2020b). An example of decisionmaking using the optimized results are introduced in this study that can aid engineers for final decisions in preliminary design stage. The Pareto frontier illustrates the tradeoff ratio of investigated objective functions that are not displayed by the metaheuristics-based results. The algorithm has been initially proposed by Hong et al.  and , utilizing Artificial Neural Networks (ANN) to generalize objective functions. Lagrange multiplier method and KKT conditions are used to handle constraining conditions (Kuhn and Tucker 1951;Karush 1939). ANNs have been investigated in optimization design of RC members by Yücel et al. (Yücel et al. 2021b) and (Yücel et al. 2021a), in design of steel beams in (Ferreira et al. 2022;Nguyen, Ly, and Tran 2021;Hosseinpour, Sharifi, and Sharifi 2020;Sharifi et al. 2020).
The present study consists of the following sections. Section 2 introduces design scenarios of a doubly RC beam, while Section 3 explains a five-step optimization. This manuscript describes the ANN-based Lagrange optimizations, leading to Pareto frontier in Section 3.2. The case study of the manuscript shows the design accuracies of the Pareto frontier is validated by structural-based large datasets in Section 4. Section 5 shows a design example based on a set of optimizations minimizing several objective functions simultaneously as decision-making guidance. Finally, practical conclusions and design recommendations are drawn in Section 6.

Selection of design parameters based on design criteria of doubly RC beams
This study presents an application of ANN-based Lagrange optimization of doubly RC beams with twoend fixed boundary conditions, as shown in Figure 1, where L denotes a span length, cross-sectional dimensions are also denoted by beam width (b) and beam depth (h). ρ r;c and ρ r;t denote compressive and tensile reinforcement ratio (or rebar ratio), respectively. A beam section (b, h), and reinforcement details (ρ r;c and ρ r;t ) of a doubly RC beam must be designed to satisfy requirements in terms of strength, serviceability, and reinforcement detailing specified in the design code ACI 318-19 (2019). In this study, nominal flexural strength (ϕM n ) is required to be equal to factored moment (M u ), which is a combination of dead load (M D ) and live load (M L ) meeting strength criteria. To avoid brittle behaviour in the case of an overload, a reinforcement strain is limited at a design strength. ε rt and ε rc are reinforcement strains at a concrete strain of 0.003. Serviceability is reflected in Δ imme and Δ long , referring to immediate deflections and time-independent deflections which are limited to L=360 and L=240, respectively. It is worth noting that the cracking moment M cr , which must be less than ϕM n =1:2 to control a minimum of flexural   (2019), is investigated when calculating Δ imme . A design of a longitudinal compressive and tensile rebar in this study includes a calculation of compressive rebar ratio (ρ r;c ) and tensile rebar ratio (ρ r;t ) based on their rebar areas A r;c and A r;t , respectively. For a rectangular crosssection of doubly RC beams, a reinforcement ratio is calculated as the ratio of rebar area to cross-sectional gross area, which is equal to b � h. Rebar details such as maximum and minimum spacing are also determined, where a clear concrete cover c of 40mm, reinforcement of ϕ29, and aggregates of ϕ20 are used in this study. Table 1 shows the nine inputs for these design criteria, which include five parameters to define a beam section (RFFC A 2086786), two material properties (f 0 c ; f y ) for concrete strength and rebar yield strength, respectively, and two externally imposed moments (M D , M L ). Seven output parameters, ϕM n ; M u ; M cr ; ε rt ; ε rc ; Δ imme ; and Δ long , are calculated from nine inputs. A forward design of doubly RC beams requires determining seven outputs from nine given input parameters. Table 2 describes input and output parameters with their nomenclature.

Selection of objective functions
In this study, a design of a doubly RC beam is optimized based on three objective functions, overall construction cost and material cost (CI), beam weight (BW), and CO 2 emissions, which are selected based on engineers', contractors', and the government's interests. The overall cost or cost index (CI) of a doubly RC beam is calculated based on Korean unit cost for concrete and rebar, as defined in Table 3 (Hong et al. 2010). CO 2 emission units and the unit weight of concrete and rebar are shown in Table 3 to calculate CO 2 emissions and the weight of a doubly RC beam.

An optimization scenario
In this study, a design of a doubly RC beam in which CI; BW; and CO 2 are minimized simultaneously is performed. Beam sections (b; h; ρ r;c ; ρ r;t ) are designed for a fix-fix boundary conditions having a beam span L ¼ 9000 mm. Beam depth h is designed in a range of 400 À 1500 mm, while a beam section has an aspect ratio b=h within 0:25 À 0:8. The concrete strength and rebar yield strength used in this study are f 0 c ¼ 40 MPa and f y ¼ 550 MPa, respectively. Dead load M D and live load M L are given as 800 kN.m and 400 kN.m, respectively. The design of ductile doubly RC beams follows ACI-based requirements in terms of strength, serviceability, and rebar details, as described in Table 2. The conventional design parameters in Table 1 contain nine input parameters and ten output parameters when generating large datasets. Optimized beam parameters (b; h; ρ r;c ; ρ r;t ) are then determined on an outputside during an optimization design of a doubly RC  beam, resulting in minimized CI; BW; and CO 2 . Input parameters and 0 output parameters implemented in beam optimization will be discussed in Section 3.2.

Optimization using the optimization toolbox of Matlab
Hong et al.  described ANNbased Lagrange optimization based on a singleobjective function. The authors introduced an ANNbased Lagrange optimization and its application to cost minimization for RC beams . In the authors' previous study, an algorithm based on the weighted sum method was used to develop a unified function of objective (UFO) for RC columns (Hong et al. XXXX). This study derives a UFO for RC beams so that several objective functions can be optimized simultaneously. UFO is optimized based on MATLAB's optimization toolbox; thus, readers do not have to derive any optimization code. Engineers and users only have to provide constraint conditions of their interests. Engineers can follow five steps demonstrated in Figure 2 to optimize the unified function of objective minimizing cost, beam weight, and CO 2 emissions of doubly RC beams.
MATLAB's toolbox is used to perform UFO optimization. Section 3.2 explains the five steps in detail. A general description of a five-step optimization is introduced as follows. Initially, structural datasets are generated based on structural mechanics-based software by using randomly selected input parameters within their design ranges to obtain corresponding outputs. Structural mechanics-based software is developed based on strain-compatibility theory following a conventional design. An Artificial neural network (ANN) is established on the basis of design and training parameters such as a number of hidden layers, a number of neurons, and a training method chosen according to optimization natures. ANNs are trained by generated datasets to derive objective functions, CI, BW, and CO 2 emissions. These ANN-based objective functions can replace analytical ones, which are difficult to be derived by structural engineers analytically and explicitly because there are abounding design parameters for each type of structure. The design of RC beams should meet various requirements imposed by the design codes, architectural interests, and construction limitations. f ANN i x ð Þ, with x being an input vector based on i ¼ 1 to m is an objective function generated by ANNs, where m is a number of objective functions. Each objective function f ANN i x ð Þ is "weighted" by a fraction w i to be united in a UFO as shown in Eq. (1), where the weighted sum is simple and convenient for MOO problems. A set of optimizations in MOO problems yields a Pareto frontier. As shown in Eq. (2), the fractions of all objective functions sum up to 1. UFO is optimized using fractions ranging from 0 to 1. In this study, the FMINCON function from the MATLAB's optimization toolbox (MathWorks 2020b) is used to find saddle points for a constrained nonlinear multivariable UFO. Engineers and users provide constrained conditions and objective functions, determining optimization natures. Lagrange functions shown in Equations (3a) and (3b) (6) and Eq.
(7). Saddle points x opt are found when a convergence criterion g opt ffi 0, as shown in Eq. (8), is met.

Five steps to optimize the design of doubly RC beams based on three objective functions
This section describes five steps minimizing cost, beam weight, and CO 2 emissions simultaneously for doubly RC beams. A Pareto frontier optimizing a doubly RC beam design based on multiple-objective functions is obtained by five steps.
Step 1 describes how to establish an ANN to generate ANN-based objective functions of CI; BW, and CO 2 . A topology of an ANN used in this study is described in Figure 3.
Step 2 introduces design and an optimization scenario. Equality constraints EC x ð Þ and inequality constraints IC x ð Þ are also imposed in Step 2 based on the code requirements and other interests from architects and engineers. AI-based objective functions derived in Step 1, EC x ð Þ and inequality constraints IC x ð Þ formulated in Step 2 are used to calculate ANN-based Lagrange optimization shown in Eq. (3a) for optimization of single objective function in Step 3. The optimized design parameters b; h; ρ rt ; and ρ rc are obtained on an output-side using the Newton-Raphson iteration, described in Figure 4 shows a flow of an optimization for beam designs based on the ANN-based Lagrange optimization algorithm established in Figure 2. The optimized beam    Table 4a are input parameters in the ANN model and structural mechanicsbased software to calculate corresponding design parameters ϕM n ; M u ; M cr ; ε rt ; ε rc ; Δ imme ; Δ long ; CI; CO 2 ; and BW on an outside. The obtained parameters b; h; ρ rt ;ρ rc ; L; f 0 c ; f y ; M D ; M L ; ϕM n ; M u ; M cr ; ε rt ; ε rc ; Δ imme ; Δ long ; CI; CO 2 ; and BW are optimized design parameters, which yield a Pareto frontier.
Step 4 derives UFO based on equally spaced fractions of objective functions shown in Equations (1) and (2) to globalize multiple objective functions derived from a trained ANN as shown in Step 1. UFO is used in Eq. (3b), formulating Lagrange function of multiple objective functions. Saddle (stationary) points of Lagrange function of UFO under KKT conditions are the optimized design parameters to identify a Pareto frontier. The Pareto frontier based on UFO, which is obtained by UFO-based ANN as shown in Figure 3, is verified by structural mechanics-based software Autobeam where ϕM n ; M u ; M cr ; ε rt ; ε rc ; Δ imme ; Δ long ; CI; CO 2 ; and BW are validated.

Step 1-Establishment of an ANN
Step 1.1: Selecting random design ranges to generate large datasets Table 1 shows the nomenclatures of 19 input and output parameters generated in 100,000 datasets. Large datasets of ductile doubly RC beams are generated by randomly varying input parameters (L; b; h; ρ r;c ; ρ r;t ; f 0 c ; f y ; M D , M L ) within their design ranges described in Table  5. Corresponding outputs (ϕM n ; M u ; M cr ; ε rt ; ε rc ; Δ imme ; Δ long ; CI; BW; CO 2 ) are obtained using Autobeam, structural mechanicsbased software developed by Nguyen and Hong et al. (Nguyen and Hong 2019), which is based on a strain-compatibility-based algorithm. It is worth noting that the design criteria shown in Table 2 are included in Autobeam. To consider a wide range of beam designs, a beam length L randomly varies from 812m. Beam depth h is randomized in a range of 0:4 À 1:5m, while beam width is generated to satisfy an aspect ratio b=h of 0:25 À 0:8. Load parameters M D and M L are generated within 250 À 1500kN:m and 125 À 750kN:m, respectively. Concrete compressive strength and rebar yield strength are randomized in their ranges, such as 30 À 50MPa and 400 À 600MPa, respectively. The 100,000 datasets used to train ANNs derived in Step 1.2 are provided in Table 6, where nonnormalized datasets and datasets normalized between À 1 and 1 are shown.
Step 1.2: Selecting training parameters Training parameters implemented in training ANNs include the following: • A number of hidden layers • A number of neurons • A type of activation functions applied to all layers including the output layer • A number of required epochs and number of validated epochs • A training method, such as training on entire datasets (TED), parallel training method (PTM), chained training scheme with revised sequence (CRS) (Hong 2020).
Training parameters are selected based on the complexity of large datasets and the availability of computation hardware. An ANN in this study uses three, four, and five hidden layers, each of which contains 30, 40, and 50 neurons. Table 7 shows the training accuracies. The activation functions are tanh at hidden layers and linear at an output layer. PTM is selected in which nine input parameters are mapped to each of the 10 outputs separately. Figure 3 shows the topology of an ANN. MATLAB training toolbox (MathWorks 2020b) is used to train ANNs.

Step 2-Defining MOO problems
Step 2.1: Deriving objective functions This study investigates a design of doubly RC beams having an optimization scenario described in Section 2.3. Functions of three objectives CI; BW , and CO 2 to be minimized are derived from an ANN denoted as f ANN  Table 4a and four input variables (b; h; ρ r;c , and ρ r;t ), which are to be determined based on the Newton-Raphson iteration. g N x ð Þ is used to normalize inputs so that it has a dimension of 9 � 1. These normalized inputs are then interconnected to hidden layers through interior weight and bias matrices before being sent to an output layer. ω 1 ð Þ CI is a weight matrix to transfer an input parameters x to the first hidden layer, making its dimension 30 � 9, representing 30 neurons and 9 inputs.  Table 7.
Step 2.2: Formulation of equality and inequality constraints Six equalities and eleven inequalities are formulated in Table 8 to reflect constrained conditions based on design criteria and engineers' needs.

Step 3-Optimization based on a singleobjective function
This step is required to maximize and minimize each objective function to define their boundaries for normalization because each normalized objective function is multiplied by a fraction to derive UFO. Optimization based on a single-objective function is based on the FMINCON function in MATLAB optimization toolbox, which was introduced in Section 3.1, including equalities and inequalities formulated in Table 8. Equations (12) and (14) (9) and (11).   . Optimized design parameters on a Pareto frontier illustrated in Figure 6.
Single-objective functions f ANN CI x ð Þ; f ANN BW x ð Þ; and f ANN CO 2 x ð Þ are normalized from their corresponding minimum to maximum provided in Table 9, as shown in Eq. (15) to (17).

Step 4-Formulation of UFO
(a) Defining fractions Fractions of three objective functions w CI ; w CO 2 ; andw BW are generated based on the linspace function in MATLAB, which is a function of linear spacing to assure that they are equally distributed between 0 and 1. Eq. (18) defines 19 fractions of CI w CI ð Þ that are equally spaced between 0 and 1 so that their spacing is equal to 1=18. In Eq. (19), 19 fractions of CO 2 w CO 2 ð Þ are, then, equally spaced between 0 and w CI obtained in Eq. (18). Fractions of BW w BW ð Þ are calculated with respect to w CI and w CO 2 using Eq. (20). A total of 361 fractions of three objective functions are generated based on the formulas presented in Table 10 and Figure 4. A number of fractions is adjustable considering optimization natures.  Table 10.

Step 5-Optimization of UFO
Substituting equalities, inequalities shown in Table 8 and UFO having 361 fractions shown in Eq. (21) into Eq.
(3) provides a Lagrange function for UFO, which minimizes designs of a doubly RC beam using the FMINCON function of MATLAB. Figure 6 shows a Pareto frontier using 361 fractions optimizing designs of doubly RC beams based on the optimization scenario presented in Section 2.3. Each point on the Pareto frontier corresponding to a fraction in Figure 5 uncovers design parameters when optimizing doubly RC beams, providing engineers and decision-makers with particular trade-off ratios contributed by each of the three objective functions. An optimization accuracy is estimated based on two criteria:

Discussions on the Pareto frontier
(1) Calculation accuracies obtained using Eq. (22) are defined to indicate computational errors.
Output ANN and Output Check are obtained output parameters ϕM n ; M u ; M cr ; ε rt ; ε rc ; Δ imme ; Δ long ; CI; CO 2 ; and BW using optimized parameters b; h; ρ rt , ρ rc and prescribed parameters L; f 0 c ; f y ; M D ; and M L based on the ANN model and structural mechanics-based software Autobeam. In Figure 5(b) and Table 11, six optimized designs based on six fractions on a Pareto frontier are verified. Design points P1, P2, and P3 minimize the UFO of the multi-objective functions CI; CO 2 , and BW based  Table 4(a). As shown in Figure 7, the Pareto frontier representing minimized design parameters of a doubly RC beam predicts the lower boundary of one million structural datasets, verifying the optimization accuracies of the proposed AI-based Lagrange optimization.

Decision-making based on a Pareto frontier
As shown in Figure 6, a Pareto frontier provides optimized design parameters by simultaneously minimizing three objective functions (cost CI; CO 2 , and beam weight BW), while a particular trade-off among the three objective functions is represented  Table 13. Design parameters of two designs marked in blue color in Table 12.
by a fraction w CI : w CO 2 : w BW . The optimized design parameters and their corresponding fractions can significantly aid engineers and decision-makers in a preliminary design stage, offering them optimal designs minimizing CI; CO 2 , and BW simultaneously. This section's optimization example presents a design recommendation for a doubly RC beam. When the budget is assumed between 82,000 and 85,000 KRW/m, a doubly RC beam with constraint conditions shown in Table 4 can be designed. When a design yields the least weight and CO 2 emissions at the same time, the ranges of CO 2 emissions and beam weight corresponding to the budget between 82,000 and 85,000 KRW/m can be estimated as 0.181-0.187 t-CO 2 /m and 6.079-6.311 kN/m, respectively, as shown in Figure 8. The green dots in Figure 8 denote a Pareto based on 361 fractions obtained in Figure 5 and Table 10. A Pareto frontier shown in Figure 8 illustrates an estimated least design range of CO 2 and BW when cost ranges from 82,000 to 85,000 KRW/m. Figure 8 helps determine a design meeting the least BW of 6.079 kN/m in the range (6.079-6.311 kN/m) that corresponds to the minimized CI= 84,662.2 KRW/m and CO 2 = 0.187 t-CO 2 /m when minimizing beam weight is preferable. A fraction of the design shown in Figure 8 is 0.11:0.49:0.40, indicating that minimizing CI (CI= 84,662.2 KRW/m) contributes 11% to the optimization based on multi-objective functions, UFO, while contributions of minimizing of CO 2 (CO 2 = 0.187 t-CO 2 /m) and BW (BW ¼ 6:079 kN/m) are 49% and 40% to the unified optimization, respectively. As shown in Table 12, engineers without any computational limitations can choose one of design diverse options based on the trade-off ratios or fractions among the three objective functions corresponding to a project preference, allowing them to calculate minimized CO 2 and BW for a pre-determined budget CI and vice versa. For example, engineers can choose the design having the fractions of three objective functions w CI : w CO 2 : w BW of 0.28:0.36:0.36 and 0.55:0.10:0.35 as shown in Table 13. An optimized design can now be selected based on multi-objective functions, assisting engineers to perform efficient designs.

Conclusions and contributions
This study proposed an optimization of multi-objective functions by which engineers and policymakers can make decisions. To enable simultaneous optimization, ANN-trained objective functions are globalized into one function UFO. A five-step optimization is described in detail and can be applied to help engineers design diverse structural systems. Optimization examples based on the multi-objective functions will assist engineers and decision-makers in the preliminary design stage. The following are some conclusions drawn from this study.
(1) This study formulates the unified functions of objectives (UFO) where multiple objective functions are globalized into one function. Saddle points of the Lagrange function of the UFO are found where multiple objective functions are optimized simultaneously. Engineers, contractors, and government can simply establish objective functions according to their own interests.
(2) A Pareto frontier is obtained from multiobjective functions, UFO, and confirmed by yielding a lower boundary which coincides with one obtained from a large structural dataset, ensuring the proposed optimization's reliability. The accuracy of optimized designs provided by a five-step optimization based on ANN-based Lagrange optimization is verified by computation errors as large as 5.5%.
(3) Design requirements imposed by structural codes are reflected by equality and equality constraints while RC beams are optimized based on unified functions of objectives (UFO). The multiple design variables are now obtained using Newton-Raphson iteration which solves huge differential equations representing Jacobi equations of UFO.
(5) It is now possible to find holistic design parameters optimizing multi-objective functions with identifying fractions contributed by each objective function. The proposed ANN-based Pareto frontier offers the trade-off ratios that are not well presented when using the metaheuristics-based MOO algorithm.
(6) The optimized design parameters and their corresponding fractions can significantly aid engineers and decision-makers in a preliminary design stage, offering them optimal designs minimizing CI; CO 2 , and BW simultaneously. An optimized design can now be selected based on multiobjective functions that can be implemented in achieving efficient designs.

Notes on contributors
Dr. Won-Kee Hong is a Professor of Architectural Engineering at Kyung Hee University. Dr. Hong received his Master's and Ph.D. degrees from UCLA, and he worked for Englelkirk and Hart, Inc. (USA), Nihhon Sekkei (Japan) and Samsung Engineering and Construction Company (Korea) before joining Kyung Hee University (Korea). He also has professional engineering licenses from both Korea and the USA. Dr. Hong has more than 30 years of professional experience in structural engineering. His research interests include a new approach to construction technologies based on value engineering with hybrid composite structures. He has provided many useful solutions to issues in current structural design and construction technologies as a result of his research that combines structural engineering with construction technologies. He is the author of numerous papers and patents both in Korea and the USA. Currently, Dr. Hong is developing new connections that can be used with various types of frames including hybrid steel-concrete precast composite frames, precast frames and steel frames. These connections would help enable the modular construction of heavy plant structures and buildings. He recently published a book titled as "Hybrid Composite Precast Systems: Numerical Investigation to Construction" (Elsevier).
Thuc Anh Le is a Master's student at Structural Engineering/ Construction Technology at Kyung Hee University, Korea.