Optimization design method for vibration control structures based on multi-level step-yielding dampers

Abstract

The present invention discloses an optimization design method for vibration control structure based on multi-level graduated yield damper. The method includes the steps of inputting an uncontrolled vibration control structure model, determining the installation position of the damper and the structure optimization prediction, determining and installing the preset parameters of the algorithm program, and performing the analysis calculation of the algorithm program independently for the load modes of small earthquakes, medium earthquakes, and large earthquakes, obtaining the damper yield load capacity, yield displacement, and hysteresis curve corresponding to the maximum fitness value in the corresponding mode, constructing the skeleton curve of the multi-level stage yield damper according to the hysteresis curves in the three modes, and then back-calculating the hysteresis curve of the multi-level stage yield damper according to the constructed skeleton curve, and then designing the subsequent multi-level stage yield damper.Compared with the original structure, the optimized structure obtained by the optimization design method of the present invention can not only maximize the performance of the damper, but also meet the control requirements of the structure inter-story displacement angle.

Description

[Technical field]

[0001] The present invention relates to the technical field of damper vibration control structure, and particularly to an optimization design method for a vibration control structure based on multi-level step-yield dampers. [Background technology]

[0002] Energy dissipation vibration control technology is a seismic technology that is currently widely used in countries around the world, and is commonly applied both at home and abroad. Currently, various types of dampers have been developed. Current conventional technology mainly focuses on the placement position of the damper, and the type of damper is mainly focused on viscous dampers, but there is currently no research on the parameters of step-yielding dampers. Therefore, in terms of applying step-yielding dampers to vibration control, maximizing the performance of step-yielding dampers and improving the vibration control effect of the damper while increasing economic efficiency will greatly contribute to vibration control of structures in actual construction. Summary of the Invention [Problem to be solved by the invention]

[0003] SUMMARY OF THE PRESENT EMBODIMENT The objective of the present invention is to provide an optimization design method for a vibration control structure based on a multi-level graduated yield damper, in order to solve the problems presented in the background art above. [Means for solving the problem]

[0004] In order to achieve the above objectives, the present invention provides an optimization design method for a vibration control structure based on a multi-level step-yield damper, S1, inputting an uncontrolled vibration control structure model into SAP2000 finite element software; S2, determining the installation position of the damper and the structural optimization prediction; S3, determining and setting preset parameters of the algorithm program; S4, setting a small earthquake load mode; S5, enter the algorithm program, perform analytical calculation, and output the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value in the small earthquake load mode; S6, establishing a moderate earthquake load mode; S7, enter the algorithm program, perform analysis calculation, and output the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value in the moderate earthquake load mode; S8, setting a large earthquake load mode; S9, enter the algorithm program, perform analytical calculations, and output the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value in the large earthquake load mode; S10, constructing a skeleton curve of a multi-level step-yield damper according to the output hysteresis curves of the three modes; S11, back-calculating the hysteresis curve of the multi-level step-yield damper by the constructed skeleton curve; S12, designing and installing a multi-level graduated yield damper.

[0005] In a preferred embodiment, in step S3, the preset parameters of the algorithm program include the initial population number, population random range interval, binary coding length, seismic waves and seismic modes of the vibration control structure model, and structural optimization forecast.

[0006] In a preferred embodiment, the step of entering an algorithmic program and performing analytical calculations comprises: S101, setting initial data of an algorithm program; S102, randomly selecting a damper yield load capacity and stiffness value for a fixed number of pairs within a preset range; S103, inputting an initial data file; S104, using python (registered trademark) software to read damper yield load capacity and stiffness value data; S105, changing the parameters of the damper installed by SAP2000 software; S106, performing analysis using software SAP2000, outputting structural analysis results, and obtaining result data files; S107, reading the result data using matlab software, and substituting the result data into a fitness function to calculate the fitness; S108, analyzing the vibration control result, outputting the damper yield load capacity and yield displacement corresponding to the maximum fitness value in the current population, and judging whether the number of iterations has been reached; S109, if the number of iterations has not been reached, a new initial population is obtained through selection, crossover, and mutation, and the process is repeated again; S110, if the number of iterations is reached, terminating the calculation and outputting the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value calculated by the algorithm program.

