Development of hybrid performance-based optimization algorithm for structures equipped with vibration damper devices

Nowadays, various types of vibration damping systems are being implemented in different buildings to diminish seismic effects on structures. However, engineers are faced with the challenging task of developing an optimum design for structures utilizing a proper type of damping device based on new techniques such as the performance-based design method. Therefore, this research was aimed at developing a multi-objective optimization algorithm by hybridizing the particle swarm optimization (PSO) and gravitational search algorithm (GSA) to obtain an optimum design for structures equipped with vibration damper devices based on the performance-based design method. Then, the developed hybrid algorithm (PSOGSA) would be capable of optimizing the damping system simultaneously with the optimized details of the structural sections, including the steel rebars, by satisfying all the design criteria. For this purpose, a special process for the design of structures equipped with vibration damper devices according to the performance-based design method was developed by considering of a wide range of vibration damping systems. The proposed PSOGSA optimization framework was then implemented to design a 12-storey reinforced concrete structure equipped with different types of dampers to minimize the structural weight while satisfying all the prescribed performance-based design acceptance criteria. The results indicated that the proposed optimization method was able to successfully optimize the details of the structural members as well as the type and properties of the damper, which significantly improved the structural response in terms of the formation of plastic hinges and the structural movements.


Introduction
The rapid growth of urbanization has resulted in a historic increase in the number of structures and high-rise buildings. Excessive vibrations in buildings are expected to be induced by external forces, such as earthquakes and winds. The induced vibrations may lead to structural damage and the unsatisfactory performance of structures. Induced vibrations can cause significant inconveniences, if not, casualties. Therefore, to enhance the safety and functionality of structures, control of vibrations in structural systems has long been considered by structural designers, and serves as a major source of technological competitiveness [1].
Many vibration control technologies, including damping, vibration isolation, control of excitation forces, and vibration absorbers, have been adopted to decrease the damage and amend the performance of structures. Each system has its own limitations and advantages, and the choice of a particular control system is usually decided by taking into consideration several factors, such as effectiveness, convenience, and life cycle cost [2].
The use of seismic control systems has increased, but choosing the best damper and installing it into a building are very important for reducing vibrations in structures when subjected to seismic loading. The controlling devices reduce damage significantly by increasing the safety and serviceability of the structure, and preventing the building from collapsing during an earthquake. Therefore, much research is being carried out to find the best solution [3]. Some technologies have been developed and adopted to control excessive vibrations, to mitigate their impact on the structural response, and to keep them within sustainable limits during unpredictable events, such as earthquakes [4]. To enhance their performance and ensure their functionality, civil infrastructures require reliable and efficient mechanisms to provide structures with the adequate capacity to maintain their structural integrity under such loadings. Therefore, seismic control systems are needed to enhance the structural behaviour, such as by increasing the energydissipation capacity of civil infrastructures [5].
Cha et al. [6] investigated multiple optimization methods to achieve several targets in which expense and reliability were concurrently optimized to obtain multiple hazardstricken output levels. The direct performance-based design method developed by Cha et al. [7] provides multiple control design layouts for various efficiency levels with different risk levels, using multi-target optimization approaches. Based on the findings, the system stability method determines the storey with the highest total value of near drifts, enabling the seismic fragility of multi-storey buildings to be more precisely measured.
The Liapunov equation was used by García [8] to calculate the mean-square response of the random process and to improve the parameters of a pendulum tuned mass damper (PTMD). The results revealed that the performance of the optimized PTMD is improved when the structure has a linear hysteresis loop. The enhanced PTMD minimizes the mean value of the standard deviation of the principal structural displacement. Furthermore, for a variety of dynamic qualities, the optimized device using the suggested approach decreases vibrations, limits the development of plasticity, and protects the safety of the basic structure on soil. Xian et al. [9] provided an efficient optimization approach for large-scale non-linear viscous dampers that considered both the random seismic excitations and unpredictable damper parameters. To optimize the structure, the method of moving asymptotes (MMA) was applied. Under random seismic excitations, the suggested optimization framework was effectively applied to the stochastic optimum design of uncertain viscous dampers for a long-span suspension bridge. The results implied that, in addition to random seismic excitations, the uncertainty of damper settings must be addressed.
Moghaddam et al. [10,11] proposed an optimum strength distribution method for the seismic design of tall buildings, and thereafter, Hajirasouliha et al. [12] implemented this method to develop a new performance-based optimization technique for RC frames with friction wall dampers [13] to reduce the computation cost [14].
The performance-based earthquake engineering (PBEE) philosophy is aimed at improving seismic risk decisionmaking using scientifically sound assessment and design procedures, which provide stakeholders with informed decision-making options. The first wave of PBEE procedures, such as FEMA-273 1997 [15], represented a significant advance towards standardized decisions. Relations would be formed between structural reaction indices (interstorey drifts, inelastic deformations, and member forces) and performance-orientated definitions, such as immediate occupancy, life safety, and collapse prevention. The performance-based plastic design (PBPD) technique for the design of humpbacked braced frames (BRBFs) was proposed by Sahoo [16]. All the BRBFs tested in this analysis achieved the output goals for yield mechanisms and target drift speeds. As the PBPD is a straightforward design process, no iterations have been done to achieve the performance goals of the BRBF. Gaxiola-Camacho et al. [17] suggested a new reliability appraisal method for filling this expertise gap by providing an alternative to the existing design requirements for the protection of life. Giannakouras and Zeris [18] evaluated the direct displacement-based seismic design (DDBD) provisions for modelling such irregular RC structures and their seismic performance on non-linear static and dynamic analyses. The approach was viewed favourably by the engineering community through major changes that were suggested simultaneously. While the scope of its implementation is constantly expanding, for situations such as frames with a retrograde irregularity no specific guidelines are usually given, arguing that no fundamental changes to the method are expected in such situations. Guidelines are given to ensure that the DDBD approach is used to produce objective performance. Zhang and Tian [19] proposed an optimal seismic design framework for multi-storey reinforced concrete moment frames based on simplified results. The suggested solution minimizes building costs and integrates plastic rotation, and optionally, inter-storey drift as optimization constraints. The developed optimization method consists of two steps, deciding the feasible border region in the field of strength and rigidity, and maximizing the area of the material used.
Based on an extensive review of the literature, as discussed above, the main challenge was identified as the optimization of a structural design that is equipped with vibration damper devices according to the performance-based design method. Since the design of structures based on this method must satisfy a few acceptance criteria for the seismic response of structures, the approach employed is to select the best vibration dissipation type and properties simultaneously, while optimizing the structural members (beams and columns), which is a very complicated task using the conventional optimization techniques.
Therefore, this study attempted to mitigate the extreme effects of earthquake excitations on structures by improving the performance of supplemental vibration dissipation systems, and by optimizing the structural design. Since the optimization process, according to the performance-based design method, must consider several effective parameters regarding the seismic response of structures in terms of structural nodal displacements, structural member forces, plastic hinges, and structural construction costs (building weights), a multi-objective optimization algorithm was developed to optimize all the related parameters at the same time. Since, in the literature the particle swarm optimization approach was unable to search whole parts of the problem (random selection) and resulted in a reduced optimization phase, therefore, in this study the particle swarm optimization (PSO) method was hybridized with the gravitational search algorithm (GSA) to optimize the performance of earthquake energy-dissipation systems simultaneously while optimizing the structural characteristics, including those of the beam and column sections and also the steel reinforcement details.

