140 likes | 284 Views
Probabilistic Analysis using FEA. A. Petrella. What is Probabilistic Analysis. All input parameters have some uncertainty What is the uncertainty in outcome metrics? How sensitive are outcomes to different inputs
E N D
Probabilistic Analysis using FEA A. Petrella
What is Probabilistic Analysis • All input parameters have some uncertainty • What is the uncertainty in outcome metrics? • How sensitive are outcomes to different inputs • Which inputs are most important and how can we design for a specific probability of performance?
What is Probabilistic Analysis Outcome Probabilities & Sensitivities Validated Deterministic Model Model Input Uncertainties Probability Tissue Properties Performance Metric External Loads Response and Failure Prediction Device Placement Sensitivity Factors
Probabilistic Methods • Monte Carlo (MC) is the simplest prob method… input distributions randomly sampled to form trials • MC is robust and will always converge, but this usually requires many thousands of trials • It may be impractical to perform 1000’s of trials with an FE model that requires hours for one solution • There are more advanced methods that require fewer trials and many modern programs implement these methods… e.g., ANSYS uses DOE + Response Surface
Prob… an example with Excel Random variables, normally distributed h = 400 ± 20 mm b = 100 ± 5 mm P = 1000 ± 50 N E = 200 ± 10 GPa P h L = 2400 mm b
Standard Normal Distribution CDF PDF m = 0 s = 1
Standard Normal Distribution • Normal (m=0, s=1) • Standard normal variate • (Note: Halder uses S) • All normal distributions can be simply transformed to the standard normal distribution
Back to the Beam Example… 500 MC To get the 10% lower and 90% upper bounds… Use Excel functions: “large()” and “small()”
Beam Example in ANSYS • ANSYS uses the term…“Sig Sigma Analysis”…this is most likely marketing since 6s is popular in industry • Prob trials are taken from a response surface (quadratic polynomial regression) built on a results from a DOE • This is how ANSYS avoids 1000’s of trials required for a brute force MC
Beam Example in ANSYS - Sensitivity Sensitivity factors are the components of a unit vector in the direction of the function gradient…(i.e., stress = f(h,b,P,E)) …then sqrt(sum(si2)) = 1 sh sb sP sE sh sb sP sE
How does Prob Compare? • Provides information on sensitivities similar to DOE and Response Surface methods • Prob provides insight into how uncertainty in your input parameters will affect outcome metrics • Allows you to design for probability of specific outcomes… e.g., 90% upper bound