[0007] In a preferred embodiment, in step S101, the initial data of the algorithm program includes the number of algorithm iterations, the stiffness and yield load capacity range of the damper, the genetic selection rate, the genetic mutation rate, and the genetic crossover rate.

[0008] In a preferred embodiment, in step S107, the ratio of the absolute value of the difference between the maximum inter-story shear force and the top story acceleration after the vibration control structure optimization and the corresponding values ​​of the original structure to the corresponding values ​​of the original structure is calculated as a percentage, and then the square root of the sum of squares is taken as a fitness calculation formula,

[0009]

number

[0010] Here, V max is the maximum inter-story shear force of the optimized structure (kN), V is the maximum inter-story shear force of the original structure (kN), and V i is the inter-story shear force (kN) of each story of the structure, and v2 is the maximum top story acceleration of the optimized structure (mm / s 2), and v1 is the acceleration of the top floor of the original structure (mm / s 2 ) and fitness is the fitness function.

[0011] In a preferred embodiment, in step S107, the ratio of the absolute value of the difference between the maximum inter-story shear force and the top story acceleration after the vibration control structure optimization and the corresponding values ​​of the original structure to the corresponding values ​​of the original structure is calculated as a percentage, and the ratio is used as a fitness calculation formula;

[0012]

number

[0013] Here, V max is the maximum inter-story shear force of the optimized structure (kN), V is the maximum inter-story shear force of the original structure (kN), and V i is the inter-story shear force (kN) of each story of the structure, and v2 is the maximum top story acceleration of the optimized structure (mm / s 2 ), and v1 is the acceleration of the top floor of the original structure (mm / s 2 ), where fitness is a fitness function and α, β are weighting coefficients.

[0014] In a preferred embodiment, in step S107, a ratio of an absolute value of a difference between the maximum inter-story displacement angle after the vibration-damping structure optimization and a corresponding value of the original structure to the corresponding value of the original structure is calculated as a percentage to obtain a fitness calculation formula;

[0015]

number

[0016] Here, θ max is the maximum inter-story displacement angle of the optimized structure, θ is the maximum inter-story displacement angle of the original structure, and fitness is the fitness function.

[0017] In a preferred embodiment, in step S107, the ratio of the absolute value of the difference between the maximum inter-story shear force and the top story acceleration after the vibration control structure optimization and the corresponding values ​​of the original structure to the corresponding values ​​of the original structure is calculated as a percentage, and the ratio is used as a fitness calculation formula;

[0018]

number

[0019] Here, θ max is the maximum story drift angle of the optimized structure, θ is the maximum story drift angle of the original structure, and a max is the maximum acceleration of the optimized structure, a is the maximum inter-story acceleration of the original structure, and u max is the maximum interlayer displacement of the optimized structure, u is the maximum interlayer displacement of the original structure, α, β, and χ are weighting coefficients, and fitness is the fitness function.

[0020] In a preferred embodiment, the structural analysis results output in step S106 include the inter-story shear force of each structural layer, the top story acceleration, the structural inter-story displacement angle, the inter-story acceleration, the maximum inter-story displacement, the damper yield load capacity, the damper displacement, and the hysteresis curve.

[0021] In a preferred embodiment, the population random range interval includes upper and lower limits of the damper preset stiffness from the stiffness of the damper that meets the structural maximum story drift angle required by the code to 25% of the uncontrolled structural maximum story stiffness, and upper and lower limits of the damper preset load capacity are set from 0 to 25% of the structural maximum story shear force.

