彈塑性力學第十章

2023/2/61第十章彈性力學的能量原理§10-1幾個基本概念和術語§10-2虛功方程§10-3功的互等定理

§10-4虛位移原理和最小勢能原理§10-5虛應力原理和最小余能原理§10-6基于能量原理的近似解法2023/2/62第十章彈性力學的能量原理

彈性力學的解法之一為彈性力學邊值問題求解體系——靜力法。在前面各章中就圍繞平面問題、扭轉(zhuǎn)問題和空間軸對稱問題進行了具體分析和研究。2023/2/63第十章彈性力學的能量原理

彈性力學問題的解法還有另一種解法：以能量形來建立彈性力學求解方程——能量法（從數(shù)學意義上說也可認為變分法）。本章主要介紹幾個基本能量原理以及基于能量原理的近似解法。在介紹能量原理以前，先介紹幾個基本概念和術語。2023/2/64§10-1幾個基本概念和術語1.1應變能U和應變余能Uc：dijijij應變能

U在第四章中已定義過：應變能密度2023/2/65——彈性關系

如果將幾何關系引入應變能,U、W

為位移的函數(shù)。應變余能（類似應變能）定義§10-1幾個基本概念和術語2023/2/66應變余能密度——單位體積的應變余能Wc與積分路徑無關，只與終止狀態(tài)和初始狀態(tài)有關。

Wc=ijij

為全微分§10-1幾個基本概念和術語dijijijdij2023/2/67

——逆彈性關系且

W+Wc=ijij

§10-1幾個基本概念和術語dijijijdij2023/2/68

但§10-1幾個基本概念和術語材料為線彈性時2023/2/69各向同性線性材料的應力應變關系§10-1幾個基本概念和術語將幾何關系引入上式U=U(ui)應變能是位移的函數(shù)2023/2/610

代入Uc表達式各向同性線性材料的應力應變關系§10-1幾個基本概念和術語2023/2/611§10-1幾個基本概念和術語應變能、應變余能的計算舉例xolP圖示等截面桿，承受軸向荷載P作用。桿截面面積為

A,材料應力—應變關系分別為（1）=E，(2)=E1/2.試計算外力功T、應變能U和應變余能Uc。解：（1）=E2023/2/612§10-1幾個基本概念和術語xolP

T=U=Uc=Pl/2l=Pl/(EA)P=N=lEA/l,U=l2EA/(2l),Uc=P2l/(2EA),（2）=E1/2T=U=2023/2/613§10-1幾個基本概念和術語xolP

T=U2023/2/614§10-1幾個基本概念和術語xolPUc=Pl–U=Pl-U2023/2/615§10-1幾個基本概念和術語作業(yè)：圖示結(jié)構(gòu)各桿等截面桿，截面面積為A,結(jié)點C承受荷載P作用,材料應力—應變關系分別為（1）=E，(2)=E1/2。試計算結(jié)構(gòu)的應變能U和應變余能Uc。lPCBAx

ylC’2023/2/6161.2可能位移ui(k)和可能應變ij(k)：可能應變ij(k)：由ui(k)通過幾何方程導出的

可能位移ui(k)：在V內(nèi)連續(xù)且可微，在

su上滿足：§10-1幾個基本概念和術語2023/2/6171.3可能應力ij(k)：ij,j(k)+fi=0

(a)在s上滿足

(b)滿足式(a)、（b)——滿足靜力方程可能應力

ij(k)：在V內(nèi)滿足

§10-1幾個基本概念和術語2023/2/618

兩種可能位移ui(k1)和ui(k2)之差稱為虛位移ui，而由兩種可能位移狀態(tài)對應的可能應變

ij(k1)

、ij(k2)之差稱為虛應變ij

。1.4虛位移ui和虛應變ij

：ui=0

在su上齊次位移邊界條件。ij=(ui,j+uj,i)/2

在V內(nèi)§10-1幾個基本概念和術語2023/2/6191.5虛應力ij

：在s

上:njij=0；ij=ij(k1)-ij(k2)

