CN116626752B - Ground vibration rotation component solving method based on field surface deformation rate - Google Patents
Ground vibration rotation component solving method based on field surface deformation rate Download PDFInfo
- Publication number
- CN116626752B CN116626752B CN202310674306.4A CN202310674306A CN116626752B CN 116626752 B CN116626752 B CN 116626752B CN 202310674306 A CN202310674306 A CN 202310674306A CN 116626752 B CN116626752 B CN 116626752B
- Authority
- CN
- China
- Prior art keywords
- field
- displacement
- time course
- translational
- simulation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 12
- 238000006073 displacement reaction Methods 0.000 claims abstract description 81
- 238000004088 simulation Methods 0.000 claims abstract description 43
- 230000001133 acceleration Effects 0.000 claims abstract description 33
- 239000013598 vector Substances 0.000 claims abstract description 20
- 238000012417 linear regression Methods 0.000 claims description 4
- 238000000611 regression analysis Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 2
- 230000004069 differentiation Effects 0.000 abstract 1
- 230000010354 integration Effects 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/24—Recording seismic data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention belongs to the field of ground vibration simulation in civil engineering disciplines, and provides a ground vibration rotation component solving method based on a ground surface deformation rate. Determining position coordinate information of each simulation point in the field; acquiring known multipoint translational seismic data of the surface of the site; combining the position coordinates of each simulation point with the earthquake motion translational displacement data to obtain a translational displacement field of the field, and determining an equation expression of a three-dimensional fitting curved surface of the displacement field at each moment; according to the displacement field curved surface equation, the normal vector corresponding to the displacement field curved surface tangent plane at each simulation point at any adjacent moment is obtained, and projected into three planes, and the included angle of the two normal vectors in the corresponding projection planes is obtained; the ratio of the included angle of the normal vectors of two adjacent moments in the projection plane to the corresponding time interval is used as the deformation rate of the field in the corresponding rotation direction, the rotation component speed time course of the earthquake motion is obtained, and the rotation component acceleration time course and the displacement time course at each simulation point of the field are obtained through differentiation and integration processing.
Description
Technical Field
The invention relates to the field of ground vibration simulation in civil engineering disciplines, in particular to a ground vibration rotation component solving method based on a ground surface deformation rate.
Background
China belongs to one of the most frequent countries of the world, and the earthquake disasters historically occur cause huge casualties and economic losses. Reasonable seismic inputs are the basis for accurately calculating the seismic response of an engineering structure, and thus, related studies on the field of seismic operations have been receiving extensive attention from students. During the seismic process, the seismic waves propagate from the source to the surface of the ground, and the surface generates three translational components and three rotational components (two rocking components and one torsion component) at the same time. In the past practical earthquake damage of the earthquake, examples of great deformation and serious damage to civil engineering structures caused by various earthquake vibration rotation components, such as torsion of a monument tower body, damage of high-rise buildings and the like, occur. These actual seismic disaster results indicate that the impact of the seismic rotational component on the structure is not negligible; in particular, for some swing and torsional matrix type structures, the seismic rotation component plays a decisive role in the seismic response of the structure.
Because the precision level of the seismic rotation observation instrument is low, a great deal of research is concentrated on the translational component of the seismic vibration, and the research on the rotational component is ignored once. At the end of the 19 th century, with the development of sensor technology, rotating earthquake motion acquisition equipment such as optical fibers, ring laser gyroscopes and the like are continuously developed, and rotating earthquake motion data can be obtained by observing and recording through a dense array. However, since the rotation component observation technique is not mature, the application is not wide enough, and the limitation of the observation field cannot be directly applied to engineering analysis. Some researchers have begun to try to derive the rotational component by theoretical derivation of the translational component of the seismic motion, most of which are based on linear elastic theory, using a number of assumptions, such as regarding the ground as a rigid body plane, and reacting the ground to linear elasticity, and thus the rotational component of the seismic motion obtained by these methods lacks accuracy.
At present, since the observation technology of the earthquake motion rotating station is not mature, the actual measurement record of the earthquake motion rotating component is very deficient, and the main stream calculation method has a plurality of defects due to the use of a large number of assumptions. Therefore, research on a calculation method of the ground vibration rotation component of the field surface is developed, and the research is very important for reasonably and accurately calculating the dynamic response of the structure under the action of the ground vibration.
Disclosure of Invention
The invention aims to provide a ground vibration rotation component solving method based on the field surface deformation rate, which provides a ground vibration input foundation for scientific and reasonable earthquake resistance analysis and design of a large building structure.
The technical scheme of the invention is as follows: firstly, determining position coordinate information at each simulation point in a field in a Cartesian coordinate system; secondly, collecting and processing the known multipoint translational earthquake motion time course data on the surface of the field; then, combining the position coordinates and translational displacement data of each simulation point in the field to obtain translational displacement fields of the field at each moment in the seismic event, and analyzing and determining a reasonable deformation form equation of the displacement field curved surface at each moment aiming at the specific field; then, according to a field deformation equation, obtaining normal vectors corresponding to tangential planes of displacement fields at any adjacent moment in the earthquake duration of each simulation point of the field, respectively projecting the obtained normal vectors into the XOZ, YOZ and XOY planes, and obtaining the included angles of the two normal vectors in the corresponding projection planes; finally, taking the ratio of the included angles in the three projection planes to the time difference of the corresponding moments as the change rate of the field in the corresponding direction, and regarding the deformation rate of the field as the speed representation of the corresponding vibration rotation component, and repeatedly performing the calculation on every two moments in the earthquake time course to obtain the speed time course of the vibration rotation component; differentiating and integrating the velocity schedule can obtain acceleration schedules and displacement schedules of three rotational components (two wobble components, one torsion component) at each simulation point of the field.
A method for solving the ground vibration rotation component based on the field surface deformation rate comprises the following steps:
s1, determining position coordinates of each simulation point on the surface of a field in a Cartesian coordinate system;
s2, collecting known multipoint translational earthquake motion time course data of the surface of the treatment site;
s3, at any moment in the whole time course of the earthquake motion, translational earthquake motion data (translational earthquake motion time course can be derived from actual observation data of the dense array) at all simulation points of the surface of the field and position coordinates are combined to obtain translational displacement field data of the field, and regression analysis is carried out on the displacement field data to determine a fitting form of a displacement field curved surface;
s4, determining a displacement field curved surface Z at any moment ti of earthquake motion in a multiple linear regression mode ti The (X, Y) equation is:
Z ti (X,Y)＝f ti (X,Y)
wherein f ti The specific form of (X, Y) is determined by multiple linear regression from the translational displacement field of the surface of the field;
the equation of the deformation form of the displacement field is converted into a mathematical three-dimensional form, and X, Y and Z are expressed by three principal axis coordinates:
F ti (X,Y,Z)＝f ti (X,Y)-Z
the position coordinates at any one of the simulation points in the field are known as (X 0 ,Y 0 ,Z 0 ) Solving the normal vector corresponding to the simulation point at the tangential plane of the displacement field curved surface through a field displacement field deformation equationThe expression is as follows:
wherein F is ti,X (X 0 ,Y 0 ,Z 0 )，F ti,Y (X 0 ,Y 0 ,Z 0 )，F ti,Z (X 0 ,Y 0 ,Z 0 ) Respectively F ti (X, Y, Z) at point (X 0 ,Y 0 ,Z 0 ) Derivatives with respect to X, Y, Z;
at the earthquake moment ti+deltat, aiming at the data of the displacement field of the surface of the field, the curved surface deformation equation of the displacement field at the moment ti+deltat is re-calculated as follows:
F ti+Δt (X,Y,Z)＝f ti+Δt (X,Y)-Z
obtaining the point (X) from the displacement field equation 0 ,Y 0 ,Z 0 ) Tangential plane at curved surface of displacement field at time ti+delta t and corresponding normal vectorThe expression is as follows:
wherein F is ti+Δt,X (X 0 ,Y 0 ,Z 0 )，F ti+Δt,Y (X 0 ,Y 0 ,Z 0 )，F ti+Δt,Z (X 0 ,Y 0 ,Z 0 ) Respectively F ti+Δt (X, Y, Z) at point (X 0 ,Y 0 ,Z 0 ) Derivatives with respect to X, Y, Z;
s5, the normal vectors corresponding to the two earthquake moments are obtainedAnd normal vector->Projected into the planes XOZ, YOZ, XOY in the Cartesian coordinate system, respectively, to find +.>And->Included angles in the corresponding projection planes are:
will beThe ratio of the included angle in the corresponding projection plane to the time difference delta t is taken as the deformation rate of each simulation point in the field in the X, Y and Z directions, and the change of the angular displacement of the field surface at adjacent moments represents the action of the ground vibration rotation component and comprises two swinging components and one torsion component; the deformation rate of the field in the X, Y and Z directions is used as the speed representation of the corresponding rotation component, namely:
repeating the steps in the earthquake motion time history, and obtaining the speed time courses of the earthquake motion rotating components at different positions of the surface of the field by solving the deformation rate of each simulation point of the field in each section of continuous time interval;
s6, differentiating and integrating the obtained seismic rotation component speed time courses at different simulation points in a time domain respectively to obtain acceleration time courses and displacement time courses of the seismic rotation components at each simulation point.
According to the whole flow, the earthquake motion component (two swinging directions and one torsion direction) at each simulation point can be obtained by solving the translational motion components (X, Y and Z directions) at a plurality of simulation points on the surface of the field.
The invention has the beneficial effects that: the earthquake motion time course of the surface rotation component of the site required by the dynamic response analysis of the civil engineering structure under the earthquake action is obtained through a manual solving method, and an important load input foundation is provided for the multidimensional and multipoint earthquake-resistant analysis and design of the large-scale structure.
Drawings
FIG. 1 is a basic logic idea of the proposed method of the present invention;
FIG. 2 is a schematic diagram of solving normal vectors corresponding to simulation points at the tangent plane of a displacement field curved surface at two adjacent moments according to the present invention; FIG. 2 (a) is a schematic diagram showing the normal vector at the tangent plane of the curved surface of the displacement field at the earthquake motion time ti; FIG. 2 (b) is a schematic diagram of the corresponding normal vector at the tangent plane of the curved surface of the displacement field at the earthquake motion time ti+Δt;
FIG. 3 is a schematic diagram of the deformation rate calculation according to the present invention, wherein the normal vector of the simulation point is projected into each plane to calculate the included angle;
FIG. 4 is a schematic view of planar positions of target simulation points of an array site according to an embodiment of the present invention;
FIG. 5 is a graph of known three-way translational seismic acceleration time courses at the surfaces of four station sites, line1-a1, line6-a2, line11-a3, line16-a4, in an array site in accordance with an embodiment of the present invention; FIG. 5 (a) is a known x-direction translational seismic acceleration time course at the surface of a line1-a1 station site; FIG. 5 (b) is a known y-translational seismic acceleration time course at the surface of the line1-a1 station site; FIG. 5 (c) is a known z-translational seismic acceleration time course at the surface of a line1-a1 station site; FIG. 5 (d) is a known x-direction translational seismic acceleration time course at the surface of a line6-a2 station site; FIG. 5 (e) is a known y-translational seismic acceleration time course at the surface of a line6-a2 station site; FIG. 5 (f) is a known z-translational seismic acceleration time course at the surface of a line6-a2 station site; FIG. 5 (g) is a known x-direction translational seismic acceleration time course at the surface of the line11-a3 station site; FIG. 5 (h) is a known y-translational seismic acceleration time course at the surface of the line11-a3 station site; FIG. 5 (i) is a known z-translational seismic acceleration time course at the surface of a line11-a3 station site; FIG. 5 (j) is a known x-direction translational seismic acceleration time course at the surface of the line16-a4 station site; FIG. 5 (k) is a known y-translational seismic acceleration time course at the surface of the line16-a4 station site; FIG. 5 (l) is a known z-translational seismic acceleration time course at the surface of a line16-a4 station site;
FIG. 6 is a known three-way translational seismic displacement schedule at the surfaces of four station sites, line1-a1, line6-a2, line11-a3, line16-a4, in an array site in accordance with an embodiment of the present invention; FIG. 6 (a) is a known x-direction translational seismic displacement time course at the surface of a line1-a1 station site; FIG. 6 (b) is a known y-translational seismic displacement time course at the surface of a line1-a1 station site; FIG. 6 (c) is a known z-translational seismic displacement time course at the surface of a line1-a1 station site; FIG. 6 (d) is a known x-direction translational seismic displacement time course at the surface of a line6-a2 station site; FIG. 6 (e) is a known y-translational seismic displacement time course at the surface of a line6-a2 station site; FIG. 6 (f) is a known z-translational seismic displacement time course at the surface of a line6-a2 station site; FIG. 6 (g) is a known x-direction translational seismic displacement time course at the surface of the line11-a3 station site; FIG. 6 (h) is a known y-translational seismic displacement time course at the surface of the line11-a3 station site; FIG. 6 (i) is a known z-translational seismic displacement time course at the surface of a line11-a3 station site; FIG. 6 (j) is a known x-direction translational seismic displacement time course at the surface of the line16-a4 station site; FIG. 6 (k) is a known y-translational seismic displacement time course at the surface of the line16-a4 station site; FIG. 6 (l) is a known z-translational seismic displacement time course at the surface of the line16-a4 station site;
FIG. 7 is a three-way rotational seismic acceleration time course obtained by solving at the surfaces of four station sites, line1-a1, line6-a2, line11-a3, line16-a4, in an array site according to an embodiment of the present invention; FIG. 7 (a) is an x-direction rotational seismic acceleration time course solved at the line1-a1 station site surface; FIG. 7 (b) is a y-direction rotational seismic acceleration time course solved at the line1-a1 station site surface; FIG. 7 (c) is a z-direction rotational seismic acceleration time course solved at the line1-a1 station site surface; FIG. 7 (d) is the x-direction rotational seismic acceleration time course solved at the line6-a2 station site surface; FIG. 7 (e) is a y-direction rotational seismic acceleration time course solved at the line6-a2 station site surface; FIG. 7 (f) is a z-direction rotational seismic acceleration time course solved at the line6-a2 station site surface; FIG. 7 (g) is the x-direction rotational seismic acceleration time course solved at the line11-a3 station site surface; FIG. 7 (h) is the y-direction rotational seismic acceleration time course solved at the line11-a3 station site surface; FIG. 7 (i) is a z-direction rotational seismic acceleration time course solved at the line11-a3 station site surface; FIG. 7 (j) is the x-direction rotational seismic acceleration time course solved at the line16-a4 station site surface; FIG. 7 (k) is the y-direction rotational seismic acceleration time course solved at the line16-a4 station site surface; FIG. 7 (l) is a z-direction rotational seismic acceleration time course solved at the line16-a4 station site surface;
FIG. 8 is a three-way rotational seismic displacement time course obtained by solving at the surfaces of four station sites, line1-a1, line6-a2, line11-a3, line16-a4, in an array site according to an embodiment of the present invention; FIG. 8 (a) is an x-direction rotational seismic displacement time course solved at the line1-a1 station site surface; FIG. 8 (b) is a y-direction rotational seismic displacement time course solved at the line1-a1 station site surface; FIG. 8 (c) is a z-direction rotational seismic displacement time course solved at the line1-a1 station site surface; FIG. 8 (d) is an x-direction rotational seismic displacement time course solved at the line6-a2 station site surface; FIG. 8 (e) is a y-direction rotational seismic displacement time course solved at the line6-a2 station site surface; FIG. 8 (f) is a z-direction rotational seismic displacement time course solved at the line6-a2 station site surface; FIG. 8 (g) is an x-direction rotational seismic displacement time course solved at the line11-a3 station site surface; FIG. 8 (h) is a y-direction rotational seismic displacement time course solved at the line11-a3 station site surface; FIG. 8 (i) is a z-direction rotational seismic displacement time course solved at the line11-a3 station site surface; FIG. 8 (j) is an x-direction rotational seismic displacement time course solved at the line16-a4 station site surface; FIG. 8 (k) is a y-direction rotational seismic displacement time course solved at the line16-a4 station site surface; FIG. 8 (l) is a time course of the resolved z-direction rotational seismic displacement at the line16-a4 station site surface.
Detailed Description
The following describes the embodiments of the present invention further with reference to the drawings and technical schemes. The example of the array site shown in fig. 4 is selected for specific illustration and description, all the stations in the site are selected as target simulation points, and the time course of the translational components of the multipoint earthquake motion at each station is known, so as to simulate the time course of the rotational earthquake motion at all the stations of the site. The basic idea of the invention is shown in fig. 1, and the specific embodiment is as follows:
(1) The radius of the array field is R from inside to outside 1 ＝30m、R 2 ＝35m、R 3 ＝45m、R 4 ＝60m、R 5 In total, 101 station simulation points are arranged in the field at 80m, are distributed in 20 directions of the field, and are sequentially line1 to line20, and 5 station simulation points are arranged in each direction, and are sequentially a1 to a5. The site center station line1-a0 is taken as the origin of coordinates, so that the position coordinates of the simulation points of each station can be determined.
(2) The known multi-point translational earthquake motion data at each simulation point of the surface of the array site is analyzed and processed, translational displacement time-course data at each station of the surface of the site is collected, the translational earthquake motion sampling frequency is 100Hz, the time interval of the earthquake motion is 0.01s, and the duration is 20s in the embodiment. The known three-way translational earthquake motion acceleration time course and displacement time course at four stations of line1-a1, line6-a2, line11-a3 and line16-a4 are selected as examples for illustration, and are shown in fig. 5 and 6.
(3) The collected earthquake motion translation displacement data and position coordinates of each simulation point of the surface of the site are combined to establish a three-dimensional translation displacement field of the site where the array is located, and the three-dimensional curved surface form of the displacement field of the array site can be determined by carrying out regression analysis on the curved surface of the displacement field, so that the three-dimensional curved surface form of the displacement field of the array site can be better simulated by adopting a paraboloid.
(4) According to the field displacement field data at each moment in the earthquake motion time history, determining a specific expression of a displacement field parabolic equation at each moment, thereby obtaining normal vectors corresponding to tangential planes of a displacement field fitting curved surface at each simulation point at any adjacent moment in the earthquake motion, respectively projecting the normal vectors obtained at two moments into XOZ, YOZ and XOY planes, and obtaining an included angle of the two normal vectors in the corresponding projection planes.
(5) And taking the ratio of the included angle of the two normal vectors in the projection plane to the time interval as the rotational deformation rate of each simulation point of the field in the corresponding direction, and solving the rotational deformation rate of each simulation point of the surface of the field at all times of earthquake so as to obtain the speed time course of the earthquake motion component at all simulation points.
(6) And differentiating and integrating the obtained seismic vibration rotation component speed time course in the time domain, so as to further obtain the rotation component acceleration time course and displacement time course at each simulation point of the surface of the field. Because of the large number of stations in the field, the rotational seismic acceleration time course and displacement time course obtained at four simulation points of line1-a1, line6-a2, line11-a3, line16-a4 are selected as examples for illustration, as shown in fig. 7 and 8.
Claims (1)
1. The method for solving the ground vibration rotation component based on the field surface deformation rate is characterized by comprising the following steps of:
s1, determining position coordinates of each simulation point on the surface of a field in a Cartesian coordinate system;
s2, collecting known multipoint translational earthquake motion time course data of the surface of the treatment site;
s3, at any moment in the whole time course of the earthquake motion, combining translational earthquake motion displacement data and position coordinates of all simulation points on the surface of the field to obtain translational displacement field data of the field, and carrying out regression analysis on the displacement field data to determine a fitting form of a displacement field curved surface;
s4, determining a displacement field curved surface Z at any moment ti of earthquake motion in a multiple linear regression mode ti The (X, Y) equation is:
Z ti (X,Y)＝f ti (X,Y)
wherein f ti The specific form of (X, Y) is determined by multiple linear regression from the translational displacement field of the surface of the field;
the equation of the deformation form of the displacement field is converted into a mathematical three-dimensional form, and X, Y and Z are expressed by three principal axis coordinates:
F ti (X,Y,Z)＝f ti (X,Y)-Z
the position coordinates at any one of the simulation points in the field are known as (X 0 ,Y 0 ,Z 0 ) Solving the normal vector corresponding to the simulation point at the tangential plane of the displacement field curved surface through a field displacement field deformation equationThe expression is as follows:
wherein F is ti,X (X 0 ,Y 0 ,Z 0 )，F ti,Y (X 0 ,Y 0 ,Z 0 )，F ti,Z (X 0 ,Y 0 ,Z 0 ) Respectively F ti (X, Y, Z) at point (X 0 ,Y 0 ,Z 0 ) Derivatives with respect to X, Y, Z;
at the earthquake moment ti+deltat, aiming at the data of the displacement field of the surface of the field, the curved surface deformation equation of the displacement field at the moment ti+deltat is re-calculated as follows:
F ti+Δt (X,Y,Z)＝f ti+Δt (X,Y)-Z
obtaining the point (X) from the displacement field equation 0 ,Y 0 ,Z 0 ) Tangential plane at curved surface of displacement field at time ti+delta t and corresponding normal vectorThe expression is as follows:
wherein F is ti+Δt,X (X 0 ,Y 0 ,Z 0 )，F ti+Δt,Y (X 0 ,Y 0 ,Z 0 )，F ti+Δt,Z (X 0 ,Y 0 ,Z 0 ) Respectively F ti+Δt (X, Y, Z) at point (X 0 ,Y 0 ,Z 0 ) Derivatives with respect to X, Y, Z;
s5, the normal vectors corresponding to the two earthquake moments are obtainedAnd normal vector->Projected into the planes XOZ, YOZ, XOY in the Cartesian coordinate system, respectively, to find +.>And->Included angles in the corresponding projection planes are:
will beThe ratio of the included angle in the corresponding projection plane to the time difference delta t is taken as the deformation rate of each simulation point in the field in the X, Y and Z directions, and the change of the angular displacement of the field surface at adjacent moments represents the action of the ground vibration rotation component and comprises two swinging components and one torsion component; the deformation rate of the field in the X, Y and Z directions is taken as a pairThe speed of the component to be rotated is characterized by:
repeating the steps in the earthquake motion time history, and obtaining the speed time courses of the earthquake motion rotating components at different positions of the surface of the field by solving the deformation rate of each simulation point of the field in each section of continuous time interval;
s6, differentiating and integrating the obtained seismic rotation component speed time courses at different simulation points in a time domain respectively to obtain acceleration time courses and displacement time courses of the seismic rotation components at each simulation point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310674306.4A CN116626752B (en) | 2023-06-08 | 2023-06-08 | Ground vibration rotation component solving method based on field surface deformation rate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310674306.4A CN116626752B (en) | 2023-06-08 | 2023-06-08 | Ground vibration rotation component solving method based on field surface deformation rate |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116626752A CN116626752A (en) | 2023-08-22 |
| CN116626752B true CN116626752B (en) | 2023-10-24 |
Family
ID=87641772
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310674306.4A Active CN116626752B (en) | 2023-06-08 | 2023-06-08 | Ground vibration rotation component solving method based on field surface deformation rate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116626752B (en) |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000121743A (en) * | 1998-10-14 | 2000-04-28 | Osaka Gas Co Ltd | Evaluating method for seismic shake distribution |
| WO2014209141A2 (en) * | 2014-10-30 | 2014-12-31 | Instytut Technik Innowacyjnych Emag | Method and system for assessing a risk of high-energy earth bursts generated by underground mining |
| CN109683195A (en) * | 2019-01-09 | 2019-04-26 | 武汉理工大学 | A kind of earthquake motion rocking component acquisition methods |
| CN109931878A (en) * | 2018-07-13 | 2019-06-25 | 上海海事大学 | A kind of building curtain wall seismic deformation monitoring method based on digital speckle label |
| CN110389379A (en) * | 2019-07-12 | 2019-10-29 | 中国地震局地球物理研究所 | The near-fault ground motion Acceleration time course approximating method of ground permanent displacement can be characterized |
| CN111272366A (en) * | 2020-03-02 | 2020-06-12 | 东南大学 | Bridge displacement high-precision measurement method based on multi-sensor data fusion |
| CN111458744A (en) * | 2020-04-09 | 2020-07-28 | 西南交通大学 | Spatial rotation seismic oscillation simulation method |
| CN112698406A (en) * | 2020-12-22 | 2021-04-23 | 大连理工大学 | Method for simulating seawater velocity potential under seismic excitation of ocean field |
| CN113009566A (en) * | 2021-03-19 | 2021-06-22 | 中国地震局地球物理研究所 | In-situ seismic motion prediction model and construction method thereof |
| CN113484911A (en) * | 2021-04-30 | 2021-10-08 | 石家庄学院 | Earthquake displacement reaction calculation method based on earthquake motion evolution power spectrum |
| CN113608262A (en) * | 2021-08-11 | 2021-11-05 | 中国地质大学（北京） | Seismic data processing method and device for calculating rotation component by using translation component |
| CN115469360A (en) * | 2022-09-13 | 2022-12-13 | 中铁二院工程集团有限责任公司 | Cross-fault seismic motion field simulation method based on finite fault and equivalent velocity pulse |
| CN115630282A (en) * | 2022-11-14 | 2023-01-20 | 同济大学 | Method for effectively monitoring interlayer displacement of high-rise building based on earthquake background noise |
| JP2023012586A (en) * | 2021-07-14 | 2023-01-26 | 大成建設株式会社 | Method and system for estimating size of earthquake vibrations |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102787676B (en) * | 2011-05-20 | 2014-07-16 | 青岛理工大学 | Multi-objective performance-based aseismic design method of engineering structure |
-
2023
- 2023-06-08 CN CN202310674306.4A patent/CN116626752B/en active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000121743A (en) * | 1998-10-14 | 2000-04-28 | Osaka Gas Co Ltd | Evaluating method for seismic shake distribution |
| WO2014209141A2 (en) * | 2014-10-30 | 2014-12-31 | Instytut Technik Innowacyjnych Emag | Method and system for assessing a risk of high-energy earth bursts generated by underground mining |
| CN109931878A (en) * | 2018-07-13 | 2019-06-25 | 上海海事大学 | A kind of building curtain wall seismic deformation monitoring method based on digital speckle label |
| CN109683195A (en) * | 2019-01-09 | 2019-04-26 | 武汉理工大学 | A kind of earthquake motion rocking component acquisition methods |
| CN110389379A (en) * | 2019-07-12 | 2019-10-29 | 中国地震局地球物理研究所 | The near-fault ground motion Acceleration time course approximating method of ground permanent displacement can be characterized |
| CN111272366A (en) * | 2020-03-02 | 2020-06-12 | 东南大学 | Bridge displacement high-precision measurement method based on multi-sensor data fusion |
| CN111458744A (en) * | 2020-04-09 | 2020-07-28 | 西南交通大学 | Spatial rotation seismic oscillation simulation method |
| CN112698406A (en) * | 2020-12-22 | 2021-04-23 | 大连理工大学 | Method for simulating seawater velocity potential under seismic excitation of ocean field |
| CN113009566A (en) * | 2021-03-19 | 2021-06-22 | 中国地震局地球物理研究所 | In-situ seismic motion prediction model and construction method thereof |
| CN113484911A (en) * | 2021-04-30 | 2021-10-08 | 石家庄学院 | Earthquake displacement reaction calculation method based on earthquake motion evolution power spectrum |
| JP2023012586A (en) * | 2021-07-14 | 2023-01-26 | 大成建設株式会社 | Method and system for estimating size of earthquake vibrations |
| CN113608262A (en) * | 2021-08-11 | 2021-11-05 | 中国地质大学（北京） | Seismic data processing method and device for calculating rotation component by using translation component |
| CN115469360A (en) * | 2022-09-13 | 2022-12-13 | 中铁二院工程集团有限责任公司 | Cross-fault seismic motion field simulation method based on finite fault and equivalent velocity pulse |
| CN115630282A (en) * | 2022-11-14 | 2023-01-20 | 同济大学 | Method for effectively monitoring interlayer displacement of high-rise building based on earthquake background noise |
Non-Patent Citations (3)
| Title |
|---|
| Seismic response analyses of transmission towers under multidimensional ground motions with rocking and torsion components;Jie Zhang等;《Journal of Aerospace Engineering》;第33卷(第6期);全文 * |
| SMART-1台阵地震动激励下输电塔线体系反应分析;田利等;《工程力学》;第30卷;全文 * |
| 提取转动地震动的曲面拟合法及结构多维响应分析;张杰;《中国博士学位论文全文数据库工程科技Ⅱ辑》(2022年第03期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116626752A (en) | 2023-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lister et al. | Fabric development in shear zones: theoretical controls and observed phenomena | |
| CN110132515B (en) | Time-course-level iteration real-time hybrid test method based on model updating | |
| CN112632837B (en) | Method for determining longitudinal earthquake-resistant numerical value of underground structure | |
| CN102230785A (en) | Indoor 3D (3-dimensional) dimension measurement method | |
| CN105824987A (en) | Wind field characteristic statistical distributing model building method based on genetic algorithm | |
| CN105954796A (en) | Method and apparatus for determining hypocenter position of micro earthquake | |
| CN114004104A (en) | CORS site selection method based on checkerboard test | |
| CN116626752B (en) | Ground vibration rotation component solving method based on field surface deformation rate | |
| CN109164416A (en) | A kind of sound localization method of three planes, five yuan of microphone arrays | |
| CN107796643A (en) | A kind of model-free quick damnification recognition method theoretical based on statistical moment | |
| Liu et al. | Longitudinal integral response deformation method for the seismic analysis of a tunnel structure | |
| Ngeljaratan et al. | Novel digital image correlation instrumentation for large-scale shake table tests | |
| CN116299740A (en) | Space domain gravity multi-parameter analysis forward modeling method of rotary rectangular prism | |
| CN111462319B (en) | Three-dimensional numerical simulation method for basin-mountain coupling site vibration | |
| Ma et al. | Application and research of 3D laser scanning technology in steel structure installation and deformation monitoring | |
| McGinnis et al. | Application of multiple digital image correlation sensors in earthquake engineering | |
| CN103353286A (en) | Method for detecting coordinates of cross node of two cylindrical rod members | |
| CN105891780A (en) | Indoor scene positioning method and indoor scene positioning device based on ultrasonic array information | |
| CN218097857U (en) | Test structure for acquisition system | |
| Zhou et al. | A preliminary research establishing the present-time intraplate blocks movement model on the Chinese mainland based on GPS data | |
| Ding et al. | Study on the natural vibration characteristics of single-point diamond fly cutting machine tool based on transfer matrix method for multibody systems | |
| Han et al. | Application of phase three dimensional laser scanner in high altitude large volume irregular structure | |
| CN115096534B (en) | Compliance surface identification method based on non-reference point partition test | |
| CN114781432B (en) | Displacement resolving method based on multi-source information fusion and trend removal fluctuation analysis | |
| Huo et al. | A study of impact point detecting method based on seismic signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |