Search Results

Piston design is a challenging engineering problem which involves complex physics and requires satisfying multiple performance objectives. Uncertainty in piston operating conditions and variability in piston design variables are inevitable and must be accounted for. The piston assembly can be a major source of engine mechanical friction and cold start noise, if not designed properly. In this paper, an analytical piston model is used in a deterministic and probabilistic (reliability-based) multi-objective design optimization process to obtain an optimal piston design. The model predicts piston performance in terms of scuffing, friction and noise, In order to keep the computational cost low, efficient and accurate metamodels of the piston performance metrics are used. The Pareto set of all optimal solutions is calculated allowing the designer to choose the “best” solution according to trade-offs among the multiple objectives.

Using the total probability theorem, we propose a method to calculate the failure rate of a linear vibratory system with random parameters excited by stationary Gaussian processes. The response of such a system is non-stationary because of the randomness of the input parameters. A space-filling design, such as optimal symmetric Latin hypercube sampling or maximin, is first used to sample the input parameter space. For each design point, the output process is stationary and Gaussian. We present two approaches to calculate the corresponding conditional probability of failure. A Kriging metamodel is then created between the input parameters and the output conditional probabilities allowing us to estimate the conditional probabilities for any set of input parameters. The total probability theorem is finally applied to calculate the time-dependent probability of failure and the failure rate of the dynamic system. The proposed method is demonstrated using a vibratory system.

An approach for Probabilistic Reanalysis (PRA) of a system is presented. PRA calculates very efficiently the system reliability or the average value of an attribute of a design for many probability distributions of the input variables, by performing a single Monte Carlo simulation. In addition, PRA calculates the sensitivity derivatives of the reliability to the parameters of the probability distributions. The approach is useful for analysis problems where reliability bounds need to be calculated because the probability distribution of the input variables is uncertain or for design problems where the design variables are random. The accuracy and efficiency of PRA is demonstrated on vibration analysis of a car and on system reliability-based optimization (RBDO) of an internal combustion engine.

An efficient reliability estimation method is presented for engineering systems with multiple failure regions and potentially multiple most probable points. The method can handle implicit, nonlinear limit-state functions, with correlated or non-correlated random variables, which can be described by any probabilistic distribution. It uses a combination of approximate or “accurate-on-demand,” global and local metamodels which serve as indicators to determine the failure and safe regions. Samples close to limit states define transition regions between safe and failure domains. A clustering technique identifies all transition regions which can be in general disjoint, and local metamodels of the actual limit states are generated for each transition region.

Quality and performance are two important customer requirements in vehicle design. Driveline clunk negatively affects the perceived quality and must be therefore, minimized. This is usually achieved using engine torque management, which is part of engine calibration. During a tip-in event, the engine torque rate of rise is limited until all the driveline lash is taken up. However, the engine torque rise, and its rate can negatively affect the vehicle throttle response. Therefore, the engine torque management must be balanced against throttle response. In practice, the engine torque rate of rise is calibrated manually. This paper describes a methodology for calibrating the engine torque in order to minimize the clunk disturbance, while still meeting throttle response constraints. A set of predetermined engine torque profiles are calibrated in a vehicle and the transmission turbine speed is measured for each profile. The latter is used to quantify the clunk disturbance.

This paper explores the trade-off between reliability-based design and robustness for an automotive under-hood thermal system using the iSIGHT-FD environment. The interaction between the engine cooling system and the heating, ventilating, and air-conditioning (HVAC) system is described. The engine cooling system performance is modeled using Flowmaster and a metamodel is developed in iSIGHT. The actual HVAC system performance is characterized using test bench data. A design of experiment procedure determines the dominant factors and the statistics of the HVAC performance is obtained using Monte Carlo simulation (MCS). The MCS results are used to build an overall system response metamodel in order to reduce the computational effort. A multi-objective optimization in iSIGHT maximizes the system mean performance and simultaneously minimizes its standard deviation subject to probabilistic constraints.

An efficient approach for series system reliability-based design optimization (RBDO) is presented. The key idea is to apportion optimally the system reliability among the failure modes by considering the target values of the failure probabilities of the modes as design variables. Critical failure modes that contribute the most to the overall system reliability are identified. This paper proposes a computationally efficient, system RBDO approach using a single-loop method where the searches for the optimum design and for the most probable failure points proceed simultaneously. Specifically, at each iteration the optimizer uses approximated most probable failure points from the previous iteration to search for the optimum. A second-order Ditlevsen upper bound is used for the joint failure probability of failure modes. Also, an easy to implement active strategy set is employed to improve algorithmic stability.

