--- title: Haptic Feedback Calibration url: https://www.emergentmind.com/topics/haptic-feedback-calibration type: topic --- # Haptic Feedback Calibration Searching arXiv for recent papers on haptic feedback calibration and related calibration methodologies. Haptic feedback calibration comprises the procedures that make sensed state, actuator output, coordinate systems, and perceived sensation consistent with an intended interaction. In the cited literature, calibration appears in encountered-type robotic displays, sim-to-real manipulation, friction-modulating touchscreens, wearable vibrotactile systems, electrotactile and EMS interfaces, and deformable virtual environments. Its concrete targets vary by platform: spatial co-location between virtual and physical contact, accurate transformation between sensor and robot frames, equalization of actuator output, compensation of mechanical and electrical transfer functions, stabilization of latency-sensitive control loops, and user-specific matching of perceived intensity across sites or modalities [2309.16768][2507.08572][2504.21477]. ## 1. Calibration objectives across haptic system classes A first class of calibration problems is geometric. In encountered-type displays, the central requirement is that a real actuator reach the same locus and time as a virtual contact. The Tracking Calibrated Robot (TCR) formulates this as the consistency condition that when virtual contact occurs, $\|x_f^v-x_o^v\|=0$, the real system should also satisfy $\|x_f^r-x_o^r\|=0$, and more generally that $\|x_f^v-x_o^v\|-\|x_f^r-x_o^r\|=0$ throughout the approach [2309.16768]. A related but distinct geometric problem appears in robotic sim-to-real transfer, where the simulator provides $x_s=f(q_s)\in\mathbb{R}^3$ while reality provides only touchscreen contact coordinates $x_r=(x,y)\in\mathbb{R}^2$ with $z_r=0$ at contact; calibration then becomes a learned transformation $C_S=T(C_R)$ that corrects the simulated target before inverse kinematics [2507.08572]. A second class is transduction calibration. Surface haptics require friction, latency, and spatial uniformity to be measured physically rather than inferred from commanded input. One framework defines the friction coefficient as $\mu=F_T/F_N$, the friction range as $\Delta\mu=\mu_H-\mu_L$, and end-to-end latency as $\Delta T=t_2-t_1$ [2308.15190]. HapTable similarly calibrates electromechanical vibrotactile actuation by measuring frequency response functions over an 84-point grid and calibrates electrostatic friction by voltage-driven normal-force modulation [2103.16510]. A third class is perceptual calibration. Wearable and electrical interfaces must account for strong inter-user variability. SensoPatch therefore separates normalization of 25 piezoresistive sensors from equalization of perceived vibration across body sites such as the upper arm, shoulder, and lower back [2409.19155]. The wearable haptic bracelet study shows that normal and shear cues cannot be equalized globally by either displacement or force, and instead require a user-specific point of subjective equality (PSE) [1911.02104]. The survey of electrical stimulation haptics generalizes this logic: thresholds, dynamic ranges, and comfortable operating regions depend on site, electrode geometry, waveform, age, and sex, so calibration must be per user and often per electrode [2504.21477]. ## 2. Spatial and geometric alignment In robotic and mixed-reality systems, calibration commonly begins with rigid registration. TCR estimates a transformation between the VR world frame and the robot base frame from paired samples of controller positions and robot end-effector positions. With centered point clouds $A\in\mathbb{R}^{3\times N}$ and $B\in\mathbb{R}^{3\times N}$, the paper solves the orthogonal Procrustes problem $$ R=\arg\min \|\Omega A-B\|_F \quad \text{subject to } \Omega^T\Omega=I, $$ using the SVD of $M=BA^T$, $M\to U\Sigma V^T$, and $R=UV^T$. Translation is then computed as $t=\mu_B-R\mu_A$, yielding $$ T_{BV}=\begin{bmatrix}R & t\\0 & 1\end{bmatrix}. $$ Using 486 paired samples, the reported mean squared alignment error is $0.00485$ m, ანუ $4.85$ mm [2309.16768]. The same structural idea appears in planar contact tasks, but with a different geometry. For touchscreen-based sim-to-real transfer, three reachable edge points on the physical screen define a best-fit screen plane in simulation. With points $P_0$, $P_x$, and $P_y$, the basis is built as $e_x=(P_x-P_0)/\|P_x-P_0\|$, $e_y=\mathrm{proj}_{plane}(P_y-P_0)/\|\mathrm{proj}_{plane}(P_y-P_0)\|$, and $e_z=e_x\times e_y$, giving $$ R_S^B=[e_x\ e_y\ e_z],\qquad t_S^B=P_0. $$ Any screen point then maps by $[X_B\ Y_B\ Z_B]^T=R_S^B[X_S\ Y_S\ 0]^T+t_S^B$ [2507.08572]. This suggests that planar haptic calibration often reduces a 3D problem to a lower-dimensional contact manifold plus a residual correction model. Tool-based systems add frame transformations at the wrench level. In robotic surgery training, the wrist sensor frame $S$, tool-tip frame $T$, robot base $B$, and haptic-device frame $H$ are linked by adjoint mappings, with force-only rendering implemented by successive rotations, $$ {}^Tf=R_{T_S}\,{}^Sf_{ext},\qquad {}^Bf=R_{B_T}\,{}^Tf,\qquad {}^Hf=R_{H_B}\,{}^Bf. $$ The TCP at the tip is calibrated explicitly because every subsequent compensation and rendering step depends on it [2604.27385]. SmartBelt shows a geometry-free alternative. Instead of measuring microphone locations around the waist, it interpolates a 360-by-28 lookup table of expected TDoAs from eight calibration recordings and then derives the personalized motor directions from TDoA zero-crossings of adjacent microphone pairs. This yields a per-user set of motor angles $\theta_i$ that adapts to waist size and belt deformation without explicit geometric measurement [2202.13974]. ## 3. Sensor, actuator, and transfer-function calibration Once geometry is fixed, calibration shifts to the physical relation between command and delivered haptic output. In teleoperated construction, accelerometer outputs are calibrated by bias, scale, and axis-alignment correction, $$ \mathbf{a}_{tool}(t)=\mathbf{R}\,\mathbf{S}\,\big(\mathbf{a}_{meas}(t)-\mathbf{b}\big), $$ and actuator placement is optimized using the handle-to-source energy ratio $$ E_{\mathrm{ratio}}=\frac{E_x+E_y+E_z}{E_x^*+E_y^*+E_z^*}. $$ The paper reports that the lower sensor location on the tool gave the highest SNR for contact and that the actuator orientation perpendicular to the upper part of the handle gave the highest $E_{\mathrm{ratio}}$ [2302.00741]. Wearable string-based devices use a simpler but explicit actuation model. The lightweight fingernail device converts motor torque to string tension through $\tau_m=K_t i$ and $T\approx\tau_m/r_{eff}$, then uses an empirically fitted linear current–force relation in software. Bench calibration found a linear current–force relation up to approximately $1.36$ N, with the shortfall from the theoretical $1.57$ N attributed to internal friction and the increased effective winding radius due to string diameter [2506.21417]. Surface haptics often require full frequency-domain characterization. HapTable measures displacement FRFs over a $7\times12$ grid under a $0$ to $625$ Hz sine sweep and represents the local surface response as $$ H(\omega,r)=\frac{X(\omega,r)}{U(\omega)}. $$ For multi-actuator synthesis, the linear model is $$ X(\omega,r)\approx \sum_{i=1}^{N} H_i(\omega,r)U_i(\omega), $$ and, in a general multichannel implementation, actuator inputs would be obtained from a regularized least-squares solution. HapTable instead uses precomputed lookup tables that select actuator-frequency combinations producing maximal displacement contrast between source and destination loci [2103.16510]. Electrovibration calibration is more explicitly model-based. One recent touchscreen study models electrovibration-induced friction as a first-order low-pass mapping from message voltage to friction, $$ H_{f-e}(s)=\frac{F_f(s)}{V_m(s)}=\frac{K}{1+\tau s}, $$ with a speed-dependent cutoff $$ f_c(\nu)=385.68+13.811\nu \quad \text{Hz}, $$ valid for $\nu$ in mm/s, and average gain $K\approx0.0123$ N/V at 150 Vpp [2505.11162]. The input waveform is amplitude-modulated to undo the $V^2$ nonlinearity, $$ V_i(t)=\sqrt{V_m(t)+|\min(V_m(t))|}\cos(2\omega_c t), $$ with a carrier of 7 kHz [2505.11162]. In force-sensing surgical systems, calibration begins from the factory sensor model $w_{meas}=Cv+b$, followed by bias capture, gravity compensation, and moving-average smoothing of the compensated force estimate [2604.27385]. The common pattern across these systems is that calibration is rarely a single step: it couples static parameter identification, placement optimization, and dynamic filtering. ## 4. Psychophysical equalization and user-specific mapping Perceptual calibration is required when the same physical command does not yield the same felt intensity across users, sites, or modalities. The bracelet study provides a direct example. Using a staircase-based method of adjustments with a two-up one-down rule, it found that for normal references of 1 mm, 2 mm, and 3 mm, the corresponding shear PSEs were 1.7 mm, 2.7 mm, and 3.7 mm, while a 3 mm shear reference corresponded to a normal PSE of 1.9 mm [1911.02104]. The group-average displacement mapping is therefore summarized as $$ d_{shear}\approx d_{normal}+0.7\ \text{mm}, $$ but the paper explicitly reports substantial inter-user variability and concludes that a single global mapping is insufficient [1911.02104]. SensoPatch extends this logic to modular vibrotactile feedback. The glove side uses region averaging and thresholding, while the actuator side may be equalized psychophysically by fitting a logistic psychometric function $$ p(x)=\gamma+(1-\gamma-\lambda)\frac{1}{1+\exp(-(x-x_0)/s)} $$ to perceived intensity or detection responses for each motor and body site [2409.19155]. The paper’s experiments found intensity discrimination accuracies of 76% on the upper arm, 76% on the shoulder, and 73% on the lower back; single-motor location discrimination was 90% on the upper arm, 88% on the shoulder, and 81% on the lower back; two-motor simultaneous discrimination was best on the shoulder at 66% [2409.19155]. This suggests that calibration is not merely parameter fitting but also site selection. Electrical stimulation makes user-specific calibration even more explicit. The survey emphasizes threshold estimation, dynamic-range mapping, charge and charge-density constraints, and per-electrode calibration. Core formulas are $$ Q=I\cdot PW,\qquad D=\frac{Q}{A_{elec}},\qquad I_{stim}=I_{rh}\left(1+\frac{t_c}{PW}\right). $$ Typical cutaneous ranges are current amplitudes often below 1 mA at fingertips and operational ranges of 0–6 mA, pulse widths of 50–400 $\mu$s, and frequencies of 10–300 Hz; EMS typically uses longer pulse widths of 0.2–1.0 ms and frequencies of 20–100 Hz [2504.21477]. Because thresholds increase with age, differ by site, and drift with impedance, calibration must be closed-loop or periodically repeated [2504.21477]. Psychophysical calibration can also be task-specific rather than purely sensory. In digital musical instruments, TorqueTuner resets the current knob angle as the zero reference on each mode change, then parameters such as spring stiffness $k_s$, damping $b_s$, detent amplitude $A_d$, and detent spacing $\Delta\theta=45^\circ$ are tuned against subjective criteria including comfort, flexibility, ease of control, and helpfulness [2405.10502]. The reported association between musical background and preferred haptic mode—77.8% of wind instrument players preferring Spring mode and 62.5% of string players preferring Detent mode—indicates that calibration may legitimately incorporate user history rather than only biomechanics [2405.10502]. ## 5. Dynamic compensation, rendering laws, and latency Calibration is incomplete without a rendering law that remains stable under motion, delay, and model mismatch. TCR uses distance matching: the robot end-effector is commanded so that $$ d_v(t)=\|x_f^v(t)-x_o^v(t)\|,\qquad d_r(t)=\|x_f^r(t)-x_o^r(t)\|,\qquad d_r(t)=d_v(t). $$ With the front-of-user constraint, the command becomes $$ x_o^{r,cmd}(t)=x_f^r(t)+\alpha(t)\hat v,\qquad \alpha(t)=d_v(t), $$ or, equivalently, $$ x_o^{r,cmd}(t)=R x_f^v(t)+t+d_v(t)\hat v. $$ The dot-product condition keeps the end-effector in front of the user for safety [2309.16768]. In sim-to-real transfer, the calibrated mapping is inserted upstream of inverse kinematics: $$ C_S=T(x_d),\qquad q_s^*=IK(C_S). $$ The paper compares an affine-plus-interpolation baseline $M1$, a partially nonlinear network $M2$, and a fully nonlinear network $M3$ that predicts $(x_s,y_s,z_s)$ directly [2507.08572]. The reported cross-validation for $M3$ is MSE $<0.001$ in normalized units, and online correction is obtained by replacing the nominal simulated target with the calibrated one before execution [2507.08572]. Electrovibration work makes the compensation law explicit. Given the forward model $H_{f-e}(s,\nu)=K/(1+s/\omega_o(\nu))$, the proposed inverse is $$ H_{f-e}^{-1}(s,\nu)=\frac{1+s/\omega_o(\nu)}{K}, $$ so that the message voltage for a desired friction spectrum is $$ V_m^*(s,\nu)=\frac{F_f^*(s)}{K}\left(1+\frac{s}{\omega_o(\nu)}\right). $$ This is a speed-aware equalization strategy rather than a purely empirical lookup [2505.11162]. In robotic surgery training, force rendering uses a direction-preserving saturating nonlinearity, $$ {}^Hf_{scaled}=\left(\frac{{}^Hf}{\|{}^Hf\|}\right)\tanh\left(\frac{\|{}^Hf\|}{f_{scale}}\right)f_{max}, $$ with $f_{max}=3$ N and $f_{scale}=7$ N [2604.27385]. The stated rationale is high sensitivity in the sub-1–2 N regime and smooth saturation below the Touch device capability of approximately 3.3 N [2604.27385]. A different stability strategy appears in deformable simulation. The mesh-deformation framework decouples a 1000 Hz haptic thread from a 70–900 Hz visual/physics thread and relies on continuous collision detection, continuous penalty forces, and a local damping kernel $$ G_d(x)= \begin{cases} \frac{1}{1+k_1 r}, & r