Reliability-based shape optimization of structures undergoing fluid–structure interaction phenomena

Fluid–structure interaction phenomena are often roughly approximated when the stochastic nature of a system is considered in the design optimization process, leading to potentially significant epistemic uncertainty. In this paper, after reviewing the state-of-the-art methods in robust and reliability-based design optimization of problems undergoing fluid–structure interaction phenomena, a computational framework is presented that integrates a high-fidelity aeroelastic model into reliability-based design optimization. The design optimization problem is formulated pursuant to the reliability index and performance measure approaches. The system reliability is evaluated by a first-order reliability analysis method. The steady-state aeroelastic problem is described by a three-field formulation and solved by a staggered procedure, coupling a potentially detailed structural finite element model and a finite volume discretization of the Euler flow. The design and imperfection sensitivities are computed by evaluating the analytically derived direct and adjoint coupled aeroelastic sensitivity equations. The computational framework is verified by the optimization of three-dimensional wing structures. The lift-to-drag ratio is maximized, subject to stress constraints, by varying shape, thickness, and material properties. Uncertainties in structural parameters, including design parameters, operating conditions, and modeling uncertainties are considered. The results demonstrate the need for reliability-based optimization methods, for the design of structures undergoing fluid–structure interaction phenomena, and the applicability of the proposed framework to realistic design problems. Comparing the optimization results for different levels of uncertainty shows the importance of accounting for uncertainties in a quantitative manner.     

Efficient strategies for reliability-based design optimization of variable stiffness composite structures

This study investigates efficient design optimization frameworks for composite structures with uncertainties related to material properties and loading. The integration of two decoupled reliability-based design optimization methodologies with a decoupled discrete material optimization is proposed to determine material and fiber orientation for three-dimensional composite structures. First, a deterministic and decoupled discrete material optimization is used for baseline comparison. The objective is to minimize the cost of composite structures with the design variables comprising of the piecewise patch orientations and material properties of the fiber reinforced composites. The reliability-based design optimization includes a hybrid method, and also the sequential optimization and reliability assessment method. In the sequential optimization and reliability assessment method, the inverse reliability analysis is evaluated using a stochastic response surface method and a first order reliability approach. Comparing the methods based on the optimal material and fiber orientations, the uncertainties in loads and material properties lead to different optimal layouts compared to the deterministic solutions. The numerical results also reveal that the hybrid method applied in reliability based designs results in negligible additional computational cost.     

Reliability-based structural optimization of frame structures for multiple failure criteria using topology optimization techniques

Topology optimization methods using discrete elements such as frame elements can provide useful insights into the underlying mechanics principles of products; however, the majority of such optimizations are performed under deterministic conditions. To avoid performance reductions due to later-stage environmental changes, variations of several design parameters are considered during the topology optimization. This paper concerns a reliability-based topology optimization method for frame structures that considers uncertainties in applied loads and nonstructural mass at the early conceptual design stage. The effects that multiple criteria, namely, stiffness and eigenfrequency, have upon system reliability are evaluated by regarding them as a series system, where mode reliabilities can be evaluated using first-order reliability methods. Through numerical calculations, reliability-based topology designs of typical two- or three-dimensional frames are obtained. The importance of considering uncertainties is then demonstrated by comparing the results obtained by the proposed method with deterministic optimal designs.     

Computational strategies for reliability-based structural optimization of aeroelastic limit cycle oscillations

Gradient-based optimization, via the adjoint method, is needed to realistically enable the reliability-based design of a nonlinear unsteady aeroelastic system with many random and/or deterministic design variables. The adjoint derivatives of a time-marched system entail a cumbersome reverse-time integration, and so a time-periodic spectral element scheme is used here to efficiently capture the gradients of the limit cycle oscillations. Further reductions in the computational cost of the monolithic-time adjoint vector are obtained with proper orthogonal decomposition, which projects the large system onto a reduced basis. Design reliability is computed with the first order reliability method, which provides an estimate of the failure probability without resorting to sampling-based approaches (infeasible for large systems). Analytical gradients are needed to obtain the most probable point (in the random variable space), as well as the reliability design derivatives. These computational strategies are utilized to locate the optimal thickness distribution of a cantilevered wing operating beyond its flutter point in supersonic flow (via piston theory). Specifically, the wing mass is minimized under both deterministic and non-deterministic limit cycle oscillation amplitude constraints, with both structural and flow uncertainties considered in the latter.     

