一階有限元方法研究以及圓柱渦振問題的數(shù)值模擬
本文選題:有限元 + 分區(qū)伽遼金; 參考:《湖南大學(xué)》2016年博士論文
【摘要】:本文在總結(jié)以往學(xué)者在彈性力學(xué)及流體力學(xué)有限元方法的基礎(chǔ)之上,提出了基于分區(qū)伽遼金方程的變分方法。針對(duì)常規(guī)有限元方法導(dǎo)致彈性力學(xué)應(yīng)力精度下降的問題,提出了一階有限元改進(jìn)算法。針對(duì)不可壓Navier-Stokes方程求解的困難,提出了不含壓力項(xiàng)的一階流體動(dòng)力學(xué)方程系統(tǒng)。并對(duì)高雷諾數(shù)下高質(zhì)量比圓柱渦致振動(dòng)現(xiàn)象進(jìn)行了數(shù)值模擬和分析。本文主要包括以下一些工作:(1)提出了分區(qū)伽遼金方程,放松了分區(qū)交界面上位移、應(yīng)力連續(xù)的條件。分析了弱形式降低連續(xù)性的現(xiàn)象。建立了基于分區(qū)伽遼金方程(分區(qū)加權(quán)殘數(shù)法)的求解體系,為構(gòu)造各種單元,特別是擬協(xié)調(diào)元、雜交元提供了理論基礎(chǔ)。對(duì)于彈性力學(xué)問題,在統(tǒng)一的構(gòu)架下,基于分區(qū)伽遼金方程,導(dǎo)出了分區(qū)弱形式、分區(qū)廣義虛功方程和分區(qū)變分原理。分析了積分形式解的組成模式。提出了選取權(quán)函數(shù)要滿足的條件;诜謪^(qū)伽遼金方程的變分方法為其后的一階有限元方法提供了理論依據(jù)。(2)將一階有限元方法應(yīng)用到二維、三維彈性力學(xué)問題中。首先推導(dǎo)了彈性力學(xué)問題的一階弱形式,然后使用Free Fem++軟件完成有限元編程工作。并且采用典型算例來(lái)比較一階算法和常規(guī)有限元算法的精度。通過數(shù)值算例,驗(yàn)證了一階有限元方法使應(yīng)力精度與位移精度同階,應(yīng)力的精度得到了提高。同時(shí)數(shù)值計(jì)算的結(jié)果還表明采用一階算法,達(dá)到相同的應(yīng)力精度,比常規(guī)有限元法花費(fèi)的時(shí)間要少。一階有限元算法為有限元應(yīng)力精度的提升提供了一個(gè)新的思路。(3)介紹了不可壓Navier-Stokes方程,介紹了特征線分裂算法,并編程將其用于二維圓柱繞流的計(jì)算模擬。并將一階有限元解法應(yīng)用到不可壓粘性流體計(jì)算中。對(duì)于不可壓粘性流動(dòng),提出了不含壓力項(xiàng)的一階流體動(dòng)力學(xué)方程系統(tǒng)。基于有限元方法,對(duì)應(yīng)力和速度采用同階插值,對(duì)兩平行平板間的穩(wěn)態(tài)粘性流動(dòng)及二維非定常圓柱繞流進(jìn)行了數(shù)值計(jì)算,并分別和精確解以及標(biāo)準(zhǔn)測(cè)試算例進(jìn)行對(duì)比。采用有限元方法對(duì)提出的不含壓力項(xiàng)的一階流體動(dòng)力學(xué)方程系統(tǒng)進(jìn)行求解,可以對(duì)應(yīng)力和速度采用同階插值,從而可以避免Navier-Stokes方程求解過程中反復(fù)使用速度導(dǎo)數(shù)而導(dǎo)致精度下降的問題。(4)研究了高雷諾數(shù)下高質(zhì)量比圓柱的渦致振動(dòng)問題。對(duì)所使用的流體SST湍流模型理論以及一些相關(guān)的參數(shù)做了詳細(xì)的介紹。采用基于伽遼金最小二乘有限元法的CFD計(jì)算軟件,通過使用質(zhì)量-彈簧-阻尼系統(tǒng)以及SST湍流模型進(jìn)行了流固耦合數(shù)值模擬。流固耦合分析結(jié)果與文獻(xiàn)試驗(yàn)結(jié)果基本吻合,驗(yàn)證了采用的流固耦合計(jì)算方法的正確性。從而可以為該類問題的數(shù)值模擬提供參考。數(shù)值模擬的結(jié)果表明:在計(jì)算邊界層由于強(qiáng)烈逆壓梯度而引起的分離流動(dòng)問題時(shí),采用SST湍流模型是比較合適的;并且,在高雷諾數(shù)下,高質(zhì)量比圓柱會(huì)出現(xiàn)高幅分支現(xiàn)象,從而為該類工程問題的研究提供參考。
[Abstract]:In this paper , based on the elastic mechanics and hydromechanics finite element method , this paper presents a variational method based on the partition - Galerkin equation . A first - order finite element method is proposed to solve the problem of decreasing the stress precision of the two - dimensional cylinder .
【學(xué)位授予單位】:湖南大學(xué)
【學(xué)位級(jí)別】:博士
【學(xué)位授予年份】:2016
【分類號(hào)】:O302
【相似文獻(xiàn)】
相關(guān)期刊論文 前10條
1 劉鳴放;車穎濤;;理工院校相關(guān)專業(yè)增設(shè)有限元方法選修課程的可行性探討[J];高等函授學(xué)報(bào)(自然科學(xué)版);2010年02期
2 劉鳴放;車穎濤;;理工類部分本科專業(yè)增設(shè)《有限元方法》選修課程的可行性探討[J];商丘職業(yè)技術(shù)學(xué)院學(xué)報(bào);2010年02期
3 張斐然;;大學(xué)理工類本科專業(yè)增設(shè)《有限元方法》選修課之探討[J];商丘師范學(xué)院學(xué)報(bào);2011年12期
4 袁益讓;一類退化非線性拋物型方程組的變網(wǎng)格有限元方法[J];科學(xué)通報(bào);1985年15期
5 趙登虎,李志敏;有限元方法中網(wǎng)格編碼的優(yōu)化問題[J];工科數(shù)學(xué);2001年02期
6 劉震,李起升,白永強(qiáng);有限元方法的保結(jié)構(gòu)算法(英文)[J];河南科學(xué);2004年05期
7 陳樂生;;第一講 什么是有限元方法[J];木工機(jī)床;2006年02期
8 李宏;魏小溪;;奇異非線性拋物方程的時(shí)空有限元方法[J];高等學(xué)校計(jì)算數(shù)學(xué)學(xué)報(bào);2007年01期
9 孫繼華;趙洪賢;韓曉華;董欣;孟令華;李慶卓;黃緒萍;杜石巖;;基于有限元方法的凹槽超聲檢測(cè)[J];計(jì)測(cè)技術(shù);2009年06期
10 智晉寧;;有限元方法課程教學(xué)改革與實(shí)踐[J];安徽工業(yè)大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版);2010年05期
相關(guān)會(huì)議論文 前10條
1 許鶴華;周蒂;;非連續(xù)有限元方法的發(fā)展及其在地球科學(xué)中的應(yīng)用[A];第七屆全國(guó)數(shù)學(xué)地質(zhì)與地學(xué)信息學(xué)術(shù)會(huì)議論文摘要匯編[C];2004年
2 徐方遷;何世堂;;厚金屬柵力學(xué)負(fù)載貢獻(xiàn)反射系數(shù)的有限元方法[A];中國(guó)聲學(xué)學(xué)會(huì)2005年青年學(xué)術(shù)會(huì)議[CYCA'05]論文集[C];2005年
3 許鶴華;;連續(xù)時(shí)間有限元方法在求解非穩(wěn)態(tài)熱傳導(dǎo)的應(yīng)用[A];2000年中國(guó)地球物理學(xué)會(huì)年刊——中國(guó)地球物理學(xué)會(huì)第十六屆年會(huì)論文集[C];2000年
4 陳文;陳林;傅卓佳;;河海大學(xué)“工程與科學(xué)數(shù)值模擬軟件”的研究與開發(fā)[A];慶祝中國(guó)力學(xué)學(xué)會(huì)成立50周年暨中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)’2007論文摘要集(下)[C];2007年
5 曹雄;晉長(zhǎng)秋;;兩種有限元方法能量守恒分析[A];中國(guó)工程物理研究院科技年報(bào)(2000)[C];2000年
6 陳銳敏;;求解電磁位場(chǎng)的高階曲邊有限元方法[A];1987年全國(guó)微波會(huì)議論文集(上)[C];1987年
7 申文;馮西橋;;細(xì)胞粘附的有限元模擬[A];損傷、斷裂與微納米力學(xué)進(jìn)展:損傷、斷裂與微納米力學(xué)研討會(huì)論文集[C];2009年
8 蔚喜軍;符鴻源;常謙順;;用有限元方法求解雙曲守恒律[A];中國(guó)工程物理研究院科技年報(bào)(1998)[C];1998年
9 龍丹冰;劉西拉;;特大增量步算法在二維連續(xù)體分析上的拓展[A];中國(guó)計(jì)算力學(xué)大會(huì)'2010(CCCM2010)暨第八屆南方計(jì)算力學(xué)學(xué)術(shù)會(huì)議(SCCM8)論文集[C];2010年
10 隋永楓;;陀螺系統(tǒng)時(shí)間有限元的內(nèi)點(diǎn)法[A];中國(guó)計(jì)算力學(xué)大會(huì)'2010(CCCM2010)暨第八屆南方計(jì)算力學(xué)學(xué)術(shù)會(huì)議(SCCM8)論文集[C];2010年
相關(guān)博士學(xué)位論文 前10條
1 周振華;自適應(yīng)連續(xù)內(nèi)罰有限元方法和自適應(yīng)多罰間斷Galerkin方法[D];南京大學(xué);2014年
2 周少玲;非牛頓流體模型的最小二乘有限元方法[D];上海大學(xué);2015年
3 宋飛;間斷、組合多尺度有限元方法的分析與計(jì)算[D];南京大學(xué);2016年
4 紀(jì)海峰;求解帶有界面和不規(guī)則區(qū)域的橢圓問題的浸入界面有限元方法[D];南京師范大學(xué);2015年
5 錢凌志;流體及流體耦合問題的有限元方法研究[D];南京師范大學(xué);2016年
6 卞正寧;一階有限元方法研究以及圓柱渦振問題的數(shù)值模擬[D];湖南大學(xué);2016年
7 何斯日古楞;發(fā)展型方程的混合間斷時(shí)空有限元方法[D];內(nèi)蒙古大學(xué);2011年
8 王春梅;橢圓型偏微分方程的弱有限元方法研究[D];南京師范大學(xué);2014年
9 王奇生;幾類初邊值問題重疊型非匹配網(wǎng)格的有限元方法及收斂性分析[D];湘潭大學(xué);2007年
10 郭會(huì);幾類發(fā)展方程的最小二乘有限元方法[D];山東大學(xué);2006年
相關(guān)碩士學(xué)位論文 前10條
1 王懷志;航天器典型結(jié)構(gòu)中高頻動(dòng)力學(xué)環(huán)境預(yù)示的能量有限元方法[D];哈爾濱工業(yè)大學(xué);2015年
2 徐宇;基于有限元方法的心臟力學(xué)仿真[D];哈爾濱工業(yè)大學(xué);2015年
3 張雨晴;兩類流體力學(xué)方程組的兩重變分尺度有限元方法[D];溫州大學(xué);2015年
4 李生濤;PBX變形破壞的宏細(xì)觀數(shù)值模擬[D];北京理工大學(xué);2015年
5 錢雪;二維Sobolev方程的局部間斷Galerkin有限元方法[D];南京大學(xué);2014年
6 宋航;多尺度橢圓問題的粗細(xì)網(wǎng)格有限元方法[D];南京大學(xué);2013年
7 張煒;多尺度橢圓問題的間斷Petrov-Galerkin有限元方法[D];南京大學(xué);2014年
8 張瓊潔;橢圓界面問題近似的非匹配界面罰有限元方法[D];南京大學(xué);2014年
9 付海博;基于有限元方法電法測(cè)井模型的數(shù)值仿真及應(yīng)用[D];電子科技大學(xué);2015年
10 彭聰;時(shí)域有限元方法在仿真微波無(wú)源器件中的應(yīng)用[D];電子科技大學(xué);2015年
,本文編號(hào):2000324
本文鏈接:http://www.sikaile.net/shoufeilunwen/jckxbs/2000324.html