§10-1幾個基本概念和術語在V內(nèi)：ij,j=0滿足齊次靜力方程。2023/2/620§10-2虛功方程2.1虛功方程Su

S

在給定體力、面力和約束情況下，如果找到兩種狀態(tài):第一種狀態(tài)：

在給定的體力fi和面力，已知（找到）可能應力狀態(tài)ij(k1)在V內(nèi)：ij(k1)+fi=0

;在s=s:

2023/2/621第二種狀態(tài)：

彈性體處于可能變形狀態(tài)ui(k2)

、ij(k2)

則第一種狀態(tài)外力在第二種狀態(tài)可能位移作的外力虛功等于第一種狀態(tài)可能應力在第二種狀態(tài)可能應變上作的虛變形功。

——虛功原理§10-2虛功方程在s=su:

2023/2/6222.2虛功方程的證明：§10-2虛功方程2023/2/623§10-2虛功方程2023/2/624代入虛功方程左端，得并注意虛功方程未涉及本構(gòu)關系，所有在各種材料性質(zhì)虛功方程成立。則We=Wi§10-2虛功方程2023/2/625虛功方程雖然對兩種不相干的可能狀態(tài)成立，但一般應用是一種為真實狀態(tài)，另一種為虛設可能狀態(tài)（虛設狀態(tài)）。qP=1§10-2虛功方程2023/2/626§10-3功的互等定理將虛功方程用于線彈性體可導出功的互等定理。同一彈性體處于兩種真實狀態(tài)。第一種狀態(tài)：滿足所有方程。

第二種狀態(tài)：

滿足所有方程。2023/2/627§10-3功的互等定理根據(jù)虛功方程第一種狀態(tài)的外力在第二種狀態(tài)的相應彈性位移上做功第二種狀態(tài)外力在第一種狀態(tài)的相應彈性位移上做功2023/2/628§10-3功的互等定理對于線彈性體本構(gòu)關系

W12=W212023/2/629§10-3功的互等定理功的互等定理優(yōu)點：可以避免求解物體內(nèi)的應力、應變和位移場的復雜過程，而直接從整體變形的角度來處理問題。PPb第一種狀態(tài)的外力在第二種狀態(tài)的相應彈性位移上所做的功等于第二種狀態(tài)外力在第一種狀態(tài)的相應彈性位移上所做的功。第一狀態(tài)：一對力P作用在直桿的垂直方向，局部效應，在桿兩端點伸長？2023/2/630§10-3功的互等定理

第二狀態(tài)：讓一對力Q作用同一桿兩端點，很易求得一對力Q引起桿橫向縮短。對兩種狀態(tài)應用功的互等定理P=Q

Q第二狀態(tài)引起的

易求:QQx2023/2/631§10-3功的互等定理QQx2023/2/632

§10-4虛位移原理和最小勢能原理4.1虛位移原理

運用虛功原理，但一種狀態(tài)為與真實外力平衡的狀態(tài)，ij、fi、

、;而第二狀態(tài)為可能變形狀態(tài)，為真實狀態(tài)位移的變分:

ui=0

在su上

ui、ij=(ui,j+uj,i)/2

在V內(nèi)虛設狀態(tài)2023/2/633§10-4虛位移原理和最小勢能原理根據(jù)虛功方程，真實的外力與應力狀態(tài)在虛設的齊次可能位移上做功彈性體應力與外力處于平衡狀態(tài)，對于任意虛設的齊次微小位移及應變，則外力在虛位移上做的虛功等于應力在虛應變上做的虛功

——虛位移方程。2023/2/634§10-4虛位移原理和最小勢能原理代入原虛位移方程將虛位移方程重新改寫2023/2/635§10-4虛位移原理和最小勢能原理代入原虛位移方程虛位移方程為平衡方程和力的邊界條件的積分形式。虛位移原理舉例2023/2/636§10-4虛位移原理和最小勢能原理虛位移原理舉例

