泛函極值及變分法

第二章泛函極值及變分法(補充內(nèi)容)變分的基本概念泛函和變分泛函是一種廣義的函數(shù)，是指對于某一類函數(shù)｛y(%)｝中的每一個函數(shù)y(%)，變量J有一值與之對應，或者說數(shù)J對應于函數(shù)y(%)的關系成立，則我們稱變量J是函數(shù)y(%)的泛函，記為J[y(%)]。例1：如果表示兩固定端點A(%A,yA)，B(%B，yB)間的曲線長度J(圖2.1.1)，則由微積分相關知識容易得到：(2.1.1)J=I%B\,1+(dy/d%)2(2.1.1)%A顯然，對于不同的曲線y(%)，對應于不同的長度J，即J是函數(shù)y(%)的函數(shù)，J=J[y(%)]。圖2.1.1兩點間任一曲線的長度例2：歷史上著名的變分問題之 最速降線問題，如果2.1.2所示。設在不同鉛垂線上的兩點P1與P2連接成某一曲線，質(zhì)點P在重力作用下沿曲線由點P1自由滑落到點P2,這里不考慮摩擦作用影響，希望得到質(zhì)點沿什么樣的曲線滑落所需時間最短。圖2.1.2圖2.1.2最速降線問題選取一個表示曲線的函數(shù)y(%)，設質(zhì)點從P1到P2沿曲線y=y(%)運動，則其運動速度為:dtd%dsv=dtd%dt其中，S表示曲線的弧長，t表示時間，于是:,,"+y'2,

