有限差分法彈性力學(xué)

有限差分法彈性力學(xué)第1頁，課件共30頁，創(chuàng)作于2023年2月

工程中常用的數(shù)值解法有有限單元法和差分法。

有限單元法

是以有限個單元的集合體來代替連續(xù)體，屬于物理上的近似。

差分法

是把彈性力學(xué)的基本方程和邊界條件（一般均為微分方程）近似地改用差分方程（代數(shù)方程）來表示，把求解微分方程的問題改換成為求解代數(shù)方程的問題，屬于數(shù)學(xué)上的近似。第2頁，課件共30頁，創(chuàng)作于2023年2月

第一節(jié)差分方程

第二節(jié)應(yīng)力函數(shù)的差分解

第三節(jié)深梁應(yīng)力函數(shù)的差分解第3頁，課件共30頁，創(chuàng)作于2023年2月第一節(jié)差分方程

差分法是沿用已久的一種數(shù)值解法。隨著計算機(jī)的普及和相應(yīng)的軟件發(fā)展，此法成為解彈性力學(xué)問題的一種有效的方法。第4頁，課件共30頁，創(chuàng)作于2023年2月

我們在彈性體上，用相隔等間距h而平行于坐標(biāo)軸的兩組平行線織成正方形網(wǎng)格，Δx=Δy=h，如圖。

設(shè)f=f(x,y)為彈性體內(nèi)的某一個連續(xù)函數(shù)。該函數(shù)在平行于x軸的一根網(wǎng)線上，如在３-０-１上,它只隨x坐標(biāo)的改變而變化。在鄰近結(jié)點(diǎn)０處，函數(shù)f可展為泰勒級數(shù)如下：第5頁，課件共30頁，創(chuàng)作于2023年2月

我們將只考慮離開結(jié)點(diǎn)０充分近的那些結(jié)點(diǎn)，即（x-x0）充分小。于是可不計（x-x0）的三次及更高次冪的各項，則上式簡寫為：在結(jié)點(diǎn)３，x=x0-h,在結(jié)點(diǎn)1,x=x0+h，代入(b)得：第6頁，課件共30頁，創(chuàng)作于2023年2月聯(lián)立(c),(d),解得差分公式：

同理，在網(wǎng)線４-０-２上可得到差分公式第7頁，課件共30頁，創(chuàng)作于2023年2月

差分公式(１-１)及(１-３)是以相隔2h的兩結(jié)點(diǎn)處的函數(shù)值來表示中間結(jié)點(diǎn)處的一階導(dǎo)數(shù)值，可稱為中點(diǎn)導(dǎo)數(shù)公式。

以相鄰三結(jié)點(diǎn)處的函數(shù)值來表示一個端點(diǎn)處的一階導(dǎo)數(shù)值，可稱為端點(diǎn)導(dǎo)數(shù)公式。

應(yīng)當(dāng)指出：中點(diǎn)導(dǎo)數(shù)公式與端點(diǎn)導(dǎo)數(shù)公式相比，精度較高。因為前者反映了結(jié)點(diǎn)兩邊的函數(shù)變化，而后者卻只反映了結(jié)點(diǎn)一邊的函數(shù)變化。因此，我們總是盡可能應(yīng)用前者，而只有在無法應(yīng)用前者時才不得不應(yīng)用后者。第8頁，課件共30頁，創(chuàng)作于2023年2月

以上(１-１)~(１-４)是基本差分公式,從而可導(dǎo)出其它的差分公式如下：第9頁，課件共30頁，創(chuàng)作于2023年2月第二節(jié)應(yīng)力函數(shù)的差分解

當(dāng)不計體力時，我們已把彈性力學(xué)平面問題歸結(jié)為在給定邊界條件下求解雙調(diào)和方程的問題。用差分法解平面問題，就應(yīng)先將雙調(diào)和方程變換為差分方程，而后求解之。第10頁，課件共30頁，創(chuàng)作于2023年2月

一旦求得彈性體全部節(jié)點(diǎn)的φ值后，就可按應(yīng)力分量差分公式（對節(jié)點(diǎn)0）算得彈性體各節(jié)點(diǎn)的應(yīng)力。第11頁，課件共30頁，創(chuàng)作于2023年2月

可見，用差分法解平面問題，共有兩大任務(wù)：一、建立差分方程將（1-6~8）代入雙調(diào)和方程

對于彈性體邊界以內(nèi)的每一結(jié)點(diǎn)，都可以建立這樣一個差分方程。整理即得第12頁，課件共30頁，創(chuàng)作于2023年2月

二、聯(lián)立求解這些線性代數(shù)方程，就能求得各內(nèi)結(jié)點(diǎn)處的值。

為了求得邊界上各結(jié)點(diǎn)處的φ值，須要應(yīng)用應(yīng)力邊界條件，即：

一般建立和求解差分方程，在數(shù)學(xué)上不會遇到很大困難。但是，當(dāng)對于邊界內(nèi)一行的（距邊界為h的）結(jié)點(diǎn)，建立的差分方程還將涉及邊界上各結(jié)點(diǎn)處的φ值，并包含邊界外一行的虛結(jié)點(diǎn)處的φ值。第13頁，課件共30頁，創(chuàng)作于2023年2月代入上式，即得：

l1=cos(N,x)=cosα=dy/ds,l2=cos(N,y)=sinα=-dx/ds,于是，式（a）可改寫為：由右圖可見，第14頁，課件共30頁，創(chuàng)作于2023年2月關(guān)于邊界上任一點(diǎn)處由此得：

的值，可將（b）式從A點(diǎn)到B點(diǎn)對s積分得到：第15頁，課件共30頁，創(chuàng)作于2023年2月

將此式亦從A點(diǎn)到B點(diǎn)沿s進(jìn)行積分，就得到邊界上任一點(diǎn)B處的φ值。為此利用分部積分法，得：

由高等數(shù)學(xué)可知，第16頁，課件共30頁，創(chuàng)作于2023年2月將式(b),(c)代入，整理得：

由前知，把應(yīng)力函數(shù)加上一個線性函數(shù)，并不影響應(yīng)力。因此，可設(shè)想把應(yīng)力函數(shù)加上a+bx+cy，然后調(diào)整a，b，c三個數(shù)值，使得由式(d)及式(c)可見，設(shè)即可根據(jù)面力分量及導(dǎo)數(shù)求得為已知，第17頁，課件共30頁，創(chuàng)作于2023年2月

從圖易看出，式（2－3）右邊的積分式表示A與B之間的x方向的面力之和；式（2－4）右邊的積分式表示A與B之間的y方向的面力之和；式（2－5）右邊的積分式表示A與B之間的面力對于B點(diǎn)的矩。于是式(d)，式(c)即簡化為：第18頁，課件共30頁，創(chuàng)作于2023年2月

至此，我們解決了怎樣計算邊界上各結(jié)點(diǎn)的值的問題。

至于邊界外一行虛結(jié)點(diǎn)處的值，則可用邊界上結(jié)點(diǎn)處的或值和邊界內(nèi)一行相應(yīng)結(jié)點(diǎn)處的值來表示。例如，對于圖1中的虛結(jié)點(diǎn)14，因為有

第19頁，課件共30頁，創(chuàng)作于2023年2月所以有

當(dāng)求出全部結(jié)點(diǎn)上的φ值以后，我們就可按應(yīng)力分量的差分公式（2－1）計算應(yīng)力分量。

用差分法解彈性平面問題時，可按下列步驟進(jìn)行：（1）在邊界上任意選定一個結(jié)點(diǎn)作為基點(diǎn)A，取第20頁，課件共30頁，創(chuàng)作于2023年2月

然后由面力的矩及面力之和算出邊界上所有各結(jié)點(diǎn)處φ

的值，以及所必需的一些及值，即垂直于邊界方向的導(dǎo)數(shù)值。

（2）應(yīng)用公式（2－6），將邊界外一行虛結(jié)點(diǎn)處的φ值用邊界內(nèi)的相應(yīng)結(jié)點(diǎn)處的φ值來表示。

（3）對邊界內(nèi)的各結(jié)點(diǎn)建立差分方程（2－2），聯(lián)立求解這些結(jié)點(diǎn)處的值。

（4）按照公式（2－6），算出邊界外一行的各虛結(jié)點(diǎn)處的φ值。第21頁，課件共30頁，創(chuàng)作于2023年2月（5）按照公式（2－1）計算應(yīng)力的分量。說明：

1．以上是針對單連體導(dǎo)出的結(jié)果。對于多連體，情況就不象這樣簡單。

2.

如果一部分邊界是曲線的，或是不與坐標(biāo)軸正交，則邊界附近將出現(xiàn)不規(guī)則的內(nèi)結(jié)點(diǎn)。對于這樣的結(jié)點(diǎn)，差分方程（2－2）必須加以修正。第22頁，課件共30頁，創(chuàng)作于2023年2月第三節(jié)例深梁的應(yīng)力函數(shù)差分解

現(xiàn)以如圖所示的混凝土深梁為例，應(yīng)用應(yīng)力函數(shù)的差分解求出應(yīng)力分量。已知混凝土深梁上邊受有均布向下的鉛直荷載q，并由下角點(diǎn)處的反力維持平衡。第23頁，課件共30頁，創(chuàng)作于2023年2月(1）計算邊界上各結(jié)點(diǎn)的φ

、和值。取A為基點(diǎn)，且由上面公式所得的計算結(jié)果見下表。第24頁，課件共30頁，創(chuàng)作于2023年2月（2）計算邊界外以行各虛結(jié)點(diǎn)處的值。由式（2-6）及前表可得第25頁，課件共30頁，創(chuàng)作于2023年2月(3）邊界內(nèi)各結(jié)點(diǎn)的差分方程,由式（2-2）可知第26頁，課件共30頁，創(chuàng)作于2023年2月聯(lián)立求解上式，可得（以qh2為單位)

（4）計算結(jié)點(diǎn)外一行各結(jié)點(diǎn)處的值。由(a)、(b)、(c)可得（5）計算應(yīng)力。對于結(jié)點(diǎn)M，由式（1）可知

第27頁，課件共30頁，創(chuàng)作于2023年2月同理可得

沿著梁的中線MA，的變化如下圖和右圖所示。

第28頁，課件共30頁，創(chuàng)作于2023年2月差分方法的優(yōu)缺點(diǎn)優(yōu)點(diǎn)：

1.差分法是解微分方程邊值問題和彈性力學(xué)問題的有效方法；

2.差分法簡便易行；

3.對于某些結(jié)構(gòu)，為了更精確地分析局部的應(yīng)力狀態(tài)，可以用差分法進(jìn)行分析；缺點(diǎn)：

1.對于曲線邊界和斜邊界等產(chǎn)生的不等間距網(wǎng)格的