Technical overview of the equivalent static loads method for non-linear static response structural optimization

Linear static response structural optimization has been developed fairly well by using the finite element method for linear static analysis. However, development is extremely slow for structural optimization where a non linear static analysis technique is required. Optimization methods using equivalent static loads (ESLs) have been proposed to solve various structural optimization disciplines. The disciplines include linear dynamic response optimization, structural optimization for multi-body dynamic systems, structural optimization for flexible multi-body dynamic systems, nonlinear static response optimization and nonlinear dynamic response optimization. The ESL is defined as the static load that generates the same displacement field by an analysis which is not linear static. An analysis that is not linear static is carried out to evaluate the displacement field. ESLs are evaluated from the displacement field, linear static response optimization is performed by using the ESLs, and the design is updated. This process proceeds in a cyclic manner. A variety of problems have been solved by the ESLs methods. In this paper, the methods are completely overviewed. Various case studies are demonstrated and future research of the methods is discussed.     

An efficient procedure for structural reliability-based robust design optimization

This paper focuses on the development of an optimization tool with the aim to obtain robust and reliable designs in short computational time. The robustness measures considered here are the expected value and standard deviation of the performance function involved in the optimization problem. When using these robustness measures combined, the search of optimal design appears as a robust multiobjective optimization (RMO) problem. Reliable design addresses uncertainties to restrict the structural probability of failure. The mathematical formulation for the reliability based robust design optimization (RBRDO) problem is obtained by adding a reliability based constraint into the RMO problem. As both, statistics calculations and the reliability analysis could be very costly, approximation technique based on reduced-order modeling (ROM) is also incorporated in our procedure. The selected ROM is the proper orthogonal decomposition (POD) method, with the aim to produce fast outputs considering structural non-linear behavior. Moreover, to obtain RBRDO designs with reduced CPU time we propose others developments to be added in the integrated tool. They are: Probabilistic Collocation Method (PCM) to evaluate the statistics of the structural responses and, also, an approximated reliability constraints procedure based on the Performance Measure Approach (PMA) for reliability constraint assessment. Finally, Normal-Boundary Intersection (NBI) or Normalized Normal-Constraint (NNC) multiobjective optimization techniques are employed to obtain fast and even spread Pareto robust designs. To illustrate the application of the proposed tool, optimization studies are conducted for a linear (benchmark) and nonlinear trusses problems. The nonlinear example consider different loads level, exploring the material plasticity. The integrated tool prove to be very effective reducing the computational time by up to five orders of magnitude, when compared to the solutions obtained via classical standard approaches.     

‘Integral’ finite elements for designing the strong frames of highly maneuverable airplanes

The experience of using classical finite element programs in the computer-aided design of airframes is discussed. A set of features in the classical finite element method is mentioned which hinder the implementation of programs for automated structural optimization with respect to strength requirements. Principles simplifying implementation of such programs are proposed, including the principle of two-dimensional elements being identical to natural structural elements. The approximating functions are introduced on the basis of exact solutions for a beam and an anisotropic plate subjected to concentrated and distributed loads, these solutions being subsequently used to describe displacements of these elements. This approach is exemplified by test problems and is demonstrated in application to the highly maneuverable aircraft strong frame design.     

Elastoplastic truss design using a displacement based optimization

A displacement based optimization (DBO) method is applied to truss design problems with material nonlinearities, to explore feasibility and verify efficiency of the method as compared to traditional structural optimization. Minimum-weight truss sizing problems with various path-independent elastoplastic laws, including elastic perfectly plastic, linear strain hardening, and Ramberg–Osgood models, are investigated. This type of material nonlinearity allows us to naturally extend the linear elastic truss sizing in the DBO setting to nonlinear problems. To implement the methodology a computer program that uses the commercially available optimizer DOT by VR&D and IMSL Linear Programming solver by Visual Numerics is developed. Several test problems are successfully solved using the DBO approach and solutions are compared to those available in the literature, demonstrating significant reduction of computational time in comparison to the traditional structural optimization method. In particular, the DBO approach is found to be suitable for truss topology design since the method allows member areas to have cross-sectional areas equal to zero exactly.     

Limit state reliability optimization accounting for geometric effects

This paper proposes a unified formulation of reliability-based limit state optimization problems for discrete structures. Two approaches are presented: neglecting and accounting for nonlinear geometric effects at the limit state. For discrete structures, the classical limit state problem may be defined as a linear program, whereas the limit state problem accounting for geometric effects is considered as a two-level linear program. A portal frame example is considered in detail. Results of deterministic- and reliability-based optimization are compared. In both cases the impact on structural safety due to neglecting geometric effects is discussed.     

Limit state reliability optimization accounting for geometric effects

This paper proposes a unified formulation of reliability-based limit state optimization problems for discrete structures. Two approaches are presented: neglecting and accounting for nonlinear geometric effects at the limit state. For discrete structures, the classical limit state problem may be defined as a linear program, whereas the limit state problem accounting for geometric effects is considered as a two-level linear program. A portal frame example is considered in detail. Results of deterministic- and reliability-based optimization are compared. In both cases the impact on structural safety due to neglecting geometric effects is discussed.     

