Seismic dynamics simulation method and system for gravity dam based on time domain deconvolution algorithm

Abstract

The present invention discloses a gravity dam seismic dynamic simulation method and system based on a time-domain deconvolution algorithm, and belongs to the technical field of seismic analysis of gravity dams. The existing seismic analysis schemes for gravity dams do not take into account the coupling between the gravity dam, reservoir water, and dam foundation, and cannot accurately describe the structural response of the gravity dam in an earthquake, which affects the accuracy of the seismic analysis of the gravity dam. The seismic dynamic simulation method for a gravity dam based on a time-domain deconvolution algorithm of the present invention fully considers the compressibility of reservoir water, dam foundation damping, and the interaction between the gravity dam, dam foundation, and reservoir water, and introduces time-domain deconvolution to calculate seismic effects, and then performs dynamic time-history calculations for the seismic water conditions to which the gravity dam is subjected, thereby being able to accurately describe the structural response of the gravity dam in an earthquake, facilitating efficient and reliable seismic simulation analysis of the gravity dam, thereby effectively improving the accuracy of the seismic simulation calculation results of the gravity dam.

Description

Technical Field

[0001] The invention relates to a gravity dam seismic dynamics simulation method and system based on a time domain deconvolution algorithm, and belongs to the technical field of gravity dam seismic analysis. Background Art

[0002] As an important part of the water conservancy project hub, the gravity dam plays a huge engineering benefit in regulating the spatial and temporal distribution of water resources. Looking at the concrete gravity dam cases that have been subjected to strong earthquakes at home and abroad, there are weak parts of the gravity dam that are earthquake-resistant. When an earthquake occurs, the key parts of the gravity dam may crack. Once it is severely damaged by a strong earthquake, it will also cause the dam body to become unstable and the bank slope to collapse. Not only will the project itself be severely damaged, but more importantly, it will lead to unimaginable secondary disasters. Therefore, it is very important to analyze the seismic safety of the gravity dam, do a good job in the seismic design of the concrete dam, and propose feasible reinforcement measures.

[0003] Chinese literature (Li Bo, Xin Quancai. Analysis of seismic dynamic response of gravity dam considering reservoir water compressibility [J]. People's Yellow River, 2018 (8): 116-119.) To study the influence of reservoir water compressibility on the dynamic characteristics of the dam body, a gravity dam was taken as an example to establish four different finite element models of reservoir water action. The natural frequency of the dam body, the peak stress and displacement of the key parts of the dam body, and the dynamic water pressure distribution on the dam surface obtained by comparing the traditional added mass model, the half-reduced added mass model, and the model ignoring the compressibility of reservoir water showed that: the added mass model exaggerates the effect of water to a certain extent, resulting in the largest decrease in the fundamental frequency of the dam body, increasing the peak stress at the heel of the dam, and significantly amplifying the dynamic water pressure in the middle and upper parts of the dam body, but weakening the displacement of the dam top to a certain extent; the compressible reservoir water model can better reflect the interaction between reservoir water and the dam body, so the compressibility of reservoir water should not be ignored in the calculation.

[0004] However, the above scheme only considers the coupling effect of dam body and reservoir water, and does not consider the influence of dam foundation. However, at the actual earthquake site, the earthquake first acts on the dam foundation, and then the dam foundation drives the movement and deformation of the dam body and surrounding geographical structures, resulting in the response of the water body. The moving water will react on the dam surface, thereby affecting the structural response of the gravity dam. Therefore, if the mutual coupling effect between the gravity dam and reservoir water and dam foundation structure is not considered, it will affect the accuracy of the seismic dynamic simulation of the gravity dam.

[0005] Furthermore, the above scheme is based on the design response spectrum in the "Code for Seismic Design of Hydraulic Structures in Hydropower Projects" and artificially synthesizes seismic waves. Although the response spectrum method takes into account the resonance effect caused by the dynamic characteristics of the structure, it still regards the seismic inertia force as static force, and cannot reflect the important influence of the duration of seismic motion on the degree of structural damage. It can only obtain the maximum value of the structural vibration response and only conducts static analysis on the gravity dam. It cannot accurately describe the structural response of the gravity dam in an earthquake, which affects the accuracy of the seismic analysis of the gravity dam.

[0006] The information disclosed in this Background Art is only for understanding the background of the inventive concept and therefore it may include information that does not constitute the prior art. Summary of the invention

[0007] In response to the above problem or one of the above problems, an object of the present invention is to provide a gravity dam seismic dynamic simulation method and system based on a time-domain deconvolution algorithm. By creating a dam body generation simulation model, a dam body modal analysis model, a seismic action calculation model, and a fluid-solid coupling response calculation model, the compressibility of reservoir water, the dam foundation damping, and the interaction between the gravity dam, the dam foundation and the reservoir water are fully considered, and time-domain deconvolution is introduced to calculate the seismic action, and then the modal analysis of the overall structure of the gravity dam, the dam foundation and the reservoir water is performed. The dynamic time history calculation is performed for the seismic water conditions to which the gravity dam is subjected. Therefore, the structural response of the gravity dam in an earthquake can be accurately described, and efficient and reliable seismic simulation analysis of the gravity dam is facilitated, thereby effectively improving the accuracy of the seismic simulation calculation results of the gravity dam, which has important guiding significance for the structural seismic design, operation and maintenance, and reinforcement engineering of the gravity dam.

[0008] In response to the above problem or one of the above problems, the second object of the present invention is to provide a gravity dam seismic dynamic simulation method and system based on a time domain deconvolution algorithm, which takes the seismic action calculated based on the time domain deconvolution algorithm as input, considers the coupling of gravity dam-reservoir water-dam foundation, and uses simulation means to perform dynamic time-history analysis and calculation, thereby effectively improving the accuracy of the seismic calculation of the gravity dam.

[0009] To achieve one of the above purposes, the first technical solution of the present invention is:

[0010] The seismic dynamics simulation method of gravity dam based on time domain deconvolution algorithm includes the following contents:

[0011] By generating a simulation model of the dam body created in advance, a dam body simulation object for earthquake-resistant dynamics simulation is established according to the gravity dam design data and geological information, which at least includes the gravity dam, reservoir water and dam foundation;

[0012] Using the previously created dam body modal analysis model, taking into account the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption effect of the reservoir water boundary, the modal analysis of the dam body simulation object including the gravity dam, reservoir water and dam foundation is performed to calculate the natural frequency and structural modal vibration shape of the dam body simulation object;

[0013] Based on the earthquake action calculation model created in advance and the geological information around the gravity dam, the earthquake acceleration time history is determined, and the earthquake action input data is calculated through the time domain deconvolution algorithm;

[0014] According to the fluid-solid coupling response calculation model created in advance, the seismic action input data is applied to the dam body simulation object, and the dynamic time history analysis is carried out according to the natural frequency and structural modal vibration shape of the dam body simulation object to obtain the seismic response data of the dam body simulation object, and complete the seismic dynamic simulation of the gravity dam based on the time domain deconvolution algorithm.

[0015] The present invention creates a dam body generation simulation model, a dam body modal analysis model, an earthquake action calculation model, and a fluid-solid coupling response calculation model, fully considers the compressibility of reservoir water, dam foundation damping, and the interaction between the gravity dam, the dam foundation and reservoir water, introduces time domain deconvolution to calculate the earthquake action, and then performs modal analysis on the overall structure of the gravity dam, the dam foundation and reservoir water, and performs dynamic time history calculation on the seismic water conditions experienced by the gravity dam. Therefore, the structural response of the gravity dam in an earthquake can be accurately described, and efficient and reliable earthquake-resistant simulation analysis of the gravity dam is facilitated, thereby effectively improving the accuracy of the earthquake-resistant simulation calculation results of the gravity dam, and having important guiding significance for the structural earthquake-resistant design, operation and maintenance, and reinforcement engineering of the gravity dam.

[0016] As the preferred technical measures:

[0017] The method of establishing a dam body simulation object for seismic dynamics simulation by generating a simulation model of the dam body created in advance and according to the gravity dam design data and geological information is as follows:

[0018] Based on the gravity dam design engineering drawings and geological survey reports, establish two-dimensional plane simulation objects of the gravity dam, reservoir water and dam foundation;

[0019] Then the two-dimensional plane simulation object is divided into unit grids to obtain the first dam body grid;

[0020] Based on the location information of the boundary between reservoir water, gravity dam and dam foundation, set the number of encryption layers and the size of the encrypted area grid;

[0021] According to the number of encryption layers and the mesh size of the encryption area, the first dam body mesh is encrypted to obtain the second dam body mesh;

[0022] Based on the interface information between reservoir water and bedrock, additional grid nodes are set on the basis of the second dam body grid to characterize the properties of the reflected wave, thereby obtaining the third dam body grid;

[0023] Based on the third dam body grid, the interface group, free boundary group and absorbing boundary group of the gravity dam are coupled to distinguish the propagation characteristics of seismic waves, so as to obtain the dam body simulation object for fluid and structure coupling calculation.

[0024] As the preferred technical measures:

[0025] Using the previously created dam body modal analysis model, considering the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption effect of the reservoir water boundary, the modal analysis of the dam body simulation object including the gravity dam, reservoir water and dam foundation is performed. The method for calculating the natural frequency and structural modal vibration shape of the dam body simulation object is as follows:

[0026] Based on the structural material properties of the gravity dam, set the structural density, Young's modulus and Poisson's ratio of the dam simulation object;

[0027] According to the structural density, Young's modulus and Poisson's ratio, the mass matrix and stiffness matrix of the dam simulation object are established;

[0028] Calculate the hydrodynamic pressure field of the fluid domain based on the medium characteristics and wave impedance data of wave propagation;

[0029] Based on the hydrodynamic pressure field of the fluid and the wave transmission and reflection characteristics of the interface, and taking into account the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water and the absorption effect of the reservoir water boundary, the boundary conditions and constraints of the bottom of the dam foundation and the upstream and downstream surfaces are determined;

[0030] According to the boundary conditions and constraints, the damping matrix and constraint excitation vector of the dam simulation object are established;

[0031] According to the mass matrix, stiffness matrix, damping matrix and constraint excitation vector of the dam simulation object, the linear system of the gravity dam is established;

[0032] The modal solution of the gravity dam linear system is performed to obtain the natural frequency and structural mode shape of the dam simulation object.

[0033] As the preferred technical measures:

[0034] According to the medium characteristics and wave impedance data of wave propagation, the method for calculating the fluid hydrodynamic pressure field in the fluid domain is as follows:

[0035] The Helmholtz equation is established based on the medium characteristics and wave impedance of wave propagation;

[0036] Solve the Helmholtz equation to obtain the fluid wave vector;

[0037] Based on the fluid characteristics of reservoir water being uniform, compressible, irrotational and inviscid, a wave equation satisfied by water pressure is established;

[0038] According to the wave equation, construct the propagation expression of the fluid domain;

[0039] Substitute the fluid wave vector into the propagation expression to obtain the water pressure data in the fluid domain;

[0040] The water pressure data is subjected to inverse Fourier transform to obtain the fluid dynamic water pressure field in the time domain.

[0041] As the preferred technical measures:

[0042] Based on the hydrodynamic pressure field of the fluid and the wave transmission and reflection characteristics of the interface, and considering the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption effect of the reservoir water boundary, the method for determining the boundary conditions and constraints of the bottom of the dam foundation and the upstream and downstream surfaces is as follows:

[0043] Based on the hydrodynamic pressure field of the fluid and the wave transmission and reflection characteristics of the interface, the first node boundary condition on the interface between the dam body and the reservoir water is set;

[0044] Considering the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption of the reservoir water boundary, the influence of the flexible bedrock or rigid bedrock on the boundary conditions of the first node is obtained;

[0045] Based on the impact information, the reservoir water absorption coefficient is set to represent the pressure gradient;

[0046] On the free surface of the reservoir water, ignoring the surface wave effect of the reservoir water, set the second node boundary condition;

[0047] At the upstream position of the reservoir water far from the dam surface, that is, at 2.5 times the gravity dam elevation, a stationary boundary is imposed to obtain the third node boundary condition;

[0048] The first node boundary condition, the second node boundary condition, the third node boundary condition and the reservoir water absorption coefficient are coupled to obtain the boundary conditions and constraints of the bottom and upstream and downstream surfaces of the dam foundation.

[0049] As the preferred technical measures:

[0050] Based on the earthquake action calculation model created in advance, the earthquake acceleration time history is determined according to the geological information around the gravity dam, and the method of calculating the earthquake action input data through the time domain deconvolution algorithm is as follows:

[0051] According to the earthquake information of the actual geographical location of the gravity dam, the free-field earthquake acceleration time history function required for the dynamic time history analysis is determined, which includes the free-field horizontal earthquake acceleration time history function and the free-field vertical earthquake acceleration time history function;

[0052] The free-field horizontal seismic acceleration time history function and the free-field vertical seismic acceleration time history function are integrated by the time-domain deconvolution algorithm to obtain the free-field horizontal seismic velocity time history function and the free-field vertical seismic velocity time history function;

[0053] According to the free-field horizontal seismic velocity time history function and the free-field vertical seismic velocity time history function, the horizontal velocity time history and the vertical velocity time history of the earthquake's horizontal wave and vertical wave at the bottom of the dam foundation are calculated;

[0054] Based on the horizontal velocity time history, vertical velocity time history, equivalent shear and compression wave velocities in the node force layer, and considering the viscous absorption characteristics of the dam foundation boundary, the node force calculation equations are constructed;

[0055] Based on the node force calculation equations, the node forces in the vertical direction and the horizontal direction caused by the incident seismic wave are calculated, that is, the seismic action input data.

[0056] As the preferred technical measures:

[0057] According to the previously created fluid-solid coupling response calculation model, the dynamic time history analysis is carried out to obtain the seismic response data of the dam simulation object as follows:

[0058] Apply earthquake input data to the dam simulation object, and perform dynamic time history analysis based on the natural frequency and structural modal vibration shape of the dam simulation object to obtain the displacement field of the dam simulation object at several moments;

[0059] Based on the displacement field, calculate the earthquake response acceleration, velocity, displacement and stress data of the dam simulation object at different times;

[0060] According to the earthquake response acceleration, velocity, displacement and stress data, the impact of earthquake on the structural stability of the dam simulation object is evaluated to obtain the seismic response data.

[0061] As the preferred technical measures:

[0062] The method to obtain the displacement field of the dam simulation object at several moments is as follows:

[0063] According to the propagation characteristics of the interface position of the dam simulation object, the displacement potential and displacement potential are established;

[0064] Based on the displacement potential and displacement potential, as well as the natural frequency and structural mode vibration shape, the weak form of the fluid-solid coupling equation is obtained, which includes the solid domain calculation formula and the fluid domain calculation formula;

[0065] Using the solid domain calculation formula and the fluid domain calculation formula, a linear system consisting of symmetric matrices is obtained;

[0066] On the basis of the linear system, impedance boundaries are set to obtain a set of fluid-solid coupling discrete equations with a time format;

[0067] Decompose the earthquake action input data into the force on the dam foundation close to the upstream surface, the force on the bottom of the dam foundation and the force on the dam foundation close to the downstream surface;

[0068] The structural state at the initial moment is obtained, the calculation time step is selected, and the force of the dam foundation close to the upstream surface, the force at the bottom of the dam foundation, and the force close to the downstream surface of the dam foundation are input into the fluid-solid coupling discrete equation group for solution to obtain the displacement field of the dam simulation object at each moment.

[0069] As the preferred technical measures:

[0070] According to the earthquake response acceleration, velocity, displacement and stress data, the impact of earthquake on the structural stability of the dam simulation object is evaluated and the method for obtaining the seismic response data is as follows:

[0071] Calculate the stability coefficient of the dam simulation object based on the earthquake response acceleration, velocity, displacement and stress data;

[0072] Based on the stress data, the calculation results of the stress tensor of each node of the dam simulation object are obtained;

[0073] According to the calculation results of the stress tensor, three stress results are obtained at each node position, including horizontal stress component, vertical stress component and shear stress component;

[0074] Construct a stress tensor matrix according to horizontal stress components, vertical stress components and shear stress components;

[0075] Calculate the stress tensor matrix to obtain the eigenvalue of the stress tensor and the principal stress at each position of the dam simulation object;

[0076] Compare the principal stress with the tensile strength and compressive strength of the concrete material to obtain the comparison results;

[0077] The stability coefficients and comparison results are summarized to obtain seismic response data, which is used to evaluate the impact of earthquakes on the structural stability of the dam simulation object.

[0078] To achieve one of the above purposes, the second technical solution of the present invention is:

[0079] The gravity dam seismic dynamics simulation system based on time domain deconvolution algorithm includes:

[0080] one or more processors;

[0081] A storage device for storing one or more programs;

[0082] When the one or more programs are executed by the one or more processors, the one or more processors implement the above-mentioned gravity dam seismic dynamics simulation method based on the time domain deconvolution algorithm.

[0083] To achieve one of the above purposes, the second technical solution of the present invention is:

[0084] Compared with the prior art solutions, the present invention has the following beneficial effects:

[0085] The present invention creates a dam body generation simulation model, a dam body modal analysis model, an earthquake action calculation model, and a fluid-solid coupling response calculation model, fully considers the compressibility of reservoir water, dam foundation damping, and the interaction between the gravity dam, the dam foundation and reservoir water, introduces time domain deconvolution to calculate the earthquake action, and then performs modal analysis on the overall structure of the gravity dam, the dam foundation and reservoir water, and performs dynamic time history calculation on the seismic water conditions experienced by the gravity dam. Therefore, the structural response of the gravity dam in an earthquake can be accurately described, and efficient and reliable earthquake-resistant simulation analysis of the gravity dam is facilitated, thereby effectively improving the accuracy of the earthquake-resistant simulation calculation results of the gravity dam, and having important guiding significance for the structural earthquake-resistant design, operation and maintenance, and reinforcement engineering of the gravity dam.

[0086] Furthermore, the present invention uses seismic effects based on time domain deconvolution calculations as input, considers the coupling of gravity dam-reservoir water-dam foundation, and uses simulation methods to perform dynamic time-history analysis and calculations, which can effectively solve the problem that the existing gravity dam seismic calculations have too many assumptions and the calculation results are not comprehensive and accurate enough. BRIEF DESCRIPTION OF THE DRAWINGS

[0087] Figure 1 It is a schematic diagram of a first flow chart of the seismic dynamics simulation method of a gravity dam according to the present invention;

[0088] Figure 2 It is a second flow chart of the gravity dam seismic dynamics simulation method of the present invention;

[0089] Figure 3 A schematic diagram of a two-dimensional simulation grid of a gravity dam according to the present invention;

[0090] Figure 4 A schematic diagram of the boundary conditions of the fluid-solid interface of the present invention. DETAILED DESCRIPTION

[0091] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0092] On the contrary, the present invention covers any substitution, modification, equivalent method and scheme made on the essence and scope of the present invention as set by the claims. Further, in order to make the public have a better understanding of the present invention, some specific details are described in detail in the detailed description of the present invention below. Those skilled in the art can fully understand the present invention without the description of these details.

[0093] Unless otherwise specified, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art to which the present invention pertains. The terms used herein are only for the purpose of describing specific embodiments and are not intended to limit the present invention. The term "or / and" as used herein includes any and all combinations of one or more of the associated listed items.

[0094] like Figure 1 As shown, the first specific embodiment of the gravity dam seismic dynamics simulation method based on the time domain deconvolution algorithm of the present invention is:

[0095] The seismic dynamics simulation method of gravity dam based on time domain deconvolution algorithm includes the following contents:

[0096] By generating a simulation model of the dam body created in advance, a dam body simulation object for earthquake-resistant dynamics simulation is established according to the gravity dam design data and geological information, which at least includes the gravity dam, reservoir water and dam foundation;

[0097] Using the previously created dam body modal analysis model, taking into account the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption effect of the reservoir water boundary, the modal analysis of the dam body simulation object including the gravity dam, reservoir water and dam foundation is performed to calculate the natural frequency and structural modal vibration shape of the dam body simulation object;

[0098] Based on the earthquake action calculation model created in advance and the geological information around the gravity dam, the earthquake acceleration time history is determined, and the earthquake action input data is calculated through the time domain deconvolution algorithm;

[0099] According to the fluid-solid coupling response calculation model created in advance, the seismic action input data is applied to the dam body simulation object, and the dynamic time history analysis is carried out according to the natural frequency and structural modal vibration shape of the dam body simulation object to obtain the seismic response data of the dam body simulation object, and complete the seismic dynamic simulation of the gravity dam based on the time domain deconvolution algorithm.

