4. Sources of uncertainty in airbag correlation
Airbag correlation is the process of aligning a CAE model with physical sled or barrier test results so that the simulation can be trusted for restraint-system design decisions. In practice, correlation is never a single deterministic comparison. Even when the same vehicle, the same airbag module, and the same crash pulse are used, repeated tests scatter around a mean response. This scatter is the uncertainty budget that any credible correlation study must quantify. Four dominant contributors are material variability, manufacturing tolerances, sensor drift, and environmental conditions. Understanding how each enters the model, how large its effect typically is, and how to separate it from genuine model-form error is the core skill of this chapter.
Material Variability
Airbag fabrics, coatings, and inflator propellant behave differently from lot to lot. The dominant fabric parameters are areal density, yarn tenacity, weave construction, and the silicone or polyurethane coating weight. For a typical 470 dtex nylon 6.6 fabric, areal density might be specified at 235 ± 10 g/m², and biaxial tensile stiffness can vary by 8–15 % between production lots. Because the bag's internal pressure and deployment kinematics depend on fabric permeability and stiffness, this variability propagates directly into the occupant response.
A practical way to bound the effect is a one-factor-at-a-time (OFAT) sensitivity sweep. Hold all other inputs at nominal, then run the model at the low and high material bounds:
# Example OFAT sweep for fabric areal density
nominal_density = 235 g/m2
low_density = 225 g/m2 # -4.3 %
high_density = 245 g/m2 # +4.3 %
# Observed effect on 50th percentile HIC36 (illustrative)
# low_density -> HIC36 = 612
# nominal -> HIC36 = 640
# high_density -> HIC36 = 671
# Sensitivity = (671 - 612) / (2 * 640) = 4.6 % per 10 g/m2
The numbers above are illustrative, not measured. The point is the method: convert each material parameter into a normalized sensitivity coefficient so that parameters with different units can be ranked. Inflator propellant burn rate is usually the single largest material contributor, often producing 5–10 % variation in bag peak pressure and 3–8 % in occupant HIC. Fabric permeability is second, and coating weight third.
Manufacturing Tolerances
Manufacturing tolerances are geometric and assembly-related rather than constitutive. They include bag panel cut dimensions, seam stitch density and thread tension, vent-hole diameter and location, inflator mounting clocking, and module-to-steering-wheel gap. A vent hole specified at Ø 25 ± 0.5 mm changes effective vent area by roughly ± 4 %, which shifts bag pressure decay and can move chest deflection by 2–5 mm.
Unlike material variability, manufacturing tolerances are often asymmetric and correlated. Stitch density and thread tension tend to vary together because they come from the same sewing machine setup. A correlation study that treats them as independent will underestimate the combined effect. The recommended practice is to build a tolerance stack-up table and then run a small Design of Experiments (DOE) rather than a full factorial.
| Parameter | Nominal | Tolerance | Typical effect on HIC36 |
|---|---|---|---|
| Vent hole diameter | 25 mm | ± 0.5 mm | ± 2 % |
| Seam stitch density | 45 stitches/100 mm | ± 3 | ± 1.5 % |
| Inflator clocking | 0° | ± 5° | ± 3 % |
| Module gap | 2.0 mm | ± 0.5 mm | ± 1 % |
Sensor Drift
Sensor drift is often overlooked because it is a measurement-system error, not a physical variability. Accelerometers, load cells, and pressure transducers all exhibit zero offset, sensitivity shift, and thermal drift. A 1 % sensitivity error on a chest accelerometer integrates twice into displacement and can produce a 5–10 mm error in chest deflection over a 100 ms event. This is comparable to the physical effects being studied, so it cannot be ignored.
Mitigation follows a calibration discipline. Before each test series, perform an end-to-end shunt calibration and record the zero-offset. During correlation, treat the measured signal as y_meas = y_true + b + s·y_true, where b is offset and s is sensitivity error. If the same sensor chain is used for all tests, the drift is systematic and will bias the correlation rather than add scatter. The practical fix is to include a calibration uncertainty term in the correlation metric, or to re-calibrate between test groups.
Environmental Conditions
Environmental conditions affect both the physical test and the material state. The most influential are ambient temperature and humidity at the time of deployment, and the soak temperature of the inflator and fabric. Propellant burn rate is strongly temperature-dependent: a typical gas generator may show a 15–25 % change in peak pressure between −35 °C and +85 °C. Fabric stiffness and permeability also shift with moisture content, though the effect is smaller, typically 3–7 % over a 10–90 % RH range.
For correlation, the test environment must be recorded and matched in the model. If the physical test was conducted at 23 °C and 50 % RH but the model uses room-temperature material cards, the correlation will absorb the environmental difference as apparent model error. The correct approach is to define a reference condition, then apply temperature and humidity scaling factors to the material inputs before comparison.
Building the Uncertainty Budget
The four sources above do not act independently in a real correlation study. A defensible workflow is:
- Classify each source as random (lot-to-lot material, tolerance stack) or systematic (sensor calibration, environmental mismatch).
- Quantify each source with a sensitivity coefficient from OFAT or DOE runs, expressed in the same output unit (e.g., % HIC36).
- Combine random sources in root-sum-square (RSS) if they are independent, or use a covariance term if they are correlated.
- Subtract systematic bias from the correlation error before judging model quality.
- Report the residual as the model-form uncertainty that remains after all known sources are accounted for.
For example, if material variability contributes ± 4 %, tolerance ± 3 %, and sensor drift ± 2 % (all random and independent), the RSS combined uncertainty is √(4² + 3² + 2²) = 5.4 %. If the observed correlation error is 7 %, then roughly 5.4 % is explained by input uncertainty and only about 1.6 % remains as genuine model error. Without this decomposition, the engineer might wrongly conclude the model is poor and over-tune it, degrading its predictive capability.
The practical takeaway is that airbag correlation is a statistical exercise, not a single overlay plot. By explicitly budgeting material variability, manufacturing tolerances, sensor drift, and environmental conditions, the engineer can state not just whether the model matches, but how much of the mismatch is attributable to known uncertainty and how much is real model-form error worth fixing.