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.

It is challenging to perform probabilistic analysis and design of large-scale structures because probabilistic analysis requires repeated finite element analyses of large models and each analysis is expensive. This paper presents a methodology for probabilistic analysis and reliability based design optimization of large scale structures that consists of two re-analysis methods; one for estimating the deterministic vibratory response and another for estimating the probability of the response exceeding a certain level. The deterministic re-analysis method can analyze efficiently large-scale finite element models consisting of tens or hundreds of thousand degrees of freedom and large numbers of design variables that vary in a wide range. The probabilistic re-analysis method calculates very efficiently the system reliability for many probability distributions of the design variables by performing a single Monte Carlo simulation.

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.

Early in the engineering design cycle, it is difficult to quantify product reliability due to insufficient data or information to model uncertainties. Probability theory can not be therefore, used. Design decisions are usually based on fuzzy information which is imprecise and incomplete. Various design methods such as Possibility-Based Design Optimization (PBDO) and Evidence-Based Design Optimization (EBDO) have been developed to systematically treat design with non-probabilistic uncertainties. In practical engineering applications, information regarding the uncertain variables and parameters may exist in the form of sample points, and uncertainties with sufficient and insufficient information may exist simultaneously. Most of the existing optimal design methods under uncertainty can not handle this form of incomplete information. They have to either discard some valuable information or postulate the existence of additional information.

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.

In this paper the development of statistical metamodels and statistical fast running models is presented first. They are utilized for propagating uncertainties in a multi-discipline design optimization process. Two main types of uncertainty can be considered in this manner: uncertainty due to variability in design variables or in random parameters; uncertainty due to the utilization of metamodels instead of the actual simulation models during the optimization process. The value of the new developments and their engagement in multi-discipline design optimization is demonstrated through a case study. An underwater vehicle is designed under four different disciplines, namely, noise radiation, self-noise due to TBL excitation, dynamic response due to propulsion impact loads, and response to an underwater detonation.

The cathode gas diffusion layer (GDL) is an important component of polymer electrolyte membrane (PEM) fuel cell. Its design parameters, including thickness, porosity and permeability, significantly affect the reactant transport and water management, thus impacting the fuel cell performance. This paper presents an optimization study of the GDL design parameters with the objective of maximizing the current density under a given voltage. A two-dimensional single-phase PEM fuel cell model is used. A multivariable optimization problem is formed to maximize the current density at the cathode under a given electrode voltage with respect to the GDL parameters. In order to reduce the computational effort and find the global optimum among the potential multiple optima, a global metamodel of the actual CFD-based fuel cell simulation, is adaptively generated using radial basis function approximations.

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 presents a validation case study for an Energy Finite Element Analysis (EFEA) formulation through comparison to test data. The EFEA comprises a simulation tool for computing the structural response of a complex structure and the amount of the radiated power. The EFEA formulation presented in this paper can account for periodic stiffeners, for partial fluid loading effects on the outer part of the structure, and for internal compartments filled with heavy fluid. In order to validate these modeling capabilities of the EFEA two 1/8th scale structures representing an advanced double hull design and a conventional hull design of a surface ship are analyzed. Results for the structural vibration induced on the outer bottom part of the structure are compared to available test data. The excitation is applied at two different locations of the deck structure. Good correlation is observed between the numerical results and the test data.

In applications of the Energy Finite Element Analysis (EFEA) there is an increasing need for developing comprehensive models with a large number of elements which include both structural and interior fluid elements, while certain parts of the structure are considered to be exposed to an external fluid loading. In order to accommodate efficient computations when using simulation models with a large number of elements, joints, and domains, a substructuring computational capability has been developed. The new algorithm is based on dividing the EFEA model into substructures with internal and interface degrees of freedom. The system of equations for each substructure is assembled and solved separately and the information is condensed to the interface degrees of freedom. The condensed systems of equations from each substructure are assembled in a reduced global system of equations. Once the global system of equations has been solved the solution for each substructure is pursued.

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.

The Energy Boundary Element Analysis (EBEA) has been utilized in the past for computing the exterior acoustic field at high frequencies (above ∼400Hz) around vehicle structures and numerical results have been compared successfully to test data [1, 2 and 3]. The Energy Finite Element Analysis (EFEA) has been developed for computing the structural vibration of complex structures at high frequencies and validations have been presented in previous publications [4, 5]. In this paper the EBEA is utilized for computing the acoustic field around a vehicle structure due to external acoustic noise sources. The computed exterior acoustic field comprises the excitation for the EFEA analysis. Appropriate loading functions have been developed for representing the exterior acoustic loading in the EFEA simulations, and a formulation has been developed for considering the acoustic treatment applied on the interior side of structural panels.

A reliability analysis method is presented for time-dependent systems under uncertainty. A level-crossing problem is considered where the system fails if its maximum response exceeds a specified threshold. The proposed method uses a double-loop optimization algorithm. The inner loop calculates the maximum response in time for a given set of random variables, and transforms a time-dependent problem into a time-independent one. A time integration method is used to calculate the response at discrete times. For each sample function of the response random process, the maximum response is found using a global-local search method consisting of a genetic algorithm (GA), and a gradient-based optimizer. This dynamic response usually exhibits multiple peaks and crosses the allowable response level to form a set of complex limit states, which lead to multiple most probable points (MPPs).

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.

Utilizing simulation technology is important for designing a structure with increased survivability to a load from an explosion. The pressure wave from the blast and the fragments hitting the structure must be simulated in such an analysis. Commercial software can be utilized through the development of appropriate interfaces for performing such computations. In this paper an approach is presented for combining commercially available Eulerian and Lagrangian solvers for performing blast event simulations. A capability has been developed for automatically creating the Eulerian finite element given the finite element model for the structure. The effect of moisture in the soil properties is considered during the generation of the soil - explosive - air model used by the Eulerian solver. Tracers are defined in the Eulerian model for all structural finite elements which are on the outer part of the structure and are subjected to the load from the blast.