第7章 逐步法——對一般動力荷載的反應(yīng)

1、高等結(jié)構(gòu)動力學(xué)第七章第七章高等結(jié)構(gòu)動力學(xué)對一般動力荷載的反應(yīng)逐步法高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué) 分析承受任意動力荷載的線性結(jié)構(gòu)，分析承受任意動力荷載的線性結(jié)構(gòu)，Duhamel積積分或頻域分析，提供了最方便的解法。分或頻域分析，提供了最方便的解法。 這兩種方法的推導(dǎo)過程中都使用了疊加原理，只這兩種方法的推導(dǎo)過程中都使用了疊加原理，只能適用于線性體系，即反應(yīng)過程中體系的特性保持不能適用于線性體系，即反應(yīng)過程中體系的特性保持不變。變。 另一方面，有許多種重要的結(jié)構(gòu)動力學(xué)問題，體另一方面，有許多種重要的結(jié)構(gòu)動力學(xué)問題，體系不能視作線性的。如：足以引起嚴(yán)重破壞的地震運系不能視作線性的。如：足以引起嚴(yán)重破

2、壞的地震運動下的建筑物反應(yīng)等等。因此，還需要發(fā)展適用于非動下的建筑物反應(yīng)等等。因此，還需要發(fā)展適用于非線性體系的其它分析方法。線性體系的其它分析方法。高等結(jié)構(gòu)動力學(xué)動力反應(yīng)分析的方法動力反應(yīng)分析的方法1 、疊加法、疊加法線性體系，即反應(yīng)過程中體系的特性保線性體系，即反應(yīng)過程中體系的特性保持不變；持不變；2、逐步法、逐步法體系不能視作線性的，要發(fā)展適用于非體系不能視作線性的，要發(fā)展適用于非線性體系的方法。線性體系的方法。高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué)0( )pp(7-1)(7-1)0mvcvkvp(7-2)(7-2)( )( )( )hpvvv(7-3)(7-3)( )

3、exp()(cossin)hDDvAB 021( )()pacvpkk(7-4)(7-4)0123( )exp()cos exp()sinDDvAAAA (7-5)(7-5)高等結(jié)構(gòu)動力學(xué)其中其中00232vA12A200AvA30221()DAvA同樣地，可獲得時間步長期間的速度為同樣地，可獲得時間步長期間的速度為13223( )() exp() cos (A+A )exp()sinDDDDvAAA (7-6)(7-6)高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué)圖圖E7-2 分段精確計算的反應(yīng)分段精確計算的反應(yīng)高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué)0000mvcvkvp00001()vpcvkvm01

4、101 21 2 vvvvvvhh；1 21 2010012211()()vvvvvvvhhh（7-77-7）（7-87-8）高等結(jié)構(gòu)動力學(xué)010121(2)vvvvh（7-97-9）21010002()hvvvpcvkvm2100001()2hvpcvkvvvm（7-107-10）1102vvvh1102vvhv高等結(jié)構(gòu)動力學(xué)2100000()2hvvhvpcvkvm（7-117-11）10011()2vvvvh10102()vvvvh（7-127-12）高等結(jié)構(gòu)動力學(xué) 另一類一般性的逐步進行動力反應(yīng)分析的數(shù)值方法是，對每一時間步，從初始到最終條件應(yīng)用積分前進一步。這個基本概念可用如下式子表

5、示：100100( )( )hhvvvdvvvd （7-13a7-13a）（7-13b7-13b）高等結(jié)構(gòu)動力學(xué) 最終速度和位移依據(jù)這些值的初始值加一個積分表達式。速度的變化依賴于加速度歷程的積分，而位移的變化依賴于相應(yīng)的速度積分。 為了進行這類分析，首先需要假設(shè)在時間步的持續(xù)時間內(nèi)加速度是如何變化的；加速度的假設(shè)也控制了速度的變化，因而可以由這一步向前獲得下一時間步。高等結(jié)構(gòu)動力學(xué)Euler-Gauss方法 最簡單是假設(shè)加速度在時間步持續(xù)時間內(nèi)為常數(shù)，結(jié)果是在持續(xù)時間內(nèi)速度為線性，位移為二次曲線 著名的Euler-Gauss方法。高等結(jié)構(gòu)動力學(xué) 列式特性的樣式，假設(shè)常量加速度是由初值及步長持

6、續(xù)時間內(nèi)所獲得的最終加速度的平均。在此圖中也顯示了速度和位移的表達式，它們是對此加速度在這步持續(xù)時間內(nèi)任意時刻 由逐次積分所獲得的，把 代入這些表達式而獲得最終速度和位移。h高等結(jié)構(gòu)動力學(xué)圖7-3 基于常平均加速度的運動h0t1t210001()4hvvv hvv1001()2hvvvv1v 0v 0v 0v011()2avvvv001( )()2vvvv20001( )()4vvvvvh 加速度（常數(shù)） 速度（線性） 位移（二次的）高等結(jié)構(gòu)動力學(xué)為了對任意步開始這種分析，首先需要計算初始加速度0v 時刻0tt式（7-7）所示的動力平衡表達式獲得。另外，最終加速度1v 需要應(yīng)用隱式列式，它的值

7、可以由迭代獲得。對，這可以解1v 它的值可以由迭代獲得。對開始時用任意假設(shè)的值，再用圖7-3所列式（a）和（b）得到1v 和1v的值。然后，用與式（7-7）相當(dāng)?shù)谋磉_式從動力平衡方程計算1t時刻1v 值的一個改進，由此再導(dǎo)得速度和1v 和1v 和位移1v 1v的改進值，最后，迭代收斂于這時間步最終加速度的一個固定的值，這個過程可以前進一步到下一時間步。高等結(jié)構(gòu)動力學(xué)迭代的列式11111()vpcvkvm1v 開始時用任意假設(shè)的值 ， 再用圖7-3所列式（a）210001()4hvvv hvv1001()2hvvvv高等結(jié)構(gòu)動力學(xué) 常平均加速度法的主要優(yōu)點：是無條件穩(wěn)定的。也就是說，從一步到下一

8、步不管時間步長選得如何長，誤差不會放大。因此時間步長的選擇只需要考慮所定義動力激勵和結(jié)構(gòu)的振動反映特性。高等結(jié)構(gòu)動力學(xué)Newmark法 一種更一般的逐步列式是由Newmark提出的，前面的方法可以作為它的特殊情況。但是也可以在其他一些形式下應(yīng)用。在Newmark列式中，對最終速度和位移的基本積分式（7-13）如下所示：高等結(jié)構(gòu)動力學(xué)1001(1)vvhvhv22100011()2vvhvh vh v（7-14a7-14a）（7-14b7-14b）由式（7-14a）顯然可見，系數(shù) 提供了在初始和最終加速度改變影響之間的線性變化的權(quán)重，類似地，系數(shù)提供了在這些初始和最終加速度對位移改變貢獻的權(quán)重。

9、高等結(jié)構(gòu)動力學(xué) 從該列式性能的研究發(fā)現(xiàn)，系數(shù)控制了由這個逐步法導(dǎo)致的人工阻尼量；如果=1/2，方法是無人工阻尼的，因此這個值被推薦用于標(biāo)準(zhǔn)的單自由度分析。 在式（7-14a）和式（7-14b）中令系數(shù)=1/2和=1/4，此時可以看到， Newmark列式直接退化為圖7-3所示最終速度和位移的表達式。因此， Newmark=1/4法也可以歸諸于常平均加速度法。高等結(jié)構(gòu)動力學(xué) 另一方面，如果取作1/6（用=1/2），最終速度和位移的表達式成為1001()2hvvvv221000136hhvvv hvv（7-15a7-15a）（7-15b7-15b）高等結(jié)構(gòu)動力學(xué) 這些結(jié)果也可以如圖7-4所示，由假