[0022] Compared with the prior art, the beneficial effects of the present invention are as follows: the optimization design method of the present invention independently performs calculation analysis for small earthquake load mode, medium earthquake load mode, and large earthquake load mode, respectively, obtains the damper yield load capacity, yield displacement, and hysteresis curve corresponding to the maximum value of the fitness in the corresponding mode, constructs the skeleton curve of the multi-level stage yield damper according to the hysteresis curve in the three modes, and reverse-calculates the hysteresis curve of the multi-level stage yield damper according to the constructed skeleton curve, and then carries out the subsequent design of the multi-level stage yield damper. The optimized structure obtained by the optimization design method of the present invention can not only maximize the performance of the damper compared to the original structure, but also meet the control requirements of the structure inter-story displacement angle. Under the action of different intensity earthquakes, the vibration control structure optimization design method based on the multi-level stage yield damper proposed in the present invention all have good structural vibration control effect, are suitable for popularization and application to other displacement type stage yield dampers, and have important reference value and engineering application value in the optimization of architectural structural vibration control design. [Brief description of the drawings]

[0023] [Figure 1] 1 is a flowchart of an optimization design method according to a preferred embodiment of the present invention. [Diagram 2] 1 is a flowchart of an analytical calculation by an algorithm program according to a preferred embodiment of the present invention. [Diagram 3] FIG. 4 is a schematic diagram of an installation position of a damper according to an embodiment of the present invention. [Figure 4] FIG. 2 is a schematic diagram of a seismic wave according to an embodiment of the present invention. [Figure 5A] 13 is a corresponding hysteresis curve for a small earthquake loading mode according to an embodiment of the present invention. [Figure 5B] 13 is a corresponding hysteresis curve for a moderate earthquake loading mode according to an embodiment of the present invention. [Figure 5C] 13 is a corresponding hysteresis curve for a large earthquake loading mode according to an embodiment of the present invention. [Figure 6]4 is a schematic diagram of the damper yield load capacity and yield displacement output of a stage yield damper according to one embodiment of the present invention. FIG. [Figure 7] FIG. 13 is a schematic diagram of story displacement angle data before and after optimization design according to an embodiment of the present invention. [Figure 8] FIG. 13 is a schematic diagram of inter-story acceleration data before and after optimization design according to an embodiment of the present invention. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0024] The technical solutions in the embodiments of the present invention are clearly and completely described below, and all other embodiments obtained by those skilled in the art without any creative work are within the scope of protection of the present invention.

[0025] Example 1 As shown in FIGS. 1 and 2, the optimization design method for a vibration control structure based on a multi-level graduated yield damper according to a preferred embodiment of the present invention includes the following steps: Step S1, an uncontrolled vibration control structure model is input into the SAP2000 finite element software, and the uncontrolled vibration control structure model is an original structure model without a damper. Step S2: determine the installation position of the damper and the structural optimization prediction, including but not limited to the prediction of meeting the reference inter-story displacement angle limit value and inter-story acceleration optimization rate. Step S3, determine and set preset parameters of algorithm program. Specifically, the preset parameters of algorithm program include initial population number, population random range interval, binary coding length, seismic wave and seismic mode of vibration control structure model, and structural optimization prediction. Step S4: Set the small earthquake load mode. Step S5: Enter into the algorithm program to perform analytical calculations, and output the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value in the small earthquake load mode. Step S6: Set the moderate earthquake load mode. Step S7: Enter into the algorithm program to perform analytical calculations, and output the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value in the moderate earthquake load mode. Step S8: Set the large earthquake load mode. Step S9: Enter into the algorithm program to perform analytical calculations and output the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value in the large earthquake load mode.

[0026] Specifically, as shown in FIG. 2, the steps of entering the algorithm program and performing the analysis calculation in steps S5, S7, and S9 are as follows: step S101 of setting the initial data of the algorithm program; step S102 of randomly selecting the damper yield load capacity and stiffness value of a fixed number of pairs within a preset range; step S103 of inputting the initial data file; step S104 of reading the damper yield load capacity and stiffness value data using the python software; step S105 of changing the parameters of the damper set by the SAP2000 software; step S106 of inputting the parameters of the damper into the software SAP2000 to execute the analysis, outputting the data of the structural analysis result, and storing the result data. step S106 of obtaining a data file; step S107 of reading the result data using MATLAB software and substituting the result data into a fitness function to calculate the fitness; step S108 of analyzing the vibration control result, outputting the damper yield load-bearing capacity and yield displacement corresponding to the maximum fitness value in the current population, and judging whether the number of iterations has been reached; step S109 of obtaining a new initial population through selection, crossover and mutation and repeating again if the number of iterations has not been reached; and step S110 of terminating the calculation and outputting the damper yield load-bearing capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value calculated by the algorithm program if the number of iterations has been reached.

