CN115438546A - High-rise industrial factory building vibration comfort degree calculation method based on pedestrian, equipment and automobile load excitation - Google Patents
High-rise industrial factory building vibration comfort degree calculation method based on pedestrian, equipment and automobile load excitation Download PDFInfo
- Publication number
- CN115438546A CN115438546A CN202211082391.7A CN202211082391A CN115438546A CN 115438546 A CN115438546 A CN 115438546A CN 202211082391 A CN202211082391 A CN 202211082391A CN 115438546 A CN115438546 A CN 115438546A
- Authority
- CN
- China
- Prior art keywords
- finite element
- vibration
- factory building
- pedestrian
- load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for calculating the vibration comfort degree of a high-rise industrial factory building under the excitation of pedestrian, equipment and automobile loads, which comprises the following steps: 1. establishing a finite element model of the overall structure of the high-rise industrial factory building; 2. the simulation of the moving load is realized by utilizing a finite element model, a loading plane is an XY plane, a relation between a loading path and time is given by defining a loading area and a loading speed, and the moving load is realized by assigning pressure intensity; 3. the method comprises the steps that vibration response of a high-rise industrial factory building under pedestrian load is tested and measured on the spot, and a vibration acceleration time-course curve of a test floor is obtained; 4. adjusting a finite element model of the overall structure of the high-rise industrial factory building; 5. establishing transmission rules of vibration responses of the floors of the pedestrian, equipment and automobile under the excitation of three types of pedestrian, equipment and automobile loads at different positions of the same floor and different floors through the adjusted finite element model simulation; 6. evaluating comfort level: and obtaining the comfort evaluation of the finite element model simulated at different positions of the same floor and different floors.
Description
Technical Field
The invention discloses a vibration comfort degree calculation method, in particular to a high-rise industrial factory building vibration comfort degree calculation method based on pedestrian, equipment and automobile load excitation.
Background
Comfort refers to the overall evaluation of the satisfaction that a person experiences from both physiological and psychological aspects of an objective environment.
In the prior art, in the field of buildings, in a common vibration comfort degree calculation method under artificial load, vibration source loads are simple harmonic loads, a calculation model is a single-degree-of-freedom motion equation, and finally, a floor vibration peak value acceleration value is obtained through simplification and is used as an index for judging the floor vibration comfort degree. The computational model of prior art can be fine solution people cause the vibration comfort level problem of load floor down, nevertheless along with falling to the ground of "industry upstairs" policy, the appearance of more and more high-rise industrial factory building, if: high and new medicine industry garden factory building (mainly relate to equipment vibration), new energy automobile garden factory building (mainly relate to car vibration), therefore, the excitation source that these different factory buildings exist has very big difference with single artificial load excitation, and there is the limitation in the application of current theoretical calculation method. Because the vibration excitation which can appear in actual use is lacked before the factory building is put into use formally, the comfort of the factory building can not be verified by detection methods such as in-situ test and the like; the prior art is analyzing the comfort of the floor elements, while the comfort analysis of the overall structure is still blank.
Disclosure of Invention
Aiming at the defect that the calculation method of the vibration comfort degree under the excitation of the single man-induced load in the prior art is difficult to meet the analysis requirement of the whole structure comfort degree under the excitation of loads of multiple excitation sources, the invention provides the calculation method of the vibration comfort degree of the high-rise industrial factory building under the excitation of loads of pedestrians, equipment and automobiles.
The technical scheme adopted by the invention for solving the technical problems is as follows: a high-rise industrial factory building vibration comfort degree calculation method based on pedestrian, equipment and automobile load excitation comprises the following steps:
s1, modeling: establishing a finite element model of the overall structure of the high-rise industrial factory building;
s2, simulation of moving load: the simulation of the moving load is realized by utilizing a finite element model, a loading plane is an XY plane, a relation between a loading path and time is given by defining a loading area and a loading speed, and the moving load is realized by assigning pressure intensity;
s3, testing the pedestrian load on the spot: the method comprises the steps that vibration response of a high-rise industrial factory building under pedestrian load is tested and measured on the spot, and a vibration acceleration time-course curve of a test floor is obtained;
s4, adjusting a finite element model: adjusting a finite element model of the overall structure of the high-rise industrial factory building according to the test result of the step S3;
s5, establishing a transfer rule: establishing transmission rules of vibration responses of the three types of floors, namely pedestrian, equipment and automobile loads, on different positions of the same floor and different floors through the adjusted finite element model simulation;
step S6, comfort evaluation: based on the load transfer rule of the high-rise industrial factory building under the action of different loads, a finite element model is obtained to simulate and establish the evaluation of the comfort degree of the floor at different positions of the same floor and different floors by the vibration response of the floor under the excitation of the pedestrian load, the equipment load and the automobile load.
The technical scheme adopted by the invention for solving the technical problem further comprises the following steps:
in the step S1, a finite element model of the overall structure of the high-rise industrial factory building is established by using professional finite element analysis software ABAQUS.
When the ABAQUS is used for establishing the finite element model of the integral structure of the high-rise industrial factory building, wherein the steel pipe concrete is used for modeling, the steel pipe adopts a shell unit, the material property assignment selects from bottom to top, the concrete adopts a solid unit, the material property assignment selects from a neutral axis to two ends, the contact between the steel pipe and the concrete is set to be hard contact and friction contact, and the penalty is set to be 0.4; the steel beams and the structural columns adopt beam units in line units; the floor slab and the shear wall adopt solid units; the steel tube and concrete constraint in the steel tube concrete is set as a shell unit and an entity unit constraint, and the rest of the steel tube concrete involved in contact are set as binding constraints; the bottommost boundary condition is consolidation.
And in the step S2, carrying out a loading path test on the floor slab according to a rectangular short edge midspan mode.
In the step S3, the vibration test sensor is disposed in the middle of the entire floor slab during field measurement, and if the floor slab is set in a partition, the vibration test sensor is set in the middle of each floor slab in each partial partition.
In the step S3, during the field measurement, the test path is set along the transverse midspan of the floor slab, and the path and the motion parameters of the field measurement are the same as the parameters set in the step S2.
In the step S4, the same load is set according to the same walking path, the same sensor position is selected to extract the floor vibration acceleration, the difference between the actual measurement result and the finite element simulation result is analyzed, and if the error is less than 10%, the error is considered to be free; if the error is more than or equal to 10%, adjusting the constitutive relation, the contact setting, the constraint setting, the material attribute and the stress surface relation of the finite element model until obtaining a finite element numerical simulation result which is identical with the actual measurement result, and then considering that the established finite element model has reliability.
The beneficial effects of the invention are: the method is used for establishing a load model based on a large amount of actual measurement and numerical simulation data, considers man-induced loads, equipment loads and automobile loads, includes several common load types of high-rise industrial plants, and has better applicability; the method of combining actual measurement with numerical simulation is adopted, and the vibration comfort degree of the structure is analyzed more accurately at lower cost; the invention also considers the transmission rule of the vibration excitation source and divides the structure into horizontal and vertical functional subareas from the integral angle; the invention establishes a refined model of the high-rise industrial factory building, and can calculate the vibration response of the high-rise industrial factory building more accurately and reasonably, thereby comprehensively analyzing the vibration comfort degree of the whole structure.
The method for calculating the vibration comfort degree of the high-rise industrial factory building is provided based on the man-caused load, the equipment load and the automobile load, overcomes the defect that the excitation source of the high-rise industrial factory building is only the man-caused load at present under a single working condition, takes the equipment load and the automobile load which are required to be considered by the high-rise industrial factory building into consideration at present, and enlarges the applicability of the method for calculating the vibration comfort degree of the high-rise industrial factory building; because vibration excitation in the actual use process cannot be generated before the plant is put into use formally, the plant comfort can not be verified by detection methods such as in-situ tests and the like, the comfort of the plant structure can be obtained more accurately at lower cost by adopting the method of combining the actual measurement with the numerical simulation; by adopting the method, the comfort condition can be analyzed from the whole structure, and the horizontal and vertical functional partitions of the structure can be more reasonably planned; the simulation numerical calculation adopted by the method can accurately obtain the vibration comfort level index of the high-rise industrial factory building, so that whether the comfort level limit value is met or not is judged.
The invention will be further described with reference to the accompanying drawings and specific embodiments.
Drawings
FIG. 1 is a schematic diagram of a finite element structural model according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a finite element model moving load DLOAD subroutine code written in an embodiment of the present invention;
FIG. 3 is a schematic diagram of a sensor layout area according to an embodiment of the present invention;
FIG. 4 is a time-course curve diagram of floor vibration data of area A under actual measurement of pedestrian load in the embodiment of the invention;
FIG. 5 is a time-course curve diagram of floor vibration data of a B area under actual measurement of pedestrian load in the embodiment of the invention;
FIG. 6 is a time-course curve diagram of floor vibration data of area A under finite element simulation pedestrian load in the embodiment of the invention;
FIG. 7 is a time-course curve diagram of floor vibration data of a B area under finite element simulation pedestrian load in the embodiment of the invention;
FIG. 8 is a time-course curve diagram of floor vibration data of the area A under three loads of finite element simulation in the embodiment of the invention;
FIG. 9 is a time-course curve diagram of B-zone floor vibration data under three loads of finite element simulation in the embodiment of the invention;
FIG. 10 is a time-course curve diagram of vibration data of a next floor area A under three loads of finite element simulation in the embodiment of the invention;
FIG. 11 is a flow chart illustrating an embodiment of the present invention.
Fig. 12 is a screenshot of the program for setting the "load module" in the embodiment of the present invention.
Fig. 13 is a screenshot of the program for setting the amplitude curve option in the embodiment of the present invention.
Fig. 14 is a screenshot of the program for setting the parameter of the amplitude curve according to the embodiment of the present invention.
Detailed Description
The present embodiment is a preferred embodiment of the present invention, and other principles and basic structures that are the same as or similar to the present embodiment are within the scope of the present invention.
The invention aims to provide a method for calculating the vibration comfort level of a high-rise industrial factory building under the action of various load excitations (pedestrian, equipment and automobile loads), which makes up the defects of the current building vibration comfort level analysis on the equipment load and the automobile load excitation, enables a numerical simulation result to be more accurate by combining an actual measurement method with a finite element simulation method, and fills the gap of the whole structure vibration comfort level analysis in the current field.
The invention relates to a method for calculating the vibration comfort degree of a high-rise industrial factory building under the excitation of pedestrian, equipment and automobile loads, which comprises the following steps of:
s1, modeling: and establishing a finite element model of the overall structure of the high-rise industrial factory building, wherein in the embodiment, a professional finite element analysis software ABAQUS is used for establishing the finite element model of the overall structure of the high-rise industrial factory building.
The finite element model is a calculation model for performing a finite element analysis, which provides all necessary raw data for the finite element calculation. The process of establishing a finite element model is called finite element modeling, which is the key of the whole finite element analysis process, whether the model is reasonable or not directly influences the precision of a calculation result, the length of calculation time, the size of storage capacity and the completeness of the calculation process, the task of the finite element modeling is to abstract an actual problem or a design scheme into a finite element model which can provide all input data for numerical calculation, and the model quantitatively reflects the characteristics of various aspects of the geometry, materials, loads, constraints and the like of an analysis object.
In this example, when ABAQUS is used to build a finite element model of the overall structure of a high-rise industrial factory building (see FIG. 1), wherein concrete-filled steel tubes are modeled, steel tubes are modeled by "shell elements" (material property assignment choice "from bottom to top"), concrete is modeled by "solid elements" (material property assignment choice "from neutral axis to both ends"), contact between steel tubes and concrete is set as "hard contact" and "frictional contact", (in this example, both contact methods are "hard contact" refers to a contact type in which the contact surface between concrete and steel tubes is vertical, "frictional contact" refers to a contact type in which the contact surface between concrete and steel tubes is tangential, both contacts exist, similarly, when a box is dragged on the ground, the ground is hard contact against the gravity of the box, and the ground is frictional contact with the box, and "penalty" is set as 0.4 "(in this example, the penalty" is a special call in this software, and means the coefficient of friction, because the magnitude of the coefficient of friction between steel tubes and concrete determines the degree of frictional contact, and the general "value of the steel tube concrete" is set as 0.4); the steel beams and the structural columns adopt beam units in line units (as the modeling of components in ABAQUS software adopted in the embodiment has three unit forms, namely solid units, line units (the line units are divided into the beam units and truss units) and shell units, the modeling units selected by the components are only explained here, the assignment is certain, the concrete is endowed with three material properties of density, elasticity and concrete damage plasticity, the steel pipes are endowed with three material properties of density, elasticity and plasticity, and the assignment can be carried out according to a conventional finite element model); the floor and the shear wall adopt solid units (the concrete endows three material properties of density, elasticity and concrete damage plasticity), the configured steel bars endow three material properties of density, elasticity and plasticity, and the assignment is carried out according to a conventional finite element model); the steel tube and concrete constraint in the steel tube concrete is set as a shell unit and an entity unit constraint, and the rest of the steel tube concrete involved in contact are set as binding constraints; the bottommost boundary condition is consolidation.
S2, simulation of moving load: in this embodiment, the DLOAD subprogram is written by using Visual Studio based on the secondary development interface of ABAQUS to realize simulation of the mobile load, and other manners may be adopted in specific implementation.
In this embodiment, when a DLOAD subroutine is written by using Visual Studio based on an ABAQUS secondary development interface to realize simulation of a moving load (see fig. 2), a loading plane is an XY plane, a relational expression of a loading path and time is given by defining a loading area and a loading speed, and the moving load is realized by assigning pressure.
In this embodiment, a specific moving load subroutine is taken as an example, and the specific steps are as follows:
in this embodiment, the move load subroutine includes: the first 5 lines of code are DLOAD default frames;
the code in line 6 is given x and y position variables, v speed variable and s path variable;
the code in line 7 gives a speed value, which in this embodiment is the average speed of normal walking of 1.5m/s (the reason why the value is 1500 instead of 1.5 is that the unit system is mm selected according to the ABAQUS modeling), and in specific implementation, other values can be selected;
line 8 codes are journey = speed multiplied by time;
the association between the values with codes x and y in lines 9 and 10 and the finite element model, specifically the lower IF cycle;
the code in lines 11-20 is an IF loop, specifically, the area of mobile load loading is given first, in this embodiment IF the y-coordinate is between 400-500mm or the y-coordinate is between 900-1000mm (the foot face width is assumed to be 100mm, so 400-500mm is one foot and 900-1000mm is the other foot, where only two feet are set to walk in these two lines), the x-coordinate is between s and s-420 (mentioned in the foregoing code,
v＝1500
s＝v*TIME(2)
v is 1.5m/s 2 And t time is an independent variable in the code) (in the embodiment, the length of the human foot is 42cm, and other normal symbols can be selected during the specific implementationThe length of a fit foot), the pressure is 24KPa (the pressure is the pressure to the ground when a normal person walks, 0.24 is determined by mm and N according to the unit system selected by ABAQUS modeling, and other values can be selected during specific implementation); if the x-coordinate or the y-coordinate is not within this range, the pressure is 0.
In this embodiment, the simulated floor slab is subjected to a loading path test in a manner that the short sides of the simulated floor slab span a center, that is, the path is arranged parallel to the long sides of the rectangle and at the center line position of the short sides (since the real building includes the simulated building, neither of which is a strict rectangle, and the rectangle in the present invention is a shape similar to the rectangle).
S3, testing the pedestrian load on the spot: and (3) testing the vibration response of the high-rise industrial factory building under the pedestrian load on the spot, and measuring the vibration response on the spot to obtain a vibration acceleration time-course curve of the test floor as shown in figures 4 and 5.
In this embodiment, go to and measure on the spot, the vibration test sensor is arranged and is seen in fig. 3, the vibration test sensor sets up the principle: arranged in the middle of the whole floor, so that the deflection is larger and the vibration response is theoretically maximum, if the floor is arranged in a subarea, the middle of each part of the floor is taken in each subarea.
In this embodiment, the floor slab is divided into two areas, namely, a floor slab and a floor slab B, between the two areas (in this embodiment, since the floor slab is formed by two areas, namely, the two areas a and B, the floor slab is divided into two areas a and B, the current two areas are the middle of the two partial floor slabs, the deflection is the largest theoretically, the vibration response is the largest, the two positions are selected, and the two areas are selected for testing whether the comfort level is satisfied, namely, whether the vibration response is satisfied, so that the two areas are selected at the most unfavorable position, namely, the position with the largest deflection, and the two areas must be selected from this angle.
In the present embodiment, the test is performed by a pedestrian in the field, and the walking speed of the pedestrian is about 1.5m/s in the test. The time course curve of the floor vibration acceleration tested and recorded in the field is shown in figures 4 and 5.
S4, adjusting a finite element model: and adjusting the finite element model of the overall structure of the high-rise industrial factory building by taking the test results of the figure 4 and the figure 5 as the guide.
In this embodiment, the test results shown in fig. 4 and 5 are taken as guidance, the same load is set according to the same walking path, the same sensor position is selected to extract the floor vibration acceleration, the difference between the actual measurement result and the finite element simulation result is analyzed (in this embodiment, the peak acceleration obtained by the actual measurement and the finite element simulation is analyzed, and how large the error is compared directly), if the error is less than 10%, the error is considered to be free, and if the error is greater than or equal to 10%, the constitutive relation, the contact setting, the constraint setting, the material property and the stress surface relation of the finite element model are adjusted, and the adjustment can be performed according to the conventional adjustment mode of the finite element model (for example: 1% of the actual measured results are obtained each time, and the established finite element model is considered to have reliability through comparison of multiple measurements and actual measurements, wherein the adjustment is a series of complex processes, for example, the stress-strain relationship (constitutive relationship) of concrete and steel pipe materials of input software can directly influence the finite element simulation result, the contact arrangement methods among beam plates, beam columns and the like are multiple, the plate beams (four sides are consolidated or two sides are consolidated and two sides are hinged and the like), the constraint arrangement among beam columns (whether the consolidation or the hinging is carried out) is multiple, and the rigidity of the composite floor slab can not be directly obtained, the rigidity of the composite floor slab can only be converted through theoretical calculation and then considered to be equivalent to the concrete slab, because the finite element model cannot be completely matched with the actual structure, the adjustment is close to the actual measured results according to the actual measured results, namely, the error is less than 10%, which is a common method for model adjustment, and the established finite element model is considered to have reliability through comparison of the multiple measurements and the actual measurements and the simulation, (because of the limitation of test conditions (it is difficult to have perfect conditions like theoretical analysis, and there are environmental interference or human error, machine error, etc.), finite element simulation is unlikely to be the same as the actual result, so the comparison between the two is that the allowable error exists, and when the error is within 10%, the finite element model is considered to have reliability), the result simulated by the model has qualitative significance.
S5, establishing a transfer rule: and establishing the transmission rule of the vibration response of the floor under the excitation of three types of pedestrian, equipment and automobile loads on the same layer and different layers through the adjusted finite element model simulation.
In this embodiment, when the vibration response transmission rule of the floor under the excitation of three types of excitation of pedestrian, equipment and automobile load is simulated by the finite element model, the setting of pedestrian and automobile load (the code in the step S2 can be said to be a code of moving load, so the code is universal, and only the ranges of v (speed), F (pressure intensity, determined by mass) and x and y in the code need to be changed according to specific conditions, the length and width of the foot in the above embodiment can be changed into the length and width of the tire of the vehicle) can be realized by changing the speed and the length and width of the tire through the DLOAD subprogram; the simulation of the equipment vibration is realized by setting the load as a dynamic load, giving a vibration amplitude value and specifically setting according to the vibration frequency of different machines.
The following embodiment provides a specific setting mode, in this embodiment, in the "load module" of the software, setting "concentrated force", as shown in fig. 12, cf3 is Z-axis direction, 10 refers to force 10KN downward, and "Amplitude" is obtained by setting a periodically varying Amplitude curve (Periodic), as shown in fig. 13, specifically setting fig. 14, and "time span" is obtained by starting from a certain analysis step or all analyses, in this embodiment, for all analyses, "Circular frequency" is circle frequency, and the lower two rows are starting time and initial Amplitude, setting to 0, and a and B of the table are set to 5 and 0, and different machines may be set differently.
Fig. 8 is a time-course curve of vibration acceleration of a floor slab in an area a under the action of the area a by three excitations, fig. 9 is a time-course curve of vibration acceleration of a floor slab in an area B under the action of the area a by three excitations, and fig. 10 is a time-course curve of vibration acceleration of a floor slab in an area a under the action of the area a by three excitations. Drawing (A)8, comparing with the figure 9, the transmission rule of the vibration response in the same layer can be obtained, and the transmission rule is related to the distance between the areas A and B and the structural material and the type; FIG. 8 is a graph comparing the transmission law of vibration response at different layers with FIG. 10, and the transmission law is related to the layer height, the floor structure material and the type (comparing the acceleration peak value, FIG. 8 is the A area floor acceleration time course curve of the vibration source at the area A, the peak value is 39.35mm/s 2 FIG. 9 is the time curve of the acceleration of the floor slab in the area B with the peak value of 6.91mm/s 2 The energy loss of the vibration transmitted from the area A to the area B can be seen, a relational expression of the vibration along with the change of the distance can be obtained by simulating the floor slabs with different distances for many times, and the rule can reasonably perform horizontal plane function partition, for example, an office and a production operation area can be arranged accordingly, so that the vibration of the production operation area is transmitted to the office for a few times; the same method is adopted to research the vibration attenuation of the upper layer and the lower layer, the influence of the vibration of the upper layer on each part of the lower layer can be obtained, a relational expression of the vibration changing along with the height of the layer can be obtained, and the method can also be used for functional division setting in the vertical direction, for example, the vibration of the 4 layers is transmitted to the 3 layers to attenuate the remaining half and is transmitted to the 1 layer to be remained a few, and then the 1 layer can be considered to be completely set as an office; these two parts of work require multiple tests to be compared with finite element simulations, which are not yet fully validated and proposed, and are considered feasible. )
Step S6, comfort evaluation: based on comfort evaluation of high-rise industrial factory building under the different load effect, in this embodiment, based on the finite element model that has been perfected, the whole vibration condition of high-rise industrial factory building under the effect of three kinds of load of simulation can be got whether the holistic comfort level of structure satisfies standard requirement and the whole level of structure to and vertical functional partitioning suggestion, if: the whole industrial plant is suitable for placing machine equipment, working areas and transport vehicle passages.
In this embodiment, according to the above-mentioned perfect finite element model, under the effect of various loads, more accurate evaluation can be obtained from a plurality of angle simulations based on comfort level standard, including:
(1) The horizontal and vertical vibration peak acceleration of the whole structure;
(2) The position with larger vibration response of each floor;
(3) The vibration transmission response rule of the same layer;
(4) The vibration response transmission rules of different layers;
whether the whole structure meets the standard comfort level can be evaluated by (1); the vibration sensitive position of the floor can be obtained by the step (2), so that the horizontal function partition is carried out or proper vibration stopping measures are taken; the reduction rule of the vibration response on the same layer can be obtained through the step (3), so that the function partition can be better carried out in the horizontal direction; the reduction rule of the vibration response in the adjacent layers can be obtained through the step (4), so that the functional partition can be better performed in the vertical direction.
The method can take the automobile load and the equipment load into consideration, and expand the applicability of the vibration comfort degree calculation method so as to supplement the singleness of the artificial load of the excitation source in the prior art; the method combining actual measurement with numerical simulation is adopted to analyze accurately with lower cost, has economy, and can perfectly solve the problems that the prior art lacks vibration excitation which can appear in actual use before the factory building is put into use formally and can not carry out in-situ test and other detection methods to verify the comfort of the factory building; according to the invention, by establishing the integral numerical model, the horizontal and vertical vibration transmission rules of the integral structure can be better analyzed, so that reference is provided for division of the horizontal and vertical functional partitions of the integral structure, the problem of vibration comfort is avoided, and the blank of analyzing the comfort of the integral structure is filled.
The method is used for establishing a load model based on a large amount of actual measurement and numerical simulation data, considers man-induced loads, equipment loads and automobile loads, includes several common load types of high-rise industrial plants, and has better applicability; the method of combining actual measurement with numerical simulation is adopted, so that the vibration comfort degree of the structure is more accurately analyzed at lower cost; the invention also considers the transmission rule of the vibration excitation source and divides the structure into horizontal and vertical functional subareas from the integral angle; the invention establishes a refined model of the high-rise industrial factory building, and can calculate the vibration response of the high-rise industrial factory building more accurately and reasonably, thereby comprehensively analyzing the vibration comfort degree of the whole structure.
The method for calculating the vibration comfort level of the high-rise industrial factory building is provided based on the man-induced load, the equipment load and the automobile load, overcomes the defect that the excitation source of the high-rise industrial factory building is only the man-induced load at present under the single working condition, considers the equipment load and the automobile load which are considered by the high-rise industrial factory building at present, and enlarges the applicability of the method for calculating the vibration comfort level of the high-rise industrial factory building; because vibration excitation in the actual use process cannot be generated before the factory building is put into use formally, the comfort of the factory building cannot be verified by detection methods such as an in-situ test and the like, the comfort of the factory building structure can be obtained more accurately at lower cost by adopting the method of combining actual measurement with numerical simulation of the invention; by adopting the method, the comfort condition can be analyzed from the whole structure, and the horizontal and vertical functional partitions of the structure can be more reasonably planned; the simulation numerical calculation adopted by the method can accurately obtain the vibration comfort level index of the high-rise industrial factory building, so that whether the comfort level limit value is met or not is judged.
Claims (7)
1. A method for calculating the vibration comfort degree of a high-rise industrial factory building based on pedestrian, equipment and automobile load excitation is characterized by comprising the following steps: the method comprises the following steps:
s1, modeling: establishing a finite element model of the overall structure of the high-rise industrial factory building;
s2, simulation of the moving load: the simulation of the moving load is realized by utilizing a finite element model, a loading plane is an XY plane, a relational expression of a loading path and time is given by defining a loading area and speed, and the moving load is realized by assigning pressure intensity;
s3, testing the pedestrian load on the spot: the method comprises the steps that vibration response of a high-rise industrial factory building under pedestrian load is tested and measured on the spot, and a vibration acceleration time-course curve of a test floor is obtained;
s4, adjusting a finite element model: adjusting a finite element model of the overall structure of the high-rise industrial factory building according to the test result of the step S3;
s5, establishing a transfer rule: establishing transmission rules of vibration responses of the floors of the pedestrian, equipment and automobile under the excitation of three types of pedestrian, equipment and automobile loads at different positions of the same floor and different floors through the adjusted finite element model simulation;
step S6, comfort evaluation: based on the load transfer rule of the high-rise industrial factory building under the action of different loads, a finite element model is obtained to simulate and establish the evaluation of the vibration response of the floor under the excitation of three types of pedestrian load, equipment load and automobile load at different positions of the same floor and different floors.
2. The method for calculating the vibration comfort level of the high-rise industrial factory building under the excitation of pedestrian, equipment and automobile loads according to claim 1, which is characterized in that: in the step S1, a finite element model of the overall structure of the high-rise industrial factory building is established by using professional finite element analysis software ABAQUS.
3. The method for calculating the vibration comfort level of the high-rise industrial factory building under the excitation of pedestrian, equipment and automobile loads according to claim 2, wherein the method comprises the following steps: when the ABAQUS is used for establishing the finite element model of the integral structure of the high-rise industrial factory building, wherein the steel pipe concrete is used for modeling, the steel pipe adopts a shell unit, the material property assignment selects from bottom to top, the concrete adopts a solid unit, the material property assignment selects from a neutral axis to two ends, the contact between the steel pipe and the concrete is set to be hard contact and friction contact, and the penalty is set to be 0.4; the steel beams and the structural columns adopt beam units in line units; the floor slab and the shear wall adopt solid units; the steel tube and concrete constraint in the steel tube concrete is set as a shell unit and an entity unit constraint, and the rest of the steel tube concrete involved in contact are set as binding constraints; the bottommost boundary condition is consolidation.
4. The method for calculating the vibration comfort level of the high-rise industrial factory building under the excitation of the pedestrian load, the equipment load and the automobile load according to claim 1 is characterized in that: and in the step S2, carrying out a loading path test on the floor slab according to a rectangular short edge midspan mode.
5. The method for calculating the vibration comfort level of the high-rise industrial factory building under the excitation of pedestrian, equipment and automobile loads according to claim 1, which is characterized in that: in the step S3, the vibration test sensor is arranged in the middle of the whole floor slab during field measurement, and if the floor slab is arranged in a subarea, the vibration test sensor is arranged in the middle of each part of the subarea where the floor slab is arranged.
6. The method for calculating the vibration comfort level of the high-rise industrial factory building under the excitation of pedestrian, equipment and automobile loads according to claim 1, which is characterized in that: in the step S3, during the field measurement, the test path is set along the transverse midspan of the floor slab, and the path and the motion parameters of the field measurement are the same as the parameters set in the step S2.
7. The method for calculating the vibration comfort level of the high-rise industrial factory building under the excitation of pedestrian, equipment and automobile loads according to claim 1, which is characterized in that: in the step S4, the same load is set according to the same walking path, the same sensor position is selected to extract the floor vibration acceleration, the difference between the actual measurement result and the finite element simulation result is analyzed, and if the error is less than 10%, the error is considered to be free; and if the error is more than or equal to 10%, adjusting the constitutive relation, the contact setting, the constraint setting, the material attribute and the stress surface relation of the finite element model until a finite element numerical simulation result which is identical with the actual measurement result is obtained, and determining that the established finite element model has reliability.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211082391.7A CN115438546A (en) | 2022-09-06 | 2022-09-06 | High-rise industrial factory building vibration comfort degree calculation method based on pedestrian, equipment and automobile load excitation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211082391.7A CN115438546A (en) | 2022-09-06 | 2022-09-06 | High-rise industrial factory building vibration comfort degree calculation method based on pedestrian, equipment and automobile load excitation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115438546A true CN115438546A (en) | 2022-12-06 |
Family
ID=84247561
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211082391.7A Pending CN115438546A (en) | 2022-09-06 | 2022-09-06 | High-rise industrial factory building vibration comfort degree calculation method based on pedestrian, equipment and automobile load excitation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115438546A (en) |
-
2022
- 2022-09-06 CN CN202211082391.7A patent/CN115438546A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109783961A (en) | A kind of analog detecting method and simulated testing system of the automobile frame fatigue life based on CAE | |
| Lavan et al. | Seismic behavior of viscously damped yielding frames under structural and damping uncertainties | |
| Park et al. | A model updating method with strain measurement from impact test for the safety of steel frame structures | |
| CN116738556B (en) | Modeling method, device and equipment for construction wall and storage medium | |
| Ponsi et al. | A multi-objective optimization approach for FE model updating based on a selection criterion of the preferred Pareto-optimal solution | |
| CN110399661B (en) | Discrete-continuous coupling-based steel bridge deck pavement interlayer shear test simulation method | |
| Taher et al. | Earthquake input and state estimation for buildings using absolute floor accelerations | |
| CN107016194A (en) | Bridge cable system deterioration method for analyzing performance based on multi-scale finite Meta Model | |
| CN117236053A (en) | Test system for measuring displacement of prototype in packaging material in falling process through sensor | |
| CN115422649A (en) | Method for predicting corrosion fatigue life of hanging structure and related equipment | |
| CN115438546A (en) | High-rise industrial factory building vibration comfort degree calculation method based on pedestrian, equipment and automobile load excitation | |
| Nguyen | A deep learning platform for evaluating energy loss parameter in engineering structures | |
| CN103528775B (en) | A kind of structural health detection method based on response sensitivity | |
| CN105912465A (en) | Test method for guidance and navigation software of satellite-bone attitude and orbit control | |
| Raphael et al. | Sensor data interpretation in bridge monitoring—A case study | |
| CN108376196A (en) | It is a kind of to improve existing foot bridge comfort level method | |
| Wąsik et al. | Model‐Based Dynamic Structural Modification of Machine Tools | |
| Yuen | Structural identification of Tsing Ma bridge | |
| Froli et al. | Experimental static and dynamic tests on a large-scale free-form Voronoi grid shell mock-up in comparison with finite-element method results | |
| Ghaffary et al. | Uncertainty quantification and probabilistic performance-based assessment of nonlinear behavior of moment resisting frames under design and extreme wind loads | |
| US20200151293A1 (en) | Systems and methods for the analysis of structural components | |
| Aktas et al. | Computer-Aided Mold Design Algorithm for Precast Concrete Elements. | |
| Ruediger et al. | Dealing with uncertainties in manufacturing process simulations | |
| CN117606960B (en) | Method and system for scale test for simulating tunnel defect evolution mechanism in operation period | |
| Feltrin et al. | Seismic Reliability Maps of Code-Compliant Italian Reinforced-Concrete Bare and Infilled Frame Buildings |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |