Three-dimensional displacement high-precision monitoring method and system for shock insulation support

Abstract

The invention relates to the technical field of data processing, in particular to a three-dimensional displacement high-precision monitoring method and system for a seismic isolation support, and the method comprises the steps: obtaining acceleration data of the seismic isolation support in different axial directions and a spectrogram of each axial direction within a preset time period; according to the amplitude distribution condition of each frequency in each axial spectrogram, screening out a seismic main frequency, and according to the amplitude proportion condition of the seismic main frequency in each axial spectrogram, obtaining a relative energy value in each axial direction; correcting the acceleration data of each axial direction according to the correlation between the acceleration data of each axial direction and the acceleration data of other axial directions in combination with the distribution condition and the relative energy value of the acceleration data of each axial direction to obtain corrected data of each axial direction; and determining the three-dimensional displacement data of the shock insulation support based on the integral result of each axial correction data. According to the invention, other interference components can be effectively inhibited, and the monitoring precision of the seismic isolation support displacement data is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of data processing, in particular to a three-dimensional displacement high-precision monitoring method and system for an isolation bearing. BACKGROUND

[0002] An isolation bearing is a device installed between the foundation and the upper structure of a building, designed to absorb and reduce the energy transmitted by seismic waves, thereby reducing the impact of earthquakes on the building. By accurately monitoring the three-dimensional displacement of the isolation bearing under seismic waves, the deformation of the isolation bearing can be evaluated in real time, and its working condition can be judged, and potential damage or failure risks can be found in time. In addition, monitoring the three-dimensional displacement of the isolation bearing under seismic waves also helps to optimize the isolation design and verify the isolation effect, providing reliable basis for subsequent engineering improvement.

[0003] Generally, the three-dimensional displacement high-precision monitoring of the isolation bearing during the earthquake is carried out by acceleration integration method. Acceleration sensors are often used to monitor the acceleration changes of the isolation bearing under the action of earthquakes in real time. These sensors can accurately capture the acceleration data of the bearing in X, Y and Z axial directions. By performing second-order integration on the acceleration data, the displacement of the bearing in the three axial directions can be calculated, thereby providing a scientific basis for seismic design optimization and improving the seismic capacity and safety of the building.

[0004] However, when using acceleration integration method to determine the three-dimensional displacement of the isolation bearing under earthquake, it may also be disturbed by other non-seismic signals, which may come from equipment noise, temperature change or other environmental factors. Although these signals are not related to seismic waves, they will also affect the detection results of acceleration sensors, thereby introducing errors and making the three-dimensional displacement detection results of the isolation bearing less accurate. SUMMARY

[0005] In order to solve the technical problem that the three-dimensional displacement of the isolation bearing under earthquake determined by the existing acceleration integration method is affected by other factors, resulting in poor accuracy of the three-dimensional displacement detection results of the isolation bearing, the purpose of the present application is to provide a three-dimensional displacement high-precision monitoring method and system for an isolation bearing, and the technical solution adopted is as follows: In a first aspect, the present application provides a three-dimensional displacement high-precision monitoring method for an isolation bearing, comprising: obtaining acceleration data of the isolation bearing in different axial directions and frequency spectrum of each axial direction within a preset time period; screening out the main frequency of the earthquake according to the amplitude distribution of each frequency in the frequency spectrum of each axial direction, and obtaining the relative energy value of each axial direction according to the amplitude proportion of the main frequency of the earthquake in the frequency spectrum of each axial direction; According to the correlation between the acceleration data of each axis and the acceleration data of other axes, the distribution of the acceleration data of each axis and the relative energy value, the acceleration data of each axis is corrected to obtain the correction data of each axis; Based on the integral result of the correction data of each axis, the three-dimensional displacement data of the isolation bearing is determined.

[0006] Preferably, the main frequency of the earthquake is screened according to the amplitude value distribution of each frequency in the frequency spectrum of each axis, specifically including: According to the amplitude value distribution of each frequency and the adjacent frequency in the frequency spectrum of each axis, the distribution characteristic value of each frequency in the frequency spectrum of each axis is obtained in combination with the value of each frequency; According to the distribution characteristic value and the amplitude value distribution of each frequency in all axes, the main frequency index of the earthquake of each frequency is obtained. The frequency corresponding to the maximum value of the main frequency index of the earthquake is taken as the main frequency of the earthquake.

[0007] Preferably, the distribution characteristic value of each frequency in the frequency spectrum of each axis is obtained according to the amplitude value distribution of each frequency and the adjacent frequency in the frequency spectrum of each axis in combination with the value of each frequency, specifically including: For the frequency spectrum of any one axis, any one frequency is taken as a selected frequency; All frequencies are arranged in a preset order according to the value, and the reference frequency adjacent to the selected frequency is obtained in the arrangement order of the frequency; Based on the amplitude value of the selected frequency and the amplitude value of the reference frequency, a amplitude comparison coefficient is determined; based on the value of the selected frequency, a frequency characteristic coefficient is determined; The product of the amplitude comparison coefficient and the frequency characteristic coefficient is taken as the distribution characteristic value of the selected frequency in the frequency spectrum of the arbitrary axis.

