Anti-Roll Bar Modeling for NVH and Vehicle
Total Page:16
File Type:pdf, Size:1020Kb
Anti-Roll Bar Model for NVH and Vehicle Dynamics Analyses Anti-Roll Bar Model for NVH and Vehicle Dynamics Analyses Tobolar,Anti-Roll Jakub and Bar Leitner, Modeling Martin and Heckmann, for NVH Andreas and Vehicle Dynamics Analyses 99 Jakub Tobolárˇ1 Martin Leitner1 Andreas Heckmann1 1German Aerospace Center (DLR), Institute of System Dynamics and Control, Wessling, [email protected] Abstract A The latest extension of the DLR FlexibleBodies Library concerns the field of automotive applications, namely the anti-roll bar. For the particular purposes of NVH and ve- hicle dynamics, the anti-roll bar module provides two ap- propriate levels of detail, both being based upon the beam preprocessor. In this paper, the procedure on preparing the models and their application for particular automotive related analyses is presented. Keywords: anti-roll bar, vehicle chassis, flexible body, beam model, finite element C S Figure 1. Vehicle axle with an anti-roll bar (color emphasized, 1 Introduction courtesy of Wikimedia Commons). Whenever an automotive suspension is excited in vertical direction due to road irregularities or driving maneuvers (Heckmann et al., 2006) provides capabilities to incorpo- in an asymmetrical way, i.e. differently on the right and rate data that originate from FE models in Modelica mod- the left side of the vehicle, the roll motion of the car body els. Thus, a tool chain to perform vehicle dynamics sim- is stimulated. This concerns – in common case – the com- ulation including the structural characteristics of anti-roll fort and driving experience of the car passengers. In limit bars is in principle available. conditions’ situations, such a roll motion can influence the In common design tasks, driving maneuvers or noise, road-holding forces in a way that vehicle’s driving safety vibration and harshness (NVH) scenarios are first ana- is affected significantly. Consequently, it is advantageous lyzed in multibody simulations. Then, the FE method is to introduce an additional suspension component in partic- used to redesign the anti-roll bar in order to improve its ular tailored to influence the dynamical roll motion char- characteristics. Subsequently, a FE to multibody interface acteristics independently from the layout of the vertical has to be used to prepare the new FE data for the DLR suspension. This so-called anti-roll bar (also called stabi- FlexibleBodies Library and, finally, the vehicle dynamics lizer or anti-sway bar), see e.g. (Rill, 2012) or (Heißing simulation has to be invoked again in order to assess the and Ersoy, 2011), connects the suspensions on the right modification. This tool chain or loop, respectively, is in- and the left side of the vehicle’s axle by a cranked bar that convenient and makes it difficult to set up computational acts as a torsional spring, see Figure1. optimization procedures. Therefore, the design of the anti-roll bar is mainly tar- The given background motivates the introduction of a geted on its torsional stiffness, but also has to comply with new modeling capability called AntiRollBar into the DLR the available space at the underfloor and must allow for at- FlexibleBodies Library. In the present paper, a principle of tachments to the vehicle body and to both vertical suspen- the flexible body modeling and of the beam theory behind sions. These requirements quite often result in the anti-roll the AntiRollBar model is given in Sections 2.1 and 2.2, bar to be a geometrical complex structural element that is respectively. In Section3, a framework of the AntiRollBar prone for dynamical vibrations. and its parametrization is discussed. Section4 presents In daily practice, the Finite Element (FE) method first simulation experiments provided. turned out to be the adequate tool to design the geometri- cal and the structural properties of anti-roll bars. However, the driving behavior of vehicles, to which the anti-roll bar 2 Theoretical Background significantly contributes, is commonly developed using 2.1 Flexible Bodies Theory multibody simulation – generally utilizing the MultiBody package of the Modelica Standard Library in the Modelica The mechanical description of flexible bodies in multi- community. In addition, the DLR FlexibleBodies Library body systems is based on the floating frame of reference DOI Proceedings of the 13th International Modelica Conference 99 10.3384/ecp1915799 March 4-6, 2019, Regensburg, Germany Anti-Roll Bar Model for NVH and Vehicle Dynamics Analyses where the following quantities and symbols appear: m body mass, I3 3 × 3 identity matrix, dCM(q) position of center of mass, J(q) inertia tensor, Ct (q) inertia coupling matrix (translational), Cr(q) inertia coupling matrix (rotational), hw (w;q;q˙) gyroscopic and centripetal forces, he external forces, Me structural mass matrix, Figure 2. Vector chain of the floating frame of reference. K e structural stiffness matrix, De structural damping matrix. 1 approach, i.e. the absolute position r =r(c;t) of a specific In the context of the anti-roll bar modeling, the struc- body particle is subdivided into three parts: tural mass, stiffness and damping masses are gained as the result of a FE preprocessing step whose background is • the position vector rR = rR(t) to the body’s reference given in the following section. Note, that the FE prepro- frame, cessing is implemented internally so that the user does not need to switch to a different modeling tool. • the initial position of the body particle within the body’s reference frame, i.e. the Lagrange coordinate c 6= c(t), 2.2 FE Beam Theory • and the elastic displacement u(c;t) that is approxi- mated by a Taylor expansion, here limited to first The structural models used in multibody analysis are order terms, with space-dependent mode shapes usually obtained from FE analysis and subsequently re- F(c) 2 3;n and time-dependent modal amplitudes R duced by e.g. modal decomposition approaches. For q(t) 2 n, cf. (Wallrapp, 1994): R the anti-roll bar structural models, a simple, classical fi- nite element beam formulation is employed. Therein, r = rR +c +u ; u = F q : (1) the three-dimensional problem is split into a two- dimensional, cross-sectional analysis and a subsequent, one-dimensional analysis along the beam’s reference axis. All terms in equation (1) are resolved w.r.t. the body’s Solving the two-dimensional problem simply involves floating frame of reference (R). That’s why the angular ve- integration of the material properties (Young’s modu- locity of the reference frame w R have to be taken into ac- lus, shear modulus and density) over the specified cross- count when the kinematic quantities velocity v and accel- sectional geometry. With the resulting cross-sectional eration aR of a particle are derived, see (Heckmann et al., stiffness and inertia resultants, the corresponding consti- 2006). The decomposition in equation (1) makes it possi- tutive matrix C can be built. Along with a given strain ble to superimpose a large nonlinear overall motion of the field, one arrives at the description for the force and mo- reference frame with small elastic deformations. ment distributions along the beam axis. The kinematic quantities are inserted into Jourdain’s principle of virtual power. Subsequently, the equations The one-dimensional analysis is based on (Bazoune of motion of an unconstrained flexible body are formu- et al., 2003), where an adjustable Timoshenko beam el- lated neglecting deflection terms of higher than first order ement was implemented, that uses linear shape functions (Wallrapp, 1994, (38)): for longitudinal displacements and torsional deformation. In order to describe the bending deformation, a cubic 2 32 3 ansatz function is used in the lateral displacements and mII sym: a 3 R corresponding rotational fields. The unknown coefficients md˜ J w˙ = 4 CM 54 R5 can be solved using the description for the total slopes in- C C M q¨ t r e cluding a constant transversal shear, the force and moment 2 3 0 equilibrium equations and the discrete boundary condi- = hw − 4 0 5 + he; (2) tions at both ends of the beam. These equations can then K e q +De q˙ be partitioned by discerning between displacement field contribution and discrete boundary condition excitation. 1Both vectors and matrices are written in bold symbols, whereby In a parametric space from [0,1], the shape functions that vectors are of lower case letters and matrices of upper case letters. are needed for the shear displacements and bending rota- 100 Proceedings of the 13th International Modelica Conference DOI March 4-6, 2019, Regensburg, Germany 10.3384/ecp1915799 Anti-Roll Bar Model for NVH and Vehicle Dynamics Analyses tions read: few degrees of freedom (DoF). Therein, all DoF are par- 1 titioned into an "analysis" set (index a) and an "omitted" N1 = (1 − 3x 2 + 2x 3); set (index o) and a constraint w.r.t. the omitted degrees of bs + 1 Q freedom being force free is applied. The resulting trans- 2 1 2 3 Nbs = (3x − 2x ); formation matrix between the analysis set and the original, 1 + Q full DoF set reads, 1 1 2 3 2 Nbb = (x − 2x + x + 1=2(2x − x )Q); 1 + Q xa I = −1 xa : (5) l xo −K K oa N2 = (−x 2 + x 3 + 1=2x 2Q); oa bb 1 + Q In order for the mass and stiffness matrix to be compat- Q N1 = (1 − x); ible with the generalized equations of motion, a normal ss + 1 Q modes analysis is performed afterwards. The resulting 2 Q Nss = x; eigenfrequencies and eigenvectors form a modal solution 1 + Q set that is used to reduce the structural degrees of free- 1 Ql dom to the user-specified number of frequencies and mode Nsb = − 1=2x; 1 + Q shapes retained.