Optimization of geometrically non‐linear symmetric systems with coincident critical points

A new formulation is presented for optimum design of an elastic symmetric structure for specified non‐linear buckling load factor. It is shown that the method of sensitivity analysis of bifurcation load factor developed in Ohsaki and Vetani (Int. J. Numer. Methods Engng. 1996; 39 : 1707–1720) can also be applied for the case where the structure reaches coincident critical points including a limit point. Based on the method of sensitivity analysis, an algorithm is presented for finding optimum designs for specified coincident critical points. The well‐known danger of designing a structure that exhibits coincident buckling is discussed in detail. It is shown in the examples of trusses that the structural volume may be successfully reduced as a result of optimization even if the reduction of the maximum load factor due to possible asymmetric initial imperfection is considered. Copyright © 2000 John Wiley & Sons, Ltd.     

Mesh‐free radial point interpolation method for the buckling analysis of Mindlin plates subjected to in‐plane point loads

In this paper, a mesh‐free approach is employed for buckling analysis of Mindlin plates that are subjected to in‐plane point loads. The radial point interpolation method (RPIM) is used to approximate displacements based on nodes. Variational forms of the system equations are established. Two‐step solution procedures are implemented. The non‐uniform pre‐stress distribution of plate is first obtained using the RPIM based on a two‐dimensional (2D) elastic plane stress problem. This predetermined non‐uniform pre‐stress distribution is then used to compute buckling loads of plate using the RPIM based on Mindlin's plate assumption. The RPIM can easily handle any number and location of nodes in the plate domain for a desired computational accuracy without major difficulties in solving the initial stresses and buckling loads. Numerical examples considered here include circular and rectangular Mindlin plates that are subjected to in‐plane uniform and point loads with different aspect ratios and boundary conditions. The present results are validated against the available analytical and numerical solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Efficient non‐linear reanalysis of skeletal structures using combined approximations

Non‐linear reanalysis of large‐scale structures usually involves much computational effort, because the set of non‐linear equations must be solved repeatedly during the solution process. Various approximations that are often used for linear reanalysis are not sufficiently accurate for non‐linear problems. In this study, solution procedures based on the combined approximations approach are developed and compared in terms of efficiency and accuracy. Various path‐independent non‐linear analysis and reanalysis problems are considered, including material non‐linearity, geometric non‐linearity and buckling analysis. Numerical examples demonstrate the effectiveness of the procedures presented. It is shown that in various cases accurate results can be achieved efficiently. Copyright © 2007 John Wiley & Sons, Ltd.     

Optimum design of cross-sectional profiles of pultruded box beams with high ultimate strength

Pultruded box beams under bending may be subjected to local buckling which causes premature failure of the beams. As such it is important to design pultruded box beams with high local buckling resistance to increase their ultimate strength. This paper presents an optimum design approach for cross-sectional profiles of pultruded box beams of (approximately) the same mass with emphasis on accomplishment of high local buckling resistance through finite element analysis. Five different sectional profiles have been designed by stiffening a simple box, and finite element analysis is used to study linear buckling and postbuckling of the beams. Results for the critical loads of linear buckling and local buckling judged by stress variation, stresses and deformations in postbuckling are presented. The computational results show one of the proposed sectional profiles does not develop local buckling and produces much smaller stresses and deformations within the load range of interest.     

Post buckling analysis of shear deformable shallow shells by the boundary element method

This paper presents four boundary element formulations for post buckling analysis of shear deformable shallow shells. The main differences between the formulations rely on the way non‐linear terms are treated and on the number of degrees of freedom in the domain. Boundary integral equations are obtained by coupling boundary element formulation of shear deformable plate and two‐dimensional plane stress elasticity. Four different sets of non‐linear integral equations are presented. Some domain integrals are treated directly with domain discretization whereas others are dealt indirectly with the dual reciprocity method. Each set of non‐linear boundary integral equations are solved using an incremental approach, where loads and prescribed boundary conditions are applied in small but finite increments. The resulting systems of equations are solved using a purely incremental technique and the Newton–Raphson technique with the Arc length method. Finally, the effect of imperfections (obtained from a linear buckling analysis) on the post‐buckling behaviour of axially compressed shallow shells is investigated. Results of several benchmark examples are compared with the published work and good agreement is obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

Second‐order analysis and design of steel structures allowing for member and frame imperfections

