Pseudo Mass-Spring-Dampers

• Pseudo mass-spring-dampers are reduced-order models that represent complex vibration and resonator systems using effective mass, damping, and stiffness parameters derived via asymptotic reduction or system identification.
• They capture nonlinear, amplitude-dependent effects such as inertia amplification and effective negative stiffness, enabling accurate reproduction of frequency response, damping, and transmissibility.
• These models support practical applications in metamaterials, sea state estimation, and vibration isolation through state-space embedding, experimental validation, and design optimization.

Pseudo mass-spring-dampers are reduced-order dynamical models in which complex vibration, wave-interaction, or resonator systems are represented by effective or “pseudo” mass ($M$), stiffness ($K$), and damping ($C$) parameters. These parameters are typically constructed via asymptotic reduction, energy-equivalence, or system identification such that the resulting ordinary differential equation or state-space system reproduces a subset of the critical dynamic properties—typically frequency response, damping, or transmissibility—of the underlying multi-degree-of-freedom or nonlinear system. In modern research, pseudo mass-spring-damper models unify applications spanning locally resonant acoustic metamaterials, nonlinear damped resonators with inertia amplification, advanced vibration isolators such as KDamper implementations, and data-driven estimators for inverse problems in seakeeping. This article gives a technical, method-centric summary, including derivations, design guidance, and validation results from key works (Damme et al., 2021, Tiwari et al., 27 Nov 2025, Sapountzakis et al., 2017).

1. Mathematical Formulation of Pseudo Mass-Spring-Damper Models

Pseudo mass-spring-damper models generalize the canonical second-order ODE: $m \ddot{x} + c\dot{x} + kx = F(t)$ by allowing $m$, $c$, and $k$ to be explicit functions of the system’s state, excitation amplitude, or parametric configuration, or by deriving them as effective quantities emerging from nonlinear, nonlocal, or coupled subsystems.

1.1. Nonlinear Amplitude-Dependent Form (Inertia-Amplification Resonators)

In locally resonant structures exploiting inertia amplification (IA), the governing equation adopts the form (Damme et al., 2021): $$F(t) = \tilde{M}(U) \, \ddot{U} + \tilde{C}(U) \, \dot{U}|\dot{U}| + \tilde{K} U$$ where $\tilde{M}(U)$ and $\tilde{C}(U)$ encode amplitude-dependent mass and quadratic damping: $$\tilde{M}(\epsilon) = M + m\left[1 + \frac{(1-\epsilon)^2}{A(\epsilon)}\right], \quad \tilde{C}(\epsilon) = \frac{m}{2L}\left[\frac{1-\epsilon}{A(\epsilon)} + \frac{(1-\epsilon)^3}{A(\epsilon)^2}\right]$$ with $\epsilon = U/(2L)$.

1.2. Energy-Equivalence in Wave-Excited Systems

For ship seakeeping, a pseudo mass-spring-damper model is constructed by normalization and parameter identification to match the vessel’s heave/pitch response amplitude operator (RAO) (Tiwari et al., 27 Nov 2025): $M \ddot{x}(t) + C \dot{x}(t) + x(t) = p(t)$ with $M(\eta)$, $C(\eta)$ explicitly mapped from vessel geometric and hydrodynamic parameters.

1.3. Effective Stiffness in Negative-Stiffness Isolators

In KDamper systems, an auxiliary negative-stiffness element is realized using preloaded Belleville disc springs, leading to a reduced form (Sapountzakis et al., 2017): $$m \ddot{x} + c_{\rm eff}(\dot{x}) + k_{\rm eff}(x) x = F(t)$$ where $k_{\rm eff}(x)$ includes the incremental negative stiffness derived as the derivative of a nonlinear force-displacement law.

2. Nonlinear and Amplitude-Dependent Effects

The pseudo mass-spring-damper paradigm is designed to capture effects—especially nonlinearity and parameter variability—not possible with linear models.

• Inertia-Amplification Resonators: Nonlinear kinematic couplings yield mass and damping coefficients strongly dependent on vibration amplitude ($A$), producing amplitude-dependent damping ($c_{\rm eff}(A)$) and frequency shift [$2106.02576$].
• KDamper Systems: Effective negative stiffness softens or hardens with excursion $u(x)$, imparting strong nonlinearity in transmissibility and isolation characteristics [$1705.05622$].
• Sea State Estimators: Pseudo parameters can be recursively updated to incorporate time-varying environmental forcing, parametric drift, and input stochasticity [$2511.21997$].

3. Analytical Reduction and Parameter Identification

Pseudo mass-spring-damper constructs originate via formal analytical reduction:

• Lagrangian Derivation: For IA resonators, full nonlinear kinematics are reduced using generalized coordinates, yielding closed-form expressions for $\tilde{M}(\epsilon)$, $\tilde{C}(\epsilon)$, and validating the amplitude-dependent behavior through energy analysis [$2106.02576$].
• Energy/RAO Matching: In wave-structure problems, equivalence is reached by matching the frequency response (RAO) using Nelson’s hydrodynamic formulas, mapping physical vessel parameters to pseudo $M(\eta)$, $C(\eta)$ [$2511.21997$].
• Effective Stiffness Linearization: In KDamper, the two-mass equations of motion collapse to a reduced pseudo system by eliminating the auxiliary DoF, with negative stiffness computed directly by differentiating the nonlinear Belleville spring law [$1705.05622$].

