CN108763611A - A kind of wing structure random eigenvalue analysis method based on probabilistic density evolution - Google Patents
A kind of wing structure random eigenvalue analysis method based on probabilistic density evolution Download PDFInfo
- Publication number
- CN108763611A CN108763611A CN201810261867.0A CN201810261867A CN108763611A CN 108763611 A CN108763611 A CN 108763611A CN 201810261867 A CN201810261867 A CN 201810261867A CN 108763611 A CN108763611 A CN 108763611A
- Authority
- CN
- China
- Prior art keywords
- equation
- eigenvalue
- formula
- probability density
- characteristic value
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Complex Calculations (AREA)
Abstract
The wing structure random eigenvalue analysis method based on probabilistic density evolution that the invention discloses a kind of, belongs to field of structural design.It fully considers uncertain factor present in large and complex structure, under conditions of sample information abundance, quantification characterization is carried out to uncertain parameter using random device.The finite element equation for establishing free vibration of structures is converted into generalized eigenvalue problem by modal coordinate, and is derived to Eigenvalue Sensitivity.On this basis, characteristic value probability density evolution equation is established, Eigenvalue Sensitivity is introduced into equation, obtains the characteristic value probability density evolution equation based on sensitivity analysis.It is canonical form by introducing new variables by characteristic value probability density evolution equation abbreviation, reducing format using finite difference method and total variance solves characteristic value probability density function.It is preferable that numerical result shows that the characteristic value probability density function that the method for the present invention obtains coincide with monte carlo method, and can significantly reduce and calculate the time.
Description
Technical field
The present invention relates to field of structural design, more particularly to a kind of wing structure random character based on probabilistic density evolution
It is worth analysis method.
Background technology
Structural natural frequencies and characteristic value can provide the information of characterization structural dynamic characteristics, the especially intrinsic frequency of low order
Rate and characteristic value are usually used to the kinetic stability of evaluation structure.Meanwhile considering the structure optimization of intrinsic frequency and characteristic value
It is also very promising in field of structural design.Therefore, structural natural frequencies and characteristic value are generally viewed as in structure design
Important indicator parameter.Structural Eigenvalue Problem has just had received widespread attention since the eighties in last century, up to the present proposes
A variety of methods for solving Structural Eigenvalues, and be successfully applied to Structural Reanalysis, sensitivity analysis and Modifying model etc.
In Practical Project problem.Under the scope of certainty, structural natural frequencies and characteristic value can include mass matrix by that will solve
It is obtained with the generalized eigenvalue equation of stiffness matrix.But there are multi-source uncertain factors in practical structures, such as:(1) material
The uncertainty for expecting parameter makes engineering material since the factors such as manufacturing environment, technical conditions, the multiphase feature of material influence
Elasticity modulus, Poisson's ratio, mass density have uncertainty;(2) uncertainty of geometric dimension, due to manufacture and installation error
Make geometrical scale such as thickness, cross-sectional area etc. that there is uncertainty, the presence of these uncertain factors to make structure feature
Value problem solving has difficulties.
For uncertain caused by material and geometric parameter, generally using random device to the uncertain factor amount of progress
Change characterization, utilizes the distribution characteristics of stochastic variable characterising parameter.To the Structural Eigenvalue Solve problems containing stochastic variable, existing solution
Certainly approach can be divided into three categories:(1) Monte-carlo Simulation Method;(2) response surface method of deploying;(3) Moment equation method.However, covering
Special Carlow analogy method needs to carry out a large amount of sample point analysis, and it is larger to expend computing resource;Response surface method of deploying and square side
Journey method can only handle small uncertain problem.The Matrix perturbation and direct method of analysis of variance proposed in recent years can also solve
Structure random parameters problem, but the above method can only obtain the first-order reliability method information of characteristic value.It is special when needing to obtain structure
When value indicative probability density function, there is presently no effective analysis methods to solve, and limits complicated knot to a certain extent
The development of structure stability analysis technology.In conclusion there is an urgent need for developing one kind, can quickly, accurately to solve structure random parameters general
The new method of rate density function, to overcome conventional method to calculate drawback of long duration, that precision is low.
Invention content
The technical problem to be solved in the present invention is:For conventional process structure random parameters problem method computational efficiency
The problems such as low, characteristic value probability density function is difficult to obtain proposes that a kind of wing structure based on probabilistic density evolution is special at random
Value indicative analysis method.This method fully considers uncertain factor present in large and complex structure, using random device to not true
Determine parameter and carries out quantification characterization.Under the conditions of certainty, free vibration of structures equation is established, is converted by modal coordinate
For generalized eigenvalue problem, and numerical solution is carried out to Eigenvalue Sensitivity.On this basis, characteristic value probability density is established
EVOLUTION EQUATION, and Eigenvalue Sensitivity is introduced into constructed equation, it is close to obtain the characteristic value probability based on sensitivity analysis
Spend EVOLUTION EQUATION.It is canonical form by introducing new variables by characteristic value probability density evolution equation abbreviation, using finite difference
Method and total variance, which reduce format, can obtain characteristic value probability density function.
The present invention solve the technical solution that uses of above-mentioned technical problem for:A kind of wing structure based on probabilistic density evolution
Random eigenvalue analysis method, includes the following steps:
Step (1), the finite element equation for establishing free vibration of structures:
In formula, M is mass matrix, and K is stiffness matrix, and q is generalized coordinates,For generalized acceleration;
Step (2) introduces modal vector u, utilizesFree vibration finite element equation can be turned
Turn to generalized eigenvalue equation:
K { u }-λ M { u }=0 (2)
Step (3), under the conditions of certainty, Structural Eigenvalue λ can be obtained by following determinant equation:
| K- λ M |=0 (3)
Step (4), Structural Eigenvalue λ can be expressed as about the sensitivity of structural parameters:
In formula, λiFor the i-th rank characteristic value, uiFor the i-th rank feature vector, b is structural parameters;
Step (5) establishes Structural Eigenvalue probability density evolution equation, is expressed as form:
In formula, Θ=(Θ1,...,Θs) it is that s ties up random uncertain parameter, p λ Θ (λ, θ, b) are that the joint of (λ, Θ) is general
Rate density function, N are characterized value quantity;When only one characteristic value (when taking N=1) of consideration, formula (5) can be rewritten as:
Step (6), in the domain of variation Ω of uncertain parameter Θ, equably take NtotalA sample point, is denoted as:
Θq(q=1 ..., Ntotal), and it is N that domain of variation Ω, which is divided,totalA subdomain, is denoted as Ωq(q=1 ...,
Ntotal);
Step (7), by equation (6) in subdomain ΩqInterior integral can obtain:
Step (8), by exchange integral and sequence of solving the derivation, can be by formula (7) abbreviation:
In formula,For the probability density function corresponding to q-th of sample point;
Eigenvalue Sensitivity is introduced into equation (8) by step (9), can obtain the probability of feature based value sensitivity
Density evolution equation:
Step (10) introduces new parameter z=λ (Θq, b) and b, substituting into equation (9) can obtain:
Step (11) determines that primary condition is:
In formula, δ is Dirac function,
Equation (10) is rewritten as following form by step (12):
In formula, a=λ (Θq,b);
Step (13) can obtain following difference scheme using finite difference method and total variance reduction format:
In formula,zm=m Δs z (m=0, ± 1 ...), bk=k Δs b (k=0,1 ...),
For current limiter,WithIt is represented by:
Step (14), will be in NtotalIt is calculated at a sample pointSummation, can obtain:
Step (15) takes bk=1, then the expression formula of Structural Eigenvalue probability density function can be obtained:
Wherein, in the step (4), the derivation of Eigenvalue Sensitivity is:
The equation of step (2) is expressed as component form, i.e.,:
K{ui}-λiM{ui}=0 (17)
In formula, λiFor the i-th rank characteristic value, uiFor the i-th rank feature vector.Local derviation is asked to structural parameters b at equation (17) both ends
Number, can obtain following relationship:
In equation (18) both sides while premultiplicationAnd it carries out transposition and converts and can obtain:
Wherein, in the step (7), the concrete mode of exchange integral and sequence of solving the derivation is:
Wherein, in the step (11), can be by primary condition is discrete:
In formula, zm=m Δs z (m=0, ± 1 ...), Δ z indicate the size of mesh opening in the directions z;
Wherein, in the step (13), the difference scheme condition of convergence to be met is:
The beneficial effects of the invention are as follows:
The present invention proposes the Probability Density Evolution Method based on sensitivity, can divide structure random parameters
Analysis obtains Structural Eigenvalue probability density function.When the mean value and standard deviation of given input parameter, what the method for the present invention obtained
The characteristic value probability density function that characteristic value probability density function is obtained with monte carlo method coincide preferably, has higher number
It is worth precision.It, can be in the premise item for fully ensuring that computational accuracy when being analyzed complicated wing structure using the method for the present invention
Under part, the required repeated sample point analysis of traditional Monte Carlo Analogue Method is avoided, is greatly reduced and calculates the time, improves and calculates
Efficiency.The method of the present invention is of great significance for shortening wing structure design cycle, lift structure stability.
Description of the drawings
Fig. 1 is the finite element model of wing structure;
Fig. 2 is the first rank characteristic value probability density function;
Fig. 3 is second-order characteristic value probability density function;
Fig. 4 is third rank characteristic value probability density function;
Fig. 5 is fourth order characteristic value probability density function;
Fig. 6 is the 5th rank characteristic value probability density function;
Fig. 7 is the method implementation process of the present invention.
Specific implementation mode
Hereinafter reference will be made to the drawings, and the design example of the present invention is described in detail.It should be appreciated that selected example only for
Illustrate the present invention, rather than limits the scope of the invention.
(1) using HIRENASD wing structures as object, structural finite element model is established, as shown in Figure 1;
(2) wing structure unit material property parameters are given, as shown in table 1;
1 wing structure material properties parameter of table
In table 1, E is elasticity modulus, and G is shear model, and μ is Poisson's ratio, and ρ is density;
(3) it is uncertain parameter, Normal Distribution, mean value and standard deviation such as table 2 to choose elastic modulus E and density p
It is shown;
The mean value and standard deviation of 2 material properties parameter of table
(4) constant interval of E and ρ [+6 σ of μ 6 σ, μ] is divided into 20 subintervals, then the E and ρ of subinterval boundary are:
In this way, sample point (the E formedi,ρj) share 441;
(5) q-th of sample point (E is takeni,ρj)q, according to (Ei,ρj)qMaterial properties are set, are obtained using business software Nastran
It takes to being applied to sample point (Ei,ρj)qFirst five rank eigenvalue λ1~λ5;
(6) it establishes respectively and is directed to Structural Eigenvalue λ1~λ5Probability density evolution equation:
(7) primary condition can discrete be:
In formula,
(8) finite difference scheme is set as:
In formula,zm=m Δs z (m=0, ± 1 ...), bk=k Δs b (k=0,1 ...),
(9) current limiterIt is set as:
(10) using by finite difference method, meeting the condition of convergenceUnder the premise of, it can be with
It obtains corresponding to sample point (Ei,ρj)qProbability density function
(11) step (5)~(10) are repeated, the corresponding probability density function of all 441 sample points is calculated
Being summed can obtain:
(12) b is takenk=1, then it can calculate structure first five rank eigenvalue λ1~λ5Probability density function expression formula be:
(13) Monte-carlo Simulation Method is utilized, the sample point of N=1000 Normal Distribution is taken, by each sample
This character pair value λ1~λ5Calculating, obtain λ1~λ5Probability density function;
(14) both the above method obtains λ1~λ5Probability density function as shown in figures 2-6, result illustrates the present invention in figure
The result that method obtains is coincide preferable with monte carlo method result;
(15) according to λ1~λ5Probability density function calculates separately mean value (μ) and standard deviation (σ), and the results are shown in Table 3:
3 result of calculation of table compares
From table 3 it can be seen that the max calculation error of two methods is no more than 4%, illustrate that the method for the present invention precision is preferable;
(16) two methods calculating total time-consuming is respectively:TThe method of the present invention=327338.70s, TMonte carlo method=614400.23s.
Time comparing result show the method for the present invention can be reduced by about 50% calculating take, to significantly improve wing structure feature
It is worth the computational efficiency of probability density function.
In conclusion the present invention proposes a kind of structure random eigenvalue analysis method based on probabilistic density evolution.First
Free vibration of structures equation is established, generalized eigenvalue problem is converted by modal coordinate.On this basis, feature is established
It is worth probability density evolution equation, and Eigenvalue Sensitivity is introduced into constructed equation, obtains the spy based on sensitivity analysis
Value indicative probability density evolution equation, characteristic value probability density letter can be obtained by reducing format using finite difference method and total variance
Number.Numerical result shows the characteristic value that the characteristic value probability density function that the method for the present invention obtains is obtained with monte carlo method
Probability density function coincide preferable, and can significantly reduce and calculate the time, to for wing structure Analysis of Vibration Characteristic and
Stability Design provides new approaches.
The specific steps that the above is only the present invention, are not limited in any way protection scope of the present invention, expansible to answer
For Aircraft structural design field, any technical scheme formed by adopting equivalent transformation or equivalent replacement, all falls within this hair
Within bright rights protection scope.
Claims (5)
1. a kind of wing structure random eigenvalue analysis method based on probabilistic density evolution, it is characterised in that:Realize step such as
Under:
Step (1), the finite element equation for establishing wing structure free vibration:
In formula, M is mass matrix, and K is stiffness matrix, and q is generalized coordinates,For generalized acceleration;
Step (2) introduces modal vector u, utilizesFree vibration finite element equation can be converted to extensively
Adopted eigenvalue equation:
K { u }-λ M { u }=0 (2)
Step (3), under the conditions of certainty, Structural Eigenvalue λ can be obtained by following determinant equation:
| K- λ M |=0 (3)
Step (4), Structural Eigenvalue λ can be expressed as about the sensitivity of structural parameters:
In formula, λiFor the i-th rank characteristic value, uiFor the i-th rank feature vector, b is structural parameters;
Step (5) establishes Structural Eigenvalue probability density evolution equation, is expressed as form:
In formula, Θ=(Θ1,...,Θs) it is that s ties up random uncertain parameter, pλΘ(λ, θ, b) is the joint probability density of (λ, Θ)
Function, N are characterized value quantity;When only one characteristic value (when taking N=1) of consideration, formula (5) can be rewritten as:
Step (6), in the domain of variation Ω of uncertain parameter Θ, equably take NtotalA sample point, is denoted asAnd Ω points are N by domain of variationtotalA subdomain, is denoted as Ωq(q=1 ..., Ntotal);
Step (7), by equation (6) in subdomain ΩqInterior integral can obtain:
Step (8), by exchange integral and sequence of solving the derivation, can be by formula (7) abbreviation:
In formula,For the probability density function corresponding to q-th of sample point;
Eigenvalue Sensitivity is introduced into equation (8) by step (9), can obtain the probability density of feature based value sensitivity
EVOLUTION EQUATION:
Step (10) introduces new parameter z=λ (Θq, b) and b, substituting into equation (9) can obtain:
Step (11) determines that primary condition is:
In formula, δ is Dirac function,
Equation (10) is rewritten as following form by step (12):
In formula, a=λ (Θq,b);
Step (13) can obtain following difference scheme using finite difference method and total variance reduction format:
In formula,zm=m Δs z (m=0, ± 1 ...), bk=k Δs b (k=0,1 ...),For current limiter,WithIt is represented by:
Step (14), will be in NtotalIt is calculated at a sample pointSummation, can obtain:
Step (15) takes bk=1, then the expression formula of Structural Eigenvalue probability density function can be obtained:
2. a kind of wing structure random eigenvalue analysis method based on probabilistic density evolution according to claim 1,
It is characterized in that:In the step (4), the derivation of Eigenvalue Sensitivity is:
The equation of step (2) is expressed as component form, i.e.,:
K{ui}-λiM{ui}=0 (17)
In formula, λiFor the i-th rank characteristic value, uiFor the i-th rank feature vector, partial derivative is asked to structural parameters b at equation (17) both ends,
It can obtain following relationship:
In equation (18) both sides while premultiplicationAnd transplant and can obtain:
3. a kind of wing structure random eigenvalue analysis method based on probabilistic density evolution according to claim 1,
It is characterized in that:In the step (7), the concrete mode of exchange integral and sequence of solving the derivation is:
4. a kind of wing structure random eigenvalue analysis method based on probabilistic density evolution according to claim 1,
It is characterized in that:In the step (11), can be by primary condition is discrete:
In formula, zm=m Δs z (m=0, ± 1 ...), Δ z indicate the size of mesh opening in the directions z.
5. a kind of wing structure random eigenvalue analysis method based on probabilistic density evolution according to claim 1,
It is characterized in that:In the step (13), the difference scheme condition of convergence to be met is:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810261867.0A CN108763611B (en) | 2018-03-28 | 2018-03-28 | Wing structure random eigenvalue analysis method based on probability density evolution |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810261867.0A CN108763611B (en) | 2018-03-28 | 2018-03-28 | Wing structure random eigenvalue analysis method based on probability density evolution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108763611A true CN108763611A (en) | 2018-11-06 |
CN108763611B CN108763611B (en) | 2022-03-15 |
Family
ID=63980444
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810261867.0A Active CN108763611B (en) | 2018-03-28 | 2018-03-28 | Wing structure random eigenvalue analysis method based on probability density evolution |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108763611B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112800650A (en) * | 2021-01-27 | 2021-05-14 | 重庆大学 | Structural time-varying reliability analysis method considering normal distribution discrete initial condition |
CN114186395A (en) * | 2021-11-23 | 2022-03-15 | 大连理工大学 | Engineering structure system reliability rapid calculation method considering multiple randomness |
CN114186447A (en) * | 2021-11-11 | 2022-03-15 | 大连理工大学 | Non-invasive random finite element method for random analysis of plate structure |
CN114429060A (en) * | 2021-12-02 | 2022-05-03 | 中国兵器科学研究院宁波分院 | Method for assessing structure dislocation failure and service life prediction in fatigue vibration |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050072234A1 (en) * | 2003-05-20 | 2005-04-07 | Weidong Zhu | System and method for detecting structural damage |
CN105183703A (en) * | 2015-10-27 | 2015-12-23 | 北京航空航天大学 | Complex mode random eigenvalue direct variance calculation method based on matrix perturbation theory |
CN105718634A (en) * | 2016-01-15 | 2016-06-29 | 北京航空航天大学 | Airfoil robust optimization design method based on non-probability interval analysis model |
CN107633118A (en) * | 2017-08-28 | 2018-01-26 | 西北工业大学 | Wing structure Parameter Sensitivity Analysis method based on random interval mixed model |
-
2018
- 2018-03-28 CN CN201810261867.0A patent/CN108763611B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050072234A1 (en) * | 2003-05-20 | 2005-04-07 | Weidong Zhu | System and method for detecting structural damage |
CN105183703A (en) * | 2015-10-27 | 2015-12-23 | 北京航空航天大学 | Complex mode random eigenvalue direct variance calculation method based on matrix perturbation theory |
CN105718634A (en) * | 2016-01-15 | 2016-06-29 | 北京航空航天大学 | Airfoil robust optimization design method based on non-probability interval analysis model |
CN107633118A (en) * | 2017-08-28 | 2018-01-26 | 西北工业大学 | Wing structure Parameter Sensitivity Analysis method based on random interval mixed model |
Non-Patent Citations (3)
Title |
---|
QIU ZHIPING,QIU HECHEN: "A direct-variance-analysis method for generalized stochastic eigenvalue problem based on matrix perturbation theory", 《SCIENCE CHINA TECHNOLOGICAL SCIENCES》 * |
YUTING DAI,CHAO YANG: "Smolyak-Grid-Based Flutter Analysis with the Stochastic Aerodynamic Uncertainty", 《DISCRETE DYNAMICS IN NATURE AND SOCIETY》 * |
王小兵 等: "基于随机因子法的随机结构动力特性分析研究", 《振动与冲击》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112800650A (en) * | 2021-01-27 | 2021-05-14 | 重庆大学 | Structural time-varying reliability analysis method considering normal distribution discrete initial condition |
CN112800650B (en) * | 2021-01-27 | 2023-05-23 | 重庆大学 | Structural time-varying reliability analysis method considering normal distribution discrete initial conditions |
CN114186447A (en) * | 2021-11-11 | 2022-03-15 | 大连理工大学 | Non-invasive random finite element method for random analysis of plate structure |
CN114186447B (en) * | 2021-11-11 | 2024-09-17 | 大连理工大学 | Non-invasive random finite element method for random analysis of plate structure |
CN114186395A (en) * | 2021-11-23 | 2022-03-15 | 大连理工大学 | Engineering structure system reliability rapid calculation method considering multiple randomness |
CN114186395B (en) * | 2021-11-23 | 2024-08-02 | 大连理工大学 | Engineering structure system reliability quick calculation method considering multiple randomness |
CN114429060A (en) * | 2021-12-02 | 2022-05-03 | 中国兵器科学研究院宁波分院 | Method for assessing structure dislocation failure and service life prediction in fatigue vibration |
CN114429060B (en) * | 2021-12-02 | 2022-12-27 | 中国兵器科学研究院宁波分院 | Method for examining structure dislocation failure and service life prediction in fatigue vibration |
Also Published As
Publication number | Publication date |
---|---|
CN108763611B (en) | 2022-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108763611A (en) | A kind of wing structure random eigenvalue analysis method based on probabilistic density evolution | |
CN105843073B (en) | A kind of wing structure aeroelastic stability analysis method not knowing depression of order based on aerodynamic force | |
CN104899380B (en) | A kind of Reliability of Slope Stability Sensitivity Analysis based on Monte Carlo simulation | |
CN107133406B (en) | Rapid search method for static voltage stability domain boundary of power system | |
Preece et al. | Probabilistic small-disturbance stability assessment of uncertain power systems using efficient estimation methods | |
CN107123994B (en) | Linear solving method of interval reactive power optimization model | |
WO2020248228A1 (en) | Computing node load prediction method in a hadoop platform | |
Matarazzo et al. | STRIDE for structural identification using expectation maximization: Iterative output-only method for modal identification | |
Slivinski et al. | A hybrid particle–ensemble Kalman filter for Lagrangian data assimilation | |
US8180619B2 (en) | System and method for digital effects analysis | |
CN107563550A (en) | A kind of Optimal Configuration Method of the real-time distributed state estimation of power distribution network based on PMU and PMU | |
CN102946253A (en) | Data sampling method and data sampling system, as well as method and system applying data sampling method and data sampling system to parameter identification | |
Leiteritz et al. | How to avoid trivial solutions in physics-informed neural networks | |
Sajadi et al. | Identification of transient stability boundaries for power systems with multidimensional uncertainties using index-specific parametric space | |
CN103885867A (en) | Online evaluation method of performance of analog circuit | |
Liu et al. | Parameter calibration in wake effect simulation model with stochastic gradient descent and stratified sampling | |
CN112365044A (en) | Tunnel face failure probability prediction method based on k nearest neighbor algorithm and support vector machine | |
Yao et al. | Reduced-order modeling for limit-cycle oscillation using recurrent artificial neural network | |
CN108205713A (en) | A kind of region wind power prediction error distribution determination method and device | |
CN109446557A (en) | A kind of random aeroelastic system method for analyzing stability based on probabilistic density evolution | |
Simiu et al. | Wind load factors for use in the wind tunnel procedure | |
CN110222098A (en) | Electric power high amount of traffic abnormality detection based on flow data clustering algorithm | |
CN111859303B (en) | Soil humidity fusion method and system based on dynamic Bayesian average | |
CN115687854B (en) | High-precision soil sample parameter measuring method and system thereof | |
Liu et al. | Model safety standard of railway signaling system for assessing the conformity of safety critical software based weighted factors analysis |
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 |