Conceptual design of aeroelastic structures by topology optimization

A topology optimization methodology is presented for the conceptual design of aeroelastic structures accounting for the fluid–structure interaction. The geometrical layout of the internal structure, such as the layout of stiffeners in a wing, is optimized by material topology optimization. The topology of the wet surface, that is, the fluid–structure interface, is not varied. The key components of the proposed methodology are a Sequential Augmented Lagrangian method for solving the resulting large-scale parameter optimization problem, a staggered procedure for computing the steady-state solution of the underlying nonlinear aeroelastic analysis problem, and an analytical adjoint method for evaluating the coupled aeroelastic sensitivities. The fluid–structure interaction problem is modeled by a three-field formulation that couples the structural displacements, the flow field, and the motion of the fluid mesh. The structural response is simulated by a three-dimensional finite element method, and the aerodynamic loads are predicted by a three-dimensional finite volume discretization of a nonlinear Euler flow. The proposed methodology is illustrated by the conceptual design of wing structures. The optimization results show the significant influence of the design dependency of the loads on the optimal layout of flexible structures when compared with results that assume a constant aerodynamic load.     

A reliability-based multidisciplinary design optimization procedure based on combined probability and evidence theory

To address the reliability-based multidisciplinary design optimization (RBMDO) problem under mixed aleatory and epistemic uncertainties, an RBMDO procedure is proposed in this paper based on combined probability and evidence theory. The existing deterministic multistage-multilevel multidisciplinary design optimization (MDO) procedure MDF-CSSO, which combines the multiple discipline feasible (MDF) procedure and the concurrent subspace optimization (CSSO) procedure to mimic the general conceptual design process, is used as the basic framework. In the first stage, the surrogate based MDF is used to quickly identify the promising reliable regions. In the second stage, the surrogate based CSSO is used to organize the disciplinary optimization and system coordination, which allows the disciplinary specialists to investigate and optimize the design with the corresponding high-fidelity models independently and concurrently. In these two stages, the reliability-based optimization both in the system level and the disciplinary level are computationally expensive as it entails nested optimization and uncertainty analysis. To alleviate the computational burden, the sequential optimization and mixed uncertainty analysis (SOMUA) method is used to decompose the traditional double-level reliability-based optimization problem into separate deterministic optimization and mixed uncertainty analysis sub-problems, which are solved sequentially and iteratively until convergence is achieved. By integrating SOMUA into MDF-CSSO, the Mixed Uncertainty based RBMDO procedure MUMDF-CSSO is developed. The effectiveness of the proposed procedure is testified with one simple numerical example and one MDO benchmark test problem, followed by some conclusion remarks.     

Component and system reliability-based topology optimization using a single-loop method

We perform reliability-based topology optimization by combining reliability analysis and material distribution topology design methods to design linear elastic structures subject to random inputs, such as random loadings. Both component reliability and system reliability are considered. In component reliability, we satisfy numerous probabilistic constraints which quantify the failure of different events. In system reliability, we satisfy a single probabilistic constraint which encompasses the component events. We adopt the first-order reliability method to approximate the component reliabilities and the inclusion-exclusion rule to approximate the system reliability. To solve the probabilistic optimization problem, we use a variant of the single loop method, which eliminates the need for an inner reliability analysis loop. The proposed method is amenable to implementation with existing deterministic topology optimization software, and hence suitable for practical applications. Designs obtained from component and system reliability-based topology optimization are compared to those obtained from traditional deterministic topology optimization and validated via Monte Carlo simulation.     

Optimum values of structural safety factors for a predefined reliability level with extension to multiple limit states

