Traveling wave effect analysis method, device, equipment and storage medium for long-span bridges

Abstract

The present invention discloses a method, device, equipment and storage medium for analyzing the traveling wave effect of a long-span bridge, and relates to the technical field of bridge seismic analysis. The method comprises determining the mass damping coefficient and the stiffness damping coefficient according to the circular frequency of the vibration mode of the long-span bridge and the damping ratio of the long-span bridge; using the mass damping coefficient and the stiffness damping coefficient to correct the first time-history curve of the seismic wave, and obtaining the second time-history curve of the seismic wave. The problem in the prior art that only the mass damping coefficient can be corrected, but the mass damping coefficient and the stiffness damping coefficient cannot be corrected at the same time is solved. The second time-history curve of the seismic wave is corrected for the apparent wave velocity, and the third time-history curve of the seismic wave is obtained. The problem in the prior art that the equivalent shear wave velocity is used as the apparent wave velocity, resulting in distortion of the calculation and analysis results, and the inability to correctly evaluate the seismic safety of the structure is solved.

Description

Technical Field

[0001] The present invention relates to the technical field of bridge seismic analysis, and in particular to a method, device, equipment and storage medium for analyzing traveling wave effects of a large-span bridge. Background Art

[0002] Long-span bridges generally have large spans, and seismic design specifications generally require the calculation of multi-point excitation effects of structures. Multi-point excitation effects are mainly achieved by inputting different seismic waves at different supporting points of long-span bridges. When a certain seismic wave is used as a reference and only phase adjustment is performed at different supporting points without spectrum adjustment to consider the spatial variation characteristics of ground motion, it is called the traveling wave effect method, which is the most important influencing factor in the analysis of multi-point excitation effects of long-span bridges.

[0003] The traveling wave effect mainly refers to the reduction of coherence characteristics caused by the time difference of seismic waves from the source to the two measuring points due to different propagation paths. The traveling wave effect method assumes that seismic waves propagate along the ground at a certain speed, the waveform remains unchanged, only the time lag occurs, the attenuation of the amplitude is not considered, and the propagation speed of seismic waves is regarded as a constant.

[0004] In the prior art, when the Rayleigh damping model is used to calculate the structural dynamic response, the error increases with the increase of the mass damping coefficient and the stiffness damping coefficient. There is a problem that only the mass damping coefficient can be corrected, but the mass damping coefficient and the stiffness damping coefficient cannot be corrected at the same time. Summary of the invention

[0005] In view of the defects existing in the prior art, the purpose of the present invention is to provide a method, device, equipment and storage medium for analyzing the traveling wave effect of a long-span bridge, which can solve the problem that when the Rayleigh damping model is used to calculate the dynamic response of the structure, the error increases with the increase of the mass damping coefficient and the stiffness damping coefficient. There is a problem that only the mass damping coefficient can be corrected, and the mass damping coefficient and the stiffness damping coefficient cannot be corrected at the same time.

[0006] In order to achieve the above purpose, the technical solution adopted by the present invention is:

[0007] In the first aspect, this scheme provides a method for analyzing the traveling wave effect of a long-span bridge, including:

[0008] According to the circular frequency of the vibration mode of the long-span bridge and the damping ratio of the long-span bridge, the mass damping coefficient and the stiffness damping coefficient are determined;

[0009] The first time history curve of the seismic wave is corrected by using the mass damping coefficient and the stiffness damping coefficient to obtain the second time history curve of the seismic wave.

[0010] In some optional schemes, according to the formula: Determine the mass damping coefficient and stiffness damping coefficient according to the formula: Obtaining a second time history curve of seismic waves;

[0011] Among them, α is the mass damping coefficient, β is the stiffness damping coefficient, ω i is the circular frequency of the i-th vibration mode of the long-span bridge, ω j is the circular frequency of the jth vibration mode of the long-span bridge, ξ is the damping ratio of the long-span bridge, is the acceleration time history of the second time history curve of the seismic wave, is the acceleration time history curve of the first time history curve, is the velocity time history curve of the first time history curve, u g It is the displacement time history curve of the first time history curve of the seismic wave.

[0012] In some optional solutions, before the first time history curve of the seismic wave is corrected by using the mass damping coefficient and the stiffness damping coefficient, the method further includes:

[0013] The original time history curve of the seismic wave is corrected for zero line drift to obtain the first time history curve of the seismic wave.

[0014] In some optional solutions, performing zero line drift correction on the original time history curve of the seismic wave to obtain the first time history curve of the seismic wave includes:

[0015] Using a Butterworth high-pass filter to perform low-frequency noise reduction on the original acceleration time-history curve to obtain an acceleration time-history curve of the first time-history curve;

[0016] Integrate the acceleration time history curve of the first time history curve once, and subtract the average value of the pre-earthquake part, set the initial velocity to zero, and obtain the velocity time history curve of the first time history curve;

[0017] The velocity time history curve of the first time history curve is integrated once to obtain the displacement time history curve.

[0018] In some optional schemes, after obtaining the second time history curve of the seismic wave, the method further includes: performing apparent wave velocity correction on the second time history curve of the seismic wave to obtain a third time history curve of the seismic wave.

[0019] In some optional solutions, performing apparent wave velocity correction on the second time history curve of the seismic wave to obtain the third time history curve of the seismic wave includes:

[0020] According to the thickness and shear wave velocity of each soil layer, the representative value of the shear wave velocity from the earthquake source to the surface is obtained;

[0021] According to the representative value of shear wave velocity and epicenter distance, the equivalent soil layer apparent wave velocity is obtained;

[0022] The third time history curve of the seismic wave is obtained according to the equivalent soil layer apparent wave velocity, the span of the long-span bridge and the second time history curve of the seismic wave.

[0023] In some optional schemes, according to the formula: Get the representative value of shear wave velocity according to the formula: Obtain the equivalent soil layer apparent wave velocity according to the formula: Obtain the third time history curve of seismic waves;

[0024] in, is the representative value of shear wave velocity, D i is the thickness of the i-th soil layer, i = 1~n, n is the number of soil layers, D is the focal depth, v si Shear wave velocity of the i-th soil layer, is the apparent wave velocity of the equivalent soil layer, S is the epicenter distance, It is the acceleration time history curve of the third time history curve, t is time, and d is the span of the long-span bridge.

[0025] In a second aspect, the present invention provides a device for analyzing traveling wave effects of a long-span bridge, comprising:

[0026] A damping coefficient acquisition module, which is used to determine the mass damping coefficient and the stiffness damping coefficient according to the circular frequency of the vibration mode of the long-span bridge and the damping ratio of the long-span bridge;

[0027] The time history curve acquisition module is used to modify the first time history curve of the seismic wave by using the mass damping coefficient and the stiffness damping coefficient to acquire the second time history curve of the seismic wave.

[0028] In a third aspect, the present invention provides a computer device, comprising a processor, a memory, and a computer program stored in the memory and executable by the processor, wherein when the computer program is executed by the processor, the steps of any of the above-mentioned methods for analyzing the traveling wave effect of a large-span bridge are implemented.

[0029] In a fourth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the steps of any of the above-mentioned methods for analyzing traveling wave effects on large-span bridges are implemented.

[0030] Compared with the prior art, the advantages of the present invention are: this solution determines the mass damping coefficient and the stiffness damping coefficient according to the circular frequency of the vibration mode of the long-span bridge and the damping ratio of the long-span bridge; the mass damping coefficient and the stiffness damping coefficient are used to correct the first time history curve of the seismic wave to obtain the second time history curve of the seismic wave. This solves the problem in the prior art that when the Rayleigh damping model is used to calculate the dynamic response of the structure, the error increases with the increase of the mass damping coefficient and the stiffness damping coefficient. There is a problem that only the mass damping coefficient can be corrected, and the mass damping coefficient and the stiffness damping coefficient cannot be corrected at the same time. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0032] Figure 1 Schematic diagram of a flow chart of a method for analyzing traveling wave effects of a long-span bridge in an embodiment of the present invention;

[0033] Figure 2 A schematic diagram of a multi-layer seismic wave propagation path in an embodiment of the present invention;

[0034] Figure 3 A schematic block diagram of the structure of a computer device in an embodiment of the present invention;

[0035] Figure 4 The figure is a flow chart of a method for analyzing traveling wave effects on large-span bridges taking into account correction of apparent wave velocity in an embodiment of the present invention. DETAILED DESCRIPTION

[0036] In order to make the purpose, technical solution and advantages of the embodiments of the present application clearer, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of this application.