[0008] Preferably, the main frequency index of the earthquake of each frequency is obtained according to the distribution characteristic value and the amplitude value distribution of each frequency in all axes, specifically including: The product of the amplitude value and the distribution characteristic value of the selected frequency in the frequency spectrum of each axis is taken as the main frequency index of the selected frequency in each axis; the cumulative sum of the main frequency index of the selected frequency in all axes is taken as the main frequency index of the selected frequency of the earthquake.

[0009] Preferably, the acceleration data of each axis is corrected according to the correlation between the acceleration data of each axis and the acceleration data of other axes, the distribution of the acceleration data of each axis and the relative energy value, to obtain the correction data of each axis, specifically including: According to the comprehensive result of the correlation between the acceleration data of each axial direction and the acceleration data of each other axial direction, a data usability factor of each axial direction is obtained; According to the data usability factor of each axial direction and the balance of the acceleration data of the same axial direction, a relative distribution of the data of each axial direction is analyzed, and a seismic wave effective factor of each axial direction is obtained in combination with the difference between the relative distribution of the data and the relative energy value; The acceleration data of each axial direction is corrected by using the seismic wave effective factor of each axial direction, and corrected data of each axial direction is obtained.

[0010] Preferably, the data usability factor of each axial direction is obtained according to the comprehensive result of the correlation between the acceleration data of each axial direction and the acceleration data of each other axial direction, and specifically includes: Any one axial direction is taken as a selected axial direction, and each axial direction other than the selected axial direction is taken as a reference axial direction; A correlation coefficient between the acceleration data corresponding to the selected axial direction and each reference axial direction is calculated, and the average of the correlation coefficients between the selected axial direction and all reference axial directions is taken as the data usability factor of the selected axial direction.

[0011] Preferably, the relative distribution of the data of each axial direction is analyzed according to the data usability factor of each axial direction and the balance of the acceleration data of the same axial direction, and the seismic wave effective factor of each axial direction is obtained in combination with the difference between the relative distribution of the data and the relative energy value, and specifically includes: The product of the average of all acceleration data of the selected axial direction and the data usability factor is taken as a response characteristic index of the selected axial direction, and the ratio between the response characteristic index of the selected axial direction and the cumulative sum of the response characteristic indexes of all axial directions is taken as a relative contribution degree of the selected axial direction. The difference between the relative contribution degree of the selected axial direction and the relative energy value is negatively correlated and normalized to obtain the seismic wave effective factor of the selected axial direction.

[0012] Preferably, the acceleration data of each axial direction is corrected by using the seismic wave effective factor of each axial direction to obtain the corrected data of each axial direction, and specifically includes: The product of the seismic wave effective factor of the selected axial direction and each acceleration data of the selected axial direction is taken as the corresponding correction data of each acceleration data of the selected axial direction.

[0013] Preferably, the relative energy value of each axial direction is obtained according to the amplitude proportion of the seismic main frequency of each axial direction, and specifically includes: The ratio between the amplitude of the seismic main frequency of each axial direction and the cumulative sum of the amplitudes of the seismic main frequencies of all axial directions is taken as the relative energy value of each axial direction.

[0014] In a second aspect, the present application provides a high-precision three-dimensional displacement monitoring system for an isolation bearing, comprising a memory, a processor, and a computer program stored in the memory and running on the processor, which, when executed by the processor, implements the steps of a high-precision three-dimensional displacement monitoring method for an isolation bearing.

[0015] The embodiments of the present application have at least the following beneficial effects: The present application first acquires acceleration data of each axis and a corresponding frequency spectrum, and extracts frequency components of each axis. Then, the first aspect screens the main frequency performance of each axis that best matches the seismic wave distribution according to the characteristics of the frequency distribution, and determines the main frequency of the seismic wave. The second aspect can reflect the energy proportion of the seismic wave propagation in each axis according to the frequency amplitude proportion of each axis at the seismic main frequency. Further, not only the correlation of acceleration data between different axes is considered, but also the consistency between the distribution of acceleration data of each axis and the relative energy value, so that the process of correcting the acceleration data of each axis fully considers the characteristics of the seismic wave, and the higher the consistency between the acceleration data and the characteristics of the seismic wave, the higher the data retention. Finally, the corrected acceleration data will be used for the acceleration integration method to obtain the high-precision displacement value of each axis of the isolation bearing under the action of the earthquake. The correction process of the acceleration data of each axis can strengthen the components of the acceleration data of each axis related to the seismic wave, while suppressing noise or interference signals unrelated to the seismic wave, so as to obtain more accurate acceleration data under the influence of the seismic wave, and improve the precision of the displacement data obtained by the final integration method. BRIEF DESCRIPTION OF DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, and the advantages thereof, a brief introduction will be given to the drawings needed in the embodiments or prior art description. Obviously, the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can be obtained without creative labor based on these drawings.