[0027] In step S101, the initial data of the algorithm program includes the number of algorithm iterations, the stiffness and yield load capacity range of the damper, the genetic selection rate, the genetic mutation rate, and the genetic crossover rate.

[0028] In step S107, the fitness is obtained by one of the following four calculation formulas, but is not limited to these, and the fitness function can be adjusted as necessary.

[0029] 1. Calculate the absolute value of the difference between the maximum inter-story shear force and the top story acceleration after vibration control structure optimization and the corresponding values ​​of the original structure, and the ratio to the corresponding values ​​of the original structure in percentage. Then, take the square root of the sum of squares to obtain the fitness calculation formula.

[0030]

number

[0031] Here, V max is the maximum inter-story shear force of the optimized structure (kN), V is the maximum inter-story shear force of the original structure (kN), and V i is the inter-story shear force (kN) of each story of the structure, and v2 is the maximum top story acceleration of the optimized structure (mm / s 2 ), and v1 is the acceleration of the top floor of the original structure (mm / s 2 ) and fitness is the fitness function.

[0032] The fitness calculation using formula (1) is applied when it is necessary to comprehensively consider the optimization rate of the top story acceleration and maximum inter-story shear force of the structure and control the seismic response of the structure.

[0033] 2. The fitness calculation formula is calculated by calculating the absolute value of the difference between the maximum inter-story shear force and the top story acceleration after vibration control structure optimization and the corresponding values ​​of the original structure, and the ratio of these values ​​to the corresponding values ​​of the original structure, expressed as a percentage.

[0034]

number

[0035] Here, V max is the maximum inter-story shear force of the optimized structure (kN), V is the maximum inter-story shear force of the original structure (kN), and V iis the inter-story shear force (kN) of each story of the structure, and v2 is the maximum top story acceleration of the optimized structure (mm / s 2 ), and v1 is the acceleration of the top floor of the original structure (mm / s 2 ), where fitness is a fitness function and α, β are weighting coefficients.

[0036] The fitness calculation using equation (2) is applied when controlling seismic response by optimizing the inter-story shear force and top story acceleration at different ratios.

[0037] 3. The fitness calculation formula is calculated by calculating the ratio of the absolute value of the difference between the maximum story displacement angle after the vibration control structure optimization and the corresponding value of the original structure to the corresponding value of the original structure in percentage;

[0038]

number

[0039] Here, θ max is the maximum inter-story displacement angle of the optimized structure, θ is the maximum inter-story displacement angle of the original structure, and fitness is the fitness function.

[0040] The fitness calculation according to Equation (3) is applied when controlling the inter-story displacement of the structure by considering the displacement angle optimization rate of the structure, so as to meet the requirement of the reference structural displacement angle and control the structural deformation.

[0041] 4. The fitness calculation formula is calculated by calculating the absolute value of the difference between the maximum inter-story shear force and the top story acceleration after vibration control structure optimization and the corresponding values ​​of the original structure, and the ratio of these values ​​to the corresponding values ​​of the original structure, expressed as a percentage.

[0042]

number

[0043] Here, θ max is the maximum story drift angle of the optimized structure, θ is the maximum story drift angle of the original structure, and a maxis the maximum acceleration of the optimized structure, a is the maximum inter-story acceleration of the original structure, and u max is the maximum interlayer displacement of the optimized structure, u is the maximum interlayer displacement of the original structure, α, β, and χ are weighting coefficients, and fitness is the fitness function.

[0044] The fitness calculation according to Equation (4) is applied when the optimization rates of the structure inter-story displacement angle, the layer angular acceleration, and the maximum inter-story displacement at different ratios are comprehensively considered to avoid the unification of the objective function.