4. State-Space Modeling, Estimation, and Control

Pseudo mass-spring-damper models are well-suited for state-space embedding and model-based estimation/control.

• Sea State Estimation: The vessel-wave system is cast as a discrete-time state-space model, augmenting unknown parameters and excitation into the state, and employing a square-root Cubature Kalman Filter (SRCKF) for Bayesian estimation. Process noise covariance is derived analytically, and estimator uncertainty is rigorously lower-bounded using the posterior Cramér–Rao lower bound (PCRLB) [$2511.21997$].
• Dampers in Structural Control: The amplitude-dependent damping inherent in IA resonators allows for passive self-regulation of vibration attenuation, as demonstrated in experiments where the −3 dB bandwidth and phase slope of frequency response functions exhibit strong amplitude dependence [$2106.02576$].
• Design Optimization: For KDamper implementations, classical TMD tuning (Den Hartog) is extended to include negative stiffness and mass ratio, so double-peak equalization and broad-band minima are achieved at minimal auxiliary mass [$1705.05622$].

5. Experimental Validation and Prototyping

Extensive experimental work confirms the efficacy and predictive accuracy of pseudo mass-spring-damper models.

System	Model Form	Experimental Finding
IA resonators (Damme et al., 2021)	Nonlinear, amplitude-dependent $\tilde{M}(U),\tilde{C}(U)$	Amplitude-dependent peak broadening and frequency shift
Sea-state estimator (Tiwari et al., 27 Nov 2025)	Linear pseudo MSD with parameter adaptation	$<$2.5% error in $H_s$, $<$4% in $T_z$ in simulation
KDamper (Sapountzakis et al., 2017)	Effective negative stiffness via disc springs	Transmissibility peaks driven to near zero at low mass ratio

• Metamaterials and Resonators: Apparent dynamic mass amplification and amplitude-dependent loss are observed experimentally, confirming analytical predictions for IA-based pseudo dampers [$2106.02576$].
• Sea State Estimation: Monte Carlo simulation and high-fidelity panel-method data validate the ability of pseudo MSD-based estimators to recover significant wave height ($H_s$) and zero-crossing period ($T_z$) without prior knowledge of full transfer functions [$2511.21997$].
• KDamper Prototyping: Belleville-based KDamper reduces main mass vibrations by an order of magnitude with negligible tuning mass, confirmed across both displacement- and velocity-perturbed cases [$1705.05622$].

6. Design Guidelines and Parameter Selection

Design of pseudo mass-spring-damper systems proceeds by specifying target performance, then mapping specifications to structural and pseudo-parameter choices.

• Inertia-Amplification (IA):
  • Choose target resonance $f_0$ and “base” mass $M$.
  • Set auxiliary mass $m$ and geometric IA factor $\alpha=1/\sin^2\theta_0$.
  • Compute stiffness $k = (2\pi f_0)^2(M+m/\sin^2\theta_0)$.
  • Linearized damping $$c_0 = \frac{m}{2L}\frac{\cos^2\theta_0}{\sin^4\theta_0}$$; amplitude-dependent $c_{\rm eff}(A)=\frac{8 c_0 \omega_0 A}{\pi}$.
  • Keep strain-parameter $\epsilon=A/(2L) \lesssim 0.05$ for linear approximation fidelity [$2106.02576$].
• Sea-State Model:
  • Parameters $M(\eta)$ and $C(\eta)$ are estimated recursively, enabling model adaptation to real-time operational conditions, without fixed transfer functions [$2511.21997$].
• KDamper:
  • Select low mass ratio $\mu = m_p/m \sim 0.01$–$0.05$.
  • Negative stiffness ratio $\kappa = |k_n|/k_s$ tuned close to static stability threshold $\kappa_{\max}(\mu)$.
  • Dashpot $c_p$ is computed from “equal-peak” criterion analogous to TMD optimization [$1705.05622$].

7. Applications and Performance Boundaries

Pseudo mass-spring-damper models support a diverse range of applications:

• Metamaterial and phononic unit cells: Apparent dynamic mass and tunable bandgaps at low frequency via IA mechanisms [$2106.02576$].
• Onboard sea state measurement: Real-time estimation with limited sensor information and no a priori vessel transfer functions [$2511.21997$].
• Ultra-low-frequency vibration isolation: KDamper implementations achieving broadband isolation with minimal auxiliary mass and negative-stiffness elements [$1705.05622$].

Limitations include stability concerns for negative-stiffness elements (critical $\kappa$), model validity for large amplitudes or non-linearity, and the restriction of sea-state estimation to heave/pitch with unidirectional long-crested waves. In all cases, experimental and simulation results converge closely with theory, supporting wide adoption of the pseudo mass-spring-damper modeling framework.