[0037] The embodiments of the present invention are further described in detail below with reference to the accompanying drawings.

[0038] like Figure 1 As shown, in the first aspect, the present invention provides a method for analyzing traveling wave effect of a long-span bridge, comprising:

[0039] S1: Determine the mass damping coefficient and stiffness damping coefficient based on the circular frequency of the vibration mode of the long-span bridge and the damping ratio of the long-span bridge.

[0040] S2: Correcting the first time history curve of the seismic wave by using the mass damping coefficient and the stiffness damping coefficient to obtain the second time history curve of the seismic wave.

[0041] In this embodiment, according to the formula: Determine the mass damping coefficient and stiffness damping coefficient according to the formula: Obtaining a second time history curve of seismic waves;

[0042] Among them, α is the mass damping coefficient, β is the stiffness damping coefficient, ω i is the circular frequency of the i-th vibration mode of the long-span bridge, ω j is the circular frequency of the jth vibration mode of the long-span bridge, ξ is the damping ratio of the long-span bridge, is the acceleration time history of the second time history curve of the seismic wave, is the acceleration time history curve of the first time history curve, is the velocity time history curve of the first time history curve, u g It is the displacement time history curve of the first time history curve of the seismic wave.

[0043] In this embodiment, the damping ratio of the long-span bridge is 0.02 when the main beam is a steel structure, 0.05 when the main beam is a concrete structure, and 0.03 when the main beam is a steel-concrete composite structure.

[0044] In some optional embodiments, before the first time history curve of the seismic wave is corrected by using the mass damping coefficient and the stiffness damping coefficient, the method further includes:

[0045] S0: Perform zero line drift correction on the original time history curve of the seismic wave to obtain the first time history curve of the seismic wave.

[0046] In this embodiment, the first time-history curve may be an original time-history curve.

[0047] Step S0 specifically includes:

[0048] S01: Use a Butterworth high-pass filter to perform low-frequency noise reduction on the original acceleration time-history curve to obtain an acceleration time-history curve of the first time-history curve.

[0049] In this embodiment, the minimum pass frequency of the Butterworth high-pass filter is the minimum value between the first-order natural frequency of the long-span bridge and 0.05 Hz.

[0050] S02: Integrate the acceleration time history curve of the first time history curve once, and subtract the average value of the pre-earthquake part, set the initial velocity to zero, and obtain the velocity time history curve of the first time history curve.

[0051] In this embodiment, the pre-earthquake portion is 0 to 20 seconds.

[0052] S03: Integrate the velocity time history curve of the first time history curve once to obtain a displacement time history curve.

[0053] like Figure 2 The figure shows the propagation path of seismic waves in multiple soil layers. When seismic waves propagate to the interface between different soil layers, refraction occurs.

[0054] like Figure 4 As shown, in some optional embodiments, after obtaining the second time history curve of the seismic wave, the method further includes:

[0055] S3: Perform apparent wave velocity correction on the second time history curve of the seismic wave to obtain the third time history curve of the seismic wave.

[0056] Step S3 specifically includes:

[0057] S31: According to the soil thickness and shear wave velocity of each soil layer, a representative value of the shear wave velocity from the earthquake source to the surface is obtained.

[0058] S32: Obtain the equivalent soil layer apparent wave velocity based on the representative value of shear wave velocity and epicenter distance.

[0059] S33: Obtain a third time history curve of the seismic wave according to the equivalent soil layer apparent wave velocity, the span of the long-span bridge and the second time history curve of the seismic wave.

[0060] In this embodiment, according to the formula: Get the representative value of shear wave velocity according to the formula: Obtain the equivalent soil layer apparent wave velocity according to the formula: Obtain the third time history curve of seismic waves;

[0061] in, is the representative value of shear wave velocity, D i is the thickness of the i-th soil layer, i = 1~n, n is the number of soil layers, D is the focal depth, v si Shear wave velocity of the i-th soil layer, is the apparent wave velocity of the equivalent soil layer, S is the epicenter distance, is the acceleration time history curve of the third time history curve, t is time, and d is the span of the long-span bridge. The soil thickness and shear wave velocity are obtained through geotechnical drilling and indoor dynamic triaxial tests.

[0062] In this embodiment, when the calculated value of the equivalent soil layer apparent wave velocity is lower than 500 m / s, it is taken as 500 m / s, and when the calculated value of the equivalent soil layer apparent wave velocity is higher than 3000 m / s, it is taken as 3000 m / s.

