Haptic shared control systems that support drivers by means of added torques on the steering wheel are often tuned heuristically. To allow for more systematic design, this paper focuses on the Four Design Choices Architecture (FDCA) and systematically analyzes its tuning with an offline simulation model for the driver’s control behavior and neuromuscular system. These analyses indicated that within the FDCA architecture the Level of Haptic Support (LoHS), which is a feedforward channel supporting negotiation of upcoming curves, is a main contributor to joint system performance. In a driving simulator experiment, the adaptation to and acceptance of different LoHS levels was investigated. Driver acceptance was found to increase with increasing LoHS values up to 1. Objective metrics, including torque conflict (70% reduction), steering effort (81% reduction), steering wheel reversal rate, and lateral deviation all improved, indicating that with the FDCA a high LoHS is both acceptable and, in fact, preferred.

Haptic Shared Control (HSC) systems offer a means to support human drivers in the transition to fully-automated driving. Matching HSC systems settings with drivers’ time-varying neuromuscular system (NMS) dynamics requires real-time HSC adaptations. This paper presents an experimental validation of a previously proposed method for predicting drivers’ time-varying neuromuscular admittance using an ‘average’ grip force scheduled linear parameter varying (LPV) model. The quality of LPV model predictions is compared to that of Recursive Least Squares (RLS) fits of an admittance model on the same data. Ten participants performed steering wheel manipulation tasks with steering wheel perturbations that needed to be kept within a certain displacement boundary by adapting their grip force. Time-invariant (TI) and time-varying (TV) boundary levels were used to, respectively, construct and validate the LPV model. Results show that the average relation between admittance and grip force that underlies the current LPV method varies too much between TV and TI tasks, hampering accurate admittance predictions. Compared to the quality-of-fit of 80-90% obtained with RLS on the TV data, the LPV model’s predictions are insufficiently accurate and do not exceed 55% on average. An approach that enables individual instead of average LPV models to be constructed directly from TV experiment data needs to be pursued for HSC implementations.

A pursuit-tracking manual control model is introduced that includes an observer-like internal model to predict human detection of a change in controlled element dynamics. The internal model’s innovation signal, the difference between the observed and expected system response, is studied for its capacity to drive the detection of a change. The model’s performance is tested for different crossover frequencies, remnant power ratios, observer gains, and detection threshold settings, through Monte Carlo analysis of simulated pursuit-tracking tasks where the controlled element transitions from single to double integrator dynamics. The model shows highly accurate detection performance for a wide range in the observer gain, with a true positive rate of approximately 1 and a false positive rate of approximately 0.02. The high true and low false positive rates, combined with average detection times that match experimental human-in-the-loop data, show the observer model’s potential for accurately predicting human detection of a change in controlled element dynamics.

InceptionTime neural network models were trained to detect distractions in manual control tasks with pursuit and preview displays. Training and test data were collected in an experiment where ten participants were deliberately distracted from the primary control task using the Surrogate Reference Task. Overall, distractions are easier to detect in pursuit tasks, with test accuracies of around 80% and 60% for pursuit and preview data, respectively. With preview, human controllers see the future target trajectory, which enables them to mitigate distraction effects. Unexpectedly, data with longer distractions from ‘hard’ secondary tasks are more difficult to classify than ‘easy’ distractions; an effect attributed to differences in human behavior between the training and test data collection conditions. These results show clear opportunities for neural network models to detect distractions, in real-time, for increasing safety of human-operated vehicles.

