Building multi-stage frequency identification method, subsystem and storage medium

Abstract

The application discloses a building multi-order frequency identification method, a subsystem and a storage medium, and belongs to the technical field of building monitoring. The application obtains physical parameters of a subsystem installed on the top of a target building, and obtains vibration data of the subsystem in the vibration process of the target building under the action of general environment in a preset time period; determines contact force data of the subsystem and the target building and structural dynamic characteristic parameters of the target building according to the physical parameters of the subsystem and the vibration data; establishes an equivalent relationship among the contact force data, the structural dynamic characteristic parameters and characteristic parameters corresponding to multi-order frequencies of the target building according to the kinematics principle of a coupling system; and determines the multi-order frequencies according to the equivalent relationship, the contact force data and the dynamic characteristic parameters. The application realizes indirect identification of building multi-order frequencies, and improves the efficiency of building monitoring.

Description

Technical Field

[0001] This invention relates to the field of building monitoring, and in particular to a method, subsystem, and storage medium for identifying multi-order frequencies in buildings. Background Technology

[0002] With continuous economic development, my country has seen a surge in high-rise and super high-rise buildings, necessitating dynamic performance monitoring to ensure their safe operation. Currently, multi-sensor monitoring methods are commonly used, deploying numerous sensors along the floors of high-rise and super high-rise buildings to monitor their dynamic characteristics in real time. However, as building heights and the number of high-rise and super high-rise buildings continue to increase, the required number of sensors has skyrocketed, resulting in high costs and low efficiency for large-scale urban monitoring of high-rise and super high-rise buildings.

[0003] The above content is only used to help understand the technical solution of the present invention and does not represent an admission that the above content is prior art. Summary of the Invention

[0004] The main objective of this invention is to provide a method, subsystem, and storage medium for identifying multi-order frequencies in buildings. This aims to achieve indirect identification of multi-order frequencies in buildings, thereby improving the efficiency of building monitoring.

[0005] To achieve the above objectives, the present invention provides a method for identifying multi-order frequencies of buildings, the method comprising the following steps:

[0006] The physical parameters of the subsystem installed on the top of the target building are obtained, and the vibration data of the subsystem during the vibration process of the target building under normal environmental conditions within a preset time period are obtained.

[0007] The contact force data between the subsystem and the target building and the structural dynamic characteristic parameters of the target building are determined based on the physical parameters and the vibration data.

[0008] Based on the kinematics principle of coupled systems, establish the equivalent relationship between the contact force data, the structural dynamic characteristic parameters, and the characteristic parameters corresponding to the multi-order frequencies of the target building;

[0009] The multi-order frequencies are determined based on the equivalence relationship, the contact force data, and the structural dynamic characteristic parameters.

[0010] Optionally, the subsystem comprises, from bottom to top: a damper, a spring, and a counterweight, and the step of obtaining the physical parameters of the subsystem installed on the top of the target building includes:

[0011] The damping coefficient of the damper, the stiffness coefficient of the spring, and the mass of the counterweight are obtained. The physical parameters include the damping coefficient, the stiffness coefficient, and the mass.

[0012] Optionally, the vibration data includes: first displacement data and first velocity data at the end of the subsystem closer to the target building, and second displacement data and second velocity data at the other end of the subsystem farther from the target building. The step of determining the contact force data between the subsystem and the target building based on the physical parameters and the vibration data includes:

[0013] Determine the displacement difference between the first displacement data and the second displacement data, and determine the velocity difference between the first velocity data and the second velocity data;

[0014] The first elastic force between the subsystem and the target building is calculated based on the displacement difference data and the stiffness coefficient, and the first damping force between the subsystem and the target building is calculated based on the velocity difference data and the damping coefficient.

[0015] The contact force data is determined based on the first elastic force and the first damping force.

[0016] Optionally, the vibration data includes: first displacement data and first velocity data at the end of the subsystem closest to the target building, and second displacement data and second velocity data at the other end of the subsystem furthest from the target building. The step of determining the structural dynamic characteristic parameters of the target building based on the physical parameters and the vibration data includes:

[0017] The second elastic force of the subsystem is determined based on the second displacement data and the stiffness coefficient, and the second damping force of the subsystem is determined based on the second velocity data and the damping coefficient.

[0018] The structural dynamic characteristic parameters are calculated based on the contact force data, the second elastic force, and the second damping force.