Conventional linear analysis is deficient in handling the design of slender frames since the effective length and other non‐linear effects are difficult to assess accurately. Some proposed non‐linear analyses cannot be directly employed for practical design since they are unable to re‐produce the buckling curves of the basic structural element, a simple column under axial force, by a single element per member. This paper describes an advanced element, using the same physical significance as the advanced and second‐order analysis proposed by a number of researchers and the BS5950(2000) and AS4100 (1995), for practical design of slender steel frames. The proposed element captures the physical behaviour of a structural member that the buckling strength of the member can be predicted using a single element per member and without assuming any effective length. Copyright © 2004 John Wiley & Sons, Ltd.     

Non‐linear analysis of structures using high performance hybrid elements

In this work, we describe the formulation and implementation for stress‐based hybrid elements for conducting non‐linear analysis of elastic structures. The motivation behind developing these elements is that they should be as simple to use as standard displacement‐based isoparametric brick elements, but at the same time, be relatively immune to the shortcoming that these elements suffer from, namely, ‘locking’ problems which occur when they are used to model plate/shell geometries, almost incompressible materials or when the elements are distorted, and so on. The formulation is based on a two‐field mixed variational principle. Numerical examples are presented to demonstrate the excellent performance of the proposed elements on a variety of challenging problems involving very large deformations, buckling, mesh distortions, almost incompressible materials, etc. Copyright © 2006 John Wiley & Sons, Ltd.     

A robust non‐linear mixed hybrid quadrilateral shell element

In the paper a non‐linear quadrilateral shell element for the analysis of thin structures is presented. The variational formulation is based on a Hu–Washizu functional with independent displacement, stress and strain fields. The interpolation matrices for the mid‐surface displacements and rotations as well as for the stress resultants and strains are specified. Restrictions on the interpolation functions concerning fulfillment of the patch test and stability are derived. The developed mixed hybrid shell element possesses the correct rank and fulfills the in‐plane and bending patch test. Using Newton's method the finite element approximation of the stationary condition is iteratively solved. Our formulation can accommodate arbitrary non‐linear material models for finite deformations. In the examples we present results for isotropic plasticity at finite rotations and small strains as well as bifurcation problems and post‐buckling response. The essential feature of the new element is the robustness in the equilibrium iterations. It allows very large load steps in comparison to other element formulations. Copyright © 2005 John Wiley & Sons, Ltd.     

Tracing the equilibrium path beyond compound critical points

Procedures to handle critical points with coincident or nearly coincident buckling loads in a geometrically non‐linear continuation process are presented. The proposed branch switching algorithm is based on a Liapunov–Schmidt–Koiter‐type asymptotic reduction which is used for the prediction step into the post‐buckling range. A critical review of the existing methodology for the numerical treatment of post‐bifurcation branches at a multiple bifurcation point is also given. In the review the main emphasis is given to robustness, but also implementational issues and computational requirements are addressed. Some modifications to the existing algorithms are proposed in order to improve their robustness or efficiency. Numerical examples of a flat plate and a thin‐walled structure are shown. Copyright © 1999 John Wiley & Sons, Ltd.     

A refined non‐linear non‐conforming triangular plate/shell element

A refined non‐conforming triangular plate/shell element for geometric non‐linear analysis of plates/shells using the total Lagrangian/updated Lagrangian approach is constructed in this paper based on the refined non‐conforming element method for geometric non‐linear analysis. The Allman's triangular plane element with vertex degrees of freedom and the refined triangular plate‐bending element RT9 are used to construct the present element. Numerical examples demonstrate that the accuracy of the new element is quite high in the geometric non‐linear analysis of plates/shells. Copyright © 2003 John Wiley & Sons, Ltd.     

Buckling of thin‐walled cylindrical shells under axial compression

Lightweight thin‐walled cylindrical shells subjected to external loads are prone to buckling rather than strength failure. The buckling of an axially compressed shell is studied using analytical, numerical and semi‐empirical models. An analytical model is developed using the classical shell small deflection theory. A semi‐empirical model is obtained by employing experimental correction factors based on the available test data in the theoretical model. Numerical model is built using ANSYS finite element analysis code for the same shell. The comparison reveals that the analytical and numerical linear model results match closely with each other but are higher than the empirical values. To investigate this discrepancy, non‐linear buckling analyses with large deflection effect and geometric imperfections are carried out. These analyses show that the effects of non‐linearity and geometric imperfections are responsible for the mismatch between theoretical and experimental results. The effect of shell thickness, radius and length variation on buckling load and buckling mode has also been studied. Copyright © 2009 John Wiley & Sons, Ltd.     

