A time-frequency domain joint inversion-based isolation bearing displacement trajectory reconstruction method

Abstract

The application discloses a kind of based on time-frequency domain joint inversion isolation bearing displacement trajectory reconstruction method, comprising: with the original vibration response of isolation bearing as input, according to the energy migration order of isolation bearing in reciprocating shear motion loading section, unloading section and residual stay section, inversion obtains bearing hysteresis response fragment chain;With bearing hysteresis response fragment chain as input, time domain integration inversion and frequency domain phase inversion are carried out in each bearing hysteresis response fragment, and according to the reverse connection relationship between loading section and unloading section, identify false displacement component and eliminate, retain real displacement component, generate bearing real displacement fragment;According to the connection order of bearing hysteresis response fragment chain, adjacent bearing real displacement fragment is locked and residual displacement is inherited, and the reconstruction result of isolation bearing displacement trajectory is obtained.The application can distinguish real displacement component and drift component, realize the reconstruction of the displacement trajectory of the isolation bearing with multiple hysteresis cycles.

Description

Technical Field

[0001] This invention relates to the field of seismic isolation monitoring technology, and in particular to a method for reconstructing the displacement trajectory of seismic isolation supports based on joint time-frequency domain inversion. Background Technology

[0002] Seismic isolation technology is a crucial technique in the seismic design of building structures, bridge structures, and important engineering facilities. By installing seismic isolation bearings between the superstructure and foundation, it enables the structure to generate controlled relative displacement under seismic loading, thereby reducing the transmission of seismic energy to the superstructure. During an earthquake, seismic isolation bearings typically undergo continuous hysteretic motion processes, including loaded slip, unloaded rebound, and post-earthquake residual residence. Their displacement trajectories directly reflect the energy dissipation performance, recovery capacity, and post-earthquake residual deformation state of the seismic isolation bearings. Therefore, accurately obtaining the displacement trajectories of seismic isolation bearings is of great significance for the safety assessment of seismically isolated structures, post-earthquake damage determination, and the operation and maintenance management of seismic isolation bearings.

[0003] Currently, displacement monitoring of seismic isolation bearings mainly employs direct displacement measurement using displacement gauges or displacement inversion methods based on acceleration signal integration. Direct displacement measurement suffers from limitations such as limited installation space, high equipment costs, and insufficient long-term operational reliability. In contrast, displacement inversion based on acceleration signals is widely used due to its ease of sensor deployment and broad applicability. However, acceleration signals are susceptible to sensor zero drift, baseline offset, and environmental noise during integration, causing integration errors to accumulate over time and ultimately resulting in a noticeable low-frequency drift in the displacement trajectory. Therefore, existing technologies typically employ high-pass filtering, baseline correction, or frequency domain filtering to eliminate this drift.

[0004] However, seismic isolation bearings differ significantly from ordinary vibrating structures, exhibiting typical hysteretic motion characteristics in their displacement response. During a complete seismic response, a seismic isolation bearing undergoes multiple motion stages, including loading, unloading, and residual dwell phases, with clear directional transformation and residual displacement transfer relationships between these stages. Existing technologies mostly treat the entire monitoring process as a continuous vibration signal for unified processing, failing to fully utilize the hysteretic motion characteristics of the seismic isolation bearing itself to constrain the displacement inversion results. Therefore, when eliminating low-frequency drift, the actual low-frequency slip displacement of the seismic isolation bearing is easily misjudged as a drift component and filtered out. Simultaneously, the actual residual displacement after the earthquake is easily incorrectly zeroed out during baseline correction, resulting in inconsistencies between the reconstructed displacement trajectory and the actual motion process of the seismic isolation bearing. Summary of the Invention

[0005] This invention provides a method for reconstructing the displacement trajectory of seismic isolation bearings based on joint inversion in the time and frequency domains. By utilizing the hysteretic motion laws of the seismic isolation bearings during loading, unloading, and residual dwell processes, it can distinguish between the true displacement components and the drift components, maintain the continuous transmission of residual displacement, and realize the reconstruction of the displacement trajectory of seismic isolation bearings with multiple hysteresis cycles.

[0006] Firstly, a method for reconstructing the displacement trajectory of seismic isolation bearings based on joint time-frequency domain inversion is provided, including the following steps:

[0007] S1. Taking the original vibration response of the seismic isolation bearing as input, and based on the energy migration order of the loading segment, unloading segment and residual dwell segment of the seismic isolation bearing in the reciprocating shear motion, the original vibration response is inverted into a seismic isolation bearing hysteresis response segment chain. The bearing hysteresis response segment chain is used to characterize the continuous motion process of the seismic isolation bearing from the loading segment, unloading segment to the residual dwell segment.

[0008] S2. Using the support hysteresis response segment chain as input, perform time-domain integral inversion and frequency-domain phase inversion in each support hysteresis response segment, and identify and remove the pseudo displacement components formed by sensor zero drift according to the reverse connection relationship between the loading segment and the unloading segment, retain the real displacement components formed by the shear deformation of the seismic isolation support, and generate the real displacement segment of the support.

[0009] S3. Using the actual displacement segment of the support as input, according to the sequence of the support hysteresis response segment chain, perform hysteresis endpoint locking and residual displacement inheritance on the actual displacement segments of adjacent supports to obtain the displacement trajectory reconstruction result of the seismic isolation support.

[0010] Optionally, the original vibration response is used as input, and a short-time energy sequence is calculated based on the original vibration response to obtain the energy migration envelope. On the energy migration envelope, the continuous rising edge of energy, the continuous falling edge of energy, and the low-energy segment where the energy is stably lower than a preset threshold are detected, which correspond to the energy injection, energy release, and residual residence state after energy dissipation of the seismic isolation bearing in the reciprocating shear motion, respectively.

[0011] Optionally, the period along which the energy continuously rises is identified as the loading segment, which characterizes the process of the seismic isolation bearing sliding under force in the shear direction and the displacement increasing in the same direction; the period along which the energy continuously falls is identified as the unloading segment, which characterizes the process of the seismic isolation bearing rotating in the sliding direction and the displacement rebounding in the opposite direction; the low-energy segment is identified as the residual dwell segment, which characterizes the process of the seismic isolation bearing maintaining residual displacement and the force almost disappearing after completing one shear rebound.

[0012] Optionally, the loading segment, unloading segment, and residual dwell segment are sequentially spliced ​​together according to the time sequence of the start time of the loading segment, unloading segment, and residual dwell segment to form a support hysteresis response segment chain consisting of alternating loading segment, unloading segment, and residual dwell segment; the temporal connection of the support hysteresis response segment chain directly maps to the entire continuous motion process of the seismic isolation support from force slippage and directional rotation to residual dwell.