[0019] Optionally, the feature parameters include the similarity matrix corresponding to the multi-order frequencies, the equivalence relationship is the first set of equations between the contact force data, the structural dynamic feature parameters, and the similarity matrix, the target matrix includes the multi-order frequencies, and the step of determining the multi-order frequencies based on the equivalence relationship, the contact force data, and the structural dynamic feature parameters includes:

[0020] Using the contact force data and the structural dynamic characteristic parameters as known quantities in the first set of equations, the first set of equations is solved to obtain the similarity matrix.

[0021] Calculate at least one eigenvalue of the similarity matrix;

[0022] The multi-order frequency is determined based on the at least one characteristic value.

[0023] Optionally, the step of establishing the equivalent relationship between the contact force data, the structural dynamic characteristic parameters, and the characteristic parameters corresponding to the multi-order frequencies of the target building based on the kinematic principles of the coupled system includes:

[0024] A mechanical analysis is performed on the target building, and the kinematic equations of the target building are determined based on the first preset degree of freedom data of the target building.

[0025] Based on the kinematic principles of the coupled system, the kinematic equations and the contact force equations between the subsystem and the target building are transformed into the first set of equations.

[0026] Optionally, after determining the multi-order frequency based on the at least one feature value, the method further includes:

[0027] The multi-order frequency is determined as a multi-order verification frequency based on the second set of equations, the contact force data, and the structural dynamic characteristic parameters. The second set of equations is determined based on the kinematic equations corresponding to the second preset degree of freedom data and the kinematic principle of the coupled system. The second preset degree of freedom data is greater than the first preset degree of freedom data.

[0028] When the multi-order frequency meets the preset verification conditions, the multi-order frequency is determined to be correct;

[0029] When the multi-order frequency does not meet the preset verification conditions, it is determined that the multi-order frequency is incorrect;

[0030] The multi-order frequency includes multiple first sub-frequencys corresponding to different orders, the multi-order verification frequency includes multiple second sub-frequencys corresponding to different orders, and the preset verification condition is that the first sub-frequency and the second sub-frequency corresponding to the same order are equal.

[0031] Optionally, after determining the multi-order frequency based on the at least one feature value, the method further includes:

[0032] When the stiffness coefficient is less than the preset stiffness coefficient, or when the damping coefficient is less than the preset damping coefficient, the multi-order frequency is determined to be incorrect.

[0033] When the stiffness coefficient is greater than or equal to the preset stiffness coefficient and the damping coefficient is greater than or equal to the preset damping coefficient, the step of determining the multi-order frequency as the multi-order verification frequency based on the second set of equations, the contact force data, and the structural dynamic characteristic parameters is executed.

[0034] Furthermore, to achieve the above objectives, the present invention also provides a subsystem installed on the top of a target building. The subsystem includes a memory, a processor, and a building multi-frequency identification program stored in the memory and executable on the processor. The building multi-frequency identification program is configured to implement the steps of the building multi-frequency identification method described above.

[0035] In addition, to achieve the above objectives, the present invention also provides a storage medium storing a building multi-frequency identification program, wherein when the building multi-frequency identification program is executed by a processor, it implements the steps of the building multi-frequency identification method described above.

[0036] This invention proposes a method for identifying multiple frequencies of a building, which, compared to current multi-sensor monitoring methods, eliminates the need to install sensors on each floor. The method includes: acquiring physical parameters of a subsystem installed on the top of a target building, and acquiring vibration data generated by the subsystem as the target building vibrates within a preset time period; determining contact force data between the subsystem and the target building, and structural dynamic characteristic parameters of the target building, based on the physical parameters and the vibration data; and processing the contact force data and the structural dynamic characteristic parameters according to the kinematics principle of coupled systems to obtain the multiple frequencies of the target building. Compared to current multi-sensor monitoring methods, this application avoids directly monitoring the vibration of the target building by monitoring vibration data generated by the vibration intensity of the target building. Furthermore, by processing the vibration data and physical parameters using the kinematics principle of coupled systems, the frequency of the target building can be calculated. This allows the frequency of the target building to be determined simply by monitoring the motion of the subsystem, reducing the number of sensors required on each floor and improving the efficiency of detecting multiple frequencies of the target building. Attached Figure Description

[0037] Figure 1 This is a schematic diagram of the structure of a subsystem of the hardware operating environment involved in the embodiments of the present invention;

[0038] Figure 2 This is a flowchart illustrating the first embodiment of the multi-order frequency identification method for buildings according to the present invention.

[0039] Figure 3 This is a flowchart illustrating the third embodiment of the multi-order frequency identification method for buildings according to the present invention.