dt= d%v設重力加速度為g，則y=京。因為P1和P2點的橫坐標分別為11到％,那么質(zhì)點從P1到P2所用時間便為：J[y(1)]=j12；%1+y dx%<2g[y(11)-y(1)“211+[y'(1)]2廠d1 (2.1.2)11[2g[y(11)-y(1)]J則最速降線問題對應于泛函J[y(1)]取最小值?；仡櫤瘮?shù)的微分：對于函數(shù)的微分有兩種定義：一種是通常的定義，即函數(shù)的增量：Ay=y(1+A1)-y(1)=A(1)A1+p(1,A1) (2.1.3)其中A(1)與A1無關，且有A1-0時p(1,A1)-0,于是就稱函數(shù)y(1)是可微的，其線性部分稱為函數(shù)的微分dy=A(1)A1=y'(1)A1，函數(shù)的微分就是函數(shù)增量的主部。函數(shù)微分的另外一種定義：通過引入一小參數(shù)￡，對y(1+8A1)關于￡求導數(shù)，并令《一0的途徑得到，即：6y(二 [0=y，(1+sA1)A110=y，(1)A1=dy (2.1.4)上式說明y(1+￡A1)在￡=0處關于￡的導數(shù)就是函數(shù)y(1)在1處的微分。相應地，在泛函J[y(1)]中，變量函數(shù)y(1)的增量在其很小時稱為變分，用6y(1)或6y表示，指y(1)與它相接近的y1(1)的差，即：3y(1)=y(1)-y1(1)。泛函的變分也有類似的兩個定義：對于函數(shù)y(1)的變分6y(1)所引起的泛函的增量為AJ=J[y(1)+5y(1)]-J[y(1)]，當5y(1)-0時泛函增量的線性主部就稱為泛曲在函數(shù)y(1)處的變分，記為6J，即：5J=J[y(1)+5y(1)]-J[y(1必-0=4y(1)，5y(1)] (2.1.5)其中L[y(1),6y(1)]是泛函增量的線性主部，而且其對于變分6y(1)是線性的。另一種定義：拉格朗日的泛函變分定義為：

泛函變分是J[y（%）+￡5y（x）]對g的導數(shù)在￡=0時的值，即:—a.5J=—J[y(x)+g5y(x)] =L[y(x),5y(x)]洗 ￡f0首先，我們進行泛函:J=J[y(x)]=jx2F(x,y(x),y'(x))dxx1的變分。此泛函的增量可以用Taylor展式表示為:二jx2|aFAy+空Ay，+±xj行 行'C1a2F/,、- (Ay)2+2a2 7a2f, a2f . AyAy+ (Ay)2+a二jx2|aFAy+空Ay，+±xj行 行'C1a2F/,、- (Ay)2+2a2 7a2f, a2f . AyAy+ (Ay)2+ayay' a(y，”dx(2.1.8)當Ayf0，上式積分中的前兩項是增量的線性主部，后面的項為高階無窮小量。根據(jù)變分的定義，該泛函的變分為:5J=5J=jx2x1fafsafsX

5y+5ydx〔ay'ay，I(2.1.9)（2.1.9）也稱為泛函j的一階變分，而（2.1.8）式的后三項為二階變分，記作8J，即:52j=jx2x1a2f&、.a252j=jx2x1a2f&、.a2fss,a2f^J--(5y)2+t^75y5y+、八(5y)2dxay2 ayay a(y”也可以通過拉格朗日泛函變分的定義，得到:5J=——J[y(x)+s5y(x)]0=jx2x1=；（I$y+Fy助此結果與（2.1.9）是相同的。a—F(x,y+g5y,y+￡&y)asdx￡f0類似地，如果泛函的值決定于兩個函數(shù)，并且這些函數(shù)是兩個變量的函數(shù)，如:J=J\u(x,y),v(x,y)]=jF(x,y,u,v,u,u,v,v)ds其變分為:af, af, af, af, af,——5u+ 5v+——5u+——5u+ 5vx y xaf+ar6vydsy依此類推，不難得到多個多元函數(shù)的變分。此處，泛函的變分滿足下面的一些運算規(guī)律：6{j[y(%)]+J[y(x)]}=3J[y(x)]+6J[y(x)] (2.1.14a)1 2 1 26{J[y(%)].J[y(%)]}=6J[y(%)]J[y(%)]+J.[y(%)]-6J[y(%)] (2.1.14b)J[y(%)]] 6J[y(%)1J」y(%)]—J[y(%)16J」y(%)] 0114、(3)6,t>=—1 2-z1- 1-： 2 (2.1.14c)[J[y(%)]J J[y(%)]J222(4)6{J[y(%)]}"=n{J[y(%)]}，t.6J[y(%)] (2.1.14d)2.1.2泛函的極值和變分問題本節(jié)將討論泛函的極值和變分。微積分知識：函數(shù)取極值的必要條件(但不是充分條件)：對于一個連續(xù)可導函數(shù)，如果其在定義域的某(些)點函數(shù)有極值，那么這個函數(shù)的一階導數(shù)在這(些)點等于零，這個(些)點就是函數(shù)的極值點或駐點。對于泛函的極值問題，也有類似的結論，即泛函取極值的必要條件是其一階變分6J=0。簡要證明:假設函數(shù)y(%)是泛函J所定義的函數(shù)集合中的任一函數(shù)，這里不妨設泛函J[y(%)]在函數(shù)y(%)處有極大值，那么對于任一實變量a，必有:J[y(%)]>J[y(%)+a6y(%)]令f(a)=J[y(%)+a6y(%)]，則有:f(a)L°=J[y(%)]”[y(%)+a6y(%)]=f(a)上式表示f(a)在a=0處有極大值，根據(jù)函數(shù)取極值的必要條件:df(a)?ddf(a)?da 3df(a)?daQ=o=0，得至U：dJ[y(%)+a6y(%)]a=0=6J=0(2.1.17)由此就得到泛函取極大值的必要條件是其一階變分為零。同樣的方法可以證明，泛函取極小值的必要條件也是其一階變分為零。泛函實現(xiàn)局部極大或極小值的充要條件：泛函實現(xiàn)局部極大或極小值的充要條件與函數(shù)取極值的充要條件類似，除了其一階變分為零外，還需要考察二階變分的情況:1)若泛函J[y(%)]在y(%)處取局部極大值，其充分必要條件為:5J\y(%)]=0, 82J\y(%也0 (2.1.18)2)若泛函J[y(%)]在y(%)處取局部極小值，其充分必要條件為：5J[y(%)]=0, 52J[y(%)]>0 (2.1.19)通常，我們將求泛函極值的問題稱為變分問題。變分法的基本預備定理：如果函數(shù)F(%)在線段(%1,%2)上連續(xù)，且對于只滿足某些一般條件的任意選取的函數(shù)5y(%)，有：』%2F(%)5y(%)d%=0 (2.1.20)%1則在線段(%p%2)上有：F(%)=0 (2.1.21)這里5y(%)滿足的一般條件為：①一般或若干階可微；②在(%p%2)的端點外為0；③5y(%)|<￡或5y(%)|<￡和5yr(%)|<￡等。對于多變量問題，也有類似的變分定理。二維：函數(shù)F(%，y)在(％，y)平面S內(nèi)連續(xù)，設5u(%,y)在S的邊界上為零，5u|<￡,5u|<￡,5uj<￡且滿足連續(xù)性以及一階或若干階的可微性，對于這樣選取的5u(%,y)，若有:』sF(%,y)5u(%,y)d%dy=0 (2,1.22)則在區(qū)域S內(nèi)有：F(%,y)=0 (2.1.23)現(xiàn)在我們來研究最簡單的泛函：J=J[y(%)]=)％2F(%,y(%),yr(%))d% (2.1.24)%1的極值問題。其中F為%，y和y，的函數(shù)，且F(%,y,y')是三階可微的。確定泛函極值的曲線y=y(%)的邊界是固定不變的，且有：y(%1)=y^y(%2)=y2 (2.1.25)

采用拉格朗日法來求其泛函變分，有:(2.1.26)J[y+￡3y]=Jx2F(x,y+￡3y,y'+￡§y')dxx1(2.1.26)令：y=y+￡3y,y，=y，+￡3y1令：（利用復合函數(shù)求導法則）—J[y+—J[y+￡3y]=fd d一F(x,y+￡3y,y，+￡3y，)3y+ F(%,y+￡3y,y'+￡3y')3y'dxxi dydyr(2.1.27)令￡T0，則:3J= J[y+￡3y] =Jx2O￡ ￡t0 3J= J[y+￡3y] =Jx2O￡ ￡t0 x1OF OFOyy+Oy'其中:OF O OF O,機=OyF(x,y,y),OyF(x,y,y)（利用分部積分）fb.jfb.fbfb.udv=Jd(uv)一vdu=uvJ-Jvdu.a a aaaJ;Idxlay7)3ydx3ydx廠上式中利用到了固定邊界條件3y（x1）=3y（x2）=0。最后，可得到變分的極值條件:3J=fx2x13ydx=0根據(jù)變分法預備定理，得到:d(OF、上式中，關于x的導數(shù)為全導數(shù)，「「(，八dEaEdx0EdyOEdy'

E=E(x,y',y'),——= + -+ —、 dx 0xdxOydxOy'dx,d(OF、(2.1.32)O2FdxO2FdyO2Fdy'

d(OF、(2.1.32)OxOy'dx OyOy'dxOy'2dx（2.1.31）式即為著名的歐拉方程，是歐拉于1744年得到的，也稱為歐拉一拉格朗日方程。a）F不顯含y時（2.1.31）式經(jīng)過一次積分得一階微分方程:a）F不顯含y時（2.1.31）式經(jīng)過一次積分得一階微分方程:Fy=Fy=CC是積分常數(shù)。b）F不顯含x時,（2.1.31）式做如下變化:由（2.1.32）式:Ftiy"+F,y，+Fb）F不顯含x時,（2.1.31）式做如下變化:由（2.1.32）式:Ftiy"+F,y，+Fyy，—Fy=0.因此:經(jīng)過一次積分得:y〃「y'Ry'Fy,_F=C(2.1.34)C是積分常數(shù)。例2.1.2（續(xù)）：最速降線問題11+y'2—, dx=0.\：2gy1+y'2解：F(x,y,y')=—<2gy不顯含x，按式（2.1.34），有:y，.2y'1+y'2y，2y，21+y，2v2gyv2gy(1+y，2)v；2gy此式簡化為:y(1+y/2)=c，其中c=-1—>01 12gc2引入?yún)?shù)y,=ctg。，則有：y=——y=——c^-—1+Ctg20csin20=1Cc(1-cos20).21由于包=ctg0dx因而：dx=d^-ctg因而：dx=d^-ctg02csin0cos0d0

1 ctg02jsin20d0=c.(1-cos20)d0.積分后得：x=c(0-2-sin20)+c2=c2-(20-sin20)+c2.由初始條件：y(0)=0,知：c2=0.于是最速降線問題的解為：c -x=T(20-sin20).《cy=才(1-cos20).2其中c1由邊界條件y(x1)=y1來確定。再令：①=20,R=—12就得到：x=R(①-sin①)y=R(1-cos①)從解析幾何知，上述方程是擺線的參數(shù)方程，因此最速降線是半徑為R的圓沿x軸轉(zhuǎn)動時圓周一點所描出的曲線中的一段。例2.1.3求泛函J[y]=f(y'2-2ycosx)dx在條件y(0)=0,y(兀)=0下的極值曲線。0解：F(x,y,y))=y'2-2ycosx.貝U：F=—2cos%,F=2y，—(二)=2y〃.y y' dxSy，由空——(至)=0，對應的歐拉方程為:Sy dxSy'-2cosx-2y"=0ny"+cosx=0., ，—一一2代入邊界條件得：c「；,c2=-1.???y=cosx+-x-1.兀較復雜的泛函的歐拉方程可以仿照上述方法導出。a)對于取決于一個自變量的幾個函數(shù)的泛函J[y1(x),y2(x),...yn(x)]=fbF(x;y1,y2...yn;yJy；...y「)dx.a泛函J[y1(x),y2(x),...yn(x)]的變分問題對應于下列歐拉方程組:(2.1.35)退--(遲)=0, (i=1,2...n)(2.1.35)Sy dxSy，i i打例2.1.4：求泛函J[y(x);z(x)]=)2(y&+z3+2yz)dx在邊界條件：0y(0)=0, y(—)=1, z(0)=0, z(—)=-1下的極值曲線。解：F(x;y,z;y',z')=y'2+z'2+2yz、，"一z=0則有歐拉方程組：yz”-y=0消去z，得方程：y⑷-y=0由此解出：y=cex+ce-x+ccosx+csinx再由z=y”得：z=c1ex+c2e-x-c3cosx-c4sinx利用邊界條件：c=c=c=0,c=112 3 4因而極值曲線為：y=sinx,z=一sinxb)對于泛函取決于y(x)及其n階導數(shù)的情況"b" ，”J=F(x,y,y',...,y(n))dx.

(2.1.36)其歐拉方程：Fy-—^y'+-^—Fy"+...+(-1)n——Fy(n)=-0dxdx2 dxn(2.1.36)它的通解含有2n個任意帶數(shù)，它們由2n個邊界條件來確定。… 一門A 例2.1.5：J[y]=一一py〃2+pydx,(mp為常數(shù))-1\2 )邊界條件：y(-1)=0,y'(-1)=0,y(1)=0,y'(1)=0.求泛函J[y]的極值曲線。解：相應的歐拉方程為：股(4)+p=0由此得至U：y=———x4+cx3+c^x2+cx+c24日 1 2 3 4利用邊界條件得到：y=--p-(x2-12)2.24日c)對于泛函取決于多元函數(shù)的情況J=BIF(x,y,z;u;u,u,u)dxdydz=0V其對應的歐拉方程為：(2.1.37)例2.1,6(2.1.37)例2.1,6：設泛函J[u]=』)Qdxdy解：歐拉方程：alcFu解：歐拉方程：alcFu=0,Fu=2u,—Fu=2u.Fu=2u,萬-Fu=2u則:a則:a2u+a2u=0ax2ay2例2.1.7：例2.1.7：設泛函J[u]』Q其對應的歐拉方程是泊松方程:dxdy.10

10a2M a2u布+蚤=fx，y).可動邊界的變分問題，變分問題中的邊界條件所謂可動邊界是指極值曲線（或曲面）的兩個端點或其中一個端點（或邊界）并不通過預先給定的點（或邊界）。在我們前面所討論的泛函J」%2F（%,y,y，）d%在邊界條件y（%）=八和y（%2）=y2下的極值問%1題，此時的固定邊界條件稱為幾何邊界條件或稱為強加邊界條件。所謂“強加”，是指這些邊界條件是在變分問題中預先強加上去的。如果我們在求泛函J=52F（%,y,y）d%的極值問題時，端點%1和％2的值均不給定。泛函%1取極值的必要條件依然是：dyd%(aydyd%(ay，/aF%=%--2 ay'5y%=%,=0(2.1.42)與固定邊界變分不同的是，這里的8y在端點處并不總是為零，可以為任意的。這樣，由極值條件5J=0除了可得到歐拉方程（2.1.31）外，還有:運%=%運%=%=0,ay 1%=%=0

2(2.1.43)上式由變分得出的條件稱為自然邊界條件?？梢钥闯鲞@樣的邊界條件不是預先給定的，而是從變分原理的5J=0自動導出，它是保證極值存在而必須滿足的條件。在力學問題中，無約束時變分原理將自動補充邊界處所缺的力學邊界條件，因而自然邊界條件往往表現(xiàn)為力學邊界條件。a）每個函數(shù)端點分別在直線％=。和%=b上.泛函j的極值函數(shù)除了要滿足歐拉方程外，還應滿足：Fj=F,| =0. （2.1.44）b）更一般情況，如果泛函J[y]」bF（%,y,y'）d%中,函數(shù)y的端點（。％）與（by）分別在曲線a叭%,y）=0與甲（%,y）=0上移動，泛函J的極值函數(shù)除滿足歐拉方程外，還要滿足橫截條件（transversalitycondition）o11(2.1.45)F-y'Fx(2.1.45)F-y'Fx x=bF,yx=b這里a與b本身是待定參數(shù)。泛函的條件極值有些變分問題，容許函數(shù)有時還會受到附加約束條件的限制，這就是條件極值問題。對于這種極值問題，可用類似于處理多元函數(shù)的條件極值的Lagrange乘數(shù)法，把范函條件極值問題轉(zhuǎn)化為無條件極值問題。定理(Lagrange)：略這個方法還可以推廣到等周問題，即有如下定理：歐拉定理：略例2.1.8：等周問題在平面上，給定長度為l的所有封閉光滑的曲線中，求一條曲線，使它所圍成區(qū)域的面積A最大。設所求曲線上的參數(shù)方程為:且x(10)=x(11),y(10)=y(ti).約束條件：J%qx2+y2dt=l(常數(shù))10格林公式：若函數(shù)X(x,y)和F(x,y)及其一階偏導數(shù)在閉區(qū)域D上連續(xù)，則有:12

dxy=EXd;+Ydy.L令y=%,X=-y,其中L是區(qū)域D的邊界，且積分沿L的正方向（即逆時針方向）。由格林公式，曲線l所圍成的面積:t0t0A=2J;dy-yd;=2于是等周問題可歸結為求泛函:J=』t1(%y-y;)dt在等周條件或（1）下的極大值。作輔助泛函:I=2Jt11;y-y;+入《;2+y2)dtt0由:H由:H--Hy==0ydtH——H%'=0其對應的歐拉方程組為:y--dt-y+?q;2+y2J積分后得:C 入；=2ci=2ci、：；2+y2V%2+y2整理后:這是圓族方程，令:13入=cost2 (0WtW2n)y一cy一c12代入等周條件，得:e=J2''—sin21+—cos21dt=J2無—dt=九no4 4 o2、2即：X=-、2于是G-c1+(y一c)=—2 1 12n利用邊界條件,x(10)=x(11)y(t0)=y(t1)可定出J,c2,故所求極值曲線是一個圓。142.2力學中的變分原理在力學中，我們有各種各樣的原理，諸如能量守恒原理、動量守恒原理、達朗伯原理、虛位移原理、哈密頓原理等等。而作為古典力學基礎的著名的牛頓運動定律實質(zhì)上也是原理。原理可以被分為兩類，即非變分的原理和變分的原理。非變分的原理直接研究真實的運動；而變分原理則不然，它不是專注于實際的運動，而是考察一定約束條件下所容許的一切可能的運動，從中挑選出實際實現(xiàn)的一種真實運動來。如果說非變分的原理提供的是各種各樣普通的函數(shù)關系，那么變分原理應該是考察相應于各種運動狀態(tài)的某些特征量（泛函）并取極值（通常對應于真實運動），這便是我們所熟知的變分的含義。由此可以看出，變分原理是在縱觀全局的基礎上更一般地來論述運動的，較之非變分的原理進行了更多的概括與抽象。這樣說并不是貶低非變分的原理的重要性，事實上，很多變分的原理和非變分的原理在一定條件下都是可以互相推導或是等價的，只是各種原理的表述方式不同，因而在不同場合下應用時方便程度不同罷了。力學原理又可分為微分形式的表述和積分形式的表述。前者適用于運動的每一瞬時以及任意局部點，而后者適用于有限的時間間隔以及有限區(qū)域內(nèi)。在力學的諸多原理中，虛功原理是最基本的，其他的若干原理可從它得到。下面，我們首先介紹虛功原理。虛功原理虛功原理亦稱虛位移原理。在分析力學中，由質(zhì)點系組成的力學體系的虛功原理是熟知的。對于一個由N個質(zhì)點組成的質(zhì)點系而言，如果考慮的是靜平衡問題，則有分析力學的虛功原理：2NF5r=0 （2.2.1）iii=1其中F（（i=1，2,…，N）是作用在質(zhì)點系上的給定力，包括非理想的約束力等；。r.（i=1，2,…，3N）是質(zhì)點系滿足約束的任意一組無限小虛位移矢量。進一步，如果作用在質(zhì)點系上的諸力均是有勢的，亦即對于諸力F存在勢函數(shù)V，使得F=-aV/6r（i=1，2,…，N），I I I則上述的虛功原理可轉(zhuǎn)化為最小勢能原理。在靜止的平衡力學系統(tǒng)的所有容許位移中，真實的位移使勢能的變分為零，即。V=0。虛功原理指出，系統(tǒng)平衡時的位置是指系統(tǒng)可能有的一切位置（對應各種虛功值）中的這樣一種位置，此時作用力所作虛功之和為零。這樣，從系統(tǒng)可能有的一切運動狀態(tài)中確實挑選出了平衡這樣一種實際實現(xiàn)的運動狀態(tài)。作為泛函的虛功取極值虛功為零）時對應著真實15

的運動(平衡狀態(tài))。下面我們給出彈性連續(xù)體的虛功原理表述。設彈性體在體力fx，fy，fz以及表面為Fx，F(xiàn)的運動(平衡狀態(tài))。下面我們給出彈性連續(xù)體的虛功原理表述。設彈性體在體力fx，fy，fz以及表面為Fx，F(xiàn)y，F(xiàn)z作用下處于平衡。以｛，「，*\，7，7｝表示任一點處的應力分量，則在彈性體內(nèi)有平衡方程:So 。工 St——+ xy+—xrSx Sy SzSt So St——xy+ +——yzSx Sy SzSt St So―xr+ yz+ zSx Sy Sz+fx=0以及在應力邊界r,上，滿足力學邊界條件:X=Fx，Y=F，Z=Fz其中：X=lo+mt+ntz，Y=lt+mo+ntz，Z=ltz+mtz+noz這里｛l，m，n｝表示彈性體表面上一點的外法線方向余弦。我們假定彈性體平衡時的真實的位移為｛u，v，w｝，從這個平衡位置對物體施加一組任意的無限小虛位移｛腦，入，Jw｝，于是便有:-JJJVSo Sr Stx+ xy+xzSx Sy Sz（St So StJJh-F1u+Y-F}v+（Z-F^hw11.Sr Sr So x三+苫=+fjW1v+rt x y z（2.2.5）其中dv，ds分別表示彈性體的體積元和面積元。這里虛位移的選擇應滿足另一部分位移邊界ru上的幾何條件，即:（2.2.6）5u=0,5v=0,5（2.2.6）利用高斯公式，并經(jīng)過分部積分，（2.2.5）可進一步化簡為:JJJC3s+o3s+o3s+t3y+t3y+t3y]v（2.2.7（2.2.7）-JJJf3u+f3V+f3wdiv-JJF3u+F3v+F3w2s=0yxyz 「txyz3sxShu3syS3VSyS3sxShu3syS3VSyS3w3s= ,3yxy遮+迎Sy Sx,3yyzS3vS3w

+ 3yxzS3wS3u

+ （2.2.8）為虛應變，(2.2.7)即為彈性體虛功原理。按照彈性力學定義:為虛應變，(2.2.7)即為彈性體虛功原理。按照彈性力學定義:16u=1川Q￡+o￡+o￡+tY+ty+ty>dv （2.2.9）2Vxxyyzzxyxyyzyzxzxz稱為彈性體的變形能。因此（2.2.7）中的第一項即為虛變形能，第二項（取正號）為體積力所作的虛功，第三項（取正號）為表面力所作的虛功。由于上述過程是從平衡位置施以虛變形，故虛功簡單地表示為力與虛位移之乘積，并無因子1/2,這是虛功有別于真實功的主要特點。將式（2.2.7）進行移項，不難看出：在任一虛位移過程中，外力作的總虛功等于彈性體的總虛變形能。上述推導說明虛功原理是物體在外力作用下并滿足一定的幾何邊界條件而處于平衡的必要條件。相反的推導過程，我們完全可以利用虛位移的｛6u，6v，6w｝的任意性而得到力學平衡方程（2.2.2）以及力學（自然）邊界條件（2.2.3）。這說明虛功原理同時也是彈性體平衡及力學邊界條件的充分條件。虛功原理是彈性力學中的變分原理的基礎，其在有限元法中也具有極其重要的應用價值。盡管我們是從彈性平衡的角度給出了虛功方程，但是一般說來，虛功原理具有普遍意義，它可以適用于一切結構，不論材料是線性還是非線性，也不論物體的變形是彈性或非彈性。2.2.2最小勢能原理上節(jié)所介紹的虛功原理對于任何應力-應變關系的結構均成立，不論是彈性或是非彈性Y,YY,Y,Y ｝分別表示彈性體內(nèi)一點xyyzxz必存在一個正定的狀態(tài)函數(shù)令t,O,O,T,t,T卜口L,￡,￡,的應力和應變分量，在小變形情形下，U0=U0Cx，￡y,￡z,Yxy,Yyz，Yxz）使得：dUdU0=od￡+od￡+od￡z+tdy+tdyz+tdyz（2.2.10）這里，U0稱為應變能函數(shù)或應變能密度。狀態(tài)函數(shù)U0是單值的函數(shù)，因而dU0是全微分，有：。U。U。。U。U。U。U5U5UO=——0-,O=——0-,O=——0-,T=——0-,T=——0-,T=——0x%y洗z洗xy 5yyz 5yxz 5yxyzxyyzxz這樣，虛功原理（2.2.7）就變?yōu)?（2.2.11）1711Jf5w+f5v+f8w1v+JJFdu+FBv+F8w1sxyz ^txyzau.sY

+——o8ydvay au.sY

+——o8ydvay xzJxz（2.2.12）=I——o8s+——o8s+——o8s+——o8y+——o8yyIasx as y as zay xyay yzx y z xy yz=JJJV8U0dv=8JJJVU0dv=8U上式中U=fflVU0dv為彈性應變能。進一步，如果作用于彈形體上的體力和表面力均為有勢力，即存在勢函數(shù)①（u,v,w）和中（u，v，w），使得：（2.2.13ab）-8①=f8u+f8v+f8w,一8%=F8u+F8v+F8w（2.2.13ab）從而：JJJV§8u+f8v+f8w1v+Jj8u8u+F8v+F8w4s=-8JJJ①dv-8II乎ds=-8 ①dv+JJ+ds=-8VV rt V rt這里，V表示外力勢能，則（2.2.12）亦表示為:（5（U+V）=0如果定義n=U+V為系統(tǒng)的總勢能，故6n=0。此式或（2.2.15）稱為勢能駐值原理，即在滿足已知幾何邊界條件的一切容許位移u，v，w中，真實的位移使得系統(tǒng)總勢能（泛函）取極值。2.2.3虛余能原理前面兩節(jié)中所介紹的虛功原理及其應用于彈性連續(xù)體而得到的最小勢能原理，都是以位移作為未知函數(shù)的，位移一旦求得，根據(jù)幾何關系式和應力應變關系式不難得到相應的應變和應力分量。但是，在很多工程實際問題中往往也需要直接以應力作為待求的未知函數(shù)，尤其在以應力為目標的近似解法中，如果依舊沿用先求位移而后通過微分求應變再得到應力的方法，勢必會影響應力解的精度。實際應用的需要自然應運而生了相應的虛余功原理及最小余能原理。下面的討論仍以彈性連續(xù)體為例，而且其應力應變可呈各種關系。設彈性體在已知體力以及給定的邊界條件下處于平衡，｛u，v，w｝和L,s,s,y,y,y｝分別表示彈性體xyz.xyyz.xz.內(nèi)一點處的位移分量和應變分量。因此，在彈性體內(nèi)，有:au av aw八s =0,s =0,s =0xax yay z az187x「=0,=0,7xz」=0(2.2.16)以及在邊界上:u一u=0,v—v=0,(2.2.17)rR步，我們設平衡時的應力狀態(tài)為｛%，Gy，%，Txy，7｝，并假定物體從這個平衡狀態(tài)接受一組任意的、無限小的虛應力（應力的變分）｛Jj空，J,，抗yz，”3于是有:IIIV￡--Vu|3q+￡7xy-It+VX3Txy+VwV|3oVzJZ3Tyz+丫Ifb—ubF+(v-vbF+&-wbFIs3Txz迅ru(2.2.18)其中｛鞏，6Fy,JFz｝是表面力相對應于虛應力的虛變化。新的應力分量應該不違背彈性連續(xù)體的平衡方程和力學邊界條件，如下:乜+5c)+?(+3txyVy)+且Vy+3t)+—(t+3t)+f=0xy xy VZ xz xz x++3q)+—C+3t)+f=0y y Vz yz yz y(2.2.19)以及:—(t+3t)+—C+3t)+—(q+3o)+f=0VxXXZ Vyyz yz VzZzzF+3F=lb+3o)+mTxx xx xyF+3F=lT+3t +m&yyxyxy(F+3F=lk+3t)+miZZ XZ, XZ,+3t)+n?+3t)xy xzxz+3o4nt+3t)yy)yz yz.+3tAnG+3o)(2.2.20)由彈性體平衡方程（2.2.2）和面力表達關系式（2.2.4），可以得到:以及:V3o 13TV3tX+VxV3T ^^+VxV3t xr+Vx *y+VyV3o +VyV3T y^+Vy y^=0VzV3。 nZ=0Vz(2.2.21)193F-13。+m3r+n3rxxxyxz3F-13r+m3。+n3ryxyyyz3F-13r+m3r+n3。zxzyzz<（（2.2.22）利用高斯積分公式并對(2.2.18)進行化簡，我們可以得到:11Jc3。+￡3。+￡3。+y3T+y3T+y3T1vyyzzxyxyyzyzxzxz+川V（S3。 S3r S3r'（S3r S3。 S3r）Sy（S3r S3r S3。'——x+——+——zSy+m3r+n3r）+（3r+m3。+n3rV+（3r+m3r+n3。rx xyxz xy+JJ Cu-u）5F+0-VIF+W-WIF-dsyzyzru（2.2.23）注意到方程（2.2.21）、（2.2.22），最后得到：JJJ（3。+￡3。+￡3。+y3r+y3r+y3r1vJ〉xyyzzxyxyyzyzxzxz卜3F+V3F+W3Fdss（2.2.24）ru x上式即為彈性體的虛余能原理，其與(2.2.7)表示的虛功原理形成互補形式。上式的左邊代表彈性體的總虛余能，右端代表面力的變分在實際位移上所做的功。2.2.4最小余能原理在小變形情形下，彈性力學的一般理論指出，必定存在一個正定的狀態(tài)函數(shù)*。,。，。，r,r,r0xyzxyyzxzdU*-￡d。+￡ d。+￡ d。+y dr +y dr +y dr0xxyyzz（2.2.25）xyxyyzyzxzxz這里，U*稱為余能函數(shù)或余能密度。狀態(tài)函數(shù)U*是單值的函數(shù)，因而dU*是全微分，有x S。xSU-0-,S。ySU——0-S。zxySrxySryzxz.SU——0-Srxz（2.2.26）從而虛余功原理變?yōu)镴JJ包:3。+吆VISa xS。 yy+——S。+一zSrxyJJJ包:3。+吆VISa xS。 yy+——S。+一zSrxy Sr yzxy(5其中，SU+一Srru x3U0dv-JJruxxzxzJ（2.2.27）（U*+丫*）=0（2.2.28）U*-11JVU0dV為彈性體的余能，V*-JJ4f+VF+WF1s為外力余能。如果定ru xy20義n*（u*+v*）為系統(tǒng)的總余能，故有8n*=0。式（2.2.28）表示余能的極值原理，事實上，這時余能泛函）為極小值，故得最小余能原理:在滿足平衡方程和應力邊界條件的所有各組應力分量的函數(shù)中，真實的一組應力分量應使系統(tǒng)的余能（泛函）取極小值。此處限于篇幅，我們不再給出其證明過程。2.2.5連續(xù)介質(zhì)的哈密頓原理前面幾節(jié)中，我們僅介紹了彈性體系的靜力平衡問題及其原理，即體系在平衡時所取的一真實狀態(tài)，以區(qū)別與任何其他可能的一切狀態(tài)。而在工程問題中，還會涉及到考慮時間變量的動力學問題，不同時刻對應于不同的狀態(tài)。就數(shù)學本質(zhì)而言，靜力問題和動力問題沒有原則區(qū)別，只是僅僅增加了自變量的個數(shù)（即在空間坐標自變量的基礎上增加了時間變量），但其物理意義的差別是明顯的，即從靜力平衡過渡到了動力學問題。當牛頓建立了以三大定律及萬有引力定律為基礎的力學理論后，無數(shù)的自然現(xiàn)象都得到了定量的說明，事情似乎很完善了。后來拉格朗日在18世紀提出了一個變分原理，從這個變分原理出發(fā)，能夠十分方便地解決許多力學問題，并且由此還可以推導出力學中的很多定律。他還創(chuàng)立了拉格朗日運動方程，其比牛頓的運動方程適用的范圍更廣，而且用起來更為便捷。此后，哈密頓（Hamilton）發(fā)展了拉格朗日的理論，于1834年提出了有名的哈密頓原理。本節(jié)中，我們將引入對應于動力學問題的哈密頓原理。首先介紹離散質(zhì)點系統(tǒng)的哈密頓原理，然后將其推廣得出彈性連續(xù)體的形式。設具有N個質(zhì)點的系統(tǒng)相對于慣性參考系的位移由矢量r「r2,…，N給出，根據(jù)質(zhì)點系的達朗伯（D’Alembert）原理有：（2.2.29）工（F—mrK（2.2.29）iiiii=1上式中mi為第i個質(zhì)點的質(zhì)量，F(xiàn)i為作用于第i個質(zhì)點上的力。我們來考察動能的變分3T=33T=3（1g）—N^mr2

2ii

'i=1 ）g..."mr3r=iiii=1dfy. ）gy一乙mr3r一乙mr3rA/ iii iii?i=1 /i=1（2.2.30）結合（2.2.29），有:3t+3W=3t+3W=（2.2.31）——mr3rTOC\o"1-5"\h\z\o"CurrentDocument"Af iii?i=1 /這里，T=1gmr2表示質(zhì)點系的動能，W=gFr為外力所作的功。對（2.2.31）在任意2ii iii=1 i=121兩個時刻11和12間關于時間進行積分，得到:)mr5riii)(2.2.32)因為系統(tǒng)在時刻11和t2的位置狀態(tài)可以認為是給定的，便有5r|t=t=5r因而（2.2.32）進一步成為：JJ2（5T+5Wht=0ti如果所有的外力均為有勢力，即F=-aV/ar.（i=1，2,…，N），則6W=-6V，(2.2.33)故上式還可寫成：5Jt2(T—VLt=0ti(2.2.34)記L=T-V，稱為拉格朗日函數(shù)，對于保守系統(tǒng)，積分運算和變分運算是可以交換的，即有5J'2Ldt=0。ti所以，哈密頓原理可以敘述如下：對于有勢力作用下的完整質(zhì)點系而言，在由時刻11狀態(tài)到時刻12狀態(tài)的所有可能的運動中，實際實現(xiàn)的運動使得積分表示的泛函:J=Jt2Ldtti取極值。這里J有時也被稱作哈密頓作用量。(2.2.35)引入廣義坐標｛q1，%，…，qN｝，則L=L（q,q，…,q;q,q，…,q;t），進而哈密頓原理

12 N 1 2N1 2N（2.2.34）對應的（保守系統(tǒng)）拉格朗日方稱為d(al)al~-dtla(i)aq'n， n(n=1，2,…，N)(2.2.36)雖然，對于所假定的系統(tǒng)哈密頓原理和拉格朗日方程是等價的，但前者可適用于具有無窮多自由度的系統(tǒng)，因而從這個意義上講，哈密頓原理的適用性更廣泛。接下來，我們將哈密頓原理從離散質(zhì)點系推廣到彈性連續(xù)系統(tǒng)。將慣性力加入到彈性連續(xù)系統(tǒng)的虛功原理中，即:JJJC5s+o5s+o5s+t5y+t5y+t5yVxxyyzzxyxyyzyzxzxzfff(d2u d2v d2w^]=-JJJp 5u+ 5V+ 5wdv+V[d12 d12 d12 }+JJF5u+F5v+F5w1srtxyzJJJ(f5u+f5v+f5wlvVxyz(2.2.37)22如前所述，如果外力勢記為匕應變能記為U,即:一5V=111f5u+f5v+f5w%v+11F5u+F5v+F5w%s(2.2.38ab)Vxyz rtxy(2.2.38ab)5U=JJJb5s+o5s+o5s+t5y+t5y+t5ydvvVxxyyzzxyxyyzyzxzxz進而有:_111p_111pV(蘇5u+瓦5v+標5wJdv_5U+□=0(2.2.39)對上式在從時刻11和12間關于時間進行積分，并采用分部積分等算法，可以得到:1t1t2TLt1v「I2+”2+W1v為彈性連續(xù)體的動能。(2.2.40)令n=U+V為系統(tǒng)的總勢能，以及L=T-n為拉朗日函數(shù)，（2.2.40）進一步可表示成:t2Lt2Ldt]二0(2.2.41)此時，哈密頓原理敘述為：彈性連續(xù)體從時刻11狀態(tài)到時刻12狀態(tài)的所有可能的運動（包括彈性體的形變）中，實際實現(xiàn)的運動使拉格朗日函數(shù)在這段時間內(nèi)對時間的積分取極值。這里以梁的振動問題為例來說明哈密頓原理的具體運用。不難寫出梁的動能和應變能分別為:，T=11Lp2，T=11Lp202dx)2dx(2.2.42)其中，L表示梁的長度為梁的單位長度的質(zhì)量，w（x，t）為梁的撓度。根據(jù)哈密頓原理:1t21t2T-U■dLt1 」r32w丫(3x2Jdxdt(2.2.43)3w35w 32w325wdxdt=0dxdt=03131 3x23x2因為有6w（x，11）=6w（x，12）=0，于是:t1二1dxt1二1dxdt=1121Lt1 0I3r3w,

p1p——5w"5w[dxdt312 I啜5jdx_t1(2.2.44)對于（2.2.43）中的第二項關于應變能的變分，進一步化簡后得到:23r S2wS25w'EI V Sx2Sx2,dxdtSrS2wS5w—— z-SxVSx2SxS2rS2w—■—Sx2VSx2(2.2.45)5wdxdt+EIjt2」 11d.x2d.xSxVSx2Jx-Ldt最后得到:-j12jL11 0S2wstS4w、——+EI一St25wdxdtjlrSw丫2+Ip—5wIdx-EI0VStJ11根據(jù)變分法預備定理，S2wS5wSrS2w(2.2.46)得梁的自由振動方程為:S2w S4w八p +EI 0St2 Sx4(2.2.47)以及相應的自然邊界條件:(1)w=0或-=0

St(當t=11及t=12)(2.2.48a)(2)（3）S2wSx2S3wSx3-0或w=0（當x=0及x=L）(當x=0及x=L)(2.2.48b)(2.2.48c)其中（1）相當于梁的初始速度給定為零或者位移為零；（2）相當于給定梁在兩端部的彎矩為零或者轉(zhuǎn)角為零；而（3）相當于給定梁在兩端部的剪力為零或者位移為零。242.3變分法的近似解法數(shù)學物理中的變分原理建立了各種類型的微分方程邊值問題與泛函取駐值的等價關系。變分問題的古典解法是通過解歐拉方程來解變分問題，然而，因為求解微分方程往往并不容易，所以這個方法并不能達到預期的結果，這就要求我們必須直接尋求針對于變分問題的新方法，即直接方法。變分學的直接方法是指不通過解歐拉方程而直接近似地求解變分問題的方法。這種方法最先大量用于求解彈性力學問題，隨著電子計算機的廣泛使用和計算方法的發(fā)展，現(xiàn)今變分學的直接方法已有多種，它們的應用范圍也越來越廣。變分法的近似解法一立茲法及其應用立茲法是變分問題直接解法中最重要的一種，其基本思想是用選定的函數(shù)序列的有限線性組合逼近變分問題的極值曲線?，F(xiàn)用一個簡單的變分問題：(2.3.1)(2.3.2)JL]=廠1F(x,y,y')dx(2.3.1)(2.3.2)x/ 八0/ 八y(x)=0,y(x)=0o0 1來說明Ritz法的解題步驟。在此邊條是兩端固定的特殊情況，不失一般性，當邊條是非齊次的，即:x一xx一x而二者不同時為零時，作函數(shù)代換：z=y(x) y--1——yx一x1x一x01 0 1 0便有：z(x0)=%(q)=0解題步驟：(1)取定一相對完備函數(shù)列u(x),u(x), ,u(x),并使其中每一個都滿足邊條(2.3.2)。該序列選取對下一步計算復雜程度有很大影響。(2)將線性組合：yk(x)=,]6)+C2u2(x)+ +Cu^(x) (2.3.3)視作(2.3.1)的近似解，將(2.3.3)式代入(2.3.1)-得關于C,C,,C的函數(shù)，1 2 kJ[c,C,,c]=Jx1Fx,^kCvt^(x),￡cu(x)dx (2.3.4)% i=1 i=1(3)求(2.3.4)的極值，由方程組:25

aj=0(i=1,2,k)，解得：C,C,,Cac 12ki??????代入(2.3.3)，便得到問題的近似解。憂k(X)的選取:(1)對沒有約束條件的問題，可取函數(shù)系:xi(i=1,2,)或者:1,cos兀x,sin兀x,cos2兀(1)對沒有約束條件的問題，可取函數(shù)系:xi(i=1,2,)或者:1,cos兀x,sin兀x,cos2兀x,sin2兀x, ,cosk兀x,sink兀x(k=1,2,)(2)若要求邊條：y(0)=y(D=。，可取函數(shù)系為:3-X(1-X),n-1,2 ，3,3X,3X2, ,3Xn若要求邊條：y(0)=y(1)=y'(0)=y"(1)=0，可取函數(shù)系為:3,3x,3x2, ,3xn，3=x2(1-x)2,n-1,2(4)對于積分形式的泛函，即PDE的邊值問題，常取多項式為函數(shù)系:a)當「為矩形時，可取函數(shù)系為:3,3X,3y,3X2,3町,3y2,3%3,3%2y,3町2,3戶,342-X2)Q2-y2)b)當r為圓X2+y2=a2時，可取函數(shù)系為:3,3x,3y,3x2,3xy,3y2,3x3,3x2y,3xy2,3y3,2.3.2變分法的近似解法一一伽遼金法(Galerkin)及其應用在實際問題中，不是所有邊值問題都存在相應的泛函，伽遼金法是更廣泛一類微分方程邊值問題的近似方法。(繞過找不到合適泛函的困難)伽遼金法精度較高，計算量不大，應用較廣泛。解題步驟：(1)取定一相對完備函數(shù)列U](x),u2(x),,uk(x),并使其中每一個都滿足邊條。(2)將線性組合：yk(x)=C1u1(X)+C2u2(X)++Cyk(x)作為微分方程的近似解，其中C1,C2,,Ck待定，顯然它滿足邊界條件。 …26(3)由』Q—f)j(Qdx=0,(i=1,2k)解出Cj代回⑶便得近似解。V立茲法和伽遼金法的局限性： …(1)對于邊界形狀比較復雜(如多邊形)的區(qū)域，想找到合適的滿足邊界條件的相對完備的函數(shù)系是困難的。即使勉強湊成，也需要相當高次的多項式或某些函數(shù)的特殊結構，需較高技巧性。(2)被積函數(shù)一般都是高次，計算量大，若積分區(qū)域稍許復雜些，則難以計算。注意：伽遼金法只適用于齊次邊界條件，對于非齊次邊界條件，可通過適當?shù)淖兞看鷵Q化為齊次邊界條件。例2.3.1：用Ritz法求變分問題：J[y]=J1(y'2-y2+4xy)Zx< 0〔y(0)=0,y(1)=0的近似解。解：選相對完備的函數(shù)序列："x)=(1-x)xk(k=1,2,).其中每個函數(shù)均滿足邊條；作線性組合：y(x)工CU(x)?！筴 iii=1(1)先選取y1=c1x(1―x)代入泛函J[y]得：4]-4]」(2(12%)2C2X2(1X)2+CX2(1Jxjyy」=JC]」=Jcj(1—2x)2—cjx2(1—x)2+qx2(1—x)dx令dJc2=0dC]得：J11c(1—2x)2—2[cx(1—x)—2x]x(1—x)d^x=00 1 1由此：—c=——，得：c=——51 3 1 9得一次近似解：y1(x)=-5x(1一x)(2)再取y2(x)=qx(1—x)+c2x2(1—x).由類似上述運算得到：27

/、…、/71 7、y#)=-2x(1-x)(礪+41%).現(xiàn)討論邛x),y2(%)的近似程度：對該變分問題，其對應的Euler方程為:F=y，2-y2+4xy即d(dF)八=0aydxl辦')￡(I=-2y〃.af a￡(I=-2y〃.——=-2y+4x,——=2y，,ay ay')ny"+y+2x=0其精確解為：y(x)=2x-至”.sinl, 1234 . , ,在x=1,2,3,4四點比較值如表:5555xyy1y20.2-0.0722-0.0888-0.07240.4-0.1256-0.1333-0.12510.6-0.1420-0.13330.14150.8-0.1050-0.0888-0.1053由表可見，y1（x）相當粗糙，但y2（x）的誤差已經(jīng)很小。例2.3.2：用Ritz法求解問題：28

28對完備系可取作：Q=(a2-%2)(b2-y2),%=q%2,Q=Qy2,一次近似解為：U1(x,y)=C1(a2-x2)(b2-y2)代入泛函得：J[c]=JaJb4C2[x2(b2-y2)2+y2(a2-x2)2]-c(a2-x2)(b2-y2)\xdy1 -a-b1 1積分后，令業(yè)=0，得:

d%5c= 14(a2+b2)得：u1(x，y)=5(a2-x2)(b2-y2)4(a2+b2)例233：用伽遼金法求邊值問題Jx2y"+xy'+(x2-1)y=01y(1)=1,y(2)=2解：因邊界條件非齊次，所以先作變量代換，把邊界條件化為齊次。令y=z+x，則上述邊值問題化為：x2z"+xz'+(x2-1)z+x3=0取Q(x)=(x-1)(2-x)，一次近似解為：Z]=4]①](x)=a1(x-1)(2-x)=aJ-2+3x-x2)z；=(3-2x)a1,z；=-2a1x2(—2aj+x(3—2x)a1+(x2—1)aj—2+3x—x2)+x3=(2-5x2+3x3-x4)a1+x3J2[(2-5x2+3x3-x4)a1+x3](x-1)(2-x)dx=0w311 3八252積分得 a+—=0,a= 42015 1311252Z1—(x-1)(2-xZ1所以原邊值問題的一次近似解為:252y=—(x-1)(2-x)+x例2.3.4：用Galerkin方法求Poisson方程邊值問題:29

S2u S2uAu=——+——=-2的近似解。解：相對完備系可取作w,wx2,wy2,wx4,wx2y2,wy4,其中：w=(a2一x2)(b2一y2)一次近似解為:U(x,y)=c(a2-x2)(b2-y2)[-2ci(a2一x2+b2一y2)](a2一x2)(b2一y2)dxdy=0積分后得代數(shù)方程:32—a332—a3b3=09-45-a3b3(a2+b2)q+4(a2+b2)因此，一階近似解為:4(因此，一階近似解為:4(a2+b2)5(a2—x2)(b2—y2)u= 30泛函的極大值和極小值問題如果函數(shù)y(x)在x=x0附近的任意點上的值都不大(小)于y(x0),也即dy=y(x)-y(x0)<0(>0)時，則稱函數(shù)y(x)在x=x0上達到極大(極小)，而且在x=x0上，有。 ° °dy=0 (2-1)對于泛函n[y(x)],也有類似的定義。如果泛函n[y(x)]在任何一條與y=y0(x)接近的曲線上的值不大(或不小)于n[y0(x)]，也就是，如果6□=R[y(x)]-R[y0(x)]<0(或>0)時，則稱泛函n[y(x)]在曲線y=y(x)上達到極大值(或極小值)，而且在y=y(x)00上，有6口二0 (2-2)在這里，對于泛函的極值概念有進一步說明的必要，凡說到泛函的極大(或極小)值，主要是說泛函的相對的極大(或極小)值，也就是說，從互相接近的許多曲線來研究一個最大(或最小)的泛函值，但是曲線的接近有不同的接近度。因此，在泛函的極大極小的定義里，還應說明這些曲線有幾階的接近度。如同一般函數(shù)極大(極小)討論一樣，如果泛函在y=y0(x)曲線上有強極大(極小)值，不僅對于那些既是函數(shù)接近而且導數(shù)也接近的y(x)而言是極大(極小)值，而且對于那些只是函數(shù)接近但導數(shù)不接近的y(x)而言，也是極大(極小)值，所以泛函在y=y0(x)曲線上是強極大(極小)值時，也必在y=y0(x)上是弱極大(極小)值。反之，則不然，即泛函在y=y0(x)曲線上有弱極大(極小)值時，不一定是強極大(極小)值，因為有可能對于那些只是函數(shù)接近但導數(shù)不接近的y(x)而言，有一個比函數(shù)與導數(shù)都接近的y(x)所求的極大(極小)更大(小)的極大(極小)值存在。所以弱極大(極小)，不能滿足強極大(極小)的要求。這一概念可以推廣到包含多個函數(shù)的泛函中去。求解泛函極值的歐拉方程變分法的早期工作是如何將泛函駐值問題轉(zhuǎn)化為微分方程問題。當把泛函的駐值問題轉(zhuǎn)化為微分方程時，第一步工作就結束了，下一步是如何求解這一微分方程。這種求解方法在實際應用上碰到很大的困難。自從里茲提出直接求泛函極值的近似法(里茲法)以后，人們才認識到直接從泛函極值出發(fā)，而避免從微分方程式出發(fā)更為有效與方便，這樣的處理方法可以充分利用電子計算機的作用。于是人們研究的目標有所轉(zhuǎn)移，即把原來從泛函駐值問題化為微分方程問題，轉(zhuǎn)變?yōu)榘盐⒎址匠虇栴}轉(zhuǎn)變?yōu)槎x一個泛函，而成為泛函求駐值的問題。對于前一種問題由歐拉、拉格朗日等已建立了一套比較成熟、比較系統(tǒng)的方法，而對于后一類問題，雖然正在大力進行工作，但尚不成熟。目前用的多的方法，還是根據(jù)微分方程物理和工程背景，采取嘗試和核對的方法，即先試猜一個泛函的極值和駐值問題，然后再核對一下，看它是否與原來的微分方程問題等價。這種方法在以后的變分原理中將經(jīng)常用到?，F(xiàn)在研究最簡單泛函(2-3)式的極值問題所得到的歐拉方程，其中能確定泛函極值曲線y=y(x)的邊界是固定不變的，而且有y(xj=yJy(x/=y2，函數(shù)F(x,y,y')將認為是三階可微的。 1口二廠2F[x,y(x),y'(x)]dx (2-3)x1首先讓我們用拉格朗日法來求泛函的變分n[y+s6y]=尸2F[x,y+sSy,y'+s6y']dxxi于是有31

a_ _一一a_ ..一n[y+￡8y]-a_ _一一a_ ..一n[y+￡8y]-jx2{一F[x,y+￡8y,y'+￡8y^]8y+as %aya—F[x,y+s8y,y+s8y]8y}dxaya aF aF8n--n[y+￡8y]l -卜?[丁8y+-^8y']dxas ￡-o %ay ay(2-4)aF a其中k--F(x,y,y)，ay ay而且「2竺8y'dx=「2{—［逞8y］--d(竺)8y}dx

%ay， x1dxay，」dxay，對于固定邊界條件，因為有8y(x)-8y(x)-0，所以af.——8ydx=ay，y—「2_d_(aF)8ydx

%dx%y〃y(2-5)將(2-5)式代入(2-4)式，得到變分極值條件8n=fx2［吧—A(空)］8ydx=0%1ay dxay，)根據(jù)變分法的基本預備定理，求得本題的歐拉方程為aF_d(aF)-0 -ay dxay1這里必須指出，上式中的第二項是對x的全導數(shù)，不是偏導數(shù)，且F-F(x,y,y')(2-6)(2-7)所以daF、 a2F a2Fdya2Fdy'—(——) + + —dxay'axay'axay'dxay'2dx-F"+F"y'+F"y〃、…， 、八J(2-8)町yy其中F" F" F"，都是F(x,y,y')對x,y,y'的二階偏導數(shù)。y'=?，y"-學，所xy yy yy dx dx2以歐拉方程(2-7)式也可以寫成F'-F"-F"y'-F"y"-0 (2-9)y xy yyyy這就是1744年歐拉所得的著名方程。該方程也被稱為歐拉-拉格朗日方程。(2-9)式是關于y(x)的一個二階微分方程，其積分常數(shù)有兩個］和c2，它的積分曲線y-y(x,c/c2)叫做極值曲線，只有在這族極值曲線上，泛函(2-3)式才能達到極值，積分常數(shù)是由極值曲線通過y(x1)-yjy(x2)-y2這兩個端點條件所決定的。把泛函的變分作為泛函增量的主部，也同樣得到歐拉方程(2-7)式及(2-8)式。求泛函增量主部的過程實質(zhì)上與求微分的過程非常相似。例如從(2-3)式，因為積分限是固定的(不變的),所以有8口-8「2F(x,y,y')dx-「28F(x,y,y')dx其8F是從y，y'增量引起的，其主部為 ‘a(chǎn)F aF8F(x,y,y)--8y+寺8y

ay ay于是得到(2-4)式，這和拉格朗日法得到的變分表達式是相同的。這里還應指出，(2-9)式這樣的歐拉方程，有下列四種特殊的情況，應該予以注意。(2)F(x,y,y')和x無關，即(2-10)F-F(y,y(2-10)32于是(2-9)式可以寫成上式可以簡化為Ff77〃 ， 77〃 〃—Fy上式可以簡化為Ff77〃 ， 77〃 〃—Fy—Fy-0y yy' y'y'(2-11)一次積分后—(F_^y，)=0dx 6y，(2-12)F-笠y，-c6/, 1(2-13)其中其中c1為積分常數(shù)。(2)F(x,y,y')和y無關，(2-14)代入(2-7)式，得(2-14)代入(2-7)式，得(2-15)積分得6F——二c6y，(2-16)其中其中c為積分常數(shù)。(3)F(x,y,y')和y'無關于是歐拉方程為(2-17)(2-18)它不是微分方程，不包含什么特定常數(shù)，一般情況，所討論的變分問題不存在，只在個別的情況下，當曲線(2-18)式通過固定端點時，才存在可能達到極值的曲線。(4)F(x,y于是歐拉方程為(2-17)(2-18)它不是微分方程，不包含什么特定常數(shù)，一般情況，所討論的變分問題不存在，只在個別的情況下，當曲線(2-18)式通過固定端點時，才存在可能達到極值的曲線。(4)F(x,y,y')是y'的線性函數(shù)，即FQ,y,y')=PQ,y)+QQ,y)/(2-19)于是歐拉方程為6P6Q,dQ八——+—yr——=06y 6y dx(2-20)但是dP_6Q上6Q, 十 ydy 6x 6y(2-21)所以(2-20)式可以簡化為6P6Q————=06y 6x(2-22)它也不是一個微分方程式，因為它沒有y'項，一般說來它不滿足固定端點條件，因此，變分問題根本不存在。現(xiàn)在我們將上述變分問題推廣到含有高階導數(shù)的泛函的極值問題和泛函變分得到的歐拉方程。我們研究泛函n[y(x)]二尸2F[x,y(x),y'(x),y〃(x), ,y(n)(x)]dx(2-23)x1的極值，其中泛函F被認為對于y(x)，y'(x)，y'f(x),...,y(n)(x)是n+2階可微的，并且假定，端點上有固定條件33y(n-1)(x)=y(ny(n-1)(x)=y(n-1)1 1，y(n-1)(x)=y(n-1)2 2Y (2-24)我們將假定，極值在2〃階可微af.?+ 8y(n)}dxay(n)（2-25）（2-26）（2-27）（2-28）（2-29）（2-30）梁的位能等（2-31）其中X為梁彎曲后的曲率，它和撓度w（x）的關系為d2wdx2y(x)=y，y'(x)=y'，y"(x)=y"，…，.11 11 11y(x)=y，y'(x)=y'，y"(x)=y"，…，22 22 22端點上不僅給出函數(shù)值，而且還給出直至"-1階導數(shù)的值。曲線y=y（x）上達到。用上面相同的求泛函變分方法，我們可以證明：. af af af6n=jx2｛一Sy+一Syr+一8y〃+..x1ay ay ay其中用簡略符號8y代替8y(x)，8y(k)代替8y(k)(x)=上hy(x)1dxk積分(2-25)式中的第二項可以分部積分一次，得jx2空A(8y)dx=空8yIx2-jx2—(空)8ydxx1aydxay)x： “1dx'ay77將積分(2-25)式中第三項分部積分兩次，得 1jx2空上(8y)dx二空8yrIx2--d-(空)8yIx2+jx2_dL(空)8ydxx1 ay dx2“ a/「j dx 8y〃 x； “1dx2”/最后一項經(jīng)過n次分部積分后，得xafdn af dafx2 (8y)dx= 8y(n-1)Ix2— ( )y(n-2)|x2+x1ay(n)dxn dy(n) x1 dxdy(n) %—F(-1)njx2 ()8ydxx1dxnBy(n)根據(jù)變分法的預備定理，(2-25)式為零時，得 1d,af、d2af、 /1、dnaf、_F'- (——)+——(——)+…十(-1)n一( )=0ydxdy' dx2ay" dxndy(n)這是y=y（x）的2n階微分方程式，一般稱之為泛函Q-23）式的歐拉-泊桑方程，而它的積分曲線就是所討論變分問題的解（極值曲線）。這個方程的解通常有2n個特定常數(shù)，由2n個端點條件（2-24）式?jīng)Q定的?！纠?-1】梁在橫向載荷作用下的彎曲問題，就是含有較高階導數(shù)的泛函極值問題的一個例子。設梁的抗彎剛度為EJ，兩端固定，在橫向分布載荷q（x）作用下發(fā)生彎曲變形（或稱撓度）例x），如圖2-1所示。端點固定條件為攻（0）=攻w（0）=01攻（L）二攻'（L）=0J在梁達到平衡時，其總位能達到最小值。于梁在彎曲時所貯存的彎曲能，它等于?1 .U=jLEJz2dx02dx2dwT1+(=”Idx」這里假定撓度很小，略去高次項。（2-31）式可以寫成34

1 d2wU=JLEJ(—)2dx

02 dx載荷q(x)在變形—(x)上的位能為V=—JLq(x)—(x)dx (2-32)0于是，梁所形成的總位能口為1 d2—n=U+V=JL[EJ( )2—q(x)—(x)]dx (2-33)02 dx梁的平衡條件為—(x)使總位能達到最小值，即汨=0。于是利用變分計算，并利用固定端條件(2-30)式，得d4—Sn=JL[EJ ——q(x)]S—(x)dx=0 (2-34)0 dx4利用變分法的預備定理，求得梁的平衡方程為d4—EJ——q(xd4—EJ——q(x)=0dx4這就是歐拉-泊桑方程。注意到(2-34)式在靜力學中被稱為虛位移原理，滿足端點位移約束條件的虛位移。下面討論另一種形式的泛函n(0,0,0)=JJF(0,0,0)dxdy+JG(0)ds(2-35)xy的歐拉方程。函數(shù)中0(x,y)在域R內(nèi)連續(xù)，其邊界S由S和S組成Sc°b為給定的,式中中b303y。bc(在S上)

b其中現(xiàn)在對(2-35)泛函取一次變分，得到8n=JJ[8n=JJ[R3Fc不3