System and Method for Importance Sampling Based TimeDependent Reliability Prediction

0Associated
Cases 
0Associated
Defendants 
0Accused
Products 
5Forward
Citations 
0
Petitions 
1
Assignment
First Claim
1. A system for generating a reliability prediction for components of a vehicle, the system comprising:
 sensors electrically coupled to a data acquisition system for obtaining data related to the components from a random input process; and
a data analysis system, wherein the data analysis system comprises a computer processor electrically coupled to a computer memory, and the computer memory includes programming for the computer processor to perform the steps of;
(A) retrievably storing the data in the computer memory;
(B) characterizing the random input process;
(C) determining a decorrelation length;
(D) scaling up the standard deviation of a white noise level of the data;
(E) computing a covariance matrix of an original time series and of a scaled time series;
(F) beginning evaluation of a sample function;
(G) generating a scaled up sample function to produce an inflated domain;
(H) performing at least one of running a test or running a simulation model of the vehicle;
(I) computing a scaled vehicle response at a series of time steps until a first occurrence of a failure;
(J) when the failure occurs, computing a likelihood ratio based on an original joint probability density function and a sampling joint probability density function;
(K) determining whether an estimated vehicle response is equal to or greater than a threshold response, and when the estimated vehicle response is not equal to or greater than the threshold response, incrementing the time step and returning to the step (I), and when the estimated vehicle response is equal to or greater than the threshold response;
(L) incrementing a failure counter by 1 at the current time step;
(M) determining whether the number of the sample functions has exceeded a target number of sample functions and when the target number of sample functions is not exceeded, incrementing to the next sample evaluation and returning to the step (G), and when the target number of sample functions is exceeded;
(N) computing a safe number of the sample functions;
(O) calculating a failure rate estimation; and
(P) determining whether the failure rate estimation variance exceeds a predetermined value and the scale factor is greater than a predetermined amount, and when the failure rate estimation variance exceeds a predetermined estimation variance value and the scale factor is greater than a predetermined amount, reducing the scale factor by a predetermined amount and returning to the step (D), and when the failure rate estimation variance exceeds the predetermined estimation variance value;
(Q) providing the reliability prediction to a user, and ending the method.
1 Assignment
0 Petitions
Accused Products
Abstract
A system and a method of generating a reliability prediction for components of a vehicle. The system and the method include implementing importance sampling in dynamic vehicle systems when the vehicle is subjected to timedependent random terrain input. Alternatively, simulation data may be implemented. The system and the method include determining a decorrelation length, scaling up the standard deviation of white noise, and calculation of a likelihood ratio.
17 Citations
View as Search Results
Generating cumulative wearbased indicators for vehicular components  
Patent #
US 9,530,256 B2
Filed 12/22/2015

Current Assignee
International Business Machines Corporation

Sponsoring Entity
International Business Machines Corporation

Generating cumulative wearbased indicators for vehicular components  
Patent #
US 9,286,735 B1
Filed 09/26/2014

Current Assignee
International Business Machines Corporation

Sponsoring Entity
International Business Machines Corporation

Monitoring and planning for failures of vehicular components  
Patent #
US 9,454,855 B2
Filed 09/26/2014

Current Assignee
International Business Machines Corporation

Sponsoring Entity
International Business Machines Corporation

Integrating economic considerations to develop a component replacement policy based on a cumulative wearbased indicator for a vehicular component  
Patent #
US 9,514,577 B2
Filed 09/26/2014

Current Assignee
International Business Machines Corporation

Sponsoring Entity
International Business Machines Corporation

Method and System for Detecting Transients in Power Grids  
Patent #
US 20130158903A1
Filed 11/15/2012

Current Assignee
Mitsubishi Electric Research Laboratories

Sponsoring Entity
Mitsubishi Electric Research Laboratories

Optimized Stochastic Resonance Signal Detection Method  
Patent #
US 20100169051A1
Filed 02/22/2010

Current Assignee
Syracuse University

Sponsoring Entity
Syracuse University

SYSTEM AND METHOD FOR FEEDBACKGUIDED TEST GENERATION FOR CYBERPHYSICAL SYSTEMS USING MONTECARLO  
Patent #
US 20100094611A1
Filed 10/14/2009

Current Assignee
NEC Corporation

Sponsoring Entity
NEC Corporation

Detector and method for estimating data probability in a multichannel receiver  
Patent #
US 7,813,438 B2
Filed 09/11/2008

Current Assignee
University of Utah Research Foundation

Sponsoring Entity
University of Utah Research Foundation

LOAD CALCULATING DEVICE AND LOAD CALCULATING METHOD  
Patent #
US 20080249743A1
Filed 03/26/2008

Current Assignee
Toshiba Corporation

Sponsoring Entity
Toshiba Corporation

Materialsbased failure analysis in design of electronic devices, and prediction of operating life  
Patent #
US 20080015827A1
Filed 03/13/2007

Current Assignee
Vextec Corporation

Sponsoring Entity
Vextec Corporation

System and method for deriving a hierarchical event based database optimized for analysis of complex accidents  
Patent #
US 7,805,390 B2
Filed 02/26/2007

Current Assignee
ServiceNow Incorporated

Sponsoring Entity
International Business Machines Corporation

System and method to facilitate simulation  
Patent #
US 20050288918A1
Filed 06/24/2004

Current Assignee
Chen Thomas W., Berta Eugene

Sponsoring Entity
Chen Thomas W., Berta Eugene

Method and apparatus for predicting failure in a system  
Patent #
US 20030004679A1
Filed 01/08/2002

Current Assignee
Vextec Corporation

Sponsoring Entity
Vextec Corporation

Method and apparatus for predicting failure in a system  
Patent #
US 7,006,947 B2
Filed 01/08/2002

Current Assignee
Vextec Corporation

Sponsoring Entity
Vextec Corporation

Method and system of confidence scoring  
Patent #
US 6,353,767 B1
Filed 08/25/2000

Current Assignee
General Electric Company

Sponsoring Entity
General Electric Company

Method and system for computing path dependent probabilities of attaining financial goals  
Patent #
US 7,031,935 B1
Filed 07/31/2000

Current Assignee
Riskmetrics Group Incorporated

Sponsoring Entity
Riskmetrics Group Incorporated

Estimation method and system for complex securities using lowdiscrepancy deterministic sequences  
Patent #
US 5,940,810 A
Filed 07/30/1997

Current Assignee
Trustees Of Columbia University In The City Of New York

Sponsoring Entity


20 Claims
 1. A system for generating a reliability prediction for components of a vehicle, the system comprising:
 sensors electrically coupled to a data acquisition system for obtaining data related to the components from a random input process; and
a data analysis system, wherein the data analysis system comprises a computer processor electrically coupled to a computer memory, and the computer memory includes programming for the computer processor to perform the steps of;
(A) retrievably storing the data in the computer memory;
(B) characterizing the random input process;
(C) determining a decorrelation length;
(D) scaling up the standard deviation of a white noise level of the data;
(E) computing a covariance matrix of an original time series and of a scaled time series;
(F) beginning evaluation of a sample function;
(G) generating a scaled up sample function to produce an inflated domain;
(H) performing at least one of running a test or running a simulation model of the vehicle;
(I) computing a scaled vehicle response at a series of time steps until a first occurrence of a failure;
(J) when the failure occurs, computing a likelihood ratio based on an original joint probability density function and a sampling joint probability density function;
(K) determining whether an estimated vehicle response is equal to or greater than a threshold response, and when the estimated vehicle response is not equal to or greater than the threshold response, incrementing the time step and returning to the step (I), and when the estimated vehicle response is equal to or greater than the threshold response;
(L) incrementing a failure counter by 1 at the current time step;
(M) determining whether the number of the sample functions has exceeded a target number of sample functions and when the target number of sample functions is not exceeded, incrementing to the next sample evaluation and returning to the step (G), and when the target number of sample functions is exceeded;
(N) computing a safe number of the sample functions;
(O) calculating a failure rate estimation; and
(P) determining whether the failure rate estimation variance exceeds a predetermined value and the scale factor is greater than a predetermined amount, and when the failure rate estimation variance exceeds a predetermined estimation variance value and the scale factor is greater than a predetermined amount, reducing the scale factor by a predetermined amount and returning to the step (D), and when the failure rate estimation variance exceeds the predetermined estimation variance value;
(Q) providing the reliability prediction to a user, and ending the method.  View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10)
 sensors electrically coupled to a data acquisition system for obtaining data related to the components from a random input process; and
 11. A method of generating a reliability prediction for components of a vehicle, the method comprising the steps of:
 (A) obtaining data related to the components from a random input process and retrievably storing the data in a computer memory, and via programming stored in the computer memory implementing a computer processor to perform the steps of;