[0040] Figure 4 This is a flowchart illustrating the fourth embodiment of the multi-order frequency identification method for buildings according to the present invention.

[0041] Figure 5 This is a flowchart illustrating the fifth embodiment of the multi-order frequency identification method for buildings according to the present invention.

[0042] Figure 6 This is a simplified model diagram of the target building in an embodiment of the present invention;

[0043] The realization of the objective, functional features and advantages of the present invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0044] It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.

[0045] Reference Figure 1 , Figure 1 This is a schematic diagram of the structure of a subsystem of the hardware operating environment involved in an embodiment of the present invention.

[0046] like Figure 1 As shown, the subsystem may include: a processor 1001, such as a central processing unit (CPU), a communication bus 1002, an interactive device 1003, a network interface 1004, and a memory 1005. The communication bus 1002 is used to enable communication between these components. The interactive device 1003 may include a display screen or an input unit such as a keyboard. Optionally, the interactive device 1003 may also be connected to the communication bus via a standard wired interface or a wireless interface. The network interface 1004 may optionally include a standard wired interface or a wireless interface (such as a Wi-Fi interface). The memory 1005 may be a high-speed random access memory (RAM) or a stable non-volatile memory (NVM), such as a disk drive. Optionally, the memory 1005 may also be a storage device independent of the aforementioned processor 1001.

[0047] Furthermore, the subsystem may also include a camera device for measuring first displacement data and first velocity data of the subsystem at the end closest to the target building, and second displacement data and second velocity data of the subsystem at the other end furthest from the target building.

[0048] Furthermore, the subsystem comprises, from bottom to top, a damper, a spring, and a counterweight.

[0049] Those skilled in the art will understand that Figure 1 The structure shown does not constitute a limitation on the subsystem and may include more or fewer components than shown, or combine certain components, or have different component arrangements.

[0050] like Figure 1As shown, the memory 1005, which serves as a storage medium, may include an operating system, a data storage module, a network communication module, a user interface module, and a building multi-frequency identification program.

[0051] exist Figure 1 In the subsystem shown, the network interface 1004 is mainly used for data communication with other devices; the interactive device 1003 is mainly used for data interaction with the user; the processor 1001 and the memory 1005 in the subsystem of the present invention can be set in the subsystem. The subsystem calls the building multi-order frequency identification program stored in the memory 1005 through the processor 1001 and executes the building multi-order frequency identification method provided in the embodiment of the present invention.

[0052] This invention provides a method for identifying multi-order frequencies in buildings, referring to... Figure 2 , Figure 2 This is a flowchart illustrating the first embodiment of a multi-order frequency identification method for buildings according to the present invention.

[0053] In this embodiment, the building multi-order frequency identification method includes:

[0054] Step S10: Obtain the physical parameters of the subsystem installed on the top of the target building, and obtain the vibration data of the subsystem during the vibration process of the target building under normal environmental conditions within a preset time period.

[0055] The subsystem can be installed at the center of the top of the target building. The preset time period can generally be set by the user. In this embodiment, the normal environment is a non-typhoon weather environment. The wind in the normal environment acts on the target building and the subsystem as a whole, and the vibration data of the subsystem during the vibration process of the target building under the normal environment is obtained.

[0056] Step S20: Determine the contact force data between the subsystem and the target building and the structural dynamic characteristic parameters of the target building based on the physical parameters and the vibration data;

[0057] The contact force data here refers to a pair of opposite forces between the target building and the subsystem within the preset time period. That is, the contact force exerted by the target building on the subsystem is opposite in direction to the contact force exerted by the subsystem on the target building. The structural dynamic characteristic parameters here are signals composed only of the subsystem's parameters and vibration response, but contain only information about the target building and not the subsystem. The contact force data and structural dynamic characteristic parameters of the target building are calculated based on the physical parameters.

[0058] Step S30: Establish the equivalent relationship between the contact force data, the structural dynamic characteristic parameters, and the characteristic parameters corresponding to the multi-order frequencies of the target building based on the kinematic principle of the coupled system;

[0059] In this embodiment, the mathematical relationship between the contact force data, the structural dynamic characteristic parameters, and the characteristic parameters corresponding to the multi-order frequencies of the target building is established through the kinematic principle of coupled systems. The equivalent relationship is obtained by modifying the mathematical equation set according to the subspace identification method. The subspace identification method is a modal identification method for non-time-varying linear systems based on time-domain data matrix calculation. This method uses state equations and signal equations to describe the system. The subspace identification method can be expressed as:

[0060] X k+1 =AX k +BU k +w k

[0061] Y k =CX k +DU k +v k

[0062] Where A, B, C, and D are time-invariant system matrices, and X k U is the system's state vector. k and Y k These are the system's input and output signals, w k and v k These are system error and environmental noise, respectively. Applying the subspace identification method to analyze the output and input signals of a linear system can effectively identify system matrices A, B, C, and D, and extract corresponding modal features, such as system frequency, mode shape, and damping. The subspace identification method is used to establish equivalent relationships between the contact force data, the structural dynamic characteristic parameters, and the characteristic parameters corresponding to the multi-order frequencies of the target building.

[0063] Step S40: Determine the multi-order frequency based on the equivalent relationship, the contact force data, and the structural dynamic characteristic parameters.

[0064] In this embodiment, the contact force data and the structural dynamic characteristic parameters can be added to the corresponding positions of the equivalent relationship, the matrix corresponding to the multi-order frequency can be calculated according to the equivalent relationship, and the multi-order frequency can be determined according to the matrix corresponding to the multi-order frequency.

[0065] In this embodiment, compared with the current multi-sensor monitoring method, this application avoids directly monitoring the vibration of the target building by monitoring the vibration data generated by the vibration intensity of the target building, and processes the vibration data and physical parameters by coupling system kinematics principle, thereby calculating the frequency of the target building. This allows the frequency of the target building to be determined by monitoring the motion of the subsystem, reducing the number of sensors to be set on each floor and improving the efficiency of monitoring multiple frequencies of the target building.

[0066] Furthermore, based on the first embodiment, a second embodiment of the building multi-order frequency identification method of the present invention is proposed. In this embodiment, the subsystem includes, from bottom to top, a damper, a spring, and a counterweight. The step of obtaining the physical parameters of the subsystem installed on the top of the target building includes:

[0067] The damping coefficient of the damper, the stiffness coefficient of the spring, and the mass of the counterweight are obtained. The physical parameters include the damping coefficient, the stiffness coefficient, and the mass.

[0068] In this embodiment, the physical parameters are specifically determined based on the damping coefficient of the damper in the subsystem, the stiffness coefficient of the spring, and the mass of the counterweight. The counterweight mass is generally set according to the mass of the target building; for example, the counterweight mass is set to 0.1% of the target building's mass. In other embodiments, the mass of the subsystem can be determined by measuring the stiffness coefficient of the subsystem using a stiffness coefficient test method, measuring the damping coefficient of the subsystem using a damping coefficient test method, and using a weighing device.

[0069] In this embodiment, by sequentially arranging dampers, springs, and counterweights in the subsystem from bottom to top, the physical parameters of the subsystem can be determined based on the parameters of the dampers, springs, and counterweights. By replacing the dampers, springs, and counterweights corresponding to the target parameters, the submodal characteristic parameters of the subsystem can be easily adjusted, thereby improving the efficiency of monitoring the multi-order frequencies of the target building.

[0070] Furthermore, based on the second embodiment, a third embodiment of the building multi-order frequency identification method of the present invention is proposed. In this embodiment, reference is made to... Figure 3 The vibration data includes: first displacement data and first velocity data at the end of the subsystem closest to the target building, and second displacement data and second velocity data at the other end of the subsystem furthest from the target building. The contact force data between the subsystem and the target building is determined based on the physical parameters and the vibration data, including:

[0071] Step S21: Determine the displacement difference data between the first displacement data and the second displacement data, and determine the velocity difference data between the first velocity data and the second velocity data;

[0072] In this embodiment, the displacement difference data is obtained by subtracting the first displacement data from the second displacement data. Specifically, the second displacement data and the first displacement data each include second displacement sub-data and first displacement sub-data. The first displacement sub-data is subtracted from the second displacement sub-data at the same moment to obtain multiple displacement difference sub-data, and the displacement difference data is determined based on these multiple displacement difference sub-data. Similarly, the velocity difference data is obtained by subtracting the first velocity data from the second velocity data. Specifically, the second velocity data and the first velocity data each include second velocity sub-data and first velocity sub-data. The first velocity sub-data is subtracted from the second velocity sub-data at the same moment to obtain multiple velocity difference sub-data, and the velocity difference data is determined based on these multiple velocity difference sub-data.

[0073] Step S22: Calculate the first elastic force between the subsystem and the target building based on the displacement difference data and the stiffness coefficient, and calculate the first damping force between the subsystem and the target building based on the velocity difference data and the damping coefficient;

