A reduced beam and joint concept modeling approach to optimize global vehicle body dynamics

Ideally, NVH simulations become available already in the concept phase of vehicle development. The initial computer-aided design (CAD) can then be improved (by already including countermeasures), and the feasibility to balance NVH with other performance attributes is increased. In this early design stage, when exact geometrical information is not or scarcely available, conventional virtual prototyping techniques based on detailed CAD and FE models are not directly applicable. A state-of-the-art overview of concept NVH simulation methods in vehicle industry is given.This paper then presents a “Reduced Beam and Joint Modeling” approach to analyze and optimize the global bending and torsion modes of a vehicle body. Concept modifications in the body beam-like sections and in the joints are analyzed using the body reduced modal model. This small-sized model can be used to quickly and accurately optimize the low-frequency vehicle performance. The modifications are considered with respect to the existing (predecessor) model. Equivalent beam properties are estimated from the body FE model; modifications in the beam-like sections are then implemented with beam elements from a standard FE library. The joint modifications are considered through static superelements: stiffness formulations between the end points of the joint connected to the beam layout. The validity of the approach is first demonstrated on simple example models. An industrial vehicle BIW application case is subsequently presented. A beam and joint layout is created, and used for a fast sensitivity analysis to identify suitable modifications to improve the global modes. Next, two application cases are presented. First, a fast optimization analysis is performed to optimize the global body modes. Subsequently, suitable physical modifications are identified and applied to the full FE model; it is shown that the effect of these physical modifications is accurately predicted with the fast sensitivity analysis.     

Efficient optimisation of a vehicle suspension system, using a gradient-based approximation method, Part 1: Mathematical modelling

A methodology is proposed for the efficient determination of gradient information, when performing gradient based optimisation of an off-road vehicle’s suspension system. The methodology is applied to a computationally expensive, non-linear vehicle model, that exhibits severe numerical noise. A recreational off-road vehicle is modelled in MSC.ADAMS, and coupled to MATLAB for the execution of the optimisation. The successive approximation method, Dynamic-Q, is used for the optimisation of the spring and damper characteristics. Optimisation is performed for both ride comfort and handling. The determination of the objective function value is performed using computationally expensive numerical simulations.This paper proposes a non-linear pitch-plane model, to be used for the gradient information, when optimising ride comfort. When optimising for handling, a non-linear four wheel model, that includes roll, is used. The gradients of the objective function and constraint functions are obtained through the use of central finite differences, within Dynamic-Q, via numerical simulation using the proposed simplified models. The importance of correctly scaling these simplified models is emphasised. The models are validated against experimental results. The simplified vehicle models exhibit significantly less numerical noise than the full vehicle simulation model, and solve in significantly less computational time.     

Beam element modelling of vehicle body-in-white applying artificial neural network

In this study a large knowledge base is first established through numerous designs of experiments on beam elements, based on a validated finite element model of a reference vehicle body-in-white. Then a recurrent artificial neural network is applied to extract the input/output relationship between the crash dynamic characteristics and beam element features. With such established relationship, beam element features are predicted according to expected crash dynamic characteristics. Our analyses show that the predicted beam element model enables generating essential crash dynamic characteristics for concept BIW design evaluation at a reasonable level of accuracy. Last, a data assurance criterion is developed to quantitatively validate the beam element modelling.     

动态需求下蓄冷式多温共配多目标优化模型及算法

蓄冷式多温共配是当前冷链物流领域一种较新的配送模式,可以有效地满足客户对产品多样性的需求。针对蓄冷式多温共配中客户需求的动态变化,本文从需求预测的角度出发,首先根据服务客户点的不同将车辆冗余空间的装载比例做了差异化处理,实现合理利用车辆冗余空间而又不全部满载;进而从时间满意度和订单完成率两方面度量客户满意度,建立了蓄冷式多温共配多目标模型;最后构建了较为高效的算法进行求解。算例分析表明,和其他策略相比,本文方法实用性更强。     

Nonlinear finite element analyses of concrete bridge joint systems subjected to seismic actions

Following widespread damage to bridge joints in the San Francisco region from the 1989 Loma Prieta earthquake, the necessity for establishing an alternative method for seismic design of bridge joints was identified. Recognizing that conventional joint design practice based directly on shear forces results in congested reinforcement details, which are difficult to implement in practice, a rational design procedure was sought through large-scale testing of bridge joint systems and subsequent finite element and strut-and-tie analyses. The finite element part of the study is presented in this paper, which focuses on (a) identification of compression force flow and thus the load path across the joint, (b) examination of an efficient joint force transfer model, and (c) influence of cap beam prestressing. Combining the experimental and analytical results, a joint design method has been established in which reduction of joint reinforcement was achieved by treating joint shear as part of the complete force transfer across the joint, rather than as an independent action. The proposed design approach has been validated in a laboratory test on a full-scale multiple-column bridge bent.     

An improved modeling approach for circular flexure hinges with application to geometry optimization

Within this paper, a modeling approach for flexure hinges based on the Euler-Bernoulli beam theory for beams of variable cross section is investigated in a static analysis. The proposed approach is implemented in a finite beam element routine, for which two different discretizations are discussed. The results are compared to a full scale three dimensional model. It is shown that a circular flexure hinge cannot be modeled accurately with one element. An improved model with three elements across the flexure hinge length is presented which shows excellent accordance with the reference model. A geometry optimization is realized based on the improved, low-DOF model. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Model order reduction of finite element models: improved component mode synthesis

The finite element (FE) approach constitutes an essential methodology when modelling the elastic properties of structures in various research disciplines such as structural mechanics, engine dynamics and so on. Because of increased accuracy requirements, the FE method results in discretized models, which are described by higher order ordinary differential equations, or, in FE terms, by a large number of degrees of freedom (DoF). In this regard, the application of an additional methodology, referred to as the model order reduction (MOR) or DoF condensation, is rather compulsory. Herein, a reduced dimension set of ordinary differential equations is generated, i.e. the initially large number of DoF is condensed, while aiming to keep the dynamics of the original model as intact as possible. In the commercially available FE software tools, the static and the component mode syntheses (CMS) are the only available integrated condensation methods. The latter represents the state of the art generating well-correlated reduced order models (ROMs), which can be further utilized for FE or multi-body systems simulations. Taking into consideration the information loss of the CMS, which is introduced by its part-static nature, the improved CMS (ICMS) method is proposed. Here the algorithmic scheme of the standard CMS is adopted, which is qualitatively improved by adequately considering the advantageous characteristics of another MOR approach, the so-called improved reduction system method. The ICMS results in better correlated reduced order models in comparison to all the aforementioned methods, while preserving the required structural properties of the original FE model.     

Analytical approximation of the transition density in a local volatility model

We present a simplified approach to the analytical approximation of the transition density related to a general local volatility model. The methodology is sufficiently flexible to be extended to time-dependent coefficients, multi-dimensional stochastic volatility models, degenerate parabolic PDEs related to Asian options and also to include jumps.     

Power flow analysis of beam members with multiple wave types

A power flow finite element method (PFFEM) is developed by Nefske and Sung (ASME J. Vib. Stress Reliability Des. 111, 1989) for the steady state high frequency vibration analysis of beams in which there is propagation of one wave type (bending). A finite element thermal analog presented herein for power flow analysis of beam members in which both bending and longitudinal waves are present is a further development of the PFFEM. Conducting bar elements model mechanical power flow in the beam and also power transmission, reflection and conversion at structural joints. Convection elements model internal power dissipation. Properties of joint elements are in general non-linear and depend upon the power level at the joint as well as transmission and reflection coefficients. The thermal analog is verified for particular linear cases wherein properties of the joint elements are dependent only upon the transmission and reflection coefficients.     

Comparison of the simplified and full MHD models for laminar incompressible flow past a circular cylinder

