# Robust Optimal Design of Beams Subject to Uncertain Loads

Authors: Ismail Kucuk; Sarp Adali; Ibrahim Sadek
Journal: Mathematical Problems in Engineering
 (2009)

---

## Abstract

Optimality conditions are derived for the robust optimal design of beams subject to a combination of uncertain and deterministic transverse and boundary loads using a variational min-max approach. The potential energy of the beam is maximized to compute the worst case loading and minimized to determine the optimal cross-sectional shape which results in coupled nonlinear differential equations for the unknown functions except for the case of a variable width beam. The uncertain component of the transverse load acting on the beam is not known a priori resulting in load uncertainty subject only to an norm constraint. Similarly the optimal area function is subject to a volume constraint leading to an isoperimetric variational problem. The min-max approach leads to robust optimal designs which are not susceptible to unexpected load variations as it occurs under operational conditions. The solution methodology is illustrated for the variable width beam by obtaining analytical results for several cases. The efficiency of the optimal designs is computed with respect to a uniform beam under worst case loading taking the maximum deflection as the quantity for comparison. It is observed that the optimal shapes are more than 70% efficient for the examples given in this study.


---

## Body

## 1. Introduction

Under operational conditions, a structure is usually subjected to uncertainties which may arise from fluctuation and scatter of external loads, environmental conditions, boundary conditions, and geometrical and material properties. However quite often, design uncertainties arise from incomplete knowledge and unpredictable nature of the load under operational conditions.In the present study optimality conditions are derived for the optimal robust design of a beam subject to a transverse load which has unknown and known parts. Moreover uncertain moments and/or shear forces may be acting at the endpoints of the beam. In conventional design, it is common practice to neglect the load uncertainties when analyzing a structure and assessing the structural performance on the basis of a deterministic model. To compensate for performance variability caused by load variations, a safety factor is introduced magnitude of which correlates with the level of uncertainty with higher levels leading to larger safety factors. However, the safety factors specified may be either too conservative or too small to compensate for the lack of knowledge of operational loads. Efficiency and reliability of the structure can be improved by taking the load uncertainties into consideration in the design process leading to a design which is robust under load variations. This approach is equivalent to optimizing the design for worst case loading and leads to an optimization problem for the area function and to an antioptimization problem for the load function. The final design is robust in the sense that the sensitivity of the beam to load variations is substantially reduced. This is accomplished by maximizing its potential energy over loading while minimizing it with respect to its cross-sectional shape. Mathematically these results in a min-max optimal design problem can be studied by variational calculus.Finding the worst case loading on a structure corresponds to an antioptimization problem the examples of which can be found [1, 2] where optimization under uncertain bending and buckling loads is studied. Designing robust structures to carry loads that are not known in advance is discussed in papers [3–5] where the authors proposed a minmax formulation to maximize the compliance with respect to loading and to minimize it with respect to design variables. These problems lead to a coupled optimization/antioptimization formulation where the objective is to compute the “best” design under “worst” case loading resulting in a robust design insensitive to load variations.An alternative strategy to treat the uncertainties is convex modeling in which the uncertainties belong to a convex set [6–9]. This approach allows the designer to use not the averaged results but extremal properties of the system being modeled, according to the convex set chosen. The limitation of the convex modeling is that only small variations around a nominal value of the uncertain quantity can be considered and the model becomes less accurate as the variations become larger. Other methods of taking load uncertainties in the design process can be found in [10, 11]. The main objective of these techniques is to achieve robust designs which are not susceptible to failure under unexpected loads [12–14]. Recently the optimization under dimensional uncertainty was applied to the design of microbeams [15].In the present work, a system of nonlinear differential equations is derived in terms of state, design and load variables using calculus of variations and the methods of Lagrange multipliers and slack variables. These expressions serve as the optimality and antioptimality conditions of the problem. The design variable is the cross-sectional shape of the beam and is subject to a volume constraint. The load variables comprise a combination of deterministic and uncertain transverse loads as well as uncertain moments and shear forces which may act at the boundaries. The only constraints on the unknown loads involve finite norms and an upper bound on the transverse load. The optimization method involves a minmax formulation where the objective is to minimize the compliance with respect to the cross-sectional shape and maximize it with respect to the unknown loads. The formulation ensures that the optimal designs correspond to the most unfavorable loading and, therefore, these designs are conservative for any other loading. The solution methodology is illustrated with several examples involving cases which allow the computation of closed-form solutions.