[0074] The first elastic force here refers to the elastic force generated by the deformation of the subsystem and the target building. The first damping force refers to the resistance of the subsystem and the target building to vibration. The first damping force is directly proportional to the magnitude of the velocity during vibration. Specifically, the first elastic force is obtained by multiplying the displacement difference data by the stiffness coefficient, and the first damping force is obtained by multiplying the velocity difference data by the damping coefficient.

[0075] Step S23: Determine the contact force data based on the first elastic force and the first damping force.

[0076] In this embodiment, the contact force data is determined by adding the first elastic force and the first damping force. Specifically, the formula for calculating the contact force data can be expressed as:

[0077]

[0078] Where k is the stiffness coefficient, u(t) is the second displacement data, u1(t) is the first displacement data, and c is the damping coefficient. For the second speed data Let λ(t) be the first velocity data, and λ(t) be the contact force data.

[0079] In this embodiment, the contact force data is calculated using the stiffness coefficient, the damping coefficient, the first displacement data, the first velocity data, and the second velocity data, thereby improving the efficiency of calculating the multi-order frequencies of the target building.

[0080] Furthermore, based on any one of the second and third embodiments, a fourth embodiment of the building multi-order frequency identification method of the present invention is proposed. In this embodiment, reference is made to... Figure 4 The vibration data includes: first displacement and first velocity data at the end of the subsystem closest to the target building, and second displacement and second velocity data at the other end of the subsystem furthest from the target building; first displacement and first velocity data at the contact point between the subsystem and the target building, and second displacement and second velocity data at the furthest end from the contact point; the determination of the structural dynamic characteristic parameters of the target building based on the physical parameters and the vibration data includes:

[0081] Step S201: Determine the second elastic force of the subsystem based on the second displacement data and the stiffness coefficient, and determine the second damping force of the subsystem based on the second velocity data and the damping coefficient;

[0082] The second elastic force refers to the elastic force generated by the deformation of the subsystem, and the second damping force refers to the resistance of the subsystem to vibration. Specifically, the second elastic force is obtained by multiplying the second displacement data by the stiffness coefficient, and the second damping force is obtained by multiplying the second velocity data by the damping coefficient.

[0083] Step S202: Calculate the structural dynamic characteristic parameters based on the contact force data, the second elastic force, and the second damping force.

[0084] Specifically, the dynamic characteristic parameters of the structure are calculated by subtracting the second elastic force and then the second damping force from the contact force data. The calculation formula can be expressed as:

[0085]

[0086] Wherein, Y(t) is the structural dynamic characteristic parameter, λ(t) is the contact force data, k is the stiffness coefficient, c is the damping coefficient, and u is the second displacement data. This is the second speed data.

[0087] In this embodiment, the structural dynamic characteristic parameters are determined by the contact force data, the stiffness coefficient, the damping coefficient, the second displacement data, and the second velocity data, thereby improving the accuracy of calculating the higher-order frequencies of the target building.

[0088] Furthermore, based on any of the above embodiments, a fifth embodiment of the building multi-order frequency identification method of the present invention is proposed. In this embodiment, reference is made to... Figure 5The feature parameters include the similarity matrix corresponding to the multi-order frequencies, the equivalence relationship is the first set of equations between the contact force data, the structural dynamic feature parameters, and the similarity matrix, the target matrix includes the multi-order frequencies, and the determination of the multi-order frequencies based on the equivalence relationship, the contact force data, and the structural dynamic feature parameters includes:

[0089] Step S41: Using the contact force data and the structural dynamic characteristic parameters as known quantities in the first set of equations, solve the first set of equations to obtain the similarity matrix;

[0090] Specifically, the first system of equations here is:

[0091]

[0092] Y(t) = N c X(t)

[0093] in, Let N be the system state vector, and N be the system state vector. M N c Here, N is the time-invariant system matrix, M is the target building's mass matrix, K is the target building's stiffness matrix, and C is the target building's damping matrix. Specifically, for example: in... Figure 6 In the simplified model diagram of the target building shown, m represents the mass of the subsystem. The connection between the subsystem and the top of the target building, from top to bottom, includes a spring and a damper. The mass matrix of the target building is:

[0094]

[0095] The stiffness matrix of the target building is:

[0096]

[0097] The damping matrix of the target building is:

[0098]

[0099] Furthermore, N can be represented as:

[0100]

[0101] N M It can be represented as:

[0102]

[0103] N c It can be represented as:

[0104] Nc =[-k 0 0 -c 0 0]

[0105] N M F(t) under non-typhoon weather conditions can be considered as the systematic error of the first equation system, corresponding to w in the subspace identification method. k Systematic error.

[0106] W H The positional relationship between the subsystem and the target building.

[0107] Specifically, the contact force data and the structural dynamic characteristic parameters are used as known quantities in the first set of equations, and the similarity matrix of N is solved using the Model Identification Algorithm (Multivariable Output-Error States Pace, MOESP).

[0108] Step S42: Calculate at least one eigenvalue of the similarity matrix;

[0109] Since the eigenvalues ​​of the similar matrices of matrix N are equal to the eigenvalues ​​of matrix N, the correspondence can be determined by calculating at least one eigenvalue of the similar matrices.

[0110] Step S43: Determine the multi-order frequency based on the at least one characteristic value.

[0111] The multi-order frequencies are determined based on the sorting of all the eigenvalues ​​in ascending order.

[0112] In this embodiment, the accuracy of the solved multi-order frequencies is ensured by using the contact force data and the structural dynamic characteristic parameters as known quantities in the first set of equations and solving the coherence matrix of the N matrix using the MOESP algorithm.

[0113] Furthermore, the step of establishing the equivalent relationship between the contact force data, the structural dynamic characteristic parameters, and the characteristic parameters corresponding to the multi-order frequencies of the target building based on the kinematic principles of the coupled system includes:

[0114] The kinematic equations of the target building are determined based on the kinematic principles of the coupled system and the first preset degree of freedom data of the target building.

[0115] The kinematic equations and the contact force equations between the subsystem and the target building are transformed into the first set of equations according to the subspace identification method.

[0116] The first preset degree of freedom can be determined by the user. Generally, the user determines the first preset degree of freedom by assuming the order of the unknown system. For example, if the order of the unknown system is 3, the first degree of freedom is determined to be 3. Furthermore, since the degree of freedom of the target building cannot be determined during actual monitoring, the first preset degree of freedom can generally be decided by the user, and the size of the N matrix is ​​determined based on the first preset degree of freedom. Performing a mechanical analysis on the target building yields the kinematic equations of the target building, which include the forces exerted on the target building by the general environmental forces and the forces exerted on the target building by the subsystems.

[0117]

[0118] in, Let the second derivative of U(t) be the acceleration data generated by the vibration of the target building within a preset time period. The first derivative of U(t) represents the velocity data generated by the vibration of the target building within a preset time period, and U(t) represents the displacement data generated by the vibration of the target building within the preset time period.

[0119] The contact force equation can be a formula for calculating contact force data:

[0120]

[0121] Based on the kinematic principles of the coupled system, the kinematic equations and the contact force equations between the subsystem and the target building can be transformed into the first set of equations.

[0122] In this embodiment, by setting a first preset degree of freedom and transforming the kinematic equations into the first set of equations according to the kinematic principle of the coupled system, the efficiency of solving the higher-order matrix can be improved.

[0123] Furthermore, based on the fifth embodiment described above, a sixth embodiment of the building multi-order frequency identification method of the present invention is proposed. In this embodiment, after the step of determining the multi-order frequency based on the at least one feature value, the method further includes:

[0124] The multi-order frequency is determined as a multi-order verification frequency based on the second set of equations, the contact force data, and the structural dynamic characteristic parameters. The second set of equations is determined based on the kinematic equations corresponding to the second preset degree of freedom data and the subspace identification method. The second preset degree of freedom data is greater than the first preset degree of freedom data.

[0125] When the multi-order frequency meets the preset verification conditions, the multi-order frequency is determined to be correct;

[0126] When the multi-order frequency does not meet the preset verification conditions, it is determined that the multi-order frequency is incorrect;

[0127] The multi-order frequency includes multiple first sub-frequencys corresponding to different orders, the multi-order verification frequency includes multiple second sub-frequencys corresponding to different orders, and the preset verification condition is that the first sub-frequency and the second sub-frequency corresponding to the same order are equal.

