Method for determining main control parameters of seismic signal reproduction, computer equipment and medium

Abstract

The present invention discloses a method for determining a main control parameter for reproducing a seismic signal, a computer device and a medium, and relates to the field of seismic signal reproduction. The method comprises obtaining a response signal of a target engineering structure model system, determining a characteristic period of the target engineering structure model system according to the response signal, calculating a seismic triple response spectrum of a preset seismic GM signal according to parameter information of the preset seismic GM signal, determining a sensitive zone boundary period value of the target engineering structure model system according to the seismic triple response spectrum, calculating a main control parameter discrimination index according to the characteristic period and the sensitive zone boundary period value, and determining a main control parameter used when reproducing the preset seismic GM signal according to the main control parameter discrimination index and the type of the target engineering structure model system. The present invention adds velocity and displacement parameters that are more sensitive to the low-frequency components of the signal to the main control parameters, and can improve the reproduction accuracy of medium and low-frequency seismic signals.

Description

Technical Field

[0001] The present invention relates to the technical field of seismic signal reproduction, and in particular to a method for determining a main control parameter of seismic signal reproduction, a computer device and a medium. Background Art

[0002] Earthquakes are one of the most destructive natural disasters in the world, often causing serious losses including casualties and damage to civil infrastructure. The shaking table test is a test method that can truly reproduce structural responses and ground motions in the laboratory, and can reproduce the entire process of earthquake action or artificial action in real time. The shaking table system is an important way to obtain the seismic response of engineering structures and evaluate the seismic performance of engineering structures. At present, hundreds of earthquake simulation shaking table systems have been built around the world. Accurately reproducing the time series of ground motions with different frequency characteristics is one of the core issues of the shaking table test. Among the various shaking table systems mentioned above, the electro-hydraulic shaking table that generates a strong driving force through a hydraulic actuator is currently the most commonly used. In order to improve the accuracy of the seismic signal emergence, many corresponding control algorithms and technologies have been continuously developed and applied to the electro-hydraulic shaking table system. Among them, the development from proportional integral differential (PID) control that only uses displacement signal feedback control to three-parameter control (TVC) composed of displacement, velocity and acceleration signals is one of the more critical technological innovations. It improves the damping ratio and natural frequency of the hydraulic system, thereby using displacement, velocity and acceleration parameters to control the low-frequency, medium-frequency and high-frequency bands respectively, expanding the system bandwidth.

[0003] In the process of seismic motion reproduction based on the TVC framework, the acceleration parameter is generally selected as the main control parameter, and the accuracy of the acceleration signal reproduction is used as the basis for judging the accuracy of seismic motion reproduction, while the velocity and displacement are used as auxiliary control parameters. The acceleration reference signal generates three signals of acceleration (A), velocity (V), and displacement (D) after passing through the three-parameter signal generator and enters the closed-loop control system. The inner-loop three-parameter servo control introduces acceleration and velocity feedback on the basis of traditional displacement feedback control, thereby using displacement, velocity, and acceleration parameters to control low-frequency, medium-frequency, and high-frequency bands respectively. Among them, the acceleration and displacement feedback signals are measured in real time by the acceleration sensor and the displacement sensor in the experiment, and the velocity signal is synthesized by the two in real time. The outer-loop adaptive control technology includes amplitude phase control, adaptive harmonic elimination system, adaptive reversal control, online iteration system, etc. In each frame of the earthquake simulation test, the system acceleration impedance is estimated by the system drive signal and the table control acceleration signal, and then the system drive signal is continuously updated, so that the table control acceleration signal reproduction accuracy remains at a high level. This master parameter selection and reproduction control strategy limited to acceleration has achieved good results in most shaking table tests of engineering structures. However, with the continuous deepening of engineering seismic research and the diversification of test structure types, the adaptability of this method is relatively low. The velocity pulse effect in the near-fault zone of strong earthquakes and the low-frequency amplification effect of the far-field soft site make the low-frequency components in the ground motion signal prominent (as low as 0.1Hz in extreme cases). However, some high-rise structural systems and soft soil-structure interaction systems with significant strong earthquake softening characteristics have a long characteristic period (possibly higher than 1.0s), which makes them insensitive to acceleration signals. Using acceleration as the master control parameter may lead to poor reproduction accuracy of acceleration signals in the low-frequency band (especially in the absence of iterative calculation). On the other hand, due to the performance limitations of the mechanical hydraulic system of the shaking table, in order to prevent the table displacement from exceeding the limit and causing the table and the hydraulic cylinder buffer layer to collide, high-pass filtering is required in the process of the three-parameter signal generator generating velocity and displacement parameters by integrating the acceleration master control parameter. Although the magnitude of the acceleration reference signal in the low-frequency band is very small, for most test seismic motions that have been suppressed by the time domain signal based on the test similarity ratio, this high-pass filtering will not have a significant effect on the signal reproduction accuracy; however, for seismic motion time histories with a wider frequency band, especially near-fault seismic motions with rich low-frequency components and velocity pulse effects, since the frequency of some signal components is too low, the velocity time series and displacement time series reproduced by the shaking table system will have a large deviation, resulting in low accuracy of the test results. Summary of the invention

[0004] The purpose of the present invention is to provide a method for determining main control parameters of earthquake signal reproduction, which can improve the accuracy of earthquake signal reproduction of a medium- and long-period target engineering structure model system under earthquake motion with rich low-frequency components.

[0005] To achieve the above object, the present invention provides the following solutions:

[0006] A method for determining a main control parameter for reproducing a seismic signal, the method comprising:

[0007] S1: Acquire a response signal of a target engineering structure model system, and determine a characteristic period of the target engineering structure model system according to the response signal; the response signal is a signal generated by the target engineering structure model system after a preset acceleration peak value and white noise of a preset frequency range are input into the target engineering structure model system through a vibration table; the target engineering structure model system is located on the table of the vibration table;

[0008] S2: calculating the earthquake triple response spectrum of the preset earthquake GM signal according to the parameter information of the preset earthquake GM signal, and determining the boundary period value of the sensitive area of ​​the target engineering structure model system according to the earthquake triple response spectrum; the parameter information includes displacement time history, velocity time history and acceleration time history; the earthquake triple response spectrum includes displacement response spectrum, pseudo velocity response spectrum and pseudo acceleration response spectrum; the boundary period value of the sensitive area includes the boundary period value between the acceleration sensitive area and the velocity sensitive area and the boundary period value between the velocity sensitive area and the displacement sensitive area;

[0009] S3: Calculating the main control parameter discrimination index according to the characteristic period and the sensitive area boundary period value;

[0010] S4: Determine the main control parameter used to reproduce the preset seismic GM signal according to the main control parameter discrimination index and the type of the target engineering structure model system; the main control parameter includes displacement, velocity and acceleration; the type of the target engineering structure model system includes a solid material model and a granular material model.

[0011] Optionally, S1 specifically includes:

[0012] Acquire a response signal of the target engineering structure model system; the response signal includes a first acceleration response time history and a second acceleration response time history; the first acceleration response time history is acquired by a first acceleration sensor; the second acceleration response time history is acquired by a second acceleration sensor; the first acceleration sensor is located at the top of the target engineering structure model system; the second acceleration sensor is located on the table of the vibration table;

[0013] Performing Fourier transformation on the first acceleration response time history and the second acceleration response time history respectively to obtain a first Fourier amplitude function and a second Fourier amplitude function;

[0014] The characteristic period of the target engineering structure model system is determined according to the first Fourier amplitude function and the second Fourier amplitude function.

[0015] Optionally, determining the characteristic period of the target engineering structure model system according to the first Fourier amplitude function and the second Fourier amplitude function specifically includes:

[0016] Calculating an amplitude-frequency characteristic function according to the first Fourier amplitude function and the second Fourier amplitude function;

[0017] The characteristic period of the target engineering structure model system is determined according to the frequency corresponding to the maximum value of the amplitude-frequency characteristic function.

[0018] Optionally, the calculating of the earthquake triple response spectrum of the preset earthquake GM signal according to the parameter information of the preset earthquake GM signal in S2 specifically includes:

[0019] Taking the parameter information of the preset seismic GM signal as input, the Duhamel integration algorithm is adopted to calculate the seismic triple response spectrum of the preset seismic GM signal.

[0020] Optionally, the step of determining the sensitive zone boundary period value of the target engineering structure model system according to the earthquake triple response spectrum in S2 specifically includes:

[0021] Determine the peak acceleration, peak velocity and peak displacement of the preset ground motion GM signal according to the ground motion triple response spectrum;

[0022] The sensitive area boundary period value of the target engineering structure model system is determined according to the peak acceleration, the peak velocity and the peak displacement.

[0023] Optionally, determining the sensitive area limit period value of the target engineering structure model system according to the peak acceleration, the peak velocity and the peak displacement specifically includes:

[0024] Multiplying the peak acceleration, the peak velocity, and the peak displacement by their respective preset magnification factors to obtain amplified peak acceleration, amplified peak velocity, and amplified peak displacement;

[0025] Based on the amplified peak acceleration, an amplified peak acceleration straight line is drawn in the triple earthquake response spectrum; the amplified peak acceleration straight line is perpendicular to the y-axis of the acceleration coordinate system where the pseudo acceleration response spectrum in the triple earthquake response spectrum is located;

[0026] Drawing an amplified peak velocity straight line in the triple earthquake response spectrum based on the amplified peak velocity; the amplified peak velocity straight line is perpendicular to the y-axis of the velocity coordinate system where the pseudo velocity response spectrum in the triple earthquake response spectrum is located;

[0027] An amplified peak displacement straight line is drawn in the triple earthquake response spectrum based on the amplified peak displacement; the amplified peak displacement straight line is perpendicular to the y-axis of the displacement coordinate system where the displacement response spectrum in the triple earthquake response spectrum is located; the x-axis of the velocity coordinate system is in the horizontal direction, the y-axis of the velocity coordinate system is in the vertical direction, the y-axis of the acceleration coordinate system and the y-axis of the displacement coordinate system are both at an angle of 45° to the x-axis of the velocity coordinate system, and the y-axis of the acceleration coordinate system and the y-axis of the displacement coordinate system are perpendicular to each other;

[0028] The boundary period value between the acceleration sensitive area and the velocity sensitive area is determined according to the amplified peak acceleration straight line and the amplified peak velocity straight line, and the boundary period value between the velocity sensitive area and the displacement sensitive area is determined according to the amplified peak velocity straight line and the amplified peak displacement straight line; the boundary period value between the acceleration sensitive area and the velocity sensitive area and the boundary period value between the velocity sensitive area and the displacement sensitive area constitute the boundary period value of the sensitive area.

[0029] Optionally, the main control parameter discrimination index includes a first main control parameter discrimination index and a second main control parameter discrimination index; the expression of the first main control parameter discrimination index is:

[0030]

[0031] Among them, Z1 is the first main control parameter identification index, T m is the characteristic period of the target engineering structure model system, T c It is the boundary period value between the acceleration sensitive area and the speed sensitive area;

[0032] The expression of the second main control parameter discrimination index is:

[0033]

[0034] Among them, Z2 is the second main control parameter judgment index, T d It is the boundary period value between the velocity sensitive area and the displacement sensitive area.

[0035] Optionally, S4 specifically includes:

[0036] When the type of the target engineering structure model system is a granular material model, the main control parameter used to reproduce the preset seismic GM signal is determined according to the main control parameter discrimination index and the peak acceleration of the preset seismic GM signal;

[0037] When the type of the target engineering structure model system is a solid material model, the main control parameter used to reproduce the preset ground motion GM signal is determined according to the main control parameter discrimination index.

[0038] The present invention also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the main control parameter determination method for reproducing the above-mentioned seismic signal.

[0039] The present invention also provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and the computer program is suitable for being loaded by a processor and executing the main control parameter determination method for reproducing the above-mentioned seismic signal.

[0040] According to the specific embodiment provided by the present invention, the present invention discloses the following technical effects: the present invention provides a method for determining the main control parameters for reproducing earthquake motion signals, a computer device and a medium, the method first obtains the response signal of the target engineering structure model system, determines the characteristic period of the target engineering structure model system according to the response signal, calculates the earthquake triple response spectrum of the preset earthquake motion GM signal according to the parameter information of the preset earthquake motion GM signal, determines the sensitive zone boundary period value of the target engineering structure model system according to the earthquake triple response spectrum, wherein the sensitive zone boundary period value includes the boundary period value between the acceleration sensitive zone and the velocity sensitive zone and the boundary period value between the velocity sensitive zone and the displacement sensitive zone, calculates the main control parameter discrimination index according to the characteristic period and the sensitive zone boundary period value, and finally determines the main control parameter used when reproducing the preset earthquake motion GM signal according to the main control parameter discrimination index and the type of the target engineering structure model system, wherein the main control parameter includes displacement, velocity and acceleration. The present invention adds velocity and displacement parameters that are more sensitive to the low-frequency components of the signal to the selection of the main control parameters, which can improve the reproduction accuracy of medium and low-frequency earthquake motion signals. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0042] Figure 1 A schematic flow chart of a method for determining main control parameters for earthquake signal reproduction provided by the present invention;

[0043] Figure 2A diagram showing the main steps of the method for determining the main control parameters for reproducing the earthquake signal provided by the present invention;

[0044] Figure 3 A schematic diagram of a soft soil-underground structure interaction system and sensor arrangement provided by the present invention;

[0045] Figure 4 The triple response spectrum of earthquake motion for determining the boundary period of each sensitive area provided by the present invention;

[0046] Figure 5 A schematic diagram of the dynamic shear strain of soil under strong earthquake provided by the present invention;

[0047] Figure 6 It is a structural schematic diagram of a computer device provided by this application.

[0048] Explanation of symbols:

[0049] 1000 - computer device; 1001 - processor; 1002 - communication bus; 1003 - user interface; 1004 - network interface; 1005 - memory. DETAILED DESCRIPTION

[0050] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0051] The purpose of the present invention is to provide a method for determining the main control parameters for reproducing seismic motion signals, a computer device and a medium. By adding speed and displacement parameters that are more sensitive to the low-frequency components of the signal to the selection of the main control parameters, the forced high-pass filtering required in the process of generating speed and displacement parameters by integrating the acceleration main control parameters of a three-parameter signal generator is avoided, thereby improving the reproduction accuracy of medium and low-frequency seismic motion signals.

[0052] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0053] like Figure 1 and Figure 2 As shown, the present invention provides a method for determining a main control parameter for earthquake signal reproduction, the method comprising:

[0054] S1: Acquire a response signal of a target engineering structure model system, and determine a characteristic period of the target engineering structure model system according to the response signal; the response signal is a signal generated by the target engineering structure model system after a preset acceleration peak value and white noise of a preset frequency range are input into the target engineering structure model system through a vibration table; the target engineering structure model system is located on the tabletop of the vibration table.

[0055] S2: Calculate the earthquake triple response spectrum of the preset earthquake GM signal according to the parameter information of the preset earthquake GM signal, and determine the boundary period value of the sensitive area of ​​the target engineering structure model system according to the earthquake triple response spectrum; the parameter information includes displacement time history, velocity time history and acceleration time history; the earthquake triple response spectrum includes displacement response spectrum, pseudo velocity response spectrum and pseudo acceleration response spectrum; the boundary period value of the sensitive area includes the boundary period value between the acceleration sensitive area and the velocity sensitive area and the boundary period value between the velocity sensitive area and the displacement sensitive area.

[0056] S3: Calculate the main control parameter discrimination index according to the characteristic period and the sensitive area boundary period value.

[0057] S4: Determine the main control parameter used to reproduce the preset seismic GM signal according to the main control parameter discrimination index and the type of the target engineering structure model system; the main control parameter includes displacement, velocity and acceleration; the type of the target engineering structure model system includes a solid material model and a granular material model.

[0058] Among them, S1 specifically includes:

[0059] S11: Obtain a response signal of the target engineering structure model system; the response signal includes a first acceleration response time series and a second acceleration response time series; the first acceleration response time series is acquired by a first acceleration sensor; the second acceleration response time series is acquired by a second acceleration sensor; the first acceleration sensor is located at the top of the target engineering structure model system; the second acceleration sensor is located on the table top of the vibration table.

[0060] S12: Performing Fourier transformation on the first acceleration response time history and the second acceleration response time history respectively to obtain a first Fourier amplitude function and a second Fourier amplitude function.

[0061] S13: determining a characteristic period of the target engineering structure model system according to the first Fourier amplitude function and the second Fourier amplitude function.

[0062] Determining the characteristic period of the target engineering structure model system according to the first Fourier amplitude function and the second Fourier amplitude function in S13 specifically includes:

[0063] S131: Calculate an amplitude-frequency characteristic function according to the first Fourier amplitude function and the second Fourier amplitude function.

[0064] S132: Determine a characteristic period of the target engineering structure model system according to a frequency corresponding to a maximum value of the amplitude-frequency characteristic function.

[0065] Specifically, Figure 3 As shown, a given target engineering structure model system to be subjected to a shaking table test is anchored to the surface of the shaking table, and a first acceleration sensor A1 and a second acceleration sensor A2 are installed on the target engineering structure model system and the shaking table surface, respectively, wherein the target engineering structure model system includes a soil box and an underground structure. A broadband white noise with a preset acceleration peak value of 0.05g and a preset frequency range of 0 to 50Hz is input into the target engineering structure model system through the shaking table, and the response signals of the target engineering structure model system collected by the acceleration sensors A1 and A2 are recorded respectively, and the response signals include a first acceleration response time series WNa1(t) and a second acceleration response time series WNa2(t). The preset acceleration peak value of the white noise is controlled to be 0.05g in order to avoid damage to the target engineering structure model system. The first acceleration response time series WNa1(t) and the second acceleration response time series WNa2(t) are Fourier transformed respectively to obtain a first Fourier amplitude function F a1 (f) and the second Fourier amplitude function F a2 (f). Then calculate the amplitude-frequency characteristic function A of the target engineering structure model system m (f) = F a2 (f) / F a1 (f), according to the amplitude-frequency characteristic function A m (f) The frequency corresponding to the maximum value f max Determine the characteristic period T of the model system m , whose expression is T m =1 / f max .

[0066] The step S2 of calculating the earthquake triple response spectrum of the preset earthquake GM signal according to the parameter information of the preset earthquake GM signal specifically includes:

[0067] S21: taking the parameter information of the preset seismic GM signal as input, and using the Duhamel integration algorithm to calculate the seismic triple response spectrum of the preset seismic GM signal.

[0068] Optionally, the step of determining the sensitive zone boundary period value of the target engineering structure model system according to the earthquake triple response spectrum in S2 specifically includes:

[0069] S22: determining the peak acceleration, peak velocity and peak displacement of the preset ground motion GM signal according to the ground motion triple response spectrum;

[0070] S23: Determine a sensitive area boundary period value of the target engineering structure model system according to the peak acceleration, the peak velocity and the peak displacement.

[0071] Wherein, determining the sensitive area limit period value of the target engineering structure model system according to the peak acceleration, the peak velocity and the peak displacement in step S23 specifically includes:

[0072] S231: Multiplying the peak acceleration, the peak velocity, and the peak displacement by their respective preset magnification factors to obtain amplified peak acceleration, amplified peak velocity, and amplified peak displacement.

[0073] S232: drawing an amplified peak acceleration straight line in the earthquake triple response spectrum based on the amplified peak acceleration; the amplified peak acceleration straight line is perpendicular to the y-axis of the acceleration coordinate system where the pseudo acceleration response spectrum in the earthquake triple response spectrum is located.

[0074] S233: drawing an amplified peak velocity straight line in the earthquake triple response spectrum based on the amplified peak velocity; the amplified peak velocity straight line is perpendicular to the y-axis of the velocity coordinate system where the pseudo velocity response spectrum in the earthquake triple response spectrum is located.

[0075] S234: drawing an amplified peak displacement straight line in the earthquake triple response spectrum based on the amplified peak displacement; the amplified peak displacement straight line is perpendicular to the y-axis of the displacement coordinate system in which the displacement response spectrum in the earthquake triple response spectrum is located; the x-axis of the velocity coordinate system is horizontal, the y-axis of the velocity coordinate system is vertical, the y-axis of the acceleration coordinate system and the y-axis of the displacement coordinate system are both at an angle of 45° with the x-axis of the velocity coordinate system, and the y-axis of the acceleration coordinate system and the y-axis of the displacement coordinate system are perpendicular to each other.

[0076] S235: Determine the boundary period value between the acceleration sensitive area and the velocity sensitive area according to the amplified peak acceleration straight line and the amplified peak velocity straight line, and determine the boundary period value between the velocity sensitive area and the displacement sensitive area according to the amplified peak velocity straight line and the amplified peak displacement straight line; the boundary period value between the acceleration sensitive area and the velocity sensitive area and the boundary period value between the velocity sensitive area and the displacement sensitive area constitute the boundary period value of the sensitive area.

[0077] For a given preset ground motion GM signal, its acceleration time history, velocity time history, and displacement time history are recorded as xg(t), using Duhamel integration, based on xg(t) calculates the pseudo acceleration response spectrum of the preset ground motion GM signal Pseudo velocity response spectrum and displacement response spectrum x(t). The pseudo-acceleration response spectrum, pseudo-velocity response spectrum and displacement response spectrum are respectively represented by S a (T), S v (T) and S d (T) indicates that the expression is as follows:

[0078]

[0079]

[0080]

[0081] Among them, S d is the maximum value in the displacement response spectrum x(t), S d It is the maximum displacement response value of a single-degree-of-freedom system when an earthquake action is applied to a single-degree-of-freedom structure. ω is the natural circular frequency of the system under the undamped state. ω' is the actual circular frequency of vibration under the current damping ζ. ζ is the damping ratio, t represents time, τ represents any moment under the action of earthquake, S v is the pseudo velocity response spectrum The maximum value in S a is the pseudo acceleration response spectrum In this embodiment, ζ=0.05.

[0082] Then the pseudo velocity response spectrum S is plotted in the double logarithmic coordinate system v (T), and because In S v (T) Add S to the graph a (T) and S d (T) axis, becomes a four-logarithmic coordinate system, and the triple response spectrum can be obtained. Figure 4As shown, the horizontal direction is used as the velocity coordinate system T n The vertical direction is the y-axis of the velocity coordinate system. The y-axis of the acceleration coordinate system and the displacement coordinate system are both perpendicular to the T axis of the velocity coordinate system. n The axes form an angle of 45°, and the y axes of the acceleration coordinate system and the displacement coordinate system are perpendicular to each other. The pseudo acceleration response spectrum S is plotted in the above coordinate system. a (T), pseudo velocity response spectrum S v (T) and displacement response spectrum S d (T), thereby the earthquake triple response spectrum of the preset earthquake GM signal can be obtained.

[0083] Then calculate the peak acceleration of the preset ground motion GM signal Peak speed and peak displacement u g0 , which is calculated based on the acceleration time history, velocity time history and displacement time history of the preset ground motion GM signal, that is, u g0 =max|x g (t)|, and in the triple response spectrum of earthquake motion based on u g0 Draw three thin dotted lines respectively, such as Figure 4 As shown, these three thin dotted lines are perpendicular to S v (T), S a (T) and S d (T) is the y-axis of the coordinate system. The curve in the acceleration sensitive area is the pseudo acceleration response spectrum The curve in the velocity sensitive area is the pseudo velocity response spectrum The curve in the displacement sensitive area is the displacement response spectrum x(t), u g0 is the peak displacement of the preset seismic GM signal. Multiply by the magnification factor 2.12 to get the magnified Will Multiply by the magnification factor 1.65 to get the magnified will u g0 Multiply by the magnification factor 1.39 to get the magnified u g0 Then based on the above enlarged After enlargement and the enlarged u g0 Draw the amplified peak acceleration straight line, the amplified peak velocity straight line and the amplified peak displacement straight line, which are respectively Figure 4 The thick dashed lines b-c, c-d and d-e are perpendicular to the y-axes of their respective coordinate systems, c and d are the intersections of the three thick dashed lines in the triple reaction spectrum, and then Figure 4 T corresponding to point cc As the boundary period value between the acceleration sensitive area and the velocity sensitive area, T corresponding to point d d As the boundary period value between the velocity sensitive area and the displacement sensitive area, T c With T d The T corresponding to point c and point d are n Axis coordinates.

[0084] In getting T c With T d After that, it can be used to calculate the main control parameter discrimination index. The main control parameter discrimination index in step S3 includes the first main control parameter discrimination index and the second main control parameter discrimination index; wherein, the expression of the first main control parameter discrimination index is:

[0085]

[0086] Among them, Z1 is the first main control parameter identification index, T m is the characteristic period of the target engineering structure model system, T c It is the boundary period value between the acceleration sensitive area and the speed sensitive area;

[0087] The expression of the second main control parameter judgment index is:

[0088]

[0089] Among them, Z2 is the second main control parameter judgment index, T d It is the boundary period value between the velocity sensitive area and the displacement sensitive area.

[0090] Step S4 specifically includes:

[0091] When the type of the target engineering structure model system is a granular material model, the main control parameter used to reproduce the preset seismic GM signal is determined according to the main control parameter discrimination index and the peak acceleration of the preset seismic GM signal. When the type of the target engineering structure model system is a solid material model, the main control parameter used to reproduce the preset seismic GM signal is determined according to the main control parameter discrimination index.

[0092] Specifically, for the solid material model system mainly composed of reinforced concrete, steel, organic glass and other solid materials, the influence of system stiffness degradation under strong earthquake is not considered, and the main control parameter selection method is shown in Table 1.

[0093] Table 1 Selection method of main control parameters for seismic signal reproduction in shaking table test of solid material model system

[0094] Main control parameter judgment index value range Main control parameters for seismic signal reproduction <![CDATA[Z1<1.0]]> Acceleration <![CDATA[Z1≥1.0 and Z2<1.0]]> speed <![CDATA[Z2≥1.0]]> Displacement

[0095] When Z1<1.0, the target engineering structure model system is located in the acceleration sensitive area of ​​the GM earthquake response spectrum, and the acceleration parameter is selected as the main control parameter; when Z1≥1.0 and Z2<1.0, the target engineering structure model system is located in the velocity sensitive area of ​​the GM earthquake response spectrum, and the velocity parameter is selected as the main control parameter; when Z2≥1.0, the target engineering structure model system is located in the displacement sensitive area of ​​the GM earthquake response spectrum, and the displacement parameter is selected as the main control parameter.

[0096] like Figure 5 As shown in the figure, under strong earthquake, the dynamic shear strain γ (x-axis) of the soil increases, the dynamic shear modulus G decreases (left y-axis), the stiffness of the target engineering structure model system decreases, and the instantaneous characteristic period of the target engineering structure model system is extended, so that it will change from the acceleration sensitive area to the velocity sensitive area, and even the displacement sensitive area. Therefore, for the granular material model system mainly composed of granular materials such as clay and sand, it is necessary to consider the impact of the system stiffness degradation under strong earthquake. Therefore, this embodiment uses the peak acceleration The distinction is made at 0.2g, and the method for selecting the main control parameters is shown in Table 2.

[0097]

[0098]

[0099] Table 2 Selection method of main control parameters for seismic signal reproduction in shaking table test of granular material model system

[0100] When Z1<0.8, the target engineering structure model system is located in the acceleration sensitive area of ​​the GM response spectrum of the earthquake motion, and the acceleration parameter is selected as the main control parameter; when 0.8≤Z1<1.0, the instantaneous characteristic period of the target engineering structure model system under the action of strong earthquake may move from the acceleration sensitive area of ​​the GM response spectrum of the earthquake motion to the velocity sensitive area, so the weak earthquake condition Select acceleration parameter as the main control parameter, strong earthquake condition The velocity parameter is selected as the main control parameter; when Z1≥1.0 and Z2<0.8, the target engineering structure model system is located in the velocity sensitive area of ​​the seismic GM response spectrum, and the velocity parameter is selected as the main control parameter; when 0.8≤Z2<1.0, the instantaneous characteristic period of the target engineering structure model system under strong earthquake may move from the velocity sensitive area of ​​the seismic GM response spectrum to the displacement sensitive area, so the weak earthquake condition Select the speed parameter as the main control parameter, strong earthquake condition The displacement parameter is selected as the main control parameter; when Z2≥1.0, the target engineering structure model system is located in the displacement sensitive area of ​​the seismic GM response spectrum, and the displacement parameter is selected as the main control parameter.

[0101] Based on the above content, this embodiment was carried out on MTS in the National Key Laboratory of Disaster Prevention in Civil Engineering of Tongji University. Figure 3 The shaking table test of a certain soft soil-underground structure interaction system shown is used to verify the implementation of the present invention.

[0102] The vibration table has a surface size of 4m×6m and a maximum load capacity of 70t. The vibration table can input seismic motion in two horizontal directions and has three degrees of freedom. The rated peak acceleration in the X and Y directions is 1.5g, the rated peak velocity is 100cm / s, the rated stroke is 50cm, and the frequency range of the signal output is 0.1Hz~50Hz. The target engineering structure model system consists of a shear model box (steel), model soil (a mixture of sand and sawdust), and a subway station structure model (microconcrete and galvanized iron wire). The test scale is 1:25, the cross-sectional size of the structural model is 0.836m (width)×0.495m (height), and the model box size is 4.5m (length)×3.5m (width)×3.0m (height).

[0103] Then, three seismic GM signals GM01, GM02 and GM03 were selected, and time domain suppression was performed according to the test time similarity ratio of 1:5.59. The three seismic GM signals were adjusted to acceleration peaks of 0.1g, 0.1g and 0.4g for testing, that is, a comparison of the implementation examples of three working conditions was performed. The following parameters were input through the vibration table system: the signal sampling frequency was 256Hz, the unit impulse response (IR) duration was set to 2s, and the corresponding calculated frequency resolution Δf was 0.25Hz. The FIR filter uses a Hanning window, and the high-pass filter starting frequency is 2Δf, i.e. 0.5Hz. The FIR filter filters the seismic GM signal.

[0104] After installing acceleration sensors A1 and A2 on the vibration table and the top of the target engineering structure model system, respectively, input broadband white noise with an acceleration peak of 0.05g and a frequency range of 0 to 50Hz. Then, based on step S1, the characteristic period T of the target engineering structure model system can be obtained. m =0.21s. By plotting the triple response spectrum of the GM signal of the earthquake, we can get the response spectrum of GM01. The limit frequency of the sensitive area (the limit period value of the sensitive area) is T c =0.07s, T d =2.6s; the limit frequency of the response spectrum sensitive area of ​​GM02 is T c =0.11s, T d =3.0s; the limit frequency of the response spectrum sensitive area of ​​GM03 is T c =0.29s, T d =1.2s.

[0105] Condition 01: Input earthquake motion is GM01. u g0 =0.62cm; the main control parameter discrimination indexes Z1=3.0, Z2=0.08 are calculated, and the speed is used as the main control parameter for reproduction according to Table 1-2.

[0106] Condition 02: Input earthquake motion is GM02. u g0 =2.32cm; the main control parameter discrimination indexes Z1=1.9, Z2=0.07 are calculated, and the speed is used as the main control parameter for reproduction according to Table 1-2.

[0107] Condition 03: Input earthquake motion is GM03. u g0 =15.24cm; the main control parameter discrimination indexes Z1=0.84, Z2=0.07 are calculated, and the speed is used as the main control parameter for reproduction according to Table 1-2.

[0108] In the prior art, acceleration is directly used as the main control parameter for the working conditions 01, 02 and 03.

[0109] The present invention is compared with the prior art through the following contents:

[0110] (1) The error between the acceleration, velocity, displacement signal and the original seismic signal (preset seismic GM signal) of the vibration table surface reproduced by the present invention and the prior art:

[0111] The above error is measured by the error energy index In the formula, X represents the recurrence signal sequence of the prior art, and Y represents the original ground motion signal sequence. The comparison data is shown in Table 3. The target signal is the acceleration, velocity, and displacement signal of the original ground motion.

[0112] Table 3 Error energy EE between the acceleration, velocity, and displacement signals reproduced from the shaking table and the original ground motion signals XY

[0113]

[0114] Table 3 shows that compared with the prior art, the reproduction accuracy of the acceleration signal of the present invention is lower, but the reproduction accuracy of the velocity and displacement signals is greatly improved. The reason is that the velocity parameter is more sensitive to the medium and low frequency signals. The speed as the main control parameter improves the reproduction accuracy of the medium and low frequency signals and significantly weakens the influence of the high-pass filtering.