In the field of deterministic structural optimization, the designer reduces the structural cost without taking into account uncertainties concerning materials, geometry and loading. This way, the resulting optimum solution may represent a lower level of reliability and thus a higher risk of failure. It is the objective of reliability-based design optimization (RBDO) to design structures that should be both economic and reliable. The coupling between mechanical modeling, reliability analyses and optimization methods leads to very high computational costs and weak convergence stability. Since the traditional RBDO solution is achieved by alternating between reliability and optimization iterations, the structural designers performing deterministic optimization do not consider the RBDO model as a practical tool for the design of real structures. Fortunately, a hybrid method based on simultaneous solution of the reliability and the optimization problem, has successfully reduced the computational time problem. The hybrid method allows us to satisfy a required reliability level, but the vector of variables here contains both deterministic and random variables. The hybrid RBDO problem is thus more complex than that of deterministic design. The major difficulty lies in the evaluation of the structural reliability, which is carried out by a special optimization procedure. In this paper a new methodology is presented with the aim of finding a global solution to RBDO problems without additional computing cost for the reliability evaluation. The safety factor formulation for a single limit state case has been used to efficiently reduce the computational time . This technique is fundamentally based on a study of the sensitivity of the limit state function with respect to the design variables. In order to demonstrate analytically the efficiency of this methodology, the optimality condition is then used. The efficiency of this technique is also extended to multiple limit state cases. Two numerical examples are presented at the end of the paper to demonstrate the applicability of the new methodology.     

Aeroelastic tailoring considering uncertainties in material properties

The robustness of aeroelastic design optimization with respect to uncertainties in material and structural properties is studied both numerically and experimentally. The model consists of thin orthotropic composite wings virtually without fuselage. Three different configurations with consistent geometry but varying orientation of the main stiffness axis of the material are investigated. The onset of aeroelastic instability, flutter, is predicted using finite element analysis and the doublet-lattice method for the unsteady aerodynamic forces. The numerical results are experimentally verified in a low-speed wind tunnel. The optimization problem is stated as to increase the critical air speed, above that of the bare wing by massbalancing. It is seen that the design goals are not met in the experiments due to uncertainties in the structural performance of the wings. The uncertainty in structural performance is quantified through numerous dynamic material tests. Once accounting for the uncertainties through a suggested reformulation of the optimization problem, the design goals are met also in practice. The investigation indicates that robust and reliable aeroelastic design optimization is achievable, but careful formulation of the optimization problem is essential.     

Sequential optimization and reliability assessment for multidisciplinary systems design

With higher reliability and safety requirements, reliability-based design has been increasingly applied in multidisciplinary design optimization (MDO). A direct integration of reliability-based design and MDO may present tremendous implementation and numerical difficulties. In this work, a methodology of sequential optimization and reliability assessment for MDO is proposed to improve the efficiency of reliability-based MDO. The central idea is to decouple the reliability analysis from MDO with sequential cycles of reliability analysis and deterministic MDO. The reliability analysis is based on the first-order reliability method (FORM). In the proposed method, the reliability analysis and the deterministic MDO use two MDO strategies, the multidisciplinary feasible approach and the individual disciplinary feasible approach. The effectiveness of the proposed method is illustrated with two example problems.     

Modeling physical uncertainties in dynamic stall induced fluid–structure interaction of turbine blades using arbitrary polynomial chaos

A nonlinear dynamic problem of stall induced flutter oscillation subject to physical uncertainties is analyzed using arbitrary polynomial chaos. A single-degree-of-freedom stall flutter model with torsional oscillation is considered subject to nonlinear aerodynamic loads in the dynamic stall regime and nonlinear structural stiffness. The analysis of the deterministic aeroelastic response demonstrated that the problem is sensitive to variations in structural natural frequency and structural nonlinearity. The effect of uncertainties in these parameters is studied. Arbitrary polynomial chaos is employed in which appropriate expansion polynomials are constructed based on the statistical moments of the uncertain input. The arbitrary polynomial chaos results are compared with Monte Carlo simulations.     

Design of aircraft wings subjected to gust loads: A system reliability approach

A method for system reliability-based design of aircraft wing structures is presented. A wing of a light commuter aircraft designed for gust loads according to the FAA regulations is compared with one designed by system reliability optimization. It is shown that system reliability optimization has the potential of improving dramatically the safety and efficiency of new designs. The reasons for the differences between the deterministic and reliability-based designs are explained.     

Reliability-based robust assessment for multiobjective optimization design of improving occupant restraint system performance