圖示受均布荷載q作用的等跨連續(xù)梁，EI為常數(shù)，中間支座為彈性支座。試用虛位移原理寫出梁的撓曲線方程和邊界條件。lACBqlx

z解：圖示連續(xù)梁在荷載作用下，產(chǎn)生撓曲線w(x)、內(nèi)力和支座反力。

設連續(xù)梁有虛位移

w

，C點虛位移

wC

。2023/2/637§10-4虛位移原理和最小勢能原理lACBqlx

z虛位移方程利用對稱性

在計算薄梁的內(nèi)力虛功時，只考慮梁的正應力x作的虛功。2023/2/638§10-4虛位移原理和最小勢能原理lACBqlx

z2023/2/639§10-4虛位移原理和最小勢能原理lACBqlx

z分部積分，得再積分，得2023/2/640§10-4虛位移原理和最小勢能原理lACBqlx

z虛位移方程為2023/2/641§10-4虛位移原理和最小勢能原理lACBqlx

z得平衡微分方程和力的邊界條件（由虛位移方程得到）2023/2/642§10-4虛位移原理和最小勢能原理lACBqlx

z得平衡微分方程和力的邊界條件（由虛位移方程得到）而事先要求滿足的位移邊界條件2023/2/643§10-4虛位移原理和最小勢能原理4.2最小勢能原理

1．彈性體的總勢能

(k)=U(k)(ui(k))+V(k)(ui(k))=(k)(ui(k))2023/2/644§10-4虛位移原理和最小勢能原理

(k)=U(k)(ui(k))+V(k)(ui(k))=(k)(ui(k))（1）和給定;（2）已將幾何關系引入

ij=(ui,j+uj,i)/2;（3）ui(k)為可能位移：在su上

;2023/2/645§10-4虛位移原理和最小勢能原理2．由(k)

中尋求真實位移ui(k)為可能位移,有無窮多。因此，與其對應的勢能

(k)也有無窮多。要從

(k)

中找真實位移：

（1）=0（2）引入本構(gòu)關系

真實位移應滿足的方程。2023/2/646§10-4虛位移原理和最小勢能原理取=0,得引入本構(gòu)關系

——虛位移方程

2023/2/647§10-4虛位移原理和最小勢能原理或ui(k)為可能位移，同時滿足本構(gòu)方程。而=0，表明由ui(k)導出ij(k)滿足靜力方程，所以由=0

即為真解應滿足的控制方程。

2023/2/648§10-4虛位移原理和最小勢能原理最小勢能原理的表述：

在位移滿足幾何方程和位移邊界條件的前提下，如果由位移導出的相應應力還滿足平衡微分方程和力的邊界條件，則該位移必使勢能

為駐值（極值）。如果可能位移使

的變分

=0，則該位移相應應力必滿足靜力方程。=0等價與靜力方程。2023/2/649§10-4虛位移原理和最小勢能原理作業(yè)：圖示梁受荷載作用，試利用最小勢能原理導出梁的平衡微分方程和力的邊界條件。

y

qEI

x

l

M2023/2/650§10-4虛位移原理和最小勢能原理最小勢能原理舉例

（1）已知圖示桁架各桿EA相同，材料的彈性關系為

=E，試用勢能原理求各桿內(nèi)力。lPCBAx

ylC’解：計算圖示桁架的總勢能

=U

+V

=

(uc、vc)應變能2023/2/651§10-4虛位移原理和最小勢能原理lPCBAx

ylC’則，應變能為2023/2/652§10-4虛位移原理和最小勢能原理lPCBAx

ylC’荷載勢能為

=U

+V

=由總勢能

的變分

=0

，得2023/2/653§10-4虛位移原理和最小勢能原理lPCBAx

ylC’解得2023/2/654§10-4虛位移原理和最小勢能原理lPCBAx

ylC’各桿內(nèi)力為

2023/2/655§10-4虛位移原理和最小勢能原理(2)已知圖示ab薄板（厚度t=1）無體積力作用，試用最小勢能原理求位移。xq2yq1ab解：圖示薄板受雙向壓縮作用，可猜應力和應變?yōu)槌Ａ俊?/p>