Development of special procedure for the implementation of the performance-based design method for structures equipped with damper devices
The performance-based design (PBD) method is being increasingly implemented in the design of various structures under seismic excitations. According to this method, the structural members are designed to satisfy predefined criteria for the seismic response of a structure corresponding to the performance of the structure under an expected hazard level. The general concepts of the performancebased design method and its procedures, such as the ASCE 7-10 [20] and ASCE 41-13 [21], can be implemented as seismic provisions in buildings. The performance-based design method is focused on fulfilling the performance objective of a specific type of structure to suit seismic design code requirements through the sizing and detailing of structural members. However, over the last few years, numerous supplementary structural vibration control devices have been developed to reduce seismic damage in structures. These devices can be integrated into structures to improve their structural performance without modifying the size and properties of the structural members.
Therefore, this research attempted to implement the performance-based design (PBD) method on the design of a structure equipped with a vibration damper system. For this purpose, a seismic design process for structures with vibration damper devices according to the performancebased design method was proposed, as shown in Fig. 1, and demonstrated step by step as follows: (1) Define the demand and performance level for the considered structure based on its importance factor and application. (2) Select and scale a set of appropriate seismic acceleration records for both the Design-Basis Earthquake (DBE) and Maximum Considered Earthquake levels according to the considered site condition, as demonstrated in Sects. 11.4.3, 11.4.4 and 11.4.5 of the ASCE 7. (3) Select and define the characteristics corresponding to the considered structural vibration damper devices to be implemented in the structure. (4) Conduct a preliminary design of the structure according to the conventional codes such as the ACI 318-14 under the applied load. (5) Develop a numerical model of the considered structure and the vibration control devices using the finiteelement method to carry out dynamic non-linear timehistory analyses under applied loads. (6) Define and assign the details of the deformation-controlled components according to the performancebased design method to the structural members such as beams and columns. (7) Determine multiple earthquake hazard levels based on spectral response acceleration parameters from seismic hazard maps for multiple performance levels.
The targeted response spectrum is derived from spectral acceleration based on site class effects, and then, check if the building meets the predefined performance levels such as the inter-storey drift, maximum displacement, plastic hinges, and demand-to-capacity ratio of the members. (8) Perform dynamic non-linear time-history analyses for the developed structure equipped with vibration damper devices under considered earthquake records at both the DBE and MCE levels. (9) Review and check the results of the analysis for the response of the considered structure, including the inter-storey drift, plastic hinge formation in the structural members, maximum displacement of the structure, and the demand-to-capacity ratio (member force/ section capacity) for structural members. (10) If the acceptance criteria regarding the performancebased designed method are satisfied for the current structural design and the selected vibration damper devices, then the design process has been completed. (11) If the structure fails to meet the acceptance criteria, then the design process is repeated from the third step by changing the damper properties or by considering another type of damper. (12) If the structure still fails to meet the acceptance criteria after changing the damper properties or using another type of damper, then the design process is repeated from the fourth step by resizing the structural member sections and proceeding to the next step.
More details for some of the above-mentioned steps are provided in the following sections.

Seismic acceleration record
Based on the site location and condition of the structure, a set of earthquake acceleration records were needed to select and scale to fit a continuum of 5% damped acceleration response consistent with Sect. 11.4.5 of the ASCE 7.
Probabilistic analyses of the seismic risk were performed for a degree of function (43-year return period duration as 50% over 30 years) to obtain the Overall Considered Earthquake Shakes, as defined in the ASCE 7 Parts 11.4.3 and 11.4.4, using the existing models to clarify regional seismic origins and earthquake prediction equations. According to Sect. 16.1.3 of the ASCE 7-10, a minimum of three accelerogram sets are required to evaluate the response background. Each chosen accelerogram collection must contain at least two horizontal components, and the vertical component can also be used in rare cases. The analysis needs adequate modes to contain at least 90% of the building weight for each main horizontal action direction.
The complete quadratic combination (CQC) formed blended model responses. The related response parameters, including forces and displacements, were used as the linear response parameters for the practical output calculation.

Preliminary design of the structure
The initial design of the structure in the specific seismic active region was made as the preliminary building design for concrete structures, according to the American codes ACI 318-14 [21]. The material nonlinearities for the structural components were implemented in the finite-element modelling of the considered RC structure, which reflected the inelastic behaviour of the building during strong earthquakes, as per the ASCE 41-13 [21]. The damping matrix was computed using Rayleigh damping, and the natural frequencies were determined to evaluate the time-domain reaction of the non-linear system. The Newmark time-step Fig. 1 Performance-based design process for structures with a vibration damper system integration method was used to evaluate the response of the structure under applied seismic records.

