【摘要】：無網(wǎng)格(meshless method)法作為近十年來迅速興起的一種數(shù)值分析方法,它以其新穎的數(shù)值思想和先進(jìn)的數(shù)值技術(shù),得到了學(xué)術(shù)界的初步認(rèn)可和廣泛關(guān)注。不同于有限元法,無網(wǎng)格法在構(gòu)造形函數(shù)時(shí),是建立在一系列離散的點(diǎn)上,不需要借助網(wǎng)格,消除了節(jié)點(diǎn)之間的網(wǎng)格束縛,減少了因網(wǎng)格畸變而帶來的困難,從而擺脫了有限元繁雜的單元網(wǎng)格生成過程。在分析裂紋等不連續(xù)問題時(shí),無網(wǎng)格法具有獨(dú)特的優(yōu)勢(shì),其中無網(wǎng)格伽遼金法(EFG)法,以其精度高、穩(wěn)定性好和收斂速度快等特點(diǎn)成為最具發(fā)展前景的一種,本文采用這種方法對(duì)線性問題和斷裂問題進(jìn)行了研究。(1)以線彈性力斷裂學(xué)為分析手段,總結(jié)并概括了連續(xù)介質(zhì)和非連續(xù)介質(zhì)的斷裂數(shù)值方法,闡述了國內(nèi)外無網(wǎng)格法的研究進(jìn)展。(2)基于加權(quán)殘量法基本理論,論述了無網(wǎng)格法的幾種常見的插值方法并分析了各自的優(yōu)缺點(diǎn),對(duì)插值方法中移動(dòng)最小二乘法(MLS)進(jìn)行了重點(diǎn)的闡述。在此基礎(chǔ)上,對(duì)形函數(shù)里面的矩陣存在奇異性情況,給出了解決方案,對(duì)不同種類的權(quán)函數(shù)進(jìn)行數(shù)值計(jì)算和分析。(3)根據(jù)線彈性斷裂力學(xué),從應(yīng)力強(qiáng)度因子和能量釋放率兩個(gè)方面對(duì)裂紋體的力學(xué)行為進(jìn)行描述,基于無網(wǎng)格伽遼金法(EFG)并對(duì)其進(jìn)行改進(jìn),分析了線彈性斷裂力學(xué)中應(yīng)力強(qiáng)度因子的數(shù)值求法并在程序中成功實(shí)現(xiàn)。(4)對(duì)EFG法中處理不連續(xù)問題的四大類數(shù)值方法進(jìn)行了系統(tǒng)而全面的總結(jié)和分析,這些數(shù)值方法為不連續(xù)問題的選取準(zhǔn)則提供了參考依據(jù)。本文的創(chuàng)新在于對(duì)擴(kuò)充形函數(shù)進(jìn)行分析,改進(jìn)了EFG法并建立了基于單位分解思想的無網(wǎng)格伽遼金法的近似位移函數(shù),詳細(xì)闡述程序設(shè)計(jì)中的一些關(guān)鍵問題,制定程序流程圖,并給出相應(yīng)的主程序。(5)基于MATLAB軟件,編寫了改進(jìn)無網(wǎng)格伽遼金法計(jì)算程序,計(jì)算了MLS法近似函數(shù),有效的分析了線彈性問題即一維桿件的計(jì)算和二維懸臂深梁的求解,驗(yàn)證了形函數(shù)的性質(zhì)及參數(shù)的變化規(guī)律;并通過程序計(jì)算了斷裂問題中的應(yīng)力集中現(xiàn)象及應(yīng)力強(qiáng)度因子的計(jì)算,與擴(kuò)展有限元法作對(duì)比得到了較高的近似精度,討論了權(quán)函數(shù)、節(jié)點(diǎn)分布密度和影響域擴(kuò)大系數(shù)對(duì)求解精度的影響及變化規(guī)律,通過算例驗(yàn)證本文編寫計(jì)算程序的準(zhǔn)確性。
[Abstract]:Meshless method (meshless) method, as a numerical analysis method which has been rising rapidly in the past ten years, has gained the initial recognition and widespread concern in the academic circle with its new numerical ideas and advanced numerical techniques. Grid, which eliminates the constraint of grid between nodes, reduces the difficulty caused by grid distortion, and gets rid of the complex element mesh generation process of finite element. In the analysis of the discontinuity problems such as cracks, the meshless method has unique advantages, in which the EFG method is high in precision, good in stability and fast in convergence. This method is the most promising one. In this paper, this method is used to study linear problems and fracture problems. (1) the fracture numerical method of linear elastic force is used to sum up and summarize the fracture numerical methods of continuous medium and discontinuous medium, and the research progress of the meshless method at home and abroad is expounded. (2) based on the weighted residual method In this theory, several common interpolation methods of meshless method are discussed and their advantages and disadvantages are analyzed, and the moving least square method (MLS) in the interpolation method is emphatically expounded. On this basis, the singular situation of the matrix in the shape function is given, and the solution is given, and the numerical calculation and classification of different kinds of weight functions are given. (3) according to the linear elastic fracture mechanics, the mechanical behavior of the cracked body is described from two aspects of the stress intensity factor and the energy release rate. Based on the meshless Galerkin method (EFG) and improving it, the numerical method for the stress intensity factor in the linear elastic fracture mechanics is analyzed and successfully realized in the program. (4) the treatment of the EFG method is not good. Four kinds of numerical methods for continuous problems are systematically and comprehensively summarized and analyzed. These numerical methods provide reference for the selection criteria of discontinuous problems. The innovation of this paper is to analyze the extended shape function, improve the EFG method and establish the approximate displacement function of the element free Galerkin method based on the unit decomposition idea. Several key problems in the program design are elaborated in detail, the program flow chart is formulated and the corresponding main program is given. (5) based on the MATLAB software, an improved grid free Galerkin method is compiled and the approximate function of the MLS method is calculated. The linear elastic problem, that is, the calculation of the linear elastic problem, and the solution of the two-dimensional cantilever deep beam, is effectively analyzed. The properties of the shape function and the change law of the parameters, and the calculation of the stress concentration and the stress intensity factor in the fracture problem are calculated by the program. The higher approximate accuracy is obtained by comparing with the extended finite element method. The influence of the weight function, the distribution density of the node and the enlargement coefficient of the influence domain on the solution accuracy are discussed. An example is given to verify the accuracy of the program written in this paper.
【學(xué)位授予單位】：河南理工大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2014
【分類號(hào)】：TU312.3

【相似文獻(xiàn)】

相關(guān)期刊論文 前10條

1 劉加光;陳義保;羅震;;無網(wǎng)格伽遼金法在二維結(jié)構(gòu)問題中的應(yīng)用研究[J];山東理工大學(xué)學(xué)報(bào)(自然科學(xué)版);2006年03期

2 孟廣偉;周立明;李鋒;沙麗榮;;攝動(dòng)隨機(jī)局部正交無網(wǎng)格伽遼金法[J];吉林大學(xué)學(xué)報(bào)(工學(xué)版);2010年06期

3 于玲;夏茂輝;楊紅艷;翟社霞;;隨機(jī)無網(wǎng)格伽遼金法在疲勞斷裂可靠性分析中的應(yīng)用[J];山東理工大學(xué)學(xué)報(bào)(自然科學(xué)版);2010年05期

4 孟廣偉;趙云亮;周立明;李鋒;;無網(wǎng)格伽遼金法中兩種基函數(shù)的性質(zhì)[J];吉林大學(xué)學(xué)報(bào)(工學(xué)版);2009年03期

5 張亞靜;夏茂輝;張文婧;;無網(wǎng)格伽遼金法在板彎曲問題中的應(yīng)用[J];佳木斯大學(xué)學(xué)報(bào)(自然科學(xué)版);2008年06期

6 周立明;孟廣偉;劉昕暉;溫泳;;含裂紋結(jié)構(gòu)的模糊隨機(jī)無網(wǎng)格伽遼金法研究[J];華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版);2012年01期

7 湯紅衛(wèi),杜燁,王衛(wèi)東;應(yīng)用無網(wǎng)格伽遼金法模擬疲勞裂紋擴(kuò)展問題[J];山東大學(xué)學(xué)報(bào)(工學(xué)版);2005年02期

8 孟廣偉;周立明;李鋒;彭惠芬;沙麗榮;;局部正交無網(wǎng)格伽遼金法的研究及其在含多裂紋多孔結(jié)構(gòu)中的應(yīng)用[J];計(jì)算力學(xué)學(xué)報(bào);2010年06期

9 游世輝;曾德長;陳子娟;;用Taylor展開隨機(jī)無網(wǎng)格法和免疫算法分析結(jié)構(gòu)可靠性[J];機(jī)械強(qiáng)度;2008年02期