所以，位移為x,y的線性式.采用最小勢能原理求位移2023/2/656§10-4虛位移原理和最小勢能原理xq2yq1ab位移邊界條件

x=0:u=0,y=0:v=0考慮位移邊界條件取位移：

u=A1x,v=B1y薄板的總勢能

=U

+V

2023/2/657§10-4虛位移原理和最小勢能原理xq2yq1ab平面應力問題：代入應變能表達式，得2023/2/658§10-4虛位移原理和最小勢能原理

u=A1x,v=B1y代入上式，得2023/2/659§10-4虛位移原理和最小勢能原理荷載勢能V為xq2yq1ab

=2023/2/660§10-4虛位移原理和最小勢能原理由總勢能的變分

=0

，得解得2023/2/661§10-4虛位移原理和最小勢能原理代回u=A1x,v=B1y

，得求得應力為x=-q1,y=-q2,

xy=02023/2/662§10-5虛應力原理和最小余能原理5.1虛應力原理1．虛應力方程運用虛功原理，但第一種狀態(tài)為真實變形狀態(tài)，ui和ij

、fi、、第二狀態(tài)為自平衡狀態(tài)的可能應力（或真實應力的變分）ij；滿足：

ij,j=0

在V內(nèi)

njij=0

在s上（在s

上無面力）

2023/2/663§10-5虛應力原理和最小余能原理ij在

su上產(chǎn)生

Xi（在su上有反力）

Su

S

Xi

根據(jù)虛功方程2．虛應力方程表達

彈性體的應變與位移處于相容狀態(tài)，對于任意虛設的齊次容許應力ij及位移邊界上的虛反力Xi，虛應力在應變上做的虛功等于虛反力在給定位移上做的虛功。2023/2/664

Rc

RB

RA§10-5虛應力原理和最小余能原理3．虛應力原理舉例

圖示受均布荷載q作用的等跨連續(xù)梁，EI為常數(shù)。試用虛應力原理求梁的反力和彎矩。lABCqlx

z解：連續(xù)梁在荷載q作用下，產(chǎn)生撓曲線、內(nèi)力和支座反力。根據(jù)平衡關系，知2023/2/665

RA

RB

RC

Rc

RB

RA§10-5虛應力原理和最小余能原理梁的彎矩為設連續(xù)梁有虛反力和虛彎矩：lABCqlx

z應用虛應力方程lCAlx

zB2023/2/666

RA

RB

RC

Rc

RB

RA§10-5虛應力原理和最小余能原理lABCqlx

zlCAlx

zB在

su上無位移2023/2/667

Rc

RB

RA§10-5虛應力原理和最小余能原理lABCqlx

z得解得2023/2/668§10-5虛應力原理和最小余能原理5.2最小余能原理已知變形體在體力fi、面力作用及在位移邊界上有給定位移

。定義：由可能應力狀態(tài)ij(k)表示c(k)=Uc(ij(k))+Vc(Xi(k))=c(k)(ij(k))

——變形的總余能1．變形體的總余能c(k)2023/2/669§10-5虛應力原理和最小余能原理c(k)=Uc(ij(k))+Vc(Xi(k))=c(k)(ij(k))

——變形的總余能而Xi(k)=njij(k)

在su上（位移邊界上的反力）——應變余能——位移邊界的余能

由ij(k)導出的在su邊界上的反力2023/2/670§10-5虛應力原理和最小余能原理結(jié)構(gòu)總余能c(k)由可能應力ij(k)

定義。ij(k)

滿足

在V內(nèi)：ij(k)+fi=0在s=s:

（1）

c=02．由ij(k)定義的總余能中找出真實應力ij（2）引入關系式—逆彈性關系2023/2/671§10-5虛應力原理和最小余能原理c

=0，即——虛應力方程2023/2/672§10-5虛應力原理和最小余能原理由于ij

滿足自平衡的應力狀態(tài):ij,j=0