In the majority of research on incompressible magnetohydrodynamic (MHD) flows, the simplified model with the low magnetic Reynolds number assumption has been adopted because it reduces the number of equations to be solved. However, because the effect of flow on magnetic field is also neglected, the solutions of the simplified model may be different from those of the full model. As an example, the flow of an electrically conducting fluid past a circular cylinder under a magnetic field is investigated numerically using the simplified and full models in this paper. To solve the problems, two second-order compact finite difference algorithms based on the streamfunction-velocity formulation of the simplified model and the quasi-streamfunction-velocity formulation of the full model are developed respectively.Numerical simulations are carried out over a wide range of Hartmann number for steady-state laminar problems with both models. For the full model, magnetic Reynolds number (Re_m) is chosen from 0.01 to 10. The computed results show that solutions of the simplified MHD model are not exactly the same as those of the full MHD model for this flow problem in most cases even if Re_m in the full model is very low. Only in the special case that a strong external magnetic field is exerted perpendicular to the dominant flow direction, can the simplified MHD model be regarded as an approximation of the full MHD model with low Re_m.     

The Fshape Framework for the Variability Analysis of Functional Shapes

This article introduces a full mathematical and numerical framework for treating functional shapes (or fshapes) following the landmarks of shape spaces and shape analysis. Functional shapes can be described as signal functions supported on varying geometrical supports. Analyzing variability of fshapes’ ensembles requires the modeling and quantification of joint variations in geometry and signal, which have been treated separately in previous approaches. Instead, building on the ideas of shape spaces for purely geometrical objects, we propose the extended concept of fshape bundles and define Riemannian metrics for fshape metamorphoses to model geometric-functional transformations within these bundles. We also generalize previous works on data attachment terms based on the notion of varifolds and demonstrate the utility of these distances. Based on these, we propose variational formulations of the atlas estimation problem on populations of fshapes and prove existence of solutions for the different models. The second part of the article examines thoroughly the numerical implementation of the tangential simplified metamorphosis model by detailing discrete expressions for the metrics and gradients and proposing an optimization scheme for the atlas estimation problem. We present a few results of the methodology on a synthetic dataset as well as on a population of retinal membranes with thickness maps.     

Joint optimization of fleet size and maintenance capacity in a fork-join cyclical transportation system

This article presents an asset management-oriented multi-criteria methodology for the joint estimation of a mobile equipment fleet size, and the maintenance capacity to be allocated in a productive system. Using a business-centred life-cycle perspective, we propose an integrated analytical model and evaluate it using global cost rate, availability and throughput as performance indicators. The global cost components include: (i) opportunity costs associated with lost production, (ii) vehicle idle time costs, and (iii) maintenance resources idle time costs. This multi-criteria approach allows a balanced scorecard to be built that identifies the main trade-offs in the system. The methodology uses an improved closed network queueing model approach to describe the production and maintenance areas. We test the proposed methodology using an underground mining operation case study. The decision variables are the size of a load-haul-dump fleet and specialized maintenance crew levels. Our model achieves savings of 20.6% in global cost terms with respect to a benchmark case. We also optimize the system to achieve desired targets of vehicle availability and system throughput (based on system utilization). The results show increments of 7.1% in vehicle availability and 13.5% in system throughput with respect to baseline case. For the case studied, these criteria also have a maximum, which allows for further improvement if desired. The results also show the importance of using balanced performance measures in the decision process. A multi-criteria optimization was also performed, showing the Pareto front of considered indicators. We discuss the trade-offs among different criteria, and the implications in finding balanced solutions. The proposed analytical approach is easy to implement and requires low computational effort. It also allows for an easy re-evaluation of resources when the business cycle changes and relevant exogenous factors vary.     

Modelling and simulation of gas dynamics in an exhaust pipe

The gas dynamics in an exhaust pipe is studied. In particular we focus on the warm up of the catalytic converter in very short times after the engine start. This is done by combustion a small unburnt part of the exhaust gas. This process is classically modelled by gas dynamic equations. Compared to the existing literature we improve the (one-dimensional) modelling approach using a small Mach number technique and a network ansatz for the full exhaust pipe. The final simplified model on one hand still describes the main features and on the other hand it is computationally a few orders of magnitude faster than the original model. Performing numerical simulations we compare the new model to the (classical) full model and to experimental results in the literature.     

Empirical Modeling of a Hydraulic Actuator in a Vehicle Control System