[0013] Optionally, in step S2, the current support hysteresis response segment is extracted from the support hysteresis response segment chain, the acceleration signal within the current support hysteresis response segment is obtained, and the acceleration signal within the current support hysteresis response segment is subjected to time-domain integration inversion to obtain the initial displacement trajectory including low-frequency displacement components and true displacement components.

[0014] Optionally, in step S2, based on the current support hysteresis response segment, the loading segment and unloading segment that belong to the same complete hysteresis cycle as the current support hysteresis response segment are extracted from the support hysteresis response segment chain, and the reverse connection relationship between the loading segment and the unloading segment is determined. The reverse connection relationship means that in one complete hysteresis cycle of the seismic isolation bearing, the shear displacement of the bearing in the loading segment increases along the first direction, and the shear displacement of the bearing in the unloading segment rebounds along the second direction opposite to the first direction, forming a displacement direction reversal connection between loading and unloading.

[0015] Optionally, using the reverse continuity relationship, the low-frequency displacement components in the initial displacement trajectories of the loading and unloading segments are compared for directional consistency: if a certain low-frequency displacement component maintains the same cumulative trend in the loading and unloading segments, it is determined to be a low-frequency pseudo-displacement component formed by sensor zero drift and is removed; all the retained real displacement components are reconstructed to generate the real displacement segment of the support corresponding to the support hysteresis response segment.

[0016] Optionally, the support hysteresis response segment chain includes a loading segment, an unloading segment, and a residual dwell segment arranged in a time sequence. Based on the time sequence of the loading segment, unloading segment, and residual dwell segment in the support hysteresis response segment chain, the connection type between each support real displacement segment and its preceding and following segments is determined. The connection type includes loading segment to unloading segment, unloading segment to residual dwell segment, and residual dwell segment to the next loading segment. Based on the connection type of adjacent segments in the support hysteresis response segment chain, the displacement state formed at the end of the previous support real displacement segment is transferred to the next support real displacement segment. Hysteresis endpoint locking and residual displacement inheritance are performed to maintain displacement continuity of each support real displacement segment within the hysteresis cycle and between adjacent hysteresis cycles.

[0017] Optionally, the execution of hysteresis endpoint locking and residual displacement inheritance specifically includes:

[0018] For adjacent segments where the loading segment connects to the unloading segment, the displacement at the end of the loading segment is set as the extreme displacement reference at the beginning of the unloading segment, and hysteresis endpoint locking is performed: the starting displacement of the actual displacement segment of the unloading segment is corrected with this extreme displacement reference, so that the displacement trajectory of the unloading segment continuously rebounds in the opposite direction from the maximum displacement point at the end of the loading segment, eliminating the endpoint misalignment caused by segmentation.

[0019] For adjacent segments of the unloading segment and the residual dwell segment, the displacement at the end of the unloading segment is identified as the residual displacement of this hysteresis cycle, and it is transferred to the real displacement segment of the residual dwell segment as a constant displacement reference, and residual displacement inheritance is performed.

[0020] For the adjacent segment of the residual dwell segment following the next loading segment, the inherited residual displacement is used as the starting displacement of the actual displacement segment of the next loading segment, so as to realize the continuous transfer of residual displacement between adjacent hysteresis cycles.

[0021] Optionally, all the actual displacement segments of the supports after hysteresis endpoint locking and residual displacement inheritance are spliced ​​together end to end according to the sequence of the support hysteresis response segment chain to obtain a reconstruction result of the seismic isolation support displacement trajectory that is continuous on the time axis and has no displacement jump.

[0022] Secondly, a seismic isolation bearing displacement trajectory reconstruction system based on joint time-frequency domain inversion is provided, including:

[0023] The bearing hysteresis response segment chain generation module is used to take the original vibration response of the seismic isolation bearing as input and, based on the energy migration order of the loading segment, unloading segment and residual dwell segment of the seismic isolation bearing in reciprocating shear motion, invert the original vibration response into a bearing hysteresis response segment chain. The bearing hysteresis response segment chain is used to characterize the continuous motion process of the seismic isolation bearing from the loading segment, unloading segment to the residual dwell segment.

[0024] The support real displacement segment generation module is used to take the support hysteresis response segment chain as input, perform time-domain integral inversion and frequency-domain phase inversion in each support hysteresis response segment, and identify and remove the pseudo displacement components formed by sensor zero drift according to the reverse connection relationship between the loading segment and the unloading segment, retain the real displacement components formed by the shear deformation of the seismic isolation support, and generate the support real displacement segment.

[0025] The seismic isolation bearing displacement trajectory reconstruction module is used to take the actual displacement segment of the bearing as input, and according to the successive order of the hysteresis response segment chain of the bearing, perform hysteresis endpoint locking and residual displacement inheritance on the actual displacement segments of adjacent bearings to obtain the seismic isolation bearing displacement trajectory reconstruction result.

[0026] Thirdly, an electronic device is provided, comprising:

[0027] Memory, on which computer programs / instructions are stored;

[0028] When the processor executes the computer program / instruction, it implements the steps of the seismic isolation bearing displacement trajectory reconstruction method based on time-frequency domain joint inversion as described above.

[0029] Fourthly, a computer-readable storage medium is provided, on which a computer program / instruction is stored, wherein when the computer program / instruction is executed by a processor, the steps of the seismic isolation support displacement trajectory reconstruction method based on time-frequency domain joint inversion as described above are implemented.

[0030] The present invention has the following beneficial effects:

[0031] This invention identifies the loading, unloading, and residual dwell stages based on the energy migration envelope, and constructs a chain of bearing hysteretic response segments. This decomposes the complete hysteretic motion process of the seismic isolation bearing into temporal motion units with clear physical meaning. By using this chain of bearing hysteretic response segments to characterize the entire process of the seismic isolation bearing from stress-induced slippage and directional return to residual dwell, subsequent displacement inversion is based on actual hysteretic motion laws. This enhances the consistency between the displacement trajectory reconstruction results and the actual stress state of the seismic isolation bearing, and improves the physical reliability of trajectory reconstruction under complex seismic conditions.