[0045] Furthermore, the structural analysis results output in step S106 include the inter-story shear force of each structural story, the top story acceleration, the structural inter-story displacement angle, the inter-story acceleration, the maximum inter-story displacement, the damper yield load capacity, the damper displacement, and the hysteresis curve.

[0046] Furthermore, the population random range section includes upper and lower limits of the damper preset stiffness being set from the stiffness of the damper that satisfies the structural maximum story displacement angle required by the standard to 25% of the uncontrolled structural maximum story stiffness, and upper and lower limits of the damper preset load capacity being set from 0 to 25% of the structural maximum story shear force.

[0047] Step S10: construct a skeleton curve of the multi-level step-yield damper according to the output hysteresis curves of the three modes.

[0048] Step S11: Back-calculate the hysteresis curve of the multi-level step-yield damper according to the constructed skeleton curve.

[0049] Step S12: The optimized design of the vibration control structure is completed by designing and installing a multi-level step-yield damper.

[0050] Specifically, the skeleton curve of the step yield damper is constructed by the hysteresis curves of three modes, so that the skeleton curve has various hysteresis rules, and the hysteresis rules of the step yield damper are defined according to the actual situation, and finally the step yield damper is designed. The design concept and form of the step yield damper include, but are not limited to, sleeve ring and pre-gap type.

[0051] Example 2 In a specific embodiment, a nine-layer bench-mark model is taken as an example to carry out an optimization design method for a vibration control structure based on a multi-level graduated yield damper, which specifically includes the following steps: Step S1, an uncontrolled vibration control structure model is input into the SAP2000 finite element software, and the uncontrolled vibration control structure model is an original structure model without a damper. Step S2: Determine the installation positions of the dampers and the structure optimization prediction. As shown in Fig. 3, four dampers 1 are symmetrically arranged in each of the X and Y directions. Step S3: determine and set the preset parameters of the algorithm program. Specifically, the preset parameters of the algorithm program include the following: the initial population number is set to 20 (usually set to 10-50, the larger the population number, the longer the calculation time); the population random range section is from 0 to 25% of the inter-story shear force and stiffness; the binary coding length is set to 15 (the length is determined by the population random range section); the earthquake mode includes small earthquake acceleration 70cm / s^2, medium earthquake acceleration 200cm / s^2, and large earthquake acceleration 400cm / s^2; the structural optimization forecast is that the inter-story displacement angle optimization rate is 30% and the inter-story acceleration optimization rate is 15% (general provision). Step S4: Set the small earthquake load mode. Step S5: Enter the algorithm program to perform analytical calculations, and output the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value in the small earthquake load mode. Step S6: Set the moderate earthquake load mode. Step S7: Enter into the algorithm program to perform analytical calculations, and output the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value in the moderate earthquake load mode. Step S8: Set the large earthquake load mode. Step S9: Enter into the algorithm program to perform analytical calculations and output the damper yield load capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value in the large earthquake load mode. Here, the corresponding hysteresis curves in the small earthquake, medium earthquake, and large earthquake load modes are shown in Figures 5A to 5C. The damper yield load capacity and yield displacement output results of the step-yielding damper are shown in Figure 6. S10. Construct the skeleton curve of the multi-level step-yield damper according to the output hysteresis curves in the three modes. S11, the hysteresis curve of the multi-level step-yield damper is back-calculated by the constructed skeleton curve. S12, Design and install a multi-level step-yield damper.

[0052] Entering the algorithm program in steps S5, S7, and S9 to perform analysis calculation includes step S101 of setting the initial data of the algorithm program, step S102 of randomly selecting the damper yield load capacity and stiffness value of a fixed number of pairs within a preset range, step S103 of inputting the initial data file, step S104 of reading the damper yield load capacity and stiffness value data using python software, step S105 of changing the parameters of the damper set in SAP2000 software, inputting the software SAP2000 to perform analysis, outputting the data of the structural analysis result, and obtaining the result data file. The method includes step S106, step S107 of reading the result data using matlab software, and substituting the result data into a fitness function to calculate the fitness, step S108 of analyzing the vibration control result, outputting the damper yield load-bearing capacity and yield displacement corresponding to the maximum fitness value in the current population, and judging whether the iteration number has been reached, step S109 of obtaining a new initial population through selection, crossover and mutation, and repeating again, and step S110 of terminating the calculation and outputting the damper yield load-bearing capacity, yield displacement and hysteresis curve corresponding to the maximum fitness value calculated by the algorithm program if the iteration number has been reached.