Electrovibration technology enables tactile texture rendering on capacitive touchscreens by modulating friction between the finger and the screen through electrostatic attraction forces, generated by applying an alternating voltage signal to the screen. Accurate signal calibration is essential for robust texture rendering but remains challenging due to variations in sliding speed, applied force, and individual skin mechanics, all of which unpredictably affect frictional behavior. Here, we investigate how exploration conditions affect electrovibration-induced finger friction on touchscreens and the role of skin mechanics in this process. Ten participants slid their index fingers across an electrovibration-enabled touchscreen at five sliding speeds (20∼100 mm/s) and applied force levels (0.2∼0.6 N). Contact forces and skin accelerations were measured while amplitude modulated voltage signals spanning the tactile frequency range were applied to the screen. We modeled the finger-touchscreen friction response as a first-order system and the skin mechanics as a mass–spring-damper system. Results showed that sliding speed influenced the friction response's cutoff frequency, along with the estimated finger moving mass and stiffness. For every 1 mm/s increase in speed, the cutoff frequency, the finger moving mass, and stiffness increased by 13.8 Hz, 3.23×10⁻⁵ kg, and 4.04 N/m, respectively. Correlation analysis revealed that finger stiffness had a greater impact on the cutoff frequency than moving mass. Notably, we observed a substantial inter-participant variability in both finger-display interaction and skin mechanics parameters. Finally, we developed a speed-dependent friction model to support consistent and perceptually stable electrovibration-based haptic feedback across varying user conditions.

This paper applies model-based biodynamic feedthrough (BDFT) cancellation to a touchscreen dragging task during realistic vertical (vertical) and lateral (horizontal) aircraft turbulence, to mitigate erroneous turbulence-induced inputs. One-size-fits-all (OSFA) BDFT models were used to model the influence of turbulence accelerations on finger position, achieving average quality-of-fits of 61% and 69% in the vertical and horizontal screen directions, respectively. On average, 27% of the touch input error variance was mitigated these OSFA models, with individualized models providing only a marginal improvement (+4%). The application of OSFA models identified from a condition with equally-scaled vertical and horizontal motion (adjusted intensity) to the realistic turbulence condition did not significantly affect cancellation performance, indicating that BDFT models may not need to be adaptive to varying motion intensity. However, consistent with earlier work, BDFT dynamics were found to vary between vertical and horizontal finger movements, with BDFT dynamics exhibiting lower stiffness and a higher static gain for vertical BDFT. On average, the linear BDFT-related component of touch input errors contributed 41% of the overall error variance, indicating that current linear BDFT model may need to be extended to include nonlinear effects, such as varying finger friction.

Aerodynamic stall has been a critical factor in recent aircraft crashes, leading to revised regulations for simulator-based stall prevention and recovery training. However, the updated regulations still lack an objectively defined level of accuracy for simulators' stall models that ensures effective pilot training. To help determine this required accuracy, this paper investigates how the Just Noticeable Difference (JND) thresholds for deviations in a stall model's ‘stall abruptness’ parameter translate from a passive observer setting (typical JND experiments) to an active flying scenario (realistic training task). An experiment was performed in the SIMONA Research Simulator with 16 active pilots, whose sensitivity to stall abruptness variations was measured in both these scenarios. In the passive scenario, pilots' JND thresholds were measured using a staircase procedure in a symmetric stall maneuver flown by a stall autopilot. In the active scenario, the method of constant stimuli was used to determine the JND thresholds when the pilots themselves actively flew the same stall maneuver. The average JND thresholds for the passive and active scenarios were estimated by fitting a psychometric curve to the combined responses of all participants. Overall, the passive JND thresholds for the stall abruptness parameter, with an average Weber fraction of 0.11±0.094, were lower than those measured in an earlier experiment (0.16±0.14), indicating a higher sensitivity. Furthermore, the psychometric curve of the active experiment was found to lie entirely to the right of the passive psychometric function: the active JND threshold was found to be five times higher than the passive JND threshold. Overall, this indicates a decreased sensitivity to changes in stall abruptness -- and hence a reduced demand on its modeling accuracy -- when pilots are flying a stall themselves.

While human control behavior is well-understood in continuous control tasks, little is still known about how human operators detect sudden changes in the controlled element dynamics. This paper focuses on modeling this detection phase for pursuit tracking tasks. Potential triggers for the human operator to detect changes in the controlled element dynamics were investigated via a time-varying computer simulation. Based on the results, hypotheses were generated and later tested in a single-axis pursuit tracking experiment with fifteen participants. Transitions from approximate single to approximate double integrator dynamics and vice versa were investigated, for which participants indicated if they detected the transition by pressing a button. Using the button push data, a model for each transition was developed and validated. The models work under the assumption that human operators use a threshold, a multiple of the steady-state standard deviation, on certain signals to detect transitions. The models developed for the transition from single to double integrator dynamics and vice versa are proposed to trigger on the tracking error and system output acceleration, respectively. They have an accuracy of 88.9% and 99.4%, respectively. However, a consistent underestimation of the detection lag remains a limitation of both models. Nonetheless, this research helped confirm the tracking error can be used in a model for the transition from single to double integrator dynamics, proposed a model for the opposite transition, and identified that the relationship between control inputs and the system's response as a crucial factor for the detection.