10、設(shè)在時間步持續(xù)期間加速度在 和 的初始到最后值之間線性變化來得到；因此=1/6的Newmark法也稱為線加速度法。像常平均加速度法一樣，此法在實際中也是廣為應(yīng)用的。但是與=1/4方法對比，線加速度法僅是條件穩(wěn)定的?？墒牵c二階中心差分法一樣，在單自由度體系分析中這個限制并不重要，因為要獲得動力荷載和反應(yīng)的滿意表示，必須取比這一限制更短的時間步長。1v 0v 高等結(jié)構(gòu)動力學(xué)變換到顯式公式 法的隱式列式是不方便應(yīng)用的，因為每一時間步內(nèi)為了確定此步終點加速度需要進行迭代。因此，通常被修改為顯式形式，目的是最終加速度用其他反應(yīng)量表示選擇一個基本未知量（位移較好）。高等結(jié)構(gòu)動力學(xué) 再代入圖7-3式（a）

11、中，獲得最終速度表達式為11002()vvvvh（7-16b7-16b）11000244()vvvvvhh（7-16a7-16a）根據(jù)圖7-3式（b）對最終加速度求解可得高等結(jié)構(gòu)動力學(xué)在t1時刻寫出動力平衡方程1111mvcvkvp并將式（7-16a）和（7-16b）代入上式，則可導(dǎo)得僅含時間步終點未知位移v1的表達式。經(jīng)適當(dāng)歸并同類項，此式可寫為c11ck vp（7-177-17）這是一個靜力平衡方程的形式，它包含等效剛度c224cmkkhh（7-17a7-17a）高等結(jié)構(gòu)動力學(xué)和等效荷載0010001c2244()()vvppcvmvvhhh（7-17b7-17b）在式（7-17）里下標(biāo)c

12、用以標(biāo)記常平均加速度法。 1v 1v 高等結(jié)構(gòu)動力學(xué)11111()vpcvkvm而不是從式（7-16a）求，因而保留了平衡條件。使用這個顯式公式，時間步終點位移v1可直接由式（7-17）計算，所要用到的僅是時間步開始時的數(shù)據(jù)。然后，此時刻的速度 可用式（7-16b）計算，最后，此時間步終點的加速度 由求解該時刻的動力平衡方程而得 高等結(jié)構(gòu)動力學(xué) 采用同樣的方法，使用圖7-4中的式（a）和（b），也可以類似地將線性加速度法轉(zhuǎn)換為顯式形式，這些列式的位移差別就是等效剛度，等效荷載及最終速度的表達式不同。對線加速度分析來說，等效靜力平衡方程為 11ddkvp（7-187-18）高等結(jié)構(gòu)動力學(xué)圖7-4

13、 基于線性變化加速度的運動0t1t221000136hhvvv hvv1001()2hvvvv1v 0v 0v 0v100( )()vvvvh21000( )()vvvvvhh2310000( )()26vvvvvvh 加速度（線性） 速度（二次的） 位移（三次的）h高等結(jié)構(gòu)動力學(xué)001000012636(2)(2)2dvvhppmvvcvvhhh（7-18b7-18b）當(dāng)位移v 1由式（7-18）計算時，同時刻的速度可由如下表達式給出相當(dāng)于式（7-16b）:110003()22hvvvvvh（7-18c7-18c）其中下標(biāo)d表示線加速度法。等效剛度和荷載分別為 236dcmkkhh（7-18

14、a7-18a）高等結(jié)構(gòu)動力學(xué) 線加速度法僅僅是條件穩(wěn)定，但如前面所述，對于單自由度體系分析，這一點并不重要。另一方面，假設(shè)每個步長持續(xù)時間內(nèi)加速度線性變化，要比連續(xù)用常加速度法能獲得真實特性的更好近似。 實際上，數(shù)值實驗結(jié)果也證明了線加速度法結(jié)果比用常加速度步所得結(jié)果優(yōu)越?；诖死碛桑瑢巫杂啥润w系的分析推薦使用線加速度（=1/6）法。高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué)圖圖 7-5 非線性動力體系的定義：非線性動力體系的定義：(a)基本單自由度結(jié)構(gòu)；基本單自由度結(jié)構(gòu)；(b)力的平衡；力的平衡；(c)非線性阻尼；非線性阻尼；(d)非線性剛度；非線性剛度；(e)作用荷載作用荷載(a)(b)