Optimal performance of vehicle occupant restraint system (ORS) requires an accurate assessment of occupant injury values including head, neck and chest responses, etc. To provide a feasible framework for incorporating occupant injury characteristics into the ORS design schemes, this paper presents a reliability-based robust approach for the development of the ORS. The uncertainties of design variables are addressed and the general formulations of reliable and robust design are given in the optimization process. The ORS optimization is a highly nonlinear and large scale problem. In order to save the computational cost, an optimal sampling strategy is applied to generate sample points at the stage of design of experiment (DOE). Further, to efficiently obtain a robust approximation, the support vector regression (SVR) is suggested to construct the surrogate model in the vehicle ORS design process. The multiobjective particle swarm optimization (MPSO) algorithm is used for obtaining the Pareto optimal set with emphasis on resolving conflicting requirements from some of the objectives and the Monte Carlo simulation (MCS) method is applied to perform the reliability and robustness analysis. The differences of three different Pareto fronts of the deterministic, reliable and robust multiobjective optimization designs are compared and analyzed in this study. Finally, the reliability-based robust optimization result is verified by using sled system test. The result shows that the proposed reliability-based robust optimization design is efficient in solving ORS design optimization problems.     

An efficient strategy for reliability-based multidisciplinary design optimization using BLISS

This paper presents an efficient reliability-based multidisciplinary design optimization (RBMDO) strategy. The conventional RBMDO has tri-level loops: the first level is an optimization in the deterministic space, the second one is a reliability analysis in the probabilistic space, and the third one is the multidisciplinary analysis. Since it is computationally inefficient when high-fidelity simulation methods are involved, an efficient strategy is proposed. The strategy [named probabilistic bi-level integrated system synthesis (ProBLISS)] utilizes a single-level reliability-based design optimization (RBDO) approach, in which the reliability analysis and optimization are conducted in a sequential manner by approximating limit state functions. The single-level RBDO is associated with the BLISS formulation to solve RBMDO problems. Since both the single-level RBDO and BLISS are mainly driven by approximate models, the accuracy of models can be a critical issue for convergence. The convergence of the strategy is guaranteed by employing the trust region–sequential quadratic programming framework, which validates approximation models in the trust region radius. Two multidisciplinary problems are tested to verify the strategy. ProBLISS significantly reduces the computational cost and shows stable convergence while maintaining accuracy.     

A stability investigation of a simulation- and reliability-based optimization

This study developed a reliability-based design optimization (RBDO) algorithm focusing on the ability of solving problems with nonlinear constraints or system reliability. In this case, a sampling technique is often adopted to evaluate the reliability analyses. However, simulation with an insufficient sample size often possesses statistical randomness resulting in an inaccurate sensitivity calculation. This may cause an unstable RBDO solution. The proposed approach used a set of deterministic variables, called auxiliary design points, to replace the random parameters. Thus, an RBDO is converted into a deterministic optimization (DO, α-problem). The DO and the analysis of finding the auxiliary design points (β-problem) are conducted iteratively until the solution converges. To maintain the stability of the RBDO solution with less computational cost, the proposed approach calculated the sensitivity of reliability (in the β-problem) with respect to the mean value of the pseudo-random parameters rather than the design variables. The stability of the proposed method was compared to that of the double-loop approach, and many factors, such as sample size, starting point and the parameters used in the optimization, were considered. The accuracy of the proposed method was confirmed using Monte Carlo simulation (MCS) with several linear and nonlinear numerical problems.     

Continuum topology optimization with non-probabilistic reliability constraints based on multi-ellipsoid convex model

Using a quantified measure for non-probab ilistic reliability based on the multi-ellipsoid convex model, the topology optimization of continuum structures in presence of uncertain-but-bounded parameters is investigated. The problem is formulated as a double-loop optimization one. The inner loop handles evaluation of the non-probabilistic reliability index, and the outer loop treats the optimum material distribution using the results from the inner loop for checking feasibility of the reliability constraints. For circumventing the numerical difficulties arising from its nested nature, the topology optimization problem with reliability constraints is reformulated into an equivalent one with constraints on the concerned performance. In this context, the adjoint variable schemes for sensitivity analysis with respect to uncertain variables as well as design variables are discussed. The structural optimization problem is then solved by a gradient-based algorithm using the obtained sensitivity. In the present formulation, the uncertain-but bounded uncertain variations of material properties, geometrical dimensions and loading conditions can be realistically accounted for. Numerical investigations illustrate the applicability and the validity of the present problem statement as well as the proposed numerical techniques. The computational results also reveal that non-probabilistic reliability-based topology optimization may yield more reasonable material layouts than conventional deterministic approaches. The proposed method can be regarded as an attractive supplement to the stochastic reliability-based topology optimization.     