[0032] This invention decomposes the initial displacement trajectory into multiple candidate low-frequency displacement components. Combining this with the motion law that the displacement directions inevitably reverse during the loading and unloading phases, the directional consistency of each candidate low-frequency displacement component is determined. Low-frequency displacement components that maintain a consistent cumulative trend in both the loading and unloading phases are identified as pseudo-displacement components, while those that reverse with the support's motion direction are identified as true displacement components. Therefore, drift suppression and preservation of true displacement can be achieved without relying on empirical filtering parameters. This method utilizes the hysteretic characteristics of the seismic isolation bearing itself as the criterion.

[0033] This invention, based on the temporal sequence of the bearing hysteresis response segment chain, performs hysteresis endpoint locking between the loading and unloading segments, and inherits residual displacement between the unloading and residual dwell segments. The residual displacement is then passed to the next loading segment as the starting displacement reference, ensuring that the actual residual displacement formed after each hysteresis cycle is continuously retained and accumulated in subsequent cycles. The resulting reconstructed displacement trajectory of the seismic isolation bearing is not only continuous without jumps on the time axis, but also fully reflects the displacement evolution history of the seismic isolation bearing during multiple loading, unloading, and residual dwell processes. This provides more reliable data for post-earthquake performance evaluation of seismic isolation bearings, residual deformation analysis, and structural safety status determination. Attached Figure Description

[0034] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0035] Figure 1 This is a flowchart of the seismic isolation bearing displacement trajectory reconstruction method based on time-frequency domain joint inversion provided by an embodiment of the present invention;

[0036] Figure 2 This is a schematic diagram of step S2 execution provided in an embodiment of the present invention. Detailed Implementation

[0037] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be described in detail below. Obviously, the described embodiments are merely some embodiments of this invention, and not all embodiments. Based on the embodiments of this invention, all other implementation methods obtained by those skilled in the art without creative effort are within the scope of protection of this invention.

[0038] like Figures 1-2 As shown, a method for reconstructing the displacement trajectory of seismic isolation bearings based on joint time-frequency domain inversion is presented. This invention does not employ the traditional method of simply acquiring data, extracting features, and fusing them for reconstruction. Instead, it uses the unique loading, unloading, and residual dwell time of the seismic isolation bearings as the main inversion parameters. Through hysteresis response fragment chain constraints in time-frequency domain inversion, it addresses the problems in existing technologies that easily mistake post-earthquake residual displacement for integral drift elimination and unloading return phase lag for trajectory errors. The method includes the following steps:

[0039] S1. Taking the original vibration response of the seismic isolation bearing as input, and based on the energy migration sequence of the loading segment, unloading segment, and residual dwell segment of the seismic isolation bearing in the reciprocating shear motion, the original vibration response is inverted into a hysteresis response segment chain of the bearing. The hysteresis response segment chain of the bearing is used to characterize the continuous motion process of the seismic isolation bearing from the force-induced slippage and directional rotation to the residual dwell segment (i.e. from the loading segment, unloading segment to the residual dwell segment).

[0040] S11. Obtain the original vibration response of the seismic isolation bearing, and perform sliding time window energy statistics on the original vibration response to construct the corresponding short-time energy sequence; generate the energy migration envelope based on the continuous change trend of the short-time energy sequence on the time axis.

[0041] Specifically, the original vibration response is segmented and statistically analyzed according to a preset time window, and the response energy within each time window is accumulated to obtain the short-time energy value at the corresponding moment:

[0042] ;

[0043] in, Indicates the first The short-time energy value corresponding to each time window Indicates the first Within the first time window The original vibration response value corresponding to each sampling time refers to the vibration response sampling value collected by sensors installed near the seismic isolation bearing. This invention requires time-domain integral inversion and frequency-domain phase inversion, and acceleration response values ​​are preferred. Therefore... For the first Within the first time window The vibration acceleration sampled value of the seismic isolation bearing at each sampling time. This indicates the number of sampling points contained within a single time window.

[0044] Subsequently, the continuously acquired short-time energy values ​​are arranged in chronological order to form a short-time energy sequence:

[0045] ;in, Represents a short-time energy sequence. This represents the total number of time windows in a short-time energy sequence.

[0046] To further mitigate the impact of local random fluctuations on energy trend identification, the short-time energy sequence is smoothed to obtain the energy migration envelope:

[0047] ;

[0048] in, Indicates the first The energy migration envelope value corresponding to each time window Indicates the smooth window half-width. Indicates the first The short-term energy value corresponding to each time window.

[0049] The energy migration envelope is used to reflect the overall migration trend of energy injection, energy release, and energy dissipation of the seismic isolation bearing during the entire reciprocating shear motion.

[0050] S12. Based on the energy migration envelope, identify the continuously rising energy edge, the continuously falling energy edge, and the low-energy segment, and form the loading segment, unloading segment, and residual dwell segment respectively. This includes scanning the energy migration envelope point by point, calculating the change in the energy migration envelope between adjacent time moments, and identifying different motion stages based on the direction, duration, and current energy level of the change.

[0051] First, for any given moment, calculate the energy change between the current moment and the previous moment. When the energy change is consistently positive across multiple consecutive time windows, and the corresponding energy migration envelope shows a continuous upward trend, it indicates that the external load is continuously inputting energy into the seismic isolation bearing, the shear deformation inside the bearing is accumulating, and the bearing displacement is continuously increasing along the current slip direction. At this point, this continuous growth interval is defined as the continuous energy rise edge, and the time interval corresponding to this energy rise edge is marked as the loading segment. To avoid misjudgment due to local noise, the continuous growth state must last for at least 3-4 consecutive time windows, and no significant energy drop is allowed between adjacent time windows (judged by whether the energy drop value exceeds a preset drop threshold). When the growth state ends, the start time of the growth is taken as the starting point of the loading segment, and the end time is taken as the ending point, forming a complete loading segment.

[0052] Subsequently, after the loading segment ends, the search continues along the time axis. When the energy change is consistently negative within multiple consecutive time windows and the overall energy migration envelope shows a downward trend, it indicates that the seismic isolation bearing no longer continuously absorbs external energy but begins to release the deformation energy stored in the previous period. The sliding direction of the seismic isolation bearing gradually reverses and a displacement rebound occurs. At this time, the continuous downward interval is identified as the continuous energy downward edge, and the corresponding time interval is marked as the unloading segment. To prevent misidentification caused by instantaneous fluctuations, it is preferable to require the continuous downward state to last for more than 3 to 4 consecutive time windows. When the downward trend ends, the starting time of the downward movement is taken as the starting point of the unloading segment, and the ending time of the downward movement is taken as the ending point of the unloading segment, forming a complete unloading segment.

[0053] After the unloading phase ends, the energy migration envelope continues to be monitored. When the energy migration envelope drops below the preset low energy threshold and remains in an approximately constant state for multiple consecutive time windows, it indicates that the seismic isolation bearing has completed one shear rebound process and there is no obvious energy input and energy release behavior. Only the post-earthquake residual displacement is retained. At this time, the continuous low energy state is identified as the low energy segment, and the corresponding time interval is marked as the residual dwell segment.

[0054] After completing the above identification, the loading segment, unloading segment, and residual dwell segment are associated and organized in chronological order. For each loading segment, a corresponding relationship is established with the immediately following unloading segment; for each unloading segment, a corresponding relationship is established with the immediately following residual dwell segment, forming a complete motion unit including the loading segment, unloading segment, and residual dwell segment. When a new energy rising edge appears again, a new loading segment is established after the current residual dwell segment, and the above process is repeated, ultimately forming a support hysteresis response segment chain composed of multiple complete motion units including loading segment, unloading segment, and residual dwell segment connected end to end.

[0055] The loading segment characterizes the process of the seismic isolation bearing continuously sliding under shear force with displacement increasing in the same direction; the unloading segment characterizes the process of the seismic isolation bearing releasing deformation energy and undergoing directional rotation; and the residual dwell segment characterizes the process of the seismic isolation bearing maintaining residual displacement and almost disappearing force after completing one rebound. Through this division, the complete hysteretic motion process of the seismic isolation bearing can be mapped into a structured chain of bearing hysteretic response segments, providing a temporal constraint basis for subsequent inversion of the actual displacement segments of the bearing.

[0056] More specifically:

[0057] The change in the energy migration envelope between adjacent time points is expressed as follows: ;in, Indicates the first The amount of energy change corresponding to each time window Indicates the first The energy migration envelope value corresponding to each time window This represents the energy migration envelope value corresponding to the previous time window.

[0058] S121. When multiple consecutive time windows satisfy: At that time, the corresponding time interval is identified as the continuous rising edge of energy. Among them, This indicates the energy rise threshold, used to identify whether the seismic isolation bearing has entered a state of continuous energy absorption. First, the total energy change during the current monitoring period is statistically analyzed. The absolute average level of the energy level is used to determine the energy rise threshold, which is set to 1.2 to 1.4 times the average level, preferably 1.3 times. When the seismic isolation bearing is in a static or slightly vibrating state, the energy migration envelope will also produce small random fluctuations. If the threshold is too low, noise fluctuations are easily misjudged as the loading process. When the energy change exceeds the overall average fluctuation level by about 30%, it means that the external load has begun to continuously input energy into the bearing, and the bearing has entered the real shear slip stage. Therefore, the starting position of the loading segment can be identified more accurately.

[0059] The process of continuously increasing energy input to the seismic isolation bearing along the corresponding external load is therefore identified as the loading segment; the loading segment characterizes the process of the bearing continuously being subjected to force and sliding in the shear direction, with the displacement increasing in the same direction.

[0060] S122. When multiple consecutive time windows satisfy: At that time, the corresponding time interval is identified as the continuously decreasing energy edge. This indicates the energy reduction threshold, used to identify whether the seismic isolation bearing has entered a state of continuous energy release. Preferably, it adopts... Using the same adaptive determination method, when the seismic isolation bearing rebounds, the stored deformation energy will be continuously released, and the energy migration envelope will show a continuous downward trend. The decrease caused by normal measurement noise is usually small. When the decrease exceeds the overall fluctuation level by about 30%, it can be considered that the bearing has changed from the stress sliding stage to the rebound energy release stage, ensuring the stability and accuracy of the unloading section identification.

[0061] The process of energy continuously decreasing and releasing the stored deformation energy along the corresponding seismic isolation bearing is therefore identified as the unloading section; the unloading section characterizes the process of the bearing sliding direction gradually turning back and the displacement rebounding.

[0062] S123, When the following conditions are met: Furthermore, if this state persists for more than a preset duration, the corresponding time interval will be identified as a low-energy segment. This indicates the threshold for determining the low-energy section, used to identify whether the seismic isolation bearing has entered a residual dwell state. Preferably, [the threshold is...]. The threshold is set to 5% to 15% of the maximum energy migration envelope within the current monitoring period, with 10% being preferred. After the seismic isolation bearing completes one loading and unloading cycle, its internal energy has been largely dissipated, leaving only residual displacement. At this point, the energy migration envelope typically drops to around 10% of the peak energy or even lower. If the threshold is set too high, the unloading process that has not yet completed its rebound may be misidentified as a residual dwelling segment; if the threshold is set too low, the true residual dwelling state may be missed. Therefore, setting the threshold to around 10% of the peak energy can better distinguish between the continuous motion state and the residual dwelling state, which is consistent with the actual law of the hysteretic energy dissipation process of the seismic isolation bearing.

[0063] The low-energy section indicates that the seismic isolation bearing has completed one shear rebound, and the system's input and released energy tend to disappear. Therefore, it is identified as the residual dwell section. The residual dwell section is used to characterize the state in which the bearing maintains residual displacement and the force is almost gone.

[0064] S13. Organize the loading segment, unloading segment and residual dwell segment in chronological order to generate a support hysteresis response fragment chain.

[0065] Specifically, the start and end times of each loading segment, unloading segment, and residual dwell segment are extracted: ;in, Indicates the first A response fragment, Indicates the first The start time of each response segment Indicates the first The end time of a response segment Indicates the first The type identifier for each response fragment is either a loading segment, an unloading segment, or a residual dwell segment.

[0066] Then they are arranged in ascending order of their starting times: ;in, This represents a chain of support hysteresis response segments. This indicates the total number of response fragments.

[0067] When adjacent response segments satisfy the temporal sequence of "loading segment → unloading segment → residual dwell segment", a corresponding hysteresis correlation is established. When subsequent seismic action occurs again, a new loading segment continues to form after the residual dwell segment, and further forms a new unloading segment and a residual dwell segment, thus forming a support hysteresis response segment chain composed of alternating loading segments, unloading segments and residual dwell segments.

[0068] The support hysteresis response fragment chain fully records the continuous motion process of the seismic isolation support from stress slippage and directional rotation to residual dwell, and serves as the unified input basis for subsequent time-domain integral inversion and frequency-domain phase inversion.

[0069] S2. Using the support hysteresis response segment chain as input, perform time-domain integral inversion and frequency-domain phase inversion within each support hysteresis response segment. Based on the reverse connection between the loading segment and the unloading segment, identify and remove the pseudo displacement components formed by sensor zero drift, retain the real displacement components formed by the shear deformation of the seismic isolation support, and generate the real displacement segment of the support.