在V內(nèi),njij=0

在s上c=0

為一個虛應力方程,則ij(k)為可能變形狀態(tài)，而ij(k)已滿足靜力方程，

由ij(k)導出的ij(k)、ui(k)

滿足幾何方程及位移邊界條件。2023/2/673§10-5虛應力原理和最小余能原理因此，由c=0

表明由ij(k)中找出真實應力ij。3．最小余能原理表述：

在應力滿足平衡微分方程和應力邊界條件的前提下，如果由應力導出的相應應變還滿足相容條件，則該應力必使總余能c

為極值；或可能應力使得總余能的變分c=0,則由該應力導出相應位移必滿足相容條件。2023/2/674§10-5虛應力原理和最小余能原理EIB

x

l

qEI

y

lACRBRC最小余能原理舉例

兩跨連續(xù)梁受均布荷載q作用的彎曲問題，試用最小余能原理求解。RA解:（1）確定梁的可能應力狀態(tài)：梁的彎曲問題可選取支座反力RB為廣義可能應力梁的彎矩和其它支座反力可由RB表示。2023/2/675§10-5虛應力原理和最小余能原理EIB

x

l

qEI

y

lACRBRCRA2．確定結(jié)構(gòu)的總余能c(k)2023/2/676§10-5虛應力原理和最小余能原理c(k)=Uc(RB)+Vc(RB)3．由余能的變分c=0

確定RB

2023/2/677§10-5虛應力原理和最小余能原理即積分，解得2023/2/678§10-5虛應力原理和最小余能原理代入可得梁的最后彎矩方程2023/2/679§10-5虛應力原理和最小余能原理作業(yè)：利用虛應力原理和最小余能原理，求梁的反力和彎矩。

y

qEI

x

l2023/2/680§10-5虛應力原理和最小余能原理4．真實狀態(tài)總勢能總余能c關系：——虛功方程

所以真實狀態(tài)

=-c對于真實狀態(tài)的2023/2/681§10-6基于能量原理的近似解法

在前面幾節(jié)介紹了幾個最基本的能量原理，利用能量原理求問題的解，從理論上看是明確的步驟規(guī)范。如最小勢能原理：（分兩步）1．對于給定的外力和邊界條件尋找滿足幾何方程和位移邊界條件的ui(k)

函數(shù)序列，并確定

(k)

。2．由

=0尋求真解ui，即由ui(k)中找ui的控制方程，但由于問題求解域復雜性及約束的變化，利用能量原理求解析解也是無法實際進行的。但可由能量原理可以建立尋求問題近似解的有效途徑。2023/2/682§10-6基于能量原理的近似解法6.1基于虛位移原理的近似解法虛位移原理可以用來求問題的近似解法。

在給定的體積力、邊界力和邊界位移情況下，真實應力、應變和位移狀態(tài)（它們滿足所有方程）對于任意虛位移和虛應變滿足虛位移方程。虛位移原理近似解法的步驟2023/2/683§10-6基于能量原理的近似解法（1）選取可能位移（近似解），包含若干待定系數(shù)；（滿足位移邊界方程）（2）由可能位移求可能應變（滿足幾何方程）以及應力（應力一般不是可能應力），它們包含若干待定系數(shù)；（3）設滿足齊次邊界條件的虛位移，并導出相應的虛應變；2023/2/684§10-6基于能量原理的近似解法（4）將外力和包含若干待定系數(shù)的應力對虛位移作功等于零——虛位移方程（應力近似滿足靜力方程）；（5）由虛位移方程得到確定若干待定系數(shù)的方程，并由方程解出待定系數(shù)，從而得位移的近似解。舉一個應用虛位移原理求問題的近似解的例題2023/2/685§10-6基于能量原理的近似解法圖示簡支梁受均布荷載q作用,材料應力—應變關系分別為=E，試確定梁撓曲線的近似解

。解：（1）設梁的近似解為（包含若干待定系數(shù)）：

qEI

y

x