10 介玉新;劉巖;;基于蒙特卡羅積分的無網(wǎng)格伽遼金法及其在滲流分析中的應(yīng)用[J];巖土工程學(xué)報(bào);2007年12期

相關(guān)會(huì)議論文 前7條

1 殷水平;羅迎社;余敏;;無網(wǎng)格伽遼金法在金屬流變成形中的應(yīng)用[A];慶祝中國力學(xué)學(xué)會(huì)成立50周年暨中國力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)’2007論文摘要集（下）[C];2007年

2 陳莘莘;劉應(yīng)華;岑章志;;安定分析的無網(wǎng)格伽遼金法[A];北京力學(xué)會(huì)第13屆學(xué)術(shù)年會(huì)論文集[C];2007年

3 羅迎社;殷水平;余敏;;無網(wǎng)格伽遼金法在金屬流變成形中的應(yīng)用[A];中國化學(xué)會(huì)、中國力學(xué)學(xué)會(huì)第九屆全國流變學(xué)學(xué)術(shù)會(huì)議論文摘要集[C];2008年

4 趙光明;宋順成;;具有Kronecher delta插值性質(zhì)的改進(jìn)無網(wǎng)格伽遼金法[A];中國計(jì)算力學(xué)大會(huì)'2010（CCCM2010）暨第八屆南方計(jì)算力學(xué)學(xué)術(shù)會(huì)議（SCCM8）論文集[C];2010年

5 俞集輝;賈晉;;無網(wǎng)格伽遼金法在二維與三維非均勻電磁場(chǎng)問題中的應(yīng)用研究[A];第十一屆全國電工數(shù)學(xué)學(xué)術(shù)年會(huì)論文集[C];2007年

6 彭磊;鄧彩霞;;高斯-無網(wǎng)格伽遼金法解偏微分方程[A];中國力學(xué)大會(huì)——2013論文摘要集[C];2013年

7 華蕾娜;余錫平;;無網(wǎng)格伽遼金法在波浪作用下海床動(dòng)力響應(yīng)問題中的應(yīng)用[A];第二十一屆全國水動(dòng)力學(xué)研討會(huì)暨第八屆全國水動(dòng)力學(xué)學(xué)術(shù)會(huì)議暨兩岸船舶與海洋工程水動(dòng)力學(xué)研討會(huì)文集[C];2008年

相關(guān)博士學(xué)位論文 前3條

1 王暉;含裂紋功能復(fù)合材料結(jié)構(gòu)的單位分解擴(kuò)展無網(wǎng)格伽遼金法研究[D];吉林大學(xué);2016年

2 周立明;不確定性結(jié)構(gòu)的無網(wǎng)格伽遼金法研究[D];吉林大學(xué);2011年

3 趙光明;無網(wǎng)格方法及其在沖擊動(dòng)力學(xué)中的應(yīng)用研究[D];西南交通大學(xué);2006年

相關(guān)碩士學(xué)位論文 前10條

1 司建輝;無網(wǎng)格伽遼金法及應(yīng)用研究[D];西安理工大學(xué);2004年

2 蔣博;基于無網(wǎng)格伽遼金法（EFG）在線彈性斷裂分析中的求解與應(yīng)用[D];河南理工大學(xué);2014年

3 張亞靜;無網(wǎng)格伽遼金法的應(yīng)用研究[D];燕山大學(xué);2009年

4 姬威;若干邊界條件下的無網(wǎng)格伽遼金法及其改進(jìn)[D];燕山大學(xué);2010年

5 孫慧;無網(wǎng)格伽遼金法在非均質(zhì)多孔介質(zhì)水流問題中的應(yīng)用[D];遼寧師范大學(xué);2009年

6 李海龍;無網(wǎng)格伽遼金法精度分析與本質(zhì)邊界條件的處理[D];燕山大學(xué);2010年

7 于玲;隨機(jī)無網(wǎng)格伽遼金法在若干工程問題中的應(yīng)用研究[D];燕山大學(xué);2010年

8 趙玉鳳;改進(jìn)型無網(wǎng)格伽遼金法在溫度場(chǎng)的研究及應(yīng)用[D];燕山大學(xué);2012年

9 翟允;無網(wǎng)格伽遼金法計(jì)算精度的研究[D];燕山大學(xué);2009年

10 彭磊;無網(wǎng)格伽遼金法解帶控制的微分方程[D];哈爾濱理工大學(xué);2014年

，

本文編號(hào)：2157366