Checking the acceptance criteria
Each structural component was classified as primary or secondary prior to the determination of the conditions for the acceptance of the component. Then, the internal force in the structural elements and the deformation of the structure under an applied force were assessed in compliance with the specifications of the ASCE 41-13, Sect. 7.5.
The acceptance criteria for the concrete beams and columns were defined according to Table 10-7 and 10-8 in the ASCE 41-13. Also, the acceptance criteria for the steel beams and columns are presented in 9-6, ASCE 41-13.
It is commonly appropriate to use the mean demand for the response of a structure associated with the failure mode prediction to use a more restrictive behavioural mode demand, which may have catastrophic implications.
The storey drift (Δ) design obtained from the analyses should not exceed 125% of the drift limit, as stated in the ASCE 7-10, Sect. 12.12.1. (2% of inter-storey drift). In addition, it is advisable to verify the formation of plastic hinges in the structural members in terms of the deformation ratios. Therefore, the degree and distribution of the inelastic specifications of current and new primary elements and components are specified by the demand-to-capacity ratio (DCR) for the considered structural components as: where Q UD is the force due to gravity and earthquake forces measured in line with Sect. 7.5.2 of the ASCE 41-13; Q CE is the predicted component strength or unit strength, as defined in the ASCE 41-13. The DCRs corresponding to the internal forces of the members such as the axial force, shears, and moments, shall be determined for each primary component of a structure. The essential operation of the component is the one with the highest DCR. The largest DCR in a considered storey is referred to as the essential aspect of the DCR in that storey. If there are several components in a specific storey, the element with the largest computed DCR, shall be identified as the essential component for that storey.

Failure to satisfy the acceptance criteria
If the response of the building under applied loads meets the acceptance criteria, then the accepted proposed damper type and characteristics, and details of the structural members (section size, diameter, and number for steel bars) are implemented in the design process; however, if it does not meet the acceptance criteria, then the structure is redesigned and its seismic response is checked again by selecting different characteristics for the damper device or using another type of vibration dissipation system. In this process, changing the properties of the damper is considered first, and the response is checked again against the criteria. However, if it does not pass the acceptance criteria, then consider changing the type of damper. Furthermore, if the design is accepted regarding the damper properties and type, then consider redesigning for the member size to reduce the overall weight of the building, and thus, the cost of the structure. Once the redesigned structure passes the criteria, then the performance-based design process is completed.

Optimization of structure equipped with damper device using PSOGSA
An attempt was made in this research to develop an optimization procedure for the seismic design of a structure equipped with a vibration dissipation system according to the performance-based design method. For this purpose, the multi-objective optimization method was considered to minimize many effective parameters on the response of a structure simultaneously during the optimization process.
Therefore, the PSOGSA computational method for optimizing a structure equipped with an anti-vibration system was developed by combining the particle swarm optimization (PSO) and gravitational search algorithm (GSA) to optimize the properties of the earthquake energy-dissipation system (damper devices) simultaneously, while optimizing the characteristics of the structure, including the size of the structural member sections, and the number and diameter of the steel rebars to satisfy the criteria for the performancebased design method by considering the structural displacement and force capacity of the structural elements such as the beam and columns.

Hybrid PSOGSA optimization method
The PSO method, which applies the concept of social interaction for problem-solving [22], is a robust stochastic optimization technique based on the movement and intelligence of swarms. It uses several agents (particles) that constitute a swarm moving around the search space looking for the best solution [23]. In this technique, each agent is considered as a particle, and in every given iteration, a particle will have a location and a velocity. The best position reached by the particles in an iteration is considered as the local best position, while the best position achieved across all iterations is considered as the global best position. However, the PSO is not able to search in whole domains of the problem (random selection), resulting in an optimization process that is less accurate.
The gravitational search algorithm (GSA), which employs the gravitational law and laws of motion [24], has been applied successfully for solving various non-linear functions. The GSA has a flexible and well-balanced mechanism to enhance exploration and exploitation abilities.
Therefore, based on the advantages of the PSO method and the GSA algorithm, these methods were implemented in this study to develop a hybrid PSOGSA optimization technique to overcome the above-mentioned issues related to the individual PSO method and GSA algorithm.

PSOGSA optimization procedure
As mentioned before in this study, to optimize the type and properties of the structural damper device, and the details of the structural members for minimizing the effects of a seismic load on the structure, the multi-objective PSOGSA optimization method was adopted. Accordingly, the developed computational process for the proposed hybrid particle swarm optimization method and gravitational search algorithm (PSOGSA) was outlined in a flowchart, as illustrated in Fig. 2.
The computational process for the proposed hybrid PSOGSA multi-objective optimization method for the optimization of a structure equipped with damper devices was as follows: 1. The initial population of all the agents was generated as the first step which, in this study, represented the parameters for all the considered damper devices and the characteristics of the structure, including the section size for the beams and columns, and the number and diameter of the steel bars. 2. Each agent was considered as a candidate solution.
Therefore, for each parameter, a random variable was assumed from the range within the upper and lower limits, which had already been set previously. For the beams and columns, predefined 20-section dimensions were considered. However, the diameter and number of steel reinforcements were defined within the specific range. For the characteristics of the damper devices, the upper and lower limits for the effective parameters on the damper performance were defined based on the available range of the considered devices in the market. 3. Then, during the optimization process, these values were selected randomly within the defined upper and lower limits, and based on the objective function in each iteration, the search pattern was changed within the predefined ranges to explore all the possibilities to find the optimum set.
4. Since, both the PSO and GSA methods were integrated as a single system to operate in parallel, although the particle updating for the PSO and GSA was identical, however, the velocity at which the changes were updated was different for both algorithms.
The PSOGSA method substituted the local search of the original PSO with the acceleration of the GSA, which was basically a local best and worst attentive search function. As a result, the suggested algorithm was where Vi(t + 1) is the velocity of agent i at iteration t + 1, w is a weighting function, Vi(t) is the velocity of agent i at iteration t; cj is an acceleration coefficient, rand is a random number between 0 and 1, aci(t) is the acceleration of agent i at iteration t, gbest is the best solution so far, and Xi(t) is the position of agent i at iteration t.
i) Finally, the positions of the agents were updated by ii) The process of updating the velocities and positions would be stopped when an end criterion was met.

Optimization variables
The optimization variables in this study were split into two types in relation to the damper properties and structural characteristics, as demonstrated below: A-Type and characteristics of the structural anti-vibration systems: A (2)  As demonstrated, the first variables included the type of damper device and the characteristics of each structural anti-vibration system (damper devices), where these specific parameters influenced the function of the considered damper device.
Also, the second variables were related to the structural section details, including the section size of beams and columns, and the number and diameter of the steel rebars in each section.
In each iteration, the damper type was changed. It was a randomly selected agent from all the available choices of dampers, and was mixed with the best results (best agents) obtained through the fitness function from previous generations to have the benefit of the best yet results (from previous generations), and to have a chance of searching and assessing better agents from other available options (other damper devices). Although there may be a few redundant optimization generations, it will avoid any risk of missing any possibility of searching for better agents (better results). However, by considering the capability of the PSOGSA for fast convergency, and using advanced computation technologies and facilities, a few redundant generations were not effective in the overall optimization computation time in comparison to other available techniques.

