【正文】 ( ) . 1u ? 線性加速度法假定： 按線性規(guī)律變化（如圖47 ）。 ???這里應指出：線性加速度法是一種有條件的穩(wěn)定數(shù)值積分方法，穩(wěn)定性條件為時間步長 當所取得時間步長太大時，就可能出現(xiàn)不收斂的情況。.5 .3 4 線性加速度法4 . 1 6t? ????線性加速度法假定加速度在 內(nèi)按線性變化 如圖 （a ) 所示即有：1()( ) ( ) ( ) ( )ii i i i iutu t u t t t t t t t tt ??? ? ? ? ? ? ? ??（ ） 將上式積分一次，可得2( ) ( )( ) ( ) ( ) ( ) +2iii i it t u tu t u t u t t tt??? ? ??（ ） [ 4 . 1 6 ]速度按二次拋物線圖 （b ）再積分一次，有23( ) ( ) ( ) ( )( ) ( ) ( ) ( ) + +26i i i ii i it t u t t t u tu t u t u t t tt? ? ?? ? ?（ ） 位移：三次曲線. 1 . 1 it t t? ? ?在式（4 4 3 ）和式（4 4 4 ）中，令 ，則有+1( ) = ( + ) ( ) = ( ) ( ) ( ) ( ) 2i i i i i i itu t u t t u t u t u t u t t u t ?? ? ? ? ? ?（ ） 22( ) = ( ) + ( ) ( )26i i i ittu t u t t u t u t??? ? ? ?（ ） ( ) ( ) ( )i i iu t u t u t? ? ?為便于計算，將 ， 用 表示，故有266( ) = ( ) ( ) 3 ( )i i i iu t u t u t u ttt? ? ???（ ） 33( ) = ( ) ( ) ( )2i i i itu t u t u t u ttt?? ? ???（ ） .1 itt?現(xiàn)將式（4 3 2 ）中令 ，并將以上兩式代入，得26 6 3 3( ) ( ) 3 ( ) ( ) ( ) ( )2( ) ( )i i i i i iiitm u t u t u t c u t u t u tt t ttk u t P t?? ? ? ?? ? ? ? ?? ? ? ?? ? ??? ? ? ?? ? ? ?（ ） .1將上式整理后，可得與式（4 3 9 ）相同的方程( ) ( )iiK u t P t? ? ? （ ） 式中236=K k c mt t??? ?（ ） 6( ) ( ) ( ) 3 ( ) 2 3 ( ) + ( )2i i i i i iktP t P t m u t u t c u t u tt?? ? ? ?? ? ? ? ? ?? ? ? ??? ? ? ?（ ） +1.1 ( ) .1.1.1 =iiiuti t t t???有式（4 50）求出位移增量 后，代入式（4 47）和式（4 48），便可求出速度增量和加速度增量，然后，便可按式（4 33）求出 時段末 的動力反應值。對于多自由度體系，只需將以上各式中的相應物理量改寫成矩陣或向量的形式即可。. 5 . 24 平 均 常 加 速 度 法1( ) + ( + )4 . 1 5 ,i i i i itt u t t t t u t t????????考 慮 一 單 自 由 度 體 系 ， 假 定 在 內(nèi) 質(zhì) 點 加 速 度 為 常 數(shù) ， 它等 于 質(zhì) 點 在 時 刻 的 加 速 度 和 = 時 刻 的 加 速 度的 平 均 值 圖 （ a ） 即 有+1( ) + ( )( + ) =2iiiiu t u tu t t t?，（ ） 1.1ittt t t ??? ? ?故質(zhì)點的速度在 時段內(nèi)呈線性變化[ 圖4 5 （b ）] ，并且在 時刻的速度為1( ) ( )( ) = ( ) + ( ) = ( ) +2iii i i i iu t u tu t t u t u t t t t u t t??? ? ? ? ? ?，（ ） .1t?位移在 時段按拋物線變化[ 圖4 5 （c ）] ，并且在時段末的加速度為21( ) ( )( ) = ( ) + ( ) +4iii i iu t u tu t t u t u t t t??? ? ? ?（ ） .1 tt??對于式（4 2 2 ）的單質(zhì)點運動方程，在時刻 處，則有( ) ( ) ( ) ( )m u t t c u t k u t t P t t? ? ? ? ? ? ? ? ? ?（ ） .1 .2式（4 31 ）與式（4 2） 相減，可以得到增量形式表示的振動方程為( ) ( ) ( ) ( )m u t c u t k u t P t? ? ? ? ? ? ?（ ） 其中( ) = ( ) ( )( ) = ( ) ( )( ) = ( ) ( )( ) = ( ) ( )u t u t t u tu t u t t u tu t u t t u tP t P t t P t? ? ? ??? ? ? ??? ? ???? ? ??（ ） .1將速度按式（4 2 9 ）改寫成增量形式11( ) ( )= ( ) ( ) =2( ) ( )21( ) ( ) ( )2 2 2iii i iiii i iu t u tu u t t u t tu t u tu t t t u t t u t t???? ? ? ??? ? ? ? ? ? ? ?（ ） 同理，位移增量形式為2211= ( ) ( ) + ( )24i i iu u t t u t t u t t? ? ? ? ?（ ） .1由式（4 32 ）得( ) = ( ) ( )c k Pu t u t u tm m m?? ? ? ? ?（ ） .1由式（4 3 3 ）求得加速度增量244( ) = ( ) ( ) 2 ( )iiku t u t u t u ttt? ? ???（ ） .1并回代入式（4 3 4 ），有2= 2u u ut???（ ） . 1 . 1 . 1將式（4 3 7 ）和式（4 3 8 ）代入式（4 3 2 ），得( ) ( )iiK u t P t? ? ?（ ） 式中224=K k c mt t??? ?（ ） 4( ) ( ) ( ) 2 ( ) 2 ( )i i i i ikP t P t m u t u t c u tt??? ? ? ? ? ??????（ ） .1 .1式（4 40 ）和式（4 41 ）分別稱為等效剛度和等效增量荷載。11 ,t T T???以 上 的 計 算 公 式 具 有 2 階 精 度 ， 其 穩(wěn) 定 條 件 為為 體 系 的 最 小 自 振 周 期 。對 于 多 自 由 度 體 系 ， 只 需 將 參 數(shù) 和 改 寫 成 矩 陣 或向 量 的 形 式 即 可（ ） （ ） 00 ( 0 ) , ( 0 )u u u u??對 于 初 始 條 件 為 的 單 質(zhì) 點 振 動 方 程 ， 中 心差 分 法 的 具 體 迭 代 步 驟 為 ：0 0 0 01 ( ) ,u P c u k um? ? ?1 ） 計 算 基 本 數(shù) 據(jù) 和 初 始 條 件21 0 0 02tu u tu u??? ? ? ?12 ) . 1iitt ?根 據(jù) 及 以 前 時 刻 的 運 動 ， 利 用 式 （ 4 2 7 ） 計 算 時 刻 的1iu ?位 移 。. 5 . 14 中 心 差 分 法t處 ， 單 質(zhì) 點 的 震在 時 刻 動 方 程 為( ) ( ) ( ) ( )m u t c u t k u t P t? ? ?（ ） 1i i ifbi i it t t tt u u?? ? ?如 果 采 用 等 時 間 步 長 （ = ） , 則 速 度 在 時 刻處 得 向 前 差 分 和 向 后 差 分 分 別 為11= , =fbi i i iiiu u u uuutt??????所 以 其 中 心 差 分 為1122fbi i i iiu u u uut???????同 理 ， 可 得 其 加 速 度 的 中 心 差 分 近 似 為1122i i iiu u uut??????（ ） （ ） （ ） . 1 . 1 . 1將 式 （ 4 2 4 ） 和 式 （ 4 2 5 ） 代 入 式 （ 4 2 2 ） 有1 1 1 122 +2i i i i iiiu u u u um c k u Ptt? ? ? ?? ? ? ????1 1 11i i i iit u u tu? ? ??如 果 已 知 時 刻 ， 以 前 的 位 移 ， ， 則 時 刻 的位 移 可 由 下 式 求 解 ：112 2 2222i i i im c m c mu P k u ut t t t t??? ? ? ? ? ?? ? ? ? ? ?? ? ? ? ? ?? ? ? ? ?? ? ? ? ? ?1 . 1 . 1,itm c k P?而 時 刻 的 速 度 和 加 速 度 則 可 由 式 （ 4 2 4 ） 和 式 （ 4 2 5 ）求 得 。計算精度，即截斷誤差與時間步長 的關(guān) 系，如果誤差 則稱為該方法具有N 階精度。 法、 法110( ) ,NNttEt??? ? ??????????? ? ? ????收斂性，即當積分步長 時，數(shù)值 解是否收斂于精確解。N e w ma r kWi l s o n???????????????分段解析法中心差分法平均加速度法逐步積分法線性加速度法法法N e w m a rk W ilso n?????????????????顯式方法：在每一步內(nèi)計算新的反應值僅僅依賴于前面步驟已經(jīng) 獲得的量。但這兩種方法均是基于疊加原理，要求結(jié)構(gòu)體系是線彈性的，而且當載荷沒有解析表示式（如地震荷載）或解析式太復雜時，疊加原理將不再適用，也就是說， 積分、 變換或振型疊加法此時失效，因此只能采用數(shù)值分析方法。= 0 . 0 0 0 21 .0 0 0 . 0 0 0 20 .0 6 0 6 1 .0 0ZZZZZ???????????????22解得，對于 ，有：對于 ，有：故有220 .3 8 4 9 1 .0= = 0 .7 7 4 6 5 0 .4 5 1 4 61 .0 4 8 9 2 0 .8 0 7Z???????? ???2因此，第二次近似的一、二階振型為12,0 .3 6 6 9 1 .0= = 0 .7 3 8 5 0 .4 5 1 41 .0 0 .8 0 6 9= 0 .1 9 8 = 1 .5 5 4 9??????????? ???222歸一化后 有由第二輪迭代得到的 ， 與真實解差別已經(jīng)很小了。= 0 .3 5 1ZZZZ??????22解得：，對于 ，有：對于 ，有：1 51 Z???? ?????故有110 .4 4 9 8 0 .8 4 2 1= = 0 .7 9 9 9 0 .3 8 5 91 .0 0 .6 8 4 3Z???????? ???1因此，第一近似的一、二階振型為,0 .4 4 9 8 1 .0= = 0 .7 9 9 9 0 .4 5 8 31 .0 0 .8 1 2 6??????????1歸一化后 有220 .4 4 5 5 10 .8 0 2 0 .4 5 1 31 .0 0 0 .8 0 7 2?????????? ? ??????11再對 左乘算子 得作廣義質(zhì)量振型和廣義剛度矩陣的第二近似為2 2 22 2 21 .8 4 1 7 0 . 0 0 0 2()0 . 0 0 0 2 1 .8 5 5 30 .3 6 4