(B) characterizing the random input process;
(C) determining a decorrelation length;
(D) scaling up the standard deviation of a white noise level of the data;
(E) computing a covariance matrix of an original time series and of a scaled time series;
(F) beginning evaluation of a sample function;
(G) generating a scaled up sample function to produce an inflated domain;
(H) performing at least one of running a test or running a simulation model of the vehicle;
(I) computing a scaled vehicle response at a series of time steps until a first occurrence of a failure;
(J) when the failure occurs, computing a likelihood ratio based on an original joint probability density function and a sampling joint probability density function;
(K) determining whether an estimated vehicle response is equal to or greater than a threshold response, and when the estimated vehicle response is not equal to or greater than the threshold response, incrementing the time step and returning to the step (I), and when the estimated vehicle response is equal to or greater than the threshold response;
(L) incrementing a failure counter by 1 at the current time step;
(M) determining whether the number of the sample functions has exceeded a target number of sample functions and when the target number of sample functions is not exceeded, incrementing to the next sample evaluation and returning to the step (G), and when the target number of sample functions is exceeded;
(N) computing a safe number of the sample functions;
(O) calculating a failure rate estimation; and
(P) determining whether the failure rate estimation variance exceeds a predetermined value and the scale factor is greater than a predetermined amount, and when the failure rate estimation variance exceeds a predetermined estimation variance value and the scale factor is greater than a predetermined amount, reducing the scale factor by a predetermined amount and returning to the step (D), and when the failure rate estimation variance exceeds the predetermined estimation variance value;
(Q) providing the reliability prediction to a user, and ending the method.  View Dependent Claims (12, 13, 14, 15, 16, 17, 18, 19, 20)
 (A) obtaining data related to the components from a random input process and retrievably storing the data in a computer memory, and via programming stored in the computer memory implementing a computer processor to perform the steps of;