[0070] S21. Extract the current support hysteresis response segment from the support hysteresis response segment chain, and obtain the acceleration signal within the current support hysteresis response segment. Specifically, this includes reading the start and end times of the current support hysteresis response segment from the support hysteresis response segment chain, and then extracting the acceleration sampling data within this time interval from the original vibration response, thus obtaining the acceleration signal sequence corresponding to the current support hysteresis response segment. Let the start time of the current support hysteresis response segment be... The current end time of the hysteresis response segment of the support is Then, extract from the original vibration response: The acceleration signal sequence corresponding to the current support hysteresis response segment can be obtained. Time-domain integral inversion and frequency-domain phase inversion are performed on the acceleration signal in the current support hysteresis response segment to obtain the initial displacement trajectory and frequency-domain phase information.

[0071] Specifically, obtain the acceleration signal sequence corresponding to the current support hysteresis response segment:

[0072] ;in, This represents the acceleration signal sequence corresponding to the current support hysteresis response segment. Indicates the first The acceleration values ​​corresponding to each sampling time point.

[0073] First, perform time-domain integral inversion on the acceleration signal to obtain the velocity sequence: ;in, This represents the velocity sequence obtained by integration.

[0074] Then, by performing a second integration, the initial displacement trajectory is obtained: ;in, This represents the initial displacement trajectory.

[0075] Due to the influence of zero-bias error, installation error, and environmental noise on the accelerometer, the initial displacement trajectory simultaneously contains both true displacement components and low-frequency pseudo-displacement components. Therefore, frequency domain decomposition processing is performed on the initial displacement trajectory to decompose it into multiple displacement components corresponding to different frequency ranges. This includes performing frequency domain decomposition processing on the initial displacement trajectory to separate the vibration components within different frequency ranges from the overall displacement trajectory. Preferably, a frequency domain bandwidth partitioning method is used, dividing the spectrum corresponding to the initial displacement trajectory into multiple continuous frequency bands according to frequency from low to high, with each frequency band corresponding to an independent displacement component. First, the initial displacement trajectory is transformed in the frequency domain to obtain the corresponding spectral distribution. The spectrum is then divided into multiple frequency intervals according to frequency intervals. For each frequency interval, only the frequency components within that interval are retained, and the remaining frequency components are set to zero, thus forming the corresponding frequency band spectrum. Each frequency interval corresponds to the displacement change characteristics within different frequency ranges, thereby decomposing the initial displacement trajectory into multiple displacement components corresponding to different frequency ranges.

[0076] After completing the frequency domain decomposition, the frequency range below the low-frequency cutoff frequency is extracted as candidate low-frequency bands: ;in, Indicates candidate low-frequency bands. The low-frequency cutoff frequency is preferably determined based on the principal vibration frequency corresponding to the hysteresis response segment of the support. This includes performing spectral statistics on the current hysteresis response segment of the support to determine the principal vibration frequency with the highest energy proportion. Considering that the actual shear slip motion of the seismic isolation support is mainly concentrated in the frequency range near and below the principal vibration frequency, while the zero drift, integral drift, and baseline drift of the sensor usually exhibit extremely low frequency changes far below the principal vibration frequency, the low-frequency cutoff frequency is preferably set to 20% to 40% of the principal vibration frequency, preferably 30%.

[0077] For each low-frequency component in the candidate low-frequency band, the corresponding low-frequency shift component is reconstructed. ; Indicates the first A low-frequency displacement component This indicates the low-frequency displacement component number. Specifically, it involves performing separate reconstruction processing on each frequency interval of the candidate low-frequency band. For the... For each low-frequency interval, only the spectral components corresponding to that interval are retained. Then, an inverse transform operation is performed on the spectrum. Through this process, each low-frequency interval corresponds to an independent low-frequency displacement component. Since different low-frequency displacement components correspond to displacement change characteristics at different time scales, some low-frequency displacement components may originate from the actual shear slip process of the seismic isolation bearing, while others may originate from sensor zero drift or integral drift. Finally, all the reconstructed low-frequency displacement components are arranged in frequency order to form a candidate low-frequency displacement component set. ;in, This represents the set of candidate low-frequency displacement components. This represents the total number of candidate low-frequency displacement components.

[0078] The candidate low-frequency displacement component set is used to identify the true displacement components and pseudo displacement components by combining the reverse connection relationship between the loading segment and the unloading segment.

[0079] S22. Construct the reverse connection between the loading segment and the unloading segment.

[0080] Specifically, based on the current support hysteresis response segment, the loading and unloading segments that belong to the same complete hysteresis cycle as the current support hysteresis response segment are extracted from the support hysteresis response segment chain. A complete hysteresis cycle refers to a cycle unit in the support hysteresis response segment chain, which consists of a loading segment, an unloading segment, and a residual dwell segment.

[0081] Calculate the overall direction of change of the displacement trajectories corresponding to the loading and unloading segments respectively:

[0082] ;in, Indicates the direction of displacement increase of the loaded segment. Indicates the start time of the loading segment. This indicates the end time of the loading segment.

[0083] when: When, it indicates that the displacement of the loaded segment increases along the first direction;

[0084] when: When, it indicates that the displacement of the loaded segment increases in the opposite direction.

[0085] Similarly, calculate the direction of displacement change of the unloading section: ;in, This indicates the direction of displacement change in the unloading section. Indicates the start time of the unloading segment. This indicates the end time of the unloading segment.

[0086] If the following conditions are met: If so, it is determined that there is a reverse connection between the loading segment and the unloading segment.

[0087] in, This indicates that the two segments change in opposite directions.

[0088] The reverse continuity relationship characterizes that during a complete hysteresis cycle, the displacement of the seismic isolation bearing first increases along the first direction and then rebounds along the second direction opposite to the first direction, thus forming a characteristic of reverse displacement direction continuity between loading and unloading, that is, the displacement change directions of the loading segment and the unloading segment are opposite.

[0089] S23. Identify pseudo-displacement components and true displacement components based on the reverse continuity relationship.

[0090] Specifically, each low-frequency displacement component in the candidate low-frequency displacement component set is analyzed separately.

[0091] For any low-frequency displacement component: Calculate its cumulative change trend in the loading and unloading segments respectively:

[0092] ;in, Indicates the first The cumulative change of a low-frequency displacement component in the loading segment.