[0128] In this application, a second preset degree of freedom greater than the first preset degree of freedom is set, and the multi-order frequencies under the second preset degree of freedom are calculated based on the second equation set corresponding to the second preset degree of freedom, the contact force data, and the structural dynamic characteristic parameters, and are thus the multi-order verification frequencies. The correctness of the multi-order frequencies is determined by comparing the multi-order frequencies with the multi-order verification frequencies. Specifically, the multi-order frequencies include multiple first sub-frequencies corresponding to multiple orders, and the multi-order verification frequencies include multiple second sub-frequencies corresponding to multiple orders. The preset verification condition is that the first sub-frequency and the second sub-frequency corresponding to the same order are equal. In other embodiments, a third preset degree of freedom less than the first preset degree of freedom is set, and the multi-order frequencies under the third preset degree of freedom are calculated based on the second equation set corresponding to the third preset degree of freedom, the contact force data, and the structural dynamic characteristic parameters, and are thus the multi-order verification frequencies. The correctness of the multi-order frequencies is determined by comparing the multi-order frequencies with the multi-order verification frequencies.

[0129] In this embodiment, the correctness of the multi-order verification frequency is determined by calculating the multi-order verification frequency and judging whether the multi-order frequency meets the preset verification conditions.

[0130] Furthermore, based on the sixth embodiment described above, a seventh embodiment of the building multi-order frequency identification method of the present invention is proposed. In this embodiment, after the step of determining the multi-order frequency based on the at least one feature value, the method further includes:

[0131] When the stiffness coefficient is less than the preset stiffness coefficient, or when the damping coefficient is less than the preset damping coefficient, the multi-order frequency is determined to be incorrect.

[0132] When the stiffness coefficient is greater than or equal to the preset stiffness coefficient and the damping coefficient is greater than or equal to the preset damping coefficient, the step of determining the multi-order frequency as the multi-order verification frequency based on the second set of equations, the contact force data, and the structural dynamic characteristic parameters is executed.

[0133] In this embodiment, a preset stiffness coefficient and a preset damping coefficient can be set according to the target building to be monitored. When the stiffness coefficient is less than the preset stiffness coefficient, or when the damping coefficient is less than the preset damping coefficient, the multi-order frequency is determined to be incorrect. When the stiffness coefficient is greater than or equal to the preset stiffness coefficient and the damping coefficient is greater than or equal to the preset damping coefficient, it can be determined that the multi-order frequency is a multi-order verification frequency by performing the step of determining the multi-order frequency according to the second set of equations, the contact force data, and the structural dynamic characteristic parameters.

[0134] In this embodiment, increasing the stiffness coefficient of the subsystem can improve the recognition effect of the target building's higher-order frequencies, and increasing the damping coefficient of the subsystem can improve the recognition effect of the target building's manuscript frequencies using the subspace identification method.

[0135] Furthermore, embodiments of the present invention also propose a storage medium storing a building multi-frequency identification program, wherein when the building multi-frequency identification program is executed by a processor, it implements the steps of an embodiment of the building multi-frequency identification method as described in any of the preceding claims.

[0136] It should be noted that, in this document, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or system that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or system. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or system that includes that element.

[0137] The sequence numbers of the above embodiments of the present invention are for descriptive purposes only and do not represent the superiority or inferiority of the embodiments.

[0138] Through the above description of the embodiments, those skilled in the art can clearly understand that the methods of the above embodiments can be implemented by means of software plus necessary general-purpose hardware platforms. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk) as described above, and includes several instructions to cause a terminal device (which may be a mobile phone, computer, server, or network device, etc.) to execute the methods described in the various embodiments of the present invention.

[0139] The above are merely preferred embodiments of the present invention and do not limit the scope of the patent. Any equivalent structural or procedural transformations made based on the description and drawings of the present invention, or direct or indirect applications in other related technical fields, are similarly included within the scope of patent protection of the present invention.

Claims