[0100] like Figure 2 As shown, the second specific embodiment of the gravity dam seismic dynamics simulation method based on the time domain deconvolution algorithm of the present invention is:

[0101] The seismic dynamics simulation method of a gravity dam based on a time domain deconvolution algorithm includes the following steps:

[0102] Step 1: According to the gravity dam design engineering drawings and geological survey reports, combined with the actual engineering situation, establish a dam body simulation object for calculation and a two-dimensional simulation grid for finite element analysis, including the gravity dam, water body and dam foundation, which includes the following contents:

[0103] Based on the design drawings, establish two-dimensional plane simulation objects of the gravity dam, water body and dam foundation; then divide the two-dimensional plane simulation objects into unit grids. Encrypt the grid at the junction of the water body, gravity dam and dam foundation, with no less than 3 layers of encryption and the grid size in the encrypted area less than 50cm. Set additional grid nodes at the interface between the water body and bedrock to set the reflection wave properties. At the same time, set the interface group, free boundary group and absorbing boundary group of each structure to distinguish the propagation characteristics of seismic waves, and obtain a two-dimensional simulation grid that is convenient for finite element analysis for fluid and structure coupling calculations. See Figure 3 .

[0104] Step 2: Considering the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption effect of the reservoir water boundary, a modal analysis is performed on the dam body simulation object including the gravity dam, reservoir water, and dam foundation, and the natural frequency and structural mode vibration shape under normal water storage level are calculated.

[0105] Step 3: Determine the earthquake acceleration time history based on the geographical location of the dam, calculate the earthquake action input at the bottom depth of the model dam foundation through the deconvolution algorithm, and perform dynamic time history analysis.

[0106] Step 4: Post-process the simulation analysis results to calculate the earthquake response acceleration, velocity, displacement and stress of the structure at different times, and evaluate the structural stability and concrete damage and cracking.

[0107] This embodiment uses the seismic action based on time domain deconvolution calculation as input, considers the coupling of gravity dam-reservoir water-dam foundation, and uses simulation methods to perform dynamic time history analysis and calculation, which can effectively solve the problem that the existing gravity dam seismic calculation has too many assumptions and the calculation results are not comprehensive and accurate enough.

[0108] A specific embodiment of applying the present invention to conduct earthquake-resistant dynamic simulation on a concrete gravity dam of a reservoir in the east:

[0109] The method for conducting seismic dynamic simulation on a concrete gravity dam of a reservoir in the eastern part of China includes the following steps:

[0110] Step 1: Generate a simulation model through the dam body. Based on the concrete solid gravity dam drawings of a reservoir in the east, establish a two-dimensional plane simulation object of the gravity dam, water body and dam foundation according to the normal water storage condition of the gravity dam. The dam is designed for a flood that occurs once in 1,000 years, and the normal water storage level is 100m. The crest of the gravity dam is 200m long, the crest elevation is 130m, the upstream is vertical, and the downstream dam slope is 1:0.7. A 1.0m high wave-breaking wall is set on the upstream side of the dam crest, and a 0.8m high railing is set on the downstream side. A 3×3m city gate type grouting and drainage corridor is set downstream of the dam axis. The center line of the corridor is 3.5m away from the upstream dam axis, and the earthquake intensity is VII.

[0111] The two-dimensional plane simulation object is divided into quadrilateral / triangular mixed unit meshes, with quadrilateral meshes used for reservoir water and bedrock and triangular meshes used for the dam body. The mesh is encrypted at the interface between the water body, gravity dam and dam foundation, with no less than 5 layers of encryption and the mesh size in the encrypted area less than 40cm. Additional double mesh nodes are set at the interface between the water body and bedrock to set the reflection wave properties; at the same time, the interface group, free boundary group and absorbing boundary group of each structure are set to distinguish the propagation and reflection characteristics of seismic waves in different positions in the material, so as to obtain a two-dimensional simulation mesh for fluid and structure coupling calculation, with a total of 3208 nodes.

[0112] Step 2: Use the dam body modal analysis model to perform modal analysis on the dam body simulation object including the gravity dam, reservoir water and dam foundation, which includes the following contents:

[0113] The first step is to determine the material properties of various structures according to the actual project conditions, among which the "dam body material parameters" include: density, elastic modulus, Poisson's ratio and damping; the "dam foundation material parameters" include: density, elastic modulus, Poisson's ratio and damping; the "general material parameters" include: friction coefficient and cohesion.

[0114] The second step is to set the medium characteristics and wave impedance of wave propagation. Special considerations are required for the fluid domain, and the fluid hydrodynamic pressure field Satisfies the Helmholtz equation. For all virtual pressure fields , and its calculation formula is as follows:

[0115]

[0116] in, is the Laplace operator, is the fluid domain, is a solid domain, Representing the fluid domain The borders of is the boundary normal vector, is the wave speed in the fluid, For time.

[0117] Based on the assumption that water is a uniform, compressible, irrotational, and inviscid ideal fluid, the water pressure field Satisfying the wave equation, the calculation formula is as follows:

[0118]

[0119] Water pressure field The propagation mode in the fluid satisfies the following expression:

[0120]

[0121] in, is the amplitude of the wave, is an imaginary unit, In space The wave vectors in three directions (horizontal, longitudinal, and vertical) are is the angular frequency.

[0122] Then, substituting the relevant data into the Helmholtz equation, we get the following expression:

[0123]

[0124] After inverse Fourier transform, we can get the expression in time domain:

[0125]

[0126] The third step is to set the boundary conditions of the fluid-solid interface based on the wave transmission and reflection characteristics of the fluid-solid interface of the gravity dam. Figure 4 , which includes the following:

[0127] (1) The nodes on the interface S1 between the dam body and the reservoir water follow the natural boundary, and the calculation formula is as follows:

[0128]

[0129] in, represents the normal vector of the dam-reservoir interface node, is the fluid density, is the normal acceleration.

[0130] (2) On the free surface S2 of the reservoir water, ignoring that the water body is zero potential and has no viscosity, and ignoring the surface wave effect of the water body, the necessary boundary conditions are:

[0131]

[0132] (3) At a position S3 far enough upstream from the dam surface, that is, 2.5 times the elevation of the gravity dam, a static boundary is imposed, that is, the normal acceleration is 0, and the calculation formula is as follows:

[0133]

[0134] in, Indicates the normal direction at infinity at the end of the reservoir.

[0135] (4) The boundary condition of the interface S4 between the reservoir bottom and the bedrock is affected by the type of bedrock (flexible or rigid). To express the relationship between the pressure gradient and the normal acceleration, the calculation formula is as follows:

[0136]

[0137] in, is the normal direction of the reservoir bottom surface.

[0138] The fourth step is to impose horizontal and vertical displacement constraints on the bottom and upstream and downstream surfaces of the dam foundation.

[0139] The fifth step is to assemble the mass matrix, stiffness matrix and damping matrix for modal calculation, as well as the external constraint excitation vector, which includes the following:

[0140] Construct the Rayleigh damping formula, which is expressed as follows:

[0141]

[0142] in, is the damping matrix of the structure, is the mass matrix of the structure, is the stiffness matrix of the structure, and is the coefficient, and its calculation formula is:

[0143]

[0144]

[0145] in, and Respectively Order and The angular frequency of the order, and Gravity dam structure Order and The modal damping ratio of the order; Based on the structural characteristics of the gravity dam, Take 1, Take 2.

[0146] The sixth step is to perform modal solution, which includes the following:

[0147] For a general linear system, the equations to be solved under the constrained mode are:

[0148]

[0149] in, For external incentives, is the displacement vector.

[0150] In the case of fluid-solid coupling, according to the principle of virtual work, for the solid domain, the calculation formula is as follows:

[0151]

[0152] For fluid, the calculation formula is as follows:

[0153]

[0154]

[0155] in, represents the computational domain, is the density of the solid, is the displacement of the structure, and are the virtual displacement field of the solid and the virtual pressure field of the fluid respectively.

[0156] From the above formula, the discretized fluid-solid coupling linear equations can be written as follows:

[0157]

[0158] in, and They are obtained by integrating the following bilinear equations, and the integral formula is as follows:

[0159]

[0160]

[0161] Damping Matrix It is also calculated by the bilinear equation, and its integral formula is as follows:

[0162]

[0163] After modal calculation, the modal vibration shapes and normalized modal masses of the gravity dam and dam foundation at different modal frequencies are obtained, and the normalized modal masses are accumulated to obtain the cumulative modal mass. These results can help analyze the dynamic characteristics of the gravity dam structure, and the normalized modal mass and cumulative modal mass can help determine whether the calculated frequency range is sufficient to analyze the structural response of the dam. In this embodiment, the first 10 modal results are calculated, and the cumulative modal mass of the first three modes is greater than 90%, indicating that the 10th mode meets the calculation requirements. The modal vibration shape can serve as the basis for the subsequent step three.

[0164] Step 3: Determine the earthquake acceleration time history according to the geographical location of the dam, calculate the earthquake action input at the bottom depth of the model dam foundation through the deconvolution algorithm, and perform dynamic time history analysis, which includes the following steps:

[0165] The first step is to determine the free-field seismic acceleration time history required for dynamic time history analysis based on the seismic information of the actual geographical location of the gravity dam.

[0166] The seismic input of the embodiment intercepts a typical earthquake acceleration time history in the area where the gravity dam is located, and the calculation time domain is .

[0167] In the second step, the free-field earthquake time history function is calculated by deconvolution to obtain the earthquake action input at the bottom depth of the model dam foundation.

[0168] It is assumed that the horizontal earthquake action is mainly composed of shear waves (S waves), while the vertical earthquake action is mainly composed of longitudinal waves (P waves). The free-field horizontal and vertical earthquake acceleration time history functions are expressed as and Integrating the acceleration time history function, we get the seismic velocity time history function and A viscous absorbing boundary is applied to the dam foundation boundary to absorb the waves incident vertically on the surface.

[0169] The vertical and horizontal node forces caused by the incident seismic waves are applied to the bottom of the model dam foundation:

[0170]

[0171]

[0172] Where: Indicates the grid width, is the material density, and are the velocity time histories of the S-wave and P-wave at the bottom of the dam foundation, respectively. In the first calculation iteration, it is assumed that = and ; and are the equivalent shear and compression wave velocities in the layer where the nodal forces are applied.

