Embracing nonlinearities in structural design

1
www.bristol.ac.uk/composites Figure 3: Λ 2 continuation; at the most efficient structure of Λ 3 = 0.99993, the total volume of the structure varies from 200 2 = 1) to 1000 2 = 5); (a) (b) Figure 2: Λ 3 continuation (a) Equilibrium surface of beam problem with fixed average area Λ 2 =1 in multi- parameter space. The black line represents the most efficient structure; (b) Top-down view of Figure (a); Blue, stable; Red, one eigenvalue unstable; Green, two eigenvalues unstable Embracing nonlinearities in structural design Bradley Cox 1 , Rainer Groh 1 , Daniele Avitabile 2 and Alberto Pirrera 1 Introduction Nonlinear phenomena are prevalent in structural and continuum mechanics, albeit not commonly utilised. This general reluctance is due to the lack of availability of sufficiently robust computational tools to solve these complex problems. There appears to be little question that the so-called incremental iterative methods represent by far the most popular procedures for the solution of nonlinear structural mechanics. Of these, numerical continuation provides many advantages. The noteworthy, yet often overlooked, quality is its application in optimising engineering structures. Numerical continuation is coupled to a nonlinear finite element analysis framework, where, for a benchmark nonlinear problem, arched beams are analysed using structural beam elements based on Timoshenko kinematics. Future research Numerical continuation Results and discussion Motivation Model definition Within the realm of nonlinear structural analysis it is generally understood that most structures are somewhat overdesigned. This is simply to ensure that they respond linearly to the externally applied loads. The motivation behind the current research is to enable engineers to design new structures that embrace well-behaved nonlinear deformations thus potentially leading to significantly more efficient structures, by removing redundant and unnecessary stiffness, and in turn mass. Numerical continuation is a numerical method used for computing solutions of a system of parameterised nonlinear equations. The method is independent from the source of nonlinearity - in terms of structural analysis - it can be geometric, inertial, material, or it can arise from other physical phenomena such as contact, damage, delamination, buckling driven delamination, etc. The numerical method offers a number of advantages. Most significantly, it allows for a systematic exploration of the design space, and it makes typically cumbersome and expensive nonlinear calculations, often requiring onerous user intervention, relatively straightforward. In essence, numerical continuation is built around the idea of the addition of some augmenting conditions to the equilibrium equations. These conditions, specifying subsets of equilibrium states with certain properties, will further allow the addition of the extra control parameters required [1] in any analysis. Problem formulation 1 ACCIS, Queen’s Building, University of Bristol, Bristol, BS8 1TR, UK 2 School of Mathematical Sciences, University of Nottingham, NG7 2RD, UK The finite element equilibrium equations are formally expressed as where , the internal forces, are a function of both the current displacements, , and possibly some other properties of the model, . This additional property parameter could in theory be anything from a change in geometry (area, length, height, etc.), to a change in material properties (E, G, etc.) or auxiliary load cases. Model #1: Arched beam, Λ 1 Model #2 & #3: Arched beam, Λ 2 & Λ 3 In order to illustrate the capabilities of numerical continuation a simple arched beam problem has been evaluated. Here, the vector = {Λ 1 2 3 4 } comprises of only four parameters. The parameters are defined as follows: Λ 1 Externally applied load to centre-point of arch. Λ 2 Average cross-sectional area (used to evaluate total volume). Λ 3 Quadratic distribution of volume along beam length (symmetric). Λ 4 Arch height. This model delivers classic snap-through curves as illustrated in Figure 1a. The models presented herein are based on one fundamental load-displacement curve, where the load-factor 1 is the first parameter evaluated all other parameters remain fixed. This is necessary due to the process involved in continuing the parameters along an equilibrium path. The solutions at each time-step on this fundamental load-displacement path act as an initial guess i.e. starting positions for the evaluation of other parameters. Parameters Λ 2 and Λ 3 are once again solved independently of one another, but both are related to the change in cross-sectional area along the beam length, as illustrated in the following expression: Λ 3 is physically limited to 0.99995 ≤ Λ 3 ≤ 1.0002 to avoid elements containing zero or negative volume (see Figure 1). Λ 3 =1 corresponds to a constant cross-sectional area along the beam length. Λ 2 relates directly to the average area of the beam: when divided by the overall length it is used to evaluate the overall volume of the beam and hence mass. References To expand on our understanding of more complex nonlinear structures and to undertake comprehensive imperfection sensitivity analyses. Figure 1: Beam problem schematics; (a) Free-body diagram; (b) Parameter Λ 3 with a maximum beam thickness at the centre; (c) Parameter Λ 3 with the maximum thickness at either end. [1] C. Pacoste, A. Eriksson, A. Zdunek, Parameter dependence in the critical behaviour of shell structures: a numerical approach, IASS-IACM 2000 Figure 2a illustrates an equilibrium surface in which each load-displacement curve is solved varying the distribution of volume along the length of the beam. Each curve corresponds to a fixed load. The most efficient design (Figure 2a: black curve) is found by evaluating each point on the surface to find the best ratio of load-to- displacement. This design is found at Λ 3 = 0.99993; which incidentally is the solution at which the first bifurcation point and the first limit point coincide; this is illustrated graphically in Figure 2b as the point where all three colours coincide. The design of the arch is further evaluated with the continuation along parameter Λ 2 which essentially alters the solution into a 4 th dimension (see Figure 3). In this instance the equilibrium surface, and hence design space, is expanded. Figure 3 can be used to design minimum weight structures for any given load or required displacement. The most important feature is the blue region. These solutions correspond to stable, and hence real, solutions that are physically observed.