1. A method for identifying multi-order frequencies in buildings, characterized in that, The building multi-order frequency identification method includes the following steps: The physical parameters of a subsystem installed on the top of a target building are obtained, and vibration data of the subsystem during the vibration of the target building under normal environmental conditions within a preset time period are obtained. The subsystem, from bottom to top, includes a damper, a spring, and a counterweight. The normal environment is non-typhoon weather. The physical parameters include the damping coefficient of the damper, the stiffness coefficient of the spring, and the mass of the counterweight. The vibration data includes the first displacement data and the first velocity data of the subsystem at the end closer to the target building, and the second displacement data and the second velocity data of the subsystem at the other end farther from the target building. The contact force data between the subsystem and the target building and the structural dynamic characteristic parameters of the target building are determined based on the physical parameters and the vibration data. The step of determining the structural dynamic characteristic parameters of the target building includes: determining the second elastic force of the subsystem based on the second displacement data and the stiffness coefficient, and determining the second damping force of the subsystem based on the second velocity data and the damping coefficient; calculating the structural dynamic characteristic parameters based on the contact force data, the second elastic force, and the second damping force. Based on the kinematics principle of coupled systems, an equivalent relationship is established between the contact force data, the structural dynamic characteristic parameters, and the characteristic parameters corresponding to the multi-order frequencies of the target building. The characteristic parameters include the similarity matrix corresponding to the multi-order frequencies. The equivalent relationship is the first set of equations between the contact force data, the structural dynamic characteristic parameters, and the similarity matrix. Determining the multiple frequencies based on the equivalence relation, the contact force data, and the structural dynamic characteristic parameters includes: using the contact force data and the structural dynamic characteristic parameters as known quantities in the first set of equations, solving the first set of equations to obtain the similarity matrix; calculating at least one eigenvalue of the similarity matrix; and determining the multiple frequencies based on the at least one eigenvalue. The first set of equations is as follows: ,in, Let U(t) be the system state vector, and U(t) be the displacement data of the target building caused by vibration within a preset time period. Let the first derivative of U(t) be the velocity data generated by the vibration of the target building within a preset time period, and N, N M N c Here, N is the time-invariant system matrix, N represents the frequency information of the target building, and M is the mass matrix of the target building. Let λ(t) represent the positional relationship between the subsystem and the target building, and let λ(t) represent the contact force data. The systematic error of the first set of equations under non-typhoon weather conditions is given.

2. The building multi-order frequency identification method as described in claim 1, characterized in that, The step of determining the contact force data between the subsystem and the target building based on the physical parameters and the vibration data includes: Determine the displacement difference between the first displacement data and the second displacement data, and determine the velocity difference between the first velocity data and the second velocity data; The first elastic force between the subsystem and the target building is calculated based on the displacement difference data and the stiffness coefficient, and the first damping force between the subsystem and the target building is calculated based on the velocity difference data and the damping coefficient. The contact force data is determined based on the first elastic force and the first damping force.

3. The building multi-order frequency identification method as described in claim 1, characterized in that, The step of establishing the equivalent relationship between the contact force data, the structural dynamic characteristic parameters, and the characteristic parameters corresponding to the multi-order frequencies of the target building based on the kinematic principles of the coupled system includes: The kinematic equations of the target building are determined based on the kinematic principles of the coupled system and the first preset degree-of-freedom data of the target building. The kinematic equations and the contact force equations between the subsystem and the target building are transformed into the first set of equations using the subspace identification method.

4. The building multi-frequency identification method as described in claim 3, characterized in that, After the step of determining the multi-order frequency based on the at least one feature value, the method further includes: The multi-order frequency is determined as a multi-order verification frequency based on the second set of equations, the contact force data, and the structural dynamic characteristic parameters. The second set of equations is determined based on the kinematic equations corresponding to the second preset degree of freedom data and the subspace identification method. The second preset degree of freedom data is greater than the first preset degree of freedom data. When the multi-order frequency meets the preset verification conditions, the multi-order frequency is determined to be correct; When the multi-order frequency does not meet the preset verification conditions, it is determined that the multi-order frequency is incorrect; The multi-order frequency includes multiple first sub-frequencys corresponding to different orders, the multi-order verification frequency includes multiple second sub-frequencys corresponding to different orders, and the preset verification condition is that the first sub-frequency and the second sub-frequency corresponding to the same order are equal.

5. The building multi-order frequency identification method as described in claim 4, characterized in that, After the step of determining the multi-order frequency based on the at least one feature value, the method further includes: When the stiffness coefficient is less than the preset stiffness coefficient, or when the damping coefficient is less than the preset damping coefficient, the multi-order frequency is determined to be incorrect. When the stiffness coefficient is greater than or equal to the preset stiffness coefficient and the damping coefficient is greater than or equal to the preset damping coefficient, the step of determining the multi-order frequency as the multi-order verification frequency based on the second set of equations, the contact force data, and the structural dynamic characteristic parameters is executed.

6. A subsystem, characterized in that, The subsystem is installed on the top of the target building and includes: a memory, a processor, and a building multi-frequency identification program stored in the memory and executable on the processor. The building multi-frequency identification program is configured to implement the steps of the building multi-frequency identification method as described in any one of claims 1 to 5.

7. A storage medium, characterized in that, The storage medium stores a building multi-frequency identification program, which, when executed by a processor, implements the steps of the building multi-frequency identification method as described in any one of claims 1 to 5.

Images