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 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.

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.

The Energy Finite Element Analysis (EFEA) has been utilized successfully for modeling complex structural-acoustic systems with isotropic structural material properties. In this paper, a formulation for modeling structures made out of composite materials is presented. An approach based on spectral finite element analysis is utilized first for developing the equivalent material properties for the composite material. These equivalent properties are employed in the EFEA governing differential equations for representing the composite materials and deriving the element level matrices. The power transmission characteristics at connections between members made out of non-isotropic composite material are considered for deriving suitable power transmission coefficients at junctions of interconnected members. These coefficients are utilized for computing the joint matrix that is needed to assemble the global system of EFEA equations.

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 hybrid Finite Element Analysis (hybrid FEA) has been developed for performing structure-borne computations in automotive vehicle structures [1, 2 and 3]. The hybrid FEA method combines conventional FEA with Energy FEA (EFEA). Conventional FEA models are employed for modeling the behavior of the stiff members in a system. Appropriate damping and spring or mass elements are introduced in the connections between stiff and flexible members in order to capture the presence of the flexible members during the analyses of the stiff ones. The component mode synthesis method is combined with analytical solutions for determining the driving point conductance at joints between stiff and flexible members and for defining the properties of the concentrated elements which represent the flexible members when analyzing the stiff components.

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).

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.

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.

In order to develop confidence in numerical models which are used for automotive crash simulations, results are compared with test data. Modeling assumptions are made when constructing a simulation model for a complex system, such as a vehicle. Through a thorough understanding of the modeling assumptions an appropriate set of variables can be selected and adjusted in order to improve correlation with test data. Such a process can lead to better modeling practices when constructing a simulation model. Comparisons between the time history of acceleration responses from test and simulations are the most challenging. Computing accelerations correctly is more difficult compared to computing displacements, velocities, or intrusion levels due to the second order differentiation with time. In this paper a methodology for enabling the update of a simulation model for improved correlation is presented.

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.

The hybrid FEA method combines the conventional FEA method with the energy FEA (EFEA) for computing the structural vibration in vehicle structures when the excitation is applied on the load bearing stiff structural members. Conventional FEA models are employed for modeling the behavior of the stiff members in the vehicle. In order to account for the effect of the flexible members in the FEA analysis, appropriate damping and spring/mass elements are introduced at the connections between stiff and flexible members. Computing properly the values of these damping and spring/mass elements is important for the overall accuracy of the computations. Utilizing in these computations the analytical solutions for the driving point impedance of infinite or semi-infinite members introduces significant approximations.

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.

The Energy Finite Element Analysis (EFEA) has been developed for evaluating the vibro-acoustic behavior of complex systems. In the past EFEA results have been compared successfully to measured data for Naval, automotive, and aircraft systems. The main objective of this paper is to present information about the process of developing EFEA models for two configurations of a business jet, performing analysis for computing the vibration and the interior noise induced from exterior turbulent boundary layer excitation, and discussing the correlation between test data and simulation results. The structural EFEA model is generated from an existing finite element model used for stress analysis during the aircraft design process. Structural elements used in the finite element model for representing the complete complex aircraft structure become part of the EFEA structural model.

Probabilistic studies can be prohibitively expensive because they require repeated finite element analyses of large models. Re-analysis methods have been proposed with the premise to estimate accurately the dynamic response of a structure after a baseline design has been modified, without recalculating the new response. Although these methods increase computational efficiency, they are still not efficient enough for probabilistic analysis of large-scale dynamic systems with low failure probabilities (less or equal to 10-3). This paper presents a methodology that uses deterministic and probabilistic re-analysis methods to generate sample points of the response. Subsequently, tail modeling is used to estimate the right tail of the response PDF and the probability of failure a highly reliable system. The methodology is demonstrated on probabilistic vibration analysis of a realistic vehicle FE model.

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.