Multibody System Dynamics (2023) 60:289–316 https://doi.org/10.1007/s11044-023-09879-8 R E S E A R C H Efficient simulation strategy to design a safer motorcycle Steffen Maier1 · Jörg Fehr1 Received: 24 October 2022 / Accepted: 10 January 2023 / Published online: 7 February 2023 © The Author(s) 2023 Abstract This work presents models and simulations of a numerical strategy for a time and cost- efficient virtual product development of a novel passive safety restraint concept for mo- torcycles. It combines multiple individual development tasks in an aggregated procedure. The strategy consists of three successive virtual development stages with a continuously increasing level of detail and expected fidelity in multibody and finite element simulation environments. The results show what is possible with an entirely virtual concept study— based on the clever combination of multibody dynamics and nonlinear finite elements—that investigates the structural behavior and impact dynamics of the powered two-wheeler with the safety systems and the rider’s response. The simulations show a guided and controlled trajectory and deceleration of the motorcycle rider, resulting in fewer critical biomechanical loads on the rider compared to an impact with a conventional motorcycle. The numerical research strategy outlines a novel procedure in virtual motorcycle accident research with different levels of computational effort and model complexity aimed at a step-by-step vali- dation of individual components in the future. Keywords Multibody system · Nonlinear finite elements · Crashworthiness · Passive safety · Impact biomechanics 1 Introduction A motorcycle accident is a very complex event where the rider and the motorcycle interact with many environmental factors. Due to the exposed position of the riders, the vehicle itself does not provide protection to the rider in the event of a collision with an accident opponent or a roadside structure. Instead, violent ejection of the rider from the motorcycle is a likely accident pathway when the vehicle suddenly comes to a standstill [1–3]. When striking objects in their path and the ground, the consequences are often severe or fatal for the involved riders. This results in an excessive risk that motorcycles are 25 times more deadly per kilometer traveled than passenger cars, see Fig. 1a. Accident analyses of severe or fatal injuries [4] show that the personal protective equipment worn is insufficient to mitigate the consequences of serious accidents, which motivates the first research question RQ1: Would � J. Fehr joerg.fehr@itm.uni-stuttgart.de 1 Institute of Engineering and Computational Mechanics, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany http://crossmark.crossref.org/dialog/?doi=10.1007/s11044-023-09879-8&domain=pdf mailto:joerg.fehr@itm.uni-stuttgart.de 290 S. Maier, J. Fehr it be safer to be restrained to a motorcycle? This will be investigated with a motorcycle with a novel safety concept where the rider is restrained by thigh belts, multiple airbags, and leg impact protection. The concept’s goal is to guide the rider’s impact trajectory for controlled deceleration and to prevent direct impacts on opposing vehicle structures. A safety system that effectively protects the rider must function robustly in the many pos- sible accident scenarios. The possible scenarios for accidents involving motorcycles are very diverse. Because of their complex nature, it is very difficult to find passive safety solutions that are effective in as many accident configurations as possible and that at least have neu- tral or nondetrimental effects on the accident consequences in off-design accident scenar- ios. This results in many design iterations and cannot be achieved with costly experimental methods alone. Therefore, the merits of virtual methods, like multiquery and multiparameter design studies, are needed. It includes procedures that allow for a systematic approach to the design of the proposed passive safety systems and the investigation of their protective per- formance in representative accident scenarios. To test many different solutions efficiently, models with varying levels of complexity and computing effort that capture different aspects of the accident behavior are necessary, which leads to the main subject of this work, the sec- ond research question RQ2: What virtual models and methods meet the above requirements to investigate many possible design iterations and accident configurations effectively and efficiently? For this, this work presents a virtual simulation-based product design procedure consisting of three successive development stages with a continuously increasing level of detail and expected fidelity. The novelty of this work is the modeling and simulation strategy for crashworthiness investigations of a motorcycle with a systematic model generation approach with different levels of computational effort and model complexity. It outlines a novel procedure in virtual motorcycle accident research and passive safety equipment development that is very similar to common strategies in the development process of passenger cars for occupant protection. First, this work describes the aims and concept of the novel safety concept. It briefly in- troduces the essential computational models in current crashworthiness investigation, iden- tifying and substantiating the choice for the modeling and simulation strategy. An emphasis of this work is on describing the modeling of the vehicle, passive safety systems, and rider surrogates. The results focus on a meaningful description of the operating principles and their influence on accident behavior based on the virtual models. The advantages of the individual modeling stages are discussed. 1.1 Novel safety concept for motorcycles Powered two-wheelers (PTWs), or more commonly called motorcycles, are a popular means of transportation. They are attractive as a relatively cheap, agile, and space-saving mobility solution. Annual new registrations of motorcycles and mopeds, for example, in Europe had a general positive trend from 2013 to 2019 [5, 6] and many European countries exceeded pre-pandemic levels in 2021 [7]. However, because of their low level of passive protec- tion, they carry a disproportionate risk of severe or fatal injury for their riders in accidents compared to cars. The World Health Organization estimates that every year about 375,000 users of two and three-wheelers die in road traffic accidents [8]. This is 28% of all road traffic deaths. Only fatalities of riders and passengers of four-wheeled vehicles with a share of 29% were marginally more frequent, although these are used far more. A more detailed analysis of reported accidents in Germany illustrates the immense risk. As shown in Fig. 1a by accounting for the mileage for motorcycles (excluding light PTWs, such as mofas and mopeds), the risk of a fatal accident is more than 25 times higher than for cars, buses, or Efficient simulation strategy to design a safer motorcycle 291 Fig. 1 Road-traffic fatalities per kilometer traveled in Germany [9, 10] (Note: Figures for 2020 are also available, but according to [9], the mileage for 2020 is subject to great uncertainty due to the coronavirus pandemic) Fig. 2 Safety systems of the novel safety concept for motorcycles (Note: The right airbags are not displayed). A video of the safety concept is shown in [11] trucks, individually—the trend since 1991 in Fig. 1b illustrates that the risk of a serious in- jury with 920 vs. 43.9 serious injuries per billion kilometers is likewise much higher and the safety of motorcycles has not improved nearly as much as that of passenger cars in the past decades. The proposed novel safety concept shown in Fig. 2 aims to improve rider passive safety in accidents significantly. It consists of a newly designed motorcycle body, seat belts, multi- ple surrounding airbags, foam leg impact protectors, and side impact structure. The concept envisages that in the event of an accident, the rider is restrained to the motorcycle by two thigh belts. This causes the rider’s upper body to rotate around the belt restraint. The sur- rounding airbags then decelerate the upper body rotation in a controlled manner and protect the rider from hard contact with an accident opponent, the road, or road-side structures. The foam leg impact protectors absorb the impact of the legs on the motorcycle cockpit, and the side impact structure protects the lower extremities laterally. The lateral fairing is shaped so that comfortable mounting, riding, and balancing while stationary should still be 292 S. Maier, J. Fehr Fig. 3 Categories of mathematical models for investigating crashworthiness possible—which must be proved via an ergonomic analysis—while retaining the advantages of a motorcycle as a compact means of transportation. 1.2 Computational models for crashworthiness A recent review on mathematical models to assess vehicle impact crashworthiness [12] clas- sifies existing models into the five model categories shown in Fig. 3. (i) Reduced-order mod- els include lumped parameter models and commonly replicate vehicle impacts with lumped masses in interaction through simple mechanical elements such as springs and dampers. They aim to capture the fundamental impact kinematics, often for unidirectional response with very few degrees of freedom. (ii) Multibody (MB) models in crashworthiness are usu- ally more complex systems of rigid and flexible bodies, constrained by idealized joints and loaded by internal and external force elements. In generalized coordinates, the equation of motion for a holonomic MB system is M(y, t) · ÿ = kc(y, ẏ, t) + ke(y, ẏ, t), (1) where M ∈ R f ×f is the global mass matrix, ÿ ∈ R f ×1 is the acceleration vector, kc ∈ R f ×1 is the vector of the generalized Coriolis and centrifugal forces, and ke ∈ R f ×1 is the vector of the generalized applied forces. In addition, contact interactions of the bodies are often implemented via nonlinear contact characteristics also representing the deformation of the rigid bodies. A three-dimensional system of p rigid bodies with r holonomic, rheonomic constraints has a number of degrees of freedom of f = 6p − r [13]. MB models can be used for kinematic and kinetic analysis, including the structural response of vehicles [14] and the prediction of kinematic response [15] and biomechanical loading [16] of humans. (iii) Finite element (FE) methods in crashworthiness deal with geometry, material, constraint, and contact nonlinearities, making solving such models more difficult and computationally costly. The discretized equations of motion of an FE model with a Lagrangian mesh can be written as M · q̈ = f int − f ext, (2) where M ∈R n×n is the mass matrix and q̈ ∈R n×1 are the nodal accelerations that result from the internal forces f int ∈ R n×1, representing element stresses, and external forces f ext ∈ R n×1 from boundary conditions, such as acceleration fields, constraints, or contacts [17]. Depending on the element formulation, each node has different degrees of freedom, com- monly three or six. The models accurately represent complex geometries and incorporate diverse material models in the detailed representation of vehicles and human surrogates. They are particularly applied to capture vehicle structural behavior and localized effects such as detailed injury mechanisms in human body models. (iv) Crash pulses models are Efficient simulation strategy to design a safer motorcycle 293 substructure approaches [18] where the vehicle impact response is substituted, e.g., from experiments or full vehicle simulations [19]. It reduces the degrees of freedom of high fi- delity models by strategically replacing constituents with lower fidelity physical surrogate models. (v) As its mathematical equivalent, response surface models are mathematical sur- rogate models with usually even fewer degrees of freedom. They describe relationships be- tween very few input and output variables obtained by various mathematical and statistical methods. These models aim to predict and improve the behavior of complex designs with minimized computational effort. Although not remotely comparable to crash safety of passenger car occupants, there have been efforts to model motorcycle crashworthiness. A rare reduced order model of motor- cycles is in [20]. There are several MB models of motorized two-wheelers to model car impacts [21, 22], road-side barrier impacts [23–25], and solo accidents [26]. There are also published works on detailed FE models of various motorcycles types, such as sports and sport touring motorcycles [27, 28], a large touring motorcycle [29], or a tilting three-wheeled scooter [30]. However, there is a lack of reduced models to efficiently investigate isolated aspects of motorcycle or motorcyclist crash behavior in submodel approaches. Also, there is a lack in combinations of two or more models, e.g., of different detailing and complexity, into a connected workflow. This has long been state-of-the-art in accident simulation of cars. 2 Modeling 2.1 Modeling and simulation strategy Common methods and strategies used in the vehicle development process of passenger cars for occupant protection break down the accident event into several individual problems, see, e.g., [31]. Thus, in experiments, the performance of the components, the interaction of the vehicle structure with the accident opponent, and the occupants’ behavior within the vehicle interior are often considered separately. Individual components or subassemblies are sub- jected to separate component tests to assess the functionality of the individual systems and the structural capacity of a component, substructure, and vehicle body. The interaction of the restraint systems with mechanical human surrogates, anthropometric test devices (ATDs), as vehicle occupants is optimized in sled tests. In these experiments, reinforced partial vehi- cle bodies represent the relevant vehicle’s interior to investigate the effect of the passive safety systems without destroying an entire vehicle. Full-vehicle laboratory crash tests of new products are rarely carried out during the design process but rather at the very end of the development to ultimately prove the occupant protection for vehicle approval or to evaluate occupant safety for consumer ratings. Similarly, in the virtual vehicle development process, there are multiple mathematical model representations used for individual design aspects, see Sect. 1.1. Corresponding to the experimental procedure mentioned above, the depth of reproduction in occupant safety studies is often reduced not to simulate the entire system. Thus, only the vehicle structures are considered for the reconstruction of the vehi- cle interaction and dimensioning of the deformation and contact structures. The occupant restraint interaction and biomechanical loading is simulated using crash pulses and vehicle interior models. Such a breakdown, both for experiments and simulations, is usually not possible when investigating accidents involving current conventional motorcycles, e.g., in collisions with passenger cars. Here, the rider interacts with the motorcycle, opposing vehicle structures, and road and road-side structures. A safety concept that restrains the motorcyclist to the 294 S. Maier, J. Fehr Fig. 4 Modeling and simulation strategy motorcycle and isolates it from the accident environment not only has the advantage of potentially improving passive safety, which must be demonstrated, but also allows to apply similar modeling and simulation strategies for a more systematic investigation. In Fig. 4 such a design strategy is developed. It consists of three subsequent development stages with a continuously increasing level of detail: – Stage 1: The motorcycle as well as a rider surrogate are modeled in a combined MB and FE approach in the Madymo software environment (version MADYMO 2020.1), as the authors introduced in [32]. The vehicles are MB systems with joint restraint and contact characteristics based on simulation models of full-scale crash tests of conventional motorcycles and fitted to the experience from the full FE representations (stage 3). The multiple airbags and the thigh belts are modeled with 1D and 2D FE elements using the FE capabilities of the Madymo software environment. This stage features few degrees of freedom and associated low numerical costs while capturing the essential physics of the impact. It is used to initially tune and improve the safety system. – Stage 2: The rider interaction surfaces of the motorcycle cockpit are modeled fur- ther detailed as an FE model in the LS-Dyna software environment (version LS- DYNA R9.3.1 MPP), as the authors developed and showed in [33, 34]. The modeling and parameterization of the airbags’ inflation and thigh belts’ pretensioning and load- limiting is equivalent to the MB simulations. The FE rider interaction model represents deformable cockpit surfaces that include foam impact protectors, which the authors char- acterized in [35]. To replicate the crash kinematics and vehicle intrusion, the motorcycle and car are (offline-)coupled to a crash pulse. Prescribed motions represent the multiaxial Efficient simulation strategy to design a safer motorcycle 295 Fig. 5 Kinematic structure of the motorcycle and accident opponent rigid body motions from the MB simulations where the car’s body geometry acts as re- action surfaces for the airbags. This stage allows using state-of-the-art FE HBMs as rider surrogates, e.g., to examine injuries on the tissue level, while reducing numerical costs by using prescribed vehicle motions. – Stage 3: The motorcycle, the already tuned passive safety systems equivalent to stage 1 and 2, the rider surrogate, and an accident opponent are modeled as full FE representa- tions in Ls-Dyna. The motorcycle’s structurally relevant components that determine the crash behavior are deformable. The front and rear suspension, rotating wheels, and front fork steering are modeled with kinematic joints. This stage allows modeling the vehicle response with a high degree of detail to investigate the deformation characteristics of the motorcycle itself and the structural interaction with opposing vehicles. It is used to predict the performance of the finalized design accurately. The resulting vehicle interactions are used to tune the rigid body interaction of the stage 1 model. The procedure is divided into modeling and simulation, in between design iterations and updates from model validation are carried out. As shown in this work, the development progress follows a sequence from left to right, with stage 3 as the final stage. Linking design updates for the motorcycle structure and passive safety systems to the lower stages results in design loops through stage 1 and 2 to update and maintain an identical design with uniform design parameters in all the simulation stages. However, the sequence could also be modi- fied. Thus, the vehicle kinematics from the full FE vehicle interactions could be applied to the stage 2 model, as indicated with the dashed arrow. 2.2 Vehicles 2.2.1 Multibody systems To replicate the crash dynamics of the motorcycle and an opposing vehicle, a passenger car, the main parts that determine their significant crash behavior, are modeled as rigid bodies. They are connected with kinematic joints as schematically illustrated in Fig. 5 with restraint characteristics as shown in Fig. 6. The bodies are the wheels and the vehi- cle’s body each, which are specified by their geometry with ellipsoid and facet surfaces, mass m, and inertia I . The kinematic structure replicates (i) rotating wheels, (ii) front and rear suspension, (iii) front fork steering, and (iv) front fork impact deformation. The sus- pension is coupled by nonlinear spring restraints cFW(s) and cRW(θ) with front suspension 296 S. Maier, J. Fehr Fig. 6 Joint restraint characteristics of the motorcycle and accident opponent Fig. 7 Loading and unloading characteristics Mdef(ϕ) of the telescopic front fork impact deformation deflection s and rear suspension deflection angle θ , and constant dampening restraints dFW and dRW. A rotational joint couples the front telescopic fork for frontal impact response, rep- resenting reward fork deformation. Loading and unloading response via moment Mdef(ϕ) from deflection angle ϕ, shown in Fig. 7, allows for energy dissipation through hysteresis. It represents linear response up to the elastic limit ϕe ( a©) and beyond plastic deforma- tion along Mdef,loading(ϕ) ( b©). Unloading is first parallel to the hysteresis slope ( c©) and then along Mdef,unloading(ϕ) ( d©). Reloading while unloading follows that path until reach- ing Mdef,loading(ϕ). The incorporated Madymo contact interactions characteristics are simi- lar hysteretic models with load functions of penetration depths of the contacting ellipsoid and facet surfaces. For a multibody motorcycle model with equivalent kinematic structure, see [21]. Besides, [36] describes a similar model without front fork deformation. The opposing vehicle is a 2011 Honda Accord, a four-door passenger sedan with a mass of 1,668 kg, an overall height of 141 cm, and vehicle parameters according to [37]. The given front fork deflection and contact characteristics are based on manually fitted simulations of full-scale crash tests of conventional motorcycles and a Hybrid III 50th ATD against passenger cars from Dekra [38]. For a complete discussion of the simulations of conventional motorcycles with riders, see [32]. A comparison of motion sequence and mo- torcycle acceleration sensor data for a frontal impact of a Yamaha FZS 600 Fazer impacting at 48.5 km/h at a right angle into the side of the stationary car is given in Fig. 8 and Fig. 9. The accelerations are filtered with a CFC filter (CFC for channel frequency class) with filter class 60, see [39]. The motorcycle and car trajectories and accelerations conform well, but Efficient simulation strategy to design a safer motorcycle 297 Fig. 8 MB simulation of full-scale crash test SH01.01 [38] with a conventional motorcycle Yamaha FZS 600 Fazer and a helmeted Hybrid III 50th ATD against a VW Golf II Fig. 9 Resultant linear motorcycle accelerations at multiple locations of MB simulation of full-scale crash test SH01.01 compared to experimental sensor data 298 S. Maier, J. Fehr Fig. 10 Coupled FE rider interaction model with rigid body trajectories from MB simulation similar helmet impact location at the car could not be achieved with the MB approach. This is because in the MB simulation the initial suspension compression of the car at the impact side (see snapshot at 100 ms) does not correlate. The real vehicle from the crash test has a softer door and stiffer door sill area, which likely causes the vehicle to roll to the impact side initially; see also Fig. 26. An approximation as a homogeneous body in the MB approach is not sufficient here. In a future step, feedback from stage 3 simulations could help to create a more fitting/better multibody model of the opposing vehicle. To quantitatively compare the motorcycle accelerations, the CORA (CORelation and Analysis) objective rating [40] is used. It offers an objective metric for the correlation of time history signals by combining two independent methods: a corridor rating and a cross- correlation rating. Using the stipulated default parameters for weights W , corridor widths a0 and b0, and rating exponents K of the rating method (version CORAplus Release 4.0.4 [41]), the CORA rating the agreement for the primary impact (0-300 ms) is in the very upper range of fair, very near good. The used rating parameter values are given in Fig. 9. 2.2.2 Coupled finite element and multibody model An FE model represents the rider interaction surfaces and the passive safety equipment for a further detailed step in the overall simulation strategy. 3D multiaxial rigid-body trajectories from the linear and angular motion from the MB simulations replicate the vehicle motion for crash dynamics by prescribing the outer body geometry of the motorcycle, see Fig. 10. For translational motion, the initial linear velocities and the time histories of linear acceler- ation are used, and for angular motion, the velocities are used. The safety concept envisages that when deployed, the side airbags interact with the opposing vehicle structures or the ground. Therefore, the prescribed outer surface of the accident opponent represents the re- action surface. The overall procedure corresponds to experiments or simulations of vehicle occupants in which impact pulses also reproduce the accident interaction with an accident opponent. 2.2.3 Full finite element models For the highest degree of detail, a full FE representation of the motorcycle is developed from a concept design CAD model that represents the main construction features of the motor- cycle. The CAD model comprises standard–“off the shelf”—components for motorcycles from various model libraries and many custom parts. For investigations of crashworthiness, the model aims to represent the interaction with the crash opponent, structural loading and Efficient simulation strategy to design a safer motorcycle 299 Fig. 11 FE modeling pre-processing workflow with joint coupling details deformation, and energy absorption of the motorcycle structure. As a result, the model fo- cuses on representing the crash-relevant and structural components, i.e., the front wheel, tire, and suspension assembly. Parts that are assumed not to deform because they are very stiff or outside of crash deformation, such as drivetrain components, are modeled rigid, but their dimensions, weight, and position are considered. Figure 11 illustrates the steps performed on the way to the discretized, coupled system using the telescopic front fork assembly as an example. It consists of cleaning up and simplifying the CAD geometry due to the initially relatively low model quality mentioned above and is followed by the discretization and connection of the components, as shown in section cuts of some revolute and translational joint couplings below. A unique feature of the motorcycle structure is a rigid foam crash box in the cockpit nose, which aims for controlled energy transfer to prevent rollover in a frontal impact and elevated side impact structures that protect the lower extremities laterally. The model consists of 148 parts from 377,000 elements with 323,000 nodes. It has seven kinematic joints: front wheel rotation (2), telescopic front fork suspension (2), rear-wheel rotation, front fork steering, and rear swing arm rotation. The opposing vehicle, an FE model of the 2011 Honda Accord four-door passenger sedan, is chosen from the publicly available vehicle model database of the National Highway Traffic Safety Administration (NHTSA) [37, 42]. 300 S. Maier, J. Fehr Fig. 12 Structure of 1D truss and 2D membrane hybrid mesh belt with combined pretensioning and retracting capabilities 2.3 Rider restraint 2.3.1 Thigh belt modeling The structure of the belt system is schematically illustrated in Fig. 12. The two belt sec- tions are four-node 2D shell and two-node 1D truss elements. The truss nodes have only translational nodal degrees of freedom, carrying only in-line loads. The shell nodes have translational and rotational degrees of freedom, transmitting forces and moments. The belts are connected at both ends to the attachment points to the motorcycle body (anchors). At their outer ends, nonlinear retractors with pre-tensioning and load-limiting are implemented to control the belt force f (t) by manipulating the belt length l(t) during impact. The retrac- tor removes any pre-tension or initial slack in the belt by supplying or retracting belt length. After the activation or firing, a pre-tensioning regime is determined by a function of belt pay-in vs. time Lpay−in(t) = ⎧ ⎪⎨ ⎪⎩ 0, if t < tfire − 100 15 mm ms (t − tfire), if tfire ≤ t ≤ tfire + 15 ms −100 mm, if t > tfire + 15 ms . (3) Simultaneously, the force level at the outer belt ends is limited by the load-limiter to Fload−limit = 2 kN in the connected truss element by giving out belt material. The firing of the systems is set to a fixed time of tfire =12 ms after an initial accident contact. This follows findings of sensor concepts for accident detection for motorcycles in [43], i.e., based on front fork deceleration and specifications of the Honda Gold Wing sensor system in [44], the only motorcycle to date that is equipped with a standard, although only optional, airbag introduced in 2006 [45]. The retractors with belt tensioners are currently modeled as point elements, i.e., the packaging is currently not accurately considered. As given in Fig. 13a and b for a frontal crash scenario with a motorcycle impact speed of 48.5 km/h against a stationary car (same as in Fig. 8), the force and belt length response can be divided into three consecutive phases: (i) pretensioning, (ii) belt pay-out because of load limiting, and (iii) a rigid restraint without any belt in- or output. Here, pre-tensioning Efficient simulation strategy to design a safer motorcycle 301 Fig. 13 Belt end force and length during pretensioning and load limiting in a frontal impact reduces the frontal displacement of the rider in frontal collisions through early intervention of the belt restraint and provides a constant force level, thus reducing peak forces and peak deceleration of the rider. Apart from initial belt routing, the design variables are the belt geometry, such as the width, the positions of the attachment points, the belt pay-in function, and the belt load limit. 2.3.2 Airbag modeling The three types of airbags of the safety concept are single chamber 2D and 3D designs with and without internal structure and with and without exhaust openings. For example, the front airbag is a 3D design, meaning that its geometry in its undistorted reference or design configuration encloses a volume compared to a flat 2D design. The frontal airbag has a wedge shape, no internal structures, and two exhaust openings to regulate the airbag pressure and, consequently, the rider’s deceleration when impacting. The windshield is designed to fail under contact pressure from the expanding front airbag structurally. There are multiple methods to model inflation of an airbag, see overviews in [46, 47]. These are uniform pressure concepts [48] and methods that discretize the gas flows us- ing mesh-based [49] or particle-based [50] approaches. The latter have advantages in re- producing the beginning of the airbag expansion. Therefore, these are especially impor- tant for more fundamental out-of-position investigations. In this work, the uniform pressure method (UPM) is sufficient and an appropriate choice as an investigation in the context of the early development phase. As shown in Fig. 14, the airbag is discretized with three- and four-node membrane elements, since bending stiffness is often assumed negligible for sim- ulating airbag fabrics. Expansion is determined by entering and exiting gas mass flows ṁ through the inflator, airbag exhausts, and airbag seam and fabric leakage. The mass balance complies for each time step i such as ṁi,tot = ṁi,in + ṁi,out 302 S. Maier, J. Fehr Fig. 14 Concept of the uniform pressure method (UPM) = ṁi,12 + ṁi,23,exhaust + ṁi,23,leakage. (4) The internal airbag chamber variables are calculated from scalar thermodynamic equations. At the beginning of each time step, the volume of the expanding airbag chamber V is deter- mined. In the basic principle, assuming uniform pressure p and temperature T of an ideal gas in the adiabatic airbag chamber, the pressure is expressed as pi = (γ − 1) ρi ei, (5) where γ is the isentropic coefficient cp/cv , ρ is the density, and e is the specific internal energy. The dependency for two neighboring steps is ei ei−1 = ( Vi−1 Vi )γ + ( ρi ρi−1 )γ . (6) From volumes Vi−1, Vi , and specific internal energy ei−1 of the previous and current steps, the current specific internal energy ei is computed. With Eq. (5), this leads to the pressure and ultimately to the normal force acting on the airbag fabric. To define the initial geometry of an airbag chamber, there are methods of scaling and more complex folding [51, 52]. Here, scaled initial meshes and the so-called initial metric methodology (IMM) [53] are used to model the expanding control volume. This method uses two meshes, an initial mesh and an undeformed reference or design mesh, see Fig. 15. Due to scaling, the initial mesh has highly distorted elements. An algorithm calculates in- ternal forces based on the geometric differences between the initial and the reference mesh and allows the elements to relax to their reference surface area. It circumvents considerably more complex folding procedures, which is usually a laborious and time-consuming task, especially in tight packaging, and detrimental in the initial stages of product development. Using these methods, the design variables for deployment, apart from the fabric’s imple- mented material model, include the airbag’s external and internal geometry, a mass inflow function with constant inflow temperature, and the exhaust hole area. Resulting inflation variables in a frontal impact are given in Fig. 16. Again, the operation can be divided into multiple consecutive phases: after firing (i) inflation because of gas generator input, followed by (ii) retention of airbag volume and pressure until (iii) rider impact and controlled rider deceleration, where the gas deflates from the chamber through the exhausts. Efficient simulation strategy to design a safer motorcycle 303 Fig. 15 Front airbag meshes for initial metric method (a)–(b) and inflated mesh (c) Fig. 16 Gas generator and airbag chamber variables of front airbag in a frontal impact 2.3.3 Leg impact protection modeling The 15 mm foam impact protection, shown in Fig. 17a, is modeled using eight-node (hexa- hedron) and six-node (prism) 3D solid elements (ELFORM = −2; fully integrated elective reduced solids). It provides a soft surface while achieving good impact energy absorption. The designated material, an energy absorbing high-density foam (425 g/l), has an impact- rate-dependent material behavior: on fast impacts, it responds very firm, while under slower loading, it remains soft. Ls-Dyna provides a variety of material models. To capture rate-dependent loading and hysteretic unloading of foams the material model *MAT_FU_CHANG _FOAM [54] is well suited. It is used for a wide variety of foam materials [55] and can be mapped with near test data. *MAT_FU_CHANG_FOAM requires as input data engineering stress vs. engineering strain as a function of multiple constant strain rates. The material model is characterized by a series of pendulum impact and quasi-static tests, which the authors conducted and elabo- rated in [35]. As summarized in Fig. 17b, for rate-dependent loading, the stress is approxi- 304 S. Maier, J. Fehr Fig. 17 Leg protector modeling mated by a surface interpolation of the experimental test with different impact speeds v1; for rate-dependent unloading, an approximation of the quasi-static loading vconst. and a param- eterized damage formulation from LS-Dyna is used. The approximation is a bipolynomial regression by minimizing a quadratic objective function, with linear equality and inequality constraints incorporating LS-Dyna material card requirements and knowledge about surface topology. 2.4 Rider surrogates Since the workflow includes the two important simulation software environments of MBS and FE, the choice of numerical models for representing motorcyclists is broad. Thus, MB and FE models of ATDs as artificial mechanical human surrogates and virtual only human body models (HBMs) with a wide range of complexity and computational costs can be used. These are instrumented with sensors that measure accelerations, forces, moments, and de- formations at multiple body locations. The selection for the results of this work consists of the models from Fig. 18 (excluding sex variants). The positioning of kinematically cou- pled MB models is particularly simple, fast, and, above all, user-friendly. The positioning of FE models is far more complex and time-consuming, i.e., laborious in pre-processing. For ATD models, overlapping parts and elements can occur. Here, this occurs at the thighs- to-pelvis connections of the molded hip flexion angle, the knee bends, where the shanks touch the thighs, and the ATD wrists. HBMs have complex joints and motion patterns with nonidealized degrees of freedom surrounded by deformable tissue. To position the models, transient seating simulation incorporates deformation and compliance between deformable parts. For a detailed description of the authors’ work on the positioning and contact initial- ization process of ATDs, see [33]; for HBMs, see [34]. To evaluate rider surrogate loading, injury criteria correlate sensor loads with probabil- ities of certain injuries for specific body regions. Since the severity of an injury cannot always be correlated with the maxima of the load, derived quantities, such as the head in- jury criterion (HIC) [60], have been developed. It is a normalized integral of the resultant of Efficient simulation strategy to design a safer motorcycle 305 Fig. 18 Models ((a): Madymo [56], (b): LSTC [57], (c): Project VIRTUAL [58, 59]) of anthropometric test devices and human body models as rider surrogates (partially blanked). the acceleration ares(t) = ‖a(t)‖2 measured at the head’s center of gravity, and it rates the loading by its value and duration by HIC(t2 − t1) = max t1,t2 { (t2 − t1) [ 1 t2 − t1 ∫ t2 t1 ares(t)dt ]2.5 } , (7) where ares(t) is expressed in standard gravitational acceleration g = 9.81 m/s2 and t in s. Furthermore, a criterion that represents the maximum value over a certain period of time is often used, such as the a3ms-criterion of the resultant acceleration a3ms = max t1 ( min t1≤t≤t1+3ms ares(t) ) . (8) The selection of injury criteria with recommended biomechanical thresholds considered are summarized in Table 1. It is based on a comprehensive set of injury criteria and corre- sponding biomechanical limits for motorcyclists from a literature review by [38]. 3 Simulation results In the following, results from stages 1 to 3 for a frontal collision against a stationary passen- ger car equivalent to the full-scale crash test of the conventional motorcycle, see Fig. 19, are given. 3.1 Stage 1: combined multibody and finite element model This stage is used to predict the impact dynamics of the vehicles and the rider’s interaction with the passive safety system, see Fig. 20. The frontal shape of the motorcycle with its vo- luminous cockpit fairing and overhanging nose act as a crash structure that allows for some deformation. This results in a much larger front fork deformation. This high contact point counteracts a rotational moment that would otherwise cause the motorcycle to rotate around 306 S. Maier, J. Fehr Table 1 Evaluated injury criteria with biomechanical limits from [38] for motorcyclists body region injury criterion abbreviation limit head head injury criterion HIC(36) 1000 for 36 ms resultant head acceleration a3ms 80 g over 3 ms neck neck tensile force Fz,1ms 3.3 kN over 1 ms Fz,45ms 1.1 kN over 1 ms neck compression force F−z,1ms 4 kN over 1 ms F−z,45ms 1.1 kN over 1 ms neck shear force Fxy,1ms 4 kN over 1 ms Fxy,45ms 1.1 kN over 1 ms neck rearward moment My−,max 190 Nm neck forward moment My,max 57 Nm thorax resultant thorax acceleration a3ms 80 g over 3 ms chest deflection deflmax 75 mm pelvis resultant pelvis acceleration a3ms 80 g over 3 ms femur femur axial force |Fz|max 10 kN Fig. 19 Frontal impact scenario according to ISO 13232 [61] its transverse/pitch axis. This leads to a stiffer impact response as shown by the faster ve- locity drop in Fig. 21. Also, because the accident trajectory of the motorcycle with safety concept and its rider results in less rotation and less upward motion, the energy transfer to the car is greater. As shown, the belts restrain the pelvis. Belt load limiting allows for some forward displacement of the rider. The failing attachment of the windshield allows for max- imum forward displacement of the rider. The restraint leads to an upper body rotation with the head and upper torso, consequently impacting the frontal airbag. The pressure regula- tion of the airbags through exhaust holes prevents a rebound of the rider. The overlay of head linear acceleration shows that a maximum is reached during deceleration from the belt restraint, not when impacting the airbag. Figure 21 shows the velocities of the vehicles and the main body parts of the rider for the conventional motorcycle impact (a) and the motorcycle with proposed rider restraint (b). In the conventional impact, the rider is decelerated abruptly by the pelvis impacting the tank and the helmeted head impacting the car—the restrained rider is decelerated more con- tinuously over a more extended time period. Simulations of the novel motorcycle concept without its safety equipment would not be an adequate benchmark for evaluating its ben- efits. For conventional motorcycles, it is considered advantageous that a motorcyclist does not become entangled in parts of the motorcycle body and instead separates from the mo- torcycle as soon as possible [38]. The proposed body is not designed to allow such accident kinematics. Efficient simulation strategy to design a safer motorcycle 307 Fig. 20 Impact trajectories (black) with head accelerations (red) in combined MB/FE simulations with Hy- brid III 50th ATD (Note: The right airbags are not displayed) In Fig. 22, the resulting injury criteria are plotted relative to their biomechanical limits for the real-world crash test versus the MB simulation of the restraint safety concept. In conventional accidents, the head hits the roof rail. This leads to high head accelerations, a critical compression of the cervical spine with a very high peak rearward neck moment. For the impact with safety concept, the highest sensed body loads are the head, thorax, and pelvis accelerations and neck tensile forces. These loads depend on the implemented belt load limit selected based on a trade-off between tolerable body loads and acceptable forward body displacement. Compared to the conventional crash impact body loads, most criteria are reduced, lower overall, and within their respective limits. 3.2 Stage 2: prescribed finite element rider interaction model The modeling of the impact dynamics in the FE environment is used to investigate more rider surrogate models while focusing on the interaction behavior with the safety systems. ATDs are not applicable omnidirectional. Instead, they are specific for different load directions and postures and are most often designed to depict occupants of automobiles. The stage 2 simu- lation environment is well suited to study state-of-the-art HBMs. Figure 23 shows a compar- ison of the kinematic impact response of a 5th Hybrid III ATD with a female 50th VIVA 50F HBM. Overall, the impact kinematics turn out very similar. With the HBM, it can be ob- served that the belt slips down at 225 ms. This could possibly be prevented with other belt variants. In evaluating biomechanical injuries, the capabilities of the HBMs go far beyond 308 S. Maier, J. Fehr Fig. 21 Vehicle and rider impact velocities in MB and combined MB/FE simulations Fig. 22 Injury criteria relative to the respective biomechanical limit for the conventional motorcycle (top; experimental data of SH01.01) and the motorcycle with rider restraint (bottom; stage 1 simulation) those of mechanical surrogate models and their virtual representation. Here, HBMs provide much more profound insight. They, e.g., allow observing strain-based human injury mecha- nisms. The energy balance of each part or a subset of parts in LS-Dyna is the sum Etotal(t) = Ekinetic(t) + Einternal(t) + Epotential(t) + Econtact(t) + Esliding interface(t) + Ehourglass(t). (9) Figure 24 gives the energy balance of VIVA 50M and 50F HBMs during the impact. The kinetic, internal, and potential energies represent the rider’s current energy level. The contact energy is the energy transferred to external contact partners. The sliding interface energy is the stored energy of internal contacts, and the hourglass energy is an FEM-specific numerical Efficient simulation strategy to design a safer motorcycle 309 Fig. 23 Impact responses of Hybrid 5th ATD (top) and 50F ViVA+ HBM (bottom) in coupled FE/MB simu- lations (Note: The right airbags are not displayed) Fig. 24 Rider energy balances in coupled FE/MB simulations quantity and should be minimal or at least very low. Figure 24a shows that initial energy is solely its initial kinetic energy, which is reduced to a minimum during impact. Most of it is transferred to the contact partners; some of it is converted to internal energy. The rider is decelerated through the application of forces. Anything that applies these forces can be considered part of the restraint system. Examining the contributors to the contact energy transfer in Fig. 24b, the cumulative belt contact energy transfer in the final state at 300 ms has the most significant share, followed by the airbags and the remaining cockpit surfaces. The leg impact protection has a comparably small share. 3.3 Stage 3: full finite element model The full FE representation is used to provide (i) the deformation characteristics and vehicle strength of the motorcycle itself and (ii) the structural interaction between the motorcycle 310 S. Maier, J. Fehr Fig. 25 Frontal deformation characteristics of motorcycle (top) with stacked internal energy absorption of main structural motorcycle components (bottom) from full FE simulation of rigid wall impact Fig. 26 Motorcycle to car structural interaction and intrusion behavior in full FE simulation (Note: The right airbags are not displayed) and the opposing vehicle. In Fig. 25 the distance d after an impact with an initial veloc- ity vMoto = 13.4 m/s is given. The plot of distance d against the resulting rigid wall force shows several phases of impact deformation, which can be split into compression of the front tire, followed by the collapse of the front wheel rim, and then the collapse of the front fork. The evaluation of the internal energy absorption (initial dissipation and comparably low elastic release) of the main structural components involved shows that the impact en- ergy during the subsequent deformation phase until maximum compression is dissipated by the crash box, main frame, and fairing of the motorcycle. Since these all act at a high point of contact, this counteracts pitching of the motorcycle. The structural interaction and intrusion behavior between the motorcycle and the opposing vehicle is shown through a sec- tion cut of the opposing car in Fig. 26. The front wheel of the motorcycle collides with the car’s passenger side front door sill. The motorcycle cockpit deforms the car’s door inwards, while the motorcycle’s cockpit itself is little deformed. The deformation of the telescopic Efficient simulation strategy to design a safer motorcycle 311 Table 2 Solver run times for an impact scenario with 300 ms (AMD Ryzen 9 5950X 16-Core CPU@3.4GHz) stage 1: stage 2: stage 3: combined MB/FE model coupled FE/MB model full FE model w/ Hyb III 50th ATD: w/ Hyb III 50th ATD: w/ Hyb III 50th ATD: 36 min/16 CPUs (SMP) 10 h/16 CPUs (MPP) 59.3 h/16 CPUs (MPP) w/ VIVA 50M HBM: 23.3 h/16 CPUs (MPP) SMP: shared memory processing MPP: massively parallel processing front fork is, by design, much greater than on the conventional motorcycle discussed here. This analysis shows that the very rigid sill structure promotes motorcycle pitching. Never- theless, the motorcycle does not roll over. The maximum pitching angle in this simulation is β(201 ms) = 15.1◦. 4 Discussion Models with varying levels of detail and complexity have been developed. The models have modular input decks for the restraint systems, rider surrogate models, and opposing (vehicle) structures. Custom translators are used for the exchange between the software environments, which transfer, e.g., scenario parameters and FE-meshes (particularly airbag geometries) from design iterations. Each of the model and simulation environments allows for specific investigations with different simulation runtimes, see Table 2. They have individual advan- tages and are suitable for assessing and optimizing various aspects of the accident behavior of a motorized two-wheeler. The first stage MB approach features low complexity and low numerical costs while capturing the essential physics of the collision. This represents a nu- merically efficient way to tune and improve the safety system, i.e., adapting the properties, shapes, and locations of the safety components. The third stage full FE approach aims to fully represent the interaction of the collision to accurately predict the performance of the finalized design by replicating every structural component of the vehicles in great detail. However, this is bought by a challenging model generation and significantly increased com- putational efforts. The advantages of the second stage partial FE approach using MB-vehicle interaction of the numerical strategy presented here are the successive methodological model generation and gradually increased level of detail and expected fidelity while significantly reducing computation time compared to a full FE model representation. This hybrid model is comparable to already realized approaches for occupants in cars, see, e.g., [19]. This of- fers the possibility to consider a larger variance of accident scenarios or occupant diversity or to enable very complex and numerically expensive investigations with FE human body models at reduced numerical complexity. To exploit the full range of existing mathematical modeling in the area of crash simulation, reduced order models and response-surface mod- els must also be used to optimize aspects with numerical methods at reduced computational costs. Predicting the accident outcome of such a complex interaction is a complex task. At this stage, the MB model has shown that it can predict the vehicle interaction between a con- ventional motorcycle and a passenger car. Real-world performance of the safe motorcycle 312 S. Maier, J. Fehr has not yet been demonstrated, and the models are not currently validated by experiments. Nevertheless, the work shows what is possible with a solely simulation-driven approach. This will be necessary in the future to quickly develop innovative solutions with minimal funding to investigate a new approach conceptually quickly. 5 Conclusions RQ1: In a frontal collision scenario against a stationary opponent, in this case a passenger car, the simulations of the restraint system show its effect of a guided and controlled trajec- tory and deceleration of the motorcycle rider. It consists of (i) restraining the rider’s pelvis to the motorcycle, (ii) decelerating the resulting upper and lower body motion by airbags and leg impact protectors, and (iii) preventing hard contact with opposing structures. This re- sults in fewer critical biomechanical loads on the rider. There may be other load cases where detrimental effects result from the rider restraint. To thoroughly answer the question, in- vestigations of additional scenarios that include behavior in subsequent secondary accident phases after an initial impact, including solo accidents, are necessary. RQ2: The shown numerical research strategy outlines a novel procedure in virtual motor- cycle accident research and passive safety equipment development. It is a systematic model generation approach with different levels of computational effort and model complexity. It demonstrates a meaningful combination of current state-of-the-art MB and FE simulation software environments to model the same problem with method-specific advantages and disadvantages. The strategy enabled the virtual design and dimensioning of a novel safety system with little time and resources and aimed for a step-by-step validation of individual components in the future. Although the procedure is presently applied for a motorcycle with a novel restraint safety concept, applying similar procedures in virtual research for conven- tional (powered) two-wheelers is highly desirable. Acknowledgements We thank the consortium of the Safe Motorcycle for the collaboration and fruit- ful discussions, especially Laurent Doléac, ideator and inventor of the safety system. We thank the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) for supporting this work by fund- ing: EXC2075 – 390740016 under Germany’s Excellence Strategy. We acknowledge the support by the Stuttgart Center for Simulation Science (SimTech). Author contributions S.M. developed the shown modelling and simulation framework, modelled the mo- torcycle and passive safety systems, simulated and analysed the simulation data, and prepared the text and figures of the contribution. J.F advised S.M. on methods and approaches in computational multi-body and finite element methods and passive safety and reviewed the manuscript. S.M and J.F revised the manuscript together. Funding Note Open Access funding enabled and organized by Projekt DEAL. Declarations Competing interests The authors declare no competing interests. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. http://creativecommons.org/licenses/by/4.0/ Efficient simulation strategy to design a safer motorcycle 313 References 1. Berg, F.A., Bürkle, H., Schmidts, F.: Analysis of the passive safety of motorcycles using accident investi- gations and crash tests. In: Proceedings of the 16th International Technical Conference on the Enhanced Safety of Vehicles (ESV), 98-S10-O-11, Winsor Ontario, Canada, pp. 2221–2236 (1998). https://www- esv.nhtsa.dot.gov/Proceedings/0/16th International Technical Conference on the Enhanced Safety of Vehicles Proceedings%20Volume 3 of 3.pdf#page=587 2. Grandel, J., Schapter, D.: Investigation into motorcycle, driver and passenger safety in motorcycle ac- cidents with two motorcycle riders. In: Proceedings of the 11th International Technical Conference on the Experimental Safety of Vehicles (ESV), Washington, pp. 888–900 (1987). https://www-esv.nhtsa. dot.gov/Proceedings/0/Eleventh International Technical Conference on Experimental Safety of Vehicles United States.pdf#page=905 3. Lindenmann, M., Grandel, J., Berg, F.: Collision dynamics in experimental simulations of 90◦ mo- torcycle collisions against the side of moving passenger cars. In: Proceedings of the IRCOBI Con- ference, Zurich, Switzerland, pp. 289–302 (1986). http://www.ircobi.org/wordpress/downloads/irc1986/ pdf_files/1986_21.pdf 4. ACEM: MAIDS – In-depth Investigations of Accidents Involving Powered Two Wheelers, Final Report 2.0. European Association of Motorcycle Manufacturers (ACEM), Brussels, Belgium (2009). https:// www.maids-study.eu/pdf/MAIDS2.pdf 5. ACEM: ACEM Statistical Release 2018, European Association of Motorcycle Manufactur- ers (ACEM) (2019). https://www.acem.eu/images/publiq/2019/ACEM-statistical-release---2018-figures. pdf. Accessed 23 June 2022 6. ACEM: ACEM Statistical Release 2019, European Association of Motorcycle Manufactur- ers (ACEM) (2020). https://www.acem.eu/images/publiq/2019/ACEM_statistical_release_Jan-Dec_019. pdf. Accessed 23 June 2022 7. ACEM: Motorcycle and Moped Registrations in Key European Markets Surpass 2019 Pre-Pandemic Levels, European Association of Motorcycle Manufacturers (ACEM) (2022). https://www.acem.eu/ images/publiq/2022/ACEM_-_Statistical_press_release_-_January_-_December_2021.pdf. Accessed 23 June 2022 8. WHO: WHO Global Status Report on Road Safety 2018. World Health Organization (WHO) Geneva, Switzerland (2018). https://apps.who.int/iris/bitstream/handle/10665/276462/9789241565684-eng.pdf 9. BASt: Verkehrs- und Unfalldaten, Kurzzusammenstellung der Entwicklung in Deutschland. Bunde- sanstalt für Straßenwesen (BASt), Bergisch Gladbach, Germany (2021) (in German). https://www.bast. de/DE/Publikationen/Medien/VU-Daten/VU.html. Accessed 17 May 2022 10. Statistisches Bundesamt: Verkehrsunfälle Verkehr 8(7), 2080700207004 (2021) (in German). https:// www.destatis.de/DE/Themen/Gesellschaft-Umwelt/Verkehrsunfaelle/Publikationen/Downloads- Verkehrsunfaelle/verkehrsunfaelle-jahr-2080700207004.pdf?__blob=publicationFile 11. Maier, S.: Sicheres Motorrad - Ein Projekt des Technologietransferprogramms (2020) (in Ger- man). https://www.youtube.com/watch?v=4NAfxNl9fE4. YouTube Kanal der e-mobil BW GmbH, Lan- desagentur für neue Mobilitätslösungen und Automotive Baden-Württemberg 12. Noorsumar, G., Rogovchenko, S., Robbersmyr, K.G., Vysochinskiy, D.: Mathematical models for as- sessment of vehicle crashworthiness: a review. Int. J. Crashworthiness 27(5), 1545–1559 (2021). https:// doi.org/10.1080/13588265.2021.1929760 13. Schiehlen, W., Eberhard, P.: Applied Dynamics, 1st edn. Springer, Heidelberg (2014). https://doi.org/10. 1007/978-3-319-07335-4 14. Carvalho, M., Ambrósio, J.: Identification of multibody vehicle models for crash analysis using an opti- mization methodology. Multibody Syst. Dyn. 24, 325–345 (2010). https://doi.org/10.1007/s11044-010- 9221-z 15. Shang, S., Masson, C., Llari, M., Py, M., Ferrand, Q., Arnoux, P.J., Simms, C.: The predictive capac- ity of the MADYMO ellipsoid pedestrian model for pedestrian ground contact kinematics and injury evaluation. Accid. Anal. Prev. 149, 105803 (2021). https://doi.org/10.1016/j.aap.2020.105803 16. Ambrósio, J., Carvalho, M., Milho, J., Escalante, S., Martín, R.: A validated railway vehicle interior layout with multibody dummies and finite element seats models for crash analysis. Multibody Syst. Dyn. 54(2), 179–212 (2021). https://doi.org/10.1007/s11044-021-09794-w 17. Belytschko, T., Liu, W.K., Moran, B.: Nonlinear Finite Elements for Continua and Structures. Wiley, Chichester (2000) 18. Duddeck, F., Wehrle, E.: Recent advances on surrogate modeling for robustness assessment of struc- tures with respect to crashworthiness requirements. In: Proceedings of the 10th European LS-DYNA Users Conference, Würzburg, Germany, (2015). https://www.dynalook.com/conferences/10th-european- ls-dyna-conference/8 Optimization I -Robustness/01-Duddeck-TechnicalUnivMunich-P.pdf https://www-esv.nhtsa.dot.gov/Proceedings/0/16th International Technical Conference on the Enhanced Safety of Vehicles Proceedings%20Volume 3 of 3.pdf#page=587 https://www-esv.nhtsa.dot.gov/Proceedings/0/16th International Technical Conference on the Enhanced Safety of Vehicles Proceedings%20Volume 3 of 3.pdf#page=587 https://www-esv.nhtsa.dot.gov/Proceedings/0/16th International Technical Conference on the Enhanced Safety of Vehicles Proceedings%20Volume 3 of 3.pdf#page=587 https://www-esv.nhtsa.dot.gov/Proceedings/0/Eleventh International Technical Conference on Experimental Safety of Vehicles United States.pdf#page=905 https://www-esv.nhtsa.dot.gov/Proceedings/0/Eleventh International Technical Conference on Experimental Safety of Vehicles United States.pdf#page=905 https://www-esv.nhtsa.dot.gov/Proceedings/0/Eleventh International Technical Conference on Experimental Safety of Vehicles United States.pdf#page=905 http://www.ircobi.org/wordpress/downloads/irc1986/pdf_files/1986_21.pdf http://www.ircobi.org/wordpress/downloads/irc1986/pdf_files/1986_21.pdf https://www.maids-study.eu/pdf/MAIDS2.pdf https://www.maids-study.eu/pdf/MAIDS2.pdf https://www.acem.eu/images/publiq/2019/ACEM-statistical-release---2018-figures.pdf https://www.acem.eu/images/publiq/2019/ACEM-statistical-release---2018-figures.pdf https://www.acem.eu/images/publiq/2019/ACEM_statistical_release_Jan-Dec_019.pdf https://www.acem.eu/images/publiq/2019/ACEM_statistical_release_Jan-Dec_019.pdf https://www.acem.eu/images/publiq/2022/ACEM_-_Statistical_press_release_-_January_-_December_2021.pdf https://www.acem.eu/images/publiq/2022/ACEM_-_Statistical_press_release_-_January_-_December_2021.pdf https://apps.who.int/iris/bitstream/handle/10665/276462/9789241565684-eng.pdf https://www.bast.de/DE/Publikationen/Medien/VU-Daten/VU.html https://www.bast.de/DE/Publikationen/Medien/VU-Daten/VU.html https://www.destatis.de/DE/Themen/Gesellschaft-Umwelt/Verkehrsunfaelle/Publikationen/Downloads-Verkehrsunfaelle/verkehrsunfaelle-jahr-2080700207004.pdf?__blob=publicationFile https://www.destatis.de/DE/Themen/Gesellschaft-Umwelt/Verkehrsunfaelle/Publikationen/Downloads-Verkehrsunfaelle/verkehrsunfaelle-jahr-2080700207004.pdf?__blob=publicationFile https://www.destatis.de/DE/Themen/Gesellschaft-Umwelt/Verkehrsunfaelle/Publikationen/Downloads-Verkehrsunfaelle/verkehrsunfaelle-jahr-2080700207004.pdf?__blob=publicationFile https://www.youtube.com/watch?v=4NAfxNl9fE4 https://doi.org/10.1080/13588265.2021.1929760 https://doi.org/10.1080/13588265.2021.1929760 https://doi.org/10.1007/978-3-319-07335-4 https://doi.org/10.1007/978-3-319-07335-4 https://doi.org/10.1007/s11044-010-9221-z https://doi.org/10.1007/s11044-010-9221-z https://doi.org/10.1016/j.aap.2020.105803 https://doi.org/10.1007/s11044-021-09794-w https://www.dynalook.com/conferences/10th-european-ls-dyna-conference/8 Optimization I - Robustness/01-Duddeck-TechnicalUnivMunich-P.pdf https://www.dynalook.com/conferences/10th-european-ls-dyna-conference/8 Optimization I - Robustness/01-Duddeck-TechnicalUnivMunich-P.pdf 314 S. Maier, J. Fehr 19. Reichert, R., Kan, C.D., Marzougui, D., Mahadevaiah, U., Morgan, R.G., Park, C.K., Tahan, F.: Methodologies and examples for efficient short and long duration integrated occupant-vehicle crash sim- ulation. In: Proceedings of the 13th International LS-DYNA Users Conference, Detroit (2014). https:// www.dynalook.com/conferences/13th-international-ls-dyna-conference/automotive/methodologies- and-examples-for-efficient-short-and-long-duration-integrated-occupant-vehicle-crash-simulation 20. Happian-Smith, J., Macaulay, M., Chinn, B.: Motorcycle impact simulation and practical verification. In: Proceedings of the 11th International Technical Conference on the Experimental Safety of Vehicles (ESV), Washington, pp. 858–865 (1987) 21. Canaple, B., Rungen, G.P., Markiewicz, E., Drazetic, P., Happian-Smith, J., Chinn, B.P., Cesari, D.: Im- pact model development for the reconstruction of current motorcycle accidents. Int. J. Crashworthiness 7(3), 307–320 (2002). https://www.tandfonline.com/doi/abs/10.1533/cras.2002.0221 22. Sudyoddee, H., Behr, M., Llari, M., Koetniyom, S., Carmai, J.: Investigation of motorcyclist and pillion passenger injuries using numerical simulations. IOP Conf. Ser., Mater. Sci. Eng. 501, 012009 (2019). https://doi.org/10.1088/1757-899x/501/1/012009 23. Berg, A.F., Rücker, P., Gärtner, M., König, J., Grzebieta, R., Zou, R.: Motorcycle impacts into roadside barriers-real world accident studies, crash tests, and simulations carried out in Germany and Australia. In: Proceedings of the 19th International Technical Conference on the Enhanced Safety of Vehicles (ESV), 05-0095, Washington (2005). https://www-esv.nhtsa.dot.gov/Proceedings/19/05-0095-O.pdf 24. Ibitoye, A.B., Hamouda, A.M.S., Wong, S.V., Umar, R.S.R.: Simulation of motorcycle crashes with W- beam guardrail: injury patterns and analysis. Proc. Inst. Mech. Eng., H J. Eng. Med. 223(8), 1033–1040 (2009). https://doi.org/10.1243/09544119JEIM560 25. Moradi, R., Lankarani, H.M.: A Multi-Body Modeling and Design of Experiment Investigation of a Motorcyclist Impact on Roadside Barriers at Upright and Sliding Configurations, pp. 697–705 (2011). https://doi.org/10.1115/IMECE2011-62537 26. Lechner, M.: Der Zweirad-Alleinunfall: Relation zwischen realem Unfall, Versuch u. mathematischer Simulation (in German). Dissertation, Technische Universität München (1986) 27. Mongiardini, M., Walton, B., Grzebieta, R.H., McKay, M., Menictas, C., Berg, A., Rücker, P.: Development of a motorcycle FE model for simulating impacts into roadside safety barri- ers. In: Proceedings of the First International Roadside Safety Conference, E-C220, San Fran- cisco, pp. 657–673 (2017). http://www.tars.unsw.edu.au/news/2017/Mongiardini_et_al-Motorcycle_FE_ Model-TRB_eCircular-C220_First_International_Roadside_Safety_Conference.pdf 28. Schulz, N., Silvestri Dobrovolny, C., Hurlebaus, S.: Development of a finite element model of a motorcycle. In: Proceedings of the 14th International LS-DYNA Users Conference, De- troit (2016). https://www.dynalook.com/conferences/14th-international-ls-dyna-conference/modeling/ development-of-a-finite-element-model-of-a-motorcycle 29. Namiki, H., Nakamura, T., Iijima, S.: A computer simulation for motorcycle rider injury evaluation in collision. In: Proceedings of the 19th International Technical Conference on the Enhanced Safety of Vehicles (ESV), 05-0309, Washington (2005). https://www-esv.nhtsa.dot.gov/Proceedings/19/05-0309- W.pdf 30. Barbani, D., Baldanzini, N., Pierini, M.: Development and validation of an FE model for motorcycle–car crash test simulations. Int. J. Crashworthiness 19(3), 244–263 (2014). https://doi.org/10.1080/13588265. 2013.874672 31. Franz, U., Lorenz, B., Remfrey, J., Schöneburg, R.: Integrale Sicherheit von Kraftfahrzeugen: Biomechanik - Simulation - Sicherheit im Entwicklungsprozess, 4th edn. Springer, Wiesbaden (2013). https://doi.org/10.1007/978-3-8348-2608-4 (in German) 32. Maier, S., Doléac, L., Hertneck, H., Stahlschmidt, S., Fehr, J.: Evaluation of a novel passive safety concept for motorcycles with combined multi-body and finite element simulations. In: Proceedings of the IRCOBI Conference, IRC-20-38, Munich, Germany, pp. 250–265 (2020). https://doi.org/10.18419/ opus-11899 33. Maier, S., Doléac, L., Hertneck, H., Stahlschmidt, S., Fehr, J.: Finite element simulations of motorcy- clist interaction with a novel passive safety concept for motorcycles. In: Proceedings of the IRCOBI Conference, IRC-21-17, Munich, Germany, pp. 60–78 (2021). https://doi.org/10.18419/opus-11900 34. Maier, S., Kempter, F., Kronwitter, S., Fehr, J.: Positioning and simulation of human body models on a motorcycle with a novel restraint system. In: Proceedings of the IRCOBI Conference, IRC-22-22, Porto, Portugal, pp. 82–96 (2022). https://doi.org/10.18419/opus-12377 35. Maier, S., Helbig, M., Hertneck, H., Fehr, J.: Characterisation of an energy absorbing foam for motor- cycle rider protection in LS-DYNA. In: Proceedings of the 13th European LS-DYNA Users Conference, Ulm, Germany, (2021). https://doi.org/10.18419/opus-12000 36. Portal, R.J.F., Dias, J.M.P.: Multibody models for vehicle accident reconstruction. In: III European Con- ference on Computational Mechanics. Springer Netherlands, Lisbon, Portugal (2006). https://doi.org/10. 1007/1-4020-5370-3_775 https://www.dynalook.com/conferences/13th-international-ls-dyna-conference/automotive/methodologies-and-examples-for-efficient-short-and-long-duration-integrated-occupant-vehicle-crash-simulation https://www.dynalook.com/conferences/13th-international-ls-dyna-conference/automotive/methodologies-and-examples-for-efficient-short-and-long-duration-integrated-occupant-vehicle-crash-simulation https://www.dynalook.com/conferences/13th-international-ls-dyna-conference/automotive/methodologies-and-examples-for-efficient-short-and-long-duration-integrated-occupant-vehicle-crash-simulation https://www.tandfonline.com/doi/abs/10.1533/cras.2002.0221 https://doi.org/10.1088/1757-899x/501/1/012009 https://www-esv.nhtsa.dot.gov/Proceedings/19/05-0095-O.pdf https://doi.org/10.1243/09544119JEIM560 https://doi.org/10.1115/IMECE2011-62537 http://www.tars.unsw.edu.au/news/2017/Mongiardini_et_al-Motorcycle_FE_Model-TRB_eCircular-C220_First_International_Roadside_Safety_Conference.pdf http://www.tars.unsw.edu.au/news/2017/Mongiardini_et_al-Motorcycle_FE_Model-TRB_eCircular-C220_First_International_Roadside_Safety_Conference.pdf https://www.dynalook.com/conferences/14th-international-ls-dyna-conference/modeling/development-of-a-finite-element-model-of-a-motorcycle https://www.dynalook.com/conferences/14th-international-ls-dyna-conference/modeling/development-of-a-finite-element-model-of-a-motorcycle https://www-esv.nhtsa.dot.gov/Proceedings/19/05-0309-W.pdf https://www-esv.nhtsa.dot.gov/Proceedings/19/05-0309-W.pdf https://doi.org/10.1080/13588265.2013.874672 https://doi.org/10.1080/13588265.2013.874672 https://doi.org/10.1007/978-3-8348-2608-4 https://doi.org/10.18419/opus-11899 https://doi.org/10.18419/opus-11899 https://doi.org/10.18419/opus-11900 https://doi.org/10.18419/opus-12377 https://doi.org/10.18419/opus-12000 https://doi.org/10.1007/1-4020-5370-3_775 https://doi.org/10.1007/1-4020-5370-3_775 Efficient simulation strategy to design a safer motorcycle 315 37. NHTSA: Update to Future Midsize Lightweight Vehicle Findings in Response to Manufacturer Review and IIHS Small-overlap Testing. Report No. DOT HS 812 237, National Highway Transport Safety As- sociation (NHTSA), Washington (2016). https://www.nhtsa.gov/staticfiles/rulemaking/pdf/cafe/812237_ LightWeightVehicleReport.pdf 38. Berg, F.A., Rücker, P., Bürkle, H., Mattern, R., Kallieris, D.: Prüfverfahren für die passive Sicherheit mo- torisierter Zweiräder: Bericht zum Forschungsprojekt FE 82.121/1997. Berichte der Bundesanstalt für Straßenwesen. F, Fahrzeugtechnik, vol. 49. Wirtschaftsverlag NW, Verlag für neue Wissenschaft, Bre- merhaven (2004) (in German). https://bast.opus.hbz-nrw.de/opus45-bast/frontdoor/deliver/index/docId/ 221/file/F49.pdf 39. SAE J211-1:1995: Instrumentation for Impact Test – Part 1 – Electronic Instrumentation. SAE Interna- tional, Warrendale (1995). https://doi.org/10.4271/J211/1_199503 40. Gehre, C., Gades, H., Wernicke, P.: Objective rating of signals using test and simulation responses. In: Proceedings of the 21st International Technical Conference on the Enhanced Safety of Vehicles (ESV), 09-0407, Stuttgart, Germany (2009). https://www-esv.nhtsa.dot.gov/Proceedings/21/09-0407.pdf 41. Thunert, C.: CORAplus Release 4.0.4 User’s Manual (2017). https://www.pdb-org.com/en/ 42. NHTSA: Crash Simulation Vehicle Models. National Highway Transport Safety Association (NHTSA), Washington (2021). https://www.nhtsa.gov/crash-simulation-vehicle-models. Accessed 27 June 2022 43. Engel, A.: Airbag für motorisierte Zweiräder: Bericht zum Projekt 1.8903 der fahrzeugtechnischen Forschung des Bundesministers für Verkehr. Forschungsberichte des Bundesministers für Verkehr, Bere- ich Fahrzeugtechnik, vol. 8. Bundesanstalt für Straßenwesen, Bergisch Gladbach (1992) 44. Kuroe, T., Namiki, H., Iijima, S.: Exploratory study of an airbag concept for a large touring motorcy- cle: further research second report. In: Proceedings of the 19th International Technical Conference on the Enhanced Safety of Vehicles (ESV), 05-0316, Washington (2005). https://www-esv.nhtsa.dot.gov/ Proceedings/19/05-0316-O.pdf 45. Honda: 2006 Honda Gold Wing – Features (2005). https://hondanews.com/en-US/powersports/releases/ release-58200436c0998daa0c581e004c34c130-2006-honda-gold-wing-features. Accessed 20 December 2022 46. Christ, A., Büttner, J.: Detailed passenger airbag modelling for early stage events. In: Proceedings of the 11th German LS-DYNA Forum, Ulm, Germany, (2012). https://www.dynamore.de/dynamore/de/ download/papers/ls-dyna-forum-2012/documents/passive-3-1 47. Hirth, A., Haufe, A., Olovsson, L.: Airbag simulation with LS-DYNA past – present – future. In: Proceedings of the 6th German LS-DYNA Forum, Frankenthal, Germany, pp. 57–74 (2007). https://www.dynamore.de/en/downloads/papers/dynamore/en/downloads/papers/07- forum/forum07/keynote/airbag-simulation-with-ls-dyna-past-2013-present 48. Wang, J.T., Nefske, D.J.: A new CAL3D airbag inflation model. In: Proceedings of the International Congress and Exposition, Detroit (1988). https://doi.org/10.4271/880654. SAE Technical Paper Series, 880654 49. Fokin, D., Lokhande, N., Fredriksson, L.: On airbag simulation in LS-DYNA with the use of the arbi- trary Lagrangian-Eulerian method. In: Proceedings of 4th European LS-DYNA Users Conference, Ulm, Germany, pp. 11–22 (2003). https://www.dynamore.de/de/download/papers/dynamore/en/downloads/ papers/03-conference/occupant/on-airbag-simulation-in-ls-2013dyna-with-the-use 50. Lin, C.H., Cheng, Y.P.: Evaluation of LS-DYNA corpuscular particle method – passenger airbag applications. In: Proceedings of the 15th International LS-DYNA Users Conference, Detroit (2018). https://www.dynalook.com/conferences/15th-international-ls-dyna-conference/occupant- protection/evaluation-of-ls-dyna-r-corpuscular-particle-method-2013-passenger-airbag-applications/@ @download/file/37-A_OccupantProtection_057.pdf 51. Chawla, A., Bhosale, P.V., Mukherjee, S.: Modeling of folding of passenger side airbag mesh. In: SIAT 2005. The Automotive Research Association of India, ??? (2005). https://doi.org/10.4271/2005-26-059. SAE Paper No. 2005-26-059 52. Lakshminarayan, V., Lasry, D.: Finite element simulation of driver folded air bag deployment. SAE Transact. 100, 1969–1977 (1991). https://doi.org/10.4271/912904 53. Tanavde, A.S., Khandelwal, H., Lasry, D., Ni, X., Haug, E., Schlosser, J., Balakrishnan, P.: Airbag mod- eling using initial metric methodology. In: Proceedings of the International Congress and Exposition, San Francisco (1995). https://doi.org/10.4271/950875. SAE Technical Paper Series, 950875 54. LSTC: LS-DYNA Keyword User’s Manual Volume II Material Models, LS-DYNA R11. Livermore Soft- ware Technology Corporation (LSTC), Livermore (2018) 55. Croop, B., Lobo, H.: Selecting material models for the simulation of foams in LS-DYNA. In: Proceedings of the 7th European LS-DYNA Conference, Salzburg, Austria, (2009) 56. SISS: Simcenter Madymo Model Manual, Version 2020.1. Siemens Industry Software and Ser- vices BV (SISS), Rijswijk, Netherlands (2020) https://www.nhtsa.gov/staticfiles/rulemaking/pdf/cafe/812237_LightWeightVehicleReport.pdf https://www.nhtsa.gov/staticfiles/rulemaking/pdf/cafe/812237_LightWeightVehicleReport.pdf https://bast.opus.hbz-nrw.de/opus45-bast/frontdoor/deliver/index/docId/221/file/F49.pdf https://bast.opus.hbz-nrw.de/opus45-bast/frontdoor/deliver/index/docId/221/file/F49.pdf https://doi.org/10.4271/J211/1_199503 https://www-esv.nhtsa.dot.gov/Proceedings/21/09-0407.pdf https://www.pdb-org.com/en/ https://www.nhtsa.gov/crash-simulation-vehicle-models https://www-esv.nhtsa.dot.gov/Proceedings/19/05-0316-O.pdf https://www-esv.nhtsa.dot.gov/Proceedings/19/05-0316-O.pdf https://hondanews.com/en-US/powersports/releases/release-58200436c0998daa0c581e004c34c130-2006-honda-gold-wing-features https://hondanews.com/en-US/powersports/releases/release-58200436c0998daa0c581e004c34c130-2006-honda-gold-wing-features https://www.dynamore.de/dynamore/de/download/papers/ls-dyna-forum-2012/documents/passive-3-1 https://www.dynamore.de/dynamore/de/download/papers/ls-dyna-forum-2012/documents/passive-3-1 https://www.dynamore.de/en/downloads/papers/dynamore/en/downloads/papers/07-forum/forum07/keynote/airbag-simulation-with-ls-dyna-past-2013-present https://www.dynamore.de/en/downloads/papers/dynamore/en/downloads/papers/07-forum/forum07/keynote/airbag-simulation-with-ls-dyna-past-2013-present https://doi.org/10.4271/880654 https://www.dynamore.de/de/download/papers/dynamore/en/downloads/papers/03-conference/occupant/on-airbag-simulation-in-ls-2013dyna-with-the-use https://www.dynamore.de/de/download/papers/dynamore/en/downloads/papers/03-conference/occupant/on-airbag-simulation-in-ls-2013dyna-with-the-use https://www.dynalook.com/conferences/15th-international-ls-dyna-conference/occupant-protection/evaluation-of-ls-dyna-r-corpuscular-particle-method-2013-passenger-airbag-applications/@@download/file/37-A_OccupantProtection_057.pdf https://www.dynalook.com/conferences/15th-international-ls-dyna-conference/occupant-protection/evaluation-of-ls-dyna-r-corpuscular-particle-method-2013-passenger-airbag-applications/@@download/file/37-A_OccupantProtection_057.pdf https://www.dynalook.com/conferences/15th-international-ls-dyna-conference/occupant-protection/evaluation-of-ls-dyna-r-corpuscular-particle-method-2013-passenger-airbag-applications/@@download/file/37-A_OccupantProtection_057.pdf https://doi.org/10.4271/2005-26-059 https://doi.org/10.4271/912904 https://doi.org/10.4271/950875 316 S. Maier, J. Fehr 57. LSTC: Detailed Hybrid III Frontal Impact Dummy Models, Livermore Software Technology Corpora- tion (LSTC) (2022). https://www.lstc.com/download/dummy_models 58. John, J.D., Klug, C., Kranjec, M., Svenning, E., Iraeus, J.: Hello, world! VIVA+: a human body model lineup to evaluate sex-differences in crash protection. Front. Bioeng. Biotechnol. 10, 918904 (2022). https://doi.org/10.3389/fbioe.2022.918904 59. OpenVT: VIRTUAL OpenVT Gitlab Platform Repository (2021). https://openvt.eu/fem/viva/vivaplus_ v0_2/. Accessed 9 September 2021 60. Kleinberger, M., Sun, E., Eppinger, R., Kuppa, S., Saul, R.: Development of Improved Injury Crite- ria for the Assessment of Advanced Automotive Restraint Systems, NHTSA Docket No. 98-4405-9. National Highway Transport Safety Association (NHTSA), Washington (1998). https://www.nhtsa.gov/ sites/nhtsa.gov/files/criteria.pdf 61. ISO 13232:2005: Motorcycles – Test and Analysis Procedures for Research Evaluation of Rider Crash Protective Devices Fitted to Motorcycles. International Organization for Standardization, Geneva, Switzerland (2005) Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. https://www.lstc.com/download/dummy_models https://doi.org/10.3389/fbioe.2022.918904 https://openvt.eu/fem/viva/vivaplus_v0_2/ https://openvt.eu/fem/viva/vivaplus_v0_2/ https://www.nhtsa.gov/sites/nhtsa.gov/files/criteria.pdf https://www.nhtsa.gov/sites/nhtsa.gov/files/criteria.pdf Efficient simulation strategy to design a safer motorcycle Abstract Introduction Novel safety concept for motorcycles Computational models for crashworthiness Modeling Modeling and simulation strategy Vehicles Multibody systems Coupled finite element and multibody model Full finite element models Rider restraint Thigh belt modeling Airbag modeling Leg impact protection modeling Rider surrogates Simulation results Stage1: combined multibody and finite element model Stage2: prescribed finite element rider interaction model Stage3: full finite element model Discussion Conclusions Acknowledgements References