Search Results

A Reliability-Based Design Optimization (RBDO) method for multiple failure regions is presented. The method uses a Probabilistic Re-Analysis (PRRA) approach in conjunction with an approximate global metamodel with local refinements. The latter serves as an indicator to determine the failure and safe regions. PRRA calculates very efficiently the system reliability of a design by performing a single Monte Carlo (MC) simulation. Although PRRA is based on MC simulation, it calculates “smooth” sensitivity derivatives, allowing therefore, the use of a gradient-based optimizer. An “accurate-on-demand” metamodel is used in the PRRA that allows us to handle problems with multiple disjoint failure regions and potentially multiple most-probable points (MPP). The multiple failure regions are identified by using a clustering technique. A maximin “space-filling” sampling technique is used to construct the metamodel. A vibration absorber example highlights the potential of the proposed method.

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.

Finite element analysis is a well-established methodology in structural dynamics. However, optimization and/or probabilistic studies can be prohibitively expensive because they require repeated FE analyses of large models. Various reanalysis methods have been proposed in order to calculate efficiently the dynamic response of a structure after a baseline design has been modified, without recalculating the new response. The parametric reduced-order modeling (PROM) and the combined approximation (CA) methods are two re-analysis methods, which can handle large model parameter changes in a relatively efficient manner. Although both methods are promising by themselves, they can not handle large FE models with large numbers of DOF (e.g. 100,000) with a large number of design parameters (e.g. 50), which are common in practice. In this paper, the advantages and disadvantages of the PROM and CA methods are first discussed in detail.

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.

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.

The door closing effort is a quality issue concerning both automobile designers and customers. This paper describes an Excel based mathematical model for predicting the side door closing effort in terms of the required minimum energy or velocity, to close the door from a small open position when the check-link ceases to function. A simplified but comprehensive model is developed which includes the cabin pressure (air bind), seal compression, door weight, latch effort, and hinge friction effects. The flexibility of the door and car body is ignored. Because the model simplification introduces errors, we calibrate it using measured data. Calibration is also necessary because some input parameters are difficult to obtain directly. In this work, we provide the option to calibrate the hinge model, the latch model, the seal compression model, and the air bind model. The door weight effect is geometrically exact, and does not need calibration.