15、(c)(d)(e)高等結(jié)構(gòu)動力學(xué) 考慮的結(jié)構(gòu)為圖考慮的結(jié)構(gòu)為圖7-5(a)所示的單自由度體系，體系的特性所示的單自由度體系，體系的特性m，k，c和和p(t)可以理解為可以理解為2-5節(jié)所討論的廣義量，而并不只局限于圖節(jié)所討論的廣義量，而并不只局限于圖面上所示的簡單情況。作用于體系質(zhì)量面上所示的簡單情況。作用于體系質(zhì)量m上的力如圖上的力如圖7-5(b)所示，所示，彈簧力和阻尼力的一般非線性性質(zhì)分別繪于圖彈簧力和阻尼力的一般非線性性質(zhì)分別繪于圖7-5(c)和和(d)中，任中，任意作用荷載則繪于圖意作用荷載則繪于圖7-5(e)中。中。 在任一瞬間在任一瞬間t，作用于質(zhì)量，作用于質(zhì)量m上的力系的平衡要

16、求：上的力系的平衡要求：)()()()(tptftftfSDI（7-19a7-19a） 而在暫短時間而在暫短時間t以后，平衡方程將為以后，平衡方程將為)()()()(ttpttfttfttfSDI（7-19b7-19b） 從方程從方程(7-19b)減去方程減去方程(7-19a)則可得時間間隔則可得時間間隔t的運動方程的增量形式：的運動方程的增量形式：（7-207-20）)()()()(tptftftfSDI高等結(jié)構(gòu)動力學(xué) 這個方程中的增量力可表示如下：這個方程中的增量力可表示如下：（7-21a7-21a） 其中，質(zhì)量其中，質(zhì)量m不言而喻被假設(shè)為常量，而不言而喻被假設(shè)為常量，而c(t)和和k(t

17、)項則表示項則表示與時間間隔內(nèi)速度和位移相應(yīng)的阻尼和剛度特性，分別如圖與時間間隔內(nèi)速度和位移相應(yīng)的阻尼和剛度特性，分別如圖7-5(c)和和(d)所示。實際上，因為在時間增量末端的速度和位栘將所示。實際上，因為在時間增量末端的速度和位栘將依賴于這些特性，故所示割線斜率只能用依賴于這些特性，故所示割線斜率只能用迭迭代法進行計算。代法進行計算。因因此此，通常用時間間隔起點所定義的切線斜率來代替。，通常用時間間隔起點所定義的切線斜率來代替。（7-227-22）)()()()()()()()()()()()()()()()()(tpttptptvtktfttftftvtctfttftftvmtfttft

18、fSSSDDDIII tDvddftc)(tSdvdftk)(（7-21b7-21b）（7-21c7-21c）（7-21d7-21d）高等結(jié)構(gòu)動力學(xué) 將方程將方程(7-21)代入方程代入方程(7-20)，可導(dǎo)得時間，可導(dǎo)得時間t的增量平衡方程的增量平衡方程的最終形式為：的最終形式為：（7-237-23） 在這種類型的分析中，所討論的材料特性可以包括任何的非在這種類型的分析中，所討論的材料特性可以包括任何的非線性形式。因此，沒有必要規(guī)定彈簧力線性形式。因此，沒有必要規(guī)定彈簧力fs像非線性彈性材那樣像非線性彈性材那樣僅僅依賴于位移。同樣也可以很好地說明非線性滯變材料，在僅僅依賴于位移。同樣也可以很

19、好地說明非線性滯變材料，在這種材料中，力依賴于變形的過去時程以及位移的當(dāng)前值，唯這種材料中，力依賴于變形的過去時程以及位移的當(dāng)前值，唯一的要求是剛度特性必須完全由變形的過去和目前狀態(tài)所確定。一的要求是剛度特性必須完全由變形的過去和目前狀態(tài)所確定。此外，隱含的質(zhì)量不變的假定顯然是可以改動的：質(zhì)量也可以此外，隱含的質(zhì)量不變的假定顯然是可以改動的：質(zhì)量也可以表示為隨時間變化的量。表示為隨時間變化的量。)()()()()()(tptvtktvtctvm 高等結(jié)構(gòu)動力學(xué) 有很多方法可以用來進行方程有很多方法可以用來進行方程(7-23)的數(shù)值積分。的數(shù)值積分。 引入假定：在每個時間增量內(nèi)加速度線性變化，而

20、引入假定：在每個時間增量內(nèi)加速度線性變化，而且體系的特性在這個間隔內(nèi)保持為常量。且體系的特性在這個間隔內(nèi)保持為常量。 質(zhì)量在時間間隔內(nèi)的運動用圖形的形式繪于圖質(zhì)量在時間間隔內(nèi)的運動用圖形的形式繪于圖7-6中，與所中，與所假定的加速度為線性變化的方程一起，圖中還分別繪出了相應(yīng)假定的加速度為線性變化的方程一起，圖中還分別繪出了相應(yīng)的速度二次變化、位移三次變化的圖形。計算在間隔終點的速度二次變化、位移三次變化的圖形。計算在間隔終點(t)時后的值，導(dǎo)得速度和位移的增量方程如下：時后的值，導(dǎo)得速度和位移的增量方程如下：（7-6a7-6a* *）2)()(ttvtvtv 6)(2)()(22ttvttvt

21、vtv （7-6b7-6b* *）高等結(jié)構(gòu)動力學(xué)圖圖 7-6 基于線性變化加速基于線性變化加速 度的增量運動度的增量運動高等結(jié)構(gòu)動力學(xué)高等結(jié)構(gòu)動力學(xué) ()2()23()( )( )( )2( )26iiiiiiiiitttttttttv tvvtv tvvvtv tvvvvt高等結(jié)構(gòu)動力學(xué)222( )( )( )2( )( )( )2( )( )( )26v ttv tv ttttv tv ttv tttv tv tv tv t t 取高等結(jié)構(gòu)動力學(xué) 位移增量作為分析的基本變量，解式位移增量作為分析的基本變量，解式(7-6a*)可得加速度增量，并將這個可得加速度增量，并將這個表達式代入方程表達

22、式代入方程(7-6b*)，得：，得：（7-7a7-7a* *）（7-7b7-7b* *）)(3)(6)(6)(2tvtvttvttv )(2)(3)(3)(tvttvtvttv 獲得了用位移增量表達的速度和加速度增量。獲得了用位移增量表達的速度和加速度增量。把方程把方程(7-7*)代入方程代入方程(7-23),可導(dǎo)得下列的運動方程：可導(dǎo)得下列的運動方程：2663( )( )3 ( )( )( )3 ( )( )( )( )( )2tmv tv tv tc tv tv tv tk tv tp tttt （7-87-8* *）高等結(jié)構(gòu)動力學(xué)最終把所有含已知初始條件的各項移到右邊，給出最終把所有含已

23、知初始條件的各項移到右邊，給出（7-9a7-9a* *）（7-9b7-9b* *） 方程方程(7-8*)相當(dāng)于靜力增量平衡關(guān)系，即可將荷載增量除以相當(dāng)于靜力增量平衡關(guān)系，即可將荷載增量除以剛度求得位移增量。在等效荷載和剛度項中包含慣性和阻尼的剛度求得位移增量。在等效荷載和剛度項中包含慣性和阻尼的項反映了動力特性。解方程（項反映了動力特性。解方程（7-8*），得出位移增量后，將此值），得出位移增量后，將此值代入方程代入方程(7-7b*)即可獲得速度增量。即可獲得速度增量。 下一時段的初始條件由該時段起點速度和位移值加上這些增下一時段的初始條件由該時段起點速度和位移值加上這些增量值得到。量值得到。

24、)()()(tptvtk（7-87-8* *）其中：其中：236( )( )( )k tk tc tmtt)(2)(3)()(3)(6)()(tvttvtctvtvtmtptp 高等結(jié)構(gòu)動力學(xué) 分析包含兩個重要假定：分析包含兩個重要假定： (1)加速度為線性變化；加速度為線性變化； (2)質(zhì)量、阻尼和剛度特性在時間步長內(nèi)保持常量。質(zhì)量、阻尼和剛度特性在時間步長內(nèi)保持常量。 雖然時間步長很短時誤差甚小，但這兩個假定都不是完全雖然時間步長很短時誤差甚小，但這兩個假定都不是完全正確的。因此，誤差在增量平衡關(guān)系中出現(xiàn)，并積累。為了避正確的。因此，誤差在增量平衡關(guān)系中出現(xiàn)，并積累。為了避免誤差的積累，在

25、分析的每一步中利用總的平衡條件消除誤差。免誤差的積累，在分析的每一步中利用總的平衡條件消除誤差。 只要從總外荷載中減去總阻尼力和彈性力以表示時間步長只要從總外荷載中減去總阻尼力和彈性力以表示時間步長起點的加速度就行了。起點的加速度就行了。)()()(1)(tftftpmtvSD 高等結(jié)構(gòu)動力學(xué) 對任一給定的時間增量，按如下程序進行運算：對任一給定的時間增量，按如下程序進行運算： (1)初始速度和位移值是已知的。它們或是前一增量的終點初始速度和位移值是已知的。它們或是前一增量的終點值或是問題的初始條件值；值或是問題的初始條件值； (2)利用這些值及結(jié)構(gòu)特定的非線性特性，可找出時間間隔利用這些值及

26、結(jié)構(gòu)特定的非線性特性，可找出時間間隔內(nèi)的阻尼內(nèi)的阻尼c(t)、剛度、剛度k(t)以及阻尼力以及阻尼力fD(t)和彈性力和彈性力fS(t)的當(dāng)前值，的當(dāng)前值，例如圖例如圖7-5(c)和和7-5(d)中所示；中所示； (3)初始加速度由下式給出：初始加速度由下式給出：（7-257-25）)()()(1)(tftftpmtvSD 時間時間t時平衡方程的重新排列；時平衡方程的重新排列； (4)等效荷載增量等效荷載增量 和等效剛度和等效剛度 按方程按方程(7-9*)計算；計算；)(tp)(tk高等結(jié)構(gòu)動力學(xué) (5)由方程（由方程（7-8*）可求得位移增量，從而由方程）可求得位移增量，從而由方程(7-7b

27、*)可求可求得速度增量；得速度增量； (6)最后，獲得時段終點速度和位移：最后，獲得時段終點速度和位移：（7-26a7-26a）)()()()()()(tvtvttvtvtvttv（7-26b7-26b）高等結(jié)構(gòu)動力學(xué) 當(dāng)?shù)诋?dāng)?shù)?步運算完成時，這個時段的分析結(jié)束；下一個時段的步運算完成時，這個時段的分析結(jié)束；下一個時段的分析，只需將上述整個程序重復(fù)進行即可。顯然，連續(xù)進行上分析，只需將上述整個程序重復(fù)進行即可。顯然，連續(xù)進行上述運算，就可進行任意多個時間增量的分析，這樣就能算出具述運算，就可進行任意多個時間增量的分析，這樣就能算出具有任何非線性性質(zhì)的單自由度體系的全部反應(yīng)時程。線性體系有任何非

28、線性性質(zhì)的單自由度體系的全部反應(yīng)時程。線性體系自然也可以用同樣的方法進行處理。此時阻尼和剛度特性保持自然也可以用同樣的方法進行處理。此時阻尼和剛度特性保持不變，因此分析的過程將稍為簡單一些。不變，因此分析的過程將稍為簡單一些。高等結(jié)構(gòu)動力學(xué) 逐步積分法的精度，依賴于時間增量逐步積分法的精度，依賴于時間增量t的長的長度，選取時注意三個因素：度，選取時注意三個因素： (1)作用荷載作用荷載p(t)的變化速率；的變化速率； (2)非線性阻尼和剛度特性的復(fù)雜性；非線性阻尼和剛度特性的復(fù)雜性； (3)結(jié)構(gòu)振動周期結(jié)構(gòu)振動周期T。 高等結(jié)構(gòu)動力學(xué) 為可靠地反映這些因素，時間增量必須足為可靠地反映這些因素，

29、時間增量必須足夠的短，其中最后一個因素是和體系的自由振夠的短，其中最后一個因素是和體系的自由振動特性聯(lián)系在一起的。一般來說，材料特性的動特性聯(lián)系在一起的。一般來說，材料特性的變化不是關(guān)鍵性的因素，如果發(fā)生一個重大的變化不是關(guān)鍵性的因素，如果發(fā)生一個重大的突然變化，例如一個彈塑性彈簧的屈服，這時突然變化，例如一個彈塑性彈簧的屈服，這時可以引入一個特殊再細(xì)分的時間增量來精確地可以引入一個特殊再細(xì)分的時間增量來精確地處理這個影響。同時，要估算能恰當(dāng)?shù)啬M動處理這個影響。同時，要估算能恰當(dāng)?shù)啬M動力荷載主要形狀的時間增量也并不困難。力荷載主要形狀的時間增量也并不困難。高等結(jié)構(gòu)動力學(xué) 因此，如果加載過程