Objective function of optimization
The hybrid PSOGSA optimization method is an advanced technique that focuses on the high ability of the hybrid optimization algorithm to optimize multi-objective problems. Most of the studies available in the literature regarding structural optimization considered either the weight of the structure as representing the considerable cost of construction as an objective function or only the displacement response of the structure as a structural stability optimization target.
However, this research was aimed at minimizing the weight of the structure by reducing the section area of the beams, columns, and steel bars, while considering the acceptance criteria regarding the performance-based design method, which included the structural movement and the structural member forces. Therefore, all these parameters were defined as the objective function for the developed hybrid optimization method.
An inelastic time-history analysis was conducted by implementing at least seven time-history seismic record data according to the ASCE 7-10 to determine the seismic response of the structure.
Then, the peak storey drift (positive and negative for the right and left side movements, respectively) was used to determine the critical state of the structural deformation using the objective function.
Since the structure experienced different storey drifts in various directions, therefore, the maximum movement of the structure in 3 perpendicular directions of X, Y (horizontal movement) and Z (up-down movement) for each earthquake was considered in the optimization process. However, since the vertical seismic response of the building was less than the horizontal movement, thus, the vertical displacement had no significant effect on the optimization process.
The following equations were formulated to determine the effects of displacements in all directions, and the structural weight during the optimization process: where D x , D Y are the sum of the peak inter-storey drifts in the X-and Y-directions, Δ x and Δ y are the relative displacements in storey i in the X and Y-directions, respectively, h i is the height of storey I, and n is number of storeys. Also, W t is the overall structural weight, ΣW C is the total weight of the columns, ΣW B is the total weight of the beams, and ΣW s is the total weight of the steel reinforcement.
Since the displacements in each direction were in different ranges and independently affected the optimization process, therefore, to scale the displacements in the same range, the moderation coefficients, ∅ x and ∅ y , were multiplied with the sum of the inter-storey drifts in the corresponding X and Y-directions. Similarly, a moderation coefficient of ∅ t was multiplied with the sum of the weight of the beams, columns, and steel rebars to scale both the displacements and weights to have an equal impact on the objective function. Thus, the objective function was proposed as follows:

Optimization design constraints
The optimization design constraints adopted in this research included the number of plastic hinges in all the structural members for the loading and unloading events during earthquake excitation. Thus, to consider the optimization design constraints, the penalty function was defined as follows to add it to the objective function: where P is the penalty function, and PHi represents the design constraint, which is the total number of plastic hinges in the structural members. In this equation, CPHi denotes the adjusted coefficient for constraints. A large value was prescribed for CPHi to avoid the occurrence of plastic hinges. Therefore, during the optimization process, the PSOGSA was aimed at circumventing any occurrence of plastic hinges in the structural members (as the optimization constraint) while minimizing structural movements in all three directions in all the storeys (as the objective of the optimization). This was possible through the execution of the penalty function. Then, an auxiliary function was obtained by adding the penalty function to the objective function in the following form: where Ф is the competency function or the equivalent free function.

Considered damper devices
Numerous types of structural vibration damping devices have been developed over the past few years to dissipate vibrations and minimize harmful seismic effects on structures. This equipment is installed in structures as supplementary damping systems to enhance the structural efficiency and stability, without changing the size and properties of the structural members.
In this study, various types of structural damping devices were implemented to enhance the seismic performance of the structure through the dissipation of applied vibrations to the structure. The damper devices that were considered in this study were: oil damper [25], Bingham damper [26], fluid damper, friction damper [27], viscous wall damper [28], viscoelastic (rubber) bracing damper (RBD), rubber wall damper, yielding brace framework (YBF), and volumetric rubber bracing.
The details and optimization parameters for the damper devices that were considered in this study are outlined in the following sections:

Oil damper
An oil damper works by generating a drop in pressure when a viscous fluid passes through a valve or orifice in the piston within the cylinder during movement of the piston, as shown in Fig. 3. The generated damping resistance force in the oil damper is dependent on the speed of the movement and is not greatly affected by the temperature.
The parameters that were considered for a minimum damping force of 250 kN were 110 KN/mm for the internal rigidity (stiffness) and 2000 kN for the maximum damping force with an internal rigidity of 430 KN/mm.
The damping force of the oil damper is a function of the velocity of movement, as shown in the following equation: where C m is the non-linear viscosity coefficient and V d is the speed of movement.
The upper and lower limits that were considered for the effective parameters of the oil damper are listed in Table 1.

Bingham damper
A Bingham damper is made of a silicon-based filling material, steel cylinder, piston, and rod, as shown in Fig. 4. The minimum damping resistance force that was considered was 150 kN, with a non-linear viscosity coefficient of 88.3 KN/(s/mm) 0.1 , while the maximum damping force that was considered was 2000 kN, with a non-linear viscosity coefficient of 1177.4 KN/(s/mm) 0.1 .
The damping resistance force (F d ) is expressed as where, C d is the non-linear viscosity coefficient and μ is the speed.

Friction damper
A friction damper is made up of a series of steel plates that have been specially treated to increase friction. These plates are clipped together and allowed to slip when subjected to a pre-set load, as illustrated in Fig. 5. Friction dampers are designed to prevent slipping under small vibrations such as wind conditions, and they function only when subjected to high vibrations such as a seismic event. The lower limit is generally about 130% and the upper limit is about 75% of the shear force, after which the components of the friction damper will start to fail. In this study, two friction dampers were implemented, and the parameters that were considered for both these friction dampers are listed in Table 1.

Viscous wall damper
A viscous wall damper (VWD) consists of a stiffened steel tank filled with high-viscosity fluid and one or more steel vanes that extend into the tank of viscous fluid. The VWD is installed between two rigid girders within the floors, and the tank is connected to the floor, while the vanes are connected at the above level to dissipate the inter-storey movement, as illustrated in Fig. 6. Table 1 presents all the parameters that were considered for the viscous wall damper in this study.