1 Specification
The invention described here may be made, used and licensed by and for the U.S. Government for governmental purposes without paying royalty to us.
1. Field of the Invention
The present invention generally relates to a system and method for importance sampling based timedependent reliability prediction.
2. Background Art
Conventional systems and methods for long term (i.e., timedependent) reliability prediction are typically inaccurate in some examples and computationally intensive, hence expensive, in other examples. In particular, accurate, rapid, inexpensive vehicle component long term reliability prediction can be especially problematic where the components degrade due to timedependent effects such as multiple exposures to relatively small terrain and load related forces and corrosive environment effects.
Thus, there exists a need and an opportunity for an improved system and method for long term vehicle component reliability prediction. Such an improved system and method may overcome one or more of the deficiencies of the conventional approaches.
Accordingly, the present invention may provide a system and method for importance sampling based timedependent reliability prediction.
According to the present invention, a system for generating a reliability prediction for components of a vehicle is generally provided. The system includes:
sensors electrically coupled to a data acquisition system for obtaining data related to the components from a random input process; anda data analysis system, wherein the data analysis system comprises a computer processor electrically coupled to a computer memory, and the computer memory includes programming for the computer processor to perform the steps of:(A) retrievably storing the data in the computer memory;(B) characterizing the random input process;(C) determining a decorrelation length;(D) scaling up the standard deviation of a white noise level of the data;(E) computing a covariance matrix of an original time series and of a scaled time series;(F) beginning evaluation of a sample function;(G) generating a scaled up sample function to produce an inflated domain;(H) performing at least one of running a test or running a simulation model of the vehicle;(I) computing a scaled vehicle response at a series of time steps until a first occurrence of a failure;(J) when the failure occurs, computing a likelihood ratio based on an original joint probability density function and a sampling joint probability density function;(K) determining whether an estimated vehicle response is equal to or greater than a threshold response, and when the estimated vehicle response is not equal to or greater than the threshold response, incrementing the time step and returning to the step (I), and when the estimated vehicle response is equal to or greater than the threshold response;(L) incrementing a failure counter by 1 at the current time step;(M) determining whether the number of the sample functions has exceeded a target number of sample functions and when the target number of sample functions is not exceeded, incrementing to the next sample evaluation and returning to the step (G), and when the target number of sample functions is exceeded;(N) computing a safe number of the sample functions;(O) calculating a failure rate estimation; and(P) determining whether the failure rate estimation variance exceeds a predetermined value and the scale factor is greater than a predetermined amount, and when the failure rate estimation variance exceeds a predetermined estimation variance value and the scale factor is greater than a predetermined amount, reducing the scale factor by a predetermined amount and returning to the step (D), and when the failure rate estimation variance exceeds the predetermined estimation variance value;(Q) providing the reliability prediction to a user, and ending the method.
The system wherein, the step of characterizing the random input process further comprises time series modeling of the data.
The system wherein, the step of characterizing the random input process further comprises generating an autoregressive integrated moving average (ARIMA) model of the data.
The system wherein, the step of characterizing the random input process further comprises estimating feedback parameters of the data.
The system wherein, the step of characterizing the random input process further comprises estimating a standard deviation of the white noise in the data.
The system wherein, a scaling factor in the range of 1.2 to 1.5 is implemented to inflate the standard deviation of the white noise level of the data.
The system wherein, the covariance matrix is computed via YuleWalker equations.
The system further comprising the step of storing the covariance matrix in the computer memory.
The system further comprising the step of computing a likelihood ratio.
The system further comprising the step of adding the likelihood ratio to a previous sum at the same time step.
Also according to the present invention, a method of generating a reliability prediction for components of a vehicle is provided. The method including the steps of:
(A) obtaining data related to the components from a random input process and retrievably storing the data in a computer memory, and via programming stored in the computer memory implementing a computer processor to perform the steps of:(B) characterizing the random input process;(C) determining a decorrelation length;(D) scaling up the standard deviation of a white noise level of the data;(E) computing a covariance matrix of an original time series and of a scaled time series;(F) beginning evaluation of a sample function;(G) generating a scaled up sample function to produce an inflated domain;(H) performing at least one of running a test or running a simulation model of the vehicle;(I) computing a scaled vehicle response at a series of time steps until a first occurrence of a failure;(J) when the failure occurs, computing a likelihood ratio based on an original joint probability density function and a sampling joint probability density function;(K) determining whether an estimated vehicle response is equal to or greater than a threshold response, and when the estimated vehicle response is not equal to or greater than the threshold response, incrementing the time step and returning to the step (I), and when the estimated vehicle response is equal to or greater than the threshold response;(L) incrementing a failure counter by 1 at the current time step;(M) determining whether the number of the sample functions has exceeded a target number of sample functions and when the target number of sample functions is not exceeded, incrementing to the next sample evaluation and returning to the step (G), and when the target number of sample functions is exceeded;(N) computing a safe number of the sample functions;(O) calculating a failure rate estimation; and(P) determining whether the failure rate estimation variance exceeds a predetermined value and the scale factor is greater than a predetermined amount, and when the failure rate estimation variance exceeds a predetermined estimation variance value and the scale factor is greater than a predetermined amount, reducing the scale factor by a predetermined amount and returning to the step (D), and when the failure rate estimation variance exceeds the predetermined estimation variance value;(Q) providing the reliability prediction to a user, and ending the method.
The method wherein, the step of characterizing the random input process further comprises time series modeling of the data.
The method wherein, the step of characterizing the random input process further comprises generating an autoregressive integrated moving average (ARIMA) model of the data.
The method wherein, the step of characterizing the random input process further comprises estimating feedback parameters of the data.
The method wherein, the step of characterizing the random input process further comprises estimating a standard deviation of the white noise in the data.
The method wherein, a scaling factor in the range of 1.2 to 1.5 is implemented to inflate the standard deviation of the white noise level of the data.
The method wherein, the covariance matrix is computed via YuleWalker equations.
The method further comprising the step of storing the covariance matrix in the computer memory.
The method further comprising the step of computing a likelihood ratio.
The method further comprising the step of adding the likelihood ratio to a previous sum at the same time step.
The above features, and other features and advantages of the present invention are readily apparent from the following detailed descriptions thereof when taken in connection with the accompanying drawings.
The following definitions and terminology are applied as understood by one skilled in the appropriate art.
The singular forms such as “a,” “an,” and “the” include plural references unless the context clearly indicates otherwise. For example, reference to “a material” includes reference to one or more of such materials, and “an element” includes reference to one or more of such elements.
As used herein, “substantial” and “about”, when used in reference to a quantity or amount of a material, characteristic, parameter, and the like, refer to an amount that is sufficient to provide an effect that the material or characteristic was intended to provide as understood by one skilled in the art. The amount of variation generally depends on the specific implementation. Similarly, “substantially free of” or the like refers to the lack of an identified composition, characteristic, or property. Particularly, assemblies that are identified as being “substantially free of” are either completely absent of the characteristic, or the characteristic is present only in values which are small enough that no meaningful effect on the desired results is generated.
A plurality of items, structural elements, compositional elements, materials, subassemblies, and the like may be presented in a common list or table for convenience. However, these lists or tables should be construed as though each member of the list is individually identified as a separate and unique member. As such, no individual member of such list should be considered a de facto equivalent of any other member of the same list solely based on the presentation in a common group so specifically described.
Concentrations, values, dimensions, amounts, and other quantitative data may be presented herein in a range format. One skilled in the art will understand that such range format is used for convenience and brevity and should be interpreted flexibly to include not only the numerical values explicitly recited as the limits of the range, but also to include all the individual numerical values or subranges encompassed within that range as if each numerical value and subrange is explicitly recited. For example, a size range of about 1 dimensional unit to about 100 dimensional units should be interpreted to include not only the explicitly recited limits, but also to include individual sizes such as 2 dimensional units, 3 dimensional units, 10 dimensional units, and the like; and subranges such as 10 dimensional units to 50 dimensional units, 20 dimensional units to 100 dimensional units, and the like.
With reference to the Figures, the preferred embodiments of the present invention will now be described in detail. Generally, the present invention provides an improved system and an improved method for importance sampling based timedependent reliability prediction. An example of reliability prediction for components of a vehicle that is operated on a terrain providing a random input to the vehicle is discussed below as exemplary of the present invention; however, the present invention is not limited to the example discussed. One of ordinary skill in the relevant art is assumed to have a working knowledge of conventional statistical mathematical concepts, applications, and analysis techniques, as used herein, in particular, conventional reliability computations, autoregressive integrated moving average (ARIMA) modeling, Monte Carlo simulation, importance sampling, YuleWalker equations, and the like.
Referring to
The vehicle 102 generally includes sensors 110 (e.g., load cells, accelerometers, strain gages, displacement gages, force transducers, thermocouples, profile meters, etc.) that generate data, DATA, related to the terrain, TERR, and other operating and environmental conditions to which the components of the vehicle 102 are exposed. The terrain, TERR, generally results in random inputs to the vehicle 102 (see, for example,
The data acquisition system 104 is generally electrically coupled to the sensors 110. The data acquisition system 104 generally acquires the data to be analyzed, and transmits the data, DATA, to the data analysis system 106. The data, DATA, may be transmitted wirelessly (as illustrated), via recording and subsequent downloading, or hardwire interconnection.
The data analysis system 106 generally includes a memory 120 where the data, DATA, and appropriate programming may be stored and retrieved, a processor 122 that may implement the programming stored in the memory 120 to analyze the data, DATA, that is stored in the memory 120, and an input/output (I/O) (e.g., printer, display screen, keyboard, mouse, user interface, etc.) 124. The memory 120, the processor 122, and the I/O 124 are generally electrically coupled.
The I/O 124 may provide a user ability to control the operation of the system 100 generally and, in one example, may present the reliability prediction to a user via the data analysis system 106. In other examples, the data, DATA that is processed via the data analysis system 106 may comprise historically acquired data, may comprise simulated data, and may originate from sources other than the vehicle 102 and the data acquisition system 104.
Referring to
The method 2000 may obtain (i.e., acquire, download, retrieve, etc.) data, DATA (block or step 2010). In one example, the user may measure a sample of random input terrain profile or random input load excitation via operation of the vehicle 102 on the terrain, TERR. Random input load excitation can be measured using, for example, wheel force transducers or accelerometers or other of the sensors 110.
Referring to
The method 2000 may characterize the original random input process (block or step 2020). The step 2020 comprises subblocks or substeps 2022 and 2024.
The random input process is generally characterized via timeseries modeling (the subblock or substep 2022). In one example, an autoregressive integrated moving average (ARIMA) model may be implemented. As is known to one of skill in the art, when one of the terms is zero, AR, I or MA are usually dropped. For example, an I(1) model is ARIMA(0,1,0), a MA(1) model is ARIMA(0,0,1), and so forth.
For the substep 2022, the data, (e.g., DATAa), is considered the result of a random process (e.g., as illustrated on the plot of
where μ is the temporal mean of the process, ε_{i}≡N(0, σ_{e}^{2}) is Gaussian white noise and φ_{1}, φ_{2}, . . . φ_{p }are feedback parameters. All model parameters, μ, σ_{e}^{2}, φ_{1}, φ_{2}, . . . φ_{p }are to be estimated.
Estimate the model parameters (the subblock or substep 2024). As understood by one of skill in the art, different order AR models can be generated to determine the best fit. For an AR(p) model, the variance σ_{e}^{2 }of the Gaussian white noise is determined from
where γ(0) is the variance of the random process, and ρ_{p }is the value of the autocorrelation function at time lag τ=p·Δt. Similar expressions exist for higher order AR models.
After the feedback parameters are estimated, a residual series E(t)=x(t)−ˆ(t) is formed as the difference between the actual x(t) and the estimated ˆ(t) processes and statistical tests are performed to make sure that the random variables E_{t }and E_{t+τ} are uncorrelated for every τ.
When not known, the appropriate AR model type can be identified by a user by visually inspecting the plots of the autocorrelation and the partial sample autocorrelation functions for different lags (multiples of Δt; see,
where {circumflex over (σ)} is the estimated standard deviation of the random process. In the above equation, an unbiased estimation of ρ(τ) if is replaced by (n−h)^{−1}. For convenience however, the n^{−1 }term may be implemented.
The partial autocorrelation of lag h represents the autocorrelation between X_{i }and X_{i+τ} with the linear dependence of X_{i+1 }through X_{i+τ1 }removed. The partial autocorrelation is representative of the autocorrelation between X_{i }and X_{i−τ} that is not accounted for by lags 1 to τ−1, inclusive. The partial autocorrelation is generally useful in identifying the order of an autoregressive model. For an AR(p) model, it is zero for lags greater or equal to p+1. After the order p of the model is identified, the φ's and μ are estimated either by using the YuleWalker equations or alternatively, by minimizing
Statistical tests may be performed to ensure the goodness of fit. When the data, DATA, includes results from a known environment terrain, TERR; using an AR(3) (e.g., where ARIMA is used for p=3) autoregressive timeseries model may represent the random road process. For the example shown, the following three parameters of the model
were estimated. The standard deviation σ_{ε}=0.5132 of the zeromean residual process, ε_{i }was also estimated. The AR(3) model is then expressed as
Referring to
For a stationary process, the autocorrelation function generally decays rapidly, either exponentially or by overshooting into the negative region before settling down. When the autocorrelation function decay is exponential, all of the feedback parameters are generally positive, the number of feedback parameters may be sufficient to estimate the shape of the autocorrelation function. Therefore, a decorrelation length d=p which is equal to the order of the AR(p) model may be implemented because the partial autocorrelation function generally becomes insignificant after p lags.
When the autocorrelation function overshoots the zero axis, the oscillations with increasing lag generally indicate that there is at least one negative feedback parameter. Depending upon which feedback parameter is negative, the shape of the autocorrelation function can vary (change). As such, the shape of autocorrelation function is generally not determined based only on the order of the AR(p) model.
In the case illustrated on
Scaleup the standard deviation of the white noise to generate an inflated input domain (block or step 2040). The standard deviation of the white noise ε^{s}=N(0, σ_{s}^{2}) as σ_{s}=fσ_{e }is generally scaled up such that f≈1.2 to 1.5 to generate an inflated random input excitation using the timeseries model. Generally as a first estimate, implement as the upper value, f=1.5.
Compute the covariance matrix of the input timeseries (block or step 2050). The step 2050 comprises subblocks or substeps 2052 and 2054. YuleWalker equations may be implemented to compute the covariance matrix of both original (Σ) and sampling distribution (Σ_{S}) (substep 2052) using the correlation coefficients from the equation below.
where m=1, 2 . . . k and ρ_{m }is the correlation coefficient at lag m.
Store the covariance matrix in the memory 120 of the computer or data base 106, for later retrieval (substep 2054).
The terrain or the random process input analysis has been completed. Evaluation of response of the vehicle 102 to the input is generally conducted next.
Set a first sample function evaluation, N=1 (block or step 2060).
Generate a scaledup input excitation sample function (block or step 2070). Implement the previously calculated scaledup standard deviation, σ_{s }of the timeseries model. By scalingup the road excitation, an inflated sampling distribution may be generated (i.e., produced, calculated, etc.), which generally produces a large number of firstpassage failures. The large number of firstpassage failures generally advantageously decreases the required number of samples without sacrificing accuracy when compared to conventional approaches.
The scaledup excitation sample function then becomes:
where μ is the temporal mean of the process, ε_{i}≡N(0,σ_{s}^{2})
The step 2070 is generally similar to the step 2020; however, generally implemented with the higher standard deviation, σ_{s}=fσ_{e }(from the step 2040) of the Gaussian white noise, ε_{i}≡N(0,σ_{s}^{2}) (from the step 2020) while keeping all other estimated parameters same as in the step 2020.
Conduct (e.g., run, perform, etc.) a test (see discussion in connection with
Referring to
The test process as illustrated on
Referring to
In the simulation embodiment of
Referring to
Compute the vehicle response, S^{s}(t_{i}) at every time step until the first occurrence of the failure i.e. S^{s}(t_{i})≧S_{threshold }where S_{threshold }is the maximum acceptable level of the response (block or step 2090). The vehicle response S^{s}(t_{i}) may be such as vehicle acceleration, stress or strain in the component.
Referring to
Referring to
Referring to
When a failure occurs, compute the Likelihood Ratio (block or step 2100). The step 2100 may include substeps (or subblocks) 2102 and 2104.
Likelihood Ratio (e.g., the substep 2102):
where
The joint density f_{X}(x) is calculated using the k=d+1 normal random variables of the random vector X={x_{1}, X_{i−1}, . . . , X_{i−d}} as,
where μ={μ_{i}, μ_{i−1 }. . . μ_{i−d}}={
Similarly, the sampling density is given by
where μ^{S }and Σ^{S }are the mean vector and the covariance matrix associated with the inflated random input vector.
The likelihood ratio is added to the previous sum at the given time instant (e.g., the substep 2104) as:
Determine whether the condition
is satisfied, where x_{f }is the value of an inflated response at failure (decision block or step 2110).
The safe sample functions are generally calculated from the original environment so that more safe sample functions remain in the population at later times. The condition of safe sample functions remaining in the population is generally achieved by discarding a sampling sample function only when the condition
is satisfied, where x_{f }is the value of an inflated response at failure. The response in the original environment may be approximated by scaling down the inflated response using the ratio of original and sampling standard deviations of the residual process.
When the sample function is not discarded (i.e., the NO leg of the decision block 2110), the time is incremented by one step (block or step 2112); and the likelihood ratio is again computed for the next occurrence of the failure (i.e., the step 2090 is again performed). The steps 2090, 2100, and 2110 may be repeated until the sample function is discarded or until last time step in the data, DATA, is reached (completed).
When the sample function is discarded, the next (subsequent) sample function is generally evaluated from step 2070 onwards.
When the sample function is discarded (i.e., the YES leg of the decision block 2110), increment number of failures by 1 at the given time step (block or step 2120):
Failure counter, N_{f}(t_{i})=N_{f}(t_{i})+1, where the failure counter is generally an approximation of the number of failures in the original (i.e., not scaled up) domain.
Determine whether the number of sample functions has exceeded a target number of sample function evaluations (decision block or step 2130). When the number of sample functions has not exceeded the target number of sample function evaluations (i.e., the NO leg of the decision block 2130), increment to the next sample evaluation (block or step 2132), and return to step 2070. When the number of sample functions has exceeded the target number of sample function evaluations (i.e., the YES leg of the decision block 2130), compute the safe number of sample functions N_{S }(block or step 2140).
The safe number of sample functions N_{S }is generally computed at every step 2140 by subtracting the failed number of samples from the previous safe number of sample functions:
Estimate First Passage Failure Rate (block or step 2150).
Estimated first passage failure rate,
Determine whether the variance in the estimated failure rate exceeds a predetermined value (e.g., a predetermined variance) and the scale factor is greater than a predetermined scale factor (in the example described, f>1.2) (decision block or step 2160). When the variance in the estimated failure rate exceeds the predetermined value and the scale factor is greater than the predetermined scale factor (i.e., the YES leg of the decision block 2160), reduce the scale factor by a predetermined amount (e.g., for the example described, 0.1) (block or step 2162), and return to the block 2040.
When the variance in the estimated failure rate does not exceed the predetermined value and the scale factor is greater than the predetermined scale factor (i.e., the NO leg of the decision block 2160), provide the reliability prediction to the user (block or step 2170), and end the process 2000 (block or step 2180).
The embodiment demonstrated through reliability prediction analysis via the method 2000 of the quarter vehicle example of
Referring to
Similar accuracy levels are also demonstrated for the higher vehicle threshold response of 2.65 G (see,
As is apparent then from the above detailed description, the present invention may provide an improved system 100 and an improved method 2000 for generating a reliability prediction for components of a vehicle. The method 2000 includes implementing importance sampling in dynamic vehicle systems when the vehicle (e.g., the vehicle 102) is subjected to timedependent random terrain input (e.g., the terrain, TERR).
Other example systems that may advantageously implement the method 2000, may include any appropriate timedependent random input data having a large number of data points to consider when making a prediction. Such examples may include finance, econometrics, and biomedical engineering, and the like.
Various alterations and modifications will become apparent to those skilled in the art without departing from the scope and spirit of this invention and it is understood this invention is limited only by the following claims.