[0115] (2) Dynamic response of the target engineering structure model system: Figure 3 The peak displacement of the top of the model box measured by the pull-wire displacement meter Δ soil As the response index, the comparative data are shown in Table 4.

[0116] Table 4 Peak displacement Δ of the top of the model box in the present invention and the prior art soil

[0117]

[0118] Table 4 shows that compared with the prior art, the peak displacement response of the top of the model box of the present invention is increased. The reason is that the target engineering structure model system is located in the velocity sensitive area of ​​the response spectrum, and the velocity as the main control parameter improves the reproduction accuracy of the medium and low frequency signals and significantly weakens the influence of the high-pass filtering, thereby increasing the system response and making the system response closer to the true value.

[0119] The main control parameter determination method, computer equipment and medium for seismic signal reproduction provided by the present invention add speed and displacement parameters that are more sensitive to the low-frequency components of the signal into the selection of the main control parameters, and incorporate speed and displacement into the main control parameter selection framework for seismic signal reproduction, thereby avoiding the need for forced high-pass filtering in the process of generating speed and displacement parameters by integrating the acceleration main control parameter of the three-parameter signal generator, thereby improving the reproduction accuracy of medium and low frequency seismic signals; by measuring the characteristic period of the target engineering structure model system and constructing the seismic triple response spectrum, the sensitivity of the target engineering structure model system to the acceleration, velocity and displacement parameters is keenly identified, thereby selecting the most appropriate main control parameters for seismic signal reproduction, meeting the needs of seismic signals with different frequency characteristics and engineering structures with different characteristic periods in the shaking table test, effectively improving the accuracy of the test results, and being suitable for popularization and application.

[0120] The present invention also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the main control parameter determination method for reproducing the above-mentioned seismic signal.

[0121] See also Figure 6 , Figure 6 Schematic diagram of the structure of a computer device provided by this application. Figure 6As shown, the computer device 1000 may include: a processor 1001, a network interface 1004 and a memory 1005. In addition, the computer device 1000 may also include: a user interface 1003, and at least one communication bus 1002. The communication bus 1002 is used to realize the connection and communication between these components. The user interface 1003 may include a display screen (Display) and a keyboard (Keyboard). Optionally, the user interface 1003 may also include a standard wired interface and a wireless interface. The network interface 1004 may optionally include a standard wired interface and a wireless interface (such as a WI-FI interface). The memory 1005 may be a high-speed RAM memory, or a non-volatile memory (non-volatile memory), such as at least one disk storage. The memory 1005 may optionally be at least one storage device located away from the aforementioned processor 1001. As Figure 6 As shown, the memory 1005 as a computer storage medium may include an operating system, a network communication module, a user interface module, and a device control application program.

[0122] exist Figure 6 In the computer device 1000 shown, the network interface 1004 can provide a network communication function; the user interface 1003 is mainly used to provide an input interface for the user; and the processor 1001 can be used to call the device control application stored in the memory 1005 to implement the method for determining the main control parameters of the seismic signal reproduction described in the above embodiment, which will not be repeated here.

[0123] The present invention also provides a computer-readable storage medium, which stores a computer program. The computer program is suitable for being loaded by a processor and executing the method for determining the main control parameters of the seismic signal reproduction described in the above embodiment, which will not be described in detail here.

[0124] The above program can be deployed and executed on a computer device, or deployed and executed on multiple computer devices located in one location, or executed on multiple computer devices distributed in multiple locations and interconnected through a communication network. Multiple computer devices distributed in multiple locations and interconnected through a communication network can form a blockchain network.

[0125] The computer-readable storage medium may be an internal storage unit of the computer device, such as a hard disk or memory of the computer device. The computer-readable storage medium may also be an external storage device of the computer device, such as a plug-in hard disk, a smart memory card (smart media card, SMC), a secure digital (secure digital, SD) card, a flash card (flash card), etc. equipped on the computer device. Further, the computer-readable storage medium may also include both an internal storage unit of the computer device and an external storage device. The computer-readable storage medium is used to store the computer program and other programs and data required by the computer device. The computer-readable storage medium may also be used to temporarily store data that has been output or is to be output.

[0126] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0127] The principles and implementation methods of the present invention are described in this article using specific examples. The description of the above embodiments is only used to help understand the method and core idea of ​​the present invention. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for determining the main control parameters of earthquake signal reproduction, characterized in that: The method comprises: S1: Acquire a response signal of a target engineering structure model system, and determine a characteristic period of the target engineering structure model system according to the response signal; the response signal is a signal generated by the target engineering structure model system after a preset acceleration peak value and white noise of a preset frequency range are input into the target engineering structure model system through a vibration table; the target engineering structure model system is located on the table of the vibration table; S2: calculating the earthquake triple response spectrum of the preset earthquake GM signal according to the parameter information of the preset earthquake GM signal, and determining the boundary period value of the sensitive area of ​​the target engineering structure model system according to the earthquake triple response spectrum; the parameter information includes displacement time history, velocity time history and acceleration time history; the earthquake triple response spectrum includes displacement response spectrum, pseudo velocity response spectrum and pseudo acceleration response spectrum; the boundary period value of the sensitive area includes the boundary period value between the acceleration sensitive area and the velocity sensitive area and the boundary period value between the velocity sensitive area and the displacement sensitive area; S3: Calculating the main control parameter discrimination index according to the characteristic period and the sensitive area boundary period value; S4: Determine the main control parameter used to reproduce the preset seismic GM signal according to the main control parameter discrimination index and the type of the target engineering structure model system; the main control parameter includes displacement, velocity and acceleration; the type of the target engineering structure model system includes a solid material model and a granular material model.