[0093] ;in, Indicates the first The cumulative change of a low-frequency displacement component in the unloading section.

[0094] When the following conditions are met: This indicates that the low-frequency displacement component maintains a trend of increasing or decreasing in the same direction in both the loading and unloading segments.

[0095] Since the actual motion of the seismic isolation bearing will inevitably reverse direction during loading and unloading, such low-frequency components that always accumulate in the same direction do not conform to the actual hysteretic motion law of the bearing, and are therefore judged as pseudo-displacement components formed by sensor zero drift.

[0096] When the following conditions are met: This indicates that the low-frequency displacement component exhibits opposite trends during the loading and unloading phases. This characteristic is consistent with the actual shear slip and rebound process of the seismic isolation bearing, therefore it is identified as the true displacement component.

[0097] Then, a set of real displacement components is constructed. It includes all that satisfy low-frequency displacement components And all displacement components in the frequency band above the low-frequency cutoff frequency (all displacement components in the frequency band above the low-frequency cutoff frequency are directly retained as real displacement components).

[0098] S24. Reconstruct the true displacement segment of the support using the set of true displacement components.

[0099] Specifically, all true displacement components in the true displacement component set are superimposed and reconstructed:

[0100] ;

[0101] in, This represents a segment of the actual displacement of the support. Represents the set of true displacement components The j-th true displacement component in the equation.

[0102] The obtained true displacement segments of the bearings retain the true displacement change information formed during the loading, unloading and rebound of the seismic isolation bearings, while eliminating the pseudo displacement components caused by sensor zero drift, integral accumulation error and extremely low frequency baseline drift.

[0103] Finally, the actual displacement segment of the support is used as the input for step S3, which is then used for subsequent hysteresis endpoint locking and residual displacement inheritance processing.

[0104] S3. Using the actual displacement segment of the support as input, according to the sequence of the support hysteresis response segment chain, perform hysteresis endpoint locking and residual displacement inheritance on the actual displacement segments of adjacent supports to obtain the displacement trajectory reconstruction result of the seismic isolation support.

[0105] S31. Determine the connection type between the actual displacement segments of adjacent supports based on the support hysteresis response segment chain.

[0106] Specifically, the support hysteresis response fragment chain generated by S1 and the actual displacement fragments of each support generated by S2 are obtained.

[0107] The support hysteresis response fragment chain is represented as follows: ;in, This represents a chain of support hysteresis response segments. Indicates the first Hysteresis response segment of each support This indicates the total number of support hysteresis response segments.

[0108] Each support hysteresis response segment corresponds to one support actual displacement segment: ;in, Indicates the first Each support hysteresis response segment corresponds to a segment of the actual displacement of the support.

[0109] Further reading reveals the type identifiers of each segment in the support hysteresis response segment chain: ;in, Indicates the first Type identifier for each support hysteresis response segment Indicates the loaded segment. Indicates the unloading segment. Indicates the remaining dwell period.

[0110] Then, according to the time sequence in the support hysteresis response segment chain, for any two adjacent support hysteresis response segments: Perform type combination identification.

[0111] When the following conditions are met: When the connection type is determined to be a loading segment followed by an unloading segment, it indicates that the current segment is a loading segment and the next segment is an unloading segment. That is, after the seismic isolation bearing completes the force-induced sliding, it begins to rebound. In this stage, the bearing displacement first increases continuously in a certain direction. After reaching the displacement extreme point of this hysteresis cycle, the direction of motion begins to reverse and enter the rebound process. Therefore, this connection relationship corresponds to a loading segment followed by an unloading segment. Subsequently, hysteresis endpoint locking needs to be performed to ensure that the unloading segment starts continuously from the displacement extreme point of the loading segment, avoiding displacement discontinuities or jumps.

[0112] When the following conditions are met: When the current segment is unloading segment followed by residual dwell segment, the connection type is determined to be unloading segment followed by residual dwell segment. This indicates that the current segment is the unloading segment and the next segment is the residual dwell segment. That is, after the seismic isolation bearing completes its rebound, it enters a stable dwell state. In this stage, the bearing has released most of its deformation energy, and the displacement change tends to stop, but a certain residual displacement is still retained. Therefore, this connection relationship corresponds to unloading segment followed by residual dwell segment. Subsequently, the residual displacement formed at the end of the unloading segment needs to be transferred to the residual dwell segment as its constant displacement reference.

[0113] When the following conditions are met: When the current segment is determined to be a residual dwell segment followed by a new loading segment, the continuity type is determined to be a residual dwell segment followed by a new loading segment. This means that after the seismic isolation support maintains its residual displacement state, it is subjected to external loads again, starting a new round of shear slip process. In this stage, the new loading does not start from zero displacement, but continues from the residual displacement position left over from the previous cycle. Therefore, this continuity relationship corresponds to a residual dwell segment followed by a new loading segment. Subsequently, residual displacement inheritance needs to be performed, using the residual displacement retained from the previous cycle as the starting displacement of the next loading segment to ensure the displacement continuity between multiple hysteresis cycles.

[0114] Through the above identification process, the connection relationship between all the actual displacement segments of the supports is determined, which provides a basis for the subsequent execution of hysteresis endpoint locking and residual displacement inheritance.

[0115] S32. Perform hysteresis endpoint locking on adjacent segments of the loading segment and the unloading segment.

[0116] Specifically, for adjacent support displacement segments with a connection type of loading segment followed by unloading segment, the end displacement value of the loading segment is first extracted: ;in, This indicates the displacement value at the end of the loaded segment. This represents the actual displacement segment of the support corresponding to the loading segment. This indicates the end time of the loading segment.

[0117] The displacement value at the end of the loading segment corresponds to the extreme displacement point of the seismic isolation bearing during this stress-induced sliding process, and therefore it is set as the displacement extreme value benchmark.

[0118] Then, the starting displacement value of the actual displacement segment of the unloading section is extracted: ;in, This indicates the starting displacement value of the unloading section. This represents the actual displacement segment of the support corresponding to the unloading section. Indicates the start time of the unloading segment.

[0119] Calculate the displacement offset between the two: ;in, This indicates the displacement offset between the loading segment and the unloading segment.

[0120] Then, the displacement offset is used to perform an overall translational correction on the actual displacement segment of the entire unloading section:

[0121] ;in, This represents the actual displacement segment of the unloading section after the hysteresis endpoint is locked.

[0122] Through the above processing, the actual displacement segment of the unloading segment rebounds continuously in the opposite direction starting from the extreme displacement point at the end of the loading segment, thereby eliminating the end-point misalignment phenomenon caused by segmented processing.