Improvement of semi‐analytical design sensitivities of non‐linear structures using equilibrium relations

The accuracy problem of the semi‐analytical method for shape design sensitivity analysis has been reported for linear and non‐linear structures. The source of error is the numerical differentiation of the element internal force vector, which is inherent to the semi‐analytical approach. Such errors occur for structures whose displacement field is characterized by large rigid body rotations of individual elements. This paper presents a method for the improvement of semi‐analytical sensitivities. The method is based on the element free body equilibrium conditions, and on the exact differentiation of the rigid body modes. The method is efficient, simple to code, and can be applied to linear and non‐linear structures. The numerical examples show that this approach eliminates the abnormal errors that occur in the conventional semi‐analytical method. Copyright © 2001 John Wiley & Sons, Ltd.     

Effects of approximations in analyses of beams of open thin‐walled cross‐section—part II: 3‐D non‐linear behaviour

In a companion paper, the effects of approximations in the flexural‐torsional stability analysis of beams was studied, and it was shown that a second‐order rotation matrix was sufficiently accurate for a flexural‐torsional stability analysis. However, the second‐order rotation matrix is not necessarily accurate in formulating finite element model for a 3‐D non‐linear analysis of thin‐walled beams of open cross‐section. The approximations in the second‐order rotation matrix may introduce ‘self‐straining’ due to superimposed rigid‐body motions, which may lead to physically incorrect predictions of the 3‐D non‐linear behaviour of beams. In a 3‐D non‐linear elastic–plastic analysis, numerical integration over the cross‐section is usually used to check the yield criterion and to calculate the stress increments, the stress resultants, the elastic–plastic stress–strain matrix and the tangent modulus matrix. A scheme of the arrangement of sampling points over the cross‐section that is not consistent with the strain distributions may lead to incorrect predictions of the 3‐D non‐linear elastic–plastic behaviour of beams. This paper investigates the effects of approximations on the 3‐D non‐linear analysis of beams. It is found that a finite element model for 3‐D non‐linear analysis based on the second‐order rotation matrix leads to over‐stiff predictions of the flexural‐torsional buckling and postbuckling response and to an overestimate of the maximum load‐carrying capacities of beams in some cases. To perform a correct 3‐D non‐linear analysis of beams, an accurate model of the rotations must be used. A scheme of the arrangement of sampling points over the cross‐section that is consistent with both the longitudinal normal and shear strain distributions is needed to predict the correct 3‐D non‐linear elastic–plastic behaviour of beams. Copyright © 2001 John Wiley & Sons, Ltd.     

A spatially curved‐beam element with warping and Wagner effects

This paper presents a new spatially curved‐beam element with warping and Wagner effects that can be used for the non‐linear large displacement analysis of members that are curved in space. The non‐linear behaviour of members curved in space shows that the Wagner effects are substantial in the large twist rotation analysis. Most existing finite beam element models, such as ABAQUS and ANSYS cannot predict the non‐linear large displacement response of members curved in space correctly because the Wagner effects, viz. the Wagner moment and the corresponding finite strain terms, have not been considered in these finite beam elements. As a consequence, these finite beam elements do not provide correct predictions for the out‐of‐plane buckling and postbuckling behaviour of arches as well. In this paper, the symmetric tangent stiffness matrix has been derived based on the finite rotations parameterized by the conventional displacements. The warping and Wagner effects: both the Wagner moment and the corresponding finite strain terms and their constitutive relationship, are included in the spatially curved‐beam element. Two components of the initial curvature, the initial twist and their interactions with the displacements are also considered in the spatially curved‐beam element. This ensures that the large twist rotation analysis for the members curved in space is accurate. Comparisons with existing experimental, analytical and numerical results show that the spatially curved‐beam element is accurate and efficient for the non‐linear elastic analysis of curved members, buckling and postbuckling analysis of arches, and in its ability to predict large deflections and twist rotations in more arbitrarily curved members. Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of SMA composite laminates using a multiscale modelling technique

The present work deals with the analysis of smart laminates, obtained as stacking sequence of fibre‐reinforced composite laminae and composite shape memory alloy (SMA) layers. The behaviour of composite SMA (CSMA) laminate is studied developing a full micro–macro approach. In fact, a non‐linear 4‐node mixed interpolation of tensorial components (MITC4) laminate finite element, based on the first‐order shear deformation theory, is developed. The SMA layer constitutive relationship is determined solving a non‐linear homogenization problem at each non‐linear iteration of each time step for each integration Gauss point. Some numerical applications are developed in order to investigate the influence of the CSMA on the buckling behaviour of plates and on the transversal displacement control of plates subjected to different loading conditions. Copyright © 2006 John Wiley & Sons, Ltd.     