[0017] Figure 1 is a step flow chart of a high-precision three-dimensional displacement monitoring method for an isolation bearing provided by the present application; Figure 2 is an acceleration change curve of the X-axis in a period of time provided by the present application; Figure 3 is an acceleration change curve of the Y-axis in a period of time provided by the present application; Figure 4 is an acceleration change curve of the Z-axis in a period of time provided by the present application; Figure 5is a curve diagram of acceleration data of the X axis in a preset time period provided by the application; Figure 6 is a curve diagram of acceleration data of the Y axis in a preset time period provided by the application; Figure 7 is a curve diagram of acceleration data of the Z axis in a preset time period provided by the application; Figure 8 is a corresponding frequency spectrum diagram of the X axis in a preset time period provided by the application; Figure 9 is a corresponding frequency spectrum diagram of the Y axis in a preset time period provided by the application; Figure 10 is a corresponding frequency spectrum diagram of the Z axis in a preset time period provided by the application; Figure 11 is a step flow chart of the method for obtaining correction data of each axial direction provided by the application. DETAILED DESCRIPTION

[0018] In order to further illustrate the technical means and effects taken by the application to achieve the predetermined application purpose, the following describes the specific implementation, structure, features and effects of the three-dimensional displacement high-precision monitoring method and system for the seismic isolation support according to the application in combination with the preferred embodiments and the drawings. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. In addition, the specific features, structures or characteristics in one or more embodiments can be combined in any suitable form.

[0019] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the application belongs.

[0020] The following specifically describes the specific scheme of the three-dimensional displacement high-precision monitoring method and system for the seismic isolation support provided by the application in combination with the drawings.

[0021] Please refer to Figure 1 which shows a step flow chart of the three-dimensional displacement high-precision monitoring method for the seismic isolation support provided by one embodiment of the application. The method comprises the following steps: Step S100, obtaining acceleration data of the seismic isolation support in different axial directions in a preset time period and a frequency spectrum diagram of each axial direction.

[0022] Specifically, the embodiment is directed to the acceleration data performance of three different axial isolation bearings when generating vibrations, where the three axes are X-axis, Y-axis and Z-axis respectively, and the purpose is to capture the lateral movement of the isolation bearing in the horizontal plane, the longitudinal movement of the isolation bearing in the horizontal plane and the vertical movement of the isolation bearing in the vertical direction respectively. The three form a three-dimensional orthogonal coordinate system, and through the installation of sensors, the three-dimensional acceleration changes of the bearing under the action of earthquakes can be captured in real time, laying a foundation for subsequent high-precision displacement calculation.

[0023] More specifically, the embodiment uses MEMS accelerometers to collect acceleration data. As a specific example, the installation position of the MEMS accelerometer can be: for the X-axis and Y-axis MEMS accelerometers: they can be installed on both sides of the isolation bearing (or the upper and lower parts of the bearing), to ensure that the acceleration in the horizontal plane (X-axis and Y-axis) can be monitored respectively. For the Z-axis MEMS accelerometer: it should be installed in the vertical direction of the isolation bearing (usually at the lower part or bottom surface of the bearing), for measuring the acceleration in the vertical direction. This is a known technology, and only a brief introduction is given here. In other embodiments, the implementer can set it according to the specific implementation scene.

[0024] In the embodiment, when the seismic wave acts on the bearing, the sensor will immediately respond and start capturing the acceleration changes. The system records the acceleration data measured by the sensor in real time, and generates acceleration change curves of three axes. As shown in Figure 2 , Figure 3 and Figure 4 , the acceleration changes of X-axis, Y-axis and Z-axis in a period of time are collected respectively. Further, the embodiment divides the earthquake process into several preset time periods, which helps to more accurately analyze the propagation characteristics of the seismic wave in different time periods. The time length of different preset time periods is equal.

[0025] As a specific example, as shown in Figure 5 , Figure 6 and Figure 7 , the curve diagrams of the acceleration data of X-axis, Y-axis and Z-axis in a preset time period are shown, where the time length of the preset time period is 5 seconds.

[0026] Finally, the acceleration data of the isolation bearing in each axis in the preset time period is subjected to Fourier transform to obtain the frequency spectrum corresponding to the acceleration data of each axis, as shown in Figure 8 , Figure 9 and Figure 10 . Among them, the method of Fourier transform is a known technology, and will not be introduced in detail here.

[0027] Step S200: Select the main seismic frequency based on the amplitude distribution of each frequency in the spectrum diagram of each axis, and obtain the relative energy value of each axis based on the amplitude proportion of the main seismic frequency in the spectrum diagram of each axis.

[0028] Firstly, during an earthquake, the dominant frequency components of seismic waves determine the characteristics of the waves, especially their impact on structures. When performing Fourier transforms on the time-series accelerations along the X, Y, and Z axes, the resulting spectrograms reflect the vibration intensity of different frequency components. Specifically, seismic waves tend to exhibit a lower frequency range, with low-frequency components often representing the dominant seismic frequency. Lower frequencies and sharp peaks typically correspond to areas of concentrated energy transmission. Therefore, components with lower frequencies and sharp peaks in the spectrograms generally reflect the dominant frequency of the seismic wave well. Based on this characteristic, the most consistent dominant frequency characteristics, i.e., the dominant seismic frequency, are selected by analyzing the amplitude distribution of each frequency component in the spectrograms for each axis.