Rubber bracing damper (RBD)
The developed viscoelastic (rubber) bracing damper (patent No: US20190323183 A1) consists of a steel tube cylinder with a steel inner core moving inside the cylinder. A high damping viscoelastic material, such as high damping rubber, is mounted between the inner core and outer tube, and is bound to both surfaces as energy dissipater members to absorb vibration loads, as shown in Fig. 7.
The resistance force of the device is a function of the effective stiffness and effective damping of the device, where both these parameters are obtained from the hysteresis response of the RBD damper device via an experimental test. The upper and lower ranges that were considered for these parameters for the RBD damper are presented in Table 1.

Rubber wall damper (RWD)
As shown in Fig. 8, a rubber wall damper (Patent No: US10669734B2) consists of a few parallel steel panels that are bound to each other by high damping rubber. It is positioned within floors as wall panels. These dampers are highly capable of dissipating vibration energy and reducing the effects of earthquakes on structures to protect them against severe seismic excitations.
The damping force of a rubber wall damper is calculated using the following equation: where F VWD is the damping force, C VWD is the damping coefficient of the rubber wall damper, V stands for the (12) F VWD = C VWD × V inter-storey velocity, and η is the exponential coefficient for the velocity, which is in the range of 0.5-2.0. In this study, two types of rubber wall dampers were considered, where their corresponding parameters in terms of effective stiffness and effective damping are listed in Table 1.

Volumetric compression restrainer (VCR)
The developed volumetric compression restrainer device (patent No: US10914093B2) consists of a steel tube cylinder with a steel inner core moving inside the cylinder. Several steel divider plates are welded perpendicularly to the inner core and outer cylinder in parallel with each other, and the circular tubes made by the hyperplastic rubber are positioned in the spaces between the divider plates as energy dissipater members, as shown in Fig. 9. During movements of the shaft due to the vibration of a structure, rubber cubes are compressed between the divider plates and act in volumetric compression to generate a noticeable resistant force against the movement and to restrain the vibration. Two types of VCR devices were considered in this study, and the effective stiffness and damping corresponding to each VCR device is presented in Table 1.

Comparison between considered damper devices
A comparison was made of the different damper devices that were considered in terms of their range of performance (damping, stiffness, and stroke), installation condition, and components and details, as listed in Table 2.
As can be seen, among all the damper systems that were considered, the rubber wall damper and volumetric damper devices were the strongest dampers to generate a resistant force against the applied excitation force.

Integrated synchronized computation system
In this study, the ETABS finite-element program was implemented to conduct a structural seismic analysis, and, the PSOGSA optimization algorithm was codified in the MATLAB program platform. Therefore, the MATLAB and the ETABS programs were integrated and synchronized, whereby the data on the structural details were sent to ETABS, and a non-linear analysis of the structure under an applied seismic load was then performed, after which, the retrieved structural seismic response data were sent to MATLAB for the optimization process to be continued. Then, the data were resent again to the ETABS software. This loop was repeated through the whole optimization process, perhaps up to hundreds and thousands of times, until the optimization convergency criteria were satisfied, thereby successfully ending the optimization process.
The two programs were synchronized through the creation of an executable file for the ETABS program, which was implemented using the MATLAB program during the optimization process to receive the structural data (details about the geometry of the structure, structural members, material, boundary conditions, loads, dampers, and isolator devices) from MATLAB, after which the time-history seismic analysis was run. Next, the ETABS executive program sent back the related data regarding the seismic response of the structure in terms of the nodal displacements, structural member force, and the plastic hinges that appeared in the structural members for different applied earthquake loads at different levels of performance for the structure. Then, after sending all the analysis results to MATLAB, the ETABS executive program was closed to be called again using the MATLAB program in the next iteration of optimization.
The executable file for the ETABS program was created using Visual Studio and by choosing Visual Basic as the coding language. Furthermore, the link between the executable file and ETABS was done through the API documentation provided by ETABS (CSI), as shown in Fig. 10.

Application of developed PSOGSA method for optimization of 12-storey RC structure
The developed PSOGSA method was implemented to optimize the seismic design of a 12-storey structure under considered earthquakes according to the performance-based design method. As shown in Fig. 11, the considered structure in this research was a 12-storey, 3-by-5 bay reinforced concrete The columns were positioned 7.32 m apart from each other in both directions. The building was used as a commercial building/office. The seismic importance factor (I) was considered as 1.0. The building was in San Francisco, USA and situated in a high-seismic region. The acceleration spectral response parameters were determined as SS = 1.83, s1 = 0.85, position class D and stiff soil.
The system was categorized into risk Category II, with a seismic category E as per the ASCE 7-10, based on the short reaction acceleration parameters (SDS = 1.22 g at 0.139 s) and long period (SD1 = 0.85 g at 1 s). The columns and beam members were designed and detailed according to the ACI 318-14, using a compressive strength of f′c = 48 MPa for the concrete, and a yield strength of fy = 414 MPa for the steel reinforcement (ASTM Gr. 60).
The details of the beam-column sections are listed in Table 3 and depicted in Fig. 12.
Various cross-sections of the beam and column were defined as listed in Table 4 for use during the optimization process to minimize the weight of the structure and preserve the desired structural response under the applied loads.
For the optimization process, seven different ground motions were selected. The ground movement search criteria are displayed in Table 5, and the selected earthquakes, as shown in Table 6, were applied to the design according to the ASCE 7-10, for different intensities and periods. All the applied time-history data were scaled to match the target response rate.

Results and discussion
This section presents the results of the optimized design for the considered 12-storey structure equipped with damper devices, according to the performance-based design method, using the developed PSOGSA optimization technique.

Fitness function during the optimization process
Variations in the fitness function during the optimization process for the considered 12-storey structure using the PSOGSA method are shown in Fig. 13. The fitness function represented the formation of plastic hinges, the maximum displacements in the two orthogonal directions, and the maximum inter-storey drift for each time-history seismic analysis. The figure shows that the fitness function decreased and became stable within the first 100 iterations of the optimization process, and the convergency of optimization was successfully achieved by minimizing the fitness function in 269 iterations.