[0173] The velocity function at the bottom of the dam foundation can be expressed by fast Fourier transform to express its amplitude and phase, and its calculation formula is as follows:

[0174]

[0175]

[0176] In the formula Represent the S wave and P wave phases respectively.

[0177] The corresponding velocity response of S-wave and P-wave on the surface at this time for:

[0178]

[0179]

[0180] In the formula Represent the response phases of S wave and P wave respectively.

[0181] S-wave and P-wave transfer functions and for:

[0182]

[0183]

[0184] The effective gain function is:

[0185]

[0186]

[0187] The effective phase shift is:

[0188]

[0189]

[0190] A new velocity function for the dam base bottom is generated according to the following equation:

[0191]

[0192]

[0193] In the formula , represents the phase of S wave and P wave, and is the amplitude correction factor and can be determined as follows:

[0194]

[0195]

[0196] Then, an inverse fast Fourier transform is performed to obtain the new incident seismic wave velocity time history function from the above steps, and then the new equivalent nodal force is calculated and , acting on the bottom of the dam foundation, calculate the earthquake velocity time history at the surface at this time, and compare it with , Make a comparison.

[0197] Repeat the above steps until the calculated seismic velocity time history at the surface is consistent with , The match is good. and This is the required equivalent nodal force, applied to the bottom of the model dam foundation as the input of the earthquake action.

[0198] Step 3: Select the calculation time step , and the time interpolation format, the initial time is known The structural state can be calculated The vibration state at the time of discretization is The equilibrium equations are solved one by one on the modal basis to obtain the displacement field of the gravity dam at each moment.

[0199] More specifically, in order to more accurately describe the propagation characteristics of the material interface and the subsequent time-history calculation solution, the displacement potential is introduced here and displacement potential The concept is calculated as follows:

[0200]

[0201]

[0202] Then the weak form of the fluid-solid coupling equation is obtained. For the solid domain, the calculation formula is as follows:

[0203]

[0204] For the fluid domain, the calculation formula is as follows:

[0205]

[0206]

[0207] From the above formula, we can get a linear system composed of symmetric matrices, and its calculation formula is as follows:

[0208]

[0209]

[0210]

[0211]

[0212] in: is the fluid mass matrix.

[0213] Introducing impedance After that, the pressure on the boundary can be written as:

[0214]

[0215] in is the normal velocity of the fluid particles flowing out of the interface.

[0216] Therefore, substituting into the fluid motion equation, we can get the following calculation formula:

[0217]

[0218] Written in the form of virtual displacement, the calculation formula is as follows:

[0219]

[0220] Therefore, the impedance boundary is set, and its calculation formula is as follows:

[0221]

[0222] Finally, the fluid-solid coupling discrete equations with time format are obtained, which specifically include the following:

[0223] Consider the computational domain , where the wave equation is calculated as follows:

[0224]

[0225] In addition, the four boundary conditions of the fluid domain have been obtained in the third step. Then, based on the four boundary conditions, the two sides of the wave equation are multiplied by the imaginary displacement potential. , and integrate, and then apply Green's formula, we can get the following calculation formula:

[0226]

[0227] Will Substituting the expression of the pressure field on into the above formula, we get the following expression:

[0228]

[0229] Then we can get the following expression:

[0230]

[0231] in It's the border part of Boundary properties other than can be described using free surface boundaries and rigidly connected boundaries.

[0232] Put the above formula in After the integration, add the two equations together and you can write it in matrix form. The expression is as follows:

[0233]

[0234]

[0235] in: is the fluid displacement, is the structural displacement.

[0236] The equation is discretized in time, and finally the fluid-solid coupling discrete equation group that needs to be solved is obtained, and its specific expression is as follows:

[0237]

[0238] In this way, The fluid-solid coupling discrete equations that need to be solved at all times are composed of symmetric matrices, which can save calculation time to a great extent.

[0239] Furthermore, the displacement of the dam body and the dam foundation in the horizontal and vertical directions and the relevant data of the change over time can be solved. It can be seen from the maximum horizontal displacement cloud map of the gravity dam that the maximum horizontal displacement of the gravity dam during the earthquake occurs at the dam top, with a maximum value of 2.2 cm. At the same time, the maximum horizontal displacement calculated by the present invention is 27% smaller than the calculated value of the existing pseudo-static method, so the calculation result of the scheme of the present invention is more accurate.

[0240] Step 4: Post-process the displacement results according to the material properties, and derive the velocity response and acceleration response of the structure at each moment. The method for evaluating the impact of earthquake on structural stability is as follows:

[0241] Step 1: Evaluate the stability of the structural response of the gravity dam during an earthquake. According to the shear strength calculation formula in the "Specification for the Design of Concrete Gravity Dams" SL319-2018, calculate the anti-sliding stability and obtain the change in the stability coefficient at each moment. The time step is The minimum stability coefficient is 4.184, which appears in , indicating that the dam is able to resist sliding and is stable under earthquake excitation. The minimum stability coefficient given by the method not using the present invention is 3.76, which is smaller than the calculation result of the present invention, indicating that the fluid-solid coupling algorithm of the present invention solves the problem of overly conservative calculations before, and can obtain accurate results that are more in line with the actual situation and can better describe the real physical process.

[0242] Step 2: For concrete structures such as gravity dams, damage may be caused by a variety of causes, such as tensile cracking, compression failure, and stiffness degradation. Therefore, it is necessary to evaluate the stress level of the structure after dynamic time history analysis. Combined with the material properties, the stress tensor at each position is calculated. For the plane dam simulation object, the stress tensor at each node is There are 4 different non-zero components, and its matrix form can be expressed as follows:

[0243]

[0244] in: represents the normal stress, is the shear stress component. Based on the calculation results of the stress tensor, the principal stress results at each node are derived.

[0245] By solving the eigenvalue of the stress tensor composed of four stress components, the principal stress at each position of the gravity dam is obtained, and the principal stress is compared with the tensile and compressive strength of the concrete material to evaluate the risk of structural damage and ensure that the simulation results are accurate and reliable. The maximum compressive stress of the concrete calculated by the method of the present invention is 10.9MPa, which is less than the design value of the compressive strength of C30 concrete 14.3MPa, and the maximum tensile stress is 0.21MPa, which is less than the design value of the tensile strength of C30 concrete 1.43MPa; while the maximum compressive stress of the concrete obtained without using the method of the present invention is 12.9MPa, and the maximum tensile stress is 0.55MPa, both of which are greater than the calculation results of the present invention.

