Automatic Position Information of Web-Openings of Building Using Minimized Strain Energy Topology Optimization

This study presents a new engineering practice and idea that material topology optimization results may be utilized to optimally decide the positions of web-openings of structuralmembers in a building structure.Material topology optimization utilizes element densities as design parameters, that is, nominal constructionalmaterial, and then optimalmaterial distributions of densities between voids (0) and solids (1) in a given design domain represent the determination of topology and shape. That means that regions with element density values become occupied by solids in a design domain, while there are only void phases in regions where no density values exist. Therefore, the void regions of topology optimization results may provide design information that decides appropriate depositions of web-opening in structure. Numerical examples demonstrate the efficiency of the present methodological design information using optimization techniques to automatically resolve the building design of proper deposition of web-openings.


Introduction
Currently, a web-opening system which is used in modern architectural buildings and civil structures has been developed in order to obtain functional aspects of structural design such as planned layouts and special efficiency of deposition of facilities.This system can reduce buried spaces which result from sparing expenses of construction per floor in facilities, especially high-rise buildings, by cutting down the amount of material.The structural system of webopenings was fundamentally introduced by Bower [1] in 1968 who investigated ultimate strength of beam structures.Since then, Ward [2] has introduced structural behavioral comparisons between composite and noncomposite beams with web-openings.Improvement problems of bending and shear stiffness of composite beams with web-openings have been treated by Redwood and Demirdjian [3].In Korea, the determination of proper depositions and numbers and shapes of web-openings were investigated, and immediately bending and shear response analyses were carried out.Then their experiments have been carried out in steel and concrete structures by Eon et al. [4] in 1985, Koo [5] in 1998, and Lee et al. [6] in 2003.
According to abovementioned researches, it was verified that web-openings in structures can become substantially problematic in structural safety such as a decline of member stiffness.Since the openings are mainly located in web parts, the decline of shear stiffness especially becomes more obvious than that of bending stiffness.For example, a Vierendeel effect is that an occurrence of a bending moment around webopenings results in a reduction of shear resistance of crosssections; furthermore, a local torsion buckling may occur in web parts of members due to neighboring web-openings.
In general, when web-openings with varied shapes and sizes are placed in members of a given structure, structural stability of the members can be investigated through strain energy, that is, stiffness of structure.Since the web-openings which have an influence on stiffness of structure have varied shapes, numbers, sizes, and so on, optimization techniques that yield maximal stiffness of structure under defined design conditions can be employed for appropriate web-opening's deposition in structure.Although shapes, numbers, and sizes 2 Advances in Materials Science and Engineering of web-openings have to be also considered in the structural design, optimal deposition of web-openings into members is a main interest in this study because shapes, numbers, and sizes must be treated secondly after the determination of depositions of web-openings.
In this study, a topology optimization technique introduced by Bendsøe and Kikuchi [7] in 1988 is utilized in order to decide optimal depositions of web-openings.Topology optimization of structures yields optimal topology as well as optimal shape for global structural systems.In discretization of continuous design domain, a density is defined as a material property of each element, that is, an optimization design parameter.Therefore optimal shape and topology are represented by optimal density distribution contours which have maximal stiffness of structure.
Until now, there have been no determinant and official criteria for structural design with web-openings domestically, and the deposition of web-openings has been designed according to engineers' experiences and usual practices.Thus the objective of the present work is to verify whether applications of the topology optimization approach provide a new structural design method for proper web-opening's depositions.In this method, minimal strain energy or maximal stiffness is defined as objective function, which is satisfied with volume or mass constraints.Under the optimization conditions, an appropriate case with the greatest stiffness of feasible cases of deposition of web-openings may be designed according to optimal results.In advance, feasible parts of web-opening's depositions into void phases have to be sought through results of optimal density distributions of structure without web-openings.

Strain Energy Minimization Based Topology Optimization Problem
2.1.Optimization Problem.In design domain Ω  ⊆ R  ( = 2) which dominates linear elastostatic structures, topology optimization problems are defined as follows: Subject to: ∫ where (1) denotes an objective function , that is, minimal strain energy or maximal stiffness., C, and u are, respectively, strains, material tensors, and displacements.Equations (2) and (3) are optimization constraints.Equation ( 2) is equilibrium.u, b, and t are virtual displacements, body forces, and traction forces, respectively.Equation ( 3) is a volume constraint, and  0 is the limit of feasible volumes in the design domain.