[0029] As a specific example, the method for selecting the main seismic frequency based on the amplitude distribution of each frequency in the spectrum of each axis can be implemented by steps S201 to S203.

[0030] Step S201: Based on the amplitude distribution of each frequency and its adjacent frequencies in the spectrum diagram of each axis, and combined with the value of each frequency, obtain the distribution characteristic value of each frequency in the spectrum diagram of each axis.

[0031] Specifically, the feature analysis process is exactly the same for the X, Y, and Z axes. This embodiment uses the spectrum of any one axis as an example for illustration. For example, taking the spectrum of the X-axis as an example, any frequency in the spectrum corresponding to the X-axis is selected as the chosen frequency.

[0032] The first step is to arrange all frequencies according to the preset order of their values, and then obtain the reference frequency adjacent to the selected frequency from the frequency arrangement.

[0033] like Figure 8 As shown, the peaks in the downward-facing spectrum along the X-axis are sharp at lower frequencies. Based on this characteristic, the amplitude distribution between each frequency and its two adjacent frequencies is analyzed. More specifically, in this embodiment, the preset order is ascending, meaning all frequencies in the downward-facing spectrum along the X-axis are arranged in ascending order of frequency value. In this arrangement, the frequency adjacent to the left of the selected frequency and the frequency adjacent to the right of the selected frequency are both recorded as adjacent reference frequencies.

[0034] It should be noted that, for the left and right adjacent two frequencies cannot be obtained at the same time, this feature analysis is not performed.

[0035] In the second step, based on the amplitude of the selected frequency and the amplitude of the reference frequency, the amplitude comparison coefficient is determined; based on the value of the selected frequency corresponding to the negative correlation coefficient, the frequency feature coefficient is determined; and the product of the amplitude comparison coefficient and the frequency feature coefficient is taken as the distribution characteristic value of the selected frequency in the frequency spectrum of the arbitrary axis.

[0036] As a specific example, taking the i-th frequency in the X-axis frequency spectrum as the selected frequency, the distribution characteristic value of the selected frequency in the X-axis frequency spectrum can be expressed by the formula: wherein, represents the distribution characteristic value of the selected frequency in the X-axis frequency spectrum, i represents the i-th frequency in the frequency spectrum, represents the amplitude corresponding to the selected frequency in the X-axis frequency spectrum, that is, the amplitude corresponding to the i-th frequency in the X-axis frequency spectrum; represents the amplitude corresponding to a reference frequency in the X-axis frequency spectrum, that is, the amplitude corresponding to the i-1-th frequency in the X-axis frequency spectrum; represents the amplitude corresponding to another reference frequency in the X-axis frequency spectrum, that is, the amplitude corresponding to the i+1-th frequency in the X-axis frequency spectrum; represents the frequency value corresponding to the selected frequency in the X-axis frequency spectrum, for example, 5 Hz. represents the exponential function with the natural constant e as the base.

[0037] is the amplitude comparison coefficient, which reflects the amplitude comparison result between the adjacent frequencies of the selected frequency in the frequency spectrum. The greater the value, the greater the amplitude of the selected frequency component, and the smaller the amplitude of the surrounding frequencies, indicating that the degree of existence of the sharp peak at the position of the selected frequency in the frequency spectrum is greater. is the frequency feature coefficient, which is processed by the negative correlation of the frequency value, so that the low frequency feature can be screened out.

[0038] When the amplitude comparison coefficient and the frequency feature coefficient have large values, it indicates that the amplitude distribution characteristic of the selected frequency conforms to the degree of the seismic wave feature, and further indicates that the selected frequency belongs to the main frequency of the seismic wave.

[0039] It should be understood that the amplitude corresponding to the frequency is the value of the vertical coordinate in the frequency spectrum, and the value of the frequency is the value of the horizontal coordinate in the frequency spectrum.

[0040] Step S202, according to the distribution characteristic value and the amplitude distribution of each frequency in all axial corresponding, the seismic main frequency index of each frequency is obtained.

[0041] In the propagation process of the seismic wave, the main frequency of the seismic wave usually has an impact on the vibration of multiple axes. If the spectral amplitude of the three axes is high at a certain frequency and meets the seismic main frequency, it means that the frequency component has significant energy transmission in three directions. This indicates that the frequency component is highly consistent with the main frequency of the seismic wave and may be the main frequency component of the seismic wave.

[0042] Specifically, the seismic wave usually has a certain main frequency in the propagation process and can cause strong response in all directions. If the frequency component is observed to have a high amplitude in three axes on the basis of meeting the seismic main frequency, it means that it has a great contribution to the propagation of the seismic wave in all directions, so it can be considered that the frequency as a whole meets the seismic main frequency.

[0043] Based on this feature, the product of the amplitude of the selected frequency in the spectrum of each axis and the distribution characteristic value is taken as the main frequency factor of each axis at the selected frequency; and the cumulative sum of the main frequency factors of all axes at the selected frequency is taken as the seismic main frequency index of the selected frequency.