[0063] The third time history curve finally obtained is used in the seismic performance analysis of long-span bridges.

[0064] In summary, the present invention determines the mass damping coefficient and the stiffness damping coefficient according to the circular frequency of the vibration mode of the long-span bridge and the damping ratio of the long-span bridge; uses the mass damping coefficient and the stiffness damping coefficient to correct the first time history curve of the seismic wave to obtain the second time history curve of the seismic wave. This solves the problem in the prior art that when the Rayleigh damping model is used to calculate the dynamic response of the structure, the error increases with the increase of the mass damping coefficient and the stiffness damping coefficient. There is a problem that only the mass damping coefficient can be corrected, and the mass damping coefficient and the stiffness damping coefficient cannot be corrected at the same time.

[0065] Through zero-line drift correction, the problem of non-zero velocity and displacement time series obtained by integration at the end point due to many technical reasons such as low-frequency instrument noise, low-frequency environmental noise, initial acceleration value and initial velocity value, and human operation error is avoided.

[0066] By correcting the apparent wave velocity, the problem in the prior art that the equivalent shear wave velocity is used as the apparent wave velocity, resulting in distortion of the calculation and analysis results and inability to correctly evaluate the seismic safety of the structure is solved.

[0067] In a second aspect, the present invention provides a traveling wave effect analysis device for a long-span bridge, comprising:

[0068] A damping coefficient acquisition module, which is used to determine the mass damping coefficient and the stiffness damping coefficient according to the circular frequency of the vibration mode of the long-span bridge and the damping ratio of the long-span bridge;

[0069] The time history curve acquisition module is used to modify the first time history curve of the seismic wave by using the mass damping coefficient and the stiffness damping coefficient to acquire the second time history curve of the seismic wave.

[0070] It should be noted that those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described devices and modules and units can refer to the corresponding processes in the aforementioned embodiments and will not be repeated here.

[0071] The apparatus provided in the above embodiment may be implemented in the form of a computer program. The computer program may be Figure 3 Runs on the computer device shown.

[0072] See also Figure 3 , Figure 3A schematic block diagram of the structure of a computer device provided in an embodiment of the present application. The computer device may be a terminal.

[0073] like Figure 3 As shown, the computer device includes a processor, a memory, and a network interface connected via a system bus, wherein the memory may include a non-volatile storage medium and an internal memory.

[0074] The non-volatile storage medium can store an operating system and a computer program. The computer program includes program instructions, and when the program instructions are executed, the processor can execute any one of the long-span bridge traveling wave effect analysis methods.

[0075] The processor is used to provide computing and control capabilities and support the operation of the entire computer equipment.

[0076] The internal memory provides an environment for the operation of the computer program in the non-volatile storage medium. When the computer program is executed by the processor, the processor can execute any method for analyzing the traveling wave effect of a large-span bridge.

[0077] The network interface is used for network communication, such as sending assigned tasks, etc. Those skilled in the art will understand that Figure 3 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0078] It should be understood that the processor may be a central processing unit (CPU), and the processor may also be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Among them, the general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0079] In one embodiment, the processor is used to run a computer program stored in a memory to implement the steps of the above-mentioned method for analyzing the traveling wave effect of a long-span bridge.

[0080] An embodiment of the present application further provides a computer-readable storage medium, on which a computer program is stored. The computer program includes program instructions. The method implemented when the program instructions are executed can refer to the various embodiments of the present application.

[0081] The computer-readable storage medium may be an internal storage unit of the computer device described in the foregoing embodiment, 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 (SmartMedia Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), etc., equipped on the computer device.

[0082] In the description of the present application, it should be noted that the terms "upper", "lower", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present application and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present application. Unless otherwise clearly specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be an indirect connection through an intermediate medium, or it can be a connection between the two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the present application can be understood according to the specific circumstances.

[0083] It should be noted that, in this application, relational terms such as "first" and "second" are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "including a..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0084] The above description is only a specific implementation of the present application, so that those skilled in the art can understand or implement the present application. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present application. Therefore, the present application will not be limited to the embodiments shown herein, but will conform to the widest range consistent with the principles and novel features applied for herein.

Claims