Weight of optimized structure
In this section, the details of the beams and columns, and the selected optimum vibration damper type and the specifications of the optimized structure by the developed PSOGSA algorithm are presented in Table 7. As can been seen from the results, the sizes of the sections for both the beams and columns were reduced in the optimized structure compared to the bare frame. Also, among all the considered vibration damping systems, the RBD1 device exhibited the best performance in dissipating the effect of the applied vibrations and satisfying the acceptance criteria for the performancebased design method. Therefore, it was selected as the optimum damper device for all the storeys during the optimization process. The location and configuration of the installed RBD1 damper devices in the optimized 12-storey structure are shown in Fig. 14. Accordingly, the total weight of the concrete and steel materials for all the beam-column sections in the various floors of the considered 12-storey bare frame and the optimized structure are listed in Table 8. As can been from the results, the steel material that was used for the steel reinforcement was reduced by around 13% in storeys 4 to 12 in the optimized structure compared to the bare frame. However, there was a slight increase in the steel bars for the first 3 storeys as more strength was required for the implementation of the vibration damper devices.
Also, the results revealed that the total amount of concrete material for the optimized structure was noticeably reduced in the range of 47% to 34% for storeys 1 to 12, thereby indicating the efficiency of the developed PSOGSA optimization method in minimizing the overall weight of the structure from 4.98 tons to 2.85 tons (42% reduction) and reducing the construction cost of the structure through the implementation of an optimum vibration damping system that satisfied all the required acceptance criteria regarding the seismic response of the structure according to the performance-based design method. Although the optimization process resulted in small sections for the beams and columns, and a reduction of 13% in steel reinforcement for  the 4 th to 12 th storeys, however, relatively large section sizes for the beams and columns for the first-to-third floors were determined. Also, the optimization algorithm calculated a slightly higher amount of reinforcement (steel bars) for the first three storeys in comparison to the initial design. This was because the implementation of damper devices in the structure to control and minimize the movements of the structure through the dissipation of induced structural vibrations caused the generation of noticeable shear forces in the   sections of beams and columns, thereby transferring all the shear forces to the base of the structure, and causing the first few storeys to carry more accumulated forces compared to the higher storeys. Therefore, the optimization output for the structure furnished with damper devices was more obvious for the members in the middle and higher storeys compared to the few first storeys.

Displacement results
The data on seven earthquake records were considered, , the MCER ground motions should consist of not less than seven earthquake acceleration data. Therefore, using seven or more motions, the mean results can be used for the evaluation of the seismic response of a structure. Therefore, in this study, the seven earthquake acceleration records on the Design-Basis Earthquake (DBE) and Maximum Considered Earthquake (MCE) were selected based on the site location of the structure, and the conditions that were scaled to fit a continuum range of 5% damped acceleration response, consistent with Sect. 11.4.5 of the ASCE 7, 2010. Accordingly, the considered structure was subjected to the seven earthquake acceleration records at both the DBE and MCE levels, and a non-linear time-history analysis was conducted for each seismic record. Since the structure was oscillating sideways, therefore, movements to the right and left were considered as positive and negative displacements, respectively. To design the structure for critical conditions, the positive maximum displacement (Dmaxp) and the negative maximum displacement (Dmaxn) of the structural time-history response for all seven earthquakes at the DBE and MCE levels were considered for the bare frame and optimized frame, as shown in Fig. 15. Also, the mean for the positive and negative time-history displacements for each seismic record was calculated for both the bare frame and the optimized frame, as listed in Table 9, and shown in Fig. 16.
Eventually, the mean displacement and variance for the structural response under the applied seven earthquakes at the DBE and MCE levels for the bare frame and optimized frame were calculated, and are presented in Table 10. The results revealed that the implementation of the structural design using the developed optimization method was able to successfully reduce the seismic response of the structure by about 41.3% for the Design-Basis Earthquake, and 41.8% for the Maximum Considered Earthquake in comparison to the bare frame.

Results of storey drift
The overall maximum inter-storey drift along the height of the structure due to the application of the seven types of earthquake records at the DBE and MCE levels for the initial and optimized designs of the structure equipped with damper devices is shown in Fig. 17.
According to the ASCE 7-10, Sect. 12.12.1, the allowable storey drift is 2%. Therefore, as can be seen from the graphs, the inter-storey drift for the MCE level exceeded the allowable movement in both the X and Y-directions in the preliminary designed structure based on the ACI 318-14. However, in the optimized design of the structure, the interstorey drift was within the range for the same MCE level.  Overall, the results indicated that the inter-storey drift in the optimized design structure equipped with a vibration dissipation system was highly reduced in the range of 34% to 47% for both the DBE and MCE records in comparison to the inter-storey drift for the preliminary designed structure, as presented in Table 11. Therefore, it was deduced that the developed PSOGSA method effectively reduced the interstorey drift of the structure using the optimum design of the structure and the type and properties of the optimum structural damping system.

3D interaction yield surface of P-M2-M3 hinges for concrete structure
The moment-rotation curve for a P-M2-M3 hinge is a three-dimensional backbone relationship which is used to describe the post-yield behaviour of a beam-column element subjected to combined axial forces and biaxial bending moments around the local axis in 2-2 and 3-3. Therefore, this interaction curve is implemented to capture the actual non-linear behaviour of RC beams or columns subjected to external forces, and to identify the occurrence of any plastic hinges when the internal force of the members is beyond their capacity (outside of the yield surface). Figure 18 shows the 3D interaction yield surface and internal force position (force generated during application) in the corner column at the first storey of the preliminary designed structure and optimized design structure for all the various combinations of P, M2 and M3 (axial force, moment around 2-2 axes, moment around 3-3 axes) during the optimization process.
As can be seen from the results, the optimization process managed to shift 100% of the points representative of the internal force of the structural members from the outside of the three-dimensional P-M2-M3 yield surface (plastic or yielding zone) to the inside of the curve (elastic zone) by implementing the optimized design and utilizing the optimum damper devices simultaneously with a reduction in the total weight of the structure (Fig. 18).
Similarly, as can be seen in Fig. 18, about 76%, 59%, and 66% of the hinges in the same corner column (on the first floor) were shifted from out of the yield surface for the M2-M3, P-M2 and P-M3, respectively to the elastic zone (inside of the curve) during the optimization of the structure equipped with damper devices.
Therefore, these results confirmed that the developed PSOGSA algorithm was able to effectively optimize the seismic design of the structure by utilizing a vibration damping system to minimize damage in the structural members due to an applied earthquake load by satisfying the acceptance criteria for the performance-based design method.

Results on occurrence of plastic hinges
The overall number of occurrences of plastic hinges in the structural members and also the percentage of reduction in the occurrence of plastic hinges during the application of the seven types of DBE and DBE accelerations that were considered in the various optimization iterations are shown in Fig. 20, while the details are listed in Table 12 for the bare frame and optimized structure.
As can be seen from the results of the bare frame (structure before optimization and without vibration dissipation system) in Fig. 19a, and in Table 12 (for the DBE and MCE), only a few plastic hinges were observed in the columns in some of the applied earthquake records at the DBE level for the bare frame in the performance of the IO-LS and LS-CP. However, the number of plastic hinges that were formed in the columns of the bare frame for the various applied MCE records was in range of 10 to 22 plastic hinges for the performance of the LS-CP and 20 to 78 plastic hinges for the performance of the LS-IO. Therefore, it was revealed that the bare frame was partially damaged at the DBE level, but for the considered MCE, the structure was highly vulnerable in the LS-CP and almost collapsed in the LS-IO.
However, the results presented in Table 12 showed that there were no plastic hinges in the structural members of the optimized structure at both the DBE and MCE levels for the performance of all the various IO-LS and LS-CP of the structure, and as can be seen in Fig. 19b, the overall percentage of reduction in the formation of plastic hinges in the structure was 100% for the seven earthquake records that were applied at the DBE and MCE levels, respectively. Therefore, it was revealed that the developed POSGSA optimization method was able to successfully optimize the seismic design of the 12-storey structure throughout the optimization iterations by utilizing an appropriate vibration damper device to prevent the formation of any plastic hinges in the structural members for both the MCE and DBE seismic excitations.
Additionally, a non-linear time-history analysis was carried out for the MCE records using predicted structural strengths to survive trembling, as described in Section C1.3.1 of the ASCE 7-10, with decreasing likelihood (up

Inter-story Drift
Optimized Design Structure Preliminary Design Structure Fig. 17 Inter-storey drift along height for the preliminary design structure and optimized design structure for DBE and MCE records from 10% of complete and partial collapse). Then, according to the code statement, the optimized structure would be on the safe side as there would be no failure in any structural members. Therefore, it could be concluded that the developed PSOGSA algorithm successfully optimized the design of the structure equipped with damper devices according to the performance-based design method to protect the structure from any seismic damage at different earthquake levels (DBE and MCE) and within the range of the operational level of the structure (LS-IO).

Base reaction for design-basis earthquake and maximum considered earthquake
The effects of the vibration dissipation devices installed in the optimized frame on the base reaction, including the shear forces and bending moments, during a seismic excitation were studied for both a Design-Basis Earthquake and Maximum Considered Earthquake.

Base reaction in the structure for a design-basis earthquake
The analysis results demonstrated that of the overall base reaction due to the applied DBE decreased when the structure was equipped with dampers, as shown in Table 13, for the bare frame and optimized structure, respectively. The results illustrated that the dampers had a noteworthy effect on the base shear forces; however, there was not much effect on the moment forces. This reduction in the maximum base shear in the horizontal X-direction was about 27%, 39%, 34%, 9%, 26%, 30%, and 3% for the DB1 to DB7 earthquakes, respectively. Also, the minimum base shear in the same direction (X), was reduced by around 48%, 23%, 37%, 25%, 32%, 18%, and 33% in the optimized structure equipped with the selected damper device for DBE1 to DBE7, respectively.
Likewise, the maximum and minimum shear forces in the horizontal Y-direction were reduced in the optimized structure in the range of 13 to 34% and 4 to 25%, respectively.
In contrast, there was not much effect on the base moment in the optimized structure equipped with damper devices and the bare frame.
The noticeable reduction in the base shear forces was due to the action of the implemented vibration damping system in the structure, which absorbed a part of the applied vibration energy and dissipated its effect in the structural response.

Base reaction in the structure for maximum considered earthquake
The results of the overall base reaction for the bare frame and optimized structure under the applied MCE are presented in Table 14, respectively. The results illustrated that the maximum base shear in the horizontal X-direction in the optimized structure was reduced by about 27%, 39%, 34%, 9.3%, 26%, 30%, and 3% for the applied MCE1 to MCE7, respectively. Similarly, the reduction for the minimum base shear in the same direction was 49%, 23%, 36%, 25%, 32%, 18%, and 32% for MCE1 to MCE7, respectively. Accordingly, the same reductions were observed for the minimum and maximum base shear in the Y-direction in the range of 13 to 34% and 4 to 25%, respectively for the applied MCE due to the implementation of the optimized design of the structure and the vibration dissipation system. However, again, there was not much effect on the moment forces.

Demand-To-Capacity ratio of the columns
The demand-to-capacity ratio of the columns for the bare frame compared to the optimized structure was evaluated by considering the effects of the vibration dissipation device installed in the frame during the Design-Basis Earthquake and the Maximum Considered Earthquake.

Demand-to-capacity ratio of the columns for design-basis earthquake
The results revealed that the overall demand-to-capacity ratio for the columns decreased in the optimized structure equipped with dampers, as shown in Table 15 and Fig. 20.
The results illustrated that the implementation of damper devices in the structure had a noteworthy effect on the capacity of the structure. This corresponded with a decrease of about 99.8%, 98.7%, 99.5%, 99.4%, 99.5%, 99.99%, and 99.99% in the demand-to-capacity ratio of the columns

Demand-to-capacity ratio of the columns for maximum considered earthquake
The results of the overall demand-to-capacity ratio of the columns in the bare frame and optimized structure are presented in Table 16 and Fig. 21.
As can be seen from the results, the utilization of damper devices in the optimized structure led to a decrease in the demand-to-capacity ratio of the columns of about 98.3%, 98.5%, 99.8%, 97.4%, 99.7%, 99.4%, and 98.3% for the considered MCE1 to MCE 7, respectively for immediate occupancy, life safety and collapse prevention operational levels.
Therefore, the results proved that the use of the optimized design with the damper devices had a noteworthy effect on the capacity of the structure.