[0123] S33. Perform residual displacement inheritance on adjacent segments of the unloading segment and the residual dwell segment.

[0124] Specifically, for adjacent support true displacement segments with a connection type of unloading segment followed by residual dwell segment, the displacement value corresponding to the end time of the unloading segment true displacement segment is extracted: ;in, This represents the residual displacement. This indicates the end time of the unloading segment.

[0125] Since the isolation bearing has completed the rebound process after the unloading phase, the residual displacement represents the permanent displacement actually retained after the end of this hysteresis cycle.

[0126] This residual displacement is then used as the constant displacement reference for the actual displacement segment of the residual dwell period:

[0127] ;in, This represents the actual displacement segment of the residual dwell period after the residual displacement is inherited.

[0128] Through the above processing, the residual displacement level of the residual dwell section is kept consistent with that of the previous unloading section, thereby truly reflecting the residual dwell state of the seismic isolation bearing after the force disappears.

[0129] S34. Perform cyclic displacement inheritance on the adjacent segments of the remaining stationary segment that are connected to the next loading segment.

[0130] Specifically, for adjacent support real displacement segments where the continuity type is a residual dwell segment followed by the next loading segment, the residual displacement formed in the previous hysteresis cycle is continued to be transferred to the next loading segment.

[0131] First, extract the residual displacement corresponding to the remaining dwell segment. Then, extract the starting displacement value of the actual displacement segment of the next loading segment: ;in, Indicates the starting displacement value of the next loading segment. This represents the actual displacement segment of the support corresponding to the next loading segment. Indicates the start time of the next loading segment.

[0132] Calculate the displacement inheritance offset: ;in, This indicates the residual displacement inherited offset.

[0133] The residual displacement is used to inherit the offset and then the actual displacement segment of the next loading segment is translated and corrected. ;in, This represents the actual displacement segment of the next loading segment after the residual displacement is inherited.

[0134] Through the above processing, the next hysteresis cycle can continue loading from the residual displacement position left over from the previous cycle, thereby maintaining the physical continuity of the displacement evolution process of the seismic isolation support.

[0135] S35. Generate the reconstructed displacement trajectory results of the seismic isolation bearing.

[0136] Specifically, all the actual displacement segments of the supports after completing the hysteresis endpoint locking and residual displacement inheritance processing are sequentially spliced ​​together according to the time order in the support hysteresis response segment chain:

[0137] ;

[0138] in, This indicates the reconstruction results of the displacement trajectory of the seismic isolation bearing. Indicates the first A segment of the actual displacement of the support after continuous correction. This indicates the total number of support hysteresis response segments. It indicates that the pieces are joined together sequentially from beginning to end in chronological order.

[0139] The final result obtained is a reconstruction of the displacement trajectory of the seismic isolation bearing that is continuous on the time axis, has no displacement jumps, and can completely reflect the multiple loading, unloading and residual dwell processes of the seismic isolation bearing.

[0140] This invention also provides a seismic isolation bearing displacement trajectory reconstruction system based on time-frequency domain joint inversion, comprising:

[0141] The bearing hysteresis response segment chain generation module is used to take the original vibration response of the seismic isolation bearing as input and, based on the energy migration order of the loading segment, unloading segment and residual dwell segment of the seismic isolation bearing in reciprocating shear motion, invert the original vibration response into a bearing hysteresis response segment chain. The bearing hysteresis response segment chain is used to characterize the continuous motion process of the seismic isolation bearing from the loading segment, unloading segment to the residual dwell segment.

[0142] The support real displacement segment generation module is used to take the support hysteresis response segment chain as input, perform time-domain integral inversion and frequency-domain phase inversion in each support hysteresis response segment, and identify and remove the pseudo displacement components formed by sensor zero drift according to the reverse connection relationship between the loading segment and the unloading segment, retain the real displacement components formed by the shear deformation of the seismic isolation support, and generate the support real displacement segment.

[0143] The seismic isolation bearing displacement trajectory reconstruction module is used to take the actual displacement segment of the bearing as input, and according to the successive order of the hysteresis response segment chain of the bearing, perform hysteresis endpoint locking and residual displacement inheritance on the actual displacement segments of adjacent bearings to obtain the seismic isolation bearing displacement trajectory reconstruction result.

[0144] This invention also provides an electronic device, comprising:

[0145] Memory, on which computer programs / instructions are stored;

[0146] When the processor executes the computer program / instruction, it implements the steps of the seismic isolation bearing displacement trajectory reconstruction method based on time-frequency domain joint inversion as described above.

[0147] This invention also provides a computer-readable storage medium storing a computer program / instructions thereon, which, when executed by a processor, implements the steps of the seismic isolation support displacement trajectory reconstruction method based on time-frequency domain joint inversion as described above.

[0148] It is understood that the same or similar parts in the above embodiments can be referred to each other, and the contents not described in detail in some embodiments can be referred to the same or similar contents in other embodiments.

[0149] It should be understood that the functional unit modules in the various embodiments of the present invention can be concentrated in one processing unit, or each unit module can exist physically separately, or two or more unit modules can be integrated into one unit module, and can be implemented in hardware or software.

[0150] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0151] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0152] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0153] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1The steps of the function specified in one or more boxes.

[0154] Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention. Those skilled in the art can make changes, modifications, substitutions and variations to the above embodiments within the scope of the present invention.

Claims

1. A method for reconstructing the displacement trajectory of seismic isolation bearings based on joint time-frequency domain inversion, characterized in that, Includes the following steps: S1. Taking the original vibration response of the seismic isolation bearing as input, and based on the energy migration order of the loading segment, unloading segment and residual dwell segment of the seismic isolation bearing in the reciprocating shear motion, the original vibration response is inverted into a hysteresis response segment chain of the bearing. The hysteresis response segment chain of the bearing is used to characterize the continuous motion process of the seismic isolation bearing from the loading segment, unloading segment to the residual dwell segment. S2. Using the support hysteresis response segment chain as input, perform time-domain integral inversion and frequency-domain phase inversion in each support hysteresis response segment, and identify and remove the pseudo displacement components formed by sensor zero drift according to the reverse connection relationship between the loading segment and the unloading segment, retain the real displacement components formed by the shear deformation of the seismic isolation support, and generate the real displacement segment of the support. S3. Using the actual displacement segment of the support as input, according to the sequence of the support hysteresis response segment chain, perform hysteresis endpoint locking and residual displacement inheritance on the actual displacement segments of adjacent supports to obtain the displacement trajectory reconstruction result of the seismic isolation support.

