CN104573249A - Time-variant ARMA model based non-stable wind speed simulation method - Google Patents
Time-variant ARMA model based non-stable wind speed simulation method Download PDFInfo
- Publication number
- CN104573249A CN104573249A CN201510023180.XA CN201510023180A CN104573249A CN 104573249 A CN104573249 A CN 104573249A CN 201510023180 A CN201510023180 A CN 201510023180A CN 104573249 A CN104573249 A CN 104573249A
- Authority
- CN
- China
- Prior art keywords
- mrow
- wind speed
- msub
- time
- varying
- 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.)
- Pending
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 42
- 241001123248 Arma Species 0.000 title claims abstract description 35
- 238000001228 spectrum Methods 0.000 claims abstract description 43
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 11
- 238000005311 autocorrelation function Methods 0.000 claims description 8
- 238000005314 correlation function Methods 0.000 abstract description 12
- 238000010586 diagram Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 4
- 241001003127 Tarma Species 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000012731 temporal analysis Methods 0.000 description 3
- 238000000700 time series analysis Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000005654 stationary process Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention discloses a time-variant ARMA model based non-stable wind speed simulation method. The time-variant ARMA model based non-stable wind speed simulation method comprises the following steps of 1 determining parameters of a target non-stable wind speed simulated wind speed model and a time-variant ARMA model; 2 dispersing approximate non-stable fluctuating wind speed into a plurality of stable fluctuating wind speeds within a short time interval, and meanwhile adopting a non-uniform modulation function to modulate a power spectrum of each section of stable fluctuating wind speeds into an evolution spectrum; 3 establishing the time-variant ARMA model of non-stable wind speed and performing non-stable wind speed simulation; 4 outputting and displaying a simulation result and comparing the power spectrum of simulation point wind speeds and respective and a self-correlation function and a cross-correlation function with a target power spectrum, a target self-correlation function and a target cross-correlation function. The time-variant ARMA model based non-stable wind speed simulation method is used for fluctuating wind speed non-stability simulation and ensuring that the frequency component of the fluctuating wind speed changes with time, namely a time-varying characteristic of the power spectrum.
Description
Technical Field
The invention relates to a method for simulating non-stationary wind speed by adopting time series analysis, in particular to a non-stationary wind speed simulation method based on a time-varying ARMA model.
Background
For buildings (structures) such as large-span space structures, large-span bridges, high-rise (super) building structures, high-rise structures (such as guyed masts, television towers, chimneys and the like), wind load is one of the control loads of the wind-resistant design of the structure. The wind resistance analysis of the structure firstly needs to obtain sample data of wind load, and the current main means for determining the wind load comprises wind tunnel test, field actual measurement, numerical simulation and the like. With the rapid development of computer technology and the deep research of people on random process numerical simulation technology, the wind speed time curve obtained by adopting a numerical simulation method can consider the arbitrariness of conditions such as fields, wind spectrum characteristics, building characteristics and the like, so that the load obtained by simulation is as close to the actual wind power of a structure as possible, can meet the arbitrariness of certain statistical characteristics, is more representative than actual records, and is widely applied to actual engineering.
The non-stationary characteristic is a phenomenon (such as turbulent flow of an atmospheric boundary layer, thunderstorm strong wind, earthquake and the like) commonly existing in various random loads in the nature, and the amplitude and the frequency of the phenomenon change along with time, so that the non-stationary characteristic of the wind is a factor which needs to be considered when the numerical simulation is carried out on pulsating wind under certain specific environments. Particularly in downburst (i.e., a very sudden and destructive strong wind that is violently hit the ground in a violent downdraft in thunderstorm weather and propagates all around the ground from the point of impact), its very strong non-stationarity is likely to generate a greater dynamic response to the structure. Numerical simulation of the current stochastic process is mainly classified into two categories: the method is a spectral representation method based on trigonometric series superposition, and is a regression method based on a linear filtering technology. These conventional wind speed simulation methods generally approximate natural wind as a stationary gaussian process that experiences each state and perform the simulation based on the existing fluctuating wind speed power spectrum. For benign weather wind with good stability on an open and flat field, the stability assumption can be basically satisfied. However, analysis of a large number of actual test data shows that many wind speed records in a complex terrain in a strong wind environment do not meet the stability requirement. Particularly, when steady wind speed is assumed, non-steady data needs to be discarded in the case of non-steady pulsating wind in a complex terrain and strong wind environment, which can cause large analysis errors, for example, turbulence intensity values can be overestimated, and further accuracy of subsequent analysis is influenced.
Time series analysis (Time series analysis) is a statistical method of dynamic data processing. The invention is based on the random process theory and the mathematical statistics method, and researches the statistical rule followed by the random data sequence to solve the practical problem. An Autoregressive moving average model (ARMA model for short) is an important method for researching time series, and is formed by mixing an Autoregressive model (AR model for short) and a moving average model (MA model for short) as a basis. The method is based on a non-stationary wind speed model, adopts a Time series autoregressive moving average model and considers the Time-Varying property of ARMA model coefficients, namely a TARMA (Time-Varying ARMA) method, and carries out numerical simulation on the random process of non-stationary pulsating wind.
Disclosure of Invention
The invention aims to provide a non-stationary wind speed simulation method based on a time-varying ARMA model, which simulates the non-stationary of a pulsating wind speed and ensures that the frequency component of the non-stationary wind speed simulation method is time-varying, namely the time-varying property of a power spectrum; and simultaneously comparing the power spectrum and the auto-correlation function of the wind speed at the simulation point with the target power spectrum and the target auto-correlation function to ensure the effectiveness of the simulation.
In order to achieve the purpose, the invention has the following conception: the method comprises the steps of dispersing the unstable pulsating wind speed into a plurality of short-time sequences which can be approximated to the stable pulsating wind speed within a short enough time interval delta t, modulating the power spectrum through a non-uniform modulation function according to the 'evolutionary spectrum' theory to obtain a time-varying power spectrum, namely an evolutionary spectrum, generating the unstable pulsating wind speed through an established TARMA time-varying model of the unstable pulsating wind speed and based on the modulated evolutionary spectrum, and generating the wind speed while considering the spatial correlation of the pulsating wind.
According to the inventive concept, the invention adopts the following technical scheme: a non-stationary wind speed simulation method based on a time-varying ARMA model is characterized by comprising the following steps:
firstly, determining parameters of a wind speed model for simulating a target non-steady wind speed and a time varying ARMA model;
secondly, approximately dispersing the non-steady pulsating wind speed into a plurality of steady pulsating wind speeds in a short time interval, and simultaneously modulating the power spectrum of each section of steady pulsating wind speed into an evolutionary spectrum by adopting a non-uniform modulation function;
thirdly, establishing a time-varying ARMA model of the unsteady fluctuating wind speed, and carrying out unsteady wind speed simulation;
and fourthly, outputting and displaying a simulation result, and comparing the power spectrum and the auto-correlation function of the wind speed at the simulation point with the target power spectrum and the target auto-correlation function.
Preferably, the first step represents the downwind non-stationary wind speed as two parts, a deterministic time-varying mean wind and a zero mean stationary pulsating wind.
Preferably, the non-uniform modulation function in the second step is expressed by the following formula:
in the formula: omega is the circular frequency; z is a radical ofjIs the height of a certain point in space vertical to the ground;the time-varying average wind speed at a certain point in space;and counting the average wind speed for the non-steady wind speed at a certain point in space.
Preferably, the wind speed model of the non-stationary wind speed simulation is represented by the following formula:
in the formula,is the time-varying mean wind, u, at space (x, y, z)*(t) is the zero mean smooth pulsating wind process.
Preferably, the time-varying ARMA model employs the following equation:
wherein U (t) is a zero-mean non-stationary random process vector, Ai(t) is a time-varying autoregressive coefficient matrix, Bj(t) is a time-varying sliding regression coefficient matrix, p is an autoregressive order, q is a sliding regression order, and X (t) is a normally distributed white noise sequence with a variance of 1.
The non-stationary wind speed simulation method based on the time-varying ARMA model has the following advantages: the invention simulates the power spectrum of the unsteady fluctuating wind speed and the target average power spectrum of each pointComparing the wind speeds and simulating the auto-cross correlation function and the target auto-cross correlation function of the wind speeds at each pointAnd comparing and matching the results. PulseThe amplitude of the dynamic wind speed is related to the time-varying average wind speed, and the larger the time-varying average wind speed is, the larger the amplitude of the pulsating wind speed is, which is consistent with the actual wind field characteristics.
Drawings
FIG. 1 is a schematic view of non-stationary wind speed simulation points along the vertical direction of the ground;
FIGS. 2(a) and 2(b) are schematic diagrams of downburst models illustrating different stages of downburst movement toward the building;
FIG. 3 is a design framework diagram of a non-stationary wind speed simulation method based on a time-varying ARMA model;
FIG. 4 is a flowchart of a non-stationary wind speed simulation method routine based on a time-varying ARMA model;
FIGS. 5(a) to 5(c) are schematic diagrams of non-stationary and pulsating wind speeds at 800s at respective simulated spatial points along the ground vertical heights 10m, 25m and 80 m; FIGS. 5(d) to 5(f) are schematic views of a non-stationary pulsating wind speed at 800s at simulated spatial points along the ground vertical height 10m, 25m and 80 m;
6(a) through 6(c) are schematic diagrams of non-uniform modulation functions for various simulated spatial points; 6(d) to 6(f) are schematic diagrams of time-varying power spectra for various simulated spatial points; FIGS. 6(g) to 6(i) are schematic diagrams comparing simulated power spectra at various simulated spatial points with target spectra;
fig. 7(a) to 7(f) are schematic diagrams comparing the auto and cross correlation function target with the target auto and cross correlation function at 10m, 25m, 80m, 10m and 25m, 25m and 80m, 10m and 80m, and 80m of each simulated spatial point along the vertical height of the ground.
Detailed Description
The simulation procedure for non-stationary pulsating wind speeds (hereinafter storm flow as an example) according to the present invention is described in further detail below with reference to the accompanying drawings.
The invention relates to a non-stationary wind speed simulation method based on a time-varying ARMA model, which comprises the following steps:
firstly, determining a wind speed model (comprising a wind speed model, a simulation space point, a terrain parameter, simulation duration and the like) of a target non-stationary wind speed simulation, a target power spectrum and each parameter (an autoregressive order and a sliding regression order) of a time-varying ARMA model. The downwind non-stationary wind speed is expressed as two parts of deterministic time-varying mean wind and zero mean stationary pulsating wind.
In the first step, the wind speed model for non-stationary wind speed simulation is expressed as the following formula (1):
in the formula,is the time-varying mean wind, u, at space (x, y, z)*(t) is the zero mean smooth pulsating wind process. From the above stationary and non-stationary wind speed models, it can be seen that if the total downwind direction record U (t) itself is strictly stationaryAnd in a random process, the time-varying mean wind U (t) is degraded into a constant mean wind in a traditional steady wind model. In fact, the definition in the non-stationary wind speed model is the simplest model to describe the non-stationary process: stationary process + trend term.
As a non-stationary wind speed simulation was performed on the three points shown in fig. 1 using the downburst model shown in fig. 2(a) and 2(b), the autoregressive order p of the time-varying ARMA model was 16, and the slip regression order q was 1. The simulation point is located along the downburst moving direction and is 3500m away from the downburst thunderstorm center. The downdraft storm flow wind speed model adopts an Oseguera and Bowles average wind speed model and a Vicroy vertical distribution model, and the maximum wind speed V in the vertical distribution wind speedmax80m/s at a height position Zmax67 m; radial maximum wind speed V at a certain height in wind speed fieldr,max47m/s, horizontal distance r from downburst centermax1000m, radial length scaling factor Rr700 m; the time-dependent change in the intensity of a thunderstorm is represented by the following formula (2):
downburst translation velocity Vo8 m/s. The upper limit cut-off frequency is 4 pi rad, N is 2^11,meanwhile, the movement of the downburst is considered, the simulation time interval delta t is 0.1s, and the simulation time duration is 800 s.
Secondly, approximately dispersing the non-stable pulsating wind speed into a plurality of stable pulsating wind speeds within a short time interval delta t, and simultaneously modulating the power spectrum of each section of stable pulsating wind speed into an evolutionary spectrum by adopting a non-uniform modulation function; the existing power spectrum (such as a Kaimal spectrum) of the pulsating wind speed is selected, the Kaimal spectrum is taken as an example below, and the power spectrum is modulated into a Kaimal evolutionary spectrum by adopting a Kaimal non-uniform modulation function.
The power spectrum using Kaimal is represented by the following formula (3):
the Kaimal non-uniform modulation function is expressed as the following formula (4):
in the formula: omega is the circular frequency; z is a radical ofjIs the height of a certain point in space vertical to the ground;the time-varying average wind speed at a certain point in space;and counting the average wind speed for the non-steady wind speed at a certain point in space.
And thirdly, establishing a time-varying ARMA model of the unsteady fluctuating wind speed, and carrying out unsteady wind speed simulation according to the frame diagram shown in FIG. 3 and the flow diagram shown in FIG. 4, namely, considering the time-varying property of an ARMA coefficient matrix and the spatial correlation of fluctuating wind, wherein the Matlab language can be adopted in the third step.
The third step comprises the following specific steps:
assume that the following equation (5) is used for the simulated M-point non-stationary fluctuating wind speed:
U(t)=[u1(t) u2(t) … uM(t)]T (5)
the formula of the time-varying ARMA model is shown in the following formula (6):
wherein U (t) is a zero-mean non-stationary random process vector, Ai(t) is a time-varying autoregressive coefficient matrix (MxM order), Bj(t) is the moment of the time-varying slip regression coefficientThe matrix (MxM order), p is an autoregressive order, q is a sliding regression order, and X (t) is a white noise sequence with a variance of 1 and normal distribution, and satisfies the following formula (7):
as can be seen from equation (7), the key to calculating the pulsating wind speed U (t) is to determine Ai(t)、Bj(t) a coefficient matrix.
When t is equal to t0When the formula (7) is represented by the following formula (8):
two sides simultaneously right ride UT(t0-k Δ t) (k ═ 1, 2.. p), and mathematically expected to yield the following formula (9):
according to the definition of the correlation function, i.e. the following formula (10):
pair (8) then two sides multiply X right simultaneouslyT(t0-l Δ t) (l ═ 1, 2.. q), and taking mathematical expectations, we obtain the following formula (11):
namely, the following formula (12):
in addition, the current output U (t) can be seen from equation (9)0) Dependent only on the current input X (t)0) And past input X (t)0τ) with future input X (t)0+ τ) is irrelevant (τ > 0), i.e., U (t)0) And X (t)0+ τ) are independent of each other, then: rux(. tau.) -. 0,. tau. > 0 and Rux(-τ)=Rxu(τ)。
When j is 0, k-j > 0, then Rxu[(k-j)Δt]=0。
Combining formula (11) with formula (12) and developing the following formula (13):
wherein,
according to the formula of Vena-Sinkian, the following formula (14) is obtained:
can find [ Ruu(jΔt)]In this case, it is also necessary to know [ R ]ux(lΔt)]The equation can be found outi(t)、Bj(t) a coefficient matrix.
When q in formula (10) is 0, the following formula (15) is obtained:
two sides simultaneously right multiply XT(t0-l Δ t) (l ═ 0,1, 2.. q) and mathematically expected, yielding the following formula (16):
pair (15) right-handed at two sidesT(t0-k Δ t) (k ═ 0,1,2,. p), mathematically expected to yield the following formula (17):
expansion type <math>
<mrow>
<msub>
<mi>R</mi>
<mi>uu</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>kΔt</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>q</mi>
</munderover>
<msub>
<mover>
<mi>A</mi>
<mo>‾</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mi>uu</mi>
</msub>
<mo>[</mo>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mi>Δt</mi>
<mo>]</mo>
<mo>,</mo>
</mrow>
</math> Obtaining the following formula (18):
from the foregoing [ R ]uu(jΔt)]And equationCan find out
And a second formula from formulae (16), (17), andcan obtain the product <math>
<mrow>
<msub>
<mi>R</mi>
<mi>uu</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</munderover>
<msub>
<mover>
<mi>A</mi>
<mo>‾</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mi>uu</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>iΔt</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>B</mi>
<mn>0</mn>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</math> Namely, the following formula (19):
at this point, the matrix at the right end of the above equation is completely solved, and Cholesky decomposition is performed on the right equation to obtain B0(t0) Then, the matrix [ R ] is obtained by substituting formula (16) in a recursion mannerux(lΔt)]. The matrix [ R ] obtaineduu(jΔt)]And [ R ]ux(lΔt)]By substituting formula (13), a TARMA time-varying model coefficient matrix A can be determinedi(t)、Bj(t), finally obtaining t0The wind speed at that moment.
And fourthly, outputting and displaying a simulation result, and comparing the power spectrum and the auto-correlation function of the wind speed at the simulation point with the target power spectrum and the target auto-correlation function as shown in figure 5 to ensure the effectiveness of the simulated wind speed.
The above steps are analyzed and verified based on a calculation program of a non-stationary wind speed simulation method which is compiled by a Matlab platform and is based on a time-varying ARMA model.
As can be seen from fig. 5(a) to 5(f), the amplitude of the pulsating wind speed is related to the magnitude of the time-varying average wind speed, and the amplitude of the pulsating wind speed is larger as the time-varying average wind speed is larger, which corresponds to the actual wind field characteristics. To further illustrate the effectiveness of the simulation, the power spectrum and the target average power spectrum are simulated for the non-stationary pulsating wind speed at each point(FIGS. 6(a) to 6 (i))Shown) were compared and the downburst Kaimal non-uniform modulation function at three different points is shown in fig. 6. And auto and cross correlation function and target auto and cross correlation functionComparison was made (shown in fig. 7(a) to 7 (f)), and it can be seen that the results are consistent. It can be seen from the cross-correlation function graph of fig. 7 that the cross-correlation of the pulsating wind speed at each point decreases with increasing distance.
Various modifications and changes may be made to the present invention by those skilled in the art. Thus, it is intended that the present invention cover the modifications and variations of this invention provided they come within the scope of the appended claims and their equivalents.
Claims (5)
1. A non-stationary wind speed simulation method based on a time-varying ARMA model is characterized by comprising the following steps:
firstly, determining parameters of a wind speed model for simulating a target non-stationary wind speed, a target power spectrum and a time-varying ARMA model, and simultaneously determining the time-varying average wind speed of each simulation space point;
secondly, modulating the target power spectrum into an evolutionary spectrum by adopting a non-uniform modulation function, and dispersing the unstable fluctuating wind speed within a plurality of short time intervals;
thirdly, establishing a time-varying ARMA model of the unsteady fluctuating wind speed, and simulating the unsteady fluctuating wind speed according to the previously determined time-varying ARMA model parameters and an evolutionary spectrum;
and fourthly, determining the final non-stationary wind speed according to the time-varying average wind speed and the generated non-stationary pulsating wind speed, outputting and displaying a simulation result, and simultaneously comparing the power spectrum and the auto-correlation function of the wind speed at the simulation point with the target power spectrum and the target auto-correlation function.
2. The time-varying ARMA model-based non-stationary wind speed simulation method according to claim 1, wherein the first step represents downwind non-stationary wind speed as two parts of deterministic time-varying mean wind and zero mean stationary pulsating wind.
3. The time-varying ARMA model-based non-stationary wind speed simulation method according to claim 1, wherein the non-uniform modulation function in the second step is represented by the following formula:
in the formula: omega is the circular frequency;is the height of a certain point in space vertical to the ground;the time-varying average wind speed at a certain point in space;and counting the average wind speed for the non-steady wind speed at a certain point in space.
4. The time-varying ARMA model-based non-stationary wind speed simulation method according to claim 1, wherein the non-stationary wind speed simulated wind speed model is represented by the following formula:
in the formula,is a spaceThe average wind of the wind is changed in time,is a zero mean smooth pulsating wind process.
5. The time-varying ARMA model-based non-stationary wind speed simulation method according to claim 1, characterized in that the time-varying ARMA model employs the following formula:
wherein,is a zero-mean non-stationary random process vector,is a time-varying matrix of auto-regressive coefficients,is a time-varying sliding regression coefficient matrix, p is the autoregressive order, q is the sliding regression order,is a variance of1. A normally distributed white noise sequence.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510023180.XA CN104573249A (en) | 2015-01-16 | 2015-01-16 | Time-variant ARMA model based non-stable wind speed simulation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510023180.XA CN104573249A (en) | 2015-01-16 | 2015-01-16 | Time-variant ARMA model based non-stable wind speed simulation method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104573249A true CN104573249A (en) | 2015-04-29 |
Family
ID=53089303
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510023180.XA Pending CN104573249A (en) | 2015-01-16 | 2015-01-16 | Time-variant ARMA model based non-stable wind speed simulation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104573249A (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104951657A (en) * | 2015-06-24 | 2015-09-30 | 上海大学 | Non-stationary air speed prediction method based on Bubble wavelet kernel |
CN105046044A (en) * | 2015-05-29 | 2015-11-11 | 上海大学 | Non-stationary wind speed prediction method based on optimum wavelet-packet transformation |
CN105138782A (en) * | 2015-09-02 | 2015-12-09 | 上海大学 | Non-stationary pulse wind speed high-precision prediction method based on EEMD-ELM |
CN105205495A (en) * | 2015-09-02 | 2015-12-30 | 上海大学 | Non-stationary fluctuating wind speed forecasting method based on EMD-ELM |
CN105469694A (en) * | 2016-01-20 | 2016-04-06 | 肖锐 | Intelligent and safe suspension type sleeve chimney |
CN107480325A (en) * | 2017-07-03 | 2017-12-15 | 河海大学 | The non-stationary non-gaussian earthquake motion time history analogy method of spatial variability |
CN110826197A (en) * | 2019-10-21 | 2020-02-21 | 西南交通大学 | Wind speed field simulation method based on improved Cholesky decomposition closed solution |
CN111368392A (en) * | 2019-12-31 | 2020-07-03 | 重庆大学 | Single-sample non-stationary wind speed simulation method based on MEMD and SRM |
CN111737794A (en) * | 2020-05-26 | 2020-10-02 | 国网天津市电力公司电力科学研究院 | Method for constructing wind characteristic model for overhead transmission line high-rise tower |
CN113468692A (en) * | 2021-07-19 | 2021-10-01 | 大连理工大学 | Three-dimensional wind field efficient simulation method based on delay effect |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20110118246A (en) * | 2010-04-23 | 2011-10-31 | 전남대학교산학협력단 | System for calculating parameter and charge rate of battery |
CN104077478A (en) * | 2014-06-26 | 2014-10-01 | 华东交通大学 | Numerical simulation method for downburst non-stationary fluctuating wind speed |
-
2015
- 2015-01-16 CN CN201510023180.XA patent/CN104573249A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20110118246A (en) * | 2010-04-23 | 2011-10-31 | 전남대학교산학협력단 | System for calculating parameter and charge rate of battery |
CN104077478A (en) * | 2014-06-26 | 2014-10-01 | 华东交通大学 | Numerical simulation method for downburst non-stationary fluctuating wind speed |
Non-Patent Citations (3)
Title |
---|
LIZHONG CHEN ET AL.: "Vector Time-Varying Autoregressive (TVAR) Modeling Nonstationary Downburst Wind Speeds", 《THE FOURTH INTERNATIONAL SYMPOSIUM ON COMPUTATIONAL WIND ENGINEORING》 * |
李春祥等: "基于ARMA模型模拟高架桥的脉动风速时程", 《振动与冲击》 * |
李锦华等: "非平稳脉动风速的数值模拟", 《振动与冲击》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105046044A (en) * | 2015-05-29 | 2015-11-11 | 上海大学 | Non-stationary wind speed prediction method based on optimum wavelet-packet transformation |
CN104951657A (en) * | 2015-06-24 | 2015-09-30 | 上海大学 | Non-stationary air speed prediction method based on Bubble wavelet kernel |
CN105138782A (en) * | 2015-09-02 | 2015-12-09 | 上海大学 | Non-stationary pulse wind speed high-precision prediction method based on EEMD-ELM |
CN105205495A (en) * | 2015-09-02 | 2015-12-30 | 上海大学 | Non-stationary fluctuating wind speed forecasting method based on EMD-ELM |
CN105469694A (en) * | 2016-01-20 | 2016-04-06 | 肖锐 | Intelligent and safe suspension type sleeve chimney |
CN107480325B (en) * | 2017-07-03 | 2019-08-20 | 河海大学 | The non-stationary non-gaussian earthquake motion time history analogy method of spatial variability |
CN107480325A (en) * | 2017-07-03 | 2017-12-15 | 河海大学 | The non-stationary non-gaussian earthquake motion time history analogy method of spatial variability |
CN110826197A (en) * | 2019-10-21 | 2020-02-21 | 西南交通大学 | Wind speed field simulation method based on improved Cholesky decomposition closed solution |
CN111368392A (en) * | 2019-12-31 | 2020-07-03 | 重庆大学 | Single-sample non-stationary wind speed simulation method based on MEMD and SRM |
CN111368392B (en) * | 2019-12-31 | 2024-04-05 | 重庆大学 | Single-sample non-stationary wind speed simulation method based on MEMD and SRM |
CN111737794A (en) * | 2020-05-26 | 2020-10-02 | 国网天津市电力公司电力科学研究院 | Method for constructing wind characteristic model for overhead transmission line high-rise tower |
CN113468692A (en) * | 2021-07-19 | 2021-10-01 | 大连理工大学 | Three-dimensional wind field efficient simulation method based on delay effect |
CN113468692B (en) * | 2021-07-19 | 2022-05-13 | 大连理工大学 | Three-dimensional wind field efficient simulation method based on delay effect |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104573249A (en) | Time-variant ARMA model based non-stable wind speed simulation method | |
Huang et al. | Turbulence and vertical fluxes in the stable atmospheric boundary layer. Part I: A large-eddy simulation study | |
Aboshosha et al. | Turbulence characterization of downbursts using LES | |
Palma et al. | Linear and nonlinear models in wind resource assessment and wind turbine micro-siting in complex terrain | |
Ren et al. | Numerical simulation and prediction of spatial wind field under complex terrain | |
Correia et al. | Simulation of correlated wind speed and power variates in wind parks | |
Wang et al. | Large-eddy simulations of ventilation for thermal comfort—a parametric study of generic urban configurations with perpendicular approaching winds | |
Aumond et al. | Including the drag effects of canopies: real case large-eddy simulation studies | |
Dellnitz et al. | Seasonal variability of the subpolar gyres in the Southern Ocean: a numerical investigation based on transfer operators | |
Gao et al. | Probability density function for wave elevation based on Gaussian mixture models | |
Wilson et al. | Calculation of winds disturbed by an array of fences | |
CN104035096B (en) | Vertical wind profile nonlinear inversion method based on Doppler weather radar | |
Kenjereš et al. | Modelling and simulations of turbulent flows in urban areas with vegetation | |
Kitsios et al. | Subgrid model with scaling laws for atmospheric simulations | |
Huang et al. | Multi-scale simulation of time-varying wind fields for Hangzhou Jiubao Bridge during Typhoon Chan-hom | |
Zhou et al. | Turbulence modeling for the stable atmospheric boundary layer and implications for wind energy | |
Song et al. | Wind resource assessment on complex terrain based on observations of a single anemometer | |
Karmakar et al. | Conditional simulation of non-Gaussian wind velocity profiles: Application to buffeting response of Vincent Thomas suspension bridge | |
Malavelle et al. | A method to represent subgrid‐scale updraft velocity in kilometer‐scale models: Implication for aerosol activation | |
Burlando et al. | A simple and efficient procedure for the numerical simulation of wind fields in complex terrain | |
Nielsen et al. | Simulation of inhomogeneous, non-stationary and non-Gaussian turbulent winds | |
CN104951798A (en) | Method for predicting non-stationary fluctuating wind speeds by aid of LSSVM (least square support vector machine) on basis of EMD (empirical mode decomposition) | |
Hong et al. | Stochastic Fourier spectrum model and probabilistic information analysis for wind speed process | |
Turbelin et al. | Wavelet cross-correlation analysis of wind speed series generated by ANN based models | |
Liu et al. | Probabilistic turbulence spectra of boundary-layer winds based on measurement at Jiangyin Bridge site |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150429 |