Reliability-based design optimization accounts for variation. However, it assumes that statistical information is available in the form of fully defined probabilistic distributions. This is not true for a variety of engineering problems where uncertainty is usually given in terms of interval ranges. In this case, interval analysis or possibility theory can be used instead of probability theory. This paper shows how possibility theory can be used in design and presents a computationally efficient sequential optimization algorithm. The algorithm handles problems with only uncertain or a combination of random and uncertain design variables and parameters. It consists of a sequence of cycles composed of a deterministic design optimization followed by a set of worst-case reliability evaluation loops. A crank-slider mechanism example demonstrates the accuracy and efficiency of the proposed sequential algorithm.

Low vibration and noise level in internal combustion engines has become an essential part of the design process. It is well known that the piston assembly can be a major source of engine mechanical friction and cold start noise, if not designed properly. The piston secondary motion and piston-bore contact pattern are critical in piston design because they affect the skirt-to-bore impact force and therefore, how the piston impact excitation energy is damped, transmitted and eventually radiated from the engine structure as noise. An analytical method is presented in this paper for simulating piston secondary dynamics and piston-bore contact for an asymmetric half piston model. The method includes several important physical attributes such as bore distortion effects due to mechanical and thermal deformation, inertia loading, piston barrelity and ovality, piston flexibility and skirt-to-bore clearance. The method accounts for piston kinematics, rigid-body dynamics and flexibility.

An analytical method is presented in this paper for simulating piston secondary dynamics and piston-bore contact for an asymmetric half piston model including elastohydrodynamic (EHD) lubrication at the bore-skirt interface. A piston EHD analysis is used based on a finite-difference formulation. The oil film is discretized using a two-dimensional mesh. For improved computational efficiency without loss of accuracy, the Reynolds’ equation is solved using a perturbation approach which utilizes an “influence zone” concept, and a successive over-relaxation solver. The analysis includes several important physical attributes such as bore distortion effects due to mechanical and thermal deformation, inertia loading and piston barrelity and ovality. A Newmark-Beta time integration scheme combined with a Newton-Raphson linearization, calculates the piston secondary motion.

High-fidelity finite element (FE) tire-snow interaction models have the advantage of better understanding the physics of the tire-snow system. They can be used to develop semi-analytical models for vehicle design as well as to design and interpret field test results. For off-terrain conditions, there is a high level of uncertainties inherent in the system. The FE models are computationally intensive even when uncertainties of the system are not taken into account. On the other hand, field tests of tire-snow interaction are very costly. In this paper, dynamic metamodels are established to interpret interaction forces from FE simulation and to predict those forces by using part of the FE data as training data and part as validation data. Two metamodels are built based upon the Krieging principle: one has principal component analysis (PCA) taken into account and the other does not.

Styrene-Butadiene Rubber (SBR) is a copolymer of butadiene and styrene. It has a wide range of applications in the automotive industry due to its high durability, resistance to abrasion, oils and oxidation. SBR applications vary from tires to vibration isolators and gaskets. SBR is also used in tuned dampers which aim to reduce and control the angular vibrations of crankshafts, acting as an isolator and energy absorber between the tune damper's hub and the inertia ring. The dynamic properties of this polymer are therefore, very important in developing an appropriate analytical model. This paper presents the results of a series of experiments performed to determine the dynamic stiffness and damping properties of SBR. The frequency, temperature and displacement dependent properties are determined in a low frequency range from 0.4 to 150 Hz, and in a mid frequency range from 150 to 550 Hz. The most interesting property of SBR is its frequency dependent behavior.

In reliability design, often, there is scarce data for constructing probabilistic models. It is particularly challenging to model uncertainty in variables when the type of their probability distribution is unknown. Moreover, it is expensive to estimate the upper and lower bounds of the reliability of a system involving such variables. A method for modeling uncertainty by using Polynomial Chaos Expansion is presented. The method requires specifying bounds for statistical summaries such as the first four moments and credible intervals. A constrained optimization problem, in which decision variables are the coefficients of the Polynomial Chaos Expansion approximation, is formulated and solved in order to estimate the minimum and maximum values of a system’s reliability. This problem is solved efficiently by employing a probabilistic re-analysis approach to approximate the system reliability as a function of the moments of the random variables.

