CN107967240A - Shatter-proof time-histories acceleration optimization algorithm based on artificial ground shaking generating process - Google Patents
Shatter-proof time-histories acceleration optimization algorithm based on artificial ground shaking generating process Download PDFInfo
- Publication number
- CN107967240A CN107967240A CN201711178550.2A CN201711178550A CN107967240A CN 107967240 A CN107967240 A CN 107967240A CN 201711178550 A CN201711178550 A CN 201711178550A CN 107967240 A CN107967240 A CN 107967240A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- spectrum
- mfrac
- time
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
- G06F17/156—Correlation function computation including computation of convolution operations using a domain transform, e.g. Fourier transform, polynomial transform, number theoretic transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
Abstract
The invention belongs to earthquake motion optimization algorithmic technique field, propose a kind of shatter-proof time-histories acceleration optimization algorithm based on artificial ground shaking generating process, and improve the precision initially shaken, it can solve in hydro-structure numerical computations, under seismic loading, nonlinear dynamical damage calculates that the time is long, slow-footed problem.This algorithm is with NB35047 2015《Hydroelectric project hydraulic structure earthquake resistant design code》Standard design response spectrum is composed for target, using the shatter-proof time-histories acceleration function of MATLAB Software Creates.Beneficial effects of the present invention:1. introducing earthquake motion efficiently produces algorithm;2. initially vibrations have preferable fitting precision;3. reduce the time of nonlinear dynamical damage, easy to analyze.
Description
Technical field
The invention belongs to earthquake motion optimization algorithmic technique field, it is proposed that a kind of based on the resistance to of artificial ground shaking generating process
Time-histories acceleration optimization algorithm is shaken, and improves the precision initially shaken, can be solved in hydro-structure numerical computations, it is non-thread
Property Dynamic time history analysis calculate that the time is long, slow-footed problem.
Background technology
In recent years, China has built large quantities of world-class high concrete dams in west area, its height of dam it is more 200 meters with
On, the west area in China is the more hair bands of earthquake, and in the event of earthquake, the possibility of dam body damage will greatly increase, therefore,
It is very necessary, the health monitoring of high concrete dam and peace to analyze degree of impairment of the concrete dam under different earthquake fatigue resistance
Full evaluation is also widely paid close attention to.
However, in traditional non-linear dynamic time-histories, analyzing the response under different earthquake fatigue resistance needs substantial amounts of adjust
Width calculates, so generation increases over time, the ever-increasing shatter-proof time-histories acceleration function of Earthquake Intensity is one and relatively attaches most importance to
The work wanted, gradually increases it is crucial that how to generate one with time increase Earthquake Intensity, and the reaction of different time course
The spectrum function good with target spectrum fitting degree.On the basis of conventional artificial's earthquake motion, there is provided one kind optimization algorithm and plan
Conjunction works well.
The present invention has both sides purpose:One there is provided a kind of earthquake motion efficiently produce algorithm and improve initially shake
The fitting precision of response spectrum;Second, in nonlinear dynamical damage, the sound under different earthquake fatigue resistance can be once calculated
Situation is answered, compared with traditional Incremental Dynamic Analysis, reduces the time calculated after amplitude modulation, while the post processing of response is also carried
Facility is supplied.
The content of the invention
The present invention provides a kind of shatter-proof time-histories acceleration based on artificial ground shaking generating process to optimize algorithm, in tradition
On the basis of synthesizing artificial ground shaking, the precision initially shaken is improved, its response spectrum is closed with the fitting of standard design response spectrum
System is good, recycles MATLAB unconstrained optimizations to obtain target earthquake motion time history to its further fitting.This algorithm can synthesize
Meet the shatter-proof time-histories acceleration function of acceleration response spectrum fitting precision requirement under different time course, can ensure in hydro-structure
It is used for input earthquake in engineering, so as to calculate the response condition under different earthquake fatigue resistance, assessment dam body is in geological process
Under damaged degree and failure probability, shorten the Nonlinear time-history analysis time, traditional Incremental Dynamic Analysis can be replaced, it is right
Hydro-structure numerical simulation has far-reaching influence.
Technical scheme:
Shatter-proof time-histories acceleration optimization algorithm based on artificial ground shaking generating process, step are as follows:
This algorithm is with NB35047-2015《Hydroelectric project hydraulic structure earthquake resistant design code》Standard design response spectrum is
Target is composed, and using the shatter-proof time-histories acceleration function of MATLAB Software Creates, flow chart is as shown in Figure 1, be further illustrated below resistance to
The detailed process of time-histories acceleration function generation is shaken, step is as follows:
(1) generation initially shaken
1) steady earthquake motion time history
By response spectrum and the approximation relation of power spectrum, according to standard design response spectrum, it is calculated in power spectrum and Fu
Leaf amplitude spectrum, formula such as (1), (2), (3), standard design response spectrum are as shown in Figure 2;
In formula:S (ω) is power spectrum;
Sa(ω) is the standard design response spectrum on frequency domain;
ξ is damping ratio;
ω is circular frequency;
P is probability coefficent, generally takes p >=0.85;
T is Strong Ground Motion Duration length;
A (i ω) is Fourier's amplitude spectrum;
Δ ω is frequency interval;
fsFor sample frequency;
Nfft counts for Fourier transformation;
Generate phase spectrum at random in [0,2 π], Fourier's amplitude spectrum is multiplied with phase spectrum, positive and negative circular frequency is combined into
One vector, after inverse Fourier transform, it is the steady earthquake ground motion acceleration time-histories x in time domain to take real part, specific formula
It is as follows:
In formula:X is the steady earthquake ground motion acceleration time-histories in time domain;
G (ω) is phase spectrum;
2) earthquake motion time history x is generated after filtering0
Using filter function, the high-frequency content of the steady earthquake motion time history x (i ω) on frequency domain is filtered out, filter function is such as
Under:
In formula, ω1、ξ1For the first filtering parameter;
ω2、ξ2For the second filtering parameter;
H1(i ω) is the first filter function;
H2(i ω) is the second filter function;
Steady earthquake ground motion acceleration time-histories x in time domain is changed to obtain Acceleration time course x on frequency domain by Fourier
(i ω), is multiplied to obtain filtered amplitude spectrum A with filter function0(i ω), then obtain by inverse Fourier transform filtered
Earthquake ground motion acceleration time-histories x0, specific formula is as follows:
X (i ω)=fft (x) (7)
A0(i ω)=x (i ω) × H1(iω)×H2(iω) (8)
x0=f-1[A0(iω)eig(ω)] (9)
In formula, x (i ω) is the Acceleration time course on frequency domain;
A0(i ω) is filtered amplitude spectrum;
x0For filtered earthquake ground motion acceleration time-histories;
3) iterative process is corrected
Filtered earthquake ground motion acceleration time-histories x0Generation calculates spectrum, and width is corrected with calculating the ratio of spectrum using response spectrum
Value spectrum, formula are as follows:
In formula, A ' (i ω) is revised amplitude spectrum;
Sa1(ω) composes to calculate;
Using revised amplitude spectrum A ' (i ω) as the initial value of next iteration, by formula (9) and (10) constantly into
Row iteration, until to calculate spectrum and Response Spectrum fitting relation good for earthquake motion, otherwise iteration always, it is good to have finally obtained fit correlation
Good initially vibrations xg′.Fig. 3 is earthquake motion xg' response spectrum and standard design response spectrum fit solution.
4) linear function
Since shatter-proof time-histories acceleration function has under different time course, target spectrum is with calculating spectrum Sac(T,tdiff) have well
The property of fitting precision, therefore for different time course tdiff, target spectrum has into the relation of multiple with standard design response spectrum, specifically
Relational expression is as follows:
In formula, Sa(T) it is the standard design response spectrum in time domain;
ttargetFor the object time;
SaT(T,tdiff) composed for target;
tdiffTo be fitted the time of response spectrum;
The target that can be calculated by formula (11) under different time course is composed.Due to the calculating spectrum and mesh under different time course
Mark spectrum needs to meet fitting precision, therefore initially shakes xg' need to multiply linear strain-hardening function f (t), during the purpose is to make different
Response spectrum under journey, into multiple proportion, has good fit correlation even if calculating and composing with standard design response spectrum with target spectrum, public
Formula such as (12), (13):
xnew'=f (t) × xg' (13)
In formula, ttFor the total time-histories elapsed time of earthquake motion;
F (t) is linear function;
xg' it is initially to shake;
xnew' it is the earthquake motion increased with time Earthquake Intensity;
Calculate fitting degree such as Fig. 4 of spectrum and target spectrum, the earthquake motion x increased with time Earthquake Intensitynew' such as Fig. 5 institutes
Show.
(2) unconstrained optimization
Above procedure, which has calculated, initially shakes xnew', spectrum is calculated as seen from Figure 4 composes the plan under long period with target
It is general to close effect, it is therefore desirable to which further optimization calculates.Algorithm proposes Unconstrained Optimization Algorithms, further makes calculating spectrum and mesh
Mark spectrum meets required precision, and unconstrained optimization includes following three parts:
1) standard design response spectrum algorithm
Algorithm with MATLAB softwares and《Hydroelectric project hydraulic structure earthquake resistant design code》In standard design response spectrum
Show that target composes algorithm under different cycles, i.e., according to (11) formula, obtain the target spectrum S under different cyclesaT(T,tdiff)。
2) earthquake motion reaction of formation spectrum algorithm
Step (1) generation earthquake motion xnew' after, linear acceleration method step_by_step integration is carried out to it, obtains the meter under different time course
Calculate spectrum Sac(T,tdiff)。
3) algorithm is optimized
It is iterated optimization to calculating spectrum and target spectrum using least square method and calculates, even if calculates the difference composed with response spectrum
Value reaches minimum.In the calculation, it is necessary to which Reusability standard design response spectrum algorithm, earthquake motion reaction of formation compose algorithm, therefore
Needing largely to calculate the time, initial acceleration Function Fitting precision will produce Optimized Iterative below certain influence,
So x of generationnew' earthquake motion is very crucial.
Optimization is needed differently to shake the acceleration function under elapsed time in the process, therefore computationally intensive, optimization
Formula is as follows:
In formula, Tmax is the cycle maximum of standard design response spectrum;
tmaxEarthquake motion always lasts;
Sac(T,tdiff) composed for the calculating under different cycles.
According to formula (14) optimisation criteria design response spectrum algorithm and earthquake motion reaction of formation spectrum algorithm, make two kinds of differences composed
Value reaches minimum value, and then adjusts earthquake motion time history, draws shatter-proof time-histories acceleration function.Reaction after unconstrained optimization
Spectrum fitting degree is shown in Fig. 7.As seen in Figure 7, fitting precision is good, calculate spectrum target compose around fluctuation less, have compared with
Good fitting character.Fig. 8 is the shatter-proof time-histories acceleration function of generation.
Beneficial effects of the present invention:
1. introducing earthquake motion efficiently produces algorithm
Traditional earthquake motion has ascent stage, steady section and descending branch three parts, and shatter-proof time-histories acceleration function introduces filter
Wave function and earthquake motion efficiently produce algorithm, meet the requirements the response spectrum under any time-histories, shatter-proof time-histories acceleration function
Generation equivalent to the peak accelerator of many earthquake motions, different Earthquake Intensities can be represented, there is deep physical significance.
2. initially vibrations have preferable fitting precision
Algorithm synthesizes 4 5s earthquake ground motion acceleration time-histories, is composed 20s earthquake motions, under identical iterations, calculates
Method compares two methods of 4 5s earthquake motion time histories, the earthquake motion of 20s earthquake motion time histories synthesis, by two kinds of response spectrums not
With the comparison under time-histories, such as Fig. 4, Fig. 6, it can be seen that 4 5s earthquake motion time histories fit within the initial journey being fitted under different time course
Spend, therefore, algorithm is used as using the 20s earthquake motions of 4 5s synthesis initially shakes xnew', improve the essence initially shaken
Degree, has preferable effect.
3. reduce the time of nonlinear dynamical damage, easy to analyze
For arch dam, the free degree is more, and loading characteristic is complicated, and site topography, complex geologic conditions, utilize software numerical value
Simulation, when carrying out Nonlinear time-history analysis, it is possible to reduce a large amount of calculating times under different amplitude modulation, and can be very good to find out
The degree of impairment of dam body under different earthquake fatigue resistance, describes the pass between different earthquake fatigue resistance and degree of injury on a macro scale
System, easily sets up the relation between each response, easy to reprocessing analysis.
Brief description of the drawings
Fig. 1 is shatter-proof time-histories acceleration function generating process flow chart.
Fig. 2 is standard design response spectrum schematic diagram.
Fig. 3 is earthquake motion xg' calculate spectrum and standard design response spectrum fitting schematic diagram.
Fig. 4 is x under different time coursenew' calculate spectrum and target spectrum fitting schematic diagram.
Fig. 5 is the earthquake motion x with time increase Earthquake Intensity increasenew' schematic diagram.
Fig. 6 is 20s earthquake motions xgCalculate spectrum and target spectrum fitting schematic diagram.
Fig. 7 is shatter-proof time-histories acceleration function xETACalculate spectrum and target spectrum fitting schematic diagram.
Fig. 8 is shatter-proof time-histories acceleration function xETASchematic diagram.
In figure:A is the cycle (s);B is acceleration response spectrum (g);C is Strong Ground Motion Duration (s);D is acceleration magnitude
(g);A is water conservancy project specification Plays design response spectrum curve;B is the response spectrum curve of 20s earthquake motions generation;C is 0-5s targets
Response spectrum curve;D is 0-10s goal response spectral curves;E is 0-15s goal response spectral curves;F sets a song to music for 0-20s goal responses
Line;G calculates spectral curve for 0-5s;H calculates spectral curve for 0-10s;I calculates spectral curve for 0-15s;J is calculated for 0-20s and set a song to music
Line;K is earthquake ground motion acceleration time-histories.
Embodiment
Below in conjunction with attached drawing and technical solution, the embodiment of this further explanation.
A kind of earthquake motion based on MATLAB softwares efficiently produces algorithm, including MATLAB software process qualities initially shake and
Unconstrained optimization problem two parts:
1. the generation initially shaken
Composed using standard design response spectrum as target, set parameter value according to specification, obtain the response spectrum such as Fig. 2
Figure.
Set the design parameter that generation initially shakes:Sample frequency fs, Strong Ground Motion Duration length t, earthquake motion
Total time-histories elapsed time tt, it is fitted the time t of response spectrumdiff.Algorithm takes 4 5s earthquake motions synthesis 20s earthquake motions xg', in life
Into earthquake motion xg' when, it is necessary to iterate, after reaching certain iterations, can generate preferably bent with Response Spectrum fitting relation
Line.According to explanation and the specific formula of (1) step is combined, has been synthesized using 0.01s as time interval, the initially vibrations of 20s
xnew', this part simply generates the acceleration function good substantially with target spectrum fitting degree, and Fig. 4 is 20s earthquake motions xnew’
Calculate spectrum and the fit correlation of standard design response spectrum.As seen from Figure 4, for long period part, vibration is more obvious, therefore
Also need to further iteration optimization.
2. unconstrained optimization iterative algorithm
Obtain based on the initially vibrations x initially shaken under generating algorithmnew' afterwards, it is necessary to its further analysis, adopt
The iteration of earthquake motion is carried out with MATLAB unconstrained optimizations, further adjusts shatter-proof time-histories acceleration function.Algorithm has been fitted 4
Response spectrum curve under time-histories, i.e. tdiffTake 5s, 10s, 15s, 20s;ttargetValue is 10s, therefore, according to formula (11), 5s
Target spectrum be 0.5 times of standard design response spectrum, 10s response spectrums are standard design response spectrum, and so on.
The design parameter of unconstrained optimization is set:The cycle maximum of T max of standard design response spectrum, earthquake motion always last
tmax, response spectrum elapsed time T;With xnew' as initially shock value, the calculating spectrum under different time course is calculated, utilizes
Formula (14) makes it minimum with goal response spectral difference value, and MATLAB unconstrained optimizations are iterated it optimization.According to formula
(14), it is necessary to x when can be seen that each suboptimizationnew' each point of earthquake motion optimizes calculating, i.e., need herein
2000 points are fitted, it is met the required precision of 4 targets spectrum, calculate it is very big, it is necessary to the plenty of time, therefore just
Beginning earthquake motion fitting degree is better, more favourable to unconstrained optimization.Shatter-proof time-histories acceleration function after fitting as shown in figure 8,
Calculate spectrum and target spectrum fit solution is as shown in Figure 7.
Claims (1)
1. it is a kind of based on artificial ground shaking generating process shatter-proof time-histories acceleration optimization algorithm, it is characterised in that this algorithm with
NB35047-2015《Hydroelectric project hydraulic structure earthquake resistant design code》Standard design response spectrum is composed for target, using MATLAB
Software Create is shatter-proof time-histories acceleration function, step are as follows:
(1) generation initially shaken
1) steady earthquake motion time history
By response spectrum and the approximation relation of power spectrum, according to standard design response spectrum, power spectrum and Fourier's width is calculated
Value spectrum, formula such as (1), (2), (3):
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>&ap;</mo>
<mfrac>
<mrow>
<mo>-</mo>
<mfrac>
<mi>&xi;</mi>
<mrow>
<mi>&pi;</mi>
<mi>&omega;</mi>
</mrow>
</mfrac>
<msup>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mi>&pi;</mi>
<mrow>
<mi>&omega;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mi>ln</mi>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>4</mn>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>&omega;</mi>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>&omega;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>n</mi>
<mi>f</mi>
<mi>f</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:S (ω) is power spectrum;
Sa(ω) is the standard design response spectrum on frequency domain;
ξ is damping ratio;
ω is circular frequency;
P is probability coefficent, takes p >=0.85;
T is Strong Ground Motion Duration length;
A (i ω) is Fourier's amplitude spectrum;
Δ ω is frequency interval;
fsFor sample frequency;
Nfft counts for Fourier transformation;
Generate phase spectrum at random in [0,2 π], Fourier's amplitude spectrum is multiplied with phase spectrum, positive and negative circular frequency is combined into one
Vector, after inverse Fourier transform, it is the steady earthquake ground motion acceleration time-histories x in time domain to take real part, and specific formula is as follows:
X=f-1[A(iω)eig(ω)] (4)
In formula:X is the steady earthquake ground motion acceleration time-histories in time domain;
G (ω) is phase spectrum;
2) earthquake motion time history x is generated after filtering0
Using filter function, the high-frequency content of the steady earthquake motion time history x (i ω) on frequency domain is filtered out, filter function is as follows:
<mrow>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>i&xi;</mi>
<mn>1</mn>
</msub>
<mfrac>
<mi>&omega;</mi>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
</mfrac>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&omega;</mi>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>i&xi;</mi>
<mn>1</mn>
</msub>
<mfrac>
<mi>&omega;</mi>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
</mfrac>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>H</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>i&xi;</mi>
<mn>2</mn>
</msub>
<mfrac>
<mi>&omega;</mi>
<msub>
<mi>&omega;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&omega;</mi>
<msub>
<mi>&omega;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>i&xi;</mi>
<mn>2</mn>
</msub>
<mfrac>
<mi>&omega;</mi>
<msub>
<mi>&omega;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, ω1、ξ1For the first filtering parameter;
ω2、ξ2For the second filtering parameter;
H1(i ω) is the first filter function;
H2(i ω) is the second filter function;
Steady earthquake ground motion acceleration time-histories x in time domain is changed to obtain Acceleration time course x (i on frequency domain by Fourier
ω), it is multiplied to obtain filtered amplitude spectrum A with filter function0(i ω), then obtain by inverse Fourier transform filtered
Oscillatory acceleration time-histories x0, specific formula is as follows:
X (i ω)=fft (x) (7)
A0(i ω)=x (i ω) × H1(iω)×H2(iω) (8)
x0=f-1[A0(iω)eig(ω)] (9)
In formula, x (i ω) is the Acceleration time course on frequency domain;
A0(i ω) is filtered amplitude spectrum;
x0For filtered earthquake ground motion acceleration time-histories;
3) iterative process is corrected
Filtered earthquake ground motion acceleration time-histories x0Generation calculates spectrum, and amplitude spectrum is corrected with calculating the ratio of spectrum using response spectrum,
Formula is as follows:
<mrow>
<msup>
<mi>A</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mfrac>
<mrow>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, A ' (i ω) is revised amplitude spectrum;
Sa1(ω) composes to calculate;
Using revised amplitude spectrum A ' (i ω) as the initial value of next iteration, constantly change by formula (9) and (10)
Generation, until to calculate spectrum and Response Spectrum fitting relation good for earthquake motion, otherwise iteration always, finally obtain fit correlation well at the beginning of
Beginning earthquake motion xg′;
4) linear function
Since shatter-proof time-histories acceleration function has under different time course, target spectrum is with calculating spectrum Sac(T,tdiff) have it is good fit
The property of precision, therefore for different time course tdiff, target spectrum has into the relation of multiple, physical relationship with standard design response spectrum
Formula is as follows:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>a</mi>
<mi>T</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>t</mi>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
<msub>
<mi>t</mi>
<mrow>
<mi>t</mi>
<mi>arg</mi>
<mi>e</mi>
<mi>t</mi>
</mrow>
</msub>
</mfrac>
<mo>&times;</mo>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, Sa(T) it is the standard design response spectrum in time domain;
ttargetFor the object time;
SaT(T,tdiff) composed for target;
tdiffTo be fitted the time of response spectrum;
The target being calculated by formula (11) under different time course is composed;Since the calculating spectrum under different time course and target spectrum need
Meet fitting precision, therefore initially shake xg' need to multiply linear strain-hardening function f (t), it is anti-under different time course the purpose is to make
It should compose with standard design response spectrum into multiple proportion, even if calculating spectrum and target spectrum has good fit correlation, formula such as (12)
(13):
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>t</mi>
<mi>t</mi>
</msub>
<msub>
<mi>t</mi>
<mrow>
<mi>t</mi>
<mi>arg</mi>
<mi>e</mi>
<mi>t</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
xnew'=f (t) × xg' (13)
In formula, ttFor the total time-histories elapsed time of earthquake motion;
F (t) is linear function;
xg' it is initially to shake;
xnew' it is the earthquake motion increased with time Earthquake Intensity;
(2) unconstrained optimization
This algorithm proposes Unconstrained Optimization Algorithms, calculating spectrum is composed with target and meets required precision, unconstrained optimization bag
Containing following three parts:
1) standard design response spectrum algorithm
This algorithm with MATLAB softwares and《Hydroelectric project hydraulic structure earthquake resistant design code》In standard design response spectrum obtain
Go out target under different cycles and compose algorithm, i.e., according to (11) formula, obtain the target spectrum S under different cyclesaT(T,tdiff), generation
Target is composed;
2) earthquake motion reaction of formation spectrum algorithm
Step (1) generation earthquake motion xnew' after, linear acceleration method step_by_step integration is carried out to it, obtains the calculating spectrum under different time course
Value Sac(T,tdiff), the calculating spectrum of generation;
3) algorithm is optimized
Calculated using least square method calculating to compose to be iterated to optimize with target spectrum, even if calculating spectrum and the difference of response spectrum reach
To minimum;In the calculation, Reusability standard design response spectrum algorithm, earthquake motion reaction of formation spectrum algorithm, it is therefore desirable to a large amount of
The calculating time, initial acceleration Function Fitting precision produces certain influence to Optimized Iterative below, so generation
xnew' earthquake motion is very crucial;
Optimization is needed differently to shake the acceleration function under elapsed time in the process, thus it is computationally intensive, optimize formula
It is as follows:
<mrow>
<mi>min</mi>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<msub>
<mi>T</mi>
<mi>max</mi>
</msub>
</munderover>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<msub>
<mi>t</mi>
<mi>max</mi>
</msub>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>a</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>T</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>T</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>t</mi>
<mi>d</mi>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, Tmax is the cycle maximum of standard design response spectrum;
tmaxEarthquake motion always lasts;
Sac(T,tdiff) composed for the calculating under different cycles;
According to formula (14) optimisation criteria design response spectrum algorithm and earthquake motion reaction of formation spectrum algorithm, reach the difference of two kinds of spectrums
To minimum value, and then earthquake motion time history is adjusted, draw shatter-proof time-histories acceleration function.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711178550.2A CN107967240B (en) | 2017-11-23 | 2017-11-23 | Earthquake-resistant time-course acceleration optimization algorithm based on artificial earthquake motion generation process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711178550.2A CN107967240B (en) | 2017-11-23 | 2017-11-23 | Earthquake-resistant time-course acceleration optimization algorithm based on artificial earthquake motion generation process |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107967240A true CN107967240A (en) | 2018-04-27 |
CN107967240B CN107967240B (en) | 2020-12-11 |
Family
ID=62000595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711178550.2A Active CN107967240B (en) | 2017-11-23 | 2017-11-23 | Earthquake-resistant time-course acceleration optimization algorithm based on artificial earthquake motion generation process |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107967240B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108764073A (en) * | 2018-05-14 | 2018-11-06 | 华南理工大学 | A kind of acceleration filter of combination spectrum energy form fitting is made an uproar and integration method |
CN109002678A (en) * | 2018-10-17 | 2018-12-14 | 中国水利水电科学研究院 | Artificial earthquake analogy method based on hydroelectric project hydraulic structure Aseismic Design |
CN109450405A (en) * | 2018-10-31 | 2019-03-08 | 山东鲁能智能技术有限公司 | A kind of combined type software filtering method and system applied in AD acquisition |
CN109765023A (en) * | 2019-03-12 | 2019-05-17 | 哈尔滨工业大学 | A method of being quickly obtained structural seismic capacity curve |
CN110163134A (en) * | 2019-05-10 | 2019-08-23 | 大连理工大学 | A kind of structural damage area recognizing method based on weighted band-wise least square |
CN111551994A (en) * | 2020-05-11 | 2020-08-18 | 中国地震局地球物理研究所 | Frequency non-stationary seismic oscillation fitting method for matching multiple targets |
CN112241024A (en) * | 2019-07-18 | 2021-01-19 | 中国石油化工股份有限公司 | Method for improving signal-to-noise ratio of seismic data, computer storage medium and system |
CN112507415A (en) * | 2020-11-04 | 2021-03-16 | 江南大学 | Three-dimensional earthquake-resistant design earthquake motion generation method combining orthogonalization and influence matrix method |
CN113217109A (en) * | 2021-06-18 | 2021-08-06 | 中国矿业大学 | Waveform completion method for rockburst mine microseismic monitoring system |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002161578A (en) * | 2000-11-24 | 2002-06-04 | Sekisui House Ltd | Earthquake resistance and wind resistance classification decision system and method, and storage media housing decision program |
CN101236256A (en) * | 2008-01-25 | 2008-08-06 | 北京工业大学 | Multi- damping ratio goal response spectrum compatible artificial earthquake wave synthesis method |
CN101302740A (en) * | 2008-06-24 | 2008-11-12 | 东南大学 | Method for optimally arranging damper for single pylon self-anchoring type rope-suspension bridge shock absorption control |
CN101697014A (en) * | 2009-10-23 | 2010-04-21 | 中国地震灾害防御中心 | Earthquake motion synthesis method and system |
CN103926621A (en) * | 2014-05-06 | 2014-07-16 | 中国地震局工程力学研究所 | Method for building seismic oscillation attenuation relation based on two-stage residual analysis |
CN104849750A (en) * | 2015-04-21 | 2015-08-19 | 大连理工大学 | Nuclear power floor spectrum artificial seismic wave fitting method based on objective wave waveform analysis |
WO2015146533A1 (en) * | 2014-03-24 | 2015-10-01 | 独立行政法人情報通信研究機構 | Acceleration sensor and active anti-vibration apparatus using same |
CN107329932A (en) * | 2017-05-08 | 2017-11-07 | 上海交通大学 | The time-frequency domain Modal Parameters Identification decomposed based on nonlinear frequency modulation component |
-
2017
- 2017-11-23 CN CN201711178550.2A patent/CN107967240B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002161578A (en) * | 2000-11-24 | 2002-06-04 | Sekisui House Ltd | Earthquake resistance and wind resistance classification decision system and method, and storage media housing decision program |
CN101236256A (en) * | 2008-01-25 | 2008-08-06 | 北京工业大学 | Multi- damping ratio goal response spectrum compatible artificial earthquake wave synthesis method |
CN101302740A (en) * | 2008-06-24 | 2008-11-12 | 东南大学 | Method for optimally arranging damper for single pylon self-anchoring type rope-suspension bridge shock absorption control |
CN101697014A (en) * | 2009-10-23 | 2010-04-21 | 中国地震灾害防御中心 | Earthquake motion synthesis method and system |
WO2015146533A1 (en) * | 2014-03-24 | 2015-10-01 | 独立行政法人情報通信研究機構 | Acceleration sensor and active anti-vibration apparatus using same |
CN103926621A (en) * | 2014-05-06 | 2014-07-16 | 中国地震局工程力学研究所 | Method for building seismic oscillation attenuation relation based on two-stage residual analysis |
CN104849750A (en) * | 2015-04-21 | 2015-08-19 | 大连理工大学 | Nuclear power floor spectrum artificial seismic wave fitting method based on objective wave waveform analysis |
CN107329932A (en) * | 2017-05-08 | 2017-11-07 | 上海交通大学 | The time-frequency domain Modal Parameters Identification decomposed based on nonlinear frequency modulation component |
Non-Patent Citations (11)
Title |
---|
MARTÍNEZ CARLOS A. 等: "Optimal design of passive viscous damping systems for buildings under seismic excitation", 《JOURNAL OF CONSTRUCTIONAL STEEL RESEARCH》 * |
SEYEDPOOR S. M. 等: "Optimum shape design of arch dams for earthquake loading using a fuzzy inference system and wavelet neural networks", 《ENGINEERING OPTIMIZATION》 * |
XU QIANG 等: "Probability analysis for the damage of gravity dam", 《ENGINEERING》 * |
何凯: "隔震结构地震响应分析及抗震性能研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
张镭于: "时—频非平稳人造地震动的小波合成及其应用", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
徐强: "混凝土重力坝动力系统的可靠性方法研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
李玲萍 等: "基于MATLAB的数值模拟地震动加速度时程研究", 《陕西电力》 * |
熊健 等: "重力坝坝头孤立坝体在地震作用下的抗滑分析", 《水利与建筑工程学报》 * |
熊欢: "高墩大跨度连续刚构桥地震响应分析研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
白久林 等: "结构抗震分析的耐震时程方法", 《地震工程与工程振动》 * |
白久林: "钢筋混凝土框架结构地震主要失效模式分析与优化", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108764073B (en) * | 2018-05-14 | 2021-09-21 | 华南理工大学 | Acceleration noise filtering and integrating method combined with spectrum energy form fitting |
CN108764073A (en) * | 2018-05-14 | 2018-11-06 | 华南理工大学 | A kind of acceleration filter of combination spectrum energy form fitting is made an uproar and integration method |
CN109002678A (en) * | 2018-10-17 | 2018-12-14 | 中国水利水电科学研究院 | Artificial earthquake analogy method based on hydroelectric project hydraulic structure Aseismic Design |
CN109002678B (en) * | 2018-10-17 | 2020-08-18 | 中国水利水电科学研究院 | Artificial earthquake simulation method based on hydroelectric engineering hydraulic structure earthquake-resistant design |
CN109450405B (en) * | 2018-10-31 | 2022-09-16 | 山东鲁软数字科技有限公司智慧能源分公司 | Combined type software filtering method and system applied to AD acquisition |
CN109450405A (en) * | 2018-10-31 | 2019-03-08 | 山东鲁能智能技术有限公司 | A kind of combined type software filtering method and system applied in AD acquisition |
CN109765023A (en) * | 2019-03-12 | 2019-05-17 | 哈尔滨工业大学 | A method of being quickly obtained structural seismic capacity curve |
CN110163134A (en) * | 2019-05-10 | 2019-08-23 | 大连理工大学 | A kind of structural damage area recognizing method based on weighted band-wise least square |
CN110163134B (en) * | 2019-05-10 | 2022-03-29 | 大连理工大学 | Structural damage region identification method based on sub-band weighted least square |
CN112241024B (en) * | 2019-07-18 | 2024-04-09 | 中国石油化工股份有限公司 | Method for improving signal-to-noise ratio of seismic data, computer storage medium and system |
CN112241024A (en) * | 2019-07-18 | 2021-01-19 | 中国石油化工股份有限公司 | Method for improving signal-to-noise ratio of seismic data, computer storage medium and system |
CN111551994A (en) * | 2020-05-11 | 2020-08-18 | 中国地震局地球物理研究所 | Frequency non-stationary seismic oscillation fitting method for matching multiple targets |
CN111551994B (en) * | 2020-05-11 | 2023-01-17 | 中国地震局地球物理研究所 | Frequency non-stationary seismic oscillation fitting method for matching multiple targets |
CN112507415A (en) * | 2020-11-04 | 2021-03-16 | 江南大学 | Three-dimensional earthquake-resistant design earthquake motion generation method combining orthogonalization and influence matrix method |
CN112507415B (en) * | 2020-11-04 | 2024-04-19 | 江南大学 | Three-way earthquake-resistant design earthquake motion generation method combining orthogonalization and influence matrix method |
CN113217109B (en) * | 2021-06-18 | 2022-06-10 | 中国矿业大学 | Waveform completion method for rockburst mine microseismic monitoring system |
CN113217109A (en) * | 2021-06-18 | 2021-08-06 | 中国矿业大学 | Waveform completion method for rockburst mine microseismic monitoring system |
Also Published As
Publication number | Publication date |
---|---|
CN107967240B (en) | 2020-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107967240A (en) | Shatter-proof time-histories acceleration optimization algorithm based on artificial ground shaking generating process | |
Mukherjee et al. | Wavelet-based generation of spectrum-compatible time-histories | |
Park et al. | Evaluation of global atmospheric solvers using extensions of the Jablonowski and Williamson baroclinic wave test case | |
CN114578417B (en) | Pulse type earthquake motion fitting method based on wavelet adjustment | |
CN103308804A (en) | Method for extracting time-frequency parameters of power quality disturbance signals on basis of fast K-S (Kaiser-S) transformation | |
Dzwinel et al. | Supermodeling in simulation of melanoma progression | |
CN106814391A (en) | Ground micro-seismic state event location method based on Fresnel zone tomographic inversion | |
CN109002678B (en) | Artificial earthquake simulation method based on hydroelectric engineering hydraulic structure earthquake-resistant design | |
Li et al. | Generate tri-directional spectra-compatible time histories using HHT method | |
CN108169794A (en) | Nearly tomography Artificial Seismic Wave and the matched method of adjustment of Code response spectrum | |
Iyengar et al. | Generation of spectrum compatible accelerograms | |
Li et al. | Simulation and generation of spectrum-compatible ground motions based on wavelet packet method | |
CN113484911B (en) | Earthquake motion evolution power spectrum-based earthquake displacement response calculation method | |
CN109446552B (en) | Time domain calculation method for structural fatigue life under multi-axis related random excitation | |
CN110428874A (en) | A method of natural gas accumulation history is formed using Pyrolysis Experiment carbon isotope | |
Chen et al. | Consistency improved random flow generation method for large eddy simulation of atmospheric boundary layer | |
CN116227242B (en) | Random irregular wave simulation method based on white noise filtering | |
Liu et al. | Numerical verification of a two-layer Boussinesq-type model for surface gravity wave evolution | |
Mukhopadhyay et al. | Wavelet-based generation of accelerogram-consistent, spectrum-compatible motions: New algorithms and short-period overestimation | |
Shrikhande et al. | On the characterisation of the phase spectrum for strong motion synthesis | |
Zi-lan et al. | Stochastic seismic response analysis of engineering site considering correlations of critical soil dynamic parameters | |
Temperton | Some experiments in dynamic initialization for a simple primitive equation model | |
CN112285770B (en) | Seismic wave phase spectrum disturbance method based on real number wavelet transformation | |
Fang et al. | Boussinesq-type equations for nonlinear evolution of wave trains | |
Alexey et al. | Mathematical modeling of rhytmocardiografic signals spectrum |
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 |