目前數(shù)值模擬中常用方法還是有限元和有限差分法,這樣的方法需要構(gòu)造網(wǎng)格或單元,給問題解決帶來了麻煩。為了克服這樣的困難,無(wú)網(wǎng)格方法逐漸發(fā)展起來,日益受到關(guān)注。該類方法節(jié)點(diǎn)布置靈活,無(wú)需網(wǎng)格信息,部分或徹底避免了網(wǎng)格的影響,有助于推進(jìn)數(shù)值模擬技術(shù)的進(jìn)步。在這些方法中,無(wú)網(wǎng)格局部Petrov-Galerkin（MLPG）方法將Petrov-Galerkin方法與無(wú)網(wǎng)格近似函數(shù)結(jié)合在一起,在形函數(shù)構(gòu)造和數(shù)值積分過程中均不需要網(wǎng)格,是純無(wú)網(wǎng)格方法。本文首先概述了無(wú)網(wǎng)格方法及MLPG法的發(fā)展歷史。介紹了做為MLPG法研究基礎(chǔ)的加權(quán)殘量法、用于構(gòu)造近似解函數(shù)的移動(dòng)最小二乘法及兩種施加本質(zhì)邊界條件的方法。其次針對(duì)常見的滲透系數(shù)為分片常數(shù)的地下水模擬問題,建立了求解該類問題的MLPG法與配點(diǎn)法耦合求解的方法。該方法依據(jù)滲透系數(shù)為分片常數(shù)的特征,基于MLPG法和配點(diǎn)法分區(qū)域建立求解方程組,在每個(gè)子區(qū)除去分界線的節(jié)點(diǎn)上建立MLPG方程組,在各子區(qū)分界線布置的節(jié)點(diǎn)上,應(yīng)用相容條件建立配點(diǎn)方程組,從而得到了求解該類問題水頭函數(shù)數(shù)值解的耦合方程組,給出了該方程組的求解算法。將計(jì)算結(jié)果與MLPG法以及邊界元法的結(jié)果... 

【文章來源】：遼寧師范大學(xué)遼寧省

【文章頁(yè)數(shù)】：45 頁(yè)

【學(xué)位級(jí)別】：碩士

【部分圖文】：

二維基函數(shù)Fig.2.1Thetwo-dimensionalbasisfunction

遼寧師范大學(xué)碩士研究生學(xué)位論文91()()()()().mhTiiiupapax,xxxxx(2.23)其中：x是計(jì)算點(diǎn)x的鄰域x內(nèi)的節(jié)點(diǎn)。12()[(),(),,()]Tmaxaxaxax，()iax為待解量，12()[(),(),,()]Tmpxpxpxpx為基函數(shù)矩陣，m為基函數(shù)的個(gè)數(shù)。式中()ipx常選用單項(xiàng)式基函數(shù)，可由Pascal三角得到，二維和三維情況見圖2.1和圖2.2。線性基和二次基寫成具體向量形式為：線性基：()(1,,),3,()Tpxxym二維空間()(1,,,),4,()Tpxxyzm三維空間二次基：22()(1,,,,,),6,()Tpxxyxxyym二維空間222()(1,,,,,,,,,),10.()Tpxxyzxxyyyzzxzm三維空間圖2.1二維基函數(shù)圖2.2三維基函數(shù)Fig.2.1Thetwo-dimensionalbasisfunctionFig.2.2Three-dimensionalbasisfunction在用多項(xiàng)式進(jìn)行整體擬合時(shí)，適當(dāng)增加多項(xiàng)式次數(shù)擬合效果會(huì)好些，但次數(shù)過高，會(huì)生成病態(tài)方程，影響擬合效果。而移動(dòng)最小二乘法采用局部擬合的方式，較低次數(shù)的多項(xiàng)式就能達(dá)到不錯(cuò)的擬合效果。為了實(shí)現(xiàn)局部擬合，在整個(gè)研究域內(nèi)的每個(gè)離散節(jié)點(diǎn)(1,2,,)IxIN處定義一個(gè)權(quán)函數(shù)()=(,)IIIwxwxx，I為節(jié)點(diǎn)Ix對(duì)應(yīng)的權(quán)函數(shù)影響域半徑，權(quán)函數(shù)Iw只在Ix附近的一個(gè)有限區(qū)域I(稱為權(quán)函數(shù)Iw的支撐域，也稱節(jié)點(diǎn)Ix的影響域)內(nèi)大于零，在該區(qū)域外均為零，即()Iwx為緊支函數(shù)，離散點(diǎn)Ix的值只與它附近有限區(qū)域I內(nèi)的點(diǎn)有關(guān)，與區(qū)域外的點(diǎn)無(wú)關(guān)。在MLS法中，為了達(dá)到最佳模擬效果，要求(2.23)式中系數(shù)()iax的選取使得局部近似函數(shù)()hux,x在計(jì)算點(diǎn)x的鄰域x(x的定義域)內(nèi)與函數(shù)u(x)的逼近誤差最校在x的定義域內(nèi)包含的節(jié)點(diǎn)Ix=x處近似函數(shù)()hux,x與精確解之差的加權(quán)平方和為：

基于MLPG法的地下水模擬問題研究183.1.2局部積分弱形式的MLPG法離散在整個(gè)研究區(qū)域內(nèi)及其邊界上布置N個(gè)節(jié)點(diǎn)jx，j1,2,,N。依據(jù)這些節(jié)點(diǎn)建立移動(dòng)最小二乘形式的近似水頭函數(shù)：1()()(),NhjjjHhxxhx(3.7)其中，x(x,y)，j為移動(dòng)最小二乘法中形函數(shù)，jh為節(jié)點(diǎn)jx對(duì)應(yīng)的水頭函數(shù)近似值。以每個(gè)節(jié)點(diǎn)ix(i1,2,,N.)為中心，建立一個(gè)圓形子域i，將(3.7)式的近似函數(shù)代入(3.6)式中，并取v為以節(jié)點(diǎn)ix為中心的檢驗(yàn)函數(shù)iv，使得代入后式子在每個(gè)節(jié)點(diǎn)ix對(duì)應(yīng)的子域i上成立，可得求解h的離散化線性方程組：111[ddd]iiiNjijijjiijjvvvvhxxyyn12dd,1,,.iiiiHvQviN(3.8)3.2非均質(zhì)穩(wěn)定流問題的MLPG法與配點(diǎn)法耦合方法3.2.1耦合方法的建立在非均質(zhì)介質(zhì)中的實(shí)際地下水計(jì)算過程中，導(dǎo)水系數(shù)T常取為分片常數(shù)，即將區(qū)域劃分為若干均質(zhì)的子區(qū)域，在每個(gè)子區(qū)內(nèi)的T為常數(shù)，不同子區(qū)的T取不同常數(shù).對(duì)每個(gè)子區(qū)來說，用子區(qū)布置的所有節(jié)點(diǎn)建立移動(dòng)最小二乘近似函數(shù)，在除去子區(qū)分界線上布置的節(jié)點(diǎn)外的其它節(jié)點(diǎn)上列出MLPG方程組，在各子區(qū)分界線布置的節(jié)點(diǎn)上利用相容條件建立配點(diǎn)方程組，最終將每個(gè)子區(qū)建立的MLPG方程組與配點(diǎn)方程組聯(lián)立，即可計(jì)算出整個(gè)區(qū)域所有節(jié)點(diǎn)上解的近似值。以兩個(gè)分區(qū)為例，如圖3.1所示圖3.1非均質(zhì)區(qū)圖Fig.3.1Heterogeneousmap依據(jù)(3.7)式在Ⅰ區(qū)和Ⅱ區(qū)分別建立水頭近似函數(shù)：

【參考文獻(xiàn)】：
期刊論文
[1]熱彈性動(dòng)力學(xué)耦合問題的插值型移動(dòng)最小二乘無(wú)網(wǎng)格法研究[J]. 王峰,鄭保敬,林皋,周宜紅,范勇.  工程力學(xué). 2019(04)
[2]無(wú)網(wǎng)格全局介點(diǎn)法[J]. 楊建軍,鄭健龍.  應(yīng)用力學(xué)學(xué)報(bào). 2017(05)
[3]求解大Pe對(duì)流擴(kuò)散問題的MLPG新方法[J]. 陳征驥,謝文禮,李增耀,吳學(xué)紅.  工程熱物理學(xué)報(bào). 2016(08)
[4]帶源參數(shù)熱傳導(dǎo)問題的基于滑動(dòng)Kriging插值的MLPG法[J]. 王峰,林皋,鄭保敬,劉俊,李建波.  力學(xué)季刊. 2013(02)
[5]MLPG混合配點(diǎn)法在形狀優(yōu)化中的應(yīng)用研究[J]. 趙亮,李書,魯大偉.  計(jì)算力學(xué)學(xué)報(bào). 2011(01)
[6]彈性力學(xué)的插值型邊界無(wú)單元法[J]. 任紅萍,程玉民,張武.  中國(guó)科學(xué):物理學(xué) 力學(xué) 天文學(xué). 2010(03)
[7]大變形問題分析的局部Petrov-Galerkin法[J]. 熊淵博,崔洪雪,龍述堯.  計(jì)算力學(xué)學(xué)報(bào). 2009(03)
[8]無(wú)網(wǎng)格局部Petrov-Galerkin方法中本質(zhì)邊界條件的處理[J]. 趙美玲,聶玉峰,袁占斌.  應(yīng)用力學(xué)學(xué)報(bào). 2006(03)
[9]彈性力學(xué)的重構(gòu)核粒子邊界無(wú)單元法[J]. 秦義校,程玉民.  物理學(xué)報(bào). 2006(07)
[10]彈塑性力學(xué)問題的無(wú)網(wǎng)格法分析[J]. 熊淵博,龍述堯,劉凱遠(yuǎn).  機(jī)械強(qiáng)度. 2004(06)

碩士論文
[1]MLPG方法在地下水穩(wěn)定流計(jì)算中的應(yīng)用[D]. 付鑫.遼寧師范大學(xué) 2014
[2]用無(wú)網(wǎng)格局部彼得羅夫—伽遼金法求解非均質(zhì)多孔介質(zhì)中的水流問題[D]. 吳志紅.遼寧師范大學(xué) 2009



本文編號(hào)：3625578