Constrained multifidelity optimization using model calibration

Multifidelity optimization approaches seek to bring higher-fidelity analyses earlier into the design process by using performance estimates from lower-fidelity models to accelerate convergence towards the optimum of a high-fidelity design problem. Current multifidelity optimization methods generally fall into two broad categories: provably convergent methods that use either the high-fidelity gradient or a high-fidelity pattern-search, and heuristic model calibration approaches, such as interpolating high-fidelity data or adding a Kriging error model to a lower-fidelity function. This paper presents a multifidelity optimization method that bridges these two ideas; our method iteratively calibrates lower-fidelity information to the high-fidelity function in order to find an optimum of the high-fidelity design problem. The algorithm developed minimizes a high-fidelity objective function subject to a high-fidelity constraint and other simple constraints. The algorithm never computes the gradient of a high-fidelity function; however, it achieves first-order optimality using sensitivity information from the calibrated low-fidelity models, which are constructed to have negligible error in a neighborhood around the solution. The method is demonstrated for aerodynamic shape optimization and shows at least an 80% reduction in the number of high-fidelity analyses compared other single-fidelity derivative-free and sequential quadratic programming methods. The method uses approximately the same number of high-fidelity analyses as a multifidelity trust-region algorithm that estimates the high-fidelity gradient using finite differences.     

Non-probabilistic reliability-based topology optimization with multidimensional parallelepiped convex model

In this paper, a new non-probabilistic reliability-based topology optimization (NRBTO) method is proposed to account for interval uncertainties considering parametric correlations. Firstly, a reliability index is defined based on a newly developed multidimensional parallelepiped (MP) convex model, and the reliability-based topology optimization problem is formulated to optimize the topology of the structure, to minimize material volume under displacement constraints. Secondly, an efficient decoupling scheme is applied to transform the double-loop NRBTO into a sequential optimization process, using the sequential optimization & reliability assessment (SORA) method associated with the performance measurement approach (PMA). Thirdly, the adjoint variable method is used to obtain the sensitivity information for both uncertain and design variables, and a gradient-based algorithm is employed to solve the optimization problem. Finally, typical numerical examples are used to demonstrate the effectiveness of the proposed topology optimization method.     

Reliability-Based Optimization Using Evolutionary Algorithms

Uncertainties in design variables and problem parameters are often inevitable and must be considered in an optimization task if reliable optimal solutions are sought. Besides a number of sampling techniques, there exist several mathematical approximations of a solution's reliability. These techniques are coupled in various ways with optimization in the classical reliability-based optimization field. This paper demonstrates how classical reliability-based concepts can be borrowed and modified and, with integrated single and multiobjective evolutionary algorithms, used to enhance their scope in handling uncertainties involved among decision variables and problem parameters. Three different optimization tasks are discussed in which classical reliability-based optimization procedures usually have difficulties, namely (1) reliability-based optimization problems having multiple local optima, (2) finding and revealing reliable solutions for different reliability indices simultaneously by means of a bi-criterion optimization approach, and (3) multiobjective optimization with uncertainty and specified system or component reliability values. Each of these optimization tasks is illustrated by solving a number of test problems and a well-studied automobile design problem. Results are also compared with a classical reliability-based methodology.     

Structural dynamic topology optimization based on dynamic reliability using equivalent static loads

An approach for reliability-based topology optimization of interval parameters structures under dynamic loads is proposed. We modify the equivalent static loads method for non linear static response structural optimization (ESLSO) to solve the dynamic reliability optimization problem. In our modified ESLSO, the equivalent static loads (ESLs) are redefined to consider the uncertainties. The new ESLs including all the uncertainties from geometric dimensions, material properties and loading conditions generate the same interval response field as dynamic loads. Based on the definition of the interval non-probabilistic reliability index, we construct the static reliability topology optimization model using ESLs. The method of moving asymptotes (MMA) is employed as the optimization problem solver. The applicability and validity of the proposed model and numerical techniques are demonstrated with three numerical examples.