1. A method for analyzing traveling wave effects of long-span bridges, characterized in that: include: According to the circular frequency of the vibration mode of the long-span bridge and the damping ratio of the long-span bridge, the mass damping coefficient and the stiffness damping coefficient are determined; The first time history curve of the seismic wave is corrected by using the mass damping coefficient and the stiffness damping coefficient to obtain the second time history curve of the seismic wave; According to the formula: Determine the mass damping coefficient and stiffness damping coefficient according to the formula: Obtaining a second time history curve of seismic waves; Among them, α is the mass damping coefficient, β is the stiffness damping coefficient, ω i is the circular frequency of the i-th vibration mode of the long-span bridge, ω j is the circular frequency of the jth vibration mode of the long-span bridge, ξ is the damping ratio of the long-span bridge, is the acceleration time history of the second time history curve of the seismic wave, is the acceleration time history curve of the first time history curve, is the velocity time history curve of the first time history curve, u g It is the displacement time history curve of the first time history curve of the seismic wave.

2. The method for analyzing traveling wave effect of a long-span bridge according to claim 1, characterized in that: Before the first time history curve of the seismic wave is corrected by using the mass damping coefficient and the stiffness damping coefficient, the method further includes: The original time history curve of the seismic wave is corrected for zero line drift to obtain the first time history curve of the seismic wave.

3. The method for analyzing traveling wave effect of a long-span bridge according to claim 2, characterized in that: The original time history curve of the seismic wave is corrected for zero line drift to obtain the first time history curve of the seismic wave, including: Using a Butterworth high-pass filter to perform low-frequency noise reduction on the original acceleration time-history curve to obtain an acceleration time-history curve of the first time-history curve; Integrate the acceleration time history curve of the first time history curve once, and subtract the average value of the pre-earthquake part, set the initial velocity to zero, and obtain the velocity time history curve of the first time history curve; The velocity time history curve of the first time history curve is integrated once to obtain the displacement time history curve.

4. The method for analyzing traveling wave effect of a long-span bridge according to claim 1, characterized in that: After obtaining the second time history curve of the seismic wave, the method further includes: performing apparent wave velocity correction on the second time history curve of the seismic wave to obtain the third time history curve of the seismic wave.

5. The method for analyzing traveling wave effect of a long-span bridge according to claim 4, characterized in that: Performing apparent wave velocity correction on the second time history curve of the seismic wave to obtain the third time history curve of the seismic wave includes: According to the thickness and shear wave velocity of each soil layer, the representative value of the shear wave velocity from the earthquake source to the surface is obtained; According to the representative value of shear wave velocity and epicenter distance, the equivalent soil layer apparent wave velocity is obtained; The third time history curve of the seismic wave is obtained according to the equivalent soil layer apparent wave velocity, the span of the long-span bridge and the second time history curve of the seismic wave.

6. The method for analyzing traveling wave effect of a long-span bridge according to claim 5, characterized in that: According to the formula: Get the representative value of shear wave velocity according to the formula: Obtain the equivalent soil layer apparent wave velocity according to the formula: Obtain the third time history curve of seismic waves; in, is the representative value of shear wave velocity, D i is the thickness of the i-th soil layer, i = 1~n, n is the number of soil layers, D is the focal depth, v si Shear wave velocity of the i-th soil layer, is the apparent wave velocity of the equivalent soil layer, S is the epicenter distance, It is the acceleration time history curve of the third time history curve, t is time, and d is the span of the long-span bridge.

7. A device for analyzing traveling wave effects of long-span bridges for implementing the method described in claim 1, characterized in that: include: A damping coefficient acquisition module, which is used to determine the mass damping coefficient and the stiffness damping coefficient according to the circular frequency of the vibration mode of the long-span bridge and the damping ratio of the long-span bridge; The time history curve acquisition module is used to modify the first time history curve of the seismic wave by using the mass damping coefficient and the stiffness damping coefficient to acquire the second time history curve of the seismic wave.

8. A computer device, characterized in that: The computer device includes a processor, a memory, and a computer program stored in the memory and executable by the processor, wherein when the computer program is executed by the processor, the steps of the method for analyzing the traveling wave effect of a large-span bridge as described in any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, wherein when the computer program is executed by a processor, the steps of the method for analyzing the traveling wave effect of a long-span bridge as claimed in any one of claims 1 to 6 are implemented.

Images