Transcript of Embracing nonlinearities in structural design

Page 1: Embracing nonlinearities in structural design

www.bristol.ac.uk/composites

Figure 3: Λ2 continuation; at the most efficient structure of Λ3 = 0.99993, the total volume of the

structure varies from 200 (Λ2 = 1) to 1000 (Λ2 = 5);

(a)

(b)

Figure 2: Λ3 continuation (a) Equilibrium surface of beam problem with fixed average area Λ2 = 1 in multi-parameter space. The black line represents the most efficient structure; (b) Top-down view of Figure (a);

Blue, stable; Red, one eigenvalue unstable; Green, two eigenvalues unstable

Embracing nonlinearities in structural designBradley Cox1, Rainer Groh1, Daniele Avitabile2 and Alberto Pirrera1

Introduction

Nonlinear phenomena are prevalent instructural and continuum mechanics, albeit notcommonly utilised. This general reluctance isdue to the lack of availability of sufficientlyrobust computational tools to solve thesecomplex problems. There appears to be littlequestion that the so-called incremental iterativemethods represent by far the most popularprocedures for the solution of nonlinearstructural mechanics. Of these, numericalcontinuation provides many advantages. Thenoteworthy, yet often overlooked, quality is itsapplication in optimising engineeringstructures.

Numerical continuation is coupled to anonlinear finite element analysis framework,where, for a benchmark nonlinear problem,arched beams are analysed using structuralbeam elements based on Timoshenkokinematics.

Future research

Numerical continuation

Results and discussion

Motivation

Model definition

Within the realm of nonlinear structuralanalysis it is generally understood that moststructures are somewhat overdesigned. This issimply to ensure that they respond linearly tothe externally applied loads. The motivationbehind the current research is to enableengineers to design new structures thatembrace well-behaved nonlinear deformationsthus potentially leading to significantly moreefficient structures, by removing redundant andunnecessary stiffness, and in turn mass.

Numerical continuation is a numerical methodused for computing solutions of a system ofparameterised nonlinear equations. Themethod is independent from the source ofnonlinearity - in terms of structural analysis - itcan be geometric, inertial, material, or it canarise from other physical phenomena such ascontact, damage, delamination, buckling drivendelamination, etc.