## 2. Problem Formulation and Design Constraints

The differential equation governing the deflection of a variable cross-section beam subject to a combination of uncertain and deterministic loads, as shown in Figure1, can be expressed in nondimensional form asFigure 1
Beam diagram with external forces.(2.1)(an(x)′′(x))′′=F(x)+G(x),forx∈(0,1)
where I(x)=a(x)n is the stiffness function, a(x)∈C2[0,1] is the cross sectional area, y(x)∈C4[0,1] is the deflection, F(x)∈C0[0,1] is the unknown uncertain loading and G(x)∈C0[0,1] is the given deterministic loading. The primes denote the derivative with respect to x∈[0,1]. In (2.1), n=1,2,3 is a constant where n corresponds to a variable width, a geometrically similar cross-section and a variable height, respectively. The beam is subject to two boundary conditions at the endpoints x=0 and x=1. The boundary conditions may also contain uncertain and deterministic parts in the form of moments and shear forces which can be expressed as(2.2)ℬ0i(y)|x=0=f0+g0,ℬ1i(y)|x=1=f1+g1,i=1,2,
where ℬ0i and ℬ1i are boundary operators at x=0 and x=1, respectively, f0 and f1 are uncertain constants, and g0 and g1 are deterministic constants. In addition to the boundary conditions, the deflection, area and load functions should satisfy certain continuity conditions which arise as a result of discontinuities in the first derivative of the uncertain load function. Let the uncertain load function F(x) act on the beam between the points s1 and s2, that is, it is applied on the interval 0<s1≤x≤s2<1. In general F(x) can be expressed in the form(2.3)F(x)={0,if0<x≤s1,F1(x)≤Fmax,ifs1≤x≤d1,Fmax,ifd1≤x≤d2,F2(x)≤Fmax,ifd2≤x≤s2,0,ifs2≤x<1,
where d1 and d2 are unknown constants to be determined from the continuity conditions(2.4)F1(d1)=Fmax,F2(d2)=Fmax.
The optimal area function is computed for each interval and should be a continuous function of x which requires that(2.5)a-(s1)=a+(s1),a-(s2)=a+(s2),a-(d1)=a+(d1),a-(d2)=a+(d2),
where a- and a+ denote the values of the area function a(x) to the left and right of the relevant point. Similarly the moment and shear force should be continuous along the beam in the absence of point loads or moments, and in particular at these points that can be expressed as
(2.6)any′′|x=ξ-=any′′|x=ξ+,(any′′)′|x=ξ-=(any′′)′|x=ξ+
where ξ- and ξ+ are the points to the left and right of ξ=s1,s2,d1, and d2. In the sections on the implementation of the solution method, it is shown that the number of unknowns resulting from the integration of differential equations matches the number of boundary and continuity conditions leading to unique solutions.Physically the stiffness of a beam has a finite value and in practice the volume of the beam is specified leading to a design constraint which can be expressed mathematically as(2.7)∫01a(x)dx=1,
which simply constraints the volume of the beam available for optimization. Similarly the uncertain loads acting on the beam are required to have finite norms and an upper bound as would be the case under operating conditions even though these loadings are not known precisely. For the uncertain load F(x), these constraints can be expressed as(2.8)∥F(x)∥Lpp=∫01F(x)pdx=1,1<p<∞,max0≤0≤1F(x)≤Fmax,
that is, the Lp norm of the uncertain transverse load is constrained to have the value 1 and the uncertain transverse load is bounded from above. Finally, the uncertain moments and shear forces acting at the endpoints are of finite magnitude in practice that can be expressed as(2.9)m0p+m1p=η,v0p+v1p=γ,
where η>0 and γ>0 are given constants.


### 2.1. Objective Functional

The design optimization can be achieved by choosing a suitable performance index for the problem which serves as an objective functional of a minmax problem. In the present case a suitable objective functional is the potential energy of the beam which is given by(2.10)PE(a,F,m,v;y)=12∫01an(y′′)2dx-∫01(F+G)ydx+12(y′(0)m(0)-y′(1)m(1)+y(0)v(0)-y(1)v(1)),
the Euler-Lagrange equation of which yields the state equation (2.1) and the natural boundary conditions. Thus the functional PE(a,F,m,v;y) has the advantage of producing the deflection function y(x) when minimized with respect to y. Another objective functional can be defined as the compliance of the beam given by
(2.11)𝒞(a,F,m,v;y)=∫01(F+G)ydx,
which is a measure of the relative stiffness of the beam under distributed loads. Substituting (2.1) into (2.11) and performing integration by parts, we obtain(2.12)∫01(F+G)ydx=∫01an(y′′)2dx=2PE(a,F,m,v;y)+2𝒞(a,F,m,v;y).
Equations (2.11) and (2.12) show that the potential energy and the beam compliance are related as
(2.13)PE(a,F,m,v;y)=-12𝒞(a,F,m,v;y).
As such the two objective functionals are closely related. In the present study the potential energy given by (2.10) is chosen as the performance index of the optimization problem.


### 2.2. Robust Optimal Design Problem

Find the cross-sectional areaa(x) of the optimal beam for given values of n,p and the deterministic loading G(x) such that the beam's potential energy PE(a,F,m,v;y) given by (2.10) is minimized subject to the volume constraint given by (2.7) under the worst case loadings with respect to the uncertain load F(x), uncertain end moments m=(m0,m1), and shear forces v=(v0,v1) where F(x),m, and s satisfy the constraints (2.8) and (2.9).The robust design problem constitutes an optimization problem with respect to the area functiona(x), and an antioptimization problem with respect to the uncertain loads F(x),m and s. Thus the solution of the following minmax problem is sought:(2.14)mina(x)maxF,m,vPE(a,F,m,v;y),
subject to constraints (2.8) and (2.9) where PE(a,F,m,v;y) is given by (2.10). The solution to the antioptimization problem yields the worst case loadings and the solution to the optimization problem yields the optimal cross-sectional area as a function of x. Solving this problem as a nested minmax problem leads to the optimal robust area function under worst case of loading and mathematically to a system of three nonlinear differential equations in three unknowns: y(x),a(x), and F(x).


## 2.1. Objective Functional

The design optimization can be achieved by choosing a suitable performance index for the problem which serves as an objective functional of a minmax problem. In the present case a suitable objective functional is the potential energy of the beam which is given by(2.10)PE(a,F,m,v;y)=12∫01an(y′′)2dx-∫01(F+G)ydx+12(y′(0)m(0)-y′(1)m(1)+y(0)v(0)-y(1)v(1)),
the Euler-Lagrange equation of which yields the state equation (2.1) and the natural boundary conditions. Thus the functional PE(a,F,m,v;y) has the advantage of producing the deflection function y(x) when minimized with respect to y. Another objective functional can be defined as the compliance of the beam given by
(2.11)𝒞(a,F,m,v;y)=∫01(F+G)ydx,
which is a measure of the relative stiffness of the beam under distributed loads. Substituting (2.1) into (2.11) and performing integration by parts, we obtain(2.12)∫01(F+G)ydx=∫01an(y′′)2dx=2PE(a,F,m,v;y)+2𝒞(a,F,m,v;y).
Equations (2.11) and (2.12) show that the potential energy and the beam compliance are related as
(2.13)PE(a,F,m,v;y)=-12𝒞(a,F,m,v;y).
As such the two objective functionals are closely related. In the present study the potential energy given by (2.10) is chosen as the performance index of the optimization problem.


## 2.2. Robust Optimal Design Problem

Find the cross-sectional areaa(x) of the optimal beam for given values of n,p and the deterministic loading G(x) such that the beam's potential energy PE(a,F,m,v;y) given by (2.10) is minimized subject to the volume constraint given by (2.7) under the worst case loadings with respect to the uncertain load F(x), uncertain end moments m=(m0,m1), and shear forces v=(v0,v1) where F(x),m, and s satisfy the constraints (2.8) and (2.9).The robust design problem constitutes an optimization problem with respect to the area functiona(x), and an antioptimization problem with respect to the uncertain loads F(x),m and s. Thus the solution of the following minmax problem is sought:(2.14)mina(x)maxF,m,vPE(a,F,m,v;y),
subject to constraints (2.8) and (2.9) where PE(a,F,m,v;y) is given by (2.10). The solution to the antioptimization problem yields the worst case loadings and the solution to the optimization problem yields the optimal cross-sectional area as a function of x. Solving this problem as a nested minmax problem leads to the optimal robust area function under worst case of loading and mathematically to a system of three nonlinear differential equations in three unknowns: y(x),a(x), and F(x).


## 3. Optimality Conditions

The derivations of the optimality and antioptimality conditions are given next. In view of the presence of several constraints on the design problem, the Lagrange multiplier technique is implemented by introducing the Lagrangian given by(3.1)L(a,F,m,v;y)=PE(a,F,m,v;y)+μ1(∫01F(x)pdx-1)+μ2(∫01a(x)dx-1)+μ3(m0p+m1p-η)+μ4(v0p+v1p-γ)+∫01μ5(x)(F(x)-h(x)2-Fmax)dx,
where μi,i=1,2,3,4, and μ5(x) are Lagrange multipliers and h(x) is a slack variable. The variation of L(a,F,m,v;y) with respect to y gives the differential equation and the boundary conditions. The variation of L(a,F,m,v;y) with respect to a(x) yields(3.2)∫01(nan-1(y′′)2-μ2)δadx=0,
where δa is arbitrary. Thus, from the fundamental theorem of calculus of variations, it follows that(3.3)nan-1(y′′)2-μ2=0.The variation ofL(a,F,m,v;y) with respect to F(x) yields(3.4)∫01(-y+μ1pFp-1+μ5(x))δF(x)dx=0,
where δF is arbitrary. Thus,(3.5)-y+μ1pFp-1+μ5(x)=0.The variation ofL(a,F,m,v;y) with respect to h(x) yields(3.6)μ5(x)h(x)=0.
This equation implies that h(x) must be equal to zero when μ5(x)≠0, and μ5(x) must be zero when h(x)≠0. As a result, μ5(x) in equation (3.5) can be discarded since it is equal to zero when the second constraint in (2.8) is not active. Differentiation of L(a,F,m,v;y) with respect to m0, and m1 as well asv0, and v1 yields(3.7)y'(0)+2pμ3m0p-1=0,-y'(1)+2pμ3m1p-1=0,y(0)+2pμ4v0p-1=0,-y(1)+2pμ4v1p-1=0.
Using equations (2.1) and (3.5) results in the following expressions for y′′ and a(x):(3.8)y′′=c11a(1-n)/2,(3.9)a(x)=c11-1∫0x(x-t)(F(t)+G(t))dt+a1x+a2,forn=1(3.10)=c12(y′′)2/(1-n),forn=2,3,
where a1,a2, and c11 are arbitrary constants and c12=c112/(n-1). Equations (3.8) and (3.9) are optimality conditions for the area function a(x). Similarly, from (3.5) and (3.6), one obtains(3.11)y={c21Fp-1,forF<Fmax,c21Fmaxp-1,forF≥Fmax,F(x)={c22y1/(p-1),fory<(Fmaxc22)p-1,Fmax,fory≥(Fmaxc22)p-1,
where c21 is an arbitrary constant and c22=c211/(1-p). Equation (3.11) is the antioptimality condition for the worst case loading F(x). Substituting a(x) and F(x) from (3.9) and (3.11), respectively, into (2.1), we obtain a fourth-order nonlinear differential equation for y(x) given by(3.12)(c12n(y′′)(1+n)/(1-n))′′=c22y1/(p-1)+G(x)fory<(Fmaxc22)p-1,(3.13)(c12n(y′′)(1+n)/(1-n))′′=Fmax+G(x)fory≥(Fmaxc22)p-1,
where n=2,3. It is noted that for the beam with a variable width, that is, n=1, the optimality condition (3.8) uncouples the functions y(x) and a(x) that lead to linear differential equations for the unknown functions. This case is treated in more detail in the examples given to illustrate the solution methodology. Alternatively, (3.12) can be expressed as an integrodifferential equation by integrating it twice, namely,(3.14)(c12n(y′′)(1+n)/(1-n))′′=∫0x(x-t)(c22y1/(p-1)+G(t))dt+c3x+c4,
where c3 and c4 are integration constants which are to be determined from the boundary conditions (2.2). In the presence of uncertain moments and shear forces, the boundary conditions are obtained by noting that(3.15)m0=c5(-y′(0))1/(p-1),m1=c5(y′(1))1/(p-1),(3.16)v0=c6(-y(0))1/(p-1),v1=c6(y(0))1/(p-1),
where c5=(2pμ3)1/(1-p) and c6=(2pμ4)1/(1-p). If the boundary x=0 is subject to uncertain and deterministic moments m(0), then from (2.2) and (3.15), it follows that(3.17)a(0)ny′′(0)-c5(-y′(0))1/(p-1)=m0d.
Similarly, if uncertain and deterministic shear force v(1) is applied at x=1, then(3.18)a(1)ny′′(1)-c6(-y′(1))1/(p-1)=v1d,
for which (2.2) and (3.16) are used. Once the solution for y(x) is computed, the optimal area function aopt(x) and the worst case loading Fworst are obtained from (3.9)–(3.11). A nonlinear differential equation in terms of F(x) can also be obtained by noting that(3.19)a(1-n)/2=(c21c11)(Fp-1)′′,n=2,3,
which follows from equations (3.10) and (3.11). Substituting (3.8) and (3.19) into (2.1), we obtain
(3.20)c7(((Fp-1)′′)2n/(1-n)((Fp-1)′′)-1)′′=F(x)+G(x),
where c7=c112n/(n-1)c21(n+1)/(1-n) is an arbitrary constant.


## 4. Analytical Solutions

The theoretical framework developed in Sections2 and 3 to solve optimal design problems under load uncertainty, in general, requires the use of numerical methods for its implementation. However, the special case of n=1 (variable width beam) and p=2 (L2 norm of the transverse load) can be studied using closed-form solutions. These solutions, in turn, can be used to illustrate the method and investigate the efficiency of the robust designs. Moreover explicit solutions for the optimal beam shape and the worst case loadings serve as benchmark results for other cases where the solutions can be obtained numerically. In this section, the problem to be solved is summarized briefly and closed-form solutions are obtained for the optimal area function and the uncertain loads.We consider a simply supported beam subject to an uncertain load,F(x),0≤x≤1, which may be acting on part of the beam and may have an upper limit as shown in Figure 1 where G(x) is the deterministic component of the transverse load. The uncertain load can be defined as
(4.1)F(x)={0,if0≤x≤s1,f(x),ifs1≤x≤s2,0,ifs2≤x≤1,
where s1 and s2 are given parameters and f(x)∈C0[s1,s2] is an unknown continuous function. In addition the beam is subjected to uncertain moments m0 and m1 at the boundaries x=0 and x=1, respectively (Figure 1).Governing equation of the deflection of the beam given by (2.1) becomes(4.2)(a(x)y(x)′′)′′=F(x)+G(x),
for n=1. The constraints given by (2.7)–(2.9) become(4.3)∫01a(x)dx=1,∥F(x)∥L22≡∫01F(x)2dx=1,max0≤x≤1F(x)≤fmax,m02+m12=η.
For a simply supported beam, the boundary conditions are given by(4.4)y(0)=0,m(0)=a(0)y′′(0)=m0,y(1)=0,m(1)=a(1)y′′(1)=-m1,
where m(0) and m(1) are moments at the boundary points x=0 and x=1. The design problem involves the minimization of the potential energy of the beam under worst case of loading and as such involves optimization with respect to the area function a(x) and antioptimization with respect to the loading functions F(x), m0 and m1 subject to the constraints (4.3). This problem can be expressed as a minmax problem,namely,(4.5)mina(x)maxF(x),mPE(a(x),F(x),m;y),
where PE is the performance index (potential energy) given by(4.6)PE(a,F,m;y)=PE(a(x),F(x),m;y)=12∫01a(x)(y′′)2dx-∫01(F(x)+G(x))ydx+12m0y′(0)-12m1y′(1),
and m denotes the vector m=(m0,m1). In (4.5), the first term is the strain energy and the second, third and fourth terms make up the potential energy of the external loadings.For the computation of the optimal area functiona(x) and the worst case loading F(x),m0 and m1 subject to the constraints (4.3), the Lagrangian given by (3.1) becomes(4.7)L(a,F,m;y)=PE(a,F,m;y)+μ1(∫01F(x)2dx-1)+μ2(∫01a(x)dx-1)+μ3(m02+m12-η).
The variation of L(a,F,m;y) with respect to y gives the differential equation (4.2) and the boundary conditions (4.4). The variation of L(a,F,m;y) with respect to a(x) yields(4.8)∫01(y′′)2δadx+μ2∫01δadx=0,
where δa is arbitrary. Thus, from the fundamental theorem of calculus of variations, it follows that
(4.9)(y′′)2+μ2=0.
Similarly, the variations of L(a,F,m;y) with respect to F and m yield(4.10)-y+4μ1f(x)=0,forx∈[s1,s2],(4.11)y′(0)+4μ3m0=0,-y′(1)+2μ3m1=0.
Thus, the optimality condition of the problem is given by(4.12)y′′=constant=β.
Similarly, the antioptimization conditions can be expressed as(4.13a)f(x)={y2μ1,fory>2μ1fmax,fmax,fory<2μ1fmax,wheres1≤x≤s2,(4.13b)m0=-y′(0)4μ3,m0=y′(1)4μ3.
Substituting the optimality and antioptimality conditions into the differential equation (4.10), we obtain(4.14)a′′(x)={G(x)β,for0≤x≤s1,f(x)+G(x)β,fors1≤x≤s2,G(x)β,fors2≤x≤1,
where f(x) is given by (4.13a). A system of linear differential equations in y(x) and a(x) given by (4.12), (4.13a), and (4.14) can be solved simultaneously. In the present case it is possible to find an analytical solution for y(x) satisfying the boundary conditions, namely,
(4.15)y=β2x(x-1),for0≤x≤1.
Similarly, the optimal area function is given by(4.16)aopt(x)={G(x)β+c1x+c2,for0≤x≤s1,x348μ1(x-2)+1βG(x)+c3x+c4,fors1≤x≤s2,G(x)β+c5x+c6,fors2≤x≤1,
when y>2μ1fmax where G(x) is the second indefinite integral of G(x) and ci,i=1,…,6 are integration constants to be determined from the boundary conditions (4.4), and continuity conditions(4.17)a-(s1)=a+(s1),a-(s2)=a+(s2),
where a- and a+ denote the area function to the left and right of the points s1 and s2, respectively. Furthermore in the absence of concentrated loads as required by the continuity of the uncertain and deterministic loads, the shear force V(x)=(a(x)y′′(x))' on the beam will also be continuous. From the optimality condition (4.12), it follows that V(x)=βa'(x) . Thus we have the further continuity conditions(4.18)V(s1)=a-′(s1)=a+′(s1),V(s2)=a-′(s2)=a+′(s2),
where a-′ and a+′ denote the derivatives of the area function to the left and right of the points s1 and s2 , respectively. The case when y(x)<2μ1fmax for x in a finite interval will be solved in the example problems. The uncertain functions f(x),m0 and m1 can be computed from equations (4.3), (4.10), (4.11) and (4.14). In particular the uncertain loading f(x),s1≤s≤s2 is given by(4.19)f(x)={β4μ1x(x-1),fors1≤x≤d1,fmax,ford1≤x≤d2,β4μ1x(x-1),ford2≤x≤s2,
where d1 and d2 are unknown locations to be determined from the continuity conditions(4.20)f-(d1)=f+(d1)=fmax,f-(d2)=f+(d2)=fmax,
where f(x)- and f(x)+ denote the uncertain load functions to the left and right of the points d1 and d2, respectively. From equations (4.3), (4.11) and (4.14), it follows that m0=m1=η/2.It is noted that the number of unknowns equals the number of equations resulting in unique solutions. This aspect the method of solution will be illustrated in the next section by applying the technique to several problems of practical interest.To assess the efficiency of the optimal designs, comparisons are made with uniform beams under uncertain loads for whicha(x)=1 for 0≤x≤1. The antioptimality condition (4.10) applies to this case also and consequently the differential equation for a uniform beam under worst case loading becomes(4.21)d4ydx4={G(x),for0≤x≤s1,y2μ1+G(x),fors1≤x≤d1,fmax+G(x),ford1≤x≤d2,y2μ1+G(x),ford2≤x≤s2,G(x),fors2≤x≤1.
The solution of the differential equation (4.21) subject to the boundary conditions (4.4) and the constraints (4.3) gives the deflection yun(x) of a uniform beam under worst case loading. The efficiency of the design can be determined by comparing the maximum deflections of the uniform and optimal beams, namely,(4.22)Ieff=ymaxyun×100%,
where Ieff is the efficiency index in percentage, yun and ymax are the maximum deflections of the uniform and optimal beams under worst case of loadings.


## 5. Applications of Method

Example 5.1 (unconstrainedF(x) with 0<s1<s2<1).
Let the beam be subjected to only the uncertain transverse loadF(x) given by equation (4.1) with 0<s1<s2<1, that is, no uncertain moments are applied on the boundaries so that m0=m1=η=0, and there is no deterministic load applied, that is, G(x)=0. Moreover it is set equal to fmax=∞. For this case the optimal area function satisfying the moment boundary conditions in equation (4.4) can be computed from equations (4.16) as
(5.1)a(x)={c1x,for0≤x≤s1,148μ1x3(x-2)+c2x+c3,fors1≤x≤s2,c4(1-x),fors2≤x≤1.
Equations (4.20) for a(x) contain six unknowns β,μ1,c1,c2,c3, and c4 which are computed from six equations for the volume and L2 norm constraints (4.3), and the continuity conditions (4.17) and (4.18). These constants in terms of s1 and s2 are given by
(5.2)μ1=140s25-140s15-116s24+116s14+124s23-124s13,c1=5L(3s13-8s12+3s2s12+6s1-8s2s1+3s22s1+6s2-8s22+3s23),c2=5(3s14-4s13+8s23-6s22-3s24)K,c3=-5s13(3s1-4)K,c4=-5(3s13-4s12+3s2s12-4s2s1+3s22s1-4s22+3s23)L,β=±P2,
where
(5.3)L=6s14-15s13+6s13s2+10s12-15s2s12+6s22s12+10s2s1-15s22s1+6s1s23+10s22-15s23+6s24,K=-6s25+6s15+15s24-15s14-10s23+10s13,P=15s25-15s15-12s24+12s14+13s23-13s13.
A numerical example is given for the case s1=0.2 and s2=0.8 for which β=-0.08584,μ1=0.003684,c1=4.479,c2=5.655,c3=-0.1538, and c4=4.479. The optimal area function a(x) and the antioptimal F(x) (worst case loading) are shown in Figure 2.In the case of a uniform beam, the worst case loading is given by f(x)=5.825x(1-x) for s1≤x≤s2 and the corresponding deflection is y=0.04292x(1-x). In this case yun=0.01452,ymax=0.01073 and the efficiency is 74% as determined by the efficiency index given by equation (4.22).Curves of optimala(x) and worst case loading F(x) are plotted against x for s1=0.2 and s2=0.8 (Example 5.1).
(a)(b)Example 5.2 (unconstrainedF(x) with s1=0,s2=1 and deterministic loading).
Let the beam subject to only the uncertain transverse loadF(x) given by equation (4.1) with s1=0,s2=1, and the deterministic load G(x)>0. Moreover, fmax=∞. Solution for the worst case loading is given by (4.19) with s1=d1=0,s2=d2=1, that is, F(x)=(β/4μ1)x(x-1) for 0≤x≤1. The differential equation for the area function takes the form
(5.4)a′′(x)=14μ1(x2-x)+1βG(x),0≤x≤1,
where the optimality condition (4.12) is used. Let G(x) be specified as G(x)=sin((π/2)x). Then the solution of (5.4) satisfying the volume constraint (4.3) and the boundary conditions (4.4) is given as
(5.5)aopt(x)=148μ1x3(x-2)-4βπ2sin(π2x)+(4βπ2+148μ1)x,
where β=-0.1467 and μ1=-0.006694. The optimal area function aopt(x) and the worst case loading F(x) are shown in Figure 3.In the case of a uniform beam, the worst case loading is given by f(x)=5.479x(1-x) for 0≤x≤1, and the corresponding deflection is y=0.07335x(1-x). In this case yun=0.02321,ymax=0.01833, and the efficiency is 79% as determined by the efficiency index given by equation (4.22).Figure 3
Curves of optimala(x) (solid line) and F(x) (broken line) are plotted against x (Example 5.2).


## 6. Conclusions

The optimality conditions were derived for robust shape design of beams subject to unknown loads with the moment of inertia related to the area function asI=an,n=1,2,3,  and the Lp,1<p<∞ norm of the transverse load subject to an equality constraint. The potential energy of the beam was specified as the performance index which was minimized with respect to area function (optimization) and maximized with respect to load function (antioptimization). Lagrange multiplier method was used to take various constraints into account and the slack variable method was used to take the upper bound on the transverse load into account. In general the optimality conditions are expressed as a system of coupled nonlinear differential equations in terms of deflection, area and load functions necessitating the use of numerical methods. However, the special case of n=1 and p=2 was solved explicitly to illustrate the method, to provide benchmark solutions and to assess the efficiency of the designs.Analytical solutions were obtained for the special case under consideration and several example problems were solved involving various cases of loadings. Numerical results were given for the optimal area and antioptimal load functions. The efficiencies of the optimal designs were computed in terms of the maximum deflections of the optimal beam and the uniform beam under least favorable loading. It was shown that for the cases studied the design efficiency can exceed 70%. Load uncertainties arise due to the unpredictable conditions occurring under operational conditions and are usually the cause of unexpected failures. This situation indicates the importance of a robust design since a design optimized for a known load is strong only for this load, but it will be weak if the loading conditions happen to change.



---

*Source: 841017-2009-12-31.xml*