[0044] Step S203, the frequency corresponding to the maximum value of the seismic main frequency index is taken as the seismic main frequency.

[0045] The seismic main frequency index of the selected frequency represents the characteristic performance degree of the selected frequency in meeting the seismic main frequency in the three axes of X, Y, and Z axes as a whole. The greater the value of the seismic main frequency index, the greater the characteristic performance degree of the corresponding frequency in meeting the seismic main frequency, and the greater the possibility that the corresponding frequency belongs to the main frequency of the seismic wave.

[0046] In the second aspect, the seismic main frequency represents the main energy component in the propagation of the seismic wave and is usually the core frequency of the seismic wave. At the seismic main frequency, the vibration components of the X, Y, and Z axes reflect the energy distribution of the seismic wave in all directions. Based on this feature, by analyzing the amplitude proportion of the seismic main frequency of each axis, the propagation of the seismic wave in each axis is evaluated, and the relative energy value of the seismic wave in each axis is quantitatively obtained.

[0047] Specifically, the ratio between the amplitude of the seismic main frequency of each axis and the cumulative sum of the amplitudes of the seismic main frequencies of all axes is taken as the relative energy value of each axis.

[0048] By analyzing the amplitude proportion of the frequency component corresponding to the main frequency of the earthquake in different axial directions, the propagation characteristics of the seismic wave in each direction can be revealed. If the amplitude proportion of the frequency component in a certain axial direction is large, it means that the propagation energy of the seismic wave in this direction is strong, and vice versa. That is, the relative energy value corresponding to each axial direction represents the energy distribution ratio of the seismic wave in each axial direction, which is a quantitative index of the inherent characteristics of the seismic wave.

[0049] In step S300, the acceleration data of each axial direction is corrected according to the correlation between the acceleration data of each axial direction and other axial directions, combined with the distribution of the acceleration data of each axial direction and the relative energy value, to obtain the corrected data of each axial direction.

[0050] To determine whether the acceleration of the three axial directions of the X-axis, Y-axis and Z-axis conforms to the nature of the seismic wave, the consistency between the acceleration and the propagation energy of the seismic wave needs to be considered on the basis of the comparative availability of the acceleration. The greater the consistency, the smaller the degree of interference of the seismic wave, and the smaller the degree of correction of the original data. The smaller the consistency, the greater the degree of influence of the effective seismic wave component of the acceleration data by other interference signals, and the greater the degree of correction of the original data.

[0051] Therefore, the first aspect analyzes the data correlation between the acceleration data in different axial directions, the second aspect analyzes the relative effectiveness of the acceleration data in each axial direction, and the third aspect analyzes the consistency between the relative effectiveness of the acceleration data in each axial direction and the relative energy distribution. The results of the three aspects of feature analysis are combined to filter the effective seismic signal and suppress non-seismic interference in the acceleration data of each axial direction to realize the correction operation of the acceleration data.

[0052] As a specific example, as shown in Figure 11 The method for obtaining the corrected data of each axial direction can be realized by steps S301 to S303.

[0053] In step S301, the data availability factor of each axial direction is obtained according to the comprehensive results of the correlation between the acceleration data of each axial direction and other axial directions.

[0054] During the propagation of seismic waves, there is a certain directionality, that is, the seismic waves will affect the vibration of different axial directions when propagating along a specific direction. In order to determine the comparative availability of the three axial acceleration data of the X-axis, Y-axis and Z-axis, the correlation between the acceleration of different axial directions is analyzed. If the Pearson correlation coefficient between the acceleration data of the X-axis, Y-axis and Z-axis tends to 1 within a certain period of time, it means that the acceleration data of the three axial directions changes almost completely synchronously, which reflects their response to the same seismic wave. Therefore, the closer the Pearson correlation coefficient between each axial direction and other axial directions is to 1, the greater the comparative availability of the acceleration change of the X-axis, Y-axis and Z-axis.

[0055] Based on this feature, when the comprehensive result of the Pearson correlation coefficient between the acceleration data of each axial direction and other axial directions is close to 1, it indicates that the feature distribution of the acceleration data of the three axial directions has high consistency, and further indicates that their response under the action of seismic waves is reliable.

[0056] Specifically, in the first step, any one axial direction is taken as a selected axial direction, and each axial direction other than the selected axial direction is taken as a reference axial direction. For example, the axial direction corresponding to the X-axis is taken as the selected axial direction, and the axial directions corresponding to the Y-axis and the Z-axis are both reference axial directions.

[0057] In the second step, the correlation coefficient between the acceleration data of the selected axial direction and each reference axial direction is calculated, and the average of the correlation coefficients between the selected axial direction and all reference axial directions is taken as the data availability factor of the selected axial direction.

[0058] As a specific example, the data availability factor of the selected axial direction can be represented by the formula: wherein represents the data availability factor of the X-axis corresponding axial direction, that is, the data availability factor of the selected axial direction, represents the set of acceleration data of the X-axis corresponding axial direction, represents the set of acceleration data of the Y-axis corresponding axial direction, represents the set of acceleration data of the Z-axis corresponding axial direction, represents the Pearson correlation coefficient between the acceleration data of the X-axis corresponding axial direction and the Y-axis corresponding axial direction, represents the Pearson correlation coefficient between the acceleration data of the X-axis corresponding axial direction and the Z-axis corresponding axial direction.