2. The method for reconstructing the displacement trajectory of seismic isolation bearings based on joint time-frequency domain inversion according to claim 1, characterized in that, In step S1, the original vibration response is used as input, and a short-time energy sequence is calculated based on the original vibration response to obtain the energy migration envelope. On the energy migration envelope, the energy continuously rising edge, the energy continuously falling edge, and the low-energy segment where the energy is stably below a preset threshold are detected, which correspond to the loading segment, unloading segment, and residual dwell segment of the seismic isolation bearing in the reciprocating shear motion, respectively.

3. The method for reconstructing the displacement trajectory of seismic isolation bearings based on joint time-frequency domain inversion according to claim 1, characterized in that, In step S1, the loading segment, unloading segment, and residual dwell segment are sequentially spliced ​​together according to the starting time of the loading segment, unloading segment, and residual dwell segment to form a support hysteresis response segment chain consisting of alternating loading segment, unloading segment, and residual dwell segment. The temporal connection of the support hysteresis response segment chain directly maps to the entire continuous motion process of the seismic isolation support from force slippage and directional rotation to residual dwell.

4. The method for reconstructing the displacement trajectory of seismic isolation bearings based on joint time-frequency domain inversion according to claim 1, characterized in that, In step S2, the current support hysteresis response segment is extracted from the support hysteresis response segment chain, the acceleration signal within the current support hysteresis response segment is obtained, and the acceleration signal within the current support hysteresis response segment is subjected to time-domain integral inversion and frequency-domain phase inversion to obtain the initial displacement trajectory including low-frequency pseudo-displacement components and real displacement components.

5. The method for reconstructing the displacement trajectory of seismic isolation bearings based on joint time-frequency domain inversion according to claim 4, characterized in that, In step S2, based on the current support hysteresis response segment, the loading segment and unloading segment that belong to the same complete hysteresis cycle as the current support hysteresis response segment in the support hysteresis response segment chain are extracted, and the reverse connection relationship between the loading segment and the unloading segment is determined. Using the reverse connection relationship, the low-frequency displacement components in the initial displacement trajectories of the loading and unloading segments are compared for directional consistency: if a certain low-frequency displacement component maintains the same cumulative trend in the loading and unloading segments, it is determined to be a low-frequency pseudo-displacement component formed by sensor zero drift and is removed. All the retained true displacement components are reconstructed to generate the true displacement segment of the support corresponding to the hysteresis response segment of the support.

6. The method for reconstructing the displacement trajectory of seismic isolation bearings based on joint time-frequency domain inversion according to claim 1, characterized in that, In step S3, the support hysteresis response segment chain includes a loading segment, an unloading segment, and a residual dwell segment arranged in time sequence. Based on the time sequence of the loading segment, unloading segment, and residual dwell segment in the support hysteresis response segment chain, the connection type between each support real displacement segment and its preceding and following segments is determined. The connection type includes loading segment to unloading segment, unloading segment to residual dwell segment, and residual dwell segment to the next loading segment. Based on the connection type of adjacent segments in the support hysteresis response segment chain, the displacement state formed at the end of the previous support real displacement segment is transferred to the next support real displacement segment. Hysteresis endpoint locking and residual displacement inheritance are performed. All support real displacement segments processed by hysteresis endpoint locking and residual displacement inheritance are spliced ​​together end to end according to the connection order of the support hysteresis response segment chain to obtain a reconstruction result of the seismic isolation support displacement trajectory that is continuous on the time axis and has no displacement jump.

7. The method for reconstructing the displacement trajectory of seismic isolation bearings based on joint time-frequency domain inversion according to claim 6, characterized in that, The execution of hysteresis endpoint locking and residual displacement inheritance specifically includes: For adjacent segments where the loading segment connects to the unloading segment, the actual displacement of the end support of the loading segment is set as the extreme displacement reference at the beginning of the unloading segment, and hysteresis endpoint locking is performed: the starting displacement of the actual displacement segment of the unloading segment is corrected with this extreme displacement reference, so that the displacement trajectory of the unloading segment continuously rebounds in the opposite direction from the maximum displacement point at the end of the loading segment, eliminating the endpoint misalignment caused by segmentation. For adjacent segments of the unloading segment and the residual dwell segment, the actual displacement of the end support of the unloading segment is identified as the residual displacement of this hysteresis cycle, and it is transferred to the actual displacement segment of the residual dwell segment as a constant displacement reference, and residual displacement inheritance is performed. For the adjacent segment of the residual dwell segment following the next loading segment, the inherited residual displacement is used as the starting displacement of the actual displacement segment of the next loading segment, so as to realize the continuous transfer of residual displacement between adjacent hysteresis cycles.

8. A seismic isolation bearing displacement trajectory reconstruction system based on joint time-frequency domain inversion, characterized in that, include: The bearing hysteresis response segment chain generation module is used to take the original vibration response of the seismic isolation bearing as input and, based on the energy migration order of the loading segment, unloading segment and residual dwell segment of the seismic isolation bearing in reciprocating shear motion, invert the original vibration response into a bearing hysteresis response segment chain. The bearing hysteresis response segment chain is used to characterize the continuous motion process of the seismic isolation bearing from the loading segment, unloading segment to the residual dwell segment. The support real displacement segment generation module is used to take the support hysteresis response segment chain as input, perform time-domain integral inversion and frequency-domain phase inversion in each support hysteresis response segment, and identify and remove the pseudo displacement components formed by sensor zero drift according to the reverse connection relationship between the loading segment and the unloading segment, retain the real displacement components formed by the shear deformation of the seismic isolation support, and generate the support real displacement segment. The seismic isolation bearing displacement trajectory reconstruction module is used to take the actual displacement segment of the bearing as input, and according to the successive order of the hysteresis response segment chain of the bearing, perform hysteresis endpoint locking and residual displacement inheritance on the actual displacement segments of adjacent bearings to obtain the seismic isolation bearing displacement trajectory reconstruction result.

9. An electronic device, characterized in that, include: Memory, on which computer programs / instructions are stored; When the processor executes the computer program / instruction, it implements the steps of the seismic isolation bearing displacement trajectory reconstruction method based on time-frequency domain joint inversion as described in any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program / instructions stored thereon, characterized in that, When the computer program / instruction is executed by the processor, it implements the steps of the seismic isolation bearing displacement trajectory reconstruction method based on time-frequency domain joint inversion as described in any one of claims 1 to 7.

Images