Probabilistic Design and Reliability | Module 1 of 10
Why Deterministic Design Is Not Enough
A nominal safety factor says what happens at one assumed point. Reliability asks how often the real population of parts, loads, strengths, and operating conditions crosses a failure boundary.
Readiness check
Tick only what you can do closed-notes before using this module for design decisions.
- Explain a factor of safety as resistance divided by load effect.
- Read a mean and standard deviation from a small data set.
- Identify a mechanical failure mode for a shaft, beam, or fastener.
- Separate a design requirement from a calculation assumption.
- Use units consistently in stress and strength calculations.
The core idea
A deterministic design check uses one load, one strength, and one geometry. A probabilistic design check treats those quantities as populations, defines failure, and estimates the probability that a real part falls into the failure region.
n = R_nominal / S_nominalg(X) = resistance(X) - load_effect(X)P_f = P[g(X) <= 0], reliability = 1 - P_fSafety factors are still useful because they give a quick margin against a selected design case. They do not, by themselves, say how likely failure is. Two shafts can share the same nominal safety factor while one has tightly controlled material and load variation and the other has broad scatter. The first can have a much lower failure probability even though the deterministic number matches. Reliability adds the missing population view: define the failure mode, describe uncertainty in the inputs, propagate those inputs through the performance function, and interpret the resulting probability together with consequence. Robustness is related but different. A robust design changes little when inputs vary. A reliable design stays on the safe side of a defined failure boundary with acceptable probability. Risk combines probability with consequence, so a rare failure with severe consequence can still require action.
The skills, taught in order
Nominal value versus population
A nominal value is a chosen representative number. A population contains the actual spread of strengths, dimensions, loads, temperatures, and use conditions.
Factor of safety versus reliability
The factor of safety compares selected values. Reliability integrates the overlap between load effects and resistance across their distributions.
Failure mode first
Reliability has no meaning until the failure event is defined: yield, fracture, excessive deflection, leakage, seizure, loss of preload, or missed function.
Robustness and risk
Robustness measures sensitivity to variation. Risk adds consequence to probability. Do not collapse these into one word.
Evidence boundary
A calculated failure probability is credible only inside the data, assumptions, and model checks that support it.
Engineering connection: use these skills to turn variability into a design decision, not just a plot.
Worked example 1: same nominal safety factor, different reliability
Two steel rods carry a nominal axial load of 40 kN and both have nominal yield resistance of 80 kN, so each has n = 2. Design A has resistance standard deviation 4 kN and load standard deviation 3 kN. Design B has resistance standard deviation 12 kN and load standard deviation 9 kN. Estimate failure probability assuming independent normal resistance and load.
- ProblemTwo steel rods carry a nominal axial load of 40 kN and both have nominal yield resistance of 80 kN, so each has n = 2. Design A has resistance standard deviation 4 kN and load standard deviation 3 kN. Design B has resistance standard deviation 12 kN and load standard deviation 9 kN. Estimate failure probability assuming independent normal resistance and load.
- AssumptionsResistance R and load L are independent normal variables. Failure is R - L <= 0. This is a first screen, not a certification calculation.
- ModelM = R - Lmu_M = mu_R - mu_Lsigma_M = sqrt(sigma_R^2 + sigma_L^2)P_f = Phi(-mu_M / sigma_M)
- SolveBoth designs have mu_M = 80 - 40 = 40 kN. Design A has sigma_M = sqrt(4^2 + 3^2) = 5 kN, so beta = 40/5 = 8 and P_f is about 6.2e-16. Design B has sigma_M = sqrt(12^2 + 9^2) = 15 kN, so beta = 40/15 = 2.67 and P_f is about 0.0038.
- CheckThe deterministic safety factor is identical, but the margin distribution is much wider for Design B. Wider scatter makes the lower tail cross zero more often.
- ConclusionThe same nominal factor of safety can hide orders of magnitude difference in failure probability.
Worked example 2: turning a safety check into a performance function
A bracket is acceptable if maximum bending stress is below yield strength. Nominal yield is 250 MPa and nominal stress is 125 MPa. Define the deterministic factor, the performance function, and the failure event.
- ProblemA bracket is acceptable if maximum bending stress is below yield strength. Nominal yield is 250 MPa and nominal stress is 125 MPa. Define the deterministic factor, the performance function, and the failure event.
- AssumptionsYield failure is the relevant first mode. Local stress concentration and fatigue are outside this first model.
- Modeln = S_y / sigmag(X) = S_y(X) - sigma(X)failure: g(X) <= 0
- SolveThe nominal factor is n = 250/125 = 2. The probabilistic model does not use n as the random output. It uses g = S_y - sigma, with failure when the random stress exceeds the random yield strength.
- CheckIf both S_y and sigma were fixed at nominal values, g = 125 MPa and failure would be impossible in the model. That is exactly why uncertainty must be represented.
- ConclusionThe deterministic check becomes a limit-state model once inputs are allowed to vary.
Python activity
Engineering question: Compare two designs with equal nominal safety factor but different scatter.
Assumptions and units: Use the units shown in the code comments. Keep the random seed fixed while learning so the result is reproducible.
from math import erf, sqrt
def normal_cdf(x):
return 0.5 * (1.0 + erf(x / sqrt(2.0)))
cases = {
"A": {"mu_r": 80.0, "sd_r": 4.0, "mu_l": 40.0, "sd_l": 3.0},
"B": {"mu_r": 80.0, "sd_r": 12.0, "mu_l": 40.0, "sd_l": 9.0},
}
for name, c in cases.items():
mu_m = c["mu_r"] - c["mu_l"]
sd_m = (c["sd_r"] ** 2 + c["sd_l"] ** 2) ** 0.5
beta = mu_m / sd_m
pf = normal_cdf(-beta)
print(name, "beta =", round(beta, 3), "Pf =", pf)
Misconceptions and diagnostics
| Mistake | Diagnostic question | Correction |
|---|---|---|
| A safety factor automatically tells me the reliability. | Ask: what scatter in load, strength, geometry, and environment did that single number include? | Correct the assumption, then rerun the calculation or rewrite the report. |
| Reliability is one minus the safety factor. | Ask: is the factor a probability? It is not. Reliability is one minus failure probability. | Correct the assumption, then rerun the calculation or rewrite the report. |
| A low calculated failure probability proves the design is safe. | Ask: does the model include the governing failure mode and credible input distributions? | Correct the assumption, then rerun the calculation or rewrite the report. |
| Uncertainty can be removed simply by collecting more samples. | Ask: is the uncertainty aleatory variation, epistemic model uncertainty, or both? | Correct the assumption, then rerun the calculation or rewrite the report. |
Practice ladder
State the difference between a nominal load and a load distribution.
Solution guidance
Write the governing equation, state units and assumptions, then compare the answer with the relevant limit state or design decision.
For R = 120 kN and L = 60 kN, compute the nominal factor of safety.
Solution guidance
Write the governing equation, state units and assumptions, then compare the answer with the relevant limit state or design decision.
A report says n = 1.8, therefore reliability = 0.44. Identify the error.
Solution guidance
Write the governing equation, state units and assumptions, then compare the answer with the relevant limit state or design decision.
Choose one mechanical component and define one failure event that could be analyzed probabilistically.
Solution guidance
Write the governing equation, state units and assumptions, then compare the answer with the relevant limit state or design decision.
Working with AI, and proving it yourself
Useful AI support
Portfolio evidence
Create one short reliability note for this module. It must include the engineering question, assumptions, input definitions, equations, calculation, checks, interpretation, and limitations.
Retrieval and spaced review
Closed notes. Answer out loud, then reveal.
What is the performance function used in this course?
g(X) = resistance(X) - load_effect(X), with failure when g(X) <= 0.
What does a factor of safety compare?
Selected deterministic values of resistance and load effect.
How are reliability and risk different?
Reliability is probability of satisfactory performance; risk combines probability and consequence.
Can two designs with the same safety factor have different failure probabilities?
Yes, if their input scatter or distribution shapes differ.
What should be defined before any reliability calculation?
The failure mode and the consequence context.
References and source discipline
This course synthesizes original MechCompass explanations from established reliability engineering sources. Confirm current editions and standards before using them for formal design approval.
| Topic in this module | Reference direction |
|---|---|
| Reliability framing and performance functions | Haldar and Mahadevan, Probability, Reliability, and Statistical Methods in Engineering Design |
| Engineering statistics and distribution checks | NIST/SEMATECH Engineering Statistics Handbook |
| Risk and design reliability context | Der Kiureghian, Structural and System Reliability |