[0059] The above calculation process reflects the comprehensive result of the correlation coefficient between the acceleration data of the X-axis and the Y-axis, and the correlation coefficient between the acceleration data of the X-axis and the Z-axis, and balancedly represents the availability and reliability of the acceleration data of the X-axis corresponding axial direction.

[0060] Step S302, according to the data available factor of each axial and the balance of acceleration data of the same axial, analyze the relative distribution of data of each axial, combine the difference between the relative distribution of data and the relative energy value, obtain the seismic wave effective factor of each axial.

[0061] Considering that the propagation direction of the seismic wave is unique, the balance of acceleration data of each axial should be consistent with the axial component of the propagation direction of the seismic wave, if the effective distribution of acceleration data of a certain axial is close to the energy distribution in the axial, it indicates that the relationship between the acceleration data in the axial and the propagation direction of the seismic wave is consistent, that is, the characteristic distribution of the acceleration data in the axial conforms to the characteristics of the seismic wave.

[0062] Specifically, the embodiment takes any one axial as an example for illustration, specifically taking the X-axis corresponding axial as an example for illustration, that is, the embodiment takes the X-axis corresponding axial as the selected axial. First, the product of the mean value of all acceleration data of the selected axial and the data available factor is taken as the response characteristic index of the selected axial; the ratio between the response characteristic index of the selected axial and the cumulative sum of the response characteristic indexes of all axials is taken as the relative contribution degree of the selected axial.

[0063] As a specific example, the relative contribution degree of the selected axial can be expressed by the formula: Among them, The relative contribution degree of the selected axial is represented by X, the X-axis corresponding axial, The mean value of all acceleration data of the X-axis corresponding axial is represented by X, The data available factor of the X-axis corresponding axial is represented by X, The mean value of all acceleration data of the W-axis corresponding axial is represented by W, The data available factor of the W-axis corresponding axial is represented by W, and the value of W includes X, Y and Z.

[0064] Reflects the intensity of X-axis acceleration data, Reflects the consideration and availability of X-axis acceleration data, The response characteristic index, the greater the value, the stronger the availability of X-axis acceleration data, the higher the synchronization between X-axis acceleration data and acceleration data of other axials, and the more likely it is caused by the same seismic source.

[0065] Reflects the cumulative result of the total effective response of the seismic isolation bearing in the current preset time period, the ratio calculation result reflects the proportion of the effective response in the X-axis corresponding axial in the overall effective response, and reflects the relative effective contribution between the acceleration data of the X-axis and the three-dimensional overall data distribution.

[0066] Secondly, the difference between the relative contribution degree and the relative energy value of the selected axial direction is negatively correlated and normalized to obtain the seismic wave effective factor of the selected axial direction.

[0067] The relative contribution degree of the selected axial direction reflects the relative effective response proportion of the acceleration data of the selected axial direction and all axial directions as a whole, and the relative energy value of the selected axial direction reflects the energy distribution proportion of the seismic wave in the selected axial direction, representing the directional characteristic performance of the seismic wave in the X-axis direction.

[0068] When the relative contribution degree and the relative energy value of the selected axial direction are closer, it means that the actual acceleration distribution in the selected axial direction and the actual energy distribution of the seismic wave in the selected axial direction are closer or more similar, and further, it means that the actual acceleration distribution in the selected axial direction conforms to the seismic wave characteristics to a greater extent, and further, it means that the effective component of the acceleration data in the selected axial direction is more and the interference degree is smaller.

[0069] When the difference between the relative contribution degree and the relative energy value of the selected axial direction is greater, it means that the difference between the actual acceleration distribution in the selected axial direction and the actual energy distribution of the seismic wave in the selected axial direction is greater, and further, it means that the effective degree of the acceleration data in the selected axial direction is smaller and the interference degree is greater.

[0070] As a specific example, the calculation process of the seismic wave effective factor of the selected axial direction can be represented by the formula: wherein, represents the seismic wave effective factor of the selected axial direction, X represents the axial direction corresponding to the X-axis, represents the relative contribution degree of the selected axial direction, represents the relative energy value of the selected axial direction, is a linear normalization function.

[0071] At this point, the seismic wave effective factor of each axial direction can be obtained according to the same method, representing the size of the characteristic distribution of the acceleration data in each axial direction conforming to the seismic wave characteristics.

[0072] Step S303, using the seismic wave effective factor of each axial direction, the acceleration data of each axial direction is corrected to obtain the corrected data of each axial direction.

[0073] In the seismic wave analysis process, in order to remove the non-seismic wave interference signal, it is necessary to remove the components irrelevant to the seismic wave from the acceleration data. The more the characteristic distribution of the acceleration data under each axis conforms to the characteristics of the seismic wave, the stronger the correlation between the acceleration data under the axis and the main energy propagation of the seismic wave, and the higher the degree of actual acceleration distribution under the axis that can be retained. The greater the difference between the characteristic distribution of the acceleration data under each axis and the characteristics of the seismic wave, the weaker the correlation between the acceleration data under the axis and the main energy propagation of the seismic wave, and the lower the degree of actual acceleration distribution under the axis that can be retained.

