基于邊界元法的膜結構流固耦合效應數(shù)值模擬.pdf_第1頁
已閱讀1頁,還剩61頁未讀 繼續(xù)免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

1、膜結構和風場之間的耦合作用研究是膜結構抗風研究的難點。本文基于小擾動位勢理論,應用非協(xié)調邊界元法(Boundary Element Method)和有限元法(Finite Element Method),采取弱耦合策略對風-膜結構進行了流固耦合分析。
   目前的CFD軟件大都基于有限體積法(Finite Volume Method),需要對整個流場進行單元劃分,特別在解決對計算精度要求較高的問題中,劃分出的單元量巨大。另外,膜

2、結構變形的幾何非線性程度大,會在計算過程中引起初始劃分的單元奇異,必須重新劃分單元才可以繼續(xù)計算。所以單元劃分耗費了人們大量精力,導致總體計算時間很慢。邊界元法僅僅在流場邊界上離散計算區(qū)域,大大減少了未知量,節(jié)省了流體域求解的機時。膜結構作為一種曲面結構,其與流體的耦合恰恰發(fā)生在流場的邊界,邊界解可以直接用于膜結構變形計算,這恰好利用了邊界元的優(yōu)勢。
   文中要計算的流場域邊界非連續(xù)光滑,有尖角和尖邊存在,故所要求解場量的法向

3、導數(shù)不連續(xù),進而難以確定邊界條件。如果直接應用普通的邊界元進行計算將難以保證計算精度,于是本文引入了非協(xié)調邊界元技術,較好地解決了這一問題,使得計算精度大大提高。而且,本文編制的邊界元計算程序,可以根據(jù)源點處邊界面單元法向量的變化情況自動識別是否需要引入非連續(xù)邊界元,并重新劃分計算域。
   對于系數(shù)矩陣不對稱的大型邊界元方程組,必須選擇準確高效的數(shù)值求解方法,否則難以保證計算精度和計算效率。本文引入了非常高效的循環(huán)型廣義極小殘

4、余迭代算法(GMRES(m)) ,它是計算數(shù)學中共軛梯度法的一種。除此之外,為了防止系數(shù)矩陣條件數(shù)較大導致方程組病態(tài),本文應用對角線法對系數(shù)矩陣進行了預處理。上述方法的應用使得邊界元法的優(yōu)勢得到了充分發(fā)揮。若干算例驗證了上述理論的正確性和程序的可靠性。
   此外,求解三維空間曲面邊界單元上場量的梯度是本文流固耦合求解過程中的一個關鍵點和難點,本文給出了詳細的推導過程。
   最后,本文采用點渦模型模擬特征湍流,以修正已

5、有的基于邊界單元法的數(shù)值求解結果。文中共計算分析了菱形馬鞍屋面在0°和90°兩種風向角情況下的結構豎向位移,并將計算結果與本課題組在湖南大學完成的氣彈模型風洞試驗結果進行對比分析。結果發(fā)現(xiàn),低應力狀態(tài)下兩者的誤差較大,而高應力狀態(tài)下兩者誤差較小。主要由于:膜結構氣彈模型風洞試驗目前還處于起步階段,材料、測量器材、模型加工精度、試驗控制設備等諸多限制使得膜結構氣彈模型風洞試驗的精度大打折扣,再加上理論計算的諸多假定,導致了兩者之間的誤差。

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
  • 5. 眾賞文庫僅提供信息存儲空間,僅對用戶上傳內容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論