A complete probabilistic model of uncertainty in probabilistic analysis and design problems is the joint probability distribution of the random variables. Often, it is impractical to estimate this joint probability distribution because the mechanism of the dependence of the variables is not completely understood. This paper proposes modeling dependence by using copulas and demonstrates their representational power. It also compares this representation with a Monte-Carlo simulation using dispersive sampling.

An accurate and efficient Monte Carlo simulation (MCS) method is developed in this paper for limit state-based reliability analysis, especially at system levels, by using a response surface approximation of the failure indicator function. The Moving Least Squares (MLS) method is used to construct the response surface of the indicator function, along with an Optimum Symmetric Latin Hypercube (OSLH) as the sampling technique. Similar to MCS, the proposed method can easily handle implicit, highly nonlinear limit-state functions, with variables of any statistical distributions and correlations. However, the efficiency of MCS can be greatly improved. The method appears to be particularly efficient for multiple limit state and multiple design point problem. A mathematical example and a practical example are used to highlight the superior accuracy and efficiency of the proposed method over traditional reliability methods.

There are two sorts of uncertainty inherent in engineering design, the random and the epistemic. Random, or stochastic, uncertainty deals with the randomness or predictability of an event. It is well understood, easily modeled using classical probability, and ideal for such uncertainties as variations in manufacturing processes or material properties. Epistemic uncertainty deals with our lack of knowledge, our lack of information, and our own and others' subjectivity concerning design parameters. Epistemic uncertainty plays a particularly important role in the early stages of engineering design, when a lack of information about nominal values of parameters is much more important than potential variations in those parameters. Design reuse, or the design of product platforms, is an example in which epistemic uncertainty can play a crucial role in early design.

This paper presents the development of surrogate models (metamodels) for evaluating the bearing performance in an internal combustion engine. The metamodels are employed for performing probabilistic analyses for the engine bearings. The metamodels are developed based on results from a simulation solver computed at a limited number of sample points, which sample the design space. An integrated system-level engine simulation model, consisting of a flexible crankshaft dynamics model and a flexible engine block model connected by a detailed hydrodynamic lubrication model, is employed in this paper for generating information necessary to construct the metamodels. An optimal symmetric latin hypercube algorithm is utilized for identifying the sampling points based on the number and the range of the variables that are considered to vary in the design space.

A fast and accurate journal bearing elastohydrodynamic analysis is presented based on a finite difference formulation. The governing equations for the oil film pressure, stiffness and damping are solved using a finite difference approach. The oil film domain is discretized using a rectangular two-dimensional finite difference mesh. In this new formulation, it is not necessary to generate a global fluidity matrix similar to a finite element based solution. The finite difference equations are solved using a successive over relaxation (SOR) algorithm. The concept of “Influence Zone,” for computing the dynamic characteristics is introduced. The SOR algorithm and the “Influence Zone” concept significantly improve the computational efficiency without loss of accuracy. The new algorithms are validated with numerical results from the literature and their numerical efficiency is demonstrated.

A comprehensive formulation is presented for the dynamics of a rotating flexible crankshaft coupled with the dynamics of an engine block through a finite difference elastohydrodynamic main bearing lubrication algorithm. The coupling is based on detailed equilibrium conditions at the bearings. The component mode synthesis is employed for modeling the crankshaft and block dynamic behavior. A specialized algorithm for coupling the rigid and flexible body dynamics of the crankshaft within the framework of the component mode synthesis has been developed. A finite difference lubrication algorithm is used for computing the oil film elastohydrodynamic characteristics. A computationally accurate and efficient mapping algorithm has been developed for transferring information between a high - density computational grid for the elastohydrodynamic bearing solver and a low - density structural grid utilized in computing the crankshaft and block structural dynamic response.

This paper presents the development of surrogate models (metamodels) for evaluating the bearing performance in an internal combustion engine without performing time consuming analyses. The metamodels are developed based on results from actual simulation solvers computed at a limited number of sample points, which sample the design space. A finite difference bearing solver is employed in this paper for generating information necessary to construct the metamodels. An optimal symmetric Latin hypercube algorithm is utilized for identifying the sampling points based on the number and the range of the variables that are considered to vary in the design space. The development of the metamodels is validated by comparing results from the metamodels with results from the actual bearing performance solver over a large number of evaluation points. Once the metamodels are established they are employed for performing probabilistic analyses.