Density Distribution Method.
In topology optimization, material characteristics of each element which are employed through discretization of continuous design domain are defined as element densities.The densities are utilized as design parameters of topology optimization.It is represented as a simple penalty form related to Young's modulus.This is regarded as a density distribution method or SIMP using design domain concepts.Since optimal solutions of SIMP obtain superiority in terms of engineering's aspect and manufacture's ability, it has been practically used for topological optimal design.The penalty relation between Young's modulus and density is written as where E  and Φ  denote Young's modulus and density of element , respectively.The penalty parameter  is used for SIMP. 0 and Φ 0 are, respectively, nominal values of Young's modulus and density.Suppose that the defined structural system follows plane stress state with isotropic materials; material tensor C  of element  can be written including (4) as where ] 0 is a nominal Poisson's ratio.According to topology optimization problems defined as (1) and ( 2), optimal solutions are material density distribution contours with maximal stiffness.The material density values consist of void phases (0, white), solid phases (1, black), and intermediate phases (0 < value < 1, gray) in domain Ω as shown in Figure 1.

Sensitivity Analysis.
Since displacement fields depend on optimization design parameter s, a total sensitivity of objective function fin terms of s is written as a partial derivative introduced by Haug et al. [17] in 1986 as follows: where ∇ ex  f and ∇  f  ∇  u denote an explicit and implicit partial derivative terms, respectively.
Suppose that, in discrete processes, body forces b, traction forces t, differential tensor L, and Jacobi matrix are independent of design parameter s, and finally the total sensitivity formulation of objective function f is simply rewritten as follows: where u  , B  , and C  are nodal displacement vector, operator matrix, and material tensor of element .

Numerical Algorithm for Determination of Optimal Depositions of Web-Openings
The numerical algorithm for optimal depositions of webopenings using topology optimization results is shown in Figure 2, which is extended in [18] and the procedures of deposition design of web-openings are as follows: A D The optimal deposition result of web-openings is used as profitable information for web-opening's design.In continuous design domain (120 cm × 20 cm), square finite elements of 120 × 20 are discretized.As topology optimization problems, an objective function is minimal strain energy (N⋅m) and a volume constraint is restricted by 30% of total volumes.Nominal values of Young's modulus of steel, Poisson's ratio of steel, and an external force are  0 = 2.1 × 10 6 kg/cm 2 , ] 0 = 0.3, and  0 = 360 N, respectively.The penalty parameter [20] for SIMP is  = 5.0.In MBBbeam, one square of 8 cm × 8 cm (finite elements of 8 × 8) is considered as conditions of web-opening.It can be found from Figure 4 that feasible regions of webopenings are presented by models of 4 kinds due to symmetry of MBB-beam.The models of web-openings are illustrated in Figure 5 and Table 1.Here, el and el denote  and  coordinates of elements, respectively.

Stage 3: Execution of Topology Optimization of Each Model and Determination of a Web-Opening's Model with
Maximal Stiffness.Figure 6 shows convergence histories of objective function of strain energy at 50 iterations in topology optimization of MBB-beam with models of web-opening.
It can be seen from Figure 6 that the convergence history of SIMP without web-opening is similar to that of SIMP with web-openings of models b and c.The fact emphasizes the necessity of web-openings in structure with respect to the economical use of materials under structural safety.
Figures 7(a The results of the initial stage are not reliable since the optimization cannot be achieved from the beginning and requires proper iterations or times in order to be satisfied with optimization conditions and obtain optimal solutions.
Note that engineers must be concerned with the final results of optimization.Therefore it can be found that model b or c is the best choice but model a is the worst one as final measurement of each model.
Figure 8 illustrates continuous intermediate density distributions by models of web-opening's depositions.It can be found that the increase of strain energy at the edges of webopening is different according to each model.

Determination of Optimal Depositions of Web-Openings in
Beam-to-Column Connection

Stage 1: Initialization of Design Conditions and Topology
Optimization Problems.As the second test, a linear elastostatic structure with beam-to-column connection [21] is considered.The geometry, boundary, and loading conditions of the analytical model are shown in Figure 9.
The design domain (160 cm × 80 cm) of dash lines is discretized as finite elements of 80 × 40.As topology optimization problems, an objective function is minimal strain energy (N⋅m) and a volume constraint is restricted by 27% of total volumes.Nominal values of Young's modulus of steel, Poisson's ratio of steel, and an external force are  0 = 2.1 × 10 6 kg/cm 2 , ] 0 = 0.3, and  0 = 360 N, respectively.The penalty parameter for SIMP is  = 5.0.In the structure with beam-to-column connection, one square of 16 cm × 16 cm (finite elements of 8 × 8) is considered as conditions of web-opening.

