CN110389379A - The near-fault ground motion Acceleration time course approximating method of ground permanent displacement can be characterized - Google Patents
The near-fault ground motion Acceleration time course approximating method of ground permanent displacement can be characterized Download PDFInfo
- Publication number
- CN110389379A CN110389379A CN201910630179.1A CN201910630179A CN110389379A CN 110389379 A CN110389379 A CN 110389379A CN 201910630179 A CN201910630179 A CN 201910630179A CN 110389379 A CN110389379 A CN 110389379A
- Authority
- CN
- China
- Prior art keywords
- time
- histories
- acceleration
- displacement
- time course
- 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.)
- Granted
Links
- 230000001133 acceleration Effects 0.000 title claims abstract description 103
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000033001 locomotion Effects 0.000 title claims abstract description 19
- 238000012937 correction Methods 0.000 claims abstract description 36
- 230000004044 response Effects 0.000 claims abstract description 32
- 238000001228 spectrum Methods 0.000 claims abstract description 31
- 239000000284 extract Substances 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- 238000005316 response function Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 6
- 238000003325 tomography Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001308 synthesis method 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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/32—Transforming one recording into another or one representation into another
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
This application involves the near-fault ground motion Acceleration time course approximating methods that one kind can characterize ground permanent displacement, comprising the following steps: forms primary earthquake Acceleration time course, the adjustment of first time baseline is carried out to initially vibration, extracts second of baseline correction parameter of seismic acceleration time-histories, the seismic acceleration time-histories adjustment based on multiple target, to assessment of target component etc..Using the approximating method of the near-fault ground motion Acceleration time course that can characterize ground permanent displacement of the application, the seismic acceleration record that this method obtains, in addition to can satisfy acceleration response spectrum, peak accelerator and the characteristic for being able to reflect out Chi-chi earthquake velocity pulse, ground permanent displacement can be characterized simultaneously, the input data of the high seismic acceleration stroke of referential is capable of providing when carrying out Aseismic Design.
Description
Technical field
This application involves the near-fault ground motion Acceleration time course approximating methods that one kind can characterize ground permanent displacement, are applicable in
In the technical field of Aseismic Design.
Background technique
Ground permanent displacement refers to the ground deformation caused in seismic process due to fault movement.As illustrated in fig. 1 and 2
Seismic acceleration, speed and the displacement time-histories of Wenchuan earthquake station 051PXZ east-west direction, corresponding permanent displacement is about 56cm.
Many violent earthquake engineering earthquakes show that the permanent big displacement near field ground is the master that building or structures finally collapse in history
Reason is wanted, therefore is carrying out the long periods knots such as large-scale fluid storage tanks, ocean platform, tunnel, super high-rise building and Loads of Long-span Bridges
When the Aseismic Design of structure, it is necessary to pay close attention to the displacement of Seismic input, the case where especially close tomography is with the presence of permanent displacement.
For specific engineering ground, designer needs to obtain the Seismic input met certain condition, under normal circumstances
It is required to meet the relevant acceleration response spectrum in place, peak accelerator etc., to the closer site of detachment layer, it is also necessary to consider earthquake
Dynamic velocity pulse characteristic and ground permanent displacement.But it is difficult to find from real seismic record while expires in practical engineering application
The multiple targets of foot especially include the seismic acceleration time-histories of ground permanent displacement as Seismic input data, it is therefore desirable to are based on
Different targets develops corresponding earthquake synthesis method, to meet the needs of engineer application.
Summary of the invention
For the application using double baseline corrections and the method for adjustment in time domain, ground permanent displacement can be characterized by establishing one kind
Near-fault ground motion Acceleration time course approximating method, the seismic acceleration record that this method obtains, in addition to can satisfy accelerate
Response spectrum, peak accelerator and the characteristic for being able to reflect out Chi-chi earthquake velocity pulse are spent, while ground can be characterized forever
Displacement long.
This application involves the near-fault ground motion Acceleration time course approximating methods that one kind can characterize ground permanent displacement, including
Following steps:
The first step forms primary earthquake Acceleration time course
Primary earthquake Acceleration time course a is generated according to known aimed acceleration response spectrum0(t), or from actual seismic remember
It is selected in record and meets preferable accelerogram as initial time-histories a with aimed acceleration response spectrum0(t);The target accelerates
Spending response spectrum is that the statistical analysis based on acceleration transducer measurement result obtains;
Second step carries out the adjustment of first time baseline to initially vibration
High-pass filtering is carried out to initially vibration using filter function, by first time baseline Acceleration time course adjusted
It is defined as a01(t);
Third step extracts second of baseline correction parameter of seismic acceleration time-histories
To Acceleration time course a01(t) it is integrated, obtains corresponding speed time-histories v0(t) and displacement time-histories d0(t), then
Extract the parameter of second of baseline correction;
4th step, the seismic acceleration time-histories adjustment based on multiple target
It is frequency is ω that goal response, which is composed discrete,1, ω2..., ωMM single-degree-of-freedom system, to each frequency control
System point is successively adjusted;
After the first round adjustment for completing M control point, saving time-histories adjusted is a2(t), and acceleration peak value is adjusted
For target value;To Acceleration time course a2(t) it is integrated to obtain corresponding speed time-histories v2(t), line then is carried out to speed time-histories
Property fitting;
Acceleration time course is finally subtracted into offset, completes second of baseline correction and the control to permanent displacement;
5th step, the assessment to target component
Using the Acceleration time course a by adjusting above2(t) corresponding acceleration response spectrum, peak accelerator, speed are calculated
Degree and displacement time-histories, evaluation meets situation to preset target component, if being unsatisfactory for target component, with a2(t) conduct
New time-histories determines new baseline correction parameter since third step, repeats the above steps, and meets all targets until obtaining
Acceleration, speed and displacement time-histories.
Preferably, in third step, the extraction process of baseline correction parameter is as follows:
(1) start the time t of offset zero shift1
From displacement time-histories d0(t) terminal tendIt takes forward at equal intervals a little as t1, i.e. t1=tend- n Δ, Δ is between the time
Every;Then straight line d is usedn=d0,n-af,nT is fitted t1To tendThe displacement time-histories of this section, utilizes fitting parameter af,nIt is oblique to calculate time shift
Rate ratio kn=af,n+1+1/af,n, work as knOccurring the point jumpy corresponding time first is t1;
(2) the corresponding time t of peak acceleratorpga
The corresponding time t of peak acceleratorpgaIt is worth the corresponding time for Acceleration time course amplitude maximum absolute value;
(3) the initial time t of permanent displacement2
Definition does not correct displacement d0(t) time finally intersected with time shaft is td0, take t3=max [tpga,td0], then
Take t3Each moment to record end is t2, i.e. t2=t3+ n Δ t, Δ t are time interval, are fitted t using quadratic function2It arrives
tendPosition, move time-histories, and subtract offset from Acceleration time course and initial time-histories is corrected, and be displaced after calculating correction
Time-histories is in t2To tendThe flatness of periodMaximum determines t2Value,
R is t in formula2To tendThe linearly dependent coefficient of correction the displacement time-histories and its fitting a straight line of period, b is t2It arrives
tendThe least square regression straight slope of the correction displacement time-histories of period, σ are displacement time-histories t after correction2To tendPeriod
Variance.
Preferably, in the 4th step, method that each frequency control point is adjusted are as follows:
If control point frequencies omegakAcceleration time course before adjustment isIt is superimposed an increment time-histories Δ a (t), makes to calculate
Response spectrum Sa(ωk) composed with goal responseBetween difference meet required precision, i.e.,
Δ a (t)=Rh (tm-t)B(t,ωk)t≤tm (1)
In formula: tmAt the time of reacting generation for single-degree-of-freedom system maximum;B(t,ωk) be when-frequency envelope curve function;h(tm-
T) it is unit impulse response function, indicates are as follows:
Wherein: ζ is damping ratio;φkFor initial phase, guarantee the unit impulse function of superposition in tm
Place obtains peak value;
In formula (1), R is amplitude regulation coefficient, by control point frequencies omegakThe calculating response spectrum S at placea(ωk) and target
Response spectrumDifference DELTA S determine, it may be assumed that
Δ S=CR (3)
Wherein:
By calculating above, can be obtained corresponding to frequencies omegakAdjustment time-histories Δ a (t), by it is primary adjustedly
Shake time-historiesIt can be written as:
The Acceleration time course that will be calculated according to formula (5)As next frequency control point ωk+1Input when
Journey is adjusted according to formula (1)-(5) process.
Preferably, in the 4th step, integral a is carried out to Acceleration time course2(t) corresponding speed time-histories v is obtained2(t), then
Linear fit is carried out to speed time-histories, is obtained:
vf(t)=v0-af·t
Then Acceleration time course is obtained in time period t by following formula2-t1Correction offset, i.e.,
Finally by Acceleration time course in t2-t1Section subtracts offset am, complete second of baseline correction and to permanent displacement
Control.
Using the approximating method of the near-fault ground motion Acceleration time course that can characterize ground permanent displacement of the application, the party
The seismic acceleration record that method obtains, in addition to can satisfy acceleration response spectrum, peak accelerator and be able to reflect out nearly tomography
The characteristic of earthquake motion velocity pulse, while ground permanent displacement can be characterized, referential is capable of providing when carrying out Aseismic Design
The input data of high seismic acceleration stroke.
Detailed description of the invention
Fig. 1 shows seismic acceleration, speed and the position of the Wenchuan earthquake station 051PXZ east-west direction before baseline correction
Move the schematic diagram of time-histories.
Fig. 2 shows seismic acceleration, speed and the position of the Wenchuan earthquake station 051PXZ east-west direction after baseline correction
Move the schematic diagram of time-histories.
Fig. 3 shows the near-fault ground motion Acceleration time course approximating method that can characterize ground permanent displacement of the application
Flow diagram.
Fig. 4 shows the schematic diagram of goal response spectrum and fit solution in embodiments herein.
Fig. 5 is shown in embodiments herein according to the acceleration of this method fitting, speed and the signal for being displaced time-histories
Figure.
Specific embodiment
For the purposes, technical schemes and advantages of the application are more clearly understood, below in conjunction with attached drawing to the application
Embodiment be described in detail.It should be noted that in the absence of conflict, embodiment and implementation in the application
Feature in mode can mutual any combination.
As shown in figure 3, which show when the near-fault ground motion acceleration that can characterize ground permanent displacement of the application
The flow diagram of journey approximating method, specifically includes the following steps:
The first step forms primary earthquake Acceleration time course
Primary earthquake Acceleration time course a is generated according to known aimed acceleration response spectrum0(t), or from actual seismic remember
It is selected in record and meets preferable accelerogram as initial time-histories a with aimed acceleration response spectrum0(t);The target accelerates
Spending response spectrum is that the statistical analysis based on acceleration transducer measurement result obtains, and acceleration transducer is set to target area
It is interior;
Second step carries out the adjustment of first time baseline to initially vibration
High-pass filtering is carried out to initially vibration using such as Butterworth quadravalence filter function, cutoff frequency can be with
It is set as such as 0.05Hz, is defined as a by first time baseline Acceleration time course adjusted01(t)。
Third step extracts second of baseline correction parameter of seismic acceleration time-histories
To Acceleration time course a01(t) it is integrated, obtains corresponding speed time-histories v0(t) and displacement time-histories d0(t), then
The parameter of second of baseline correction is extracted, to illustrate the extraction process of following baseline correction parameter shown in Fig. 1 and 2.
(1) start the time t of offset zero shift1
From displacement time-histories d0(t) terminal tendIt takes forward at equal intervals a little as t1, i.e. t1=tend- n Δ, then uses straight line
dn=d0,n-af,nT is fitted t1To tendThe displacement time-histories of this section, time interval Δ can be taken as 0.25s, utilize fitting parameter af,n
Calculate time shift slope ratio kn, i.e. kn=af,n+1+1/af,n.Work as knOccurring the point jumpy corresponding time first is t1, Fig. 1
The time t of accelerogram determining in this way is given with 21。
(2) the corresponding time t of peak acceleratorpga
The corresponding time t of peak acceleratorpgaIt is worth the corresponding time for Acceleration time course amplitude maximum absolute value.
(3) the initial time t of permanent displacement2
Definition does not correct displacement d0(t) time finally intersected with time shaft is td0, take t3=max [tpga,td0], then
Take t3Each moment to record end is t2, i.e. t2=t3+ n Δ t (time interval Δ t=0.01s), utilizes quadratic function
It is fitted t2To tendDisplacement time-histories, and subtract from Acceleration time course offset and be corrected that (specific method can to initial time-histories
Referring to following formula 6 and 7), and time-histories is displaced in t after calculating correction2To tendThe flatness of period is maximum, i.e.,Come true
Determine t2Value, r is t in formula2To tendThe linearly dependent coefficient of correction the displacement time-histories and its fitting a straight line of period, b is t2
To tendThe least square regression straight slope of the correction displacement time-histories of period, σ are displacement time-histories t after correction2To tendTime
The variance of section.
4th step, the seismic acceleration time-histories adjustment based on multiple target
(1) composing goal response discrete be frequency is ω1, ω2..., ωMM single-degree-of-freedom system;
(2) each frequency control point is successively adjusted:
If setting control point frequencies omegakAcceleration time course before adjustment isSo to make to calculate response spectrum Sa(ωk) with
Goal response spectrumBetween difference meet required precision, a stackable increment time-histories Δ a (t), i.e.,
Δ a (t)=Rh (tm-t)B(t,ωk)t≤tm (1)
In formula: tmAt the time of reacting generation for single-degree-of-freedom system maximum;B(t,ωk) be when-frequency envelope curve function;h(tm-
T) it is unit impulse response function, can indicates are as follows:
Wherein: ζ is damping ratio;φkFor initial phase, guarantee the unit impulse function of superposition in tm
Place obtains peak value.
In formula (1), R is amplitude regulation coefficient, can be by control point frequencies omegakThe calculating response spectrum S at placea(ωk) and mesh
Mark response spectrumDifference DELTA S determine, it may be assumed that
Δ S=CR (3)
Wherein:
By calculating above, can be obtained corresponding to frequencies omegakAdjustment time-histories Δ a (t), then after once adjusting
Seismic Time-historyIt can be written as:
(3) Acceleration time course that will be calculated according to formula (5)As next frequency control point ωk+1Input
Time-histories is adjusted according to formula (1)-(5) process;
(4) after the first round adjustment for completing M control point, saving time-histories adjusted is a2(t), and acceleration peak is adjusted
Value is target value;
(5) to Acceleration time course a2(t) it is integrated to obtain corresponding speed time-histories v2(t), then speed time-histories is carried out
Linear fit obtains:
vf(t)=v0-af·t (6)
Then Acceleration time course is obtained in time period t by following formula2-t1Correction offset, i.e.,
Finally by Acceleration time course in t2-t1Section subtracts offset am, complete second of baseline correction and to permanent displacement
Control.
5th step, the assessment to target component
Using the Acceleration time course a by adjusting above2(t) corresponding acceleration response spectrum, peak accelerator, speed are calculated
Degree and displacement time-histories, evaluation meets situation to preset target component, if being unsatisfactory for target component, with a2(t) conduct
New time-histories determines new baseline correction parameter since third step, repeats the above steps, until meeting all predetermined mesh
Mark.
6th step, acquisition meet acceleration, speed and the displacement time-histories of all targets.
Embodiment 1
The bridge of certain crossover fault, distance of the bridge pier away from tomography are about 2km, and when design needs to consider the close of input earthquake
The influence of tomography pulse characteristic and permanent displacement.Result of study shows that the earthquake corresponding to 2% return period of 50 Annual exceeding probability adds
Speed time-histories peak value is about 0.69g, and bridge considers that Surface Rupture Due influences, and disposably happen suddenly shift value are as follows: the upper and lower disk of tomography
Vertical fault 1.4m, level pull open 0.8m (by 60 ° or so Dip countions), the horizontal changing of the relative positions < 1.4m of dextrorotation.
When horizontal earthquake Acceleration time course is fitted, initially 1992 that acceleration transducer measures are selected in vibration
The accelerogram of Lander earthquake, using the amplitude of peak value 0.69g adjustment initial acceleration time-histories, aimed acceleration response spectrum
Mean value spectrum using tomography away from acceleration response spectrum within the scope of 5km, as shown in Figure 4.According to the method fitting originally researched and proposed
Acceleration, speed and displacement time-histories are as shown in figure 5, be wherein 5% to the worst error of goal response spectrum fitting, from t2Start to
The mean value of latter end permanent displacement about 0.8m, speed time-histories can characterize the characteristic of the big pulse of Near-field ground motion speed, meet
The requirement of related specifications.
Using the approximating method of the near-fault ground motion Acceleration time course that can characterize ground permanent displacement of the application, the party
The seismic acceleration record that method obtains, in addition to can satisfy acceleration response spectrum, peak accelerator and be able to reflect out nearly tomography
The characteristic of earthquake motion velocity pulse, while ground permanent displacement can be characterized, referential is capable of providing when carrying out Aseismic Design
The input data of high seismic acceleration stroke.
Although embodiment disclosed by the application is as above, the content is only to facilitate understanding the application and adopting
Embodiment is not limited to the application.Technical staff in any the application technical field is not departing from this
Under the premise of the disclosed spirit and scope of application, any modification and change can be made in the implementing form and in details,
But the scope of patent protection of the application, still should be subject to the scope of the claims as defined in the appended claims.
Claims (4)
1. the near-fault ground motion Acceleration time course approximating method that one kind can characterize ground permanent displacement, which is characterized in that including
Following steps:
The first step forms primary earthquake Acceleration time course
Primary earthquake Acceleration time course a is generated according to known aimed acceleration response spectrum0(t), it or from real seismic record chooses
Choosing meets preferable accelerogram as initial time-histories a with aimed acceleration response spectrum0(t);The aimed acceleration reaction
Spectrum is that the statistical analysis based on acceleration transducer measurement result obtains;
Second step carries out the adjustment of first time baseline to initially vibration
High-pass filtering is carried out to initially vibration using filter function, is defined by first time baseline Acceleration time course adjusted
For a01(t);
Third step extracts second of baseline correction parameter of seismic acceleration time-histories
To Acceleration time course a01(t) it is integrated, obtains corresponding speed time-histories v0(t) and displacement time-histories d0(t), it then extracts
The parameter of second of baseline correction;
4th step, the seismic acceleration time-histories adjustment based on multiple target
It is frequency is ω that goal response, which is composed discrete,1, ω2..., ωMM single-degree-of-freedom system, to each frequency control point
Successively it is adjusted;
After the first round adjustment for completing M control point, saving time-histories adjusted is a2(t), and acceleration peak value is adjusted as target
Value;To Acceleration time course a2(t) it is integrated to obtain corresponding speed time-histories v2(t), Linear Quasi then is carried out to speed time-histories
It closes;
Acceleration time course is finally subtracted into offset, completes second of baseline correction and the control to permanent displacement;
5th step, the assessment to target component
Using the Acceleration time course a by adjusting above2(t) calculate corresponding acceleration response spectrum, peak accelerator, speed and
It is displaced time-histories, evaluation meets situation to preset target component, if being unsatisfactory for target component, with a2(t) as new
Time-histories determines new baseline correction parameter since third step, repeats the above steps, and the acceleration of all targets is met until obtaining
Degree, speed and displacement time-histories.
2. the near-fault ground motion Acceleration time course approximating method according to claim 1 that ground permanent displacement can be characterized,
It is characterized in that, the extraction process of baseline correction parameter is as follows in third step:
(1) start the time t of offset zero shift1
From displacement time-histories d0(t) terminal tendIt takes forward at equal intervals a little as t1, i.e. t1=tend- n Δ, Δ are time interval;So
Straight line d is used afterwardsn=d0,n-af,nT is fitted t1To tendThe displacement time-histories of this section, utilizes fitting parameter af,nCalculate time shift slope ratio
kn=af,n+1+1/af,n, work as knOccurring the point jumpy corresponding time first is t1;
(2) the corresponding time t of peak acceleratorpga
The corresponding time t of peak acceleratorpgaIt is worth the corresponding time for Acceleration time course amplitude maximum absolute value;
(3) the initial time t of permanent displacement2
Definition does not correct displacement d0(t) time finally intersected with time shaft is td0, take t3=max [tpga,td0], then take t3It arrives
Each moment of record end is t2, i.e. t2=t3+ n Δ t, Δ t are time interval, are fitted t using quadratic function2To tend's
It is displaced time-histories, and subtracts offset from Acceleration time course and initial time-histories is corrected, displacement time-histories is in t after calculating correction2
To tendThe flatness of periodMaximum determines t2Value,
R is t in formula2To tendThe linearly dependent coefficient of correction the displacement time-histories and its fitting a straight line of period, b is t2To tendWhen
Between section correction displacement time-histories least square regression straight slope, σ be correction after be displaced time-histories t2To tendThe side of period
Difference.
3. the near-fault ground motion Acceleration time course fitting side according to claim 1 or 2 that ground permanent displacement can be characterized
Method, which is characterized in that in the 4th step, method that each frequency control point is adjusted are as follows:
If control point frequencies omegakAcceleration time course before adjustment isIt is superimposed an increment time-histories Δ a (t), makes to calculate and react
Compose Sa(ωk) composed with goal responseBetween difference meet required precision, i.e.,
Δ a (t)=Rh (tm-t)B(t,ωk) t≤tm (1)
In formula: tmAt the time of reacting generation for single-degree-of-freedom system maximum;B(t,ωk) be when-frequency envelope curve function;h(tm- t) be
Unit impulse response function indicates are as follows:
Wherein: ζ is damping ratio;φkFor initial phase, guarantee the unit impulse function of superposition in tmPlace takes
Obtain peak value;
In formula (1), R is amplitude regulation coefficient, by control point frequencies omegakThe calculating response spectrum S at placea(ωk) and goal response
SpectrumDifference DELTA S determine, it may be assumed that
Δ S=CR (3)
Wherein:
By calculating above, can be obtained corresponding to frequencies omegakAdjustment time-histories Δ a (t), when by primary earthquake adjusted
JourneyIt can be written as:
The Acceleration time course that will be calculated according to formula (5)As next frequency control point ωk+1Input time-histories, press
It is adjusted according to formula (1)-(5) process.
4. the near-fault ground motion Acceleration time course fitting according to claim 1 or 2 or 3 that ground permanent displacement can be characterized
Method, which is characterized in that in the 4th step, integral a is carried out to Acceleration time course2(t) corresponding speed time-histories v is obtained2(t), then
Linear fit is carried out to speed time-histories, is obtained:
vf(t)=v0-af·t
Then Acceleration time course is obtained in time period t by following formula2-t1Correction offset, i.e.,
Finally by Acceleration time course in t2-t1Section subtracts offset am, complete second of baseline correction and the control to permanent displacement
System.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910630179.1A CN110389379B (en) | 2019-07-12 | 2019-07-12 | Near-fault seismic acceleration time-course fitting method capable of representing ground permanent displacement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910630179.1A CN110389379B (en) | 2019-07-12 | 2019-07-12 | Near-fault seismic acceleration time-course fitting method capable of representing ground permanent displacement |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110389379A true CN110389379A (en) | 2019-10-29 |
CN110389379B CN110389379B (en) | 2020-09-29 |
Family
ID=68286577
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910630179.1A Active CN110389379B (en) | 2019-07-12 | 2019-07-12 | Near-fault seismic acceleration time-course fitting method capable of representing ground permanent displacement |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110389379B (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257434A (en) * | 2020-01-17 | 2020-06-09 | 大连理工大学 | Knocking device for positioning and correcting microseism of surface |
CN111551994A (en) * | 2020-05-11 | 2020-08-18 | 中国地震局地球物理研究所 | Frequency non-stationary seismic oscillation fitting method for matching multiple targets |
CN112269207A (en) * | 2020-11-06 | 2021-01-26 | 中国地震局工程力学研究所 | Strong vibration displacement time course simulation method with permanent displacement |
CN112630827A (en) * | 2020-12-07 | 2021-04-09 | 中国地震局工程力学研究所 | Earth surface acceleration peak value parameter prediction method and system |
CN114578417A (en) * | 2022-03-09 | 2022-06-03 | 北京工业大学 | Pulse type earthquake motion fitting method based on wavelet adjustment |
CN114626409A (en) * | 2022-02-21 | 2022-06-14 | 中铁第四勘察设计院集团有限公司 | Near-fault acceleration pulse identification method, storage medium and computer equipment |
CN114675325A (en) * | 2022-03-13 | 2022-06-28 | 中国地震台网中心 | Mindlin solution-based method for estimating permanent displacement of seismic surface fracture |
CN114966835A (en) * | 2022-04-24 | 2022-08-30 | 福州大学 | Near-fault seismic motion fitting method based on improved time domain superposition method |
CN116626752A (en) * | 2023-06-08 | 2023-08-22 | 大连理工大学 | Novel solving method for ground vibration rotation component based on field surface deformation rate |
CN116882036A (en) * | 2023-09-07 | 2023-10-13 | 清华大学 | Seismic response analysis processing method and device for arch dam |
CN116973972A (en) * | 2023-06-06 | 2023-10-31 | 中国地震局地球物理研究所 | Near-fault earthquake motion fitting method suitable for engineering application |
CN117607967A (en) * | 2024-01-19 | 2024-02-27 | 中国建筑西南设计研究院有限公司 | Genetic algorithm-based earthquake motion baseline correction method and electronic equipment |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090259405A1 (en) * | 2008-04-15 | 2009-10-15 | Battelle Energy Alliance, Llc | Methods, systems, and computer-readable media for generating seismic event time histories |
JP2014181915A (en) * | 2013-03-18 | 2014-09-29 | National Research Institute For Earth Science & Disaster Provention | Earthquake vibration measuring instrument, earthquake vibration measuring system, and earthquake vibration measuring method using the same |
US20140316708A1 (en) * | 2013-04-19 | 2014-10-23 | The Board Of Trustees Of The Leland Stanford Junior University | Oriented Wireless Structural Health and Seismic Monitoring |
CN109143370A (en) * | 2018-07-25 | 2019-01-04 | 中国地震局地球物理研究所 | The bearing calibration of earthquake ground motion acceleration record baseline drift |
CN109186910A (en) * | 2018-10-10 | 2019-01-11 | 东南大学 | A kind of shaking-table test control method for Cross-fault leveling seismic wave |
-
2019
- 2019-07-12 CN CN201910630179.1A patent/CN110389379B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090259405A1 (en) * | 2008-04-15 | 2009-10-15 | Battelle Energy Alliance, Llc | Methods, systems, and computer-readable media for generating seismic event time histories |
JP2014181915A (en) * | 2013-03-18 | 2014-09-29 | National Research Institute For Earth Science & Disaster Provention | Earthquake vibration measuring instrument, earthquake vibration measuring system, and earthquake vibration measuring method using the same |
US20140316708A1 (en) * | 2013-04-19 | 2014-10-23 | The Board Of Trustees Of The Leland Stanford Junior University | Oriented Wireless Structural Health and Seismic Monitoring |
CN109143370A (en) * | 2018-07-25 | 2019-01-04 | 中国地震局地球物理研究所 | The bearing calibration of earthquake ground motion acceleration record baseline drift |
CN109186910A (en) * | 2018-10-10 | 2019-01-11 | 东南大学 | A kind of shaking-table test control method for Cross-fault leveling seismic wave |
Non-Patent Citations (1)
Title |
---|
荣棉水 等: "近断层强震观测记录基线校正的优化方法", 《土木工程学报》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257434A (en) * | 2020-01-17 | 2020-06-09 | 大连理工大学 | Knocking device for positioning and correcting microseism of surface |
CN111551994B (en) * | 2020-05-11 | 2023-01-17 | 中国地震局地球物理研究所 | Frequency non-stationary seismic oscillation fitting method for matching multiple targets |
CN111551994A (en) * | 2020-05-11 | 2020-08-18 | 中国地震局地球物理研究所 | Frequency non-stationary seismic oscillation fitting method for matching multiple targets |
CN112269207A (en) * | 2020-11-06 | 2021-01-26 | 中国地震局工程力学研究所 | Strong vibration displacement time course simulation method with permanent displacement |
CN112269207B (en) * | 2020-11-06 | 2023-04-07 | 中国地震局工程力学研究所 | Strong vibration displacement time course simulation method with permanent displacement |
CN112630827A (en) * | 2020-12-07 | 2021-04-09 | 中国地震局工程力学研究所 | Earth surface acceleration peak value parameter prediction method and system |
CN112630827B (en) * | 2020-12-07 | 2024-02-02 | 中国地震局工程力学研究所 | Surface acceleration peak parameter prediction method and system |
CN114626409A (en) * | 2022-02-21 | 2022-06-14 | 中铁第四勘察设计院集团有限公司 | Near-fault acceleration pulse identification method, storage medium and computer equipment |
CN114626409B (en) * | 2022-02-21 | 2023-09-26 | 中铁第四勘察设计院集团有限公司 | Near fault acceleration pulse identification method, storage medium and computer equipment |
CN114578417B (en) * | 2022-03-09 | 2022-09-20 | 北京工业大学 | Pulse type earthquake motion fitting method based on wavelet adjustment |
CN114578417A (en) * | 2022-03-09 | 2022-06-03 | 北京工业大学 | Pulse type earthquake motion fitting method based on wavelet adjustment |
CN114675325A (en) * | 2022-03-13 | 2022-06-28 | 中国地震台网中心 | Mindlin solution-based method for estimating permanent displacement of seismic surface fracture |
CN114966835A (en) * | 2022-04-24 | 2022-08-30 | 福州大学 | Near-fault seismic motion fitting method based on improved time domain superposition method |
CN114966835B (en) * | 2022-04-24 | 2024-06-04 | 福州大学 | Near-fault earthquake motion fitting method based on improved time domain superposition method |
CN116973972B (en) * | 2023-06-06 | 2024-01-12 | 中国地震局地球物理研究所 | Near-fault earthquake motion fitting method suitable for engineering application |
CN116973972A (en) * | 2023-06-06 | 2023-10-31 | 中国地震局地球物理研究所 | Near-fault earthquake motion fitting method suitable for engineering application |
CN116626752A (en) * | 2023-06-08 | 2023-08-22 | 大连理工大学 | Novel solving method for ground vibration rotation component based on field surface deformation rate |
CN116626752B (en) * | 2023-06-08 | 2023-10-24 | 大连理工大学 | Ground vibration rotation component solving method based on field surface deformation rate |
CN116882036B (en) * | 2023-09-07 | 2023-12-19 | 清华大学 | Seismic response analysis processing method and device for arch dam |
CN116882036A (en) * | 2023-09-07 | 2023-10-13 | 清华大学 | Seismic response analysis processing method and device for arch dam |
CN117607967A (en) * | 2024-01-19 | 2024-02-27 | 中国建筑西南设计研究院有限公司 | Genetic algorithm-based earthquake motion baseline correction method and electronic equipment |
CN117607967B (en) * | 2024-01-19 | 2024-03-26 | 中国建筑西南设计研究院有限公司 | Genetic algorithm-based earthquake motion baseline correction method and electronic equipment |
Also Published As
Publication number | Publication date |
---|---|
CN110389379B (en) | 2020-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110389379A (en) | The near-fault ground motion Acceleration time course approximating method of ground permanent displacement can be characterized | |
CN105891863B (en) | It is a kind of based on highly constrained EKF localization method | |
JP2003344213A (en) | Device and method for evaluating aseismatic performance of building | |
CN111553107B (en) | Random earthquake response analysis and safety evaluation method for pile foundation of liquefiable field | |
Ren et al. | Integrated GNSS/IMU hub motion estimator for offshore wind turbine blade installation | |
CN109446715B (en) | Large-span bridge earthquake response time course analysis method | |
CN107784154B (en) | Earthquake resistance probability evaluation method based on behavior bispectrum | |
CN106643806B (en) | A kind of inertial navigation system alignment precision appraisal procedure | |
CN107101631A (en) | A kind of ship heave measuring method based on auto-adaptive filtering technique | |
Yang et al. | Structure system estimation under seismic excitation with an adaptive extended Kalman filter | |
Hu et al. | Integrated multi‐type sensor placement and response reconstruction method for high‐rise buildings under unknown seismic loading | |
Shan et al. | Interstory drift estimation of nonlinear structure using acceleration measurement with test validation | |
Qi et al. | Influence of foundation type on seismic response of low-rise structures in liquefiable soil | |
Brownjohn et al. | Real-time long-term monitoring and static and dynamic displacements of an office tower, combining RTK GPS and accelerometer data | |
Lei et al. | A Smoothing EKF‐UI‐WDF Method for Simultaneous Identification of Structural Systems and Unknown Seismic Inputs without Direct Feedthrough | |
CN114527506A (en) | Method for simplifying artificial generation of fault seismic oscillation | |
Chikh et al. | Simplified procedure for seismic demands assessment of structures | |
He et al. | Estimation of the displacement time history of high-rise building structures using limited measurement data and structural information | |
Suresh et al. | Influence of diagonal braces in RCC multi-storied frames under wind loads: A case study | |
Millen et al. | Preliminary displacement-based assessment procedure for buildings on liquified soil | |
Bruserud et al. | Current conditions in the northern North Sea | |
Magenes et al. | Some issues on seismic design and assessment of masonry buildings based on linear elastic analysis | |
Atkinson | Use of the uniform hazard spectrum in characterizing expected levels of seismic ground shaking | |
Beck et al. | Comparison between transfer function and modal minimization methods for system identification | |
Castellano et al. | Testing of the seismic isolation system for the bust of Francesco I d'Este in Modena, Italy |
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 |