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 Kasif Furkan Ozturk, Tufan Cakir and Onur Araz of Gümüşhane University examined how the frequency content of ground motions and the thickness of such an isolation layer jointly govern the dynamic response of a five-story building. Their parametric study, published in the Bulletin of Earthquake Engineering, combined analytical modeling with a carefully classified set of earthquake records, comparing the structure resting on natural soil against the same structure founded on a geotechnical seismic isolation, or GSI, layer of six different thicknesses. The results carry practical weight for designers considering this technology.

The appeal of GSI systems lies in their simplicity of concept. Rather than installing mechanical base isolators, such as lead rubber bearings or friction pendulum devices, between a structure and its foundation, engineers instead surround or underlie the foundation with a engineered soil mixture, often soil blended with granulated scrap tire rubber. This softer, more damping-rich medium partially absorbs and dissipates seismic energy before it can climb into the building, filtering the shaking at the ground level itself. Because the approach recycles petroleum-derived waste materials such as discarded tires, it also addresses a genuine gap in waste recovery, turning an environmental liability into structural protection.

To test how well this filtering survives the messy variety of real earthquakes, the researchers grouped their ground motions according to frequency content and ran dynamic analyses across the full set. The findings on this front were consistent and, in some respects, sobering. For structures founded both on natural soil and on GSI layers, earthquakes with low frequency content produced structural responses greater than the average responses calculated across all the investigated records. In other words, slow, long-period ground motions tended to be the harshest critics of the building, regardless of what was buried beneath its foundations, and designers who evaluate structures against a generic suite of motions without regard to frequency content risk underestimating the demands that the most damaging classes of events can impose.

This frequency dependence echoes a growing body of research on soil-structure interaction. Previous work by members of the same group has shown that the frequency character of ground motion shapes the seismic response of retaining walls, box culverts, and buildings of varying aspect ratio when interaction with the surrounding soil is included. The new study extends that lesson into the GSI domain, confirming that an isolation layer does not erase the influence of how the shaking oscillates through time; the layer reshapes and attenuates the transmitted motion, but the character of the input still leaves its fingerprint on the building above.

The thickness of the isolation layer emerged as the other central variable, and here the study delivers its most actionable result. Increasing GSI thickness significantly attenuates structural responses, as expected, since a thicker compliant layer forces more of the seismic energy to deform and dissipate within the mixture rather than propagate upward. But the relationship is not linear, and the benefits diminish decisively with depth. Beyond a thickness of three meters, the researchers found that adding more material does not lead to a relatively significant further reduction in the dynamic responses of the investigated structural system.

That finding translates directly into a design recommendation. For the structural system, the natural soil conditions, the rubber-soil mixtures, and the earthquakes examined in the study, the authors suggest that selecting a GSI thickness in the range of two to three meters, inclusive, may be regarded as a practical approach for mitigating dynamic responses. Thicknesses beyond that range provide only progressively limited improvements, meaning the extra excavation, material, and cost buy shrinking returns. The proposal offers engineers a defensible starting point for a technology that has historically lacked the codified design procedures available for conventional base isolation.

The broader context makes the study timely. Geotechnical seismic isolation has advanced rapidly over the past decade, propelled by laboratory characterizations of sand-rubber and gravel-rubber mixtures, centrifuge tests demonstrating their protective performance, and shaking table experiments on medium-rise buildings founded on EPS bead-sand mixtures. Researchers have explored related variants, including soil mixed with recycled tire rubber, bitumen, and high-damping polyurethane, and have applied the concept to rural residences, highway embankments, and cut-and-cover utility tunnels. Analytical design models and fragility frameworks have followed, gradually building the engineering infrastructure needed to move GSI from promising research to routine practice.

What distinguishes the Turkish team’s contribution is its systematic treatment of two factors that interact: the spectral character of the earthquake and the geometry of the isolation layer. By coupling a frequency-classified ground motion set with a six-point thickness parametric study, the researchers generated a response surface that reveals where protection is most effective and where it saturates. Their analyses also confirm that the GSI system’s benefit is best understood as partial energy dissipation before the waves reach the foundation, a framing that positions the technology as a complement to, rather than a wholesale replacement for, structural isolation systems in high-consequence applications.

The authors note that their results should support the development of vibration control strategies aimed at improving the seismic resilience of structures, and that the data are expected to be useful for a better understanding of GSI applications in general. With no external funding and no competing interests declared, the study stands as a focused analytical contribution from a group with a sustained record in soil-structure interaction research. As cities in seismically active regions search for affordable, sustainable ways to protect ordinary buildings, the message from Gümüşhane is refreshingly concrete: bury a modest, two-to-three-meter cushion of rubber-enhanced soil beneath the foundation, recycle waste tires in the process, and expect the greatest gains in shielding structures from the low-frequency shaking that matters most, while resisting the temptation to assume that more isolation layer always means more protection.

Subject of Research: Seismic performance of buildings resting on geotechnical seismic isolation layers of varying thickness under earthquakes with different frequency content

Article Title: Effects of earthquake frequency content and isolation thickness on seismic behavior of a structure system resting on geotechnical seismic isolation layer

Article References: Ozturk, K. F., Cakir, T., & Araz, O. (2026). Effects of earthquake frequency content and isolation thickness on seismic behavior of a structure system resting on geotechnical seismic isolation layer. Bulletin of Earthquake Engineering. https://doi.org/10.1007/s10518-026-02658-8

Image Credits: AI Generated

DOI: 10.1007/s10518-026-02658-8

Keywords: geotechnical seismic isolation, earthquake engineering, rubber-soil mixtures, soil-structure interaction, frequency content, isolation thickness, seismic isolation, waste tire recycling, structural dynamics, vibration control, Bulletin of Earthquake Engineering, Gümüşhane University