The identification of time-varying, adaptive behavior of a human operator in basic manual control tasks is currently still a focus area, since most methodologies only account for time-invariant system dynamics. Previous authors have proven that estimation techniques based on ARX model structures can be used to identify time-varying HO model parameters. However, ARX methods do present several problems, such as a persistent bias in the obtained estimates of the HO model poles (neuromuscular parameters) that increases due to coupled noise and system models. Therefore, in this paper a novel identification technique based on Box-Jenkins (BJ) models is proposed, to achieve a better match between the BJ estimator's inherently uncoupled system and noise models and measured HO control dynamics. The identification process was tested offline (batch-fitting) using Ordinary Least Squares and the Prediction Error Method for both ARX and BJ models, respectively, or online when Recursive Least Squares and Recursive PEM are employed. The BJ estimator has excellent potential as an identification tool due to its bias reduction capabilities, as clearly shown in batch-fitting, although the non-linear optimization processes decrease its convergence speed by 500%. An RPEM algorithm with a forgetting factor of λ = 0.99609 and a first-order remnant model incorporated in the BJ structure was tested on Monte Carlo simulation and experimental data. While the recursive BJ estimator showed the same bias-diminishing advantages also seen in batch-fitting, the non-linear RPEM estimator's results showed much slower convergence after HO behavior adaptations and frequent instabilities of the obtained parameter estimates. Hence, further research is needed for implementing a practical bias-free HO model estimator based on the BJ model structure.

As users transition from drivers to passengers in automated vehicles, they often take their eyes off the road to engage in non-driving activities. In driving simulators, visual motion is presented with scaled or without physical motion, leading to a mismatch between expected and perceived motion. Both conditions elicit motion sickness, calling for enhanced vehicle and simulator motion control strategies. Given the large differences in sickness susceptibility between individuals, effective countermeasures must address this at a personal level. This paper combines a group-averaged sensory conflict model with an individualized Accumulation Model (AM) to capture individual differences in motion sickness susceptibility across various conditions. The feasibility of this framework is verified using three datasets involving sickening conditions: (1) vehicle experiments with and without outside vision, (2) corresponding vehicle and driving simulator experiments, and (3) vehicle experiments with various non-driving-related tasks. All datasets involve passive motion, mirroring experience in automated vehicles. The preferred model (AM2) can fit individual motion sickness responses across conditions using only two individualized parameters (gain K1 and time constant T1) instead of the original five, ensuring unique parameters for each participant and generalisability across conditions. An average improvement factor of 1.7 in fitting individual motion sickness responses is achieved with the AM2 model compared to the group-averaged AM0 model. This framework demonstrates robustness by accurately modeling distinct motion and vision conditions. A Gaussian mixture model of the parameter distribution across a population is developed, which predicts motion sickness in an unseen dataset with an average RMSE of 0.47. This model reduces the need for large-scale population experiments, accelerating research and development

Providing adequate simulator motion cues for simulated upset and stall scenarios remains challenging. This paper evaluates the potential of novel optimization-based motion cueing algorithms for upset and stall simulation. An offline analysis is performed to compare three Model Predictive Control (MPC) algorithms with varying prediction horizon lengths and prediction strategies (i.e., "Oracle", "Perfect", and "Constant") against two baseline classical washout algorithm implementations from literature. The analysis is performed for a symmetric stall scenario flown with TU Delft's Cessna Citation II laboratory aircraft. Overall, the analysis shows that the objective motion cueing quality expressed in terms of the Root Mean Square Error (RMSE) improves by 29.8% (for specific forces) and 18.7% (for rotational velocities) with the "Oracle" and "Perfect" MPC implementations compared to the reference classical washout results. For the "Constant" MPC algorithm, which in fact does not include any explicit prediction across the MPC algorithm's future prediction window, only a marginal improvement in motion quality was found. Overall, these results imply that, assuming a sufficient future reference motion prediction can be achieved, optimization-based motion cueing algorithms have the potential to provide significantly better motion cueing quality compared to classical motion cueing algorithms.