30、比較簡單，時間間隔的選取主要依賴于因此，如果加載過程比較簡單，時間間隔的選取主要依賴于結(jié)構(gòu)的振動周期。這個線性加速度方法只是有條件穩(wěn)定的，如結(jié)構(gòu)的振動周期。這個線性加速度方法只是有條件穩(wěn)定的，如果時間增量大于振動周期的一半左右，則將給出擴散的解。但果時間增量大于振動周期的一半左右，則將給出擴散的解。但是，為了提供適當(dāng)?shù)木龋瑫r間增量必須比這個短的多，因此是，為了提供適當(dāng)?shù)木龋瑫r間增量必須比這個短的多，因此不穩(wěn)定是不會成為問題的。一般來說，按照經(jīng)驗，如果增量不穩(wěn)定是不會成為問題的。一般來說，按照經(jīng)驗，如果增量一一周期比周期比t/T1/10，則可獲得可靠的結(jié)果。如果對于得出的則可獲得可靠的結(jié)果。

31、如果對于得出的解有任何懷疑，則第二次分析時取時間增值的一半進行，解有任何懷疑，則第二次分析時取時間增值的一半進行，如如果在第二次分析中反應(yīng)沒有明顯的變化，果在第二次分析中反應(yīng)沒有明顯的變化，則可以認(rèn)為數(shù)值積分所產(chǎn)生的誤差是可以忽略不計的。則可以認(rèn)為數(shù)值積分所產(chǎn)生的誤差是可以忽略不計的。高等結(jié)構(gòu)動力學(xué)例題例題 E7-2 為了說明應(yīng)用上述線性加速度逐步積分法的手算方法為了說明應(yīng)用上述線性加速度逐步積分法的手算方法,我們我們來計算圖來計算圖E7-3所給出的彈塑性單自由度剛架在所給的加荷過程下的反應(yīng)所給出的彈塑性單自由度剛架在所給的加荷過程下的反應(yīng).在分析中在分析中,時間步長取為時間步長取為0.1秒秒