[0053] Specifically, the fitness function in this embodiment is calculated by calculating, after the limit value of the story displacement angle is met, the absolute value of the difference between the maximum story shear force and the top story acceleration after the vibration control structure optimization and the corresponding values ​​of the original structure, and the ratio of the absolute value to the corresponding values ​​of the original structure, expressed as a percentage, as the fitness function, in which the maximum story displacement angle optimization rate is set to be 30% or less, and the maximum story acceleration optimization rate is set to be 15% or less. The fitness calculation formula is as follows:

[0054]

number

[0055] Here, V maxis the maximum inter-story shear force of the optimized structure (kN), V is the maximum inter-story shear force of the original structure (kN), and V i is the inter-story shear force (kN) of each story of the structure, and v2 is the maximum top story acceleration of the optimized structure (mm / s 2 ), and v1 is the acceleration of the top floor of the original structure (mm / s 2 ), where fitness is a fitness function, and α and β are weighting coefficients. In this embodiment, the weighting coefficients are 1.0.

[0056] From the analytical calculation results, as shown in Figs. 7 and 8, it is found that the structural inter-story displacement angle and inter-story acceleration optimized by the optimization design method of the present invention are effectively controlled.

[0057] Although the embodiments of the present invention have been shown and described above, those skilled in the art can make various changes, modifications, substitutions, and variations to these embodiments without departing from the principles and spirit of the present invention, and the scope of the present invention is limited by the appended claims and their equivalents.

Claims