An online pilot manual control behavior identification method, based on recursive low-order time-series model estimation, is presented and validated using experimental data. Eight participants performed compensatory tracking tasks with time-varying vehicle dynamics, where, at an unpredictable moment during a run, a sudden degradation in dynamics could occur. They were instructed to push a button when they detected a change in dynamics. Two methods to automatically detect the moment when pilot adaptation occurs from online estimated parameter traces are discussed. Results show that pilots are more accurate in detecting changes than either algorithm. But when the algorithms are correct, they are often quicker to detect pilot adaptation than pilots themselves. The presented techniques have potential but need improvements.

Recent aircraft have seen the implementation of touchscreens (TSCs) on the flight deck, as they enable more intuitive and direct human-machine interactions. However, biodynamic feedthrough (BDFT), i.e., the direct transmission of the aircraft's accelerations through the pilot's body to the control inputs, is a cause for concern, preventing safe and reliable use of TSCs in turbulence. This paper describes a simulator experiment evaluating the performance of model-based mitigation of BDFT in a TSC dragging task performed in turbulence. In the experiment, a total of nine different vertical (heave) motion perturbations were tested: multisine signals resembling turbulence, stationary (Gaussian) and variable (patchy) simulated turbulence, each at three intensity levels (RMS acceleration of 0.75, 0.5, and 0.25 m/s2). For the multisine turbulence signals, on average over 87% accuracy of the identified personalized BDFT models was achieved for the high and medium turbulence levels, reducing to 74% for the low-intensity turbulence due to degraded BDFT consistency. Furthermore, BDFT models fitted to the Gaussian turbulence data were found to achieve an accuracy comparable to that observed for the multisine motion disturbances, with only 3.5% lower performance on average. As expected, for the more time-varying patchy turbulence cases, model-based BDFT cancellation was found to be 4.7% lower than for the Gaussian turbulence data. Finally, models generalizing BDFT dynamics across participants or experimental runs were found to always be outperformed by individual participant and individual trial models, giving up to 10% higher identification performance. Overall, these findings show that a model-based approach to canceling the effects of BDFT mitigation for TSCs in turbulence is promising, but that real-time identification and time-varying BDFT models will be needed to achieve consistently high mitigation performance in realistic variable turbulence.

This paper investigates the effects of motion mismatches on simulator sickness and subjective ratings of the motion. In an open-loop driving simulator experiment, participants were driven through a recorded urban drive twelve times, in which mismatches were induced by manipulating the following three aspects in motion cueing: (i) mismatches in specific vehicle axes, (ii) mismatch types (scaling, missing, and false cues), and (iii) inconsistent scaling between different motion axes. Subjects (N=52) reported simulator sickness post-hoc (after each drive), as well as continuously during each drive, a first in simulator sickness research. Furthermore, subjective post-hoc motion incongruence ratings on the quality of the motion were extracted. Results show that longitudinal motion mismatches lead to the most simulator sickness and the highest ratings, followed by mismatches in lateral motion, then yaw rate. False cues induce the most sickness, followed by missing and then scaled motion. Inconsistent scaling between the axes has no significant effect. The continuous sickness ratings support that the occurrence and severity of simulator sickness are indeed related to mismatches in simulator motion of specific maneuvers. This paper contributes to an improved understanding of the relationship between simulator motion and sickness, allowing for more targeted motion cueing strategies to prevent and reduce sickness in driving simulators. These strategies may include the appropriate selection of the simulator, the motion cueing, and the sample of participants, following the presented results.