[0074] Specifically, the product of the selected axial seismic wave effective factor and each acceleration data under the selected axis is taken as the correction data corresponding to each acceleration data under the selected axis. It should be understood that for each acceleration data under each axis within the preset time period, the seismic wave effective factor can be used for correction to obtain the corresponding corrected acceleration data, that is, the correction data.

[0075] The seismic wave effective factor reflects the effective strength of the characteristic distribution of the acceleration data under each axis conforming to the characteristics of the seismic wave. By using the product form, the components of the acceleration data under each axis related to the seismic wave can be strengthened, and the noise or interference signal irrelevant to the seismic wave can be suppressed. The correction operation on the actual acceleration data is equivalent to retaining the effective seismic wave components in the acceleration and removing other interference signals, so as to obtain more accurate acceleration data under the influence of the seismic wave.

[0076] Step S400, based on the integral results of the correction data under each axis, the three-dimensional displacement data of the seismic isolation support is determined.

[0077] Specifically, within the preset time period, the secondary integral of all correction data under each axis can obtain the displacement data under the corresponding axis. In this embodiment, for the three different axes of the X axis, the Y axis and the Z axis, the high-precision displacement data of the seismic isolation support within the preset time period can be obtained, which can accurately reflect the high-precision displacement influence of the seismic wave on the three axes of the seismic isolation support. If there are multiple preset time periods in the displacement monitoring process of the seismic isolation support, the displacement data of all preset time periods can also be accumulated to obtain the total displacement in other embodiments.

[0078] It should be noted that the method of secondary integral of acceleration is a known technology, that is, acceleration is the second derivative of displacement, which means that displacement can be obtained by twice integral of acceleration, and will not be described in detail here.

[0079] In summary, compared with the influence caused by non-seismic waves in the prior art, the influence is gradually accumulated by the acceleration integration method, so that the three-dimensional displacement monitoring of the seismic isolation support under the earthquake is inaccurate. The embodiment of the application accurately obtains the three-dimensional high-precision displacement of the seismic isolation support in the X-axis, Y-axis and Z-axis under the earthquake by analyzing the acceleration change spectrum of the seismic isolation support in the X-axis, Y-axis and Z-axis. Specifically, first, the frequency spectrum of the X-axis, Y-axis and Z-axis in each time period is analyzed by Fourier transform, and each frequency component is extracted, so as to determine the main frequency of the seismic wave. Then, the main propagation energy of the seismic wave is measured by the amplitude intensity of the X-axis, Y-axis and Z-axis under the main frequency. Subsequently, according to whether the characteristic distribution of the acceleration of the X-axis, Y-axis and Z-axis is consistent with the characteristics of the seismic wave, the acceleration data of each axis is corrected to ensure that it is consistent with the characteristics of the seismic wave and remove the non-seismic wave interference signal. Finally, the corrected acceleration data will be used for the acceleration integration method to calculate the high-precision displacement value of the seismic isolation support in the X-axis, Y-axis and Z-axis under the action of the earthquake.

[0080] The embodiment of the application also provides a three-dimensional displacement high-precision monitoring system for a seismic isolation support, which comprises a memory, a processor and a computer program stored in the memory and running on the processor, and the computer program realizes the steps of a three-dimensional displacement high-precision monitoring method for a seismic isolation support when executed by the processor. Since the three-dimensional displacement high-precision monitoring method for a seismic isolation support is described in detail, it will not be described in detail here.

[0081] The above-described embodiments are only used to illustrate the technical solutions of the present application, and not to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that they can modify the technical solutions recorded in the foregoing embodiments, or make equivalent replacements to some technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present application, and should be included in the protection scope of the present application.

Claims

1. A method for high-precision monitoring of three-dimensional displacement of an isolation bearing, characterized in that, The method comprises the following steps: obtaining acceleration data of the isolation bearing in different axial directions and a frequency spectrum of each axial direction within a preset time period; screening a main frequency of an earthquake according to an amplitude value distribution of each frequency in the frequency spectrum of each axial direction, and obtaining a relative energy value of each axial direction according to an amplitude value proportion of the main frequency of the earthquake in the frequency spectrum of each axial direction; correcting the acceleration data of each axial direction according to a correlation between the acceleration data of each axial direction and other axial directions, in combination with a distribution of the acceleration data of each axial direction and the relative energy value, to obtain corrected data of each axial direction; determining three-dimensional displacement data of the isolation bearing based on an integral result of the corrected data of each axial direction.

2. The method for high-precision monitoring of three-dimensional displacement of an isolation bearing according to claim 1, characterized in that, The main frequency of the earthquake is screened according to the amplitude value distribution of each frequency in the frequency spectrum of each axial direction, and specifically comprises: obtaining a distribution characteristic value of each frequency in the frequency spectrum of each axial direction according to an amplitude value distribution of each frequency and an adjacent frequency in combination with a value of each frequency; obtaining a main frequency index of each frequency according to the distribution characteristic value and the amplitude value distribution of each frequency corresponding to all axial directions; taking a frequency corresponding to a maximum value of the main frequency index as the main frequency of the earthquake.

