結(jié)構(gòu)力學(xué)（ⅠⅠ）第4版 課件全套 第11-13章　矩陣位移法、結(jié)構(gòu)的動(dòng)力計(jì)算、結(jié)構(gòu)的穩(wěn)定計(jì)算

第11章矩陣位移法●本章教學(xué)的基本要求：掌握用矩陣位移法計(jì)算平面桿件結(jié)構(gòu)的原理和方法。包括單元和結(jié)點(diǎn)的劃分；單元和結(jié)構(gòu)坐標(biāo)系中單元?jiǎng)偠染仃嚨男纬桑挥脝卧ㄎ幌蛄啃纬山Y(jié)構(gòu)剛度矩陣；形成結(jié)構(gòu)綜合結(jié)點(diǎn)荷載列陣；結(jié)構(gòu)剛度方程的形成及其求解；結(jié)構(gòu)各桿內(nèi)力和支座約束力的計(jì)算。掌握矩陣位移法的計(jì)算步驟。

●

本章教學(xué)內(nèi)容的重點(diǎn)：用先處理法形成結(jié)構(gòu)剛度矩陣和綜合結(jié)點(diǎn)荷載列陣?！癖菊陆虒W(xué)內(nèi)容的難點(diǎn)：用先處理法形成結(jié)構(gòu)剛度矩陣中各步驟的物理意義；單元?jiǎng)偠染仃嚭徒Y(jié)構(gòu)剛度矩陣中剛度系數(shù)的物理意義和求法；矩陣位移法與位移法之間的聯(lián)系與區(qū)別?！?/p>

本章內(nèi)容簡介:11.1概述11.2桿件結(jié)構(gòu)的離散化11.3單元坐標(biāo)系中的單元?jiǎng)偠染仃?1.4結(jié)構(gòu)坐標(biāo)系中的單元?jiǎng)偠染仃?1.5用先處理法形成結(jié)構(gòu)剛度矩陣11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣11.7求解結(jié)點(diǎn)位移和單元桿端力11.8矩陣位移法的計(jì)算步驟和算例第11章矩陣位移法11.1概述11.1.1結(jié)構(gòu)矩陣分析結(jié)構(gòu)矩陣分析以結(jié)構(gòu)力學(xué)的原理為基礎(chǔ)，采用矩陣進(jìn)行運(yùn)算，公式緊湊，形式統(tǒng)一，便于計(jì)算過程的規(guī)格化和程序化，滿足電子計(jì)算機(jī)進(jìn)行自動(dòng)化計(jì)算的要求。11.1.2結(jié)構(gòu)矩陣分析的兩類方法力法位移法傳統(tǒng)結(jié)構(gòu)分析方法結(jié)構(gòu)矩陣分析方法矩陣力法（柔度法）矩陣位移法（剛度法）結(jié)構(gòu)矩陣分析有時(shí)也稱為桿件結(jié)構(gòu)的有限單元法。11.1.3矩陣位移法的三個(gè)基本環(huán)節(jié)矩陣位移法就是以矩陣形式表達(dá)的位移法。1.結(jié)構(gòu)的離散化2.單元分析3.整體分析把結(jié)構(gòu)劃分為有限個(gè)較小的單元，各單元只在有限個(gè)結(jié)點(diǎn)處相連。分析桿單元的桿端內(nèi)力與桿端位移之間的關(guān)系，以矩陣形式表示，建立單元?jiǎng)偠确匠?。通過考慮各結(jié)點(diǎn)的變形協(xié)調(diào)條件和平衡條件，建立整個(gè)結(jié)構(gòu)的剛度方程，以求解原結(jié)構(gòu)的結(jié)點(diǎn)位移。11.1概述11.2桿件結(jié)構(gòu)的離散化11.2.1單元與結(jié)點(diǎn)的劃分和編號(hào)1.單元與結(jié)點(diǎn)的劃分

基本要求：保證單元為等截面直桿。此外，應(yīng)兼顧計(jì)算的方便和結(jié)構(gòu)的特殊形式。2.單元與結(jié)點(diǎn)的編號(hào)

單元號(hào)用①、②、③、…表示；結(jié)點(diǎn)號(hào)用1、2、3、…表示。11.2.2兩種直角坐標(biāo)系矩陣分析中，為了矢量分析的方便，需要采用兩種直角坐標(biāo)系，即結(jié)構(gòu)坐標(biāo)系和單元坐標(biāo)系。采用右手旋轉(zhuǎn)直角坐標(biāo)系，如圖所示。11.2桿件結(jié)構(gòu)的離散化11.2.3單元桿端力和桿端位移的表示方法單元桿端截面的內(nèi)力和位移分別稱為單元桿端力和桿端位移。平面剛架中的第e個(gè)單元，其始端和末端的結(jié)點(diǎn)編號(hào)分別為i和j。11.2桿件結(jié)構(gòu)的離散化1.單元坐標(biāo)系中的單元桿端力和桿端位移在單元分析時(shí)，各桿端力和桿端位移均應(yīng)按照一定的次序排列，一般規(guī)定先“始端”后“末端”。11.2桿件結(jié)構(gòu)的離散化2.結(jié)構(gòu)坐標(biāo)系中的單元桿端力和桿端位移桿端線位移與結(jié)構(gòu)坐標(biāo)系正向一致為正，桿端轉(zhuǎn)角以順時(shí)針方向?yàn)檎?1.2桿件結(jié)構(gòu)的離散化11.2.4單元坐標(biāo)轉(zhuǎn)換為了便于利用單元坐標(biāo)系中的單元桿端力和桿端位移來建立結(jié)構(gòu)坐標(biāo)系中的結(jié)構(gòu)剛度方程，有必要建立單元桿端力和桿端位移在兩種坐標(biāo)系之間的轉(zhuǎn)換關(guān)系。11.2桿件結(jié)構(gòu)的離散化寫成矩陣形式，則有（11-5）（11-6）11.2桿件結(jié)構(gòu)的離散化其中（11-7）稱為平面剛架單元坐標(biāo)轉(zhuǎn)換矩陣，它是正交矩陣，因而滿足（11-8）不同坐標(biāo)系下，單元桿端力轉(zhuǎn)換關(guān)系如下11.2桿件結(jié)構(gòu)的離散化（11-6）（11-9）（11-10）（11-11）單元桿端位移與單元桿端力一樣，同為矢量，故不同坐標(biāo)系下單元桿端位移的轉(zhuǎn)換關(guān)系與單元桿端力的相同，具體如下11.2桿件結(jié)構(gòu)的離散化11.3單元坐標(biāo)系中的單元?jiǎng)偠染仃?1.3.1一般單元的單元?jiǎng)偠确匠毯蛦卧獎(jiǎng)偠染仃噯卧獥U端力和桿端位移之間的轉(zhuǎn)換關(guān)系，稱為單元?jiǎng)偠确匠?。單元?jiǎng)偠染仃囀菞U端力與桿端位移之間的轉(zhuǎn)換矩陣。在單元坐標(biāo)系中，單元?jiǎng)偠确匠炭杀硎緸闉閱卧鴺?biāo)系中的單元?jiǎng)偠染仃嚒?/p>

式中，在線彈性小變形范圍內(nèi)，可忽略軸向變形與彎曲變形之間的相互影響，根據(jù)桿件的拉壓胡克定律和無荷載作用時(shí)的轉(zhuǎn)角位移方程，可寫出關(guān)系式

11.3.1一般單元的單元?jiǎng)偠确匠毯蛦卧獎(jiǎng)偠染仃?1.3單元坐標(biāo)系中的單元?jiǎng)偠染仃噷懗删仃囆问胶唽憺榉Q為單元坐標(biāo)系中平面剛架一般單元的單元?jiǎng)偠确匠獭?1.3單元坐標(biāo)系中的單元?jiǎng)偠染仃嚻渲蟹Q為單元坐標(biāo)系中平面剛架一般單元的單元?jiǎng)偠染仃?，簡稱單剛。是6×6階的方陣。11.3單元坐標(biāo)系中的單元?jiǎng)偠染仃?1.3.2單元?jiǎng)偠染仃嚨男再|(zhì)1.單剛是單元固有的性質(zhì)單剛中各元素只與單元的彈性模量E、橫截面面積A、慣性矩I及桿長l等有關(guān)，而與外荷載等其他因素?zé)o關(guān)。

2.單元?jiǎng)偠认禂?shù)的物理意義單元?jiǎng)偠染仃囍械拿總€(gè)元素稱為單元?jiǎng)偠认禂?shù)，代表由于單位桿端位移引起的桿端力。1）任一個(gè)元素表示第l個(gè)桿端位移分量等于1（其余位移分量等于零）時(shí)，所引起的第m個(gè)桿端力分量的值。

11.3單元坐標(biāo)系中的單元?jiǎng)偠染仃?1.3.2單元?jiǎng)偠染仃嚨男再|(zhì)2）某一列的六個(gè)元素，表示當(dāng)該列對(duì)應(yīng)桿端位移分量等于1（其余位移分量等于零）時(shí)，所引起的六個(gè)桿端力分量。3）某一行的六個(gè)元素，分別表示各個(gè)桿端位移分量分別等于1時(shí)，所引起的按該行號(hào)順序排列的那個(gè)桿端力分量的數(shù)值。單剛中第五列元素11.3單元坐標(biāo)系中的單元?jiǎng)偠染仃?.單元?jiǎng)偠染仃囀菍?duì)稱矩陣11.3.2單元?jiǎng)偠染仃嚨男再|(zhì)單元?jiǎng)偠染仃囍形挥谥鲗?duì)角線兩邊處于對(duì)稱位置上的兩個(gè)元素（簡稱副系數(shù)或副元）相等。

4.一般單元的單元?jiǎng)偠染仃囀瞧娈惥仃嚨男辛惺街档扔诹?，不存在逆矩陣?1.3單元坐標(biāo)系中的單元?jiǎng)偠染仃?.單元?jiǎng)偠仍刂档男再|(zhì)1）主系數(shù)恒大于零。2）第二列和第五列各非零元素等值反號(hào)。3）第三列和第六列相比，只需將第三行和第六行元素交換。11.3.2單元?jiǎng)偠染仃嚨男再|(zhì)如果已知單元結(jié)點(diǎn)位移，可惟一求出單元桿端力；反之，則沒有惟一解。無外界約束的一般單元，也稱自由單元。任一自由單元還可以產(chǎn)生任意的剛體位移，故某一組滿足平衡條件的桿端力可以與彈性位移和任意剛體位移組成的多組桿端位移相對(duì)應(yīng)。11.3單元坐標(biāo)系中的單元?jiǎng)偠染仃嚥捎蒙鲜鐾瑯拥姆椒ㄟ€可以得出其他各種特殊單元（如一端固定而另一端為鉸支或定向支承或自由的梁單元等）的剛度方程和剛度矩陣。在結(jié)構(gòu)矩陣分析中，我們著眼于計(jì)算過程的程序化、標(biāo)準(zhǔn)化和自動(dòng)化，因而，無論i、j兩端的約束情況如何，都可按一般單元情況處理，即采用一種標(biāo)準(zhǔn)化形式（一般單元）的剛度矩陣。這時(shí)，需要把實(shí)際鉸支或自由端未知的位移作為求解的未知位移。此外，在采用計(jì)算機(jī)進(jìn)行分析時(shí)，一般都把桿件的軸向變形影響考慮進(jìn)來，使程序編制上更為簡單，而且對(duì)于大多數(shù)情況，這樣得到的結(jié)果更為準(zhǔn)確。11.3單元坐標(biāo)系中的單元?jiǎng)偠染仃嚕ǚQ為結(jié)構(gòu)坐標(biāo)系中的單元?jiǎng)偠确匠蹋??？傻梅Q為結(jié)構(gòu)坐標(biāo)系中的單元?jiǎng)偠染仃嚕匀皇菍?duì)稱矩陣；由于仍為自由單元，故仍然是奇異矩陣。其中元素的值還與結(jié)構(gòu)坐標(biāo)系到單元坐標(biāo)系之間的夾角a有關(guān)。a以從結(jié)構(gòu)坐標(biāo)系x軸順時(shí)針旋轉(zhuǎn)到單元坐標(biāo)系軸方向?yàn)檎?1.4結(jié)構(gòu)坐標(biāo)系中的單元?jiǎng)偠染仃嚫鶕?jù)前述桿端力和桿端位移的坐標(biāo)轉(zhuǎn)換關(guān)系，以及單元坐標(biāo)系中的單剛方程，可得比對(duì)平面剛架一般單元的單元?jiǎng)偠染仃噷?duì)稱

11.4結(jié)構(gòu)坐標(biāo)系中的單元?jiǎng)偠染仃?1.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5.1結(jié)構(gòu)的整體分析結(jié)構(gòu)剛度方程可寫作

D為待求的基本未知量（結(jié)點(diǎn)位移）列陣；

F為結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣（也稱總荷載列陣）；

K為結(jié)構(gòu)剛度矩陣，簡稱總剛，是結(jié)構(gòu)各結(jié)點(diǎn)力與結(jié)點(diǎn)位移之間的轉(zhuǎn)換矩陣。1.先處理法所謂先處理法，就是一開始即考慮結(jié)構(gòu)的支承條件，把已知的支座位移排除在基本未知量之外，相應(yīng)地，結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣中也不包括支約束力。因而，所形成的結(jié)構(gòu)剛度方程階數(shù)小，不用再修正。先處理法可很方便地處理有鉸結(jié)點(diǎn)的結(jié)構(gòu)、各種性質(zhì)的支承結(jié)點(diǎn)結(jié)構(gòu)及忽略軸向變形的結(jié)構(gòu)。11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5.2先處理法和后處理法在整體分析時(shí)，必須考慮和引入支承邊界條件，而這一條件可以在形成整體剛度方程之前或之后處理，因而形成了先處理法和后處理法兩種做法。11.5.2先處理法和后處理法在整體分析時(shí)，必須考慮和引入支承邊界條件，而這一條件可以在形成整體剛度方程之前或之后處理，因而形成了先處理法和后處理法兩種做法。2.后處理法所謂后處理法，就是先不考慮支承條件，把已知的支座結(jié)點(diǎn)位移和自由結(jié)點(diǎn)的位移也一并列入結(jié)構(gòu)的結(jié)點(diǎn)位移列陣中，形成不受約束的原始剛度方程；然后，再根據(jù)結(jié)點(diǎn)的實(shí)際支承條件修正原始剛度方程，形成結(jié)構(gòu)剛度方程，以求解結(jié)點(diǎn)位移。后處理法的基本未知量數(shù)目更多，占用計(jì)算機(jī)的內(nèi)存和機(jī)時(shí)更多。一般用于結(jié)點(diǎn)多而支座約束少、考慮軸向變形的結(jié)構(gòu)。11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5.3結(jié)點(diǎn)位移分量的統(tǒng)一編碼—總碼當(dāng)采用一般單元考慮單元的軸向變形時(shí)，結(jié)構(gòu)的每個(gè)結(jié)點(diǎn)一般有三個(gè)位移（u,v,q）。采用先處理法計(jì)算，首先應(yīng)按照結(jié)構(gòu)的結(jié)點(diǎn)順序，依次對(duì)每個(gè)結(jié)點(diǎn)的未知位移分量u、v、q統(tǒng)一編碼（又稱為結(jié)點(diǎn)未知量編碼）。1）對(duì)于取作基本未知量的有效結(jié)點(diǎn)未知量，按照結(jié)點(diǎn)順序（每個(gè)結(jié)點(diǎn)又按照u、v、q

的順序），依次連續(xù)編碼即1，2，…，n。11.5用先處理法形成結(jié)構(gòu)剛度矩陣2）對(duì)于已知的支座結(jié)點(diǎn)位移分量（包括為零和非零兩種情況），編“0”號(hào)。11.5用先處理法形成結(jié)構(gòu)剛度矩陣3）對(duì)于平面剛架內(nèi)部有鉸結(jié)點(diǎn)（全鉸或半鉸結(jié)點(diǎn)）的情況，因采用一般單元來分析，桿在鉸端的轉(zhuǎn)角也是未知量，故應(yīng)將相互鉸結(jié)的桿端編以不同的結(jié)點(diǎn)號(hào)（即進(jìn)行雙編號(hào)）。11.5用先處理法形成結(jié)構(gòu)剛度矩陣4）對(duì)于無效結(jié)點(diǎn)未知量（即計(jì)算中不需考慮的結(jié)點(diǎn)未知量），編“0”號(hào)11.5用先處理法形成結(jié)構(gòu)剛度矩陣5）對(duì)于忽略軸向變形的剛架，每個(gè)剛結(jié)點(diǎn)在某一方向上的位移分量不一定都是獨(dú)立的。11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5.4單元定位向量1.定義將單元始末兩端的結(jié)點(diǎn)位移分量的統(tǒng)一編碼（即結(jié)點(diǎn)未知量編碼）按順序排列所組成的列向量稱為單元定位向量。11.5用先處理法形成結(jié)構(gòu)剛度矩陣2.單元定位向量的作用單元定位向量不但可以用來確定單元?jiǎng)偠染仃嚫髟卦诮Y(jié)構(gòu)剛度矩陣K中的位置，以及單元桿端位移向量各元素在整個(gè)結(jié)構(gòu)的結(jié)點(diǎn)位移向量Δ中的位置，而且，可以用來確定單元等效結(jié)點(diǎn)荷載（由單元上非結(jié)點(diǎn)荷載等效轉(zhuǎn)換而得到的結(jié)點(diǎn)荷載）各元素在結(jié)構(gòu)綜合結(jié)點(diǎn)荷載列陣P中的位置。3.單元定位向量的優(yōu)點(diǎn)借助于單元定位向量，可以很方便地解決計(jì)算機(jī)自動(dòng)化計(jì)算中所需的信息。11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5.5按照直接平衡法形成結(jié)構(gòu)剛度方程和結(jié)構(gòu)剛度矩陣圖a)所示在結(jié)點(diǎn)荷載作用下的平面剛架，討論結(jié)構(gòu)剛度方程的建立。該剛架有5個(gè)結(jié)點(diǎn)、8個(gè)結(jié)點(diǎn)位移未知量，在結(jié)點(diǎn)1、結(jié)點(diǎn)2（或3）和結(jié)點(diǎn)4上作用有結(jié)點(diǎn)荷載。結(jié)構(gòu)的結(jié)點(diǎn)位移列陣和結(jié)點(diǎn)荷載列陣分別為11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5用先處理法形成結(jié)構(gòu)剛度矩陣圖b)為結(jié)構(gòu)坐標(biāo)系中單元隔離體和結(jié)點(diǎn)隔離體的受力情況。

11.5用先處理法形成結(jié)構(gòu)剛度矩陣根據(jù)各結(jié)點(diǎn)的平衡條件，可得如下方程11.5用先處理法形成結(jié)構(gòu)剛度矩陣根據(jù)變形協(xié)調(diào)條件，單元桿端位移應(yīng)等于與之相連的結(jié)點(diǎn)的結(jié)點(diǎn)位移，由圖a)，可得于是，用結(jié)構(gòu)的結(jié)點(diǎn)位移表示的結(jié)構(gòu)坐標(biāo)系中的單元①至單元③的剛度方程，分別如下列各式所示。11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5用先處理法形成結(jié)構(gòu)剛度矩陣展開以上三式，就可得到用結(jié)構(gòu)的結(jié)點(diǎn)位移表示的單元桿端力。例如

11.5用先處理法形成結(jié)構(gòu)剛度矩陣將上式代入，可得同理，可得11.5用先處理法形成結(jié)構(gòu)剛度矩陣將上兩式匯總寫成矩陣形式，得11.5用先處理法形成結(jié)構(gòu)剛度矩陣上式就是圖a)所示平面剛架的結(jié)構(gòu)剛度方程。與KD=F比較可知，其左邊的8×8階矩陣就是結(jié)構(gòu)剛度矩陣K，即11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5.6按照直接剛度法形成結(jié)構(gòu)剛度矩陣結(jié)構(gòu)剛度矩陣的元素是由單剛元素按照一定規(guī)律組成的，只要確定了單剛元素在結(jié)構(gòu)剛度矩陣中的位置，就可以由各單元的單剛直接集成結(jié)構(gòu)剛度矩陣K。利用單元定位向量確定單剛元素在結(jié)構(gòu)剛度矩陣中的行碼和列碼后，直接將單剛元素送入結(jié)構(gòu)總剛中的對(duì)應(yīng)位置，這種裝配結(jié)構(gòu)總剛的方法，稱為直接剛度法。11.5用先處理法形成結(jié)構(gòu)剛度矩陣1.裝配過程1)計(jì)算單元e在結(jié)構(gòu)坐標(biāo)系中的剛度矩陣Ke，將其定位向量中各分量（始末兩端結(jié)點(diǎn)未知量編碼——簡稱“總碼”）及單元自身的結(jié)點(diǎn)位移分量編碼——簡稱“局部碼”,分別寫在Ke的上方和右側(cè)。2)按照單元定位向量中的非零分量所指定的行碼和列碼，將各單元?jiǎng)偠染仃囍械脑?，正確地疊加到結(jié)構(gòu)剛度矩陣K中去，行、列碼相同的元素則相加。這一作法稱為對(duì)號(hào)入座，同號(hào)相加。11.5用先處理法形成結(jié)構(gòu)剛度矩陣?yán)?，上列剛架各單元的定位向量分別為

11.5用先處理法形成結(jié)構(gòu)剛度矩陣按照單元定位向量的指引，分別送入8×8階的結(jié)構(gòu)總剛中“對(duì)號(hào)入座”，這樣形成的K1、K2和K3矩陣，稱為每個(gè)單元對(duì)總剛的貢獻(xiàn)矩陣。11.5用先處理法形成結(jié)構(gòu)剛度矩陣顯然，結(jié)構(gòu)總剛K與三個(gè)單元貢獻(xiàn)矩陣的關(guān)系為K=K1+K2+K3將三個(gè)單元貢獻(xiàn)矩陣中相同行、列的各元素“同號(hào)相加”，即可最后裝配總剛。11.5用先處理法形成結(jié)構(gòu)剛度矩陣2.直接剛度法的正確性說明總剛K中任一元素kij，表示結(jié)點(diǎn)位移列陣D中第j個(gè)結(jié)點(diǎn)位移分量dj=1（其余結(jié)點(diǎn)位移分量均為零）時(shí)，所引起的結(jié)點(diǎn)荷載列陣F中第i個(gè)結(jié)點(diǎn)力分量fi之值。

例如，前述例子結(jié)構(gòu)總剛的剛度系數(shù)k12，表示當(dāng)結(jié)點(diǎn)位移分量d2=1（其余結(jié)點(diǎn)位移分量為零）時(shí)，所引起的結(jié)點(diǎn)荷載列陣F中第1個(gè)結(jié)點(diǎn)力分量f1之值，如下圖所示。11.5用先處理法形成結(jié)構(gòu)剛度矩陣2.直接剛度法的正確性說明總剛K中任一元素kij，表示結(jié)點(diǎn)位移列陣D中第j個(gè)結(jié)點(diǎn)位移分量dj=1（其余結(jié)點(diǎn)位移分量均為零）時(shí)，所引起的結(jié)點(diǎn)荷載列陣F中第i個(gè)結(jié)點(diǎn)力分量fi之值。

例如，前述例子結(jié)構(gòu)總剛的剛度系數(shù)k12，表示當(dāng)結(jié)點(diǎn)位移分量d2=1（其余結(jié)點(diǎn)位移分量為零）時(shí)，所引起的結(jié)點(diǎn)荷載列陣F中第1個(gè)結(jié)點(diǎn)力分量f1之值，如下圖所示。11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5用先處理法形成結(jié)構(gòu)剛度矩陣12當(dāng)d2=1時(shí)，則應(yīng)有單元①的d5(1)=1和單元②的d2(1)=1，分別如圖b)、c)所示。再由結(jié)點(diǎn)1水平方向的受力（圖d所示）及其平衡條件∑Fx=0，可得k12=k45(1)+k12(1)，與之前推導(dǎo)的結(jié)果一致。3.結(jié)構(gòu)剛度矩陣K中的元素的組成規(guī)律如果將單元定位向量中有編號(hào)i的所有單元稱為未知量di的相關(guān)單元，將與未知量di同屬一個(gè)單元的其他未知量稱為未知量di的相關(guān)未知量，則結(jié)構(gòu)剛度矩陣K具有如下組成規(guī)律：

11.5用先處理法形成結(jié)構(gòu)剛度矩陣(1)主對(duì)角線元素（簡稱主元）kii

kii由未知量di的相關(guān)單元的單剛中的相應(yīng)主對(duì)角線元素疊加而成。

(2)非主對(duì)角線元素kij（i≠j）kij有兩種情況：若di和dj是相關(guān)未知量，則kij=kji≠0；若di和dj不是相關(guān)未知量，則kij=kji=0

。11.5用先處理法形成結(jié)構(gòu)剛度矩陣1)結(jié)構(gòu)剛度矩陣K是一個(gè)N×N階的方陣，N為結(jié)點(diǎn)的位移未知量數(shù)。2)結(jié)構(gòu)剛度矩陣K是一個(gè)對(duì)稱矩陣，K=K

T。這可由約束力互等定理證明K中對(duì)稱于主對(duì)角線的元素兩兩相等，即kij=kji。3)結(jié)構(gòu)剛度矩陣K是正定的。即|K|>0，且kii>0。

4)結(jié)構(gòu)剛度矩陣K是一個(gè)帶狀矩陣，即非零元素分布在主對(duì)角線的附近。因此，在同一單元內(nèi)的未知量編碼差值應(yīng)盡量保持最小值。在實(shí)際結(jié)構(gòu)分析中所遇到的結(jié)構(gòu)剛度矩陣一般都是大型稀疏矩陣。非零元素很少，往往只占10%左右。11.5.7結(jié)構(gòu)剛度矩陣K的性質(zhì)11.5用先處理法形成結(jié)構(gòu)剛度矩陣

【例11-1】試形成圖示剛架的結(jié)構(gòu)剛度矩陣。已知EA=4.8×106kN，EI=0.9×105kN·m2。

解：(1)結(jié)構(gòu)離散化

(2)形成結(jié)構(gòu)坐標(biāo)系中的單剛11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.5用先處理法形成結(jié)構(gòu)剛度矩陣

(3)形成結(jié)構(gòu)剛度矩陣（對(duì)號(hào)入座，同號(hào)相加）11.5用先處理法形成結(jié)構(gòu)剛度矩陣11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣11.6.1形成單元的等效結(jié)點(diǎn)荷載列陣11.6.2計(jì)算結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣F

1.形成各單元在結(jié)構(gòu)坐標(biāo)系中的等效結(jié)點(diǎn)荷載列陣1)計(jì)算在單元坐標(biāo)系中的單元固端約束力2)計(jì)算在單元坐標(biāo)系中的單元等效結(jié)點(diǎn)荷載將固端約束力反號(hào)，得到單元坐標(biāo)系中的單元等效結(jié)點(diǎn)荷載3)計(jì)算在結(jié)構(gòu)坐標(biāo)系中的單元等效結(jié)點(diǎn)荷載

即固定結(jié)點(diǎn)即放松結(jié)點(diǎn)即坐標(biāo)轉(zhuǎn)換11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣2.“對(duì)號(hào)入座，同號(hào)疊加”形成結(jié)構(gòu)等效結(jié)點(diǎn)荷載將各單元在結(jié)構(gòu)坐標(biāo)系中的等效結(jié)點(diǎn)荷載列陣元素，按其單元定位向量的指示，裝入結(jié)構(gòu)等效結(jié)點(diǎn)荷載列陣FE中。若不同單元有裝入相同位置的元素，則進(jìn)行疊加。3.形成結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣將結(jié)構(gòu)原本所受的結(jié)點(diǎn)荷載所構(gòu)成的列陣（稱為結(jié)構(gòu)直接結(jié)點(diǎn)荷載FJ），與結(jié)構(gòu)等效結(jié)點(diǎn)荷載列陣FE疊加，形成結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣。F=FJ+FE11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣【例11-2】試求圖示結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣F。解：(1)計(jì)算單元固端約束力，，11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣(2)計(jì)算單元等效結(jié)點(diǎn)荷載

單元①：11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣單元②：11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣單元③：11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣(3)利用單元定位向量形成結(jié)構(gòu)的等效結(jié)點(diǎn)荷載列陣FE按單元定位向量，“對(duì)號(hào)入座，同號(hào)相加”

11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣(4)形成結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣F11.6結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣11.7求解結(jié)點(diǎn)位移和單元桿端力結(jié)構(gòu)剛度方程（即位移法的基本方程）11.7.1求解結(jié)點(diǎn)位移D用先處理法直接形成的結(jié)構(gòu)剛度方程，是一個(gè)線性代數(shù)方程組，其結(jié)構(gòu)剛度矩陣K為對(duì)稱正定矩陣。求解此方程組，即可得到未知結(jié)點(diǎn)位移D的惟一確定解。11.7.2求單元桿端位移11.7.3求單元桿端內(nèi)力應(yīng)該注意桿端內(nèi)力符號(hào)與傳統(tǒng)結(jié)構(gòu)分析中的差異性11.7求解結(jié)點(diǎn)位移和單元桿端力按上兩式求得單元桿端力中，桿端軸力與軸正向同向?yàn)檎?，桿端剪力與軸正向同向?yàn)檎?，桿端彎矩以順時(shí)針方向?yàn)檎?。這與本書前面各章節(jié)中內(nèi)力正負(fù)號(hào)規(guī)定不盡相同。注意在繪內(nèi)力圖時(shí)，仍應(yīng)按原規(guī)定的內(nèi)力正負(fù)號(hào)進(jìn)行繪制，即軸力以受拉為正，剪力以繞它端順時(shí)針旋轉(zhuǎn)為正，彎矩不區(qū)分正負(fù)，而應(yīng)繪在桿件受拉一側(cè)。

11.7求解結(jié)點(diǎn)位移和單元桿端力11.8矩陣位移法的計(jì)算步驟和算例11.8.1先處理法的計(jì)算步驟

1）對(duì)結(jié)構(gòu)進(jìn)行離散化。建立單元坐標(biāo)系和結(jié)構(gòu)坐標(biāo)系，并對(duì)結(jié)點(diǎn)、單元及結(jié)點(diǎn)位移分量分別進(jìn)行編號(hào)。2）形成單元坐標(biāo)系中的單元?jiǎng)偠染仃嚒?）形成結(jié)構(gòu)坐標(biāo)系中的單元?jiǎng)偠染仃嚒?）按直接剛度法形成結(jié)構(gòu)剛度矩陣。5）形成結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣。6）求解未知結(jié)點(diǎn)位移列陣。7）計(jì)算桿端內(nèi)力，繪內(nèi)力圖。1.平面剛架【例11-3】試用先處理法計(jì)算圖示剛架的內(nèi)力。已知各桿EA=4.8×106kN，EI=0.9×105kN·m2。

解：(1)形成結(jié)構(gòu)剛度矩陣（參見例11-1）11.8.2舉例11.8矩陣位移法的計(jì)算步驟和算例12345612345611.8矩陣位移法的計(jì)算步驟和算例(2)計(jì)算結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣（參見例11-2）(3)解結(jié)構(gòu)剛度方程，求結(jié)點(diǎn)位移列陣11.8矩陣位移法的計(jì)算步驟和算例解得結(jié)點(diǎn)位移為11.8矩陣位移法的計(jì)算步驟和算例(4)計(jì)算各單元桿端力單元①11.8矩陣位移法的計(jì)算步驟和算例單元②11.8矩陣位移法的計(jì)算步驟和算例單元③11.8矩陣位移法的計(jì)算步驟和算例(7)根據(jù)上述計(jì)算結(jié)果作內(nèi)力圖，如下圖所示。11.8矩陣位移法的計(jì)算步驟和算例2.連續(xù)梁單元連續(xù)梁離散得出的單元（通常忽略軸向變形），只有兩端的彎矩和與轉(zhuǎn)角存在，其余位移為零，即其單元?jiǎng)偠确匠虨?1.8矩陣位移法的計(jì)算步驟和算例故T=I，且其單剛為對(duì)于連續(xù)梁，各單元坐標(biāo)軸與結(jié)構(gòu)坐標(biāo)軸方向一致，無需坐標(biāo)轉(zhuǎn)換。不計(jì)軸向變形的無結(jié)點(diǎn)線位移剛架，也可以使用本類連續(xù)梁單元進(jìn)行計(jì)算。11.8矩陣位移法的計(jì)算步驟和算例【例11-4】試用先處理法計(jì)算圖示連續(xù)梁的內(nèi)力，并作彎矩圖。忽略各桿軸向變形的影響。解：(1)結(jié)構(gòu)坐標(biāo)系、單元坐標(biāo)系、結(jié)點(diǎn)編號(hào)、單元編號(hào)及結(jié)點(diǎn)位移未知量編碼如圖所示。(2)形成整體坐標(biāo)系中的單元?jiǎng)偠染仃嚒?1.8矩陣位移法的計(jì)算步驟和算例1211.8矩陣位移法的計(jì)算步驟和算例(3)利用單元定位向量形成結(jié)構(gòu)剛度矩陣K

(4)形成結(jié)構(gòu)綜合結(jié)點(diǎn)荷載列陣P11.8矩陣位移法的計(jì)算步驟和算例再利用單元定位向量形成結(jié)構(gòu)綜合結(jié)點(diǎn)荷載列陣(5)形成結(jié)構(gòu)剛度方程，解方程求結(jié)點(diǎn)位移11.8矩陣位移法的計(jì)算步驟和算例解方程，得結(jié)點(diǎn)位移D為(6)利用單元定位向量從結(jié)點(diǎn)位移D中取出各單元桿端位移11.8矩陣位移法的計(jì)算步驟和算例11.8矩陣位移法的計(jì)算步驟和算例M圖(kN.m)(7)作彎矩圖11.8矩陣位移法的計(jì)算步驟和算例3.桁架單元對(duì)于理想桁架，各桿件只有軸向變形，在單元坐標(biāo)系中桿端位移僅有軸向位移的存在，而其它位移為零，即其單元?jiǎng)偠确匠虨閱卧鴺?biāo)系中的單剛為11.8矩陣位移法的計(jì)算步驟和算例對(duì)稱坐標(biāo)轉(zhuǎn)換矩陣為結(jié)構(gòu)坐標(biāo)系中的單剛為11.8矩陣位移法的計(jì)算步驟和算例形成平面桁架結(jié)構(gòu)剛度矩陣K的方法與平面剛架相同。對(duì)結(jié)點(diǎn)位移分量編碼應(yīng)注意，桁架單元的結(jié)點(diǎn)角位移不作為基本未知量?！纠?1-5】試用先處理法計(jì)算圖示桁架的內(nèi)力。各桿EA相同。

解：(1)結(jié)構(gòu)坐標(biāo)系、單元坐標(biāo)系、結(jié)點(diǎn)編號(hào)、單元編號(hào)及結(jié)點(diǎn)位移未知量編碼如圖所示。11.8矩陣位移法的計(jì)算步驟和算例(2)形成整體坐標(biāo)系中的單元?jiǎng)偠染仃噯卧俸蛦卧?a=p/2單元①單元③11.8矩陣位移法的計(jì)算步驟和算例單元②和單元④:a=0單元②單元④11.8矩陣位移法的計(jì)算步驟和算例單元⑤:a=p

/4單元⑥:a=3p

/411.8矩陣位移法的計(jì)算步驟和算例(3)利用單元定位向量形成結(jié)構(gòu)剛度矩陣K(4)形成結(jié)構(gòu)的綜合結(jié)點(diǎn)荷載列陣P11.8矩陣位移法的計(jì)算步驟和算例(5)形成結(jié)構(gòu)剛度方程解方程，得結(jié)點(diǎn)位移D為11.8矩陣位移法的計(jì)算步驟和算例(6)計(jì)算各單元軸力11.8矩陣位移法的計(jì)算步驟和算例11.8矩陣位移法的計(jì)算步驟和算例FN圖(kN)11.8矩陣位移法的計(jì)算步驟和算例第12章結(jié)構(gòu)的動(dòng)力計(jì)算●

本章教學(xué)的基本要求：掌握動(dòng)力分析的基本方法及體系動(dòng)力自由度數(shù)的判別方法；掌握單自由度和兩個(gè)自由度體系運(yùn)動(dòng)方程的建立方法，及其自由振動(dòng)和在簡諧荷載作用下強(qiáng)迫振動(dòng)的計(jì)算方法；了解阻尼的作用；了解多自由度體系在一般動(dòng)力荷載作用下的強(qiáng)迫振動(dòng)；了解頻率的近似計(jì)算方法。●

本章教學(xué)內(nèi)容的重點(diǎn)：動(dòng)力自由度數(shù)的判別方法；單自由度體系運(yùn)動(dòng)方程的建立；單自由度及有限自由度（重點(diǎn)是兩個(gè)自由度）體系動(dòng)力特性的計(jì)算；單自由度、有限自由度體系在簡諧荷載作用下內(nèi)力、位移的計(jì)算；阻尼對(duì)振動(dòng)的影響?！?/p>

本章教學(xué)內(nèi)容的難點(diǎn)：用剛度法和柔度法建立單自由度體系的運(yùn)動(dòng)方程；在動(dòng)力特性和動(dòng)力反應(yīng)計(jì)算中，剛度系數(shù)和柔度系數(shù)的計(jì)算；單自由度和兩個(gè)自由度體系在簡諧荷載作用下動(dòng)力反應(yīng)的計(jì)算。12.1概述12.2單自由度體系的運(yùn)動(dòng)方程12.3單自由度體系的自由振動(dòng)12.4單自由度體系的強(qiáng)迫振動(dòng)12.5阻尼對(duì)振動(dòng)的影響12.6多自由度體系的自由振動(dòng)12.7主振型的正交性12.8多自由度體系在簡諧荷載作用下的強(qiáng)迫振動(dòng)12.9多自由度體系在任意動(dòng)力荷載作用下的強(qiáng)迫振動(dòng)12.10*無限自由度體系的自由振動(dòng)12.11近似法計(jì)算自振頻率第12章結(jié)構(gòu)的動(dòng)力計(jì)算●

本章內(nèi)容簡介:12.1概述12.1.1結(jié)構(gòu)動(dòng)力計(jì)算的任務(wù)1.基本任務(wù)結(jié)構(gòu)靜力學(xué)：主要研究結(jié)構(gòu)在靜力荷載作用下的靜力反應(yīng)?！耢o內(nèi)力●靜位移（均有唯一性解答）作為結(jié)構(gòu)設(shè)計(jì)的依據(jù)結(jié)構(gòu)動(dòng)力學(xué)：主要研究結(jié)構(gòu)在動(dòng)力荷載作用下（即強(qiáng)迫振動(dòng)）的動(dòng)力反應(yīng)。●動(dòng)內(nèi)力：求最大動(dòng)內(nèi)力，進(jìn)行強(qiáng)度設(shè)計(jì)，使之滿足強(qiáng)度要求?！駝?dòng)位移:求最大動(dòng)位移，進(jìn)行剛度設(shè)計(jì),使之保證振動(dòng)情況能為人體和生產(chǎn)(產(chǎn)品)所允許?！衿渌?速度和加速度):求最大速度、最大加速度，不超過其允許值，使之保證振動(dòng)情況能為人體和生產(chǎn)(產(chǎn)品)所允許。12.1概述1.基本任務(wù)2.研究動(dòng)力反應(yīng)的前提和基礎(chǔ)以上為結(jié)構(gòu)本身的動(dòng)力特性，與荷載無關(guān)。分析結(jié)構(gòu)的自由振動(dòng)，求出：●自振頻率（2p秒內(nèi)振動(dòng)的次數(shù)）●自振周期（振動(dòng)一次所需的時(shí)間）●自振型式（阻尼：使振動(dòng)衰減的因素；阻尼常數(shù)：反映阻尼情況的基本參數(shù)）●阻尼常數(shù)（反映阻尼情況的基本參數(shù)；阻尼：使振動(dòng)衰減的因素）12.1概述

3.土木工程中常見結(jié)構(gòu)振動(dòng)計(jì)算問題

●高層建筑、高聳結(jié)構(gòu)和大跨度橋梁的風(fēng)振分析

●各類工程結(jié)構(gòu)的抗震設(shè)計(jì)

●高速行駛的車輛對(duì)橋梁結(jié)構(gòu)的振動(dòng)影響

●動(dòng)力設(shè)備基礎(chǔ)上的振動(dòng)計(jì)算和減振、隔振設(shè)計(jì)等

●多層廠房中由于動(dòng)力機(jī)器引起的樓面振動(dòng)計(jì)算

4.本課程主要介紹具有線彈性特征的桿件結(jié)構(gòu)，在確定性動(dòng)力荷載作用下的動(dòng)力計(jì)算方法；對(duì)隨機(jī)荷載作用（如地震、風(fēng)振），也將作簡要介紹。

12.1概述12.1.2結(jié)構(gòu)動(dòng)力計(jì)算的特點(diǎn)（三個(gè)方面）1.動(dòng)力荷載的特點(diǎn)

(1)靜力荷載：荷載（大小、方向、作用位置）不隨時(shí)間而變化，或隨時(shí)間極其緩慢地變化（質(zhì)點(diǎn)被近似地視為在常力作用下作勻速運(yùn)動(dòng)，適用于慣性定律，即牛頓第二定律），以致所引起的結(jié)構(gòu)質(zhì)量的加速度（）及其慣性力（）可以忽略不計(jì)。

(2)動(dòng)力荷載(也稱干擾力)：荷載（大小、方向、作用位置）隨時(shí)間明顯變化（質(zhì)點(diǎn)在變力作用下作加速運(yùn)動(dòng)），以致所引起的結(jié)構(gòu)質(zhì)量的加速度（）及其慣性力（）是不可忽略的。12.1概述

2.動(dòng)力反應(yīng)的特點(diǎn)動(dòng)力反應(yīng)與結(jié)構(gòu)本身的動(dòng)力特性有關(guān)。因此，在計(jì)算動(dòng)力反應(yīng)之前，必須先分析結(jié)構(gòu)的自由振動(dòng)，以確定結(jié)構(gòu)的動(dòng)力特性。

3.動(dòng)力計(jì)算方法的特點(diǎn)一般采用動(dòng)靜法（亦稱慣性力法），即動(dòng)力計(jì)算轉(zhuǎn)化為靜力計(jì)算根據(jù)達(dá)朗伯原理(引入附加慣性力，考慮瞬間動(dòng)平衡）12.1概述動(dòng)力計(jì)算所建立的運(yùn)動(dòng)方程為微分方程：1）單自由度體系：一個(gè)變量的二階常微分方程2）多自由度體系：多個(gè)變量的二階常微分方程組3）無限自由度體系：高階偏微分方程12.1概述對(duì)于沖擊、突加等幾種特殊動(dòng)力荷載，則常可采用沖量法求解。12.1.3動(dòng)力荷載的分類

1.周期荷載隨時(shí)間按周期變化的荷載。(1)簡諧荷載是周期荷載中最簡單和最重要的一種。其隨時(shí)間t的變化規(guī)律可用正弦函數(shù)或余弦函數(shù)表示。12.1概述

(2)非簡諧周期荷載

凡有曲柄連桿的機(jī)器（如活塞式空氣壓縮機(jī)、柴油機(jī)、鋸機(jī)等）在勻速運(yùn)轉(zhuǎn)時(shí)都會(huì)產(chǎn)生這種荷載。2.沖擊荷載在很短時(shí)間內(nèi)驟然增減的集度很大的荷載。12.1.3動(dòng)力荷載的分類12.1概述3.突加常量荷載以某一恒值突然施加于結(jié)構(gòu)上并在較長時(shí)間內(nèi)基本保持不變的荷載。12.1.3動(dòng)力荷載的分類12.1概述在將來任一時(shí)刻的數(shù)值無法事先確定的荷載。不能用數(shù)學(xué)式定義，但可采用概率論和數(shù)理統(tǒng)計(jì)的方法，從統(tǒng)計(jì)方面來進(jìn)行定義。地震、脈沖風(fēng)壓和波浪所產(chǎn)生的荷載是其典型例子。

4.隨機(jī)荷載12.1.3動(dòng)力荷載的分類12.1概述12.1.4動(dòng)力計(jì)算中體系的自由度

1.動(dòng)力自由度的定義為了完全確定體系在運(yùn)動(dòng)過程中任一時(shí)刻質(zhì)量位置所必需的獨(dú)立幾何參數(shù)的數(shù)目，稱為體系的動(dòng)力自由度（動(dòng)力分析的基本未知量是質(zhì)點(diǎn)的位移）。

2.體系動(dòng)力自由度的簡化（常用的簡化方法有下列三種）

(1)集中質(zhì)量法

集中質(zhì)量法是從物理的角度提供一個(gè)減少動(dòng)力自由度的簡化方法。該方法把連續(xù)分布的質(zhì)量（根據(jù)靜力等效原則）集中為幾個(gè)質(zhì)點(diǎn)（質(zhì)點(diǎn)，無大小、幾何點(diǎn)，但有質(zhì)量），這樣，就把無限自由度體系，簡化成有限自由度體系。12.1概述a)具有均布質(zhì)量的簡支梁b)無窮多個(gè)集中質(zhì)量

(1)集中質(zhì)量法12.1概述12.1概述

(1)集中質(zhì)量法

(2)廣義坐標(biāo)法廣義坐標(biāo)法是從數(shù)學(xué)的角度提供的一個(gè)減少動(dòng)力自由度的簡化方法。例如，具有分布質(zhì)量的簡支梁的振動(dòng)曲線（位移曲線），可近似地用三角級(jí)數(shù)表示為12.1概述式中，是一組給定的函數(shù)，稱作位移函數(shù)或形狀函數(shù)，與時(shí)間無關(guān)；是一組待定參數(shù)，稱作廣義坐標(biāo)，隨時(shí)間而變化。因此，體系在任一時(shí)刻的位置，是由廣義坐標(biāo)來確定的。注意：這里的形狀函數(shù)只要滿足位移邊界條件，所選的函數(shù)形式可以是任意的連續(xù)函數(shù)。（a）也可寫成更一般的形式式中，是從自動(dòng)滿足位移邊界條件的函數(shù)集合中任意選取的n個(gè)函數(shù)，因此，體系簡化為n個(gè)自由度體系。廣義坐標(biāo)法將應(yīng)用于后面的振型疊加法和能量法。12.1概述(b)

(3)有限單元法有限單元法可看作是廣義坐標(biāo)法的一種特殊應(yīng)用。把體系的離散化和單元的廣義坐標(biāo)法二者結(jié)合起來，就構(gòu)成了有限單元法的概念。有限單元法其具體作法是：第二，取結(jié)點(diǎn)的位移參數(shù)和，即y1，θ1和y2，θ2為廣義坐標(biāo)。第一，將結(jié)構(gòu)離散為有限個(gè)單元。12.1概述第三，分別給出與結(jié)點(diǎn)的位移參數(shù)（均為1時(shí)）相應(yīng)的形狀函數(shù)，即、、和，稱作插值函數(shù)（它們確定了指定結(jié)點(diǎn)位移之間的形狀）。

12.1概述

(3)有限單元法第四，仿照公式(b)，體系的位移曲線可用四個(gè)廣義坐標(biāo)及其相應(yīng)的四個(gè)插值函數(shù)表示為可事先給定，讓其滿足邊界條件。這樣，就把無限自由度體系簡化為四個(gè)自由度（y1，q1，y2，q2）體系。有限單元法綜合了前面集中質(zhì)量法和廣義坐標(biāo)法的某些特點(diǎn)。須強(qiáng)調(diào)的是：動(dòng)力分析中的自由度，一般是變形體體系中質(zhì)量的動(dòng)力自由度。而前面第2章幾何組成分析中的自由度，是不考慮桿件彈性變形的體系的自由度。12.1概述

(3)有限單元法3.動(dòng)力自由度的確定(1)用廣義坐標(biāo)法或有限元法將無限自由度體系簡化為有限自由度體系時(shí)，體系的自由度數(shù)等于廣義坐標(biāo)數(shù)或獨(dú)立結(jié)點(diǎn)位移數(shù)。(2)用集中質(zhì)量法簡化得到的有限自由度體系，在確定體系的自由度數(shù)目時(shí)，應(yīng)注意以下兩點(diǎn)：

1)一般受彎結(jié)構(gòu)的軸向變形忽略不計(jì)。2)動(dòng)力自由度數(shù)不一定等于集中質(zhì)量數(shù)，也與體系是否超靜定和超靜定次數(shù)無關(guān)，但往往與計(jì)算精度有關(guān)。

確定動(dòng)力自由度的方法：一般可根據(jù)定義直接確定；對(duì)于比較復(fù)雜的體系，則可用限制集中質(zhì)量運(yùn)動(dòng)的方法（即附加支桿的方法）來確定。12.1概述

1)單自由度體系12.1概述

2)多自由度體系12.1概述12.1概述

2)多自由度體系12.2單自由度體系的運(yùn)動(dòng)方程描述體系振動(dòng)時(shí)質(zhì)點(diǎn)動(dòng)位移的數(shù)學(xué)表達(dá)式，稱為動(dòng)力體系的運(yùn)動(dòng)方程（亦稱振動(dòng)方程）。單自由度體系的動(dòng)力分析能反映出振動(dòng)的基本特性，是多個(gè)自由度體系分析的基礎(chǔ)。本章只介紹微幅振動(dòng)（線性振動(dòng)）。根據(jù)達(dá)朗伯原理建立運(yùn)動(dòng)方程的方法稱為動(dòng)靜法（亦稱慣性力法）。具體作法有兩種：剛度法和柔度法。剛度法：將力寫成位移的函數(shù)，按平衡條件列出外力（包括假想作用在質(zhì)量上的慣性力和阻尼力）與結(jié)構(gòu)抗力（彈性恢復(fù)力）的動(dòng)力平衡方程（剛度方程），類似于位移法。柔度法：將位移寫成力的函數(shù)，按位移協(xié)調(diào)條件列出位移方程（柔度方程），類似于力法。12.2單自由度體系的運(yùn)動(dòng)方程

1.單自由度體系的振動(dòng)模型12.2.1按平衡條件建立運(yùn)動(dòng)方程——?jiǎng)偠确╩12.2單自由度體系的運(yùn)動(dòng)方程(1)動(dòng)力荷載：(2)彈性恢復(fù)力：(3)阻尼力：(4)慣性力：2.取質(zhì)量m為隔離體，其上有四種力作用

m12.2單自由度體系的運(yùn)動(dòng)方程

3.建立運(yùn)動(dòng)方程根據(jù)達(dá)朗伯原理，由SFx＝0，得代入，即得m(a)

【說明】為了表述簡明，以下各方程和各圖形中的及除FP(t)外各力均省去自變量(t)。12.2單自由度體系的運(yùn)動(dòng)方程【例12-1】試用剛度法建立圖示剛架受動(dòng)力荷載FP(t)作用的運(yùn)動(dòng)方程。解：(1)確定自由度（建模）：結(jié)構(gòu)的質(zhì)量m分布于剛性橫梁，只能產(chǎn)生水平位移，屬單自由度體系。

(2)確定位移參數(shù)：設(shè)剛梁在任一時(shí)刻的位移為y，向右為正。12.2單自由度體系的運(yùn)動(dòng)方程設(shè)剛架的阻尼系數(shù)為c。(3)繪隔離體受力圖：取出隔離體。圖中給出了慣性力、阻尼力和彈性恢復(fù)力。各力均設(shè)沿坐標(biāo)正向?yàn)檎?4)列運(yùn)動(dòng)方程：按動(dòng)靜法列動(dòng)力平衡方程，可得12.2單自由度體系的運(yùn)動(dòng)方程式中代入整理，可得運(yùn)動(dòng)方程mm12.2單自由度體系的運(yùn)動(dòng)方程，，，式中，剛度系數(shù)

k又稱為樓層剛度，系指上下樓面發(fā)生單位相對(duì)位移（D＝1）時(shí)，樓層中各柱剪力之和，如圖所示。m12.2單自由度體系的運(yùn)動(dòng)方程【例12-2】試用剛度法建立圖示靜定梁的運(yùn)動(dòng)方程。解：本例為單自由度體系。取a為坐標(biāo)。在某一時(shí)刻t，體系位移圖和受力如圖所示12.2單自由度體系的運(yùn)動(dòng)方程由SMA＝0，得整理后，得運(yùn)動(dòng)方程12.2單自由度體系的運(yùn)動(dòng)方程12.2單自由度體系的運(yùn)動(dòng)方程【討論】關(guān)于剛度法的三種寫法1)隔離體平衡法。切取質(zhì)量為隔離體，列寫運(yùn)動(dòng)方程。當(dāng)結(jié)構(gòu)作用于質(zhì)量的彈性恢復(fù)力FS容易求得時(shí)，宜用此法（以質(zhì)量為對(duì)象）。2)整體平衡法?？紤]結(jié)構(gòu)整體平衡，列寫運(yùn)動(dòng)方程。當(dāng)無重剛桿上有集中質(zhì)量時(shí)，宜用此法（以結(jié)構(gòu)為對(duì)象）。3)附加約束法。其概念與靜力計(jì)算中的位移法相似。12.2單自由度體系的運(yùn)動(dòng)方程首先，在質(zhì)量上沿動(dòng)力自由度方向添加附加支桿；然后，在質(zhì)量上沿動(dòng)力自由度正向標(biāo)示相關(guān)的力；最后，動(dòng)力平衡位置，體系附加約束中豎向反力F1應(yīng)等于零。12.2單自由度體系的運(yùn)動(dòng)方程12.2.2按位移協(xié)調(diào)條件建立運(yùn)動(dòng)方程――柔度法質(zhì)量m所產(chǎn)生的水平位移，可視為由動(dòng)力荷載FP(t)、慣性力FI和阻尼力FC共同作用在懸臂柱頂端所產(chǎn)生的。根據(jù)疊加原理，得mm12.2單自由度體系的運(yùn)動(dòng)方程式中，d11為柔度系數(shù)。表示在體系的質(zhì)量上沿動(dòng)力自由度方向施加單位力時(shí)，引起該質(zhì)量沿該方向所產(chǎn)生的靜位移。將阻尼力和慣性力代入上式，即得因?yàn)閱巫杂啥润w系中1/d11＝k11（k11和d11互為倒數(shù))，故有(b)12.2單自由度體系的運(yùn)動(dòng)方程【注意】當(dāng)不是直接作用在質(zhì)量及其運(yùn)動(dòng)方向上時(shí)，則式中右邊第三項(xiàng)d11FP(t)應(yīng)改為d1PFP(t)。其中，d1P表示作用時(shí)，引起質(zhì)量沿動(dòng)力自由度方向所產(chǎn)生的位移。相應(yīng)地，前式右邊項(xiàng)FP(t)應(yīng)改為稱等效動(dòng)力荷載。同時(shí)，它與由于FP(t)作用而在質(zhì)點(diǎn)處添加的附加約束上所產(chǎn)生的支反力大小相等。12.2單自由度體系的運(yùn)動(dòng)方程

【例12-3】試用柔度法建立圖示靜定剛架受動(dòng)力荷載作用的運(yùn)動(dòng)方程。解：本題為單自由度體系的振動(dòng)。取質(zhì)量m水平方向的位移y為坐標(biāo)。運(yùn)動(dòng)方程為12.2單自由度體系的運(yùn)動(dòng)方程繪出、圖如圖所示。由圖乘法得得運(yùn)動(dòng)方程圖圖，12.2單自由度體系的運(yùn)動(dòng)方程也可寫作為等效動(dòng)力荷載，即式中，12.2單自由度體系的運(yùn)動(dòng)方程12.2.4建立運(yùn)動(dòng)方程小結(jié)

1)判斷動(dòng)力自由度數(shù)目，標(biāo)出質(zhì)量未知位移正向。

2)沿所設(shè)位移正向加慣性力、阻尼力和彈性恢復(fù)力，并冠以負(fù)號(hào)。

3)根據(jù)是求柔度系數(shù)方便還是求剛度系數(shù)方便，確定是寫柔度方程還是寫剛度方程。

一般情況下，對(duì)于靜定結(jié)構(gòu)，求柔度系數(shù)更為方便；而對(duì)于超靜定結(jié)構(gòu)，則求剛度系數(shù)更為方便。

12.2單自由度體系的運(yùn)動(dòng)方程12.3單自由度體系的自由振動(dòng)

12.3.1自由振動(dòng)1.自由振動(dòng)

2.強(qiáng)迫振動(dòng)體系質(zhì)點(diǎn)在外部干擾力作用下的振動(dòng)，稱為強(qiáng)迫振動(dòng)。由于外界的干擾，質(zhì)點(diǎn)m離開靜力平衡位置，而當(dāng)干擾力消失后，由于彈性恢復(fù)力的作用，質(zhì)點(diǎn)將在靜平衡位置附近作往返運(yùn)動(dòng)。這種在運(yùn)動(dòng)過程中不受干擾力的作用，而由初位移y0或初速度v0（即）或者兩者共同作用下所引起的振動(dòng)，稱為自由振動(dòng)或固有振動(dòng)。12.3.2運(yùn)動(dòng)方程的建立及求解根據(jù)式并令得體系無阻尼自由振動(dòng)方程為可得令，，12.3單自由度體系的自由振動(dòng)

這是一個(gè)二階常系數(shù)齊次線性微分方程。其特征方程為故通解為當(dāng)t＝0時(shí)，（初位移），可求出，而（初速度），可求出。故有12.3單自由度體系的自由振動(dòng)

為將位移方程寫成更簡單的單項(xiàng)形式，引入符號(hào)a和a。使之滿足即(簡諧振動(dòng))

于是可得第一部分：單獨(dú)由y0引起，質(zhì)點(diǎn)按規(guī)律振動(dòng)；第二部分：單獨(dú)由v0引起，質(zhì)點(diǎn)按規(guī)律振動(dòng)。，12.3單自由度體系的自由振動(dòng)

可得由，12.3單自由度體系的自由振動(dòng)

進(jìn)一步說明w、a和a的物理意義（考查一個(gè)模擬的勻速圓周運(yùn)動(dòng)）質(zhì)點(diǎn)做自由振動(dòng)時(shí)其位移隨時(shí)間變化的規(guī)律（圖b），與圖a中質(zhì)點(diǎn)做勻速圓周運(yùn)動(dòng)時(shí)其豎標(biāo)的改變規(guī)律相同。12.3單自由度體系的自由振動(dòng)

w、a和α的物理意義為：w

——自振頻率或圓頻率;a——振幅（自由振動(dòng)時(shí)最大的幅度），ymax；a

——初始相位角，標(biāo)志著t＝0時(shí)質(zhì)點(diǎn)的位置。12.3單自由度體系的自由振動(dòng)

12.3.3自由振動(dòng)中位移、速度、加速度和慣性力的變化規(guī)律由位移，可得由速度，可得由加速度，可得由慣性力，可得(a)(b)(c)(d)12.3單自由度體系的自由振動(dòng)

【注二】慣性力，即永遠(yuǎn)與位移方向一致，在數(shù)值上與位移成比例，其比例系數(shù)?！咀⒁弧坑墒?c)可知，最大加速度的絕對(duì)值等于振幅a與頻率的平方ω2的乘積，將式(a)與式(c)對(duì)照，可見，即加速度與位移成比例，比例系數(shù)為ω2，但方向相反（負(fù)號(hào)），表示加速度永遠(yuǎn)指向平衡位置。(c)(a)12.3.3自由振動(dòng)中位移、速度、加速度和慣性力的變化規(guī)律12.3單自由度體系的自由振動(dòng)

12.3.4自振周期與自振頻率1.自振周期

所以，自振周期因表示體系振動(dòng)一次所需要的時(shí)間，其單位為s（秒）。12.3單自由度體系的自由振動(dòng)

2.工程頻率表示體系每秒振動(dòng)的次數(shù)，其單位為s-1（1/秒）或Hz（赫茲）。一般建筑工程用鋼為7～8次/s，鋼筋混凝土為4次/s，屬低頻；一般機(jī)器為高頻。3.自振頻率

表示體系在2π秒內(nèi)振動(dòng)的次數(shù)，因此也稱圓頻率。其單位為rad/s（弧度數(shù)/秒），也常簡寫為s-1。w是體系固有的非常重要的動(dòng)力特性。在強(qiáng)迫振動(dòng)中，當(dāng)體系的自頻ω與強(qiáng)迫干擾力的擾頻θ很接近時(shí)（0.75≤θ/w

≤1.25區(qū)段），將會(huì)產(chǎn)生共振。為避免共振，就必須使w與θ

遠(yuǎn)離。，12.3單自由度體系的自由振動(dòng)

4.T和w的一些重要性質(zhì)(1)T和w只與結(jié)構(gòu)的m和k11有關(guān)，而與外界的干擾因素?zé)o關(guān)。干擾力的大小只能影響振幅a的大小。(3)結(jié)構(gòu)的T、w是結(jié)構(gòu)動(dòng)力性能的很重要的數(shù)量標(biāo)志。兩個(gè)外表相似的結(jié)構(gòu)，如果T、w相差很大，則動(dòng)力性能相差很大；反之，兩個(gè)外表看來并不相同的結(jié)構(gòu)，如果其T、w相近，則在動(dòng)力荷載作用下其動(dòng)力性能基本一致。工程實(shí)踐中常發(fā)現(xiàn)這樣的現(xiàn)象。(2)質(zhì)量越大，則T越大，w越??；剛度越大，則T越小，w越大。要改變T、w，只有從改變結(jié)構(gòu)的質(zhì)量或剛度（改變截面、改變結(jié)構(gòu)形式和材料）著手。12.3單自由度體系的自由振動(dòng)

5.T、w的計(jì)算公式小結(jié)(1)自振周期式中，g為重力加速度；W=mg為質(zhì)點(diǎn)的重力；Dst=Wk11，表示將重力W=mg施加于振動(dòng)方向所產(chǎn)生的靜位移。12.3單自由度體系的自由振動(dòng)

(2)自振頻率(3)工程頻率12.3單自由度體系的自由振動(dòng)

【例12-4】試求圖示等截面梁的自振周期T和自振頻率w。已知E＝206GPa＝206×109N/m2，I＝245cm4=245×10-8m4。解：采用柔度法。12.3單自由度體系的自由振動(dòng)b)圖a)原結(jié)構(gòu)，【例12-5】試求圖示結(jié)構(gòu)的w。各桿EI＝C。

圖解：12.3單自由度體系的自由振動(dòng)

【例12-6】試求圖示結(jié)構(gòu)的自振頻率w。

圖解：采用柔度法，12.3單自由度體系的自由振動(dòng)

【例12-7】求圖示結(jié)構(gòu)的自振周期T。12.3單自由度體系的自由振動(dòng)

解：采用剛度法。12.3單自由度體系的自由振動(dòng)

【例12-8】試求圖示結(jié)構(gòu)的自振頻率w。已知彈簧剛度系數(shù)，梁的彎曲剛度為EI。解：本題是靜定結(jié)構(gòu)，通常用柔度法計(jì)算比較方便。d

11=d

1+d

212.3單自由度體系的自由振動(dòng)

kk【例12-9】試求圖示結(jié)構(gòu)的自振周期T。解：，12.3單自由度體系的自由振動(dòng)

12.3單自由度體系的自由振動(dòng)

12.3.5多質(zhì)點(diǎn)（包括均質(zhì)剛性桿）的單自由度體系關(guān)于直接平衡法：根據(jù)達(dá)朗伯原理，引入附加慣性力，考慮瞬間動(dòng)平衡，建立運(yùn)動(dòng)方程，并與單自由度體系的微分方程作比較，即可確定w的表達(dá)式。其關(guān)鍵是適當(dāng)選擇質(zhì)量獨(dú)立位移參數(shù)（坐標(biāo)）。這類體系各質(zhì)點(diǎn)仍是平動(dòng)，而質(zhì)量桿一般有轉(zhuǎn)動(dòng)及平動(dòng)，因此，不能再簡單地用單質(zhì)點(diǎn)平動(dòng)公式和表示。w的計(jì)算可以有以下幾種途徑：1）直接平衡法；2）等效質(zhì)量法；3）旋轉(zhuǎn)振動(dòng)公式法等。12.3單自由度體系的自由振動(dòng)

【例12-10】試求圖示梁的自振頻率w。解：在例12-2中，按動(dòng)靜法已建立了運(yùn)動(dòng)方程改寫為12.3單自由度體系的自由振動(dòng)

【例12-11】求圖示梁的自振頻率w。解：由∑MA＝0，列動(dòng)力平衡方程，得，12.3單自由度體系的自由振動(dòng)

解:

【例12-12】求圖示靜定梁的自振頻率w。12.3單自由度體系的自由振動(dòng)

取CB梁段為隔離體圖(d)，由∑MB

＝0，列動(dòng)力平衡方程亦即(a)12.3單自由度體系的自由振動(dòng)

已知，代入式(a)，整理后,得動(dòng)力平衡方程(a)對(duì)比，得12.3單自由度體系的自由振動(dòng)

12.4單自由度體系的強(qiáng)迫振動(dòng)結(jié)構(gòu)在動(dòng)力荷載（也稱干擾力）作用下的振動(dòng)稱為強(qiáng)迫振動(dòng)或受迫振動(dòng)。本節(jié)研究無阻尼的強(qiáng)迫振動(dòng)。在公式中，若不考慮阻尼，則得單自由度體系強(qiáng)迫振動(dòng)的微分方程為或?qū)懗扇鬎P(t)不是直接作用在質(zhì)點(diǎn)上時(shí)，則可將其化為直接作用在質(zhì)點(diǎn)上的等效動(dòng)力荷載FE(t)。12.4單自由度體系的強(qiáng)迫振動(dòng)等效動(dòng)力荷載的幅值式中，F(xiàn)為動(dòng)力荷載的幅值，D1P為F作用下在質(zhì)點(diǎn)振動(dòng)方向產(chǎn)生的位移。12.4單自由度體系的強(qiáng)迫振動(dòng)12.4.1簡諧荷載作用下的動(dòng)力反應(yīng)(本節(jié)重點(diǎn))q

為簡諧荷載的頻率（擾頻）F為荷載的最大值（動(dòng)力荷載幅值）1.求解(1)齊次解12.4單自由度體系的強(qiáng)迫振動(dòng)(2)特解y*設(shè)特解為于是有代入得由此得12.4單自由度體系的強(qiáng)迫振動(dòng)令yst——?jiǎng)恿奢d幅值F產(chǎn)生的位移（最大“靜”位移）

Dst——實(shí)際靜荷載（如自重W）產(chǎn)生的位移（靜位移）12.4單自由度體系的強(qiáng)迫振動(dòng)于是有故特解即12.4單自由度體系的強(qiáng)迫振動(dòng)(3)通解系數(shù)C1和C2由初始條件確定：故通解為

亦即，則得設(shè)，，12.4單自由度體系的強(qiáng)迫振動(dòng)1)前兩項(xiàng)為自由振動(dòng)部分，與初始條件y0和v0有關(guān)?！居懻摗?)第三項(xiàng)與y0和v0無關(guān)，是隨干擾力的出現(xiàn)而伴隨產(chǎn)生的，仍屬自由振動(dòng)（按自頻w振動(dòng)），稱為伴生自由振動(dòng)。3)第四項(xiàng)為純強(qiáng)迫振動(dòng)（無阻尼），按擾頻θ

振動(dòng)。在工程中有實(shí)際意義的是平穩(wěn)階段（第四項(xiàng)）的y，即式中，yd,max稱最大動(dòng)位移（即A），為強(qiáng)迫振動(dòng)的振幅，是控制設(shè)計(jì)的重要依據(jù)。12.4單自由度體系的強(qiáng)迫振動(dòng)則強(qiáng)迫振動(dòng)的振幅

所以有b

的物理意義是：表示動(dòng)位移yd,max的最大值(亦即振幅A)是最大“靜”位移yst的多少倍，故稱動(dòng)力系數(shù)。對(duì)于單自由度體系，當(dāng)在簡諧荷載作用下，且干擾力作用于質(zhì)點(diǎn)上時(shí)，結(jié)構(gòu)中內(nèi)力與質(zhì)點(diǎn)位移成比例。所以動(dòng)力系數(shù)β

既是位移的動(dòng)力系數(shù)，又是內(nèi)力的動(dòng)力系數(shù)。令動(dòng)力系數(shù)12.4單自由度體系的強(qiáng)迫振動(dòng)2.討論（關(guān)于振幅算式的分析）強(qiáng)迫振動(dòng)的振幅其中，動(dòng)力系數(shù)(1)→0，β

→1：這說明機(jī)器轉(zhuǎn)動(dòng)很慢（θ

?ω）時(shí)，干擾力接近于靜力。一般當(dāng)<1/5時(shí)，可當(dāng)作靜力計(jì)算（例如，當(dāng)＝1/5時(shí)，β＝1.04）。位移反應(yīng)譜12.4單自由度體系的強(qiáng)迫振動(dòng)(2)→∞，

β→0：以軸為漸近線。這說明機(jī)器轉(zhuǎn)動(dòng)非?？鞎r(shí)（β

?ω，高頻荷載作用于質(zhì)體），質(zhì)體基本上處于靜止?fàn)顟B(tài)，即相當(dāng)于沒有干擾力作用（自重除外）。(3)0<<1，β

為正，且β

>1，又β

隨的增大而增大。與同號(hào)，即質(zhì)點(diǎn)位移與干擾力的方向每時(shí)每刻都相同（同相位）。位移反應(yīng)譜12.4單自由度體系的強(qiáng)迫振動(dòng)(4)>1，b

為負(fù)，其絕對(duì)值隨的增大而減小。y與異號(hào)，即質(zhì)點(diǎn)位移與干擾力的方向相反（相位相差p）。位移反應(yīng)譜12.4單自由度體系的強(qiáng)迫振動(dòng)(5)→1，β→∞（無阻尼）：→∞（振幅隨時(shí)間而逐漸增大），體系發(fā)生共振。此時(shí)有：即慣性力與彈性力平衡，而沒有什么力去與實(shí)際存在的外力FP(t)平衡，因此無論振幅多大，再維持動(dòng)力平衡均不可能。當(dāng)θ=ω時(shí)位移反應(yīng)譜12.4單自由度體系的強(qiáng)迫振動(dòng)防止共振的措施：一是調(diào)整機(jī)器的轉(zhuǎn)速θ；二是改變體系的自振頻率w（改變w的思路，不外就是改變k11，即改變截面形式、結(jié)構(gòu)形式，或是改變m）。但“共振”也是可以利用的，如利用q

＝w時(shí)，結(jié)構(gòu)振幅突出大的這一特點(diǎn)，不斷改變機(jī)器（激振器）轉(zhuǎn)速θ，可以測(cè)定結(jié)構(gòu)的w。位移反應(yīng)譜12.4單自由度體系的強(qiáng)迫振動(dòng)

3.計(jì)算步驟（單自由度體系在簡諧荷載作用下的強(qiáng)迫振動(dòng)）

(1)求自振頻率

(2)求干擾力頻率式中，n為電動(dòng)機(jī)轉(zhuǎn)速n（r/min）。（1/s）

(3)求動(dòng)力系數(shù)（注意正負(fù)號(hào)）12.4單自由度體系的強(qiáng)迫振動(dòng)(4)求動(dòng)位移幅值?動(dòng)（即A）

1)先求最大“靜”位移3.計(jì)算步驟

2)再求動(dòng)位移幅值(5)求最大位移Dmax=D動(dòng)+D靜=yd,max+Dst=|A|+Dst(6)求最大內(nèi)力

Mmax＝M動(dòng)＋M靜＝Md,max＋Mst12.4單自由度體系的強(qiáng)迫振動(dòng)【方法一】動(dòng)力系數(shù)法（僅當(dāng)FP(t)直接作用在質(zhì)點(diǎn)上時(shí)）：將|β

|F作為靜力作用在體系上，按靜力法計(jì)算(圖a)。a)動(dòng)力系數(shù)法12.4單自由度體系的強(qiáng)迫振動(dòng)【方法二】幅值法：由達(dá)朗伯原理，把位移達(dá)到最大值時(shí)，所有力的幅值加上去。注意F的施加方向，即1)當(dāng)β

為正時(shí)，F(xiàn)沿質(zhì)點(diǎn)位移方向一致施加(圖c)。2)當(dāng)β

為負(fù)時(shí)，F(xiàn)沿質(zhì)點(diǎn)位移方向反向施加(圖b)。

b)幅值法(β為負(fù))c)幅值法(β為正)FF12.4單自由度體系的強(qiáng)迫振動(dòng)【例12-13】對(duì)于圖示體系，已知下列各值：m=123kg，F(xiàn)=49N(離心力)，n=1200r/min(發(fā)電機(jī)轉(zhuǎn)速)，E=2.06×1011N/m2，I=78cm4。求梁中最大動(dòng)位移A(?動(dòng))和梁中最大動(dòng)內(nèi)力Md,max(M動(dòng))。解:(1)求自振頻率ω(2)求干擾頻率θFsinqt12.4單自由度體系的強(qiáng)迫振動(dòng)(3)求動(dòng)力系數(shù)b(4)求最大動(dòng)位移A負(fù)號(hào)表示最大動(dòng)位移與FP(t)方向相反。12.4單自由度體系的強(qiáng)迫振動(dòng)(5)求最大動(dòng)內(nèi)力Md,max：采用動(dòng)力系數(shù)法，在B點(diǎn)施加|b|F，繪彎矩圖，如圖所示，圖中Md,max＝16.33N·m。16.33N.mb)Md,max圖12.4單自由度體系的強(qiáng)迫振動(dòng)【例12-14】對(duì)于圖示體系，已知：梁上的機(jī)器總重W＝30kN，機(jī)器轉(zhuǎn)速n＝350r/min，離心力幅值F=5kN，忽略梁的自重，EI＝2.0×104kN·m，試作動(dòng)力彎矩幅值圖Md,max（即M動(dòng)）和總彎矩M圖。解:結(jié)構(gòu)為單自由度體系，采用動(dòng)力系數(shù)法求解。(1)求柔度系數(shù)δ11b)圖(m)12.4單自由度體系的強(qiáng)迫振動(dòng)(2)求自振頻率w(3)求干擾頻率q(4)求動(dòng)力系數(shù)b12.4單自由度體系的強(qiáng)迫振動(dòng)(5)作動(dòng)力彎矩幅值Md,max圖

將作用于梁上D點(diǎn)，作Md,max圖(6)作總彎矩圖M將作用于梁上D點(diǎn)，作M圖。c)動(dòng)力彎矩幅值Md,max圖（kN·m)d)總彎矩M圖（kN·m）DD12.4單自由度體系的強(qiáng)迫振動(dòng)【例12-15】圖示結(jié)構(gòu)，在柱頂有馬達(dá)，試求馬達(dá)轉(zhuǎn)動(dòng)時(shí)的最大水平位移和柱端彎矩。已知馬達(dá)和結(jié)構(gòu)重量集中于柱頂。W＝20kN，馬達(dá)水平離心力的幅值F＝250N，馬達(dá)轉(zhuǎn)速n=550r/min，柱的線剛度i=EI/h=5.88×106N·m。解：(1)求自振頻率w12.4單自由度體系的強(qiáng)迫振動(dòng)(2)求干擾頻率θ(3)求動(dòng)力系數(shù)β(4)求沿水平方向的?max負(fù)號(hào)表示Dmax與FP(t)的方向相反。12.4單自由度體系的強(qiáng)迫振動(dòng)(5)求柱端彎矩（分別采用兩種解法）【解法一】讓二柱頂C和D點(diǎn)均水平移動(dòng)?＝A，按兩端固定端產(chǎn)生順時(shí)針方向側(cè)移?＝A，繪柱彎矩圖，如圖所示，即各柱端彎矩均為【解法二】采用“動(dòng)力系數(shù)法”，將向右施加于左柱頂C點(diǎn)，再用“剪力分配法”計(jì)算，其彎矩圖如圖c所示，即各柱端彎矩均為c)解法二M圖b)解法一M圖12.4單自由度體系的強(qiáng)迫振動(dòng)

【例12-16】干擾力作用于等截面懸臂桿的中點(diǎn)，求質(zhì)點(diǎn)穩(wěn)態(tài)振幅。解：【解法一】附加支桿，將該支反力反向（即等效干擾力，亦稱等效動(dòng)力荷載）作用于質(zhì)點(diǎn)。(1)計(jì)算附加支桿內(nèi)的反力幅值（力法）12.4單自由度體系的強(qiáng)迫振動(dòng)

(2)計(jì)算作用于質(zhì)點(diǎn)上的情況12.4單自由度體系的強(qiáng)迫振動(dòng)【解法二】直接建立運(yùn)動(dòng)方程求解。施加慣性力，列柔度方程。位移y相當(dāng)于作用于質(zhì)點(diǎn)的等效力產(chǎn)生的位移。在質(zhì)量所在的B點(diǎn)12.4單自由度體系的強(qiáng)迫振動(dòng)或具體計(jì)算如下：，圖圖12.4單自由度體系的強(qiáng)迫振動(dòng)【解法三】利用幅值方程求解。在簡諧振動(dòng)中，慣性力與位移變化規(guī)律相同，即同時(shí)達(dá)到最大值，可列幅值方程.由此，得

不難化為12.4單自由度體系的強(qiáng)迫振動(dòng)

【小結(jié)】情況1：干擾力直接作用在質(zhì)點(diǎn)上12.4單自由度體系的強(qiáng)迫振動(dòng)情況2：干擾力不作用在質(zhì)點(diǎn)上12.4單自由度體系的強(qiáng)迫振動(dòng)解：【例12-17】圖示跨中帶有一質(zhì)量m的無重簡支梁，動(dòng)力荷載作用在距梁端l/4處，若，試求在荷載作用下，質(zhì)量m的最大動(dòng)位移A。FE=0.6875F圖圖，a)b)c)d)12.4單自由度體系的強(qiáng)迫振動(dòng)FE=0.6875F或12.4單自由度體系的強(qiáng)迫振動(dòng)【例12-18】試列出圖示結(jié)構(gòu)體系的運(yùn)動(dòng)方程，并繪出結(jié)構(gòu)彎矩幅值圖。已知：。解:(1)質(zhì)點(diǎn)的等效力

，，，12.4單自由度體系的強(qiáng)迫振動(dòng)

(3)按幅值法作彎矩幅值圖，如圖所示。

(2)求慣性力幅值12.4單自由度體系的強(qiáng)迫振動(dòng)【例12-19】圖示簡支梁跨中有一集中質(zhì)量m，EI為常量，跨度為l，不計(jì)梁的質(zhì)量。梁右端作用干擾力偶。試作彎矩幅值圖并求梁右端角位移的幅值。設(shè)靜力平衡時(shí)梁軸線為水平直線。已知。解：(1)求最大動(dòng)位移

，，，12.4單自由度體系的強(qiáng)迫振動(dòng)

(2)求慣性力幅值

(3)求內(nèi)力幅值

(4)求12.4單自由度體系的強(qiáng)迫振動(dòng)12.4.2一般動(dòng)力荷載（任意干擾力作用）1.瞬時(shí)沖量的動(dòng)力反應(yīng)設(shè)體系在t＝0時(shí)處于靜止?fàn)顟B(tài)。在質(zhì)點(diǎn)上施加瞬時(shí)沖量。這將使體系產(chǎn)生初速度，但初位移仍為0，即y0＝0（可以證明，y0系二階微量，可略去不計(jì)）。將y0和v0代入上式就是t=0時(shí)作用瞬時(shí)沖量S所引起的動(dòng)力反應(yīng)。即得12.4單自由度體系的強(qiáng)迫振動(dòng)如果瞬時(shí)沖量S從t=t開始作用，則式中的位移反應(yīng)時(shí)間t，應(yīng)改成(t-t

)，即上式應(yīng)改為12.4單自由度體系的強(qiáng)迫振動(dòng)2.任意動(dòng)力荷載的動(dòng)力反應(yīng)（總效應(yīng)）由式得到（對(duì)于t>t）

整個(gè)加載過程可看作一系列瞬時(shí)沖量所組成。在時(shí)，作用，在微分段內(nèi)產(chǎn)生的微分沖量為12.4單自由度體系的強(qiáng)迫振動(dòng)總反應(yīng)為此式稱為杜哈梅(J.M.C.Duhamel)積分（卷積）。這是初始處于靜止?fàn)顟B(tài)的單自由度體系在任意動(dòng)力荷載作用下的位移公式。如果(在O點(diǎn))初始位移y0和初始速度v0不為0，則總位移應(yīng)為(自由振動(dòng))

(伴生自由振動(dòng)＋純強(qiáng)迫振動(dòng)) 12.4單自由度體系的強(qiáng)迫振動(dòng)【說明1】這里為什么用dt而不用dt？(自由振動(dòng))

(伴生自由振動(dòng)＋純強(qiáng)迫振動(dòng))

我們是在考察加在不同時(shí)刻t的一系列瞬時(shí)沖量對(duì)同一時(shí)刻t的位移的影響。這里位移發(fā)生的時(shí)刻t被暫時(shí)地固定起來（是指定的常數(shù)），而瞬時(shí)沖量施加的時(shí)刻t表示時(shí)間的流動(dòng)坐標(biāo)，是變量。因此，變量的微分為dt

，而非dt。12.4單自由度體系的強(qiáng)迫振動(dòng)(自由振動(dòng))

(伴生自由振動(dòng)＋純強(qiáng)迫振動(dòng))

【說明2】在杜哈梅積分中，能否把伴生自由振動(dòng)分離出來？對(duì)于簡諧荷載，可以證明，能將杜哈梅積分分解為以下兩項(xiàng)之和，即(伴生自由振動(dòng))

(純強(qiáng)迫振動(dòng))其中12.4單自由度體系的強(qiáng)迫振動(dòng)12.4.3幾種特殊荷載作用下的動(dòng)力反應(yīng)

1.突加長期荷載當(dāng)t>0時(shí)12.4單自由度體系的強(qiáng)迫振動(dòng)所以仍系周期運(yùn)動(dòng)，但不是簡諧運(yùn)動(dòng)。t>0時(shí)，質(zhì)點(diǎn)圍繞其靜平衡位置（新的基線）作簡諧運(yùn)動(dòng)，即突加荷載所引起的最大動(dòng)位移A比相應(yīng)的最大靜位移增