The core challenge the researchers set out to solve is uncertainty. In an idealized computer model, bearing stiffness and damping are tidy numbers pulled from a manufacturer’s catalog. In the real world, they are anything but. Manufacturing tolerances, assembly errors, thermal expansion mismatches between the steel shaft and its housing, and speed-dependent centrifugal forces all conspire to shift these parameters by roughly 5 to 12 percent from their nominal values. Traditional deterministic designs may look excellent on paper but can fluctuate wildly or fail outright when reality deviates from the assumed values. Worse, most prior studies have relied on theoretical catalog data rather than measurements, ignoring what engineers call epistemic uncertainty, the gap between the model and the physical machine.

Wu and Ma’s first move was to measure rather than assume. They built a high-speed motorized spindle test rig spanning an operational range from 1,000 to 24,000 rpm, discretized into 50 speed points. At each point, the spindle ran under steady-state thermal conditions while high-precision eddy current sensors sampled vibration at 10 kilohertz. Repeating the identification procedure 20 times per speed produced a dataset of 1,000 identified parameter vectors. To extract bearing stiffness and damping coefficients from these measured frequency responses, the team employed the Multi-Innovation Stochastic Gradient, or MISG, algorithm, a recursive estimation technique that accelerates convergence by stacking innovation information from several consecutive time steps rather than using a single sample at a time.

The numerical verification of this algorithm is striking in its own right. In synthetic tests where the true parameter values were known, MISG reduced the number of convergence iterations from 606 to just 75, an 8.1-fold speedup and an 87.6 percent reduction compared with the standard stochastic gradient method, while also lowering the final normalized root-mean-square error by 15 percent. Although each iteration costs about 18 percent more computation, the faster convergence cuts total computation time by 28 percent. A sensitivity study further confirmed that algorithmic identification error contributes less than 3 percent of the total output variance, meaning the scatter in the identified parameters reflects genuine physical variability, not artifacts of the estimation procedure.

With the empirical distributions in hand, the team validated them statistically before feeding them into the optimizer. The identified stiffness parameters showed coefficients of variation between roughly 6.0 and 6.4 percent, and Shapiro-Wilk tests confirmed that a multivariate normal distribution was a physically permissible model, with quantile-quantile plots yielding coefficients of determination above 0.99. Bearing damping was not treated as an independent random variable; instead, it was coupled proportionally to the identified stiffness through a constant calibrated by modal hammer testing, keeping the probabilistic model lean while preserving fidelity. The identified stiffness values deviated only about 1.5 to 3.0 percent from the numerical baseline, with the small residual attributed to installation-specific effects such as thermal expansion mismatch and bearing-seat tolerances, exactly the factors catalog values ignore.

The optimization phase then framed spindle design as a bi-objective problem: minimize both the mean peak vibration amplitude at the tool tip and its standard deviation, subject to a probabilistic constraint that the probability of vibration exceeding a critical threshold of 195 micrometers stays below a target level. That threshold was calibrated from surface roughness requirements in precision milling, where larger displacements degrade finish beyond aerospace-grade limits. Rather than relying on first-order or second-order reliability methods, which linearize the problem and can incur errors of 20 to 30 percent in nonlinear bearing-rotor dynamics, the researchers propagated uncertainty through direct Monte Carlo simulation with 1,000 samples per candidate design. The Non-dominated Sorting Genetic Algorithm II then evolved populations of designs, with each run producing a Pareto front of optimal trade-offs.

Crucially, the team ran the optimization three times, for target failure probabilities of 1, 5, and 10 percent, generating a family of Pareto fronts that reveals the marginal cost of safety. The results expose sharp diminishing returns. Tightening the constraint from 10 to 5 percent delivered a high marginal benefit ratio of 0.55 micrometers of variance reduction per percentage point of reliability, but squeezing from 5 to 1 percent dropped that ratio to 0.39 and demanded a 7.0 percent increase in bearing preload, pushing it to its physical limit of 3,000 newtons, for only a marginal robustness gain. Such aggressive preload accelerates rolling-element wear and risks friction-induced thermal problems. Weighing these trade-offs against ISO 20816-3 vibration standards, the researchers identified the 5 percent configuration as the optimal engineering equilibrium.

The optimized design itself is physically revealing. The framework selected a front-bearing preload of about 2,804 newtons and shortened the bearing support span from roughly 750 to 537 millimeters, a geometric change that raises the first natural frequency and pushes the third critical speed from about 22,000 rpm to 25,800 rpm, safely beyond the 24,000 rpm maximum operating speed. Benchmarking against conventional approaches underscored the value of the joint formulation: a pure robust design optimization achieved low variability but left the failure probability at an unacceptable 21.07 percent, while a pure reliability-based design hit a low failure probability but suffered an inflated response standard deviation of 11.80 micrometers. The integrated framework achieved failure probability of 0.35 percent while keeping variability at 9.69 micrometers and mean amplitude at 16.44 micrometers.

Physical validation sealed the case. Twenty independent replication runs at the rated speed of 24,000 rpm showed the optimized spindle reducing mean peak vibration by 16.84 percent, from 19.77 to 16.44 micrometers, and the response standard deviation by 29.78 percent, from 13.80 to 9.69 micrometers, with all improvements statistically significant in paired-sample t-tests. Across the full 1,000 to 24,000 rpm spectrum, the optimized spindle maintained stable, constrained responses, eliminating the severe resonant amplification that plagued the baseline design near its third critical speed. The variance reduction of 50.70 percent in experimental terms also surpasses the roughly 22 percent typically reported in simulation-based robust optimization studies.

Beyond the immediate numbers, the study signals a broader shift in how rotating machinery can be engineered. By grounding probabilistic models in measured reality rather than catalog assumptions, and by mapping the full trade-off surface between performance, robustness, and reliability instead of fixing a single safety target, the framework gives engineers a quantitative language for deciding how much safety margin a given application truly needs. The authors note that the approach currently assumes linear bearing behavior within the Hertzian elastic regime and depends on sufficient experimental data, pointing toward amplitude-dependent identification schemes and non-probabilistic interval methods as future extensions for heavily loaded or multi-shaft systems. For now, the message is clear: the path to quieter, more dependable high-speed spindles runs directly through the test rig, with every ounce of uncertainty measured, quantified, and designed against.

Subject of Research: Reliability-based robust optimization of high-speed motorized spindle-bearing systems using experimentally identified bearing parameter uncertainties

Article Title: Experimental identification and robust optimization of spindle–bearing systems with reliability constraints

Article References: Wu, H., & Ma, J. (2026). Experimental identification and robust optimization of spindle–bearing systems with reliability constraints. Mechanical Sciences, 17(2), 839-855. https://doi.org/10.5194/ms-17-839-2026

Image Credits: AI Generated

DOI: 10.5194/ms-17-839-2026

Keywords: motorized spindle, bearing stiffness, vibration control, robust optimization, reliability-based design, uncertainty quantification, Monte Carlo simulation, NSGA-II, multi-innovation stochastic gradient, CNC machining, Pareto optimization, rotor dynamics

Researchers Van Thuy Tran, Thi Thanh Truc Nguyen, and Huu Loc Nguyen, affiliated with Ho Chi Minh City University of Technology and Pham Van Dong University in Vietnam, developed a reliability-based framework that combines Monte Carlo simulation with Response Surface Methodology to determine exactly how much adjustment a V-belt drive actually needs. Their central insight is deceptively simple: the geometric parameters that govern a belt drive, including the belt length and the pitch diameters of the driver and driven pulleys, are never exactly what the drawings say they are. Manufacturing tolerances, assembly errors, and material deformation introduce small random deviations, and those deviations propagate through the drive geometry in ways that deterministic formulas cannot capture.

The traditional approach divides the total adjustment of the center distance into two components. The first, often called take-up, is the extra travel needed to stretch the belt to its initial tension and to re-tension it as it elongates over time through creep, wear, and permanent deformation. The second, called slack-off, is the allowance that lets the belt slip over the pulleys during installation despite dimensional deviations. Design manuals typically prescribe fixed values for both, derived from empirical rules rather than probabilistic analysis. The Vietnamese team argued that these fixed values are blind to the statistical reality of the components, and they set out to replace them with intervals defined at a specified confidence level of 99.9 percent.

To do this, the researchers first modeled the three key geometric parameters as random variables following normal distributions. The choice of the Gaussian model rests on the central limit theorem: dimensional variations arise from many independent sources, such as machining errors, measurement uncertainty, and assembly variations, so their combined effect tends toward a normal distribution. Working with standard dimensions for classical Type-B V-belts drawn from ISO 4183, ISO 4184, and JIS K 6323, the team adopted baseline uncertainty levels of plus or minus 0.3 percent for pulley diameters and plus or minus 0.5 percent for belt length. These were treated as modelling assumptions rather than prescribed tolerance limits, and each range was interpreted as roughly plus or minus three standard deviations.

Monte Carlo simulation then became the engine of the analysis. In each of millions of iterations, the method draws a random set of values for belt length and both pulley diameters, computes the resulting center distance from the standard geometric relationship, and records the outcome. Repeated thousands or millions of times, this procedure builds up a full probability distribution of the center distance and of the adjustment range, revealing not just an average value but the extreme quantiles that matter for design. However, running the full geometric model millions of times is computationally expensive, which is where Response Surface Methodology enters the picture.

