大跨度橋梁非線(xiàn)性分析的論文
摘要：本文在前人研究的基礎上提出了統一的顫振和抖振分析方法。該方法以非線(xiàn)性有限元的直接積分法為基礎，在研究中具體解決了隨機風(fēng)速場(chǎng)的模擬、耦合自激力的時(shí)域計算和統一的顫抖振時(shí)程分析流程等關(guān)鍵問(wèn)題，考慮了結構的幾何非線(xiàn)性和有效攻角效應。本文的研究糾正了過(guò)去時(shí)程分析方法不能同時(shí)處理顫振和抖振的理論缺陷。本文還通過(guò)所編制的軟件的計算實(shí)例驗證了方法的正確性和可行性。
關(guān)鍵詞：橋梁非線(xiàn)性顫振抖振時(shí)程分析
一、前言
時(shí)程分析方法是橋梁風(fēng)工程中的主要方法之一。過(guò)去的非線(xiàn)性時(shí)域分析方法都局限于抖振。其基本流程是首先模擬橋梁風(fēng)場(chǎng)的脈動(dòng)風(fēng)速時(shí)程，根據脈動(dòng)風(fēng)速計算抖振力和自激力，然后將抖振力和自激力的計算編入非線(xiàn)性有限元程序中，最后再運用這樣的程序進(jìn)行計算。在這個(gè)流程中，非線(xiàn)性有限元程序是比較成熟的，但在脈動(dòng)風(fēng)速模擬和自激力的計算方面都還存在著(zhù)對分析有重要影響的缺陷。由于時(shí)域中耦合自激力的計算比較困難，過(guò)去的時(shí)程分析中都沒(méi)有考慮耦合的自激力，因此，這樣的分析方法不能用來(lái)分析耦合顫振[2]。
本文在此對時(shí)程分析方法進(jìn)行了改進(jìn)。首先是改進(jìn)了模擬隨機風(fēng)場(chǎng)的諧波合成法，提高了模擬的效率。然后本文實(shí)現了時(shí)域中耦合自激力的計算，從而在時(shí)域中實(shí)現了比較完善的風(fēng)荷載計算。利用這樣的風(fēng)荷載，本文在時(shí)域中統一了抖振和顫振的分析方法。在時(shí)域中實(shí)現了耦合顫抖振和顫振分析。根據這一方法，本文運用可視化編程技術(shù)，編制了大跨度橋梁非線(xiàn)性顫振和抖振時(shí)程分析的有限元程序Nbuffet，并對程序進(jìn)行了驗證。最后本文對江陰長(cháng)江大橋進(jìn)行了非線(xiàn)性顫振和抖振分析，得出了一些有益的結論。
二、脈動(dòng)風(fēng)送的.模擬
要進(jìn)行抖振時(shí)程分析就必須首先模擬作用在橋梁上的脈動(dòng)風(fēng)速時(shí)程。本文采用經(jīng)作者改進(jìn)的諧波合成系列中的一種方法，大大提高了模擬效率，為在后文進(jìn)行顫振時(shí)程分析中不斷變換風(fēng)速計算節約了時(shí)間。作用在大跨度橋梁上脈動(dòng)風(fēng)速可視為一維多變量隨機過(guò)程。眾所周知，用諧波合成法模擬一維多變量隨機過(guò)程需要計算互譜密度矩陣的Cholesky分解。該分解通常采用迭代法求借，計算最大，常常影響模擬的規模的效率。本文作者利用橋梁上各點(diǎn)的互譜密度近似相等的特點(diǎn)，導出了顯式的分解公式，并且采用了FFT技術(shù)，從而極大地提高了模擬效率。
三、風(fēng)荷載計算
引起橋梁風(fēng)振的荷載可以分為靜力風(fēng)荷載、抖振力和自激力。其中靜力荷載按常規靜力三分力系數計算，抖振力常按Scanlan的準定常理論計算。
自激力的計算一直是研究得較多的課題之一。傳統頻域抖振和顫振分析方法中的自激力都采用Scanlan提出的氣動(dòng)導數的線(xiàn)性表達式。由于該表達式是頻域和時(shí)域的混合表達式，不能在時(shí)域中求解。為了在時(shí)域中順利計算耦合自激力，Lin提出了一種用單位脈沖響應函數表達的統一自激力表達式[4]。本文按Lin的理論計算耦合自激力。Lin的理論基于二自由度耦合。然而，三自由度耦合對結構振動(dòng)的影響最近也引起了一些學(xué)者的關(guān)注。雖然并非所有的自由度之間都具有耦合特性，但從理論和形式完備的角度出發(fā)，本文將Lin的理論從二自由度推廣到三自由度，成功地實(shí)現了時(shí)域內三自由度耦合自激力的計算。
用脈沖響應函數表達的自激力適合于任意形式的振動(dòng)，也適用于正余弦振動(dòng)（顫振）。根據在正余弦振動(dòng)形式下，脈沖響應函數表達的自激力與氣動(dòng)導數表達的自激力相等價(jià)的關(guān)系，Lin導出了用脈沖響應函數表達的自激力的具體表達形式。
四、統一的額報和抖報時(shí)域分析方法
在傳統的步域分析方法中，抖振和顏振是通過(guò)完全不同的方法來(lái)分析的。其中，抖振分析用的是基于隨機振動(dòng)理論的響應譜方法，顫振分析用的是與特征值問(wèn)題有關(guān)的半逆解法或復模態(tài)解法。風(fēng)振時(shí)程分析的初衷是為了解決非線(xiàn)性情況下的抖振響應計算。但是顫振分析中所需要的計算自激力的公式在抖振時(shí)程分析中都要用到，所以從理論上講，利用計算抖振時(shí)程分析的方法同樣可以在時(shí)域中計算顫振。實(shí)際上，抖振和顫振并不是完全獨立的。在任何風(fēng)速之下，橋梁都受到抖振力和自激力的作用。當風(fēng)速較低時(shí)，自激力很小，不起控制作用，橋梁的振動(dòng)就體現為抖振。當風(fēng)速增加到一定程度時(shí)，自激力逐漸發(fā)散，并控制橋梁的運動(dòng)，橋梁就發(fā)生了顫振。因此，只要正確地描述了抖振力和自激力，運用時(shí)程分析這一仿真的分析方法，就可以算出一定風(fēng)速之下橋梁的真實(shí)運動(dòng)狀態(tài)。如果表現為隨機振動(dòng)，則說(shuō)明是抖振，我們就可以得到響應時(shí)程統計指標。如果是發(fā)散振動(dòng)，就說(shuō)明橋梁發(fā)生了顫振。只要不斷進(jìn)行搜索計算，我們就能在時(shí)域中找到橋梁的顫振臨界風(fēng)速。
根據以上設想，本文設計并首次成功地實(shí)現了時(shí)域中統一的顫振和抖振分析算法。
流程中，耦合自激力的計算是個(gè)關(guān)鍵。過(guò)去的一些抖振時(shí)程分析方法中常只近似考慮非耦合的自激力。而大跨度橋梁的顫振發(fā)散大多是受耦合自激力控制的，因此，過(guò)去的抖振時(shí)程分析方法不能用于計算顫振的原因就在于此。顫振發(fā)散的判斷依據也是關(guān)鍵之一�？紤]到結構在接近顫振臨界狀態(tài)時(shí)，振動(dòng)形式逐漸從隨機振動(dòng)過(guò)渡到諧波發(fā)散振動(dòng)，其振幅將逐漸增大，相應振動(dòng)的阻尼將逐漸減小。因此，本文先通過(guò)位移時(shí)程曲線(xiàn)觀(guān)察振幅的變化規律，當結構的振動(dòng)明顯過(guò)渡為諧波振動(dòng)時(shí)，則根據計算結構的阻尼系統，當阻尼系統為負時(shí)，則認為結構進(jìn)入顫振臨界狀態(tài)。計算實(shí)例表明，這種判斷方法與其他方法計算得到的結果一致。
五、非線(xiàn)性顫振和抖振時(shí)程分析的程序設計
除了在時(shí)域中統一顫振和抖振分析方法以外，本文研究時(shí)程分析方法的目的還在于分析不同非線(xiàn)性因素對橋梁顫振和抖振響應的影響。與大跨度橋梁抖振和顫振有關(guān)的非線(xiàn)性現象主要有：
�。�1）幾何非線(xiàn)性，包括平均風(fēng)荷載引起的位移：由于大跨度橋梁相對細長(cháng)，幾何非線(xiàn)性現象不能忽視；
�。�2）有效攻角效應：由平均風(fēng)荷載引起的位移使風(fēng)對橋梁的攻角發(fā)生變化，從而使靜力三分力系統和氣動(dòng)導數發(fā)生變化，因此附加攻角對橋梁的影響不能忽視。
根據以上分析流程并考慮這些非線(xiàn)性因素，借鑒一些通用有限元程序的理論和源代碼[5]，本文作者編制了大跨度橋梁顫振和抖振分析程序Nbuffet。該程序以FortranPowerStation（FPS）4。0為平臺，采用Fortran90語(yǔ)言編程。作者運用了FPS的Windows編程技術(shù)，使Nbuffet成為一個(gè)基于Windows95／NT平臺具有豐富的交互式功能的實(shí)用程序。
由于目前非線(xiàn)性有限元技術(shù)相對比較成熟，該部分在理論上不是本文的重點(diǎn)，因此這里不再詳述。
六、實(shí)例分析
在以上理論的基礎上，本文作者編制了相應計算機程序Nbuffet。該程序考慮了結構的幾何非線(xiàn)性和氣動(dòng)非線(xiàn)性（有效攻角引起的三分力和氣動(dòng)導數等變化），以便可以考慮這些非線(xiàn)性對結構風(fēng)振行為的影響。本文作者在程序中采用魚(yú)骨架式模型建立大跨度橋梁模型，采用桿梁的切線(xiàn)剛度矩陣和Newton—Raphson方法并引入平衡迭代來(lái)處理結構幾何非線(xiàn)性。運用所編制的程序，本文分析了江陰長(cháng)江大橋主橋的非線(xiàn)性顫振和抖振行為。
江陰長(cháng)江大橋主跨1385m，是我國目前在建的跨度最大的橋梁。豐文運用Nbuffet程序，分析了該橋不同參數下的顫振和抖振響應，并與用其他方法得到的結果進(jìn)行了比較。結果顯示，本文建立的統一的顫振和抖振分析方法在理論上和實(shí)踐上都是成功的。本文所編制的Nbuffet程序也是實(shí)用可靠的。以下分別是運用傳統頻域分析方法、風(fēng)洞模型試驗和本文的方法分析得出的一些結果對比情況。限于篇幅所限。從結果對比可以看出本文的計算結果與頻域分析方法、風(fēng)洞模型試驗的結果基本吻合。本文的主要目的是建立一套時(shí)域內顫振和抖振統一分析的方法和流程。從比較結果來(lái)看，這種方法和流程是成功的。
從比較結果中還可以得到以下一些現象：
�。�1）本文豎向響應略小于風(fēng)洞試驗結果，本文的扭轉結果又略大于風(fēng)洞試驗結果�？紤]到目前的風(fēng)振試驗和分析方法體系都尚未達到比較精確的程度，這些誤差可能來(lái)源于試驗、頻域、時(shí)域三者之間的模型誤差。
�。�2）素流對該橋的顫振臨界風(fēng)速沒(méi)有影響，即考慮抖振項的參與不影響該橋額振臨界風(fēng)速。
�。�3）只有氣動(dòng)導納因素對抖振結果影響顯著(zhù)�？梢�(jiàn)，幾何非線(xiàn)性和有效三分力及有效氣動(dòng)導數對懸索橋的影響可能要到更大的跨度才能表現出來(lái)。
七、結語(yǔ)
大跨度橋梁在非線(xiàn)性情況下的顫振和抖振分析是目前橋梁風(fēng)工程研究的熱點(diǎn)之一。本文著(zhù)重提出了時(shí)域中統一的顫振和抖振方法，同時(shí)解決了脈動(dòng)風(fēng)速的高效率模擬、結構幾何非線(xiàn)性和氣動(dòng)非線(xiàn)性的處理方法。在此基礎上，本文編制了計算程序Nbuffet并用該程序分析了江陰長(cháng)江大橋非線(xiàn)性顫振和抖振響應。結果表明本文提出的方法及所編制的程序在理論和實(shí)踐上都是正確的。
在此基礎上，我們就可以在時(shí)域中增加考慮各種非線(xiàn)性因素對結構進(jìn)行分析從而尋找結構對這些因素的敏感性；我們也可以根據時(shí)程計算來(lái)進(jìn)行非線(xiàn)性的振動(dòng)控制。而這些研究工作在頻域范圍內是難以開(kāi)展的。如果與CFD技術(shù)相結合，將可望實(shí)現從參數識別到結構宏觀(guān)計算和控制的全過(guò)程分析。從而達到與風(fēng)洞試驗互為補充的目的。
應該說(shuō)，盡管以上方法和程序是成功的，但是更重要的是要利用這種方法對所關(guān)心的橋梁進(jìn)行各種用過(guò)去的方法所不能進(jìn)行的全過(guò)程參數分析，從而得到更具有普遍規律性的結論。因此，大量的實(shí)例計算和總結是必要的。
【大跨度橋梁非線(xiàn)性分析的論文】相關(guān)文章：
大跨度橋梁設計的論文11-11
公路橋梁中大跨度橋梁設計研究論文10-22
市區橋梁施工技術(shù)分析的論文07-18
關(guān)于橋梁風(fēng)險分析方法研究論文10-28
跨度空間結構研究的論文10-26
施工控制及橋梁施工技術(shù)分析論文12-01
橋梁鉆孔灌注樁的質(zhì)量分析論文12-02