Discrete Optimization with Fuzzy Constraints

The primary benefit of fuzzy systems theory is to approximate system behavior where analytic functions or numerical relations do not exist. In this paper, heuristic fuzzy rules were used with the intention of improving the performance of optimization models, introducing experiential rules acquired from experts and utilizing recommendations. The aim of this paper was to define soft constraints using an adaptive network-based fuzzy inference system (ANFIS). This newly-developed soft constraint was applied to discrete optimization for obtaining optimal solutions. Even though the computational model is based on advanced computational technologies including fuzzy logic, neural networks and discrete optimization, it can be used to solve real-world problems of great interest for design engineers. The proposed computational model was used to find the minimum weight solutions for simply-supported laterally-restrained beams.


Introduction
The theory of fuzzy sets can be used to model imprecision, ambiguity or fuzziness in the formulation of structural optimization problems.In the formulation of such problems, a major source of imprecision, or fuzziness, occurs in the evaluation of constraints.In traditional optimization algorithms, constraints are satisfied with a tolerance defined by a crisp or non-fuzzy number.In reality, in common engineering practice, this evaluation involves many sources of approximations [1].A design of structure is considered satisfactory when the several constraints are satisfied within a given predetermined tolerance.However, when an optimization algorithm satisfies the constraints precisely (for a defined small tolerance degree of numerical computations), it can miss the true optimum design within the confines of practical and realistic approximations.
Adeli [2] demonstrated that by taking into account the fuzziness and imprecision in the constraints and employing fuzzy set theory, it is possible to reduce the objective function further and substantially increase the probability of finding the actual global optimum solution.The goal of the present research, carried out by several authors, was to model the effects of fuzziness in the formulation of a genetic algorithm (GA)-based structural design optimization problem [3][4][5][6].Another objective of the authors' research was to improve the convergence and efficiency of GAs through the use of fuzzy set theory [7][8][9].Several articles have been published on the fuzzy optimization of structures [10,11], with the objective of reducing the number of iterations and the total computer processing time needed.
Uncertainty exists in almost every real-world problem.In general, uncertainty is inseparable from measurement.It emerges from a combination of the limits of measurement with instruments and unavoidable errors in measurement.In this paper, the fuzziness was considered as part of the constraints.The constraints were developed using the neuro-fuzzy technique, which was based on past experience, recommendations and measurements.Fuzzy constraints were then used in discrete optimization along with other crisp constraints.This allowed us to include experience, recommendations and experimental measurements in the optimization problem.Optimization using non-linear programming (NLP) and fuzzy constraints has been done by Jelusic [12], however without the discrete optimization approach.While NLP deals with the continuous optimization of structures, mixed-integer linear programming (MILP) performs continuous-discrete optimization, where the structural topology, discrete materials (steel) and standard dimensions (steel sections) are known.
In order to find the minimum weight solutions for simply-supported laterally-restrained beams, the appropriate deflection limit should be specified.The comparison of different design codes showed that the deflection limits are too liberal.This paper defines the soft constraint for the deflection limit based on experiential rules acquired from experts and utilizing recommendations.This newly-developed soft constraint obtained with an adaptive network-based fuzzy inference system (ANFIS) is then used in the optimization model.ANFIS can learn from examples and is fault tolerant in the sense that it is able to handle noisy or incomplete data.The expert evaluations for the deflection limit are very subjective; therefore, the data for approximation function are expected to be vague, imprecise, incomplete or even contradictory.Additionally, the proposed ANFIS techniques include fuzzy clustering (FCM), which searches for patterns in data points.The study suggests that deflection limits could be reconsidered in the future by the experts who have a prolonged or intense experience through practice.