Using a Box-Behnken experimental design with fifteen runs, the team constructed a second-order regression model that approximates the nonlinear relationship between the input parameters and the adjustment distance. The fitted model proved remarkably accurate, explaining 97.45 percent of the variance in the data, with adjusted and predicted coefficients of determination both at 95.6 percent. Variance inflation factors of exactly 1.0 confirmed that the input variables were free of multicollinearity. The surrogate model then allowed Monte Carlo sampling to run at a fraction of the original computational cost, making a reliability-based design calculation practical rather than prohibitive.

The sensitivity analysis that followed produced a clear hierarchy of influence. Belt length emerged as the dominant factor, followed by the driver pulley diameter, while the driven pulley diameter played a weaker but still non-negligible role. Local partial derivatives at the design center quantified this ranking precisely. More striking were the nonlinear effects: the driver pulley diameter showed a concave response peaking near 160 millimeters, and the driven pulley diameter exhibited an asymmetric, non-monotonic trend shaped by a strong quadratic term and its interaction with the driver pulley. Physically, these curvatures arise because changes in pulley diameters simultaneously alter the effective belt length, the wrap angle, and the contact conditions between belt and pulley, producing effects that no linear approximation can faithfully represent.

The numerical example that anchored the study involved a conveyor transmission system with a rated power of 6.02 kilowatts and an input speed of 986 revolutions per minute. A Monte Carlo simulation with five million realizations showed that the installation slack-off adjustment follows a tight Gaussian distribution centered at 22.01 millimeters with a standard deviation of only 0.037 millimeters. At the 99.9 percent confidence level, the required slack-off range spans just 21.90 to 22.12 millimeters. The take-up adjustment, governed by belt elasticity and initial tension, ranged far more widely, from 41.92 to 96.16 millimeters at the same confidence level. A sensitivity check using truncated normal and uniform distributions changed the results by less than 0.13 percent, confirming that the conclusions do not hinge on the assumed distribution type.

The comparison with conventional practice is where the findings become genuinely provocative. The Bando V-belt Design Manual recommends a mechanical slack-off allowance of 31.75 millimeters for the studied configuration, while the conventional empirical calculation yields 35.58 millimeters. The reliability-based result of roughly 22 millimeters shows that about 30 percent of the traditionally prescribed installation margin is simply unused, representing wasted mechanism travel and unnecessary bulk. Meanwhile, the take-up story runs in the opposite direction: the probabilistic upper bound of 96.16 millimeters substantially exceeds the nominal value of about 61 millimeters, suggesting that traditional methods may underestimate how much tensioning travel is truly needed once variability is accounted for. This asymmetry, in which empirical rules simultaneously overestimate one adjustment component and underestimate the other, is a previously unreported insight and a fundamental critique of deterministic design.

To make the framework usable in practice, the researchers compiled lookup tables of upper design values for both adjustment components across standardized belt lengths from 1600 to 2800 millimeters, driver pulley diameters of 140 to 180 millimeters, and transmission ratios from 1.25 to 6.3. These tables convert millions of simulated realizations into numbers an engineer can read directly, without repeating the simulation. The authors acknowledge that experimental validation remains future work, and that the uncertainty levels were modelling assumptions rather than measured tolerances, but the framework is built to absorb application-specific data as they become available. The broader significance extends beyond belts: the study demonstrates how uncertainty can be transformed from a post-analysis afterthought into a direct design variable, offering a template for reliability-informed decisions across mechanical transmission systems where geometry, tolerance, and chance quietly interact.

Subject of Research: Reliability-based determination of adjustable center distance in V-belt drives using Monte Carlo simulation and response surface methodology

Article Title: Reliability-based determination of adjustable center distance in V-belt drives using Monte Carlo simulation and response surface methodology

Article References: Tran, V. T., Nguyen, T. T. T., & Nguyen, H. L. (2026). Reliability-based determination of adjustable center distance in V-belt drives using Monte Carlo simulation and response surface methodology. Heliyon, 12(15), Article e45554. https://doi.org/10.1016/j.heliyon.2026.e45554

Image Credits: AI Generated

DOI: 10.1016/j.heliyon.2026.e45554

Keywords: V-belt drives, Monte Carlo simulation, Response Surface Methodology, reliability-based design, center distance adjustment, uncertainty quantification, Box-Behnken design, sensitivity analysis, power transmission, mechanical design, probabilistic modeling, belt tension