lv=B1x(x-l)——滿足邊界條件2023/2/686§10-6基于能量原理的近似解法（2）由梁的近似解求可能應變以及應力（應力一般不是可能應力）：（3）設滿足齊次邊界條件的虛位移，并導出對應的虛應變；

qEI

y

x

lv=B1x(x-l),=-v’’=-2B12023/2/687§10-6基于能量原理的近似解法（4）應用虛位移方程：（5）由虛位移方程解出待定系數(shù)B1，從而得位移的近似解。

qEI

y

x

l2023/2/688§10-6基于能量原理的近似解法得材料力學的精確解2023/2/689§10-6基于能量原理的近似解法作業(yè)：利用虛位移原理的近似法求梁的彎矩。xyPEIl/2l/2P2023/2/690§10-6基于能量原理的近似解法u(k)=u0+Amum,v(k)=v0+Bmvm,w(k)=w0+Cmwm,1．選取可能位移(在給定的條件下選可能位移)式中u0

、um、v0、vm、w0、wm

均為已知連續(xù)可微函數(shù)，6.2基于最小勢能原理的近似解法（Ritz法）2023/2/691§10-6基于能量原理的近似解法而um=0、vm=0、wm=0、在Su上可能位移中的

Am、Bm和Cm為待定系數(shù)。在Su上且u0

、v0、w0滿足Su

的位移邊界條件2023/2/692§10-6基于能量原理的近似解法2．結(jié)構(gòu)的總勢能

(k)及其變分

(k)

此時結(jié)構(gòu)的總勢能不是泛函了，而是Am、Bm、

Cm的函數(shù)。

(k)=