This paper presents a three-step validation approach for subjective rating predictions of driving simulator motion incongruences based on objective mismatches between reference vehicle and simulator motion. This approach relies on using high-resolution rating predictions of open-loop driving (participants being driven) for ratings of motion in closed-loop driving (participants driving themselves). A driving simulator experiment in an urban scenario is described, of which the rating data of 36 participants was recorded and analyzed. In the experiment's first phase, participants actively drove themselves (i.e., closed-loop). By recording the drives of the participants and playing these back to themselves (open-loop) in the second phase, participants experienced the same motion in both phases. Participants rated the motion after each maneuver and at the end of each drive. In the third phase they again drove open-loop, but rated the motion continuously, only possible in open-loop driving. Results show that a rating model, acquired through a different experiment, can well predict the measured continuous ratings. Second, the maximum of the measured continuous ratings correlates to both the maneuver-based (ρ =0.94) and overall (ρ =0.69) ratings, allowing for predictions of both rating types based on the continuous rating model. Third, using Bayesian statistics it is then shown that both the maneuver-based and overall ratings between the closed-loop and open-loop drives are equivalent. This allows for predictions of maneuver-based and overall ratings using the high-resolution continuous rating models. These predictions can be used as an accurate trade-off method of motion cueing settings of future closed-loop driving simulator experiments.

One of the most widely applied identification methods for stall modeling using flight test data is based on Kirchhoff’s method of flow separation. However, this approach has not lead to a satisfactory aerodynamic pitching moment model. The introduction of the so-called X-variable, representing the point of flow separation on the wing, interferes with identification of a pitch damping term, that is required for dynamic stability. In general, Kirchhoff methods lead to models that are incompatible with nominal flight envelope models. This paper presents a nonlinear unsteady model of the pitching moment using lag states of the angle of attack measurements, identified from flight test data collected with a Cessna Citation II laboratory aircraft. The model is formulated in terms of well-known stability derivatives and is a one-on-one extension of the nominal envelope model. Model regressors are selected from a large pool of candidates using Multivariate Orthogonal Function Modeling. The candidate pool is based on a newly formulated mathematical model, such that each model contribution has a clear physical interpretation. The model has good predictive abilities and results in a reduction of 55.9% in validation MSE compared to Kirchhoff based pitching moment models.

Objective
This study investigates the impact of whole-body vibrations caused by external vehicle perturbations, such as aircraft turbulence, on the perception of electrovibration displayed on touchscreens.

Background
Electrovibration is a promising technology for providing tactile feedback on future touchscreens, potentially addressing usability challenges in vehicle cockpits. However, its performance under dynamic conditions, such as whole-body vibrations caused by turbulence, remains largely unexplored.

Method
We measured the absolute detection thresholds of 24 human participants for short (0.2 s) and long (0.5 s) duration electrovibration stimuli displayed on a touchscreen. These measurements were taken in the absence and presence of two types of turbulence motion (Gaussian and Multisine) generated by a motion simulator. Concurrently, we recorded participants’ applied contact force and finger displacements.

Results
Electrovibration stimuli displayed on vehicle cockpit touchscreens were more reliably perceived with a 0.5-s duration than a 0.2-s duration, both in the presence and absence of turbulence. Both turbulence types led to increased vibration-induced finger displacements and scan speeds in the direction of turbulence, as well as higher applied forces and force fluctuation rates. Gaussian turbulence significantly elevated perception thresholds, but only for short-duration electrovibration stimuli.

Conclusion
The findings indicate that whole-body vibrations impair the perception of short-duration electrovibration stimuli, primarily due to unintentional finger movements and increased fluctuations in applied normal force.

Application
Our findings offer valuable insights for the future design of touchscreens with tactile feedback in vehicle cockpits.

To improve the safety of commercial air transport, pilots are required to train on simulators to recognize the characteristics of an impending stall and subsequently correctly recover from it. To prevent negative training, it is important that the accuracy of the used simulation models is sufficiently high. A key approach for modeling the nonlinear, unsteady aerodynamic effects during the stall is by using Kirchhoff's theory of flow separation. However, widespread difficulties exist in correctly estimating the stall-related parameters of nonlinear flow separation models from flight test data. Therefore, the research in this paper aims to increase the obtained model accuracy by making optimal use of already existing flight data via introduction of a slice-based modeling method. This is done by analyzing the change in the parameter estimate values when applying the system identification procedure to sliced partitions of simulated flight data, for both the pre-stall and post-stall phases. These partitions incrementally increase in size with time from the stall initiation. The simulation data is generated to be representative of the available flight test data, but with known ‘truth’ values for all estimated model parameters. The estimated value for each partition was compared to the true parameter setting in the simulation model used to create the data. It was also investigated whether this coincided with points of increased Fisher information in the data. Manually, an optimal window was found for each parameter for which the estimated value and truth value were equal and sufficient Fisher information was present. For the stall-related parameters the optimal window is often not more than 10 s wider than the stall. For the linear stability and control derivatives it is found that using more data generally results in a better estimate. Finally, the optimal window sizes were used for parameter estimation on the real flight test data. Even though this method represents a prototype, in more than half of the validation cases a decrease in MSE of 10% to 35% was achieved. This shows that the new slice-based modeling method is able to improve the accuracy of nonlinear stall models without the need to gather more flight data and may have applications that reach beyond the realm of stall modeling.