3. The method of claim 2, wherein the method is characterized by, The distribution characteristic value of each frequency in the frequency spectrum of each axial direction is obtained according to the amplitude value distribution of each frequency and the adjacent frequency in combination with the value of each frequency, and specifically comprises: for the frequency spectrum of any one axial direction, taking any one frequency as a selected frequency; arranging all frequencies in a preset order according to the value, and obtaining a reference frequency adjacent to the selected frequency in the arrangement order of the frequencies; determining an amplitude comparison coefficient based on the amplitude value of the selected frequency and the amplitude value of the reference frequency, and determining a frequency characteristic coefficient based on a negative correlation coefficient corresponding to the value of the selected frequency; taking a product of the amplitude comparison coefficient and the frequency characteristic coefficient as the distribution characteristic value of the selected frequency in the frequency spectrum of the any one axial direction.

4. The method for high-precision monitoring of three-dimensional displacement of an isolation bearing according to claim 3, characterized in that, The main frequency index of each frequency is obtained according to the distribution characteristic value and the amplitude value distribution of each frequency corresponding to all axial directions, and specifically comprises: taking a product of the amplitude value and the distribution characteristic value of the selected frequency in the frequency spectrum of each axial direction as a main frequency factor of the selected frequency of each axial direction, and taking a cumulative sum of the main frequency factors of all axial directions at the selected frequency as the main frequency index of the selected frequency.

5. The method of claim 1, wherein the method is used for high-precision monitoring of three-dimensional displacement of an isolation bearing. The acceleration data of each axial direction is corrected according to the correlation between the acceleration data of each axial direction and other axial directions in combination with the distribution of the acceleration data of each axial direction and the relative energy value, to obtain the corrected data of each axial direction, and specifically comprises: obtaining a data usability factor of each axial direction according to a comprehensive result of the correlation between the acceleration data of each axial direction and other axial directions; analyzing a data relative distribution of each axial direction according to the data usability factor and a balance of the acceleration data of the same axial direction, and obtaining an earthquake wave effective factor of each axial direction in combination with a difference between the data relative distribution and the relative energy value; correcting the acceleration data of each axial direction by using the earthquake wave effective factor of each axial direction to obtain the corrected data of each axial direction.

6. The method of claim 5, wherein the method is used for high-precision monitoring of three-dimensional displacement of an isolation bearing. The data available factor of each axial direction is obtained according to the correlation between the acceleration data of each axial direction and the acceleration data of every other axial direction, and specifically includes the following steps. Any axial direction is selected as a selected axial direction, and every other axial direction except the selected axial direction is selected as a reference axial direction. A correlation coefficient between the acceleration data of the selected axial direction and the acceleration data of each reference axial direction is calculated, and an average of the correlation coefficients between the selected axial direction and all reference axial directions is taken as the data available factor of the selected axial direction.

7. The method of claim 6, wherein the method is used for high-precision monitoring of three-dimensional displacement of an isolation bearing. The data relative distribution of each axial direction is analyzed according to the data available factor of each axial direction and the equalization of the acceleration data of the same axial direction, and the seismic wave effective factor of each axial direction is obtained by combining the difference between the data relative distribution and the relative energy value, and specifically includes the following steps. A product of the average of all acceleration data of the selected axial direction and the data available factor is taken as a response characteristic index of the selected axial direction, and a ratio between the response characteristic index of the selected axial direction and a cumulative sum of the response characteristic indexes of all axial directions is taken as a relative contribution degree of the selected axial direction. A difference between the relative contribution degree of the selected axial direction and the relative energy value is negatively correlated and normalized to obtain the seismic wave effective factor of the selected axial direction. 8.The method for monitoring three-dimensional displacement of the seismic isolation bearing according to claim 7, wherein, The acceleration data of each axial direction is corrected by using the seismic wave effective factor of each axial direction to obtain corrected data of each axial direction, and specifically includes the following steps. A product between the seismic wave effective factor of the selected axial direction and each acceleration data of the selected axial direction is taken as corrected data corresponding to each acceleration data of the selected axial direction. 9.The method for monitoring three-dimensional displacement of seismic isolation bearing according to claim 1, wherein, The relative energy value of each axial direction is obtained according to the amplitude proportion of the seismic main frequency of each axial direction, and specifically includes the following steps. A ratio between the amplitude of the seismic main frequency of each axial direction and a cumulative sum of the amplitudes of the seismic main frequencies of all axial directions is taken as the relative energy value of each axial direction.

10. A high-precision monitoring system for three-dimensional displacement of an isolation bearing, comprising a memory, a processor, and a computer program stored on the memory and running on the processor, characterized in that, The computer program is executed by the processor to realize the steps of the three-dimensional displacement high-precision monitoring method for the seismic isolation support according to any one of claims 1-9.

Images