Stage 2: Execution of Topology Optimization of Beamto-Column without Web-Openings and Selection of Feasible
Models of Web-Opening's Depositions.Figure 10 shows results of topology optimization without considering web-opening such as truss shapes [22].Feasible models of web-opening's deposition can be investigated and they are models a, b, c, d, e, and g as shown in Figure 11.Model f is infeasible but here it is considered for comparison with feasible models.Table 2 shows deposition coordinates of models of webopening in design domain.Here, el and el denote  and  coordinates of elements.

Stage 3: Execution of Topology Optimization of Each Model and Determination of a Web-Opening's Model with
Maximal Stiffness.The convergence histories of objective function in topology optimization of the structure with beam-to-column connection with each model of Figure 11 and Table 2 are shown in Figure 12.Figures 12(a Figure 14 shows continuous intermediate density distributions by models a∼e and g.In the structure with beam-tocolumn connection, it can be found that larger strain energy occurs in edges of beam-column connection than in parts of web-opening.
Figure 15 illustrates comparisons of intermediate density distribution between a model f and a case without webopening.It is seen that the existence of web-opening results in the increase of strain energy and the more material density in that part is required.

Conclusions
Structural optimization is a sequential and mathematical technique to achieve the maximum or minimum of objective function that is satisfied with defined constraints in structural problems.This method can optimize the design variables such as topologies, shapes, and sizes for optimal solutions.In particular topology optimization method that yields optimal topologies for the solution utilizes design variables of constant densities into finite elements and its solution is represented as optimal distributions of material densities between 0 and 1.It means that the positions with densities in design domain have to be occupied by materials for structural stiffness and there is no requirement of materials in regions where no densities exist.Therefore the void regions of topology optimization results can become design information for appropriate deposition of web-opening into which it has no material.
Until now the topology optimization technique has been used for optimal design of structures; however in this study      it is proposed that the present method is an engineering practice and idea to be utilized for decision of proper webopening's position.Numerical examples of beam-to-column and simple beam with linear elastostatic problem demonstrate efficiency of the technique for proper web-opening's deposition.
In the web-opening of MBB-beam, the best structure with maximum stiffness is type b according to converged objective function values.In the web-opening of beamto-column connection, the best structure with maximum stiffness is type e according to converged objective function values.As can be seen in the two results, structures with web-openings avoiding regions of many optimal material density distributions have a good performance with respect to structural stiffness.
In the future, the determination of varied variables such as shape, number, and size of web-opening as well as decision of appropriate disposition of web-opening would be also investigated in order to establish official criteria of structural design of web-opening.

4. 1 .
Determination of Optimal Depositions of Web-Openings in MBB-Beam 4.1.1.Stage 1: Initialization of Design Conditions and Topology Optimization Problems.As a 2-dimensional linear elastostatic problem, topology optimization of MBB-(Messerschmitt-Bölkow-Blohm-) beam[19] is carried out.In MBB-beam, the size ratio of  (length) :  (height) is 6 : 1 and a concentrated load is applied at the center side of the upper part.MBB-beam is simply supported by a roller and hinge at left and right sides, respectively.It is shown in Figure3.
Models of Web-Opening's Depositions.Topology optimization results of MBB-beam without web-openings are represented by material density distribution contours and shown in Figure 4.
) and 7(b) illustrate strain energy values of structure in an initial stage and a final stage of topology optimization, respectively.Note that results of the initial stage of the optimization, that is, model a: very good, models b and d: not good, and model c: bad, are different from those of the final stage, that is, models b and c: very good, model d: not good, and model a: bad.

Figure 7 :Figure 8 :Figure 9 :
Figure 7: Strain energy values by models of a∼d of depositions of web-opening.

Figure 10 :Figure 11 :
Figure 10: Optimal solutions of SIMP of a structure with beam-to-column connection without web-opening.

Figure 12 :
Figure 12: Convergence histories of objective function in case of models a∼g of web-opening.
b Type c Type d Type e Type f Type g Objective function (b) Final result

Figure 13 :
Figure 13: Strain energy values by models of a∼g of depositions of web-opening.

Figure 14 :
Figure 14: Continuous intermediate density distributions by models a∼e and g of web-opening's depositions.

Figure 15 :
Figure 15: Continuous intermediate density distributions between original SIMP without web-opening and a model f.

Table 1 :
Models (types) of depositions of web-opening.

Table 2 :
Models (types) of depositions of web-opening.