有限元分析的節(jié)點(diǎn)和單元(有限元分析的節(jié)點(diǎn)和單元有哪些)
有限元法是一種有效解決數(shù)學(xué)問題的解題方法。其基本求解思想是把計(jì)算域劃分為有限個(gè)互不重疊的單元,在每個(gè)單元內(nèi),選擇一些合適的節(jié)點(diǎn)作為求解函數(shù)的插值點(diǎn),單元上所作用的力等效到節(jié)點(diǎn)上,將微分方程中的變量改寫成由各變量或其導(dǎo)數(shù)的節(jié)點(diǎn)值與所選用的插值函數(shù)組成的線性表達(dá)式,就是用叉值函數(shù)來近似代替 ,借助于變分原理或加權(quán)余量法,將微分方程離散求解。
什么是有限元
有限元法是一種有效解決數(shù)學(xué)問題的解題方法。其基本求解思想是把計(jì)算域劃分為有限個(gè)互不重疊的單元,在每個(gè)單元內(nèi),選擇一些合適的節(jié)點(diǎn)作為求解函數(shù)的插值點(diǎn),單元上所作用的力等效到節(jié)點(diǎn)上,將微分方程中的變量改寫成由各變量或其導(dǎo)數(shù)的節(jié)點(diǎn)值與所選用的插值函數(shù)組成的線性表達(dá)式,就是用叉值函數(shù)來近似代替 ,借助于變分原理或加權(quán)余量法,將微分方程離散求解。
有限元中殼單元和實(shí)體單元有什么本質(zhì)區(qū)別
從本質(zhì)上來說,所有單元都是基于變分原理的,通過插值函數(shù)將節(jié)點(diǎn)位移與單元內(nèi)部應(yīng)變聯(lián)系都一起,再通過虛功原理或變分原理得到單元?jiǎng)偠扰c恢復(fù)力的計(jì)算公式。殼單元與實(shí)體元的節(jié)點(diǎn)數(shù)、節(jié)點(diǎn)自由度數(shù)、單元位移模式(插值函數(shù))不一樣,所以導(dǎo)致他們的單剛和恢復(fù)力表現(xiàn)形式不一樣,而這就是有限元實(shí)體單元和殼單元二者的本質(zhì)區(qū)別1,單元自由度不同;2,單元數(shù)量不同;3,基本理論不同,4,應(yīng)用對象也不同。1、單元自由度。有限元中實(shí)體單元只有三個(gè)自由度,而殼單元有六個(gè)自由度。
有限元方法與有限差分到底有什么區(qū)別
1.1 概念
有限差分方法(FDM)是計(jì)算機(jī)數(shù)值模擬最早采用的方法,至今仍被廣泛運(yùn)用。該方法將求解域劃分為差分網(wǎng)格,用有限個(gè)網(wǎng)格節(jié)點(diǎn)代替連續(xù)的求解域。有限差分法以Taylor級數(shù)展開等方法,把控制方程中的導(dǎo)數(shù)用網(wǎng)格節(jié)點(diǎn)上的函數(shù)值的差商代替進(jìn)行離散,從而建立以網(wǎng)格節(jié)點(diǎn)上的值為未知數(shù)的代數(shù)方程組。該方法是一種直接將微分問題變?yōu)榇鷶?shù)問題的近似數(shù)值解法,數(shù)學(xué)概念直觀,表達(dá)簡單,是發(fā)展較早且比較成熟的數(shù)值方法。
1.2 差分格式
(1)從格式的精度來劃分,有一階格式、二階格式和高階格式。
(2)從差分的空間形式來考慮,可分為中心格式和逆風(fēng)格式。
(3)考慮時(shí)間因子的影響,差分格式還可以分為顯格式、隱格式、顯隱交替格式等。
目前常見的差分格式,主要是上述幾種形式的組合,不同的組合構(gòu)成不同的差分格式。差分方法主要適用于有結(jié)構(gòu)網(wǎng)格,網(wǎng)格的步長一般根據(jù)實(shí)際地形的情況和柯朗穩(wěn)定條件來決定。
1.3 構(gòu)造差分的方法
構(gòu)造差分的方法有多種形式,目前主要采用的是泰勒級數(shù)展開方法。其基本的差分表達(dá)式主要有三種形式:一階向前差分、一階向后差分、一階中心差分和二階中心差分等,其中前兩種格式為一階計(jì)算精度,后兩種格式為二階計(jì)算精度。通過對時(shí)間和空間這幾種不同差分格式的組合,可以組合成不同的差分計(jì)算格式。
2. FEM
2.1 概述
有限元方法的基礎(chǔ)是變分原理和加權(quán)余量法,其基本求解思想是把計(jì)算域劃分為有限個(gè)互不重疊的單元,在每個(gè)單元內(nèi),選擇一些合適的節(jié)點(diǎn)作為求解函數(shù)的插值點(diǎn),將微分方程中的變量改寫成由各變量或其導(dǎo)數(shù)的節(jié)點(diǎn)值與所選用的插值函數(shù)組成的線性表達(dá)式,借助于變分原理或加權(quán)余量法,將微分方程離散求解。采用不同的權(quán)函數(shù)和插值函數(shù)形式,便構(gòu)成不同的有限元方法。
2.2 原理
有限元方法最早應(yīng)用于結(jié)構(gòu)力學(xué),后來隨著計(jì)算機(jī)的發(fā)展慢慢用于流體力學(xué)、土力學(xué)的數(shù)值模擬。在有限元方法中,把計(jì)算域離散剖分為有限個(gè)互不重疊且相互連接的單元,在每個(gè)單元內(nèi)選擇基函數(shù),用單元基函數(shù)的線形組合來逼近單元中的真解,整個(gè)計(jì)算域上總體的基函數(shù)可以看為由每個(gè)單元基函數(shù)組成的,則整個(gè)計(jì)算域內(nèi)的解可以看作是由所有單元上的近似解構(gòu)成。在河道數(shù)值模擬中,常見的有限元計(jì)算方法是由變分法和加權(quán)余量法發(fā)展而來的里茲法和伽遼金法、最小二乘法等。
根據(jù)所采用的權(quán)函數(shù)和插值函數(shù)的不同,有限元方法也分為多種計(jì)算格式。
(1)從權(quán)函數(shù)的選擇來說,有配置法、矩量法、最小二乘法和伽遼金法;
(2)從計(jì)算單元網(wǎng)格的形狀來劃分,有三角形網(wǎng)格、四邊形網(wǎng)格和多邊形網(wǎng)格;
(3)從插值函數(shù)的精度來劃分,又分為線性插值函數(shù)和高次插值函數(shù)等。
不同的組合同樣構(gòu)成不同的有限元計(jì)算格式。
對于權(quán)函數(shù),伽遼金(Galerkin)法是將權(quán)函數(shù)取為逼近函數(shù)中的基函數(shù);最小二乘法是令權(quán)函數(shù)等于余量本身,而內(nèi)積的極小值則為對代求系數(shù)的*方誤差最小;在配置法中,先在計(jì)算域內(nèi)選取N個(gè)配置點(diǎn)。令近似解在選定的N個(gè)配置點(diǎn)上嚴(yán)格滿足微分方程,即在配置點(diǎn)上令方程余量為0。插值函數(shù)一般由不同次冪的多項(xiàng)式組成,但也有采用三角函數(shù)或指數(shù)函數(shù)組成的乘積表示,但最常用的多項(xiàng)式插值函數(shù)。
有限元插值函數(shù)分為兩大類,一類只要求插值多項(xiàng)式本身在插值點(diǎn)取已知值,稱為拉格朗日(Lagrange)多項(xiàng)式插值;另一種不僅要求插值多項(xiàng)式本身,還要求它的導(dǎo)數(shù)值在插值點(diǎn)取已知值,稱為哈密特(Hermite)多項(xiàng)式插值。單元坐標(biāo)有笛卡爾直角坐標(biāo)系和無因次自然坐標(biāo),有對稱和不對稱等。常采用的無因次坐標(biāo)是一種局部坐標(biāo)系,它的定義取決于單元的幾何形狀,一維看作長度比,二維看作面積比,三維看作體積比。在二維有限元中,三角形單元應(yīng)用的最早,近來四邊形等參元的應(yīng)用也越來越廣。對于二維三角形和四邊形電源單元,常采用的插值函數(shù)為有Lagrange插值直角坐標(biāo)系中的線性插值函數(shù)及二階或更高階插值函數(shù)、面積坐標(biāo)系中的線性插值函數(shù)、二階或更高階插值函數(shù)等。
2.3 基本原理與解題步驟
對于有限元方法,其基本思路和解題步驟可歸納為:
(1)建立積分方程,根據(jù)變分原理或方程余量與權(quán)函數(shù)正交化原理,建立與微分方程初邊值問題等價(jià)的積分表達(dá)式,這是有限元法的出發(fā)點(diǎn)。
(2)區(qū)域單元剖分,根據(jù)求解區(qū)域的形狀及實(shí)際問題的物理特點(diǎn),將區(qū)域剖分為若干相互連接、不重疊的單元。區(qū)域單元?jiǎng)澐质遣捎糜邢拊椒ǖ那捌跍?zhǔn)備工作,這部分工作量比較大,除了給計(jì)算單元和節(jié)點(diǎn)進(jìn)行編號和確定相互之間的關(guān)系之外,還要表示節(jié)點(diǎn)的位置坐標(biāo),同時(shí)還需要列出自然邊界和本質(zhì)邊界的節(jié)點(diǎn)序號和相應(yīng)的邊界值。
(3)確定單元基函數(shù),根據(jù)單元中節(jié)點(diǎn)數(shù)目及對近似解精度的要求,選擇滿足一定插值條件的插值函數(shù)作為單元基函數(shù)。有限元方法中的基函數(shù)是在單元中選取的,由于各單元具有規(guī)則的幾何形狀,在選取基函數(shù)時(shí)可遵循一定的法則。
(4)單元分析:將各個(gè)單元中的求解函數(shù)用單元基函數(shù)的線性組合表達(dá)式進(jìn)行逼近;再將近似函數(shù)代入積分方程,并對單元區(qū)域進(jìn)行積分,可獲得含有待定系數(shù)(即單元中各節(jié)點(diǎn)的參數(shù)值)的代數(shù)方程組,稱為單元有限元方程。
(5)總體合成:在得出單元有限元方程之后,將區(qū)域中所有單元有限元方程按一定法則進(jìn)行累加,形成總體有限元方程。
(6)邊界條件的處理:一般邊界條件有三種形式,分為本質(zhì)邊界條件(狄里克雷邊界條件)、自然邊界條件(黎曼邊界條件)、混合邊界條件(柯西邊界條件)。對于自然邊界條件,一般在積分表達(dá)式中可自動(dòng)得到滿足。對于本質(zhì)邊界條件和混合邊界條件,需按一定法則對總體有限元方程進(jìn)行修正滿足。
(7)解有限元方程:根據(jù)邊界條件修正的總體有限元方程組,是含所有待定未知量的封閉方程組,采用適當(dāng)?shù)臄?shù)值計(jì)算方法求解,可求得各節(jié)點(diǎn)的函數(shù)值。
3. 有限體積法
有限體積法(FiniteVolumeMethod)又稱為控制體積法。其基本思路是:將計(jì)算區(qū)域劃分為一系列不重復(fù)的控制體積,并使每個(gè)網(wǎng)格點(diǎn)周圍有一個(gè)控制體積;將待解的微分方程對每一個(gè)控制體積積分,便得出一組離散方程。其中的未知數(shù)是網(wǎng)格點(diǎn)上的因變量的數(shù)值。為了求出控制體積的積分,必須假定值在網(wǎng)格點(diǎn)之間的變化規(guī)律,即假設(shè)值的分段的分布的分布剖面。從積分區(qū)域的選取方法看來,有限體積法屬于加權(quán)剩余法中的子區(qū)域法;從未知解的近似方法看來,有限體積法屬于采用局部近似的離散方法。簡言之,子區(qū)域法屬于有限體積發(fā)的基本方法。有限體積法的基本思路易于理解,并能得出直接的物理解釋。離散方程的物理意義,就是因變量在有限大小的控制體積中的守恒原理,如同微分方程表示因變量在無限小的控制體積中的守恒原理一樣。限體積法得出的離散方程,要求因變量的積分守恒對任意一組控制體積都得到滿足,對整個(gè)計(jì)算區(qū)域,自然也得到滿足。這是有限體積法吸引人的優(yōu)點(diǎn)。有一些離散方法,例如有限差分法,僅當(dāng)網(wǎng)格極其細(xì)密時(shí),離散方程才滿足積分守恒;而有限體積法即使在粗網(wǎng)格情況下,也顯示出準(zhǔn)確的積分守恒。就離散方法而言,有限體積法可視作有限單元法和有限差分法的中間物。有限單元法必須假定值在網(wǎng)格點(diǎn)之間的變化規(guī)律(既插值函數(shù)),并將其作為近似解。有限差分法只考慮網(wǎng)格點(diǎn)上的數(shù)值而不考慮值在網(wǎng)格點(diǎn)之間如何變化。有限體積法只尋求的結(jié)點(diǎn)值,這與有限差分法相類似;但有限體積法在尋求控制體積的積分時(shí),必須假定值在網(wǎng)格點(diǎn)之間的分布,這又與有限單元法相類似。在有限體積法中,插值函數(shù)只用于計(jì)算控制體積的積分,得出離散方程之后,便可忘掉插值函數(shù);如果需要的話,可以對微分方程中不同的項(xiàng)采取不同的插值函數(shù)。
4. 比較分析
有限差分法(FDM):直觀,理論成熟,精度可眩但是不規(guī)則區(qū)域處理繁瑣,雖然網(wǎng)格生成可以使FDM應(yīng)用于不規(guī)則區(qū)域,但是對區(qū)域的連續(xù)性等要求較嚴(yán)。使用FDM的好處在于易于編程,易于并行。
有限元方法(FEM):適合處理復(fù)雜區(qū)域,精度可眩缺憾在于內(nèi)存和計(jì)算量巨大。并行不如FDM和FVM直觀。不過FEM的并行是當(dāng)前和將來應(yīng)用的一個(gè)不錯(cuò)的方向。
有限容積法:適于流體計(jì)算,可以應(yīng)用于不規(guī)則網(wǎng)格,適于并行。但是精度基本上只能是二階了。FVM的優(yōu)勢正逐漸顯現(xiàn)出來,F(xiàn)VM在應(yīng)力應(yīng)變,高頻電磁場方面的特殊的優(yōu)點(diǎn)正在被人重視。
比較一下:
有限容積法和有限差分法:一個(gè)區(qū)別就是有限容積法的截差是不定的(跟取的相鄰點(diǎn)有關(guān),積分方法離散方程),而有限差分就可以直接知道截差(微分方法離散方程)。有限容積法和有限差分法最本質(zhì)的區(qū)別是,前者是根據(jù)積分方程推導(dǎo)出來的(即對每個(gè)控制體積分),后者直接根據(jù)微分方程推導(dǎo)出來,所以前者的精度不但取決于積分時(shí)的精度,還取決與對導(dǎo)數(shù)處理的精度,一般有限容積法總體的精度為二階,因?yàn)榉e分的精度限制,當(dāng)然有限容積法對于守恒型方程導(dǎo)出的離散方程可以保持守恒型;而后者直接由微分方程導(dǎo)出,不涉及積分過程,各種導(dǎo)數(shù)的微分借助Taylor展開,直接寫出離散方程,當(dāng)然不一定有守恒性,精度也和有限容積法不一樣,一般有限差分法可以使精度更高一些。
當(dāng)然二者有聯(lián)系,有時(shí)導(dǎo)出的形式一樣,但是概念上是不一樣的。
至于有限容積法和有限元相比,有限元在復(fù)雜區(qū)域的適應(yīng)性對有限容積是毫無優(yōu)勢可言的,至于有限容積的守恒性,物理概念明顯的這些特點(diǎn),有限元是沒有的。目前有限容積在精度方面與有限元法有些差距。
有限元方法比有限差分優(yōu)越的方面主要在能適應(yīng)不規(guī)則區(qū)域,但是這只是指的是傳統(tǒng)意義上的有限差分,現(xiàn)在發(fā)展的一些有限差分已經(jīng)能適應(yīng)不規(guī)則區(qū)域。對于橢圓型方程,如果區(qū)域規(guī)則,傳統(tǒng)有限差分和有限元都能解,在求解效率,這里主要指編程負(fù)責(zé)度和收斂快慢、內(nèi)存需要,肯定有限差分有優(yōu)勢。
有限元離散基本原理
對于不同物理性質(zhì)和數(shù)學(xué)模型的問題,有限元求解法的基本步驟是相同的,只是具體公式推導(dǎo)和運(yùn)算求解不同。有限元求解問題的基本步驟通常為:第一步:問題及求解域定義:根據(jù)實(shí)際問題近似確定求解域的物理性質(zhì)和幾何區(qū)域。第二步:求解域離散化:將求解域近似為具有不同有限大小和形狀且彼此相連的有限個(gè)單元組成的離散域,習(xí)慣上稱為有限元網(wǎng)絡(luò)劃分。顯然單元越?。ňW(wǎng)絡(luò)越細(xì))則離散域的近似程度越好,計(jì)算結(jié)果也越精確,但計(jì)算量及誤差都將增大,因此求解域的離散化是有限元法的核心技術(shù)之一。第三步:確定狀態(tài)變量及控制方法:一個(gè)具體的物理問題通常可以用一組包含問題狀態(tài)變量邊界條件的微分方程式表示,為適合有限元求解,通常將微分方程化為等價(jià)的泛函形式。第四步:單元推導(dǎo):對單元構(gòu)造一個(gè)適合的近似解,即推導(dǎo)有限單元的列式,其中包括選擇合理的單元坐標(biāo)系,建立單元試函數(shù),以某種方法給出單元各狀態(tài)變量的離散關(guān)系,從而形成單元矩陣(結(jié)構(gòu)力學(xué)中稱剛度陣或柔度陣)。為保證問題求解的收斂性,單元推導(dǎo)有許多原則要遵循。 對工程應(yīng)用而言,重要的是應(yīng)注意每一種單元的解題性能與約束。例如,單元形狀應(yīng)以規(guī)則為好,畸形時(shí)不僅精度低,而且有缺秩的危險(xiǎn),將導(dǎo)致無法求解。第五步:總裝求解:將單元總裝形成離散域的總矩陣方程(聯(lián)合方程組),反映對近似求解域的離散域的要求,即單元函數(shù)的連續(xù)性要滿足一定的連續(xù)條件??傃b是在相鄰單元結(jié)點(diǎn)進(jìn)行,狀態(tài)變量及其導(dǎo)數(shù)(可能的話)連續(xù)性建立在結(jié)點(diǎn)處。第六步:聯(lián)立方程組求解和結(jié)果解釋:有限元法最終導(dǎo)致聯(lián)立方程組。聯(lián)立方程組的求解可用直接法、選代法和隨機(jī)法。求解結(jié)果是單元結(jié)點(diǎn)處狀態(tài)變量的近似值。對于計(jì)算結(jié)果的質(zhì)量,將通過與設(shè)計(jì)準(zhǔn)則提供的允許值比較來評價(jià)并確定是否需要重復(fù)計(jì)算。
矩陣法和有限元法的區(qū)別
矩陣位移法只能用于分析具有已知單元節(jié)點(diǎn)力——單元節(jié)點(diǎn)位移關(guān)系的桿系結(jié)構(gòu),而不能分析非桿系的連續(xù)體結(jié)構(gòu)。兩者的區(qū)別主要表現(xiàn)為建立的基本原理不同以及求解的側(cè)重點(diǎn)不同。矩陣位移法側(cè)重于根據(jù)結(jié)構(gòu)的幾何性質(zhì)和彈性性質(zhì)建立整體剛度矩陣,根據(jù)結(jié)構(gòu)的受載情況形成整體荷載向量,求出各單元的桿端內(nèi)力,需要疊加固端力;而有限元法側(cè)重于根據(jù)已知外荷載的作用和結(jié)構(gòu)材料的性質(zhì)將結(jié)構(gòu)進(jìn)行離散化后求解,結(jié)點(diǎn)力或是結(jié)點(diǎn)應(yīng)力是直接對位移求導(dǎo)獲得的,丟失了固端力一項(xiàng),致使應(yīng)力的精度有所下降。通過對*面剛架結(jié)構(gòu)進(jìn)行求解,得出結(jié)構(gòu)內(nèi)力,而這一區(qū)別在算例結(jié)果上得到了驗(yàn)證。
什么是有限元法和有限差分法
有限元法,有限差分法和有限體積法的區(qū)別\x0d有限差分方法(FDM)是計(jì)算機(jī)數(shù)值模擬最早采用的方法,至今仍被廣泛運(yùn)用.該方法將求解域劃分為差分網(wǎng)格,用有限個(gè)網(wǎng)格節(jié)點(diǎn)代替連續(xù)的求解域.有限差分法以Taylor級數(shù)展開等方法,把控制方程中的導(dǎo)數(shù)用網(wǎng)格節(jié)點(diǎn)上的函數(shù)值的差商代替進(jìn)行離散,從而建立以網(wǎng)格節(jié)點(diǎn)上的值為未知數(shù)的代數(shù)方程組.該方法是一種直接將微分問題變?yōu)榇鷶?shù)問題的近似數(shù)值解法,數(shù)學(xué)概念直觀,表達(dá)簡單,是發(fā)展較早且比較成熟的數(shù)值方法. 對于有限差分格式,從格式的精度來劃分,有一階格式、二階格式和高階格式.從差分的空間形式來考慮,可分為中心格式和逆風(fēng)格式.考慮時(shí)間因子的影響,差分格式還可以分為顯格式、隱格式、顯隱交替格式等.目前常見的差分格式,主要是上述幾種形式的組合,不同的組合構(gòu)成不同的差分格式.差分方法主要適用于有結(jié)構(gòu)網(wǎng)格,網(wǎng)格的步長一般根據(jù)實(shí)際地形的情況和柯朗穩(wěn)定條件來決定.\x0d構(gòu)造差分的方法有多種形式,目前主要采用的是泰勒級數(shù)展開方法.其基本的差分表達(dá)式主要有三種形式:一階向前差分、一階向后差分、一階中心差分和二階中心差分等,其中前兩種格式為一階計(jì)算精度,后兩種格式為二階計(jì)算精度.通過對時(shí)間和空間這幾種不同差分格式的組合,可以組合成不同的差分計(jì)算格式.\x0d有限元方法的基礎(chǔ)是變分原理和加權(quán)余量法,其基本求解思想是把計(jì)算域劃分為有限個(gè)互不重疊的單元,在每個(gè)單元內(nèi),選擇一些合適的節(jié)點(diǎn)作為求解函數(shù)的插值點(diǎn),將微分方程中的變量改寫成由各變量或其導(dǎo)數(shù)的節(jié)點(diǎn)值與所選用的插值函數(shù)組成的線性表達(dá)式,借助于變分原理或加權(quán)余量法,將微分方程離散求解.采用不同的權(quán)函數(shù)和插值函數(shù)形式,便構(gòu)成不同的有限元方法.有限元方法最早應(yīng)用于結(jié)構(gòu)力學(xué),后來隨著計(jì)算機(jī)的發(fā)展慢慢用于流體力學(xué)的數(shù)值模擬.在有限元方法中,把計(jì)算域離散剖分為有限個(gè)互不重疊且相互連接的單元,在每個(gè)單元內(nèi)選擇基函數(shù),用單元基函數(shù)的線形組合來逼近單元中的真解,整個(gè)計(jì)算域上總體的基函數(shù)可以看為由每個(gè)單元基函數(shù)組成的,則整個(gè)計(jì)算域內(nèi)的解可以看作是由所有單元上的近似解構(gòu)成.在河道數(shù)值模擬中,常見的有限元計(jì)算方法是由變分法和加權(quán)余量法發(fā)展而來的里茲法和伽遼金法、最小二乘法等.根據(jù)所采用的權(quán)函數(shù)和插值函數(shù)的不同,有限元方法也分為多種計(jì)算格式.從權(quán)函數(shù)的選擇來說,有配置法、矩量法、最小二乘法和伽遼金法,從計(jì)算單元網(wǎng)格的形狀來劃分,有三角形網(wǎng)格、四邊形網(wǎng)格和多邊形網(wǎng)格,從插值函數(shù)的精度來劃分,又分為線性插值函數(shù)和高次插值函數(shù)等.不同的組合同樣構(gòu)成不同的有限元計(jì)算格式.對于權(quán)函數(shù),伽遼金(Galerkin)法是將權(quán)函數(shù)取為逼近函數(shù)中的基函數(shù);最小二乘法是令權(quán)函數(shù)等于余量本身,而內(nèi)積的極小值則為對代求系數(shù)的*方誤差最小;在配置法中,先在計(jì)算域內(nèi)選取N個(gè)配置點(diǎn).令近似解在選定的N個(gè)配置點(diǎn)上嚴(yán)格滿足微分方程,即在配置點(diǎn)上令方程余量為0.插值函數(shù)一般由不同次冪的多項(xiàng)式組成,但也有采用三角函數(shù)或指數(shù)函數(shù)組成的乘積表示,但最常用的多項(xiàng)式插值函數(shù).有限元插值函數(shù)分為兩大類,一類只要求插值多項(xiàng)式本身在插值點(diǎn)取已知值,稱為拉格朗日(Lagrange)多項(xiàng)式插值;另一種不僅要求插值多項(xiàng)式本身,還要求它的導(dǎo)數(shù)值在插值點(diǎn)取已知值,稱為哈密特(Hermite)多項(xiàng)式插值.單元坐標(biāo)有笛卡爾直角坐標(biāo)系和無因次自然坐標(biāo),有對稱和不對稱等.常采用的無因次坐標(biāo)是一種局部坐標(biāo)系,它的定義取決于單元的幾何形狀,一維看作長度比,二維看作面積比,三維看作體積比.在二維有限元中,三角形單元應(yīng)用的最早,近來四邊形等參元的應(yīng)用也越來越廣.對于二維三角形和四邊形電源單元,常采用的插值函數(shù)為有Lagrange插值直角坐標(biāo)系中的線性插值函數(shù)及二階或更高階插值函數(shù)、面積坐標(biāo)系中的線性插值函數(shù)、二階或更高階插值函數(shù)等.\x0d對于有限元方法,其基本思路和解題步驟可歸納為\x0d(1)建立積分方程,根據(jù)變分原理或方程余量與權(quán)函數(shù)正交化原理,建立與微分方程初邊值問題等價(jià)的積分表達(dá)式,這是有限元法的出發(fā)點(diǎn).\x0d(2)區(qū)域單元剖分,根據(jù)求解區(qū)域的形狀及實(shí)際問題的物理特點(diǎn),將區(qū)域剖分為若干相互連接、不重疊的單元.區(qū)域單元?jiǎng)澐质遣捎糜邢拊椒ǖ那捌跍?zhǔn)備工作,這部分工作量比較大,除了給計(jì)算單元和節(jié)點(diǎn)進(jìn)行編號和確定相互之間的關(guān)系之外,還要表示節(jié)點(diǎn)的位置坐標(biāo),同時(shí)還需要列出自然邊界和本質(zhì)邊界的節(jié)點(diǎn)序號和相應(yīng)的邊界值.\x0d(3)確定單元基函數(shù),根據(jù)單元中節(jié)點(diǎn)數(shù)目及對近似解精度的要求,選擇滿足一定插值條件的插值函數(shù)作為單元基函數(shù).有限元方法中的基函數(shù)是在單元中選取的,由于各單元具有規(guī)則的幾何形狀,在選取基函數(shù)時(shí)可遵循一定的法則.\x0d(4)單元分析:將各個(gè)單元中的求解函數(shù)用單元基函數(shù)的線性組合表達(dá)式進(jìn)行逼近;再將近似函數(shù)代入積分方程,并對單元區(qū)域進(jìn)行積分,可獲得含有待定系數(shù)(即單元中各節(jié)點(diǎn)的參數(shù)值)的代數(shù)方程組,稱為單元有限元方程.\x0d(5)總體合成:在得出單元有限元方程之后,將區(qū)域中所有單元有限元方程按一定法則進(jìn)行累加,形成總體有限元方程.\x0d(6)邊界條件的處理:一般邊界條件有三種形式,分為本質(zhì)邊界條件(狄里克雷邊界條件 )、自然邊界條件(黎曼邊界條件)、混合邊界條件(柯西邊界條件).對于自然邊界條件,一般在積分表達(dá)式中可自動(dòng)得到滿足.對于本質(zhì)邊界條件和混合邊界條件,需按一定法則對總體有限元方程進(jìn)行修正滿足.\x0d(7)解有限元方程:根據(jù)邊界條件修正的總體有限元方程組,是含所有待定未知量的封閉方程組,采用適當(dāng)?shù)臄?shù)值計(jì)算方法求解,可求得各節(jié)點(diǎn)的函數(shù)值.\x0d有限體積法(Finite Volume Method)又稱為控制體積法.其基本思路是:將計(jì)算區(qū)域劃分為一系列不重復(fù)的控制體積,并使每個(gè)網(wǎng)格點(diǎn)周圍有一個(gè)控制體積;將待解的微分方程對每一個(gè)控制體積積分,便得出一組離散方程.其中的未知數(shù)是網(wǎng)格點(diǎn)上的因變量的數(shù)值.為了求出控制體積的積分,必須假定值在網(wǎng)格點(diǎn)之間的變化規(guī)律,即假設(shè)值的分段的分布的分布剖面.從積分區(qū)域的選取方法看來,有限體積法屬于加權(quán)剩余法中的子區(qū)域法;從未知解的近似方法看來,有限體積法屬于采用局部近似的離散方法.簡言之,子區(qū)域法屬于有限體積發(fā)的基本方法.\x0d有限體積法的基本思路易于理解,并能得出直接的物理解釋.離散方程的物理意義,就是因變量在有限大小的控制體積中的守恒原理,如同微分方程表示因變量在無限小的控制體積中的守恒原理一樣.限體積法得出的離散方程,要求因變量的積分守恒對任意一組控制體積都得到滿足,對整個(gè)計(jì)算區(qū)域,自然也得到滿足.這是有限體積法吸引人的優(yōu)點(diǎn).有一些離散方法,例如有限差分法,僅當(dāng)網(wǎng)格極其細(xì)密時(shí),離散方程才滿足積分守恒;而有限體積法即使在粗網(wǎng)格情況下,也顯示出準(zhǔn)確的積分守恒.就離散方法而言,有限體積法可視作有限單元法和有限差分法的中間物.有限單元法必須假定值在網(wǎng)格點(diǎn)之間的變化規(guī)律(既插值函數(shù)),并將其作為近似解.有限差分法只考慮網(wǎng)格點(diǎn)上的數(shù)值而不考慮值在網(wǎng)格點(diǎn)之間如何變化.有限體積法只尋求的結(jié)點(diǎn)值,這與有限差分法相類似;但有限體積法在尋求控制體積的積分時(shí),必須假定值在網(wǎng)格點(diǎn)之間的分布,這又與有限單元法相類似.
有限元分析的節(jié)點(diǎn)和單元有限元分析的節(jié)點(diǎn)和單元(有限元分析的節(jié)點(diǎn)和單元有哪些)
發(fā)表評論
還沒有評論,來說兩句吧...