The pressure information of a hydraulic actuator plays a key role in modern vehicle control and fault diagnosis. The difficulty in measuring pressure directly naturally motivates such an indirect approach as an observer whose accuracy depends heavily upon the availability of a high fidelity model of a hydraulic actuator. Notwithstanding its success in understanding the dynamics of the first principle model of a hydraulic actuator, it is not suited to controller/observer design due to its complexity. This paper presents an alternative to the first principle modeling methodology: an empirical approach to hydraulic actuator modeling. Linear and nonlinear system identification techniques are applied to obtain low-order models of a hydraulic actuator. Experimental results show how the empirical models reproduce the key features of a hydraulic actuator with sufficient accuracy.     

Analytical elastostatic stiffness modeling of parallel manipulators considering the compliance of the link and joint

This paper discusses the analytical elastostatic stiffness modeling of parallel manipulators (PMs) considering the compliance of the link and joint. The proposed modeling is implemented in three steps: (1) the limb constraint wrenches are formulated based on screw theory; (2) the strain energy of the link and the joint is formulated using material mechanics and a mapping matrix, respectively, and the concentrated limb stiffness matrix corresponding to the constraint wrenches is obtained by summing the strain energy of the links and joints in the limb; and (3) the overall stiffness matrix is assembled based on the deformation compatibility equations. The strain energy factor index (SEFI) is adopted to describe the influence of the elastic components on the stiffness performance of the mechanism. Matrix structural analysis (MSA) using Timoshenko beam elements is applied to obtain analytical expressions for the compliance matrices of different joints through a three-step process: (1) formulate the element stiffness equation for each element; (2) extend the element stiffness equation to obtain the element contribution matrix, allowing the extended overall stiffness matrix to be obtained by summing the element contribution matrices; and (3) determine the stiffness matrices of joints by extracting the node stiffness matrix from the extended overall stiffness matrix and then releasing the degrees of freedom of twist. A comparison with MSA using Euler–Bernoulli beam elements demonstrates the superiority of using Timoshenko beam elements. The 2PRU-UPR PM is presented to illustrate the effectiveness of the proposed approach. Finally, the global SEFI and scatter matrix are used to identify the elastic component with the weakest stiffness performance, providing a new approach for effectively improving the stiffness performance of the mechanism.     

Optimal maintenance scheduling for a complex manufacturing system subject to deterioration

We address the problem of determining inspection strategy and replacement policy for a deteriorating complex multi-component manufacturing system whose state is partially observable. We develop inspection and replacement scheduling models and other simple maintenance scheduling models via employing an imperfect repair model coupled with a damage process induced by operational conditions. The system state in performance of the imperfectly repaired system is modelled using a proportional intensity model incorporating a damage process and a virtual age process caused by repair. The system is monitored at periodic times and maintenance actions are carried out in response to the observed system state. Decisions to perform imperfect repair and replacement are based on the system state and crossing of a replacement threshold. The model proposed here aims at joint determination of a cost-optimal inspection and replacement policy along with an optimal level of maintenance which result in low maintenance cost and high operational performance and reliability of the system. To demonstrate the use of the model in practical applications a numerical example is provided. Solutions to optimal system parameters are obtained and the response of the model to these parameters is examined. Finally some features of the model are demonstrated. The approach presented provides a framework so that different scenario can be explored.     

Fuzzy-hybrid modelling of an Ackerman steered electric vehicle