This erratum applies to the following published paper [1]. In Fig. 10(e) and (f) of the published version of the paper, the measured values for the t _f and T _l,f (both having values around 1 s for the considered dataset) were interchanged. This erratum includes both the published and corrected versions of Fig. 10 and the related paragraph in the paper's Results section. PUBLISHED VERSION In preview tasks, bandwidth changes yield only minor adaptations in the model parameters, see Fig. 10. Only the average look-ahead time t _f decreases slightly with bandwidth from around 1.05 to 0.9 s, Fig. 10(e) but with substantial between-participant variability, as indicated by the overlapping confidence intervals. The lower t _f may not reflect a systematic adaptation to the bandwidth, but a more subtle adaptation to minimize the errors due to the additional high-amplitude sinusoids at 2.5 and 4 rad/s, see also Fig. 8c. The general way in which participants use the available preview for control is, however, not affected by the target signal bandwidth: the target response gain (K _f ≈ 0.95, Fig. 10(d)) and lag time-constant (T _{l, f} ≈ 1.15 s, Fig. 10(f)) are approximately invariant. The estimated control dynamics in Fig. 12 show that the target trajectory is tracked almost perfectly at all frequencies below 4 rad/s, mostly because the phase lead due to τ _f allows for synchronizing the CE output with the target signal (as opposed to pursuit tasks, see Fig. 11, bottom right). Therefore, different-bandwidth target signals provide no incentive for HCs to strongly adapt their control behavior in preview tasks. (Figure Presneted) CORRECTED VERSION In preview tasks, bandwidth changes yield only minor adaptations in the model parameters, see Fig. 10. Only the average target smoothing time-constant T _{l, f} decreases slightly with bandwidth [from around 1.05 to 0.9 s, Fig. 10(f)], but with substantial between-participant variability, as indicated by the overlapping confidence intervals. The lower T _{l, f} indicates that slightly more smoothing is applied to reduce tracking of the more high-amplitude high-frequency sinusoids in the 2.5 and 4 rad/s bandwidth signals through the feedforward response, see also Fig. 8(c). The general way in which participants use the available preview for control is, however, not affected by the target signal bandwidth: the target response gain [K _f ≈ 0.95, Fig. 10(d)] and look-ahead time [τ _f ≈ 1.15 s, Fig. 10(e)] are approximately invariant. The estimated control dynamics in Fig. 12 show that the target trajectory is tracked almost perfectly at all frequencies below 4 rad/s, mostly because the phase lead due to τ _f allows for synchronizing the CE output with the target signal (as opposed to pursuit tasks, see Fig. 11, bottom right). Therefore, different-bandwidth target signals provide no incentive for HCs to strongly adapt their control behavior in preview tasks. (Figure presented).

The collective goal of the driving simulation community should be to share ideas to improve the motion cueing across driving simulators worldwide. Due to the active research and intensive usage of driving simulators over the last decades, knowledge in the field of motion cueing has been gained from experience gathered by the community, or practical experience by performing dedicated experiments. This paper discusses several points of ‘common knowledge’ in designing and evaluating motion cueing, along with their value for driving simulation. The goal of the discussion in this paper is to compare these points of common knowledge to the experiences and ideas gathered at BMW’s driving simulation center in Munich, which was opened in 2021 and hosts a fleet of fourteen driving simulators. Furthermore, we aim to bring across points of interest and outlines for future research that should be of interest to the driving simulation community. With the common goal of improving motion cueing, this contribution thus aims to improve and extend the discussion between those working and researching in the driving simulator industry.