1. A method for optimizing design of a vibration control structure based on a multi-level step-yield damper, comprising the steps of: Step S1: inputting an uncontrolled vibration-damping structural model into SAP2000 finite element software; Step S2 of determining the installation position of the damper and a structural optimization prediction; Step S3 of determining and setting preset parameters of the algorithm program; Step S4 of setting a small earthquake load mode; Step S5: enter the algorithm program, perform analytical calculations, and output the damper yield load capacity, yield displacement, and hysteresis curve corresponding to the maximum fitness value in the small earthquake load mode; Step S6 of setting a moderate earthquake load mode; Step S7: enter the algorithm program, perform analytical calculations, and output the damper yield load capacity, yield displacement, and hysteresis curve corresponding to the maximum fitness value in the moderate earthquake load mode; Step S8 of setting a large earthquake load mode; Step S9: enter the algorithm program, perform analytical calculations, and output the damper yield load capacity, yield displacement, and hysteresis curve corresponding to the maximum fitness value in the large earthquake load mode; Step S10: constructing a skeleton curve of a multi-level step-yield damper based on the output hysteresis curves in the three modes; Step S11: Back-calculating a hysteresis curve of the multi-level step-yield damper based on the constructed skeleton curve; Step S12: designing and installing a multi-level graduated yield damper; Including, The steps to enter the algorithm program and perform analytical calculations are: Step S101 of setting initial data for an algorithm program; Step S102: randomly selecting damper yield load capacity and stiffness values ​​for a fixed number of pairs within a preset range; Step S103 of inputting an initial data file; Step S104 of reading the damper yield load capacity and stiffness value data using Python (registered trademark) software; Step S105 of modifying the parameters of the damper installed in the SAP2000 software; Step S106: executing the analysis using software SAP2000, outputting the structural analysis result, and obtaining a result data file; Step S107 of reading the result data using Matlab (registered trademark) software and substituting the result data into a fitness function to calculate the fitness; Step S108: analyzing the vibration control result, outputting the damper yield load capacity and yield displacement corresponding to the maximum fitness value in the current population, and judging whether the number of iterations has been reached; If the number of iterations has not been reached, a step S109 of obtaining a new initial population through selection, crossover, and mutation, and repeating the process again; Step S110: when the number of iterations is reached, the calculation is terminated and the damper yield load capacity, yield displacement, and hysteresis curve corresponding to the maximum fitness value calculated by the algorithm program are output; Including, In step S107, the ratio of the absolute value of the difference between the maximum inter-story shear force and the top story acceleration after the vibration control structure optimization and the corresponding values ​​of the original structure to the corresponding values ​​of the original structure is calculated as a percentage, and then the square root of the sum of squares is taken to obtain a fitness calculation formula. [0010] Here, V ｍａｘ is the maximum interlaminar shear force (kN) of the optimized structure, V is the maximum interlaminar shear force (kN) of the original structure, and V ｉ is the inter-story shear force (kN) of each layer of the structure, and v ２ is the maximum top floor acceleration of the optimized structure (mm / s ２ ) and v １ is the acceleration of the top floor of the original structure (mm / s ２ ), where fitness is a fitness function, In step S107, the ratio of the absolute value of the difference between the maximum inter-story shear force and the top story acceleration after the vibration control structure optimization and the corresponding values ​​of the original structure to the corresponding values ​​of the original structure is calculated as a percentage, and the ratio is used as a fitness calculation formula. [0025] Here, V ｍａｘ is the maximum interlaminar shear force (kN) of the optimized structure, V is the maximum interlaminar shear force (kN) of the original structure, and V ｉ is the inter-story shear force (kN) of each layer of the structure, and v ２ is the maximum top floor acceleration of the optimized structure (mm / s ２ ) and v １ is the acceleration of the top floor of the original structure (mm / s ２ ), where fitness is a fitness function, and α and β are weighting coefficients. In step S107, the ratio of the absolute value of the difference between the maximum inter-story displacement angle after the vibration-damping structure optimization and the corresponding value of the original structure to the corresponding value of the original structure is calculated as a percentage, and this is used as a fitness calculation formula. [0030] Here, θ ｍａｘ is the maximum inter-story drift angle of the optimized structure, θ is the maximum inter-story drift angle of the original structure, and fitness is a fitness function. In step S107, the ratio of the absolute value of the difference between the maximum inter-story shear force and the top story acceleration after the vibration control structure optimization and the corresponding values ​​of the original structure to the corresponding values ​​of the original structure is calculated as a percentage, and the ratio is used as a fitness calculation formula. [0045] Here, θ ｍａｘ is the maximum inter-story drift angle of the optimized structure, θ is the maximum inter-story drift angle of the original structure, and a ｍａｘ is the maximum acceleration of the optimized structure, a is the maximum inter-story acceleration of the original structure, and u ｍａｘ is the maximum interlayer displacement of the optimized structure, u is the maximum interlayer displacement of the original structure, α, β, and χ are weighting coefficients, and fitness is a fitness function. A method comprising:

2. In step S3, the preset parameters of the algorithm program include: initial population number, population random range interval, binary coding length, seismic wave and seismic mode of vibration control structure model, and structure optimization prediction; The optimization design method for vibration control structure based on multi-level graduated yield damper as claimed in claim 1.

3. In step S101, the initial data of the algorithm program includes the number of algorithm iterations, the stiffness and yield load range of the damper, the genetic selection rate, the genetic mutation rate, and the genetic crossover rate. The optimization design method for vibration control structure based on multi-level graduated yield damper as claimed in claim 2.

4. The structural analysis results output in step S106 include the inter-story shear force of each structural story, the top story acceleration, the structural inter-story displacement angle, the inter-story acceleration, the maximum inter-story displacement, the damper yield load capacity, the damper displacement, and the hysteresis curve. The optimization design method for vibration control structure based on multi-level graduated yield damper as claimed in claim 1.

5. The population random range section includes that the upper and lower limits of the preset stiffness of the damper are set to 25% of the maximum inter-story stiffness of the uncontrolled structure from the stiffness of the damper that meets the maximum inter-story displacement angle of the structure required by the standard, and the upper and lower limits of the preset load-bearing capacity of the damper are set to 0 to 25% of the maximum inter-story shear force of the structure. The optimization design method for vibration control structure based on multi-level graduated yield damper as claimed in claim 2.

Images