Optimum design of geometrically nonlinear space trusses

A structural optimization algorithm is developed for geometrically nonlinear three-dimensional trusses subject to displacement, stress and cross-sectional area constraints. The method is obtained by coupling the nonlinear analysis technique with the optimality criteria approach. The nonlinear behaviour of the space truss which was required for the steps of optimality criteria method was obtained by using iterative linear analysis. In each iteration the geometric stiffness matrix is constructed for the deformed structure and compensating load vector is applied to the system in order to adjust the joint displacements. During nonlinear analysis, tension members are loaded up to yield stress and compression members are stressed until their critical limits. The overall loss of elastic stability is checked throughout the steps of algorithm. The member forces resulted at the end of nonlinear analysis are used to obtain the new values of design variables for the next cycle. Number of design examples are presented to demonstrate the application of the algorithm. It is shown that the consideration of nonlinear behaviour of the space trusses in their optimum design makes it possible to achieve further reduction in the overall weight. The other advantage of the algorithm is that it takes into account the realistic behaviour of the structure, without which an optimum design might lead to erroneous result. This is noticed in one of the design example where a tension member changed into a compression one at the end of nonlinear analysis.     

Two highly efficient second-order algorithms for training feedforward networks

We present two highly efficient second-order algorithms for the training of multilayer feedforward neural networks. The algorithms are based on iterations of the form employed in the Levenberg-Marquardt (LM) method for nonlinear least squares problems with the inclusion of an additional adaptive momentum term arising from the formulation of the training task as a constrained optimization problem. Their implementation requires minimal additional computations compared to a standard LM iteration. Simulations of large scale classical neural-network benchmarks are presented which reveal the power of the two methods to obtain solutions in difficult problems, whereas other standard second-order techniques (including LM) fail to converge.     

Nonlinear dynamic response structural optimization using equivalent static loads

It is well known that nonlinear dynamic response optimization using a conventional optimization algorithm is fairly difficult and expensive for the gradient or non-gradient based optimization methods because many nonlinear dynamic analyses are required. Therefore, it is quite difficult to find practical large scale examples with many design variables and constraints for nonlinear dynamic response structural optimization. The equivalent static loads (ESLs) method is newly proposed and applied to nonlinear dynamic response optimization. The equivalent static loads are defined as the linear static load sets which generate the same response field in linear static analysis as that from nonlinear dynamic analysis. The ESLs are made from the results of nonlinear dynamic analysis and used as external forces in linear static response optimization. Then the same response from nonlinear dynamic analysis can be considered throughout linear static response optimization. The updated design from linear response optimization is used again in nonlinear dynamic analysis and the process proceeds in a cyclic manner until the convergence criteria are satisfied. Several examples are solved to validate the method. The results are compared to those of the conventional method with sensitivity analysis using the finite difference method.     

An efficient second-order SQP method for structural topology optimization

This article presents a Sequential Quadratic Programming (SQP) solver for structural topology optimization problems named TopSQP. The implementation is based on the general SQP method proposed in Morales et al. J Numer Anal 32(2):553–579 (2010) called SQP+. The topology optimization problem is modelled using a density approach and thus, is classified as a nonconvex problem. More specifically, the SQP method is designed for the classical minimum compliance problem with a constraint on the volume of the structure. The sub-problems are defined using second-order information. They are reformulated using the specific mathematical properties of the problem to significantly improve the efficiency of the solver. The performance of the TopSQP solver is compared to the special-purpose structural optimization method, the Globally Convergent Method of Moving Asymptotes (GCMMA) and the two general nonlinear solvers IPOPT and SNOPT. Numerical experiments on a large set of benchmark problems show good performance of TopSQP in terms of number of function evaluations. In addition, the use of second-order information helps to decrease the objective function value.     

Topology optimization of nonlinear structures with damping under arbitrary dynamic loading

This study presents an extended unit load method in which the displacement of a chosen degree of freedom (DOF) in a nonlinear structure under arbitrary dynamic loading is expressed as an integration of mutual strain energy density over a continuum domain. This new integral formulation for the displacement of a chosen DOF is developed by using the virtual work principle and can be used for linear or nonlinear structural behaviours. The integral form of the displacement is then used to develop new formulations for structural topology optimization involving arbitrary dynamic loading using the moving iso-surface threshold (MIST) method. Presented are two specific topology optimization problems with two objective functions: (a) to minimize the peak of a chosen displacement; or (b) to minimize the average power spectral density (PSD) of the chosen displacement over a finite time interval. New MIST formulations and algorithms are developed for solving two damping topology optimization problems of a structure under arbitrary dynamic loading, with or without large displacements, and having cellular damping materials with multi-volume fractions. Several numerical examples are presented to demonstrate the validity and efficiency of the presented unit load method and the MIST formulations and algorithms.     

Nonlinear dynamics of viscoelastic flexible structural systems by finite element method