Physical system modelling with known parameters together with 2-D or high order look-up tables (obtained from experimental data), have been the preferred method for simulating electric vehicles. The non-linear phenomena which are present at the vehicle tyre patch and ground interface have resulted in a quantitative understanding of this phenomena. However, nowadays, there is a requirement for a deeper understanding of the vehicle sub-models which previously used look-up tables. In this paper the hybrid modelling methodology used for electric vehicle systems offers a two-stage advantage: firstly, the vehicle model retains a comprehensive analytical formulation and secondly, the ‘fuzzy’ element offers, in addition to the quantitative results, a qualitative understanding of specific vehicle sub-models. In the literature several hybrid topologies are reported, sequential, auxiliary, and embedded.In this paper, the hybrid model topology selected is auxiliary and within the same hybrid model, the first paradigm used is the vehicle dynamics together with the actuator/gearbox system. The second paradigm is the non-linear fuzzy tyre model for each wheel. In particular, conventional physical system dynamic modelling has been combined with the fuzzy logic type-II or type-III methodology. The resulting hybrid-fuzzy tyre models were estimated for a-priori number of rules from experimental data. The physical system modelling required the available vehicle parameters such as the overall mass, wheel radius and chassis dimensions. The suggested synergetic fusion of the two methods, (hybrid-fuzzy), allowed the vehicle planar trajectories to be obtained prior to the hardware development of the entire vehicle. The strength of this methodology is that it requires localised system experimental data rather than global system data. The disadvantage in obtaining global experimental data is the requirement for comprehensive testing of a vehicle prototype which is both time consuming process and requires extensive resources. In this paper the authors have proposed the use of existing experimental rigs which are available from the leading automotive manufacturers. Hence, for the ‘hybrid’ modelling, localised data sets were used. In particular, wheel-tyre experimental data were obtained from the University tyre rig experimental facilities. Tyre forces acting on the tyre patch are mainly responsible for the overall electric vehicle motion. In addition, tyre measurement rigs are a well known method for obtaining localised data thus allowing the effective simulation of more detailed mathematical models. These include, firstly, physical system modelling (conventional vehicle dynamics), secondly, fuzzy type II or III modelling (for the tyre characteristics), and thirdly, electric drive modelling within the context of electric vehicles. The proposed hybrid model synthesis has resulted in simulation results which are similar to piece-wise ‘look-up’ table solutions. In addition, the strength of the ‘hybrid’ synthesis is that the analyst has a set of rules which clearly show the reasoning behind the complex development of the vehicle tyre forces. This is due to the inherent transparency of the type II and type III methodologies. Finally, the authors discussed the reasons for selecting a type-III framework. The paper concludes with a plethora of simulation results.     

Parametric solution to the joint system identification and optimization problem

This study is concerned with the simultaneous identification and optimization of static systems. The necessity and the advantages of an integrated approach to the identification and optimization of the system model is established theoretically as well as computationally. A parametric approach to the integrated problem is proven to converge to the integrated problem solution. The general methodology of decomposition of large-scale systems is extended by implementingfeasible decomposition of the joint problem. A multilevel approach is then utilized to successfully solve example problems. Handling the system constraints via a penalty-function technique is shown to be an efficient approach when using the parametric formulation of the joint problem. Numerical results for two example problems are presented using the Univac 1108 digital computer, revealing the economic advantages and disadvantages of the integrated approach to the identification and optimization problems.The authors wish to thank Messieurs I. Lefkowitz, L. S. Lasdon, F. Gembicki, and O. B. Olagundoye for their critique, comments, and suggestions during the course of this study.     

A continuous time markov-renewal replacement model for manpower systems

A continuous time Markov-renewal model is presented that generalizes the classical Young and Almond model for manpower systems with given size. The construction is based on the associated Markov-renewal replacement process and exploits the properties of the embedded replacement chain. The joint cumulant generating function of the grade sizes is derived and an asymptotic analysis provides conditions for these to converge in distribution to a multinominal random vector exponentially fast independently of the initial distribution, both for aperiodic and periodic embedded replacement chains. A regenerative approach to the wastage process is outlined and two numerical examples from the literature on manpower planning illustrate the theory.     

可持续的中国城乡居民养老保险体系数学模型

本文从我国现行养老金政策出发,利用灰关联分析得到养老金收支的影响因素的排序,建立中国城乡居民养老金收支模型.然后利用Logistic人口预测模型估计出不同替代率下2013-2035未来23年我国养老金缺口,并利用仿真方法算出保证我国养老保险体系的可持续性的替代率的合理区间为50%-70%.进一步分析影响替代率各因素的灵敏性,并考虑延迟退休、做实个人账户以及市场化投资等因素对模型进行改进.最后以养老金替代率为切入点,对退休年龄等模型参数的调整和有关措施进行了分析与评价.