32、,它比要得到精度較高的結(jié)果所需的步長較它比要得到精度較高的結(jié)果所需的步長較長一些長一些,但為了說明這個手算方法但為了說明這個手算方法,這個步長已經(jīng)足夠了這個步長已經(jīng)足夠了. 在這個結(jié)構(gòu)里在這個結(jié)構(gòu)里,假定阻尼系數(shù)保持常量假定阻尼系數(shù)保持常量,因此非線性影響僅由于發(fā)因此非線性影響僅由于發(fā)生屈服時剛度將改變才發(fā)生的生屈服時剛度將改變才發(fā)生的.此時等效剛度可以被表達為此時等效剛度可以被表達為參見方程參見方程(8-9c) 263660.10.1kktmckt其中其中k(t)或是或是5千磅千磅/英寸或是零英寸或是零,它們分別相應(yīng)于剛架的彈性和屈服階段它們分別相應(yīng)于剛架的彈性和屈服階段,等效增量荷載由下式給出等效增量荷載由下式給出參見方程參見方程(7-9b*): 60.131mp tp tcvmcvp tvv 由方程由方程(7-7b*)給出的速度增量成為給出的速度增量成為:高等結(jié)構(gòu)動力學(xué)3030.05vvvv便于手算的表格示于表便于手算的表格示于表E7-1圖圖 E7-3 彈塑性剛架和動力