This article presents a nonlinear displacement based finite elements model to study and analyze the nonlinear dynamic response of flexible double wishbone structural vehicle suspension system considering damping effect which was not previously discussed elsewhere. Due to large deflection and moderate rotation encountered during passing over road bumps, the kinematic nonlinearity is included through von Kármán strain component. Elastic undamped as well as viscous and viscoelastic damping mechanism are considered and compared. Considering the viscoelastic damping mechanism, the viscoelastic damping mechanism is modeled based on the integral constitutive form, which is recast into an incremental form suitable for finite element implementation. Additionally, the revolute joint element is adopted to incorporate the joint flexibility in the double wishbone system. The plane frame element is adopted to model the suspension links by using Timoshenko beam theory. The developed nonlinear finite element equations of motion are solved through the incremental iterative Newmark technique. The developed procedure is verified by comparing the obtained results with analytical solution and excellent agreement is observed. The applicability of the developed procedure is demonstrated by conducting parametric studies to show the effects of the road irregularities profiles, the vehicle speed, and the material damping coefficients on the nonlinear vibrations response of the double wishbone suspension systems. The obtained results are supportive in the design and manufacturing processes of these structural systems.

     

Mixed-integer linear programming approach for global discrete sizing optimization of frame structures

This paper focuses on discrete sizing optimization of frame structures using commercial profile catalogs. The optimization problem is formulated as a mixed-integer linear programming (MILP) problem by including the equations of structural analysis as constraints. The internal forces of the members are taken as continuous state variables. Binary variables are used for choosing the member profiles from a catalog. Both the displacement and stress constraints are formulated such that for each member limit values can be imposed at predefined locations along the member. A valuable feature of the formulation, lacking in most contemporary approaches, is that global optimality of the solution is guaranteed by solving the MILP using branch-and-bound techniques. The method is applied to three design problems: a portal frame, a two-story frame with three load cases and a multiple-bay multiple-story frame. Performance profiles are determined to compare the MILP reformulation method with a genetic algorithm.     

Dynamic response topology optimization in the time domain using model reduction method

The dynamic response topology optimization problems are usually computationally expensive, so it is necessary to employ the model reduction methods to reduce computational cost. This work will investigate the effectiveness of the mode displacement method(MDM) and mode acceleration method(MAM) for time-domain response problems within the framework of density-based topology optimization. Three objective functions, the mean dynamic compliance, mean strain energy and mean squared displacement are considered. It is found that, in general cases, MDM is not suitable for time-domain response topology optimization problems due to its low accuracy of approximation, while MAM works well for problems of a wide range, and when there are many time steps, the MAM based topology optimization approach is more efficient than the direct integration based approach. So for practical applications, when the problem needs many time steps, the MAM based approach is preferred and otherwise, the direct integration based approach is suggested.     

Some general results for optimal structures

This paper considers some approaches to structural optimization that meet real design requirements to a greater extent than has previously been achieved. In all of the problems considered, the structures are statically loaded, strains and displacements are small, and buckling effects are neglected.Topology optimization problems for multi-material (both linear and non-linear) elastic trusses, in particular, bi-material ones, are considered. Several general theorems are proven. Well-known classical results are generalized to minimal mass structures with restricted stresses or stiffness values. Some general theorems for trusses of given geometry are also proven. Moreover, it is proven that in the case of restricted stresses, the minimum mass non-linear-elastic truss of a given geometry is not heavier than the corresponding optimal linear-elastic one.     

On optimal physically nonlinear trusses

Optimization problems for physically nonlinear hyperelastic trusses statically loaded by a single system of loads are investigated. Several possible formulations of the problems for a prescribed truss layout are considered. It is demonstrated that the problems lead to the design of equal absolute stress values in some rods and of active technology constraints in other rods. Some general properties of the trusses are examined. Optimal physically nonlinear hyperelastic truss layout problems under static loading are also considered. It is shown that the well-known results of Maxwell, Michell and Prager concerning linear-elastic structures are generalized for the considered problems. It is demonstrated that the results of the paper are extended to optimization problems for structures in stationary creep with an arbitrary concave stress-strain velocity relation.     

Weakly and fully coupled methods for structural optimization of flexible mechanisms

The paper concerns a detailed comparison between two optimization methods that are used to perform the structural optimization of flexible components within a multibody system (MBS) simulation. The dynamic analysis of flexible MBS is based on a nonlinear finite element formulation. The first method is a weakly coupled method, which reformulates the dynamic response optimization problem in a two-level approach. First, a rigid or flexible MBS simulation is performed, and second, each component is optimized independently using a quasi-static approach in which a series of equivalent static load (ESL) cases obtained from the MBS simulation are applied to the respective components. The second method, the fully coupled method, performs the dynamic response optimization using the time response obtained directly from the flexible MBS simulation. Here, an original procedure is proposed to evaluate the ESL from a nonlinear finite element simulation, contrasting with the floating reference frame formulation exploited in the standard ESL method. Several numerical examples are provided to support our position. It is shown that the fully coupled method is more general and accommodates all types of constraints at the price of a more complex optimization process.