2. The method for determining the main control parameters of earthquake signal reproduction according to claim 1, characterized in that: S1 specifically includes: Acquire a response signal of the target engineering structure model system; the response signal includes a first acceleration response time history and a second acceleration response time history; the first acceleration response time history is acquired by a first acceleration sensor; the second acceleration response time history is acquired by a second acceleration sensor; the first acceleration sensor is located at the top of the target engineering structure model system; the second acceleration sensor is located on the table of the vibration table; Performing Fourier transformation on the first acceleration response time history and the second acceleration response time history respectively to obtain a first Fourier amplitude function and a second Fourier amplitude function; The characteristic period of the target engineering structure model system is determined according to the first Fourier amplitude function and the second Fourier amplitude function.

3. The method for determining the main control parameters of earthquake signal reproduction according to claim 2, characterized in that: Determining the characteristic period of the target engineering structure model system according to the first Fourier amplitude function and the second Fourier amplitude function specifically includes: Calculating an amplitude-frequency characteristic function according to the first Fourier amplitude function and the second Fourier amplitude function; The characteristic period of the target engineering structure model system is determined according to the frequency corresponding to the maximum value of the amplitude-frequency characteristic function.

4. The method for determining the main control parameters of earthquake signal reproduction according to claim 1, characterized in that: Calculating the earthquake triple response spectrum of the preset earthquake GM signal according to the parameter information of the preset earthquake GM signal in S2 specifically includes: Taking the parameter information of the preset seismic GM signal as input, the Duhamel integration algorithm is adopted to calculate the seismic triple response spectrum of the preset seismic GM signal.

5. The method for determining the main control parameters of earthquake signal reproduction according to claim 1, characterized in that: Determining the sensitive zone boundary period value of the target engineering structure model system according to the earthquake triple response spectrum in S2 specifically includes: Determine the peak acceleration, peak velocity and peak displacement of the preset ground motion GM signal according to the ground motion triple response spectrum; The sensitive area boundary period value of the target engineering structure model system is determined according to the peak acceleration, the peak velocity and the peak displacement.

6. The method for determining the main control parameters of earthquake signal reproduction according to claim 5, characterized in that: Determining the sensitive area limit period value of the target engineering structure model system according to the peak acceleration, the peak velocity and the peak displacement specifically includes: Multiplying the peak acceleration, the peak velocity, and the peak displacement by their respective preset magnification factors to obtain amplified peak acceleration, amplified peak velocity, and amplified peak displacement; Based on the amplified peak acceleration, an amplified peak acceleration straight line is drawn in the triple earthquake response spectrum; the amplified peak acceleration straight line is perpendicular to the y-axis of the acceleration coordinate system where the pseudo acceleration response spectrum in the triple earthquake response spectrum is located; Drawing an amplified peak velocity straight line in the triple earthquake response spectrum based on the amplified peak velocity; the amplified peak velocity straight line is perpendicular to the y-axis of the velocity coordinate system where the pseudo velocity response spectrum in the triple earthquake response spectrum is located; An amplified peak displacement straight line is drawn in the triple earthquake response spectrum based on the amplified peak displacement; the amplified peak displacement straight line is perpendicular to the y-axis of the displacement coordinate system where the displacement response spectrum in the triple earthquake response spectrum is located; the x-axis of the velocity coordinate system is in the horizontal direction, the y-axis of the velocity coordinate system is in the vertical direction, the y-axis of the acceleration coordinate system and the y-axis of the displacement coordinate system are both at an angle of 45° to the x-axis of the velocity coordinate system, and the y-axis of the acceleration coordinate system and the y-axis of the displacement coordinate system are perpendicular to each other; The boundary period value between the acceleration sensitive area and the velocity sensitive area is determined according to the amplified peak acceleration straight line and the amplified peak velocity straight line, and the boundary period value between the velocity sensitive area and the displacement sensitive area is determined according to the amplified peak velocity straight line and the amplified peak displacement straight line; the boundary period value between the acceleration sensitive area and the velocity sensitive area and the boundary period value between the velocity sensitive area and the displacement sensitive area constitute the boundary period value of the sensitive area.

7. The method for determining the main control parameters of earthquake signal reproduction according to claim 1, characterized in that: The main control parameter discrimination index includes a first main control parameter discrimination index and a second main control parameter discrimination index; the expression of the first main control parameter discrimination index is: Among them, Z1 is the first main control parameter identification index, T m is the characteristic period of the target engineering structure model system, T c It is the boundary period value between the acceleration sensitive area and the speed sensitive area; The expression of the second main control parameter discrimination index is: Among them, Z2 is the second main control parameter identification index, T d It is the boundary period value between the velocity sensitive area and the displacement sensitive area.

8. The method for determining the main control parameters of earthquake signal reproduction according to claim 1, characterized in that: S4 specifically includes: When the type of the target engineering structure model system is a granular material model, the main control parameter used to reproduce the preset seismic GM signal is determined according to the main control parameter discrimination index and the peak acceleration of the preset seismic GM signal; When the type of the target engineering structure model system is a solid material model, the main control parameter used to reproduce the preset ground motion GM signal is determined according to the main control parameter discrimination index.

9. A computer device, characterized in that: The invention comprises a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the method according to any one of claims 1 to 8.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and the computer program is suitable for being loaded by a processor and executing the method according to any one of claims 1 to 8.

Images