[0246] Therefore, the present invention fully considers the compressibility of reservoir water, dam foundation damping and the interaction between the gravity dam, dam foundation and reservoir water, introduces time domain deconvolution to calculate the seismic effect, and uses the applied load as the input, and then conducts modal analysis on the overall structure of the gravity dam, dam foundation and reservoir water, and performs dynamic time history calculation on the seismic water conditions suffered by the gravity dam. This can solve the problem of large frequency and structural response of the results in previous seismic calculations of gravity dams, and has important guiding significance for the structural seismic design, operation and maintenance, and reinforcement engineering of gravity dams.

[0247] An embodiment of a device applying the method of the present invention:

[0248] An electronic device comprising:

[0249] one or more processors;

[0250] A storage device for storing one or more programs;

[0251] When the one or more programs are executed by the one or more processors, the one or more processors implement the above-mentioned gravity dam seismic dynamics simulation method based on the time domain deconvolution algorithm.

[0252] A computer medium embodiment using the method of the present invention:

[0253] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the above-mentioned gravity dam seismic dynamics simulation method based on a time-domain deconvolution algorithm.

[0254] Those skilled in the art should understand that the embodiments of the present application can be provided as methods, systems, and computer program products. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, optical storage, etc.) containing computer-usable program codes.

[0255] The present application is described by flowcharts or / and block diagrams of the methods, devices (systems), and computer program products of the embodiments of the present application. It should be understood that each process or / and block in the flowchart or / and block diagram and the combination of the processes or / and blocks in the flowchart or / and block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart or / and block diagram. Figure 1 Process or multiple processes and / or boxes Figure 1A device that provides the functions specified in a block or multiple blocks.

[0256] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 Process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0257] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 Process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0258] The model in this application is an object that objectively describes the morphological structure with the help of physical or virtual expressions. The object is not equal to the object and is not limited to physical and virtual. It can be a data processing function, software program, processing mode, usage method, operation method, workflow, application process, electronic hardware, circuit module, processing system and system simulation.

[0259] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field can still modify or replace the specific implementation methods of the present invention with equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be included in the scope of protection of the claims of the present invention.

Claims

1. A gravity dam seismic dynamics simulation method based on a time domain deconvolution algorithm is characterized by: Includes the following: By generating a simulation model of the dam body created in advance, a dam body simulation object for earthquake-resistant dynamics simulation is established according to the gravity dam design data and geological information, which at least includes the gravity dam, reservoir water and dam foundation; Using the previously created dam body modal analysis model, taking into account the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption effect of the reservoir water boundary, the modal analysis of the dam body simulation object including the gravity dam, reservoir water and dam foundation is performed to calculate the natural frequency and structural modal vibration shape of the dam body simulation object; Based on the earthquake action calculation model created in advance and the geological information around the gravity dam, the earthquake acceleration time history is determined, and the earthquake action input data is calculated through the time domain deconvolution algorithm; According to the fluid-solid coupling response calculation model created in advance, the earthquake input data is applied to the dam simulation object, and the dynamic time history analysis is performed according to the natural frequency and structural modal vibration shape of the dam simulation object to obtain the seismic response data of the dam simulation object, and complete the seismic dynamic simulation of the gravity dam based on the time domain deconvolution algorithm; The method to obtain the seismic response data of the dam simulation object is as follows: Apply earthquake input data to the dam simulation object, and perform dynamic time history analysis based on the natural frequency and structural modal vibration shape of the dam simulation object to obtain the displacement field of the dam simulation object at several moments; Based on the displacement field, calculate the earthquake response acceleration, velocity, displacement and stress data of the dam simulation object at different times; According to the earthquake response acceleration, velocity, displacement and stress data, the impact of earthquake on the structural stability of the dam simulation object is evaluated to obtain the seismic response data.

2. The method for seismic dynamics simulation of a gravity dam based on a time domain deconvolution algorithm according to claim 1, characterized in that: The method of establishing a dam body simulation object for seismic dynamics simulation by generating a simulation model of the dam body created in advance and according to the gravity dam design data and geological information is as follows: Based on the gravity dam design engineering drawings and geological survey reports, establish two-dimensional plane simulation objects of the gravity dam, reservoir water and dam foundation; Then the two-dimensional plane simulation object is divided into unit grids to obtain the first dam body grid; Based on the location information of the boundary between reservoir water, gravity dam and dam foundation, set the number of encryption layers and the size of the encrypted area grid; According to the number of encryption layers and the mesh size of the encryption area, the first dam body mesh is encrypted to obtain the second dam body mesh; Based on the interface information between reservoir water and bedrock, additional grid nodes are set on the basis of the second dam body grid to characterize the properties of the reflected wave, thereby obtaining the third dam body grid; Based on the third dam body grid, the interface group, free boundary group and absorbing boundary group of the gravity dam are coupled to distinguish the propagation characteristics of seismic waves, so as to obtain the dam body simulation object for fluid and structure coupling calculation.

3. The method for seismic dynamics simulation of a gravity dam based on a time domain deconvolution algorithm according to claim 1, characterized in that: Using the previously created dam body modal analysis model, considering the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption effect of the reservoir water boundary, the modal analysis of the dam body simulation object including the gravity dam, reservoir water and dam foundation is performed. The method for calculating the natural frequency and structural modal vibration shape of the dam body simulation object is as follows: Based on the structural material properties of the gravity dam, set the structural density, Young's modulus and Poisson's ratio of the dam simulation object; According to the structural density, Young's modulus and Poisson's ratio, the mass matrix and stiffness matrix of the dam simulation object are established; Calculate the hydrodynamic pressure field of the fluid domain based on the medium characteristics and wave impedance data of wave propagation; Based on the hydrodynamic pressure field of the fluid and the wave transmission and reflection characteristics of the interface, and taking into account the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water and the absorption effect of the reservoir water boundary, the boundary conditions and constraints of the bottom of the dam foundation and the upstream and downstream surfaces are determined; According to the boundary conditions and constraints, the damping matrix and constraint excitation vector of the dam simulation object are established; According to the mass matrix, stiffness matrix, damping matrix and constraint excitation vector of the dam simulation object, the linear system of the gravity dam is established; The modal solution of the gravity dam linear system is performed to obtain the natural frequency and structural mode shape of the dam simulation object.