Structural Optimization and Fuzzy Set Theory
A crisp non-fuzzy structural optimization is formulated as follows: find the vector of the design variables x such that the objective function F(x) is minimized subject to the equality and inequality constraints: where N ec is the number of equality constraints and N iec is the number of inequality constraints; g u i (x) is the upper bound on the constraint g i (x), and g l i (x) is the lower bound on the constraint g i (x).If vagueness is considered in the objective function and constraints, then the variables (x) can be obtained from a fuzzy decision D, such that the membership function µ D for the fuzzy decision D can be obtained from the intersection of the fuzzy membership functions for the objective function and constraints; see Equation (3): where µ F (x) is the membership functions for the objective function and µ g i (x) is the membership functions for the i-th inequality design constraint.From this fuzzy decision, the optimum solution (x * ) for the variable x can be obtained by using the max-min procedure [13]; see Equation ( 4): where: The max-min procedure can be solved by maximizing a scalar parameter λ (overall satisfaction parameter) [14]; see Equations ( 6)-( 9 where µ u g i (x) and µ l g i (x) are the membership functions for the upper and lower bounds of the inequality constraints µ g i (x) (Equation ( 2)), respectively.The equality constraints, in Equation (1), are not included in the fuzzy formulations because they have to be satisfied strictly.

ANFIS Architecture for the Development of Soft Constraint Functions
For a Sugeno fuzzy model, a rule set with n fuzzy "if-then" is as follows: Rule 1: If x is A 1 and y is B 1 , then: Rule i: If x is A i and y is B i , then: Rule n: If x is A n and y is B n , then: where a 1 0 , a 1 1 , a 1 2 , a i 0 , a i 1 , a i 2 , a n 0 , a n 1 , a n 2 are consequent parameters and x and y are input variables.The output of each rule is equal to the constant, and the final output is the weighted average of each rule's output.
The weights are obtained from a Gaussian membership function.
where c is the position of the center of the curve's peak and σ is the width of the curve.Parameters c and σ are premise parameters.The first membership grade of the fuzzy set (A i , B i , C i , D i ) is calculated with the following equations: where x and y are the input variables in the Gaussian membership function.After this, the product of the membership function for every rule is calculated: where w i represents the fire strength of the rule i.The ratio of the i-th rule's firing strength to the sum of all of the rule's firing strengths is defined with: In order to achieve the desired input-output mapping, the consequent and premise parameters need to be updated according to the given training data and the hybrid learning procedure.This hybrid learning procedure [15] is composed of a forward pass and backward pass.In the forward pass, the algorithm uses the least-squares method to identify the consequent parameters.In the backward pass, the errors are propagated backward, and the premise parameters are updated by gradient descent.

Discrete Optimization
Exhaustive enumeration (EE) is the simplest of the discrete optimization techniques.It evaluates an optimum solution for all combinations of the discrete variables.The best solution is obtained by scanning the list of all feasible solutions for the minimum value.The total number of evaluations is: where n d is the number of discrete variables and p i is the pre-established set of discrete values.If either n d or p i (or both) are large, it shows that much work will be required.It also shows an exponential growth in the calculations with the number of discrete variables.In a mixed optimization problem, this would involve the continuous optimum solution of a reduced model.It is not a serious problem if the mathematical model and its computer calculations are easy to implement.If the mathematical model requires extensive calculations, then some concerns may arise.Programming exhaustive enumeration is straightforward.The processing speed, large available desktop-memory and easy programming, through software like MATLAB [16], make exhaustive enumeration a very good idea today.This program is also ideal because the solution is a global optimum.The most important step is translating the mathematical model into a program code.
The number of design variables in the model is reduced by the number of discrete variables.Model reduction is involved in enumeration techniques.
The algorithm with feasibility requirements is as follows: Step 1. s * = in f , K = [0, 0, . . . , 0]For every allowable combination of (y 1 , y 2 , . . . , This application example presents how soft constraints are included in discrete optimization.ANFIS is used to integrate recommendations of deflection limits into an optimization model.

Example Design of a Simply-Supported Laterally-Restrained Beam Application
The basic design process is formed by determining the design loads acting on the structure, determining the design loads on individual elements and calculating the bending moments, shear forces and deflections of the beams.Generally, laterally-restrained beams should be checked for their ultimate limit state (ULS) and serviceability limit state (SLS).As this article is about the design of steel structural elements, the following were examined: (1) resistance of the cross-section to bending (ULS), (2) resistance to shear buckling (ULS), (3) resistance to flange-induced buckling (ULS), (4) resistance of the web to transverse forces (ULS) and ( 5) deflection (SLS).
The beam is loaded by a uniformly-distributed dead load g k and a uniformly-distributed imposed load q k , as shown in Figure 1.
The beam is loaded by a uniformly-distributed dead load gk and a uniformly-distributed imposed load qk, as shown in Figure 1.The civil engineer must evaluate the possible future levels of loading (self-weight, snow, wind), which the structure may be subjected to during its design life.Then, using hand calculations or computer methods, the loads acting on individual construction elements can be evaluated.The loads are used to calculate the shear forces, bending moments and deflections at critical sections along the construction elements.Finally, sufficient dimensions for the construction element can be defined.

Design Loads
The loads acting on a beam are divided into imposed and dead loads.For each type of loading there will be characteristic values and design values that must be estimated.In addition to this, the designer will have to determine the particular combination of loading that is likely to produce the most adverse effect on the structure in terms of shear forces, bending moments and deflections.The design loads are obtained by multiplying the characteristic loads by the partial safety factor for loads.However, before flexural members, such as beams, can be sized, the design bending moments and shear forces must be evaluated.Design shear forces and moments in beams are calculated using standard equations (Equations ( 21) and ( 22)).Having calculated the shear force and design bending moment, all that now remains to be done is to estimate the dimensions and strength of the beam required.Yield strength and section classification are used for the initial choice of the section.If the section is thick, i.e., has thick flanges and web, it can sustain the formation of a plastic hinge.On the other hand, a slender section, for example with thin flanges and web, will fail by local buckling before the yield stress can be reached.
Design action: Design bending moment: Design shear force: Strength classification: The civil engineer must evaluate the possible future levels of loading (self-weight, snow, wind), which the structure may be subjected to during its design life.Then, using hand calculations or computer methods, the loads acting on individual construction elements can be evaluated.The loads are used to calculate the shear forces, bending moments and deflections at critical sections along the construction elements.Finally, sufficient dimensions for the construction element can be defined.

Design Loads
The loads acting on a beam are divided into imposed and dead loads.For each type of loading there will be characteristic values and design values that must be estimated.In addition to this, the designer will have to determine the particular combination of loading that is likely to produce the most adverse effect on the structure in terms of shear forces, bending moments and deflections.The design loads are obtained by multiplying the characteristic loads by the partial safety factor for loads.However, before flexural members, such as beams, can be sized, the design bending moments and shear forces must be evaluated.Design shear forces and moments in beams are calculated using standard equations (Equations ( 21) and ( 22)).Having calculated the shear force and design bending moment, all that now remains to be done is to estimate the dimensions and strength of the beam required.Yield strength and section classification are used for the initial choice of the section.If the section is thick, i.e., has thick flanges and web, it can sustain the formation of a plastic hinge.On the other hand, a slender section, for example with thin flanges and web, will fail by local buckling before the yield stress can be reached.
Design action: Design bending moment: Design shear force: Strength classification: Section classification: where: and:

Resistance of Steel Cross-Sections
The structural design of steel beams primarily involves predicting the strength of their members.This requires the designer to imagine all of the ways in which each member may fail during its design life.Common modes of failure associated with beams are local buckling, shear, shear buckling, web bearing and buckling, lateral torsional buckling, bending and deflection.

Bending Moment
When shear force is absent or of a low value, the design value of the bending moment, M Ed , at each section should satisfy the following: where M c,Rd is the design resistance for bending around one principal axis, taken as follows: (a) the plastic design resistance moment of the gross section: where W pl is the plastic section modulus, for Class 1 and 2 sections only; (b) the elastic design resistance moment of the gross section: where W el,min is the minimum elastic section modulus for Class 3 sections; (c) the local buckling design resistance moment of the gross section: where W eff,min is the minimum effective section modulus for Class 4 cross-sections only; (d) the ultimate design resistance moment of the net section at bolt holes M u,Rd , if this is less than the appropriate values above.In calculating this value, fastener holes in the compression zone do not need to be considered; they would need to be if they were oversized, slotted or filled by fasteners.In the tension zone, holes do not need to be considered, provided that: The design value of the shear force V Ed at each cross-section should satisfy the following: where V c,Rd is the design shear resistance.For the plastic design, V c,Rd is taken as the design plastic shear resistance, V pl,Rd , given by: where A v is the shear area, which, for the rolled I and H sections, loaded parallel to the web, is: where: b overall breadth r root radius t f flange thickness t w web thickness h w depth of the web η conservatively taken as 1.0 A cross-sectional area Fastener holes in the web do not have to be considered in shear verification.Shear buckling resistance for unstiffened webs must additionally be considered when: For a stiffened web, shear buckling resistance will need to be considered when: where k τ is the buckling factor for shear and is given by: for a/h w < The plastic resistance moment of the section is reduced by the presence of shear force.When the design value of the shear force, V Ed , exceeds 50 percent of the plastic shear design resistance, V pl,Rd , the design resistance moment of the section, M v,Rd , should be calculated using a reduced yield strength taken as: Thus, for the rolled I and H sections, the reduced design resistance moment for the section around the major axis, M y,v,Rd , will be given by:

Shear Buckling Resistance
As noted above, the shear buckling resistance of unstiffened beam webs has to be checked when: The value of 1.0 for η for all steel grades up to and including S460 is recommended.For standard rolled beams and columns, this check is rarely necessary.However, as h w /t w is usually less than 72ε, it was not discussed in this section.

Flange-Induced Buckling
To prevent the possibility of the compression flange buckling in the plane of the web, Eurocode 3-5 [17] requires that the ratio h w /t w of the web should satisfy the following criterion: where: For loading Types (a) and (c), the web is likely to fail as a result of: (i) crushing of the web close to the flange accompanied by yielding of the flange; the combined effect is sometimes referred to as web crushing (ii) localized buckling and crushing of the web beneath the flange; the combined effect is sometimes referred to as web crippling.
For loading Type (b) the web is likely to fail as a result of: (i) web crushing (ii) buckling of the web over most of the depth of the member.
Provided that the compression flange is sufficiently restrained in the lateral direction, the design resistance of webs of beams under transverse forces can be determined in accordance with the recommendations in Eurocode 3 [17].
In Eurocode 3 [17], it is stated that the design resistance of webs to local buckling is given by: where: f yw is the yield strength of the web t w is the thickness of the web γ M1 is the partial safety factor = 1.0 L eff is the effective length of the web that resists transverse forces = χ F l y , in which χ F is the reduction factor due to local buckling.l y is the effective loaded length, appropriate to the length of the stiff bearing s s .As stated in Clause 6.3 of Eurocode 3-5 [17], s s should be taken as the distance over which the applied load is effectively distributed at a slope of 1:1, but s s ≤ h w .
Reduction factor, χ F : The reduction factor χ F is given by: where: in which: Effective load length, l y : As stated in Clause 6.5 [17] for loading Types (a) and (b), the effective load length, l y , is given by: where: and if: or if: For loading Type (c), l y is taken as the smallest value obtained from Equations ( 52) and (53), as follows: where:

Deflections
Several vertical deflections are defined in the Eurocode [18].However, the National Annex to Eurocode 3 [17] recommends that verification of vertical deflections, δ, under unfactored imposed loads should be carried out.The designer is responsible for specifying appropriate limits of vertical deflections, which should be agreed upon with the client.However, like British Standard (BS) 5950 [19], the National Annex to Eurocode 3 [17] also recommends that verifications be made on vertical deflections, δ, under unfactored imposed loads.It suggests that in the absence of other limits, the recommendations in the Eurocode may be used.
The recommendations were examined and used for the building of the model and to predict the limits of vertical deflection.Two parameters with the biggest influence on the vertical deflection limits were considered.The influence of each parameter was determined on the basis of recommendations and engineering judgment.The comparison of different design codes (Eurocode, American Institute of Timber Construction (AITC) [20]) showed, that the deflection limits are very different.The study also suggests that deflection limits could be reconsidered in the future by the designers who have a prolonged or intense experience through practice.

ANFIS for the Development of the Constraint Function
A model to limit the vertical deflection based on recommendations was developed.The ANFIS model has two inputs: applied live load LL (kN/m 2 ) and classification CLASS (-); and it has one output: LIMIT.The ANFIS-LIMIT model was proposed in order to calculate the deflection limits.MATLAB [16] and a Fuzzy Logic Toolbox were used as an interface for mathematical modeling and data handling.
One of the most important stages in the ANFIS technique is the collection of data.The training data were chosen based on the recommendations in AITC [20], Eurocode (Table 1), professional experience and past projects (Table 2).The classification used is separated into three groups.The first group is reserved for railway bridge stringers and beams used for commercial and institutional buildings with plaster ceilings.The second group is reserved for highway bridge stringers and beams used for commercial and institutional buildings without plaster ceilings.The third group is reserved for industrial roof beams.The applied load (LL) and the classification system (CLASS) were taken as input parameters; whereas, the deflection limit (LIMIT) was considered as an output parameter.The training dataset (see Table 2) can be improved by adding additional recommendations and more past experience.These values can be assigned to the parameters and deflection limit.In this model, 33 evaluations were defined for a different applied live load and classification.
For the Sugeno fuzzy model [21], a rule set with i, i I, I = {1, 2} and fuzzy "if-then" rules were defined by Equations ( 55) and (56): Rule 1: If LL is A 1 and CLASS is B 1 , then: Rule 2: If LL is A 2 and CLASS is B 2 , then: where a 1 0 , a 1 1 , a 1 2 , a 2 0 , a 2 1 , a 2 2 are consequent parameters and LL and CLASS are input variables.The calculation procedure of the ANFIS models is as follows: 1.
the membership grade of the fuzzy set the product of membership function for each rule is calculated; 3.
the ratio between the i-th rule's firing strength and the sum of all rules' firing strengths is calculated; 4.
the output of each rule is calculated; and 5.
the weighted average of each rule's output is calculated.
The first membership grade of the fuzzy set (A i , B i , C i , D i ) is calculated with Equations ( 57) and (58): where LL and CLASS are inputs to Gaussian membership functions, and the parameters c Ai , c Bi , σ Ai , σ Bi are premise parameters.In addition to this, the products between the membership functions for every rule are calculated; see Equations ( 59) and (60): where w 1 and w 2 represent the firing strength of the each rule.The weighted average of each rules' output is defined as the ratio between the i-th rule's firing strength and the sum of all of the rule's firing strengths; see Equation (61): The output of each rule is finally determined as the sum of products between the weighted average of each rule's output and the linear combination between input variables and consequent parameters: For the model, the following values were evaluated: premise parameters, consequent parameters, firing strengths and weighted averages of rules outputs.
The structure of the model is shown in Figure 2.While the nodes on the left side represent the input data, the right node stands for the output.The model includes two inputs, the applied live load LL (kN/m 2 ) and the classification CLASS (-), as well as a single output deflection limit LIMIT (-).
In a conventional fuzzy inference system, the number of rules is decided by the researcher/engineer who is familiar with the system to be modeled.There are no simple ways of determining in advance the minimum number of membership functions to achieve a desired performance level.In the present attempt, the number of membership functions assigned to each input variable was chosen empirically by examining the desired input-output data and by trial and error.For the deflection limit model, two membership functions were chosen for each input.Figure 2 shows the membership functions for LL and CLASS for the deflection model LIMIT.Note that all of the membership functions used were Gaussian membership functions, defined by Equations ( 57) and (58).determining in advance the minimum number of membership functions to achieve a desired performance level.In the present attempt, the number of membership functions assigned to each input variable was chosen empirically by examining the desired input-output data and by trial and error.For the deflection limit model, two membership functions were chosen for each input.
Figure 2 shows the membership functions for LL and CLASS for the deflection model LIMIT.Note that all of the membership functions used were Gaussian membership functions, defined by Equations ( 57) and (58).After the numbers of membership functions associated with each input were fixed, the initial values of premise parameters were set in such a way that membership functions were equally spaced along the operating range of each input variable.The LIMIT model contained two rules with two membership functions being assigned to each input variable.The total number of fitting parameters was 14, composed of eight premise parameters and six consequent parameters.These parameters were obtained by using a hybrid algorithm.MATLAB [16] was used as an interface for mathematical modeling.Premise and consequent parameters are presented in Table 3.

Fuzzy Optimization Model: Beam Implementation
In accordance with the exhaustive enumeration (EE) problem formulation, an EE optimization model for the optimization of a simply-supported beam (BEAM) was developed.Since the model BEAM is proposed to be used in the design of steel elements, all decisive design constraints were involved in the model.The model enables optimization of the system for various spans, loads and different material properties.For mathematical modeling and data input/output, a high level language, MATLAB [16], was used.The proposed optimization model includes input data, variables and the BEAM system's weight objective function, which is subjected to the structure's crisp and soft constraints; see Appendix A.
Choosing an I-beam from the list also identifies all of the design variables.This then becomes a After the numbers of membership functions associated with each input were fixed, the initial values of premise parameters were set in such a way that membership functions were equally spaced along the operating range of each input variable.The LIMIT model contained two rules with two membership functions being assigned to each input variable.The total number of fitting parameters was 14, composed of eight premise parameters and six consequent parameters.These parameters were obtained by using a hybrid algorithm.MATLAB [16] was used as an interface for mathematical modeling.Premise and consequent parameters are presented in Table 3.

Fuzzy Optimization Model: Beam Implementation
In accordance with the exhaustive enumeration (EE) problem formulation, an EE optimization model for the optimization of a simply-supported beam (BEAM) was developed.Since the model BEAM is proposed to be used in the design of steel elements, all decisive design constraints were involved in the model.The model enables optimization of the system for various spans, loads and different material properties.For mathematical modeling and data input/output, a high level language, MATLAB [16], was used.The proposed optimization model includes input data, variables and the BEAM system's weight objective function, which is subjected to the structure's crisp and soft constraints; see Appendix A.
Choosing an I-beam from the list also identifies all of the design variables.This then becomes a single variable problem, the variable being the particular item from the list of beams.
The input data (constants) represent various design data for the optimization, i.e., constants/coefficients, which are involved in the objective function and (in) equality constraints.The design data are comprised of the span length L (m), the characteristic dead g k (kN/m) and imposed q k (kN/m) loads, the yield strength of the steel f y (MPa), the modulus of elasticity E (MPa), the density of the steel gam (kg/m 3 ), bearings width ss (mm) and the allowable deflection of the beam lim (-).In addition to this, the data also include the coefficients involved in the design inequality constraints: safety factor for dead loads SFg (-), safety factor for imposed loads SFq (-), partial factor for resistance of cross-sections SFm0, partial factor for resistance of members to instability SFm1, modification factor k (-), non-dimensional slenderness lamflim (-) and the reduction factor for the relevant buckling curve ksiflim (-).
The objective variable f defines the weight of the steel beam.The aim of optimization is to find a steel beam with the minimum weight that satisfies all of the design constraints.
Design constraints that enable the ULS and SLS are satisfied by the following conditions: • Condition 1, resistance of the cross-section to bending (ULS): verified by Equation (28), by which the design bending moment M Ed (kNm) must not exceed the bending moment resistance M Rd (kNm).

•
Condition 2, resistance of the cross-section to shear (ULS): verified by Equation (33), by which the design shear force V Ed (kN) must not exceed the shear resistance V Rd (kNm).

•
Condition 3, deflection (SLS) is considered: the calculated vertical deflection of the steel beam must be less than specified by the ANFIS-LIMIT model.This is a single-variable problem, the variable being the particular item from the list of beams.In this case, there are no necessary geometrical constraints or side constraints.A complete enumeration can be performed on the selected beams from the stock list (see Appendix B), and the best beam can easily be identified.
In this example, the beam was used as a horizontal member.The objective was to design a minimum mass beam that would not fail, according to recommendations, under bending, shear and specified deflection.The length L of the beam was specified as 25 m.The steel beam was loaded by uniformly-distributed loading g k = 5 kN/m and q k = 15 kN/m, as shown below.Steel was chosen from the material of the beam.The modulus of elasticity of the steel was E = 210 GPa.The weight of the steel was ρ = 7850 kg/m 3 .The yield strength in tension was f y = 335 MPa.The specified applied load LL was 3 kN/m 2 and used a classification of 1.A safety factor of 1.35 on the permanent load and 1.50 on the variable load were assumed.The results of an exhaustive enumeration computer code showed that the optimal section is HE 1000 × 393, with a weight of 9816.4 kg.

Conclusions
The article presents how recommendations can be implemented in a discrete optimization model.Good engineering judgment should be integrated into construction design; therefore, a mathematical model was developed based on engineering judgment and past experience.For this purpose, the theory of fuzzy sets was used.The advantages of the proposed fuzzy algorithm ANFIS are acknowledged and incorporated into the model based on the imprecision and fuzziness in the code-based design constraints.The comparison of different design codes showed that the deflection limits are very different and too liberal.The expert evaluations for the deflection limit are subjective; therefore, the data obtained from experts are expected to be vague, imprecise, incomplete or even contradictory.Additionally, the proposed ANFIS techniques include fuzzy clustering (FCM), which searches for patterns in data points [22].The study also suggests that deflection limits could be reconsidered in the future by the experts who have a prolonged or intense experience through practice [23].The advantages of using the exhaustive enumeration are reduced optimum weight values and obtaining a solution that is at a global optimum.The proposed computational model was used to find minimum weight solutions for simply-supported laterally-restrained beams.For selected design variables, an optimum steel section was found based on steel sections found on the market, according to the European beams tables.The computational model is based on advanced computational technologies, including fuzzy logic, neural networks and discrete optimization.It was developed to solve real-world problems that are of great interest to design engineers.

Figure 1 .
Figure 1.Simply-supported beam and steel section.

Figure 1 .
Figure 1.Simply-supported beam and steel section.
is the area of the compression flange = b•t f f yf is the yield strength of the compression flange The factor k assumes the following values: plastic rotation utilized, i.e., Class 1 flanges: 0.3; plastic moment resistance utilized, i.e., Class 2 flanges: 0.4; elastic moment resistance utilized, i.e., Class 3 or Class 4 flanges: 0.55.4.2.6.Resistance of the Web to Transverse Forces Eurocode 3-5 [17] categorize between two types of forces applied through a flange to the web: (a) forces resisted by shear in the web (loading Types (a) and (c)).(b) forces transferred through the web directly to the other flange (loading Type (b)).

Figure 2 .
Figure 2. Fuzzy inference system implemented within the framework of adaptive networks ANFIS-LIMIT.

Figure 2 .
Figure 2. Fuzzy inference system implemented within the framework of adaptive networks ANFIS-LIMIT.

• Condition 4 ,
resistance to flange-induced buckling (ULS): to prevent the possibility of the compression flange buckling in the plane of the web.• Condition 5, Condition 6, Condition 7 and Condition 8, resistance of the web to transverse forces (ULS): to prevent the possibility of the local buckling of webs.

Table 1 .
Deflection limit according to the recommendations.

Table 2 .
Training data for the ANFIS-LIMIT model.
* 1, Railway bridge stringers, beams used for commercial and institutional buildings with plaster ceiling; 2, highway bridge stringers and beams used for commercial and institutional buildings without plaster ceiling; 3, industrial roof beams.

Table 3 .
Premise and consequent parameters of the ANFIS-LIMIT model.

Table 3 .
Premise and consequent parameters of the ANFIS-LIMIT model.

Table A1 .
Dimensions and properties of steel beams (European beams).