(k)(Am、Bm、Cm)=U(Am、Bm、Cm)+V(Am、Bm、Cm）2023/2/693§10-6基于能量原理的近似解法3.利用

=0求解方程

=0

由于Am、Bm、

Cm的增量

Am、Bm、Cm的任意性，則

=0

2023/2/694§10-6基于能量原理的近似解法需要求對于線彈性體的應變能U為待定系數(shù)的二次式，荷載勢能V

為待定系數(shù)的一次式。2023/2/695§10-6基于能量原理的近似解法2023/2/696§10-6基于能量原理的近似解法各向同性線彈性體有關Am、Bm和Cm

的線性方程組（3m個方程）2023/2/697§10-6基于能量原理的近似解法4．Ritz法在平面問題的應用

可能位移

u(k)=u0+Amum,v(k)=v0+Bmvm

u0

、v0滿足

Su

的位移邊界條件,

而

um=0、vm=0

在Su上。

平面應力和平面應變問題均不考慮位移分量w，而u、v為x、y的函數(shù)，體積力分量fz=0，面力分量2023/2/698§10-6基于能量原理的近似解法總勢能

的變分

=0

，得總勢能

=（Am,Bm）2023/2/699§10-6基于能量原理的近似解法平面應力問題，取薄板厚度

t=1。其應變能為:對于平面應變問題（取t=1）將上式中替換可得平面應變問題應變能U的表達式。2023/2/6100§10-6基于能量原理的近似解法Ritz法在平面問題舉例

設有一無限長的薄板，上下兩端固定，僅受豎向重力作用。求其位移解答。xybgo2023/2/6101§10-6基于能量原理的近似解法解：由于問題對y

軸對稱，所以推論：(1)選擇可能位移：設：xybgo2023/2/6102§10-6基于能量原理的近似解法滿足上下兩邊的邊界條件：

=U

+V

2023/2/6103§10-6基于能量原理的近似解法(2)計算應變能U

：（取

y軸兩側(cè)各1/2單位長度計算）xybgo1/21/22023/2/6104§10-6基于能量原理的近似解法(3)由

=0變分方程確定系數(shù)Bk：經(jīng)推導得：即2023/2/6105§10-6基于能量原理的近似解法或：2023/2/6106§10-6基于能量原理的近似解法得：這是k個獨立方程，可求出k個待定系數(shù)Bk。解得：2023/2/6107§10-6基于能量原理的近似解法

(4)位移解答令：處的豎向位移為：2023/2/6108§10-6基于能量原理的近似解法項數(shù)精確解0.12900.12420.12530.1250123值與級數(shù)項數(shù)的關系作業(yè)：設位移的近似解為

u=0,v=

B1y(y-b)，求其位移解答。2023/2/6109§10-6基于能量原理的近似解法最小勢能原理

Ritz法自變量自變函數(shù)u、v、w

由自變量Am、Bm

、Cm

定義的u、v、w自變量的

約束條件

幾何方程和位移

邊界條件

幾何方程和位移

邊界條件

總勢能

=（u、v、w）

泛函

=（Am、Bm、Cm）多元（3m元）函數(shù)變分等于零=0

積分方程——

靜力方程的積分形式

多元（線性）方程組

滿足

=0的解

解析解

近似解最小勢能原理與Ritz法的比較2023/2/6110§10-6基于能量原理的近似解法

在最小勢能原理中，由可能位移ui(k)

定義的總勢能(k)

，并由

=0尋求真解ui。也可改寫為靜力方程的積分形式等價于靜力方程。6.3伽遼金法（1915年）而

=0

本身表示為虛位移方程：2023/2/6111§10-6基于能量原理的近似解法伽遼金法步驟：1．設定滿足強約束條件的可能位移ui(k)

ui(k)

需要滿足的強約束條件：在Su上由u(k)、v(k)、w(k)

導出的應力ij(k)

在S

上

：可能位移ui(k)

滿足給定和條件，即滿足所有邊界條件——強約束條件。2023/2/6112§10-6基于能量原理的近似解法2．由結(jié)構(gòu)總勢能

的變分等于零導出求解方程

如果所設可能位移ui(k)

的形式與Ritz法一樣u(k)=u0+Amum,v(k)=v0+Bmvm,w(k)=w0+Cmum但可能位移ui(k)滿足強約束條件。由

=0

得這里

u=umAm,v=vmBm,w=wmCm

2023/2/6113§10-6基于能量原理的近似解法并注意Am、Bm、Cm的任意性以及

由可導出三組方程3m個方程2023/2/6114§10-6基于最小勢能原理上的近似解法3．伽遼金法在平面問題例題與Ritz法類似，不考慮w,而u、v為x、y函數(shù)。伽遼金法在平面問題的求解方程為2m個。對于平面應力問題為2023/2/6115§10-6基于能量原理的近似解法對于平面應變問題將上式中例.已知：2ab薄板（厚度t=1），無體力作用。邊界條件:在x=a和y=0：u=0,v=0;在y=b：u=0,

全部邊界為位移邊界條件，無力的邊界條件。

xyaabb2023/2/6116§10-6基于能量原理的近似解法xyaabb解：不管采用Ritz法或伽遼金法，選取的近似位移場首先要為可能位移：

u(k)=u0+Amum,v(k)=v0+Bmvm,當取m=1時，u=u0+A1u1,v=v0+B1v1,2023/2/6117§10-6基于能量原理的近似解法根據(jù)邊界條件,可選u0=0，而u1和v1在所有邊界上為零。取則2023/2/6118§10-6基于能量原理的近似解法由于全部邊界均為位移邊界，可能位移無需要求其它約束條件。可利用伽遼金法求A1、B1

。由于體力為零，則用伽遼金法的求解方程（2個）：——求出A1和B12023/2/6119§10-6基于能量原理的近似解法，

A1和B1求出后可求應力的近似解。2023/2/6120§10-6基于能量原理的近似解法作業(yè)：1.試寫出伽遼金法在梁彎曲問題的求解方程。

2.利用伽遼金法求圖示簡支梁的近似解，設梁撓度的近似解為v=

B1sinx/l。

qEI

y

x

l2023/2/6121§10-6基于能量原理的近似法6.4基于最小余能原理的近似解法c(k)=Uc(ij(k))+Vc(Xi(k))=c(k)(ij(k))