4. The method for seismic dynamics simulation of a gravity dam based on a time domain deconvolution algorithm according to claim 3, characterized in that: According to the medium characteristics and wave impedance data of wave propagation, the method for calculating the fluid hydrodynamic pressure field in the fluid domain is as follows: The Helmholtz equation is established based on the medium characteristics and wave impedance of wave propagation; Solve the Helmholtz equation to obtain the fluid wave vector; Based on the fluid characteristics of reservoir water being uniform, compressible, irrotational and inviscid, a wave equation satisfied by water pressure is established; According to the wave equation, construct the propagation expression of the fluid domain; Substitute the fluid wave vector into the propagation expression to obtain the water pressure data in the fluid domain; The water pressure data is subjected to inverse Fourier transform to obtain the fluid dynamic water pressure field in the time domain.

5. The method for seismic dynamics simulation of a gravity dam based on a time domain deconvolution algorithm according to claim 3, characterized in that: Based on the hydrodynamic pressure field of the fluid and the wave transmission and reflection characteristics of the interface, and considering the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption effect of the reservoir water boundary, the method for determining the boundary conditions and constraints of the bottom of the dam foundation and the upstream and downstream surfaces is as follows: Based on the hydrodynamic pressure field of the fluid and the wave transmission and reflection characteristics of the interface, the first node boundary condition on the interface between the dam body and the reservoir water is set; Considering the mass and damping of the dam foundation and reservoir water, the compressibility of the reservoir water, and the absorption of the reservoir water boundary, the influence of the flexible bedrock or rigid bedrock on the boundary conditions of the first node is obtained; Based on the impact information, the reservoir water absorption coefficient is set to represent the pressure gradient; On the free surface of the reservoir water, ignoring the surface wave effect of the reservoir water, set the second node boundary condition; At the upstream position of the reservoir water far from the dam surface, that is, at 2.5 times the gravity dam elevation, a stationary boundary is imposed to obtain the third node boundary condition; The first node boundary condition, the second node boundary condition, the third node boundary condition and the reservoir water absorption coefficient are coupled to obtain the boundary conditions and constraints of the bottom and upstream and downstream surfaces of the dam foundation.

6. The method for seismic dynamics simulation of a gravity dam based on a time domain deconvolution algorithm according to claim 1, characterized in that: Based on the earthquake action calculation model created in advance, the earthquake acceleration time history is determined according to the geological information around the gravity dam, and the method of calculating the earthquake action input data through the time domain deconvolution algorithm is as follows: According to the earthquake information of the actual geographical location of the gravity dam, the free-field earthquake acceleration time history function required for the dynamic time history analysis is determined, which includes the free-field horizontal earthquake acceleration time history function and the free-field vertical earthquake acceleration time history function; The free-field horizontal seismic acceleration time history function and the free-field vertical seismic acceleration time history function are integrated by the time-domain deconvolution algorithm to obtain the free-field horizontal seismic velocity time history function and the free-field vertical seismic velocity time history function; According to the free-field horizontal seismic velocity time history function and the free-field vertical seismic velocity time history function, the horizontal velocity time history and the vertical velocity time history of the earthquake's horizontal wave and vertical wave at the bottom of the dam foundation are calculated; Based on the horizontal velocity time history, vertical velocity time history, equivalent shear and compression wave velocities in the node force layer, and considering the viscous absorption characteristics of the dam foundation boundary, the node force calculation equations are constructed; Based on the node force calculation equations, the node forces in the vertical direction and the horizontal direction caused by the incident seismic wave are calculated, that is, the seismic action input data.

7. The method for seismic dynamics simulation of a gravity dam based on a time domain deconvolution algorithm according to claim 1, characterized in that: The method to obtain the displacement field of the dam simulation object at several moments is as follows: According to the propagation characteristics of the interface position of the dam simulation object, the displacement potential and displacement potential are established; Based on the displacement potential and displacement potential, as well as the natural frequency and structural mode vibration shape, the weak form of the fluid-solid coupling equation is obtained, which includes the solid domain calculation formula and the fluid domain calculation formula; Using the solid domain calculation formula and the fluid domain calculation formula, a linear system consisting of symmetric matrices is obtained; On the basis of the linear system, impedance boundaries are set to obtain a set of fluid-solid coupling discrete equations in time format; Decompose the earthquake action input data into the force on the dam foundation close to the upstream surface, the force on the bottom of the dam foundation and the force on the dam foundation close to the downstream surface; The structural state at the initial moment is obtained, the calculation time step is selected, and the force of the dam foundation close to the upstream surface, the force at the bottom of the dam foundation, and the force close to the downstream surface of the dam foundation are input into the fluid-solid coupling discrete equation group for solution to obtain the displacement field of the dam simulation object at each moment.

8. The method for seismic dynamics simulation of a gravity dam based on a time domain deconvolution algorithm according to claim 1, characterized in that: According to the earthquake response acceleration, velocity, displacement and stress data, the impact of earthquake on the structural stability of the dam simulation object is evaluated and the method for obtaining the seismic response data is as follows: Calculate the stability coefficient of the dam simulation object based on the earthquake response acceleration, velocity, displacement and stress data; Based on the stress data, the calculation results of the stress tensor of each node of the dam simulation object are obtained; According to the calculation results of the stress tensor, three stress results are obtained at each node position, including horizontal stress component, vertical stress component and shear stress component; Construct a stress tensor matrix according to horizontal stress components, vertical stress components and shear stress components; Calculate the stress tensor matrix to obtain the eigenvalue of the stress tensor and the principal stress at each position of the dam simulation object; Compare the principal stress with the tensile strength and compressive strength of the concrete material to obtain the comparison results; The stability coefficients and comparison results are summarized to obtain seismic response data, which is used to evaluate the impact of earthquakes on the structural stability of the dam simulation object.

9. The gravity dam seismic dynamics simulation system based on time domain deconvolution algorithm is characterized by: It includes: one or more processors; A storage device for storing one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors implement the seismic dynamics simulation method of a gravity dam based on a time domain deconvolution algorithm as described in any one of claims 1 to 8.

Images