The numerical method offers a number ofadvantages. Most significantly, it allows for asystematic exploration of the design space, andit makes typically cumbersome and expensivenonlinear calculations, often requiring oneroususer intervention, relatively straightforward.

In essence, numerical continuation is builtaround the idea of the addition of someaugmenting conditions to the equilibriumequations. These conditions, specifying subsetsof equilibrium states with certain properties,will further allow the addition of the extracontrol parameters required [1] in any analysis.

Problem formulation

1 ACCIS, Queen’s Building, University of Bristol,

Bristol, BS8 1TR, UK2

School of Mathematical Sciences, University of

Nottingham, NG7 2RD, UK

The finite element equilibrium equations areformally expressed as

where 𝒇, the internal forces, are a function ofboth the current displacements, 𝒅 , andpossibly some other properties of the model,𝚲. This additional property parameter could intheory be anything from a change in geometry(area, length, height, etc.), to a change inmaterial properties (E, G, etc.) or auxiliaryload cases.

Model #1: Arched beam, Λ1

Model #2 & #3: Arched beam, Λ2 & Λ3

In order to illustrate the capabilities of numericalcontinuation a simple arched beam problem hasbeen evaluated.

Here, the vector 𝚲 = {Λ1, Λ2, Λ3, Λ4} comprises ofonly four parameters. The parameters are definedas follows:• Λ1 Externally applied load to centre-point of

arch.• Λ2 Average cross-sectional area (used to

evaluate total volume).• Λ3 Quadratic distribution of volume along

beam length (symmetric).• Λ4 Arch height.

This model delivers classic snap-through curvesas illustrated in Figure 1a.

The models presented herein are based on onefundamental load-displacement curve, where theload-factor 𝜆 = Λ1 is the first parameter evaluated– all other parameters remain fixed. This isnecessary due to the process involved incontinuing the parameters along an equilibriumpath. The solutions at each time-step on thisfundamental load-displacement path act as aninitial guess i.e. starting positions for theevaluation of other parameters.

Parameters Λ2 and Λ3 are once again solvedindependently of one another, but both arerelated to the change in cross-sectional areaalong the beam length, as illustrated in thefollowing expression:

Λ3 is physically limited to 0.99995 ≤ Λ3 ≤ 1.0002 toavoid elements containing zero or negativevolume (see Figure 1). Λ3 = 1 corresponds to aconstant cross-sectional area along the beamlength. Λ2 relates directly to the average area ofthe beam: when divided by the overall length it isused to evaluate the overall volume of the beamand hence mass.

References

To expand on our understanding of more complexnonlinear structures and to undertakecomprehensive imperfection sensitivity analyses.

Figure 1: Beam problem schematics; (a) Free-body diagram; (b) Parameter Λ3 with a maximum beam thickness at the centre; (c) Parameter Λ3 with the

maximum thickness at either end.

[1] C. Pacoste, A. Eriksson, A. Zdunek, Parameter dependence in the critical behaviour of shell structures: a numerical approach, IASS-IACM 2000

Figure 2a illustrates an equilibrium surface inwhich each load-displacement curve is solvedvarying the distribution of volume along thelength of the beam. Each curve corresponds to afixed load.

The most efficient design (Figure 2a: blackcurve) is found by evaluating each point on thesurface to find the best ratio of load-to-displacement. This design is found at Λ3 = 0.99993;which incidentally is the solution at which the first

bifurcation point and the first limit pointcoincide; this is illustrated graphically in Figure2b as the point where all three colours coincide.

The design of the arch is further evaluatedwith the continuation along parameter Λ2 whichessentially alters the solution into a 4th

dimension (see Figure 3). In this instance theequilibrium surface, and hence design space, isexpanded. Figure 3 can be used to designminimum weight structures for any given load orrequired displacement. The most importantfeature is the blue region. These solutionscorrespond to stable, and hence real, solutionsthat are physically observed.