The paper describes the design and multibody dynamic analysis of a mechanically interconnected suspension, as applied to a small off-road vehicle. Interconnected suspensions use some sort of connection between the axles of a vehicle in order improve ride quality or vehicle handling. In principle, the connection may be hydraulic, pneumatic, or mechanical, but for installation in a typical passenger car, a mechanical connection would likely be impractical due to weight and complexity. In this paper, the vehicle in question is the University of Windsor SAE Baja off-road competition vehicle, and novel mechanical design is proposed. A multibody dynamic analysis is performed on the proposed design using the EoM open source multibody software developed by theUniversity ofWindsorVehicle Dynamics and Control research group in order to assess any potential performance improvements.
The goal of this research is to achieve close to real-time dynamics performance for allowing auto-pilot in-the-loop testing of unmanned ground vehicles (UGV) for urban as well as off-road scenarios. The overall vehicle dynamics performance is governed by the multibody dynamics model for the vehicle, the wheel/terrain interaction dynamics and the onboard control system. The topic of this paper is the development of computationally efficient and accurate dynamics model for ground vehicles with complex suspension dynamics. A challenge is that typical vehicle suspensions involve closed-chain loops which require expensive DAE integration techniques. In this paper, we illustrate the use the alternative constraint embedding technique to reduce the cost and improve the accuracy of the dynamics model for the vehicle.
This paper compares selected optimization-based methods for the analysis of multibody systems with redundant constraints. The following numerical schemes are examined: direct integration method, Udwadia-Kalaba formulation, two types of least-squares block solution method and Udwadia-Phohomsiri formulation. In order to compare efficiency of the algorithms, a series of simulations is performed on two exemplary McPherson struts. In the first variant, the mechanism has no redundant constraints whereas the other is overconstrained. Three constraint stabilization schemes are also compared in terms of integration errors.
This paper presents a methodology for contact detection between convex quadric surfaces using its implicit equations. With some small modifications in the equations, one can model superellipsoids, superhyperboloids of one or two sheets and supertoroids. This methodology is to be implemented on a multibody dynamics code, in order to simulate the interpenetration between mechanical systems, particularly, the simulation of collisions with motor vehicles and other road users, such as cars, motorcycles and pedestrians. The contact detection of two bodies is formulated as a convex nonlinear constrained optimization problem that is solved using two methods, an Interior Point method (IP) and a Sequential Quadratic Programming method (SQP), coded in MATLAB and FORTRAN environment, respectively. The objective function to be minimized is the distance between both surfaces. The design constraints are the implicit superquadrics surfaces equations and operations between its normal vectors and the distance itself. The contact points or the points that minimize the distance between the surfaces are the design variables. Computational efficiency can be improved by using Bounding Volumes in contact detection pre-steps. First one approximate the geometry using spheres, and then Oriented Bounding Boxes (OBB). Results show that the optimization technique suits for the accurate contact detection between objects modelled by implicit superquadric equations.
The dynamics of the turning process of a thin-walled cylinder in manufacturing is modeled using flexible multibody system theory. The obtained model is time varying due to workpiece rotation and tool feed and retarded, due to repeated cutting of the same surface. Instabilities can occur due to these consecutive cuts that must be avoided in practical application because of the detrimental effects on workpiece, tool and possibly the machine. Neglecting the small feed, the stability of the resulting periodic system with time-delay can be analyzed using the semi-discretization method. The use of an adaptronic tool holder comprising actuators and sensors to improve the dynamic stability is then investigated. Different control concepts, two collocated and two model-based, are implemented in simulation and tuned to increase the domain of stable cutting. Cutting of a moderately thin workpiece exhibits instabilities mainly due to tool vibration. In this case, the stability boundary can be significantly improved. When the instability is due to workpiece vibration, the collocated concepts fail completely. Model based concepts can still obtain some improvements, but are sensitive to modeling errors in the coupling of workpiece and tool.
The purpose of the present research relates to the sensitivity analysis of road vehicle comfort and handling performances with respect to suspension technological parameters. The envisaged suspension being of semi-active nature, this implies first to consider an hybrid modeling approach consisting of a 3D multibody model of the full car - an Audi A6 in our case - coupled with the electro-hydraulic model of the suspension dampers. Concerning parameter sensitivitie, the goal is to capture them for themselves - and not necessarily for optimization purpose - because their knowledge is of a great interest for the damper manufacturer. An important issue of the research is to consider objective functions which are based on complete time integrations along a given trajectory, the goal being - for instance - to quantify the sensitivity of the carbody rms acceleration (comfort) or of the vehicle overturning character (handling) with respect to suspension parameters. On one hand, the accuracy of the various partial derivatives computation can be greatly enhanced thanks to the symbolic capabilities of our ROBOTRAN multibody program. On the other hand, the computational efficiency of the process also takes advantage of the recursive formulation of the multibody equations of motion which must be time integrated with respect to both the generalized coordinates and their partial derivatives in case of the so-called direct method underlying sensitivity analysis.
In the present work, a tire model is derived based on geometrically exact shells. The discretization is done with the help of isoparametric quadrilateral finite elements. The interpolation is performed with bilinear Lagrangian polynomials for the midsurface as well as for the director field. As time stepping method for the resulting differential algebraic equation a backward differentiation formula is chosen. A multilayer material model for geometrically exact shells is introduced, to describe the anisotropic behavior of the tire material. To handle the interaction with a rigid road surface, a unilateral frictional contact formulation is introduced. Therein a special surface to surface contact element is developed, which rebuilds the shape of the tire.
A formulation developed at the Laboratory of Mechanical Engineering allows robust and efficient simulation of large and complex multibody systems. Simulators of cars, excavators and other systems have been developed showing that real-time simulations are possible even when facing demanding manoeuvres. Hydraulic actuators are presented in many industrial applications of multibody systems, like in the case of the heavy machinery field. When simulating the dynamics of this kind of problems that combine multibody dynamics and hydraulics, two different approaches are common: to resort to kinematically guide the variable length of the actuator, thus avoiding the need to consider the dynamics of the hydraulic system; or to perform a multi-rate integration of both subsystems if a more detailed description of the problem is required, for example, when the objective of the study is to optimize the pump control. This work addresses the inclusion of hydraulic actuators dynamics in the above-mentioned self-developed multibody formulation, thus leading to a unified approach. An academic example serves to compare the efficiency, accuracy and ease of implementation of the simplified (kinematic guidance), multi-rate and unified approaches. Such a comparison is the main contribution of the paper, as it may serve to provide guidelines on which approach to select depending on the problem characteristics.
One of the applications of tether system is in the field of satellite technology, where the mother ship and satellite equipment are connected with a cable. In order to grasp the motion of this kind of tether system in detail, the tether can be effectively modeled as flexible body and dealt by multibody dynamic analysis. In the analysis and modeling of flexible body of tether, large deformation and large displacement must be considered. Multibody dynamic analysis such as Absolute Nodal Coordinate Formulation with an introduction of the effect of damping force formulation can be used to describe the motion behavior of a flexible body. In this study, a parameter identification technique via an experimental approach is proposed in order to verify the modeling method. An example of swing-up control using the genetic algorithm control approach is performed through simulation and experiment. The validity of the model and availability of motion control based on multibody dynamics analysis are shown by comparison between numerical simulation and experiment.
The beam elements, which are widely used in the absolute nodal coordinate formulation (ANCF) can be treated as isoparametric elements, and by analogy to the classical finite element analysis (FEA) are integrated with standard, spatial Gauss-Legendre quadratures. For this reason, the shape of the ANCF beam cross section is restricted only to the shape of rectangle. In this paper, a distinct method of integration of ANCF elements based on continuum mechanics approach is presented. This method allows for efficient analysis of the ANCF beam elements with circular cross section. The integration of element vectors and matrices is performed by separation of the quadrature into the part that integrate along beam axis and the part that integrate in the beam cross section. Then, an alternative quadrature is used to integrate in the circular shape of the cross section. Since the number of integration points in the alternative quadrature corresponds to the number of points in the standard Gaussian quadrature the change in the shape of the cross section does not affects negatively the element efficiency. The presented method was verified using selected numerical tests. They show good relatively agreement with the reference results. Apart from the analysis of the beams with the circular cross section, a possibility of further modifications in the methods of the element integration is also discussed. Due to the fact that locking influence on the convergence of the element is also observed, the methods of locking elimination in the proposed elements are also considered in the paper.