Abstract
Computer simulations have become very popular in the automotive industry. In order to achieve a good conformity with field test, sophisticated vehicle models are needed. A real vehicle incorporates many complex dynamic systems, such as the drive train, the steering system and the wheel/axle suspension. On closer inspection some force elements such as shock absorbers and hydro-mounts turn out to be dynamic systems too. Modern vehicle models consist of different subsystems. Then, each subsystem may be modeled differently and can be tested independently. If some subsystems are available as a set of nested models of different complexity it will be even possible to generate overall vehicle models which are well tailored to particular applications. But, the numerical solution of coupled subsystems is not straight forward. This paper shows that the overall vehicle model can be solved very effectively by suitable interfaces and an implicit integration algorithm. The presented concept is realized in the product ve-DYNA, applied worldwide by automotive companies and suppliers.
Vehicle dynamics; vehicle model; axle modeling; drive train; multibody systems
TECHNICAL PAPERS
Vehicle modeling by subsystems
Georg Rill
GeorgRill@aol.com FH Regensburg, University of Applied Sciences Galgenbergstr. 30 D-93053 Regensburg, Germany
ABSTRACT
Computer simulations have become very popular in the automotive industry. In order to achieve a good conformity with field test, sophisticated vehicle models are needed. A real vehicle incorporates many complex dynamic systems, such as the drive train, the steering system and the wheel/axle suspension. On closer inspection some force elements such as shock absorbers and hydro-mounts turn out to be dynamic systems too. Modern vehicle models consist of different subsystems. Then, each subsystem may be modeled differently and can be tested independently. If some subsystems are available as a set of nested models of different complexity it will be even possible to generate overall vehicle models which are well tailored to particular applications. But, the numerical solution of coupled subsystems is not straight forward. This paper shows that the overall vehicle model can be solved very effectively by suitable interfaces and an implicit integration algorithm. The presented concept is realized in the product ve-DYNA, applied worldwide by automotive companies and suppliers.
Keywords: Vehicle dynamics, vehicle model, axle modeling, drive train, multibody systems
Modeling Concept
For dynamic simulation the vehicles are usually modeled by multi body systems (MBS), van der Jagt (2000). Typically, the overall vehicle model is separated into different subsystems, Rauh (2003). Fig. 1 shows the components of a passenger car model which can be used to investigate handling and ride properties. The vehicle model consists of the vehicle framework and subsystems for the steering system and the drive train.
The vehicle framework represents the kernel of the model. It at least includes the module chassis and modules for the wheel/axle suspension systems. The vehicle framework is supplemented by modules for the load, an elastically suspended engine, and passenger/seat models. A simple load module just takes the mass and inertia properties of the load into account. To describe the sloshing effects of liquid loads dynamic load models are needed, Rill and Rauh (1992). The subsystems elastically suspended engine, passenger/seat, and in heavy truck models a suspended driver's cabin can all be handled by the presented generic free body model. For standard vehicle dynamics analysis the chassis can be modeled by one rigid body. For applications where the chassis flexibility has to be taken into account a suitable flexible frame model is presented. Most wheel/axle suspension systems can be described by typical multi body system elements such as rigid bodies, links, joints and force elements, Rill (1994). Using a modified implicit Euler algorithm for solving the dynamic equations, axle suspensions with compliancies and dry friction in the damper element can be handled without any problems, Rill (2004). Due to their robustness leaf springs are still a popular choice for solid axles. They combine guidance and suspension properties which causes many problems in modeling, Fickers and Richter (1994). A leaf spring model is presented in this paper which overcomes these problems.
The steering system at least consists of the steering wheel, a flexible steering shaft, and the steering box which may also be power-assisted. Neureder (2002) has developed a very sophisticated model of the steering system which includes compliancies, dry friction, and clearance.
Tire forces and torques have a dominant influence on vehicle dynamics. The semi-empirical tire model TMeasy has mainly been developed to meet both the requirements of user-friendliness and sufficient model accuracy, Hirschberg et. al. (2002). Complex tire models such as the FTire Model provided by Gipser (1998) can be used for special applications. The module tire also includes the wheel rotation which acts as input for the drive train model. The presented drive train model is generic. It takes lockable differentials into account, and it combines front wheel, rear wheel and all wheel drive. The drive train is supplemented by a module describing the engine torque. It may be modeled quite simply by a first order differential equation or by the enhanced engine torque module en-DYNA developed by TESIS.
Road irregularities and variations in the coefficient of friction present significant impacts on the vehicle. A road model generating a two-dimensional reproducible random profile was provided by Rill (1990).
This modeling concept is realized with a MATLAB/Simulink® interface in the product ve-DYNA which also includes suitable models for the driver, TESIS.
Module Flexible Frame
Multi Body Approach to First Eigenmodes
The chassis eigenmodes of most passenger cars start at f >20Hz. Hence, for standard vehicle dynamic analysis the chassis can be modeled as one rigid body. The lower chassis stiffness of trucks and pickups results in eigenmodes starting at ¦ » 10 Hz, Fig. 2.
The first eigenmodes consist of torsion and bending of the chassis. These modes can be approximated by a multi body chassis model where the chassis is divided into three parts, Fig. 3.
Free Body Motions
The position and orientation of the reference frame xc, yc, zc which is fixed to the center body with respect to the inertial frame x0, y0, z0 is given by the rotation matrix
and the position vector
where the comma separated subscript 0 indicates that the coordinates of the vector from to are expressed in the inertial frame. The generalized coordinates roll, pitch, and yaw angle a0C, b0C, g0C as well as the coordinates x0C, y0C, z0C of the vector r0C,0 describe the free body motion of the vehicle.
Modal Coordinates
The motions of the front and rear body relative to the center body are small compared to the free body motions of the center part. Hence, the linearized rotation matrices
and the position vectors
are used to describe the orientation and position of the front and rear body relative to the center part. The vectors and denote the initial position of the front and rear body.
The generalized coordinates
describe the motions of the front and rear body relative to the center body. These motions are approximated by nM eigenmodes e1, e2, ... enM, now
where m1 , m2 , ..., mnM are modal coordinates, and EF and ER represent 6´nM matrices containing the eigenmodes.
Generalized Coordinates
The flexible chassis is modeled by 3 rigid bodies here. The orientation and the position of the bodies are described by free body motions and modal coordinates
where the 6 free body motions and the nM modal coordinates are collected in the vector yC. The dimension of depends on the number nM of modal coordinates, ny = 6 + nM.
Equations of Motion
To generate the equations of motion Jordain's Principle with generalized speeds is used. For a multi body system consisting of m rigid bodies it results in a set of two first order differential equations
The kinematic matrix K follows from the definition of the generalized speeds. The elements of the mass matrix M are given by
where mk is the mass and Qk the inertia tensor of body k. Finally, the components of the generalized forces and torques are defined by
where FAk, TAk denote the forces and torques applied to body k and , are remaining parts of the accelerations which do not depend on the derivatives of the generalized speeds.
Applied Forces and Torques
The forces and torques applied to the bodies can be written as
and
where the superscripts ext and cmp denote external and compliance forces and torques.
Applying Jordain's Principle, one part within the equations of motion describes the whole chassis motion. The compliance forces and torques are internal forces for the whole chassis and, therefore, do not show up in the corresponding parts of the generalized force vector.
If we assume that the compliance forces and torques are proportional to the motions of the front and rear body then, we will get
where cCF and cCR mark 6´6 stiffness matrices. The modal coordinate approximation Eq. 6 results in
Within Jordain's Principle the compliance forces and torques are reduced to generalized forces which are calculated by
where
are nM´nM stiffness matrices. which are defined by the modal stiffnesses cF1, cF2, ... cFnM and cR1, cR2, ... cRnM. Thus, to describe the motions of a flexible chassis only some eigenmodes and modal stiffnesses have to be provided.
Results
Depending on the vehicle layout, a flexible frame has a significant influence on the driving behavior, Fig 4. The rear axle of the considered bus is guided by four links. Here, the arrangement of the links generates a steering effect which depends on the roll angle of the rear part of the chassis and, therefore, also on the torsional stiffness of the chassis.
Module Leaf Spring
Modeling Aspects
Poor leaf spring models approximate guidance and suspension properties of the leaf spring by rigid links and separate force elements, Matschinsky (1998). The deformation of the leaf springs must be taken into account for realistic ride and handling simulations.
Within ADAMS leaf springs can be modeled with sophisticated beam-element models, ADAMS/Chassis 12.0. But, according to Fickers (1994) it is not easy to take the spring pretension into account. To model the effects of a beam, ADAMS/Solver uses a linear 6 -dimensional action-reaction force (3 translational and 3 rotational) between two markers. In order to provide adequate representation for the nonlinear cross section, usually 20 elements are used to model one leaf spring. A subsystem consisting of a solid axle and two beam-element leaf spring models would have ¦ = 6+2 * (20*6) = 246 degrees of freedom. In addition, the beam-element leaf spring model results in extremely stiff differential equations. These and the large number of degrees of freedom slow down the computing time significantly.
For real time applications the leaf springs must be modeled by a simple, but still accurate model. Fig. 5 shows a model of a solid axle with leaf spring suspension, which is typical for light truck rear axle suspension systems. There are no additional links. Hence, only the forces and torques generated by leaf spring deflections guide and suspend the axle.
The position of the axle center A and the orientation of an axle fixed reference frame xA, yA, zA are described relative to a chassis fixed frame xB, yB, zB by the displacements x, h, z and the rotation angles a, b, g which are collected in the 6´1 axle position vector
Similar to Fickers (1994) each leaf spring is modeled by five rigid bodies which are connected to each other by spherical joints, Fig. 5.
Each leaf spring is connected to the frame via the front leaf eye X. Furthermore each leaf spring is attached to the shackle at Y, and again to the frame at Z. In C the center part of each leaf spring is rigidly connected to the axle. The front eye bushings are modeled by spring/damper elements in x-, y-, and z-direction. The shackles are modeled by radial and a lateral spring/damper elements. Within each leaf spring the angles j1, y1, and j2, y2, describe the motions of part P-Q and part R-S relative to the center part. The outer parts Q-X and P-Y perform their rotations j3, y3, and j4, y4, relative to part P-Q and part R-S. As each leaf spring element is considered as a rigid rod, the roll motions can be neglected. The angles are collected in 4´1 position vectors
where y1F, y2F and y1R, y2R describe the momentary shape of the front and the rear part of the left (1) and the right (2) leaf spring.
A fully dynamic description of a solid axle with two five link leaf spring models would result in ¦ = 6+2*8=22 degrees of freedom. Compared to the beam-element model this is a really significant reduction.
But a dynamic description of the five link leaf spring model still includes some high frequent modes which will cause problems in the numerical solution of the equations of motion. As mass and inertia properties of the leaf spring model parts are small compared to the solid axle, a quasi static solution of the internal leaf spring deflection should be accurate enough within the overall vehicle model.
A quasi static solution provides the position vectors of the leaf spring parts as functions of the axle position vector, y1F = y1F (yA), y1R = y1R (yA) y2F = y2F (yA), y2R = y2R (yA). Hence, the subsystem solid axle with two leaf springs has only f=6 degrees of freedom.
Initial Shape and Pretension
At first it is assumed that the leaf spring is located in the xz-plane of the leaf spring fixed frame xL, yL, zL and its shape in the design position can be approximated by a circle which is fixed by the points X, C and Y. By dividing the arc X -Y into 5 parts of equal length the position of the links P, R, S, Q and the initial values of the angles j01, y01, j02, y02, j03, y03, j04, y04 can be calculated very easily.
In design position each leaf spring is only preloaded by a vertical load which results in zero pretension forces in the yL-direction, , and zero pretension torques around the zl-axis, , , , . In addition the torques around the xl-axis vanish, , , , .
To transfer the vertical preload F0 to the front eye bushing and the shackle, the joints P, Q, R, S must provide torques around the yL-axis, Fig. 6. The pretension forces in the front eye bushing , and in the shackle F0S, can easily be calculated from the equilibrium conditions of the five link leaf spring model,
where uYZ refers to the unit vector in the direction of the shackle, and , are the x and z components of the vector from pointing from X to Y. The pretension torques in the leaf spring joints around the yL-axis, , , , follow from
where rij , i = P,Q,R,S , j = X,Y are vectors pointing from i to j.
Compliance
The leaf spring compliance is defined in the design position by the vertical and the lateral stiffness, CV and CL. In Fig. 7a the leaf spring is approximated by a beam which is supported on both ends and is loaded in the center by the force F.
The deflection and the force F are related to each other via the stiffness c
If we transfer the beam model to the five link leaf spring model and look at the front half, Fig. 4, one will get
where a defines the length of one link, and small deflections in the xL, zL plane were assumed. The torques around the yL-axis in the joints P and Q are proportional to the deflection angles j1 and j3
The equilibrium condition results in
The leaf spring bending mode due to a single force can be approximated very well by a circular arc. Hence, the relative angle between connected links is equal, j1=j3=j and Eq. 23 can be simplified to w=3aj or or . From Eqs. 24 and 25 follows
Using Eq. 22, one finally obtains
where the beam stiffness C was replaced by the vertical leaf spring stiffness cV. Assuming symmetry, the stiffnesses in the rear joints are given by cj2 = cj1 and cj4 = cj3. The stiffnesses around the vertical axis cy1, cy2, cy3 and cy4 can be calculated in a similar way. The torsional stiffness of the leaf spring is neglected in this approach.
Actual Shape
The energy of a flexible system achieves a minimum value, E ® Min, in an equilibrium position. The energy of the five link leaf spring model is given by
where wX is the 3´1 displacement vector and cB is the 3´3 stiffness matrix of the front eye bushing, wSR, wSL are the radial and lateral shackle displacements, and cSR, cSLdenote the corresponding stiffnesses.
According to Eqs. 18 and 19, the actual shape of the leaf spring is determined by the position vectors y1 = [j1,y1,j3,y3]T and y2 = [j2,y2,j4,y4]T. If the leaf spring energy becomes a minimum, the following equations will hold
As the shackle displacements wSR, wSL do not depend on y1 and the front bushing displacement vector wX does not depend on y2 the conditions in Eq. 29 form two independent sets of nonlinear equations ¦1(y1,yA) = 0 and ¦2(y2,yA) = 0, where yA denotes the dependency of the actual position and orientation of the solid axle. These equations are solved iteratively by the Newton-Algorithm. Starting with initial guesses , one gets an improvement by solving the linear equations
Here, the Jacobians can be calculated analytically.
Leaf Spring Reaction Forces
The actual forces in the front leaf eye bushing are given by
where F0B denotes the pretension force and CB, dB are 3´3 matrices, characterizing the stiffness and damping properties of the front leaf eye bushing. The displacement vector wx in the front leaf eye bushing depend on the generalized coordinates y1 and yA which describe the actual shape of the front leaf spring part and the actual position and orientation of the solid axle. By solving Eq. 30, y1 is given as a function of yA. Hence, wx only depends on yA and its derivative can be determined by
where
A describes the velocity state of the solid axle.The radial and lateral components of the shackle forces can be calculated from
where F0S represents the pretension force, uSR, uSL are unit vectors in the radial and lateral shackle direction, and cSR, cSL, dSR, dSL are constants characterizing the stiffness and damping properties of the shackle. The shackle displacements wSR and wSL depend on the generalized coordinates y2R and yA which describe the actual shape of the rear leaf spring part and the actual position and orientation of the solid axle. Similarly to Eq. 32 the displacement velocities are given by
Finally, the shackle force reads as
Forces Applied to the Axle
The leaf springs act like generalized force elements in this approach, Fig. 8. Guidance and suspension of the solid axle is done by the resulting force
and the resulting torque
where rAB1 = rAB1(yA) , ... rAS2(yA) describe the momentary position of the front eye bushings and the shackles relative to the axle center.
As the forces in the front eye bushings FB1, FB2 and the shackle forces FS1, FS2 depend on the axle state yA, A only, the resulting force and the resulting torque are also functions of the axle state only. Since, hereby each leaf spring acts as a generalized force element, it can easily be integrated into the vehicle framework. By suppressing high frequent leaf spring eigenmodes, it is perfectly adopted to real-time application.
Bending Modes
The quasi-static approach reproduces all significant bending modes of the leaf spring, Fig.9. A leaf spring is stiffer in the lateral direction than in the vertical direction. Hence, a displacement in the front eye bushing is noticeable only on lateral leaf spring deflections.
Model Performance
The five link leaf spring model was integrated into a ve-DYNA Ford Transit vehicle model.
Using this model at the rear axle instead of a poor kinematic approach means only 85% more computer run time. Hence, real time applications are still possible. The simulation results are in good conformity to measurements, Fig. 10. The nonlinearity in the spring characteristics is caused by an additional bump stop and by the change of the shackle position during jounce and rebound. Obviously, the five link model is accurate enough.
Free Body Module
Position and Orientation
To describe the momentary state of the body E the frame xE yE, zE located in the center of gravity is used. In addition, sensor points S monitor position, velocity, and acceleration at specific body points, Fig. 11.
The frame B is fixed to the vehicle. The suspension of body E on the vehicle, frame B may consist of force elements and/or rubber mounts. The road-fixed frame 0 is considered as an inertial frame. The position of frame B with respect to the road-fixed inertial frame 0 is given by the position vector
The orientation of the frame axes is described by a rotation matrix. Three elementary rotations are put together. The sequence
results in
Hence, position and orientation of the vehicle-fixed reference frame are described by 6 generalized coordinates xB,yB,zB and aB, bB, gB.
The position and orientation of the elastically suspended body with respect to the reference frame B is defined by
and
Generalized Speeds
The velocity of the reference frame B with respect to the inertial frame 0 is given by
The velocity denoted in the inertial frame can be transformed to the reference frame
By this, the orthogonality of the rotation matrix
was already taken into consideration.
The angular velocity of the reference frame B with respect to the inertial frame may be expressed directly in reference frame B
The 6 components of n0B,B and w0B,B will be chosen as generalized speeds now. First order kinematical differential equations connect this speeds with the derivatives of the
and
where the solution of Eq. 48 is given by
The momentary state of the reference frame B is fully characterized by 6 generalized coordinates xB, yB, zB, aB, bB, gB and 6 generalized speeds n0Bx, n0By, n0Bz, w0Bx, w0By, w0Bz.
The velocity and the angular velocity of the elastically suspended body with respect to the inertial frame 0 is given by
where the derivative of the position vector and the angular velocity of the elastically suspended body follow from the Eqs. 41 and 42. They read as
and
By using the components of the velocity
and the angular velocity
as generalized speeds, Eq. 50 can be written as a set of kinematical differential equations
and
Where as the 6 generalized coordinates xE, yE, zE, aE, bE, gE describe the position and orientation of frame E relative to frame B, the 6 generalized speeds n0Ex, n0Ey, n0Ez, w0Ex, w0Ey, w0Ez are the components of the absolute velocity and angular velocity of body E.
Accelerations
The accelerations of body E with respect to the inertia frame 0 can be expressed in reference frame B. They read as
where
and
follow from the Eqs. 53 and 54.
Force Elements
If a force element is attached to the chassis at point i and to the body at point j, the momentary position of force element ij will be defined by
where
rBi,K , rEj,K are given by data, and rBE,B follows from Eq. 41. The actual length can be calculated from
and the unit vector
describes the momentary direction of the force element. If denotes the initial length of the force element, the displacement of the force element will be formed by
The displacement velocity follows from
Using Eq. 60, Eq. 61, and
Bi,K = 0 one gets
where
BE,B and wBE,B are given by the Eqs. 51 and 52. The forces Fij,B, Fji,B and the torques Tij,B, Tji,B, applied to body and chassis are determined by
and
where ¦ describes an arbitrary spring/damper characteristic.
Equations of Motion
Applying liner and angular momentum to the elastically suspended body, one obtains
and
where mE, QE,B denote mass and inertia tensor of the free body, FE,B, TE,B are the resulting forces and torques applied to the free body, and g,B is the vector of gravity expressed in the body fixed reference frame. These equations are coupled with the chassis equations of motion only by the applied forces and torques. Due to the particular choice of generalized speeds, no mass or inertia coupling terms appear.
By using this modeling technique, Seibert and Rill (1998) showed that the comfort of a passenger car is significantly influenced by the engine suspension system. The free body model can also be used to model an elastically suspended driver's cab, Rill (1993).
Subsystem Drive Train
Generic Model Structure
The subsystem drive train, Fig. 12, interacts on one hand with the engine and on the other hand with the wheels. Hence, the angular velocities of the wheels w1, , w4 and the engine or the gear output angular velocity w0 respectively are input quantities. For this reason, the calculation of the engine torque and the dynamics of the wheel rotation are performed in other subsystems. Via the tire forces and torques, the drive train is coupled with the steering system and the vehicle frame work.
The drive train model includes three lockable differentials. The angular velocities of the drive shafts wS1: front left, wS2: front right, wSF: front, wSR: rear, wS3: rear left, wS4: rear right are used as generalized coordinates.
The torque distribution of the front and rear differential is 1:1. If rF and rR are the ratios of the front and rear differential, one will get
The torque distribution of the center differential is formed by
where tF, tR denote the torques transmitted to the front and rear drive shaft, and µ is a dimensionless drive train parameter. A value of µ = 1 means front wheel drive, 0< µ <1 stands for all wheel drive, and µ = 0 is rear wheel drive. If the ratio of the center differential is given by rC then
holds.
Equation of Motion
The equation of motion for the drive train is derived from Jordain's Principle, which reads as
where Qi is the inertia of body i , i denotes the time derivative of the angular velocity, ti is the torque applied to each body, and dwidescribe the variation of the angular velocity. Applying Eq. 75 for the different parts of the drive train model results in
Collecting all terms with dwS1, dwS2, dwSF, dwSR, dwS3, dwS4 and using the abbreviations n=1-µ, , and finally leads to three blocks of differential equations
which describe the dynamics of the drive train. Due to its simple structure, an extension to a 6x6 or 8x8 drive train will be straight forward.
Drive Shaft Torques
The torques in the drive shafts are given by
and cS0, cS1, cS2, cS3, cS4, cSF, cSR denote the stiffnesses of the drive shafts. The first order differential equations can be arranged in matrix form
where
represents the vector of the angular velocities,
contains the torsional angles in the drive shafts,
is the excitation vector, and
forms a 7x6 distribution matrix.
Locking Torques
The differential locking torques are created by an enhanced dry friction model consisting of a static and a dynamic part
The dynamic parts are modeled by a torque proportional to the differential output angular velocities
where dLF, dLC, dLR are damping parameters which have to be chosen appropriately. In steady state operating conditions, the static parts , , will provide torques, even if the differential output angular velocities are equal. From the Eqs. 80, 81, and 82, one gets
By this locking torque model, the effect of dry friction inside the differentials can also be taken into account.
Numerical Solution
The equations of motion 80, 81, and 82 can be combined in a matrix differential equation
where w, Dj are given by the Eqs. 85, 86, and the mass matrix M is built by three 2x2 submatrices
where the elements of MF, MC, and MR follow from the Eqs. 80, 81, and 82. The vector of the generalized torques is written as
Because the model also includes the high frequent drive shaft vibrations the differential equations for the drive train are stiff. Hence, implicit integration algorithm should be used for the numerical solution. Vehicle dynamic equations can be solved very effectively by a modified implicit Euler algorithm, Rill (2004).
The implicit Euler-formalism for Eq. 92 and Eq. 84 results in
where h is the integration step size, and the superscripts k and k+1 indicate the states at t and t+1. Applying the Taylor-expansion to q at and wk, one gets
By using the Eqs. 84 and 96 the second term on the right side can be written as
Now, the implicit algorithm in Eq. 95 can be approximated by
which finally results in
where the partial derivatives ¶q/¶Dq and ¶q/¶q can be calculated quite easily.
Partial Derivatives
Only the dynamic locking torques , and depend on the angular velocities. Hence, one obtains
The change of q with respect to Dj leads to a 6x7 matrix
The term which is finally needed in Eq. 100 is symmetric and reads as
where the abbreviations
were used.
System Performance
Locking the differential improves the traction of a vehicle. In Fig. 13 the simulation results of a vehicle with rear wheel drive starting on a µ-split surface are shown.
At first all differentials are unlocked. The left rear wheel which is running on a low µ-plate immediately starts spinning. At t = 2.5 s the rear differential is locked. Now, the locking torque, which is generated by the drive train model forces both wheels to run with the same angular velocity.
Conclusion
Vehicle modeling by subsystems make a large variety of applications possible. The combination of simple subsystems and modules results in a vehicle model with a minimum number of data and a very good run time performance. Such "light models" can be used to develop enhanced control strategies for electronic safety devices. Depending on the focus of interest, more and more subsystems and modules may be replaced by enhanced ones. Then, sophisticated design studies or a comfort analysis are possible. If the modified implicit Euler algorithm is also applied to the critical subsystems drive train and steering system the numerical solution of the overall vehicle model still will not be time consuming.
References
[1] ADAMS/Chassis 12.0 Reference Guide.
[15] tesis: www.tesis.de.
Paper accepted: June, 2005.
Presented at XI DINAME International Symposium on Dynamic Problems of Mechanics, February 28th - March 4th, 2005, Ouro Preto. MG. Brazil.
Technical Editors: J.R.F. Arruda and D.A. Rade.
- [2] Fickers, P.; Richter, B.: Incorporating FEA-Techniques into MSA illustrated by several Rear Suspension Concepts. In: 9th European ADAMS User Conference, Frankfurt, November 21st/22nd, 1994.
- [3] Gipser, M.: Reifenmodelle für Komfort- und Schlechtwegsimulationen. In: Tagungsband zum 7. Aachener Kolloquium Fahrzeug- und Motorentechnik. IKA, RWTH Aachen und VDI 1998.
- [4] Hirschberg, W.; Rill, G.; Weinfurter, H.: User-Appropriate Tyre-Modeling for Vehicle Dynamics in Standard and Limit Situations. Vehicle System Dynamics 2002, Vol. 38, No. 2, pp. 103-125. Lisse: Swets & Zeitlinger.
- [5] van der Jagt, P.: The Road to Virtual Vehicle Prototyping; new CAE-models for accelerated vehicle dynamics development. PhD-Thesis, Tech. Univ. Eindhofen, Eindhofen 2000, ISBN 90-386-2552-9 NUGI 834.
- [6] Matschinsky, W.: Radführungen der Straßenfahrzeuge. Springer, Berlin 1998.
- [7] Neureder, U.: Untersuchungen zur Übertragung von Radlastschwankungen auf die Lenkung von Pkw mit Federbeinvorderachse und Zahnstangenlenkung. Fortschritt-Berichte VDI, Reihe 12, Nr. 518. Düsseldorf: VDI Verlag 2002.
- [8] Rauh, J.: Virtual Development of Ride and Handling Characteristics for Advanced Passenger Cars. Vehicle System Dynamics, 2003, Vol. 40, Nos. 1-3, pp. 135-155.
- [9] Rill, G.: Demands on Vehicle Modeling. In: The Dynamics of Vehicles on Road and on Tracks. Ed.: Anderson, R.J., Lisse: Swets-Zeitlinger 1990.
- [10] Rill, G.; Rauh, J.: Simulation von Tankfahrzeugen. In: Berechnung im Automobilbau, VDI-Bericht 1007. Düsseldorf: VDI-Verlag 1992.
- [11] Rill, G.: Modeling and Dynamic Optimization of Heavy Agricultural Tractors. In: 26th International Symposium on Automotive Technology and Automation (ISATA). Croydon: Automotive Automation Limited 1993.
- [12] Rill, G.: Simulation von Kraftfahrzeugen. Vieweg, Braunschweig 1994.
- [13] Rill, G.: A Modified Implicit Euler Algorithm for Solving Vehicle Dynamics Equation. To Appear In: Multibody System Dynamics, 2004.
- [14] Seibert, Th.; Rill, G.: Fahrkomfortberechnungen unter Einbeziehung der Motorschwingungen. In: Berechnung und Simulation im Fahrzeugbau, VDI-Bericht 1411. Düsseldorf: VDI-Verlag 1998.
Publication Dates
-
Publication in this collection
08 Oct 2007 -
Date of issue
Dec 2006
History
-
Accepted
June 2005 -
Received
June 2005