Validation of developed PSOGSA algorithm through other optimization methods
To validate the results of the PSOGSA optimization algorithm, it was applied for the optimization of a 12-storey RC structure, and the results in terms of a reduction in the maximum displacement response were investigated. The same process was performed using other optimization techniques such as the genetic algorithm (GA), ant colony optimization (ACO) and particle swarm optimization (PSO), and THE  results were compared to the output of the PSOGSA optimization method. Then, the optimization outputs were evaluated with respect to the fitness function, plastic hinge development in the structural members, and maximum displacements during earthquake excitations at both the DBE and MCE levels, as demonstrated in detail as follows: Figure 22 shows the decrease in the overall fitness function optimization process. Furthermore, the fitness was plotted in respect of the occurrence of plastic hinges, maximum displacement in the two applied orthogonal time-history data, the maximum inter-storey drift, and the weight of the concrete.

Fitness function
Additionally, the PSOGSA found the optimal solution after around 250 iterations in comparison with the other algorithms that found it after 900, 900, and 650 iterations for the GA (genetic algorithm), ACO (ant colony optimization), and PSO (particle swarm optimization), respectively. Therefore, much less computational power and time were required to finish the optimization process.

Plastic hinges
The reduction in the formation of plastic hinges in the structural members under earthquake excitations at the DBE and MCE levels for the optimized structure for all optimization iterations using various optimization methods is depicted in Fig. 23.
As can be seen from the graph in Fig. 23, both the PSOGSA and GA methods resulted in the greatest reduction of 88.2% in the formation of plastic hinges in the beams and columns under the DBE, while the developed PSOGSA algorithm was able to obtain the best design within just 269 iterations in comparison to the GA method with 900 iterations. Two other methods, namely the PSO and AC methods, were able to reduce the occurrence of plastic hinges to 64.7% and 41.1% within 650 and 900 iterations, respectively. These results indicated the efficiency of the developed PSOGSA method in obtaining the maximum reduction of plastic hinges in the structural members via a minimum computation time and effort in comparison to other conventional optimization methods.
The same graph corresponding to the MCE level in Fig. 23 also revealed that all the optimization methods resulted in a reduction of 73% in the occurrence of plastic hinges in the structural members. Hence, the PSOGSA algorithm produced results after only 269 iterations, thereby proving the very fast convergency of this technique in comparison to the other methods, which obtained an output after 650 and 900 iterations. Deduction percantage of total plastic hinges occurance Iteration Fig. 19 Results on the formation of plastic hinges during the optimization process

Maximum displacements
The percentage reductions in the Maximum Positive Displacement (movement to the right side) and Maximum Negative Displacement (movement to the left side) in the X and Y-directions in the various iteration steps of the different optimization techniques for the structure under the Design-Basis Earthquake (DBE) and Maximum Considered Earthquake (MCE) are listed in Table 17 and presented in Figs. 24 and 25.
As can be seen from the results, the PSOGSA and GA methods produced a maximum reduction of maximum displacements at both the DBE and MCE levels in the range of 41% ~ 43% in the X-direction and 23% to 27% in the Y-direction (perpendicular to the X-direction), while the PSO and ACO methods obtained maximum displacement reductions in the range of 38% to 43% in the X-direction and 24% to 27% in the Y-direction.
However, the results indicated that the PSOGSA method obtained the optimum point at an earlier stage after 269 iterations, while the PSO method reached the optimum point after 650 iterations, and both the GA and AC techniques had to attempt 900 iterations to complete the process.
Therefore, this proved the remarkable capability and advantages of the developed PSOGSA method to obtain an optimization convergency in a very short time compared to other similar methods, thereby leading to significant savings in computation time and cost. This advantage would be even more significant for the optimization of mega structures, high-rise buildings, and towers, which require strong and expensive computational facilities to process the mathematics and calculations for the optimization of structures furnished with damper devices.

Conclusions
The current research study developed a hybrid PSOGSA technique for the multi-objective optimization of a structure subjected to seismic excitations through the utilization of vibrational damper devices according to the performance-based design method. In this method, the details of the structural design, including the member sections and steel reinforcement, as well as the type and characteristics of the structural vibration damping system were optimized simultaneously to minimize the cost of construction and to satisfy the acceptance criteria for the performance-based design method, including structural movements and the occurrence of plastic hinges in the beams and columns. For this purpose, seven types of structural damping systems were considered, and seven earthquake records  at the DBE and MCE levels were adopted for application to the structure.
The results proved that the developed PSOGSA optimization method was able to effectively minimize the effect of earthquakes on the structure through the utilization of an optimum design and vibration damping system. This technique also greatly improved the integrity and protection of the building against the applied DBE and MCE, as demonstrated by some of the highlighted results as follows: • The outcome of this study proved that the reduction in the occurrence of plastic hinges in the optimized structure equipped with optimum damping devices was 100% for both the MCE and DBE levels. Therefore, it could be  concluded that the developed hybrid PSOGSA method effectively resulted in the optimum design of a structure equipped with the best vibration damping system that prevented any damage to the structural members, and thus, it was a more resilient structure compared to the bare frame. • The results also indicated that the weight of the optimized structure was reduced by about 42% by minimizing the sectional area of the beams and columns, and the number and diameter of the steel bars, while also by considering the minimization of the displacement of the structure under applied seismic loads in the objective function. • The displacement response of the optimized structure was reduced by about 41% and 23% in the X-and Y-directions under the applied MCE and DBE in comparison with the bare frame, which indicated the effec- tiveness of implementing the optimum damper devices to minimize movement of the structure. • The inter-storey drift in the optimized design structure with the vibration dissipation system was reduced in the range of 34% to 47% for both the DBE and MCE records in comparison to the drifts in the bare frame. • There was no formation of plastic hinges in the beam and column sections of the optimized structure, and it was revealed that the implementation of the optimum damper devices in the structure successfully diminished the effect of the DBE and even the MCE on the structure. • The application of the optimum damping devices in the optimized structure resulted in a reduction in the base shear of the building in the range of 4% to 48% for the various applied DBE and MCE levels. • A comparison of the PSOGSA technique with other optimization methods such as the genetic algorithm (GA), ant colony optimization (ACO), and particle swarm optimization (PSO) proved the high efficiency of this method in obtaining the optimum results with minimum computation time and effort. Therefore, the PSOGSA method exhibited a fast convergency rate in the range of 61% ~ 70% in comparison to the other optimization techniques that were studied.
Data availability All related data for this research are presented and available in the manuscript.

Declarations
Conflicts of interest It is declared that there is no conflict of interest in this research work.
Ethical standards There are no ethical concerns in this research work.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.