16.21 Techniques of Structural Analysis and Design Spring 2005 Unit #9 Calculus of Variations
16.21 Techniques of Structural Analysis and Design Spring 2005 Unit #9 Calculus of Variations
16.21 Techniques of Structural Analysis and Design Spring 2005 Unit #9 Calculus of Variations
(1)
v a b
Consider a function of the following form: F = F (x, u(x), u (x)) It depends on an independent variable x, another function of x (u(x)) and its derivative (u (x)). Consider the change in F , when u (therefore u ) changes: F = F (x, u + u, u + u ) F (x, u, u ) = F (x, u + v, u + v ) F (x, u, u ) expanding in Taylor series: F F 1 2F 1 2F 2 v + v + ( v ) + (v )(v ) + F 2 u u 2! u 2! uu F F = v + v + h.o.t. u u First total variation of F: F F = lim 0 F (x, u + v, u + v ) F (x, u, u ) = lim 0 F v + F v F F u = v + v = lim u 0 u u F F F = u + u u u F = F + 2
First variation of a functional: I = F (x, u(x), u (x))dx = F (x, u(x), u (x)) dx F F I = u + u dx u u Extremum of a functional u0 is the minumum of a functional if: I (u) I (u0 )u A necessary condition for a functional to attain an extremum at u0 is: I (u0 ) = 0, or dI (u0 + v, u0 + v ) =0 d =0
Note analogy with dierential calculus. Also dierence since here we require dF = 0 at = 0. d b F F I = u + u dx u u a Integrate by parts the second term to get rid of u . b F d F d F I = u + u u dx u dx u dx u a b F F d F b = udx + u u dx u u a a Require u to satisfy homogeneous displacement boundary conditions: u(b) = u(a) = 0 Then: b F d F I = udx = 0, u dx u a 4
u that satisfy the appropriate dierentiability conditions and the homoge neous essential boundary conditions. Then: F d F =0 u dx u These are the EulerLagrange equations corresponding to the variational problem of nding an extremum of the functional I . Natural and essential boundary conditions A weaker condition on u also allows to obtain the Euler equations, we just need: F
b u
= 0 u a which is satised if: u(a) = 0 and u(b) = 0 as before u(b) = 0 and
F
(a) u F (a) u F (b) u
=0
=0
= 0, or u = u0 on Su Essential boundary conditions: u
Su F Natural boundary conditions: u = 0 on S . Example: Derive Eulers equation corresponding to the total po tential energy functional = U + V of an elastic bar of length L, Youngs modulus E, area of cross section A xed at one end and subject to a load P at the other end. L EA
du
2 (u) = dx P u(L) 2 dx 0 Compute the rst variation:
EA du
du
2 = dx P u(L) 2 dx dx