Effects of approximations in analyses of beams of open thin‐walled cross‐section—part I: Flexural–torsional stability

In formulating a finite element model for the flexural–torsional stability and 3‐D non‐linear analyses of thin‐walled beams, a rotation matrix is usually used to obtain the non‐linear strain–displacement relationships. Because of the coupling between displacements, twist rotations and their derivatives, the components of the rotation matrix are both lengthy and complicated. To facilitate the formulation, approximations have been used to simplify the rotation matrix. A simplified small rotation matrix is often used in the formulation of finite element models for the flexural–torsional stability analysis of thin‐walled beams of open cross‐section. However, the approximations in the small rotation matrix may lead to the loss of some significant terms in the stability stiffness matrix. Without these terms, a finite element line model may predict the incorrect flexural–torsional buckling load of a beam. This paper investigates the effects of approximations in the elastic flexural–torsional stability analysis of thin‐walled beams, while a companion paper investigates the effects of approximations in the 3‐D non‐linear analysis. It is found that a finite element line model based on a small rotation matrix may predict incorrect elastic flexural–torsional buckling loads of beams. To perform a correct flexural–torsional stability analysis of thin‐walled beams, modification of the model is needed, or a finite element model based on a second‐order rotation matrix can be used. Copyright © 2001 John Wiley & Sons, Ltd.     

A Koiter‐Newton approach for nonlinear structural analysis

A new approach termed the Koiter‐Newton is presented for the numerical solution of a class of elastic nonlinear structural response problems. It is a combination of a reduction method inspired by Koiter's post‐buckling analysis and Newton arc‐length method so that it is accurate over the entire equilibrium path and also computationally efficient in the presence of buckling. Finite element implementation based on element independent co‐rotational formulation is used. Various numerical examples of buckling sensitive structures are presented to evaluate the performance of the method. The examples demonstrate that the method is robust and completely automatic and that it outperforms traditional path‐following techniques. This improved efficiency will open the door for the direct use of detailed nonlinear finite element models in the design optimization of next generation flight and launch vehicles. Copyright © 2013 John Wiley & Sons, Ltd.     

A non‐stationary covariance‐based Kriging method for metamodelling in engineering design

Metamodels are widely used to facilitate the analysis and optimization of engineering systems that involve computationally expensive simulations. Kriging is a metamodelling technique that is well known for its ability to build surrogate models of responses with non‐linear behaviour. However, the assumption of a stationary covariance structure underlying Kriging does not hold in situations where the level of smoothness of a response varies significantly. Although non‐stationary Gaussian process models have been studied for years in statistics and geostatistics communities, this has largely been for physical experimental data in relatively low dimensions. In this paper, the non‐stationary covariance structure is incorporated into Kriging modelling for computer simulations. To represent the non‐stationary covariance structure, we adopt a non‐linear mapping approach based on parameterized density functions. To avoid over‐parameterizing for the high dimension problems typical of engineering design, we propose a modified version of the non‐linear map approach, with a sparser, yet flexible, parameterization. The effectiveness of the proposed method is demonstrated through both mathematical and engineering examples. The robustness of the method is verified by testing multiple functions under various sampling settings. We also demonstrate that our method is effective in quantifying prediction uncertainty associated with the use of metamodels. Copyright © 2006 John Wiley & Sons, Ltd.     

A geometrically non‐linear piezoelectric solid shell element based on a mixed multi‐field variational formulation

This paper is concerned with a geometrically non‐linear solid shell element to analyse piezoelectric structures. The finite element formulation is based on a variational principle of the Hu–Washizu type and includes six independent fields: displacements, electric potential, strains, electric field, mechanical stresses and dielectric displacements. The element has eight nodes with four nodal degrees of freedoms, three displacements and the electric potential. A bilinear distribution through the thickness of the independent electric field is assumed to fulfill the electric charge conservation law in bending dominated situations exactly. The presented finite shell element is able to model arbitrary curved shell structures and incorporates a 3D‐material law. A geometrically non‐linear theory allows large deformations and includes stability problems. Linear and non‐linear numerical examples demonstrate the ability of the proposed model to analyse piezoelectric devices. Copyright © 2005 John Wiley & Sons, Ltd.