有限元分析基礎(chǔ)

武漢理工大學(xué)教務(wù)處制PAGE24有限元分析基礎(chǔ)有限元法概述在機(jī)械設(shè)計(jì)中，人們常常運(yùn)用材料力學(xué)、結(jié)構(gòu)力學(xué)等理論知識(shí)分析機(jī)械零構(gòu)件的強(qiáng)度、剛度和穩(wěn)定性問題。但對(duì)一些復(fù)雜的零構(gòu)件，這種分析常常就必須對(duì)其受力狀態(tài)和邊界條件進(jìn)行簡(jiǎn)化。否則力學(xué)分析將無法進(jìn)行。但這種簡(jiǎn)化的處理常常導(dǎo)致計(jì)算結(jié)果與實(shí)際相差甚遠(yuǎn)，有時(shí)甚至失去了分析的意義。所以過去設(shè)計(jì)經(jīng)驗(yàn)和類比占有較大比重。因?yàn)檫@個(gè)原因，人們也常常在設(shè)計(jì)中選擇較大的安全系數(shù)。如此也就造成所設(shè)計(jì)的機(jī)械結(jié)構(gòu)整體尺寸和重量偏大，而局部薄弱環(huán)節(jié)強(qiáng)度和剛度又不足的設(shè)計(jì)缺陷。近年來，數(shù)值計(jì)算機(jī)在工程分析上的成功運(yùn)用，產(chǎn)生了一門全新、高效的工程計(jì)算分析學(xué)科——有限元分析方法。該方法徹底改變了傳統(tǒng)工程分析中的做法。使計(jì)算精度和計(jì)算領(lǐng)域大大改善?！?.1有限元方法的發(fā)展歷史、現(xiàn)狀和將來歷史有限元法的起源應(yīng)追溯到上世紀(jì)40年代（20世紀(jì)40年代）。1943年R.Courant從數(shù)學(xué)的角度提出了有限元法的基本觀點(diǎn)。50年代中期在對(duì)飛機(jī)結(jié)構(gòu)的分析中，誕生了結(jié)構(gòu)分析的矩陣方法。1960年R.W.Clough在分析彈性力學(xué)平面問題時(shí)引入了“FiniteElementMethod”這一術(shù)語，從而標(biāo)志著有限元法的思想在力學(xué)分析中的廣泛推廣。60、70年代計(jì)算機(jī)技術(shù)的發(fā)展，極大地促進(jìn)了有限元法的發(fā)展。具體表現(xiàn)在：1）由彈性力學(xué)的平面問題擴(kuò)展到空間、板殼問題。2）由靜力平衡問題——穩(wěn)定性和動(dòng)力學(xué)分析問題。3）由彈性問題——彈塑性、粘彈性等問題。現(xiàn)狀現(xiàn)在有限元分析法的應(yīng)用領(lǐng)域已經(jīng)由開始時(shí)的固體力學(xué)，擴(kuò)展到流體力學(xué)、傳熱學(xué)和電磁力學(xué)等多個(gè)傳統(tǒng)的領(lǐng)域。已經(jīng)形成了一種非常成熟的數(shù)值分析計(jì)算方法。大型的商業(yè)化有限元分析軟件也是層出不窮，如：SAP系列的代表SAP2021（StructureAnalysisProgram）美國(guó)安世軟件公司的ANSYS大型綜合有限元分析軟件美國(guó)航天航空局的NASTRAN系列軟件除此以外，還有MASTER、ALGO、ABIQUES、ADINA、COSMOS等。將來有限元的發(fā)展方向最終將和CAD的發(fā)展相結(jié)合。運(yùn)用“四個(gè)化”可以概括其今后的發(fā)展趨勢(shì)。那就是：可視化、集成化、自動(dòng)化和網(wǎng)絡(luò)化?！?.2有限元法的特點(diǎn)機(jī)械零構(gòu)件的受力分析方法總體說來分為解析法和數(shù)值法兩大類。如大家學(xué)過的材料力學(xué)、結(jié)構(gòu)力學(xué)等就是經(jīng)典的解析力學(xué)分析方法。在這些解析力學(xué)方法中，彈性力學(xué)的分析方法在數(shù)學(xué)理論上是最為嚴(yán)謹(jǐn)?shù)囊环N分析方法。其解題思路是：從靜力、幾何和物理三個(gè)方面綜合考慮，建立描述彈性體的平衡、應(yīng)力、應(yīng)變和位移三者之間的微分方程，然后考慮邊界條件，從而求出微分方程的解析解。其最大的有點(diǎn)就是，嚴(yán)密精確。缺點(diǎn)就是微分方程的求解困難，很多情況下，無法求解。數(shù)值方法是一種近似的計(jì)算方法。具體又分為“有限差分法”和“有限元法”?！坝邢薏罘址ā笔菍⒌玫降奈⒎址匠屉x散成近似的差分方程。通過對(duì)一系列離散的差分方程求解，得到最終的力學(xué)問題近似解。其優(yōu)點(diǎn)就是：計(jì)算簡(jiǎn)單收斂性好。缺點(diǎn)是：計(jì)算程序無法標(biāo)準(zhǔn)化，在不能獲得整個(gè)問題的微分方程時(shí)，該方法不能運(yùn)用。由于其是將微分方程轉(zhuǎn)為差分方程，所以它是一種數(shù)學(xué)近似?！坝邢拊ā钡幕舅枷刖褪恰跋确趾蠛稀被蛘摺盎麨榱悖址e零為整”。與有限差分不同，它是在力學(xué)模型上進(jìn)行近似處理，也就是（分塊近似）。具體做法：把連續(xù)體模型轉(zhuǎn)為由有限個(gè)單元組成的離散體模型，離散體模型之間通過一些節(jié)點(diǎn)聯(lián)系。對(duì)于每一個(gè)離散體個(gè)體選擇簡(jiǎn)單的函數(shù)近似表示其中的物理變化規(guī)律（如位移等），運(yùn)用力學(xué)方程推導(dǎo)單元的平衡方程組，然后集合所有的方程組形成表征整體結(jié)構(gòu)的方程組，引入邊界條件，求取最后問題的解。優(yōu)點(diǎn)：概念清晰、易于學(xué)習(xí)理解，適用性強(qiáng)，便于電算化。缺點(diǎn)：計(jì)算精度受單元?jiǎng)澐值挠绊戄^大?！?.3有限元分析的一般過程為了能夠了解有限元分析的全貌，我們就一個(gè)簡(jiǎn)單的例子，來分析一下有限元分析的三個(gè)過程：結(jié)構(gòu)離散化、單元分析、整體分析。結(jié)構(gòu)離散化在該階段中，要完成把連續(xù)結(jié)構(gòu)的力學(xué)模型轉(zhuǎn)變?yōu)殡x散的力學(xué)模型。處理的好壞，直接影響到最后分析結(jié)果的正確與否、計(jì)算的精度和計(jì)算的效率。根據(jù)模型的傳力特性和分析的目標(biāo)，正確選擇單元類型。通常單元分為：一維單元、二維單元和三維單元。所謂一維單元就是指所求物理量?jī)H隨一個(gè)坐標(biāo)變量而變化的單元。如桁架、平面剛架和空間剛架單元。一維單元：桿單元、梁?jiǎn)卧?。二維單元：三角形單元、四邊形單元（平面類問題）三維單元：四面體單元、六面體單元等（空間問題）計(jì)算精度和計(jì)算效率：取決于單元?jiǎng)澐值男螤?、大小和分布狀況。通常單元愈多、愈密集，計(jì)算精度愈高，但計(jì)算效率愈低。有限元分析工作就是要在精度和效率兩者之間做到有機(jī)的統(tǒng)一。單元分析進(jìn)行單元分析的目的是為了到處表征單元力學(xué)特性的“單元?jiǎng)偠染仃嚒?。一般說來該過程有三種方法：直接法。虛功原理法（變分法）。加權(quán)余數(shù)法。直接法概念淺顯，易于理解物理含義。變分法需要泛函的數(shù)學(xué)知識(shí)，其推導(dǎo)過程具有嚴(yán)謹(jǐn)?shù)臄?shù)學(xué)概念。加權(quán)余數(shù)法適用于泛函不存在的應(yīng)用范圍。本教材將運(yùn)用虛功原理方程結(jié)合彈性力學(xué)和材料力學(xué)中的知識(shí)來推求幾種常見單元的單剛計(jì)算公式。現(xiàn)在先看一個(gè)簡(jiǎn)單的階梯軸的軸向拉伸問題例：如圖所示的變截面直桿，受拉力P，運(yùn)用有限元方法分析其變形。取任意單元，長(zhǎng)度為l，面積為A，單元內(nèi)任意一點(diǎn)的軸向位移和其位置坐標(biāo)成正比，即u=a0+a1x其中a0、a1為待定系數(shù)。由于桿的兩個(gè)端點(diǎn)節(jié)點(diǎn)1、2是單元上的點(diǎn)，所以它們應(yīng)該滿足上述方程。節(jié)點(diǎn)1，x1=0，∴u1=a0+a1×0=a0節(jié)點(diǎn)2，x2=l，∴u2=a0+a1×la1=(u2-u1)/l將求出的結(jié)果帶入方程并整理，就得：式中：N1、N2是形函數(shù)[N]形函數(shù)矩陣{δ}e節(jié)點(diǎn)位移向量由位移與應(yīng)變的關(guān)系知道：將上面推出的位移表達(dá)式代入，可得：上式中的[B]稱為應(yīng)變矩陣或幾何矩陣。運(yùn)用材料力學(xué)中的虎克定律，可以將應(yīng)變和應(yīng)力聯(lián)系起來。單向應(yīng)力狀態(tài)的虎克定律為××其中[S]稱為應(yīng)力矩陣。利用虛功方程可以建立力與位移之間的關(guān)系，也就是單剛方程。在后面我們將會(huì)推導(dǎo)出它的一般形式如下：式中：{F}e為單元節(jié)點(diǎn)力向量，對(duì)我們這個(gè)例子應(yīng)為[U1U2]T。[K]e為單元?jiǎng)偠染仃?。后面將推?dǎo)出它的計(jì)算公式為[D]矩陣是彈性矩陣。對(duì)于一維單元來說，就是E。所以我們這兒討論的例題：求得單剛矩陣，也就完成了單元分析?？偨Y(jié)單剛矩陣推導(dǎo)的步驟，應(yīng)該分為四步：假定單元內(nèi)位移變化的近似規(guī)律，即選擇位移模式。運(yùn)用幾何關(guān)系，推求位移與應(yīng)變的關(guān)系。應(yīng)用物理規(guī)律，把應(yīng)變與應(yīng)力聯(lián)系起來。運(yùn)用虛功方程的力與位移關(guān)系，求出單剛矩陣。單元分析是整個(gè)有限元分析的核心。不同的單元因?yàn)槠淞W(xué)特性不同，而具有不同的單元?jiǎng)偠染仃?，我們這本教材就是要學(xué)習(xí)幾種常用單元矩陣的推導(dǎo)和計(jì)算。了解各種單元的力學(xué)特性，為以后選擇單元類型打好基礎(chǔ)。整體分析由各單元?jiǎng)偠染仃嚱M集成整個(gè)結(jié)構(gòu)的總剛度矩陣。整個(gè)結(jié)構(gòu)有三個(gè)節(jié)點(diǎn)，首先將單元?jiǎng)偠染仃嚁U(kuò)充為3X3的矩陣，移動(dòng)各元素使之與單剛矩陣中的元素位置相對(duì)應(yīng)，如下：然后直接相加。把各單元的節(jié)點(diǎn)力向量組集成總的節(jié)點(diǎn)載荷向量。根據(jù)邊界條件，修改總剛度矩陣，獲得總剛方程組。邊界條件修改之前的總剛方程：修改以后（采用置“0，1”法）求解方程組，得出總的節(jié)點(diǎn)位移向量。解得的解是：有了節(jié)點(diǎn)位移，再回代到前面單元推導(dǎo)過程中的公式×和××，就可以求得每個(gè)單元的應(yīng)變和應(yīng)力了。從這個(gè)簡(jiǎn)單的例子，我們了解了有限元法求解力學(xué)問題三大步驟中的內(nèi)容，想必很多同學(xué)會(huì)說，這樣復(fù)雜，如果運(yùn)用材料力學(xué)的知識(shí)，我還來得快些。但是大家不要忘記，有限元的計(jì)算很多都是編程完成，而且現(xiàn)在很多的商業(yè)軟件都已經(jīng)完成了很多的工作。我們學(xué)習(xí)有限元主要是了解它的原理，并對(duì)常見單元的力學(xué)特性有所了解，這樣對(duì)于以后運(yùn)用有限元起到幫助作用。所以下面章節(jié)的內(nèi)容，就是圍繞這個(gè)主題展開。要達(dá)到這個(gè)目的，我們還必須學(xué)習(xí)必要的彈性力學(xué)知識(shí)。對(duì)彈性力學(xué)知識(shí)的學(xué)習(xí)，也對(duì)我們以后把握問題的本質(zhì)有幫助。

平面問題平面問題在力學(xué)研究的課題中屬?gòu)椥粤W(xué)的范疇。該類問題不僅本身具有典型性，而且在機(jī)械零構(gòu)件的分析中，也是應(yīng)用得非常廣泛。所以這類問題也稱之為經(jīng)典的力學(xué)問題。我們知道，實(shí)際的機(jī)械零構(gòu)件都是具有三維空間尺寸的物體，理應(yīng)作為三維對(duì)象處理，但是當(dāng)物體的幾何形狀和受力狀態(tài)處于某些特定的情況下，近似地簡(jiǎn)化為平面問題，不僅可以大大簡(jiǎn)化計(jì)算的工作量，而且其精度也完全能夠滿足所要求。如：直齒圓柱齒輪可在垂直與孔軸線的截平面內(nèi)作平面應(yīng)力分析就足以了解整個(gè)齒輪的受力狀態(tài)；大壩的橫斷面可作平面應(yīng)變分析來了解整個(gè)大壩受力情況等。本章是全書的重點(diǎn)，在這里不僅介紹彈性力學(xué)的基本知識(shí)。還將系統(tǒng)地講解有限元的基本概念、原理和方法。是學(xué)習(xí)以后各章節(jié)的基礎(chǔ)?！?.1外力、應(yīng)力、應(yīng)變和位移外力、應(yīng)力、應(yīng)變和位移的概念在材料力學(xué)中已經(jīng)學(xué)習(xí)過，由于這些概念在彈性力學(xué)、有限元法中具有和在材料力學(xué)中不同的規(guī)定，彈性力學(xué)中的規(guī)定和有限元法是完全相同的，所以在這里我們將按照彈性力學(xué)的習(xí)慣表達(dá)方法把他們集中的加以闡述。外力外界作用在物體上的作用力，可以分為兩大類：體積力分布在物體體積內(nèi)的力。如：重力、慣性力和磁性力等。單位體積的體力在坐標(biāo)軸上的分量X、Y、Z，稱為體力分量。符號(hào)規(guī)定為沿坐標(biāo)軸正向的為正，沿負(fù)向的為負(fù)。面力作用在物體表面的力。如輪壓、水壓等。它又可細(xì)分為集中力與分布力。面力在坐標(biāo)軸上的投影，表示為X、Y、Z。符號(hào)沿正軸為正，負(fù)軸為負(fù)。應(yīng)力彈性體受到外力作用后，內(nèi)部產(chǎn)生的抵抗變形的內(nèi)力。以彈性體中P點(diǎn)為定點(diǎn)的微單元體來考察。所謂微單元體，就是圖中PA、PB、PC的邊長(zhǎng)分別為dx、dy和dz。以下簡(jiǎn)稱這樣的微單元體為微元體。微元體每個(gè)面上的應(yīng)力都可以分解為三個(gè)應(yīng)力分量。以圖中紅面為例，分別是σx、τxy、τxz。應(yīng)力命名的規(guī)則：以應(yīng)力所在面垂直的坐標(biāo)軸為第一個(gè)下標(biāo)，應(yīng)力指向?yàn)榈诙聵?biāo)。如果下標(biāo)相同就用一個(gè)下標(biāo)表示。符號(hào)規(guī)定：正面上的應(yīng)力與坐標(biāo)軸同向?yàn)檎?，反之為?fù)。負(fù)面上的應(yīng)力與坐標(biāo)軸反向?yàn)檎?。反之為?fù)。所謂正面就是面的外法向與坐標(biāo)軸同向?yàn)檎?。反之為?fù)面。作用在兩個(gè)相互垂直的面上，并且垂直與該兩面交線的剪應(yīng)力互等。即：τxy=τyx；τyz=τzy；τxz=τzx如此以來，代表P點(diǎn)應(yīng)力狀態(tài)的應(yīng)力分量應(yīng)有6個(gè)，它們是：位移任一點(diǎn)的位置移動(dòng)用u、v、w表示它在坐標(biāo)軸上的三個(gè)投影分量。符號(hào)規(guī)定：沿坐標(biāo)軸正向?yàn)檎粗疄樨?fù)。應(yīng)變彈性體內(nèi)各點(diǎn)的位移在受力后一般是不相同的。各點(diǎn)之間距離的改變，從而使物體形狀發(fā)生變化，即所謂的變形。而物體的形狀總可以用它各部分的長(zhǎng)度和角度表示。長(zhǎng)度的改變稱為正應(yīng)變?chǔ)?，角度的改變稱為剪應(yīng)變?chǔ)谩Ｒ晕⒃w三個(gè)棱邊的線伸長(zhǎng)和角度的變化，就分別有和6個(gè)應(yīng)力分量相對(duì)應(yīng)的6個(gè)應(yīng)變分量，即：為與前面符號(hào)規(guī)定一致，這里對(duì)剪應(yīng)變的符號(hào)規(guī)定如下：正應(yīng)變伸長(zhǎng)為正，縮短為負(fù)；剪應(yīng)變使直角變小為正，變大為負(fù)?！?.2兩類平面問題前面我們講過，實(shí)際受力物體都是三維的空間物體，作用在其上的外力，通常也是一個(gè)空間力系，其應(yīng)力分量、應(yīng)變分量和位移也都是x、y、z三個(gè)變量的函數(shù)。但是當(dāng)所考察的物體具有某種特殊的形狀和特殊的受力狀態(tài)時(shí)，就可以簡(jiǎn)化為平面問題處理。彈性力學(xué)中的平面問題有兩類。平面應(yīng)力問題當(dāng)物體的長(zhǎng)度與寬度尺寸，遠(yuǎn)大于其厚度（高度）尺寸，并且僅受有沿厚度方向均勻分布的、在長(zhǎng)度和寬度平面內(nèi)的力作用時(shí)，該物體就可以簡(jiǎn)化為彈性力學(xué)中的平面應(yīng)力問題。我們分析以下其應(yīng)力特征。當(dāng)z=±t/2時(shí)，有σz=0、τzx=0、τzy=0。由于板較?。ㄏ鄬?duì)于長(zhǎng)度和寬度尺寸），外力沿板厚又是均勻分布的，根據(jù)應(yīng)力應(yīng)連續(xù)的假定（彈性力學(xué)中的基本假定），所以可以認(rèn)為，整個(gè)板的各點(diǎn)均有σz=0、τzx=0、τzy=0。如此以來，描述空間問題的6個(gè)應(yīng)力分量也就變?yōu)榱?個(gè)，即而且這些應(yīng)力分量?jī)H是x、y兩個(gè)變量的函數(shù)。平面應(yīng)變問題當(dāng)物體是一個(gè)很長(zhǎng)、很長(zhǎng)的柱形體，其橫截面沿長(zhǎng)度方向保持不變，物體承受平行于橫截面且沿長(zhǎng)度方向均勻分布的力時(shí)，該問題就可以簡(jiǎn)化為平面應(yīng)變問題處理。分析其應(yīng)力特征。假定其長(zhǎng)度方向?yàn)闊o限長(zhǎng)，那么任一橫截面都可以看作是物體的對(duì)稱面，如此則有該面上的點(diǎn)都有w=0，也就是橫截面上的所有點(diǎn)都不會(huì)發(fā)生Z方向的位移。由這一點(diǎn)可以推出也就有εz=0、τzx=0、τzy=0。和平面應(yīng)力相比較，平面應(yīng)變是εz=0，那么是否也就有σz=0呢？可能有同學(xué)想σ=Eε，當(dāng)然也就有σz=0，這是錯(cuò)誤的。平面應(yīng)變狀態(tài)下σz≠0的。雖然不等于零，但它也不是一個(gè)獨(dú)立的變量了，它由σx、σy的大小而決定。如此以來，獨(dú)立的應(yīng)力分量同平面應(yīng)力問題一樣也是3個(gè)：兩類問題的比較幾何特征平面應(yīng)力厚度＜＜長(zhǎng)度、寬度平面應(yīng)變厚度＞＞長(zhǎng)度、寬度為便于說明可講上述長(zhǎng)度看作為厚度受力特征外力都必須在其面內(nèi)且不沿厚度方向變化應(yīng)力特征平面應(yīng)力σz=0、τzx=0、τzy=0。εz≠0自由變形（無約束）平面應(yīng)變?chǔ)襷≠0但不是自變量、τzx=0、τzy=0。εz=0總以上比較可以看出，平面應(yīng)力是真正的2維（平面）應(yīng)力狀態(tài)，而平面應(yīng)變卻不是，而是3維應(yīng)力狀態(tài)，只不過σz不是獨(dú)立變量而是隨橫截面平面應(yīng)力分量而定。獨(dú)立變化的應(yīng)力分量只有3個(gè)，類似于平面應(yīng)力狀態(tài)。§2.3平衡微分方程彈性力學(xué)求解問題是從靜力學(xué)、幾何學(xué)和物理學(xué)三方面綜合考慮的。所以我們首先微元體應(yīng)該滿足的平衡條件——平衡微分方程。我們以平面問題為例推導(dǎo)，看看它應(yīng)該具有什么形式。首先對(duì)平面問題的微元體進(jìn)行受力分析圖，如左所示。物體靜力平衡的條件是：∑Fx(y)=0；∑M=0。先看∑Fx=0展開化簡(jiǎn)得同理可求得∑Fy=0滿足得條件，由∑M=0，列出方程如下：化簡(jiǎn)后得：略去微量項(xiàng)，可得：。這就是前面所將的剪應(yīng)力互等。對(duì)于平面應(yīng)變問題，微元體的前后面還有正應(yīng)力σz，不過它們是互等的。對(duì)于推導(dǎo)出來的結(jié)果，沒有任何影響。所以平面問題的平衡微分方程就是：寫成矩陣形式§2.4幾何方程考察平衡微分方程，其中具有三個(gè)未知變量σx、σy、τxy，，而只有兩個(gè)方程，方程具有無數(shù)個(gè)解。表明僅從靜力學(xué)關(guān)系無法求解該方程。我們必須從其它方面尋求幫助。彈性體在受到外力后，會(huì)發(fā)生位移和形變，從幾何上描述彈性體各點(diǎn)位移于應(yīng)變之間的關(guān)系，就是彈性力學(xué)中的又一個(gè)重要方程——幾何方程。仍然取截面的微元體ABCD，AB、CD邊長(zhǎng)為dx、dy，厚度為“1”位移u、v都是x、y的函數(shù)，即u(x,y)、v(x,y)，偏導(dǎo)數(shù)、表示位移分量u、v沿坐標(biāo)軸x的變化率，偏導(dǎo)數(shù)、表示位移分量u、v沿坐標(biāo)軸y的變化率，設(shè)A點(diǎn)的位移為u、v，那么B’點(diǎn)的位移就是：同理的D’點(diǎn)的位移分量由于α角在位移和形變很微小的情況下非常小，所以A’B’≈A’B”線段AB位移后的總伸長(zhǎng)量為A’B’-AB=A’B”-AB=uB’-uA=-u=∴，同理可得剪應(yīng)變由α、β兩個(gè)角度組成由于，所以，同理可得∴綜合以上幾何方程，并將它們寫成矩陣形式：由以上方程可以看出，當(dāng)彈性體的位移分量確定以后，由幾何方程可以完全確定應(yīng)變，反過來，已知應(yīng)變卻不能完全確定彈性體的位移。這是因?yàn)槲矬w產(chǎn)生位移的原因有兩點(diǎn)：變形產(chǎn)生的位移。因運(yùn)動(dòng)產(chǎn)生的位移。因此彈性體有位移不一定有應(yīng)變，有應(yīng)變就一定有位移?！?.5物理方程描述彈性體內(nèi)應(yīng)力與應(yīng)變關(guān)系的方程，我們稱之為物理方程，也叫材料的本構(gòu)方程。彈性力學(xué)通常研究的是各向同性材料，在三維應(yīng)力狀態(tài)下的應(yīng)力應(yīng)變關(guān)系。當(dāng)彈性體處于小變形條件下，正應(yīng)力只會(huì)引起微元體各棱邊的伸長(zhǎng)或縮短，而不會(huì)影響棱邊之間角度的變化，剪應(yīng)力只會(huì)引起角度的變化而不會(huì)引起各棱邊的伸長(zhǎng)或縮短。因此運(yùn)用力的疊加原理、單向虎克定律和材料的橫向效應(yīng)（泊松效應(yīng)），我們就可以很容易的推導(dǎo)出材料在三向應(yīng)力狀態(tài)下的虎克定律，也就是通常所說的廣義虎克定律。式中，E——材料線彈性模量G——材料剪切彈性模量μ——材料橫向收縮系數(shù)，即泊松系數(shù)。三者不是獨(dú)立的。具有以下關(guān)系：這些參數(shù)都是材料的固有屬性系數(shù)，可以通過查材料手冊(cè)獲得。如：鋼材的彈性模量E＝196～206GPa之間，通常取2.1×105MPa，μ＝0.24～0.28之間，也取為0.3進(jìn)行計(jì)算，G＝79GPa。將以上空間問題的物理方程運(yùn)用到平面問題，其形式如下：平面應(yīng)力問題的物理方程前面分析已知，平面應(yīng)力問題有σz=0、τzx=0、τzy=0。所以：從以上物理方程，也論證了我們前面說的εz≠0的結(jié)論，但由于它是由x和y方向應(yīng)力產(chǎn)生的附加無約束變形，所以通常不予以考慮。在有限元分析中更多的是運(yùn)用應(yīng)變表示的應(yīng)力關(guān)系，所以我們將上式變形一下：以上方程的矩陣表達(dá)形式為：簡(jiǎn)記為：式中：{σ}、{ε}為該問題的應(yīng)力、應(yīng)變向量。[D]為彈性矩陣。它是一個(gè)對(duì)稱矩陣，且只與材料的彈性常數(shù)有關(guān)。平面應(yīng)變問題的物理方程因?yàn)棣舲=0所以由空間物理方程的第三式得：，代入（1）、（2）式得同理變型為應(yīng)變表示應(yīng)力的形式：矩陣形式：也可簡(jiǎn)記為平面應(yīng)變問題的彈性矩陣不同于平面應(yīng)力問題的彈性矩陣，比較可以發(fā)現(xiàn)只需將平面應(yīng)力問題彈性矩陣[D]中的材料常數(shù)E換為E/(1-μ2)，μ換為μ/(1-μ)就得到了平面應(yīng)變問題的彈性矩陣。其實(shí)彈性矩陣的這種轉(zhuǎn)換方法，是彈性力學(xué)中將平面應(yīng)力結(jié)果，轉(zhuǎn)換到平面應(yīng)變問題結(jié)論的一般方法。因?yàn)樵趦煞N平面問題的描述方程中（平衡微分方程、幾何方程和物理方程），只有物理方程是不同的。§2.6邊界條件求解彈性力學(xué)問題實(shí)際就是在確定邊界條件下，求解8個(gè)基本方程（平面問題而言），以確定8個(gè)未知變量。所以從數(shù)學(xué)的角度看，就是求解偏微分方程的邊值問題。邊界條件的給出通常是各式各樣的。大體可以分為三類：第一類邊值問題給定物體的體力和面力條件，確定彈性體的應(yīng)力場(chǎng)和位移場(chǎng)。此類問題邊界以力的形式給出，所以也稱為應(yīng)力邊界條件。我們可以來考察一下應(yīng)力邊界的一般形式：是在Sσ面上給出的力的分量。平面問題如左圖所示，設(shè)陰影部分的微元體弧長(zhǎng)為ds，厚度為單元厚度“1”，其法線與X軸的夾角為θ，由陰影部分微元體的平衡條件可以推出：化簡(jiǎn)后得：此即為平面問題應(yīng)力邊界方程。第二類邊值問題給出彈性體的體力和物體表面各點(diǎn)得位移條件，確定彈性體得應(yīng)力場(chǎng)和位移場(chǎng)。由于以位移給出已知得邊界條件，所以也稱為位移邊界問題。一般得位移邊界條件為：在Su面上。第三類邊值問題給定彈性體得體力和一定邊界上得面力，其余邊界上的位移，確定其應(yīng)力場(chǎng)和位移場(chǎng)。由于邊界以力和位移兩種形式給出，所以也稱為混合邊界問題。針對(duì)不同的邊界條件，彈性力學(xué)求解的方法也有所不同?！?.7彈性力學(xué)的解題方法（解析法）應(yīng)力法由于第一類問題的邊界條件以應(yīng)力形式給出，所以以應(yīng)力作為基本的未知量的求解過程，就是人們通常所說的應(yīng)力法。由于平衡方程中有三個(gè)未知量，而只有兩個(gè)平衡微分方程，必須找出另外一個(gè)包含應(yīng)力分量的方程，才能求得方程解。考慮到彈性體變形前是一個(gè)連續(xù)體，變形后也應(yīng)是連續(xù)體的基本假設(shè)，所以要求微元體的變形一定要協(xié)調(diào)，才能使變形前、后，不會(huì)發(fā)生裂縫、重疊等現(xiàn)象。要使變形協(xié)調(diào)，就要研究幾何方程。前面介紹的平面問題幾何方程如下：分別對(duì)εx、εy、求y、x的二階偏導(dǎo)，然后相加：上式表明三個(gè)應(yīng)變分量之間應(yīng)滿足的連續(xù)性條件，我們稱之為形變協(xié)調(diào)方程（相容方程）。通過物理方程，使上述的形變協(xié)調(diào)方程換成應(yīng)力表示的形式，使之與平衡微分方程就構(gòu)成了應(yīng)力法中需要求解的方程組。具體我們來看看：利用物理方程消去相容方程中的形變分量（以平面應(yīng)力為例）利用平衡微分方程，消去上述公式中的剪應(yīng)力eq\o\ac(○,1)eq\o\ac(○,2)eq\o\ac(○,1)式對(duì)X求偏導(dǎo)，eq\o\ac(○,2)式對(duì)y求偏導(dǎo)，然后兩者相加代入相容方程，化簡(jiǎn)對(duì)于平面應(yīng)變而言，運(yùn)用前面講過的物理方程的轉(zhuǎn)換方法，只需將上式中的μ代以μ/(1-μ)就可以了。最終求解的方程組平面應(yīng)力問題平面應(yīng)變問題三個(gè)微分方程，三個(gè)未知變量，再考慮邊界條件，即可求得。如果是單連體（只具有唯一的封閉邊界）的對(duì)象，滿足了以上方程組后就是實(shí)際的解。但對(duì)于多連體（具有多個(gè)封閉邊界）的對(duì)象，中還包含有待定系數(shù)，這些待定系數(shù)會(huì)導(dǎo)致位移的解出現(xiàn)多值性。所以對(duì)于多連體的問題，還應(yīng)考慮位移的單值條件，才能最終確定。該部分的內(nèi)容可以參見徐芝綸編的《簡(jiǎn)明彈性力學(xué)教程》中圓環(huán)受均布?jí)簯?yīng)力的情況（P87）位移法位移法主要針對(duì)第二類邊界條件問題求解。解題步驟：1）改寫物理方程使之成為應(yīng)變表示應(yīng)力的形式2）應(yīng)用幾何方程表示以上得到公式中的應(yīng)變3）將它們代入平衡微分方程經(jīng)整理最后得到的位移法求解平面應(yīng)力問題方程為：兩個(gè)未知量，兩個(gè)方程，再加以邊界條件即可求得問題的解。以上介紹的解析法中，應(yīng)力法和位移法是求解彈性力學(xué)問題的基本方法。但都需要解聯(lián)立的偏微分方程組。求解過程中的數(shù)學(xué)難度，常常導(dǎo)致這種求解是無法進(jìn)行的。由于應(yīng)力法在體力為常量的情況下可以進(jìn)一步簡(jiǎn)化為求解一個(gè)單獨(dú)的微分方程的問題，所以應(yīng)力法在解析法中相對(duì)應(yīng)用較多。但即使這樣，在應(yīng)力法中，也常常采用逆解法或半逆解法?！?.8常體力情況下應(yīng)力法的簡(jiǎn)化、應(yīng)力函數(shù)及實(shí)例分析我們前面講述了彈性力學(xué)的三大方程，及應(yīng)用這三大方程的應(yīng)力法和位移法解題步驟。但是也說了要求解這些聯(lián)立的偏微分方程在數(shù)學(xué)上是存在很大難度的。很多情況下，根本無法進(jìn)行。那么彈性力學(xué)如何在實(shí)際中進(jìn)行應(yīng)用，它們和我們前面學(xué)過的材料力學(xué)區(qū)別究竟在哪里？我們將通過這一節(jié)的學(xué)習(xí)，一方面了解如何應(yīng)用這些彈性力學(xué)的方程求解問題，另一方面加深對(duì)力學(xué)概念的理解，建立力學(xué)分析問題的直觀感覺，為建立有限元模型打好基礎(chǔ)。我們知道在大多數(shù)情況下，我們分析的對(duì)象，體力是常數(shù)，它不隨x、y坐標(biāo)變化。如此以來，前面講解的第三個(gè)方程（應(yīng)力表示的相容方程），就可以簡(jiǎn)化為了：簡(jiǎn)記為：以上方程稱為拉普拉斯微分方程，數(shù)學(xué)上也稱之為調(diào)和方程，滿足調(diào)和方程的函數(shù)稱之為調(diào)和函數(shù)，及這里的。是拉普拉斯算子。這樣以來常體力情況下的應(yīng)力法方程就是：以上方程都不含有材料常數(shù)E、μ，所以平面應(yīng)力和平面應(yīng)變兩類問題具有相同的方程，這表明：在單連體問題中，只要邊界相同、受同樣的分布外力，應(yīng)力分布與材料無關(guān)；也與是平面應(yīng)力還是平面應(yīng)變的狀態(tài)無關(guān)。以上結(jié)論的意義：eq\o\ac(○,1)彈性力學(xué)平面解答的應(yīng)用范圍加寬。eq\o\ac(○,2)為實(shí)驗(yàn)應(yīng)力分析提供了理論依據(jù)（光彈實(shí)驗(yàn)）下面我們考察平衡方程：其解由其次微分方程的通解，加上任意一組特解組成。特解我們可以很容易找到。如：所以現(xiàn)在關(guān)鍵是找其次方程的通解。由第一個(gè)方程，可得：，由數(shù)學(xué)微分理論，該式是一個(gè)函數(shù)全微分的充要條件。所謂全微分就是有一個(gè)函數(shù)且同理由第二式可得：由剪應(yīng)力公式又知存在一個(gè)函數(shù)φ，可以使∴故：由于應(yīng)力與函數(shù)φ存在這樣的關(guān)系，因此函數(shù)φ即是應(yīng)力函數(shù)。我們用應(yīng)力函數(shù)來表示相容方程：上式表明φ是重調(diào)和函數(shù)。前面講過在彈性力學(xué)中，常用逆解法和半逆解法。所謂逆解法就是設(shè)定各種滿足相容方程的應(yīng)力函數(shù)，運(yùn)用σx、σy與φ的關(guān)系，求得應(yīng)力分量，再考察其滿足何種邊界條件，從而知曉這樣的應(yīng)力函數(shù)可以解決什么問題。所謂半逆解法就是根據(jù)彈性體的邊界形狀和受力關(guān)系，設(shè)定部分應(yīng)力分量為何種形式的函數(shù)，從而確定應(yīng)力函數(shù)φ，在運(yùn)用應(yīng)力函數(shù)求出所有的應(yīng)力分量，根據(jù)邊界條件確定應(yīng)力分量應(yīng)具有的最終形式。下面我們來看一個(gè)半逆解法的例子。運(yùn)用逆解法求簡(jiǎn)支梁受均布載荷的應(yīng)力分布。由材料力學(xué)知，彎曲應(yīng)力主要由彎矩M引起，剪應(yīng)力由剪力引起，而擠壓應(yīng)力由分布載荷q引起。現(xiàn)在q為不隨x變化的常量。因此我們?cè)O(shè)σy不隨x坐標(biāo)變化，即，因此我們對(duì)x積分：上式中f1(y)、f2(y)是待定函數(shù)。由于應(yīng)力函數(shù)必須滿足相容方程，所以：eq\o\ac(○,1)代入到式eq\o\ac(○,1)中考察上式可以看出它是一個(gè)x的二次方程，所以一般情況下只有兩個(gè)根。也就是說只有兩個(gè)位置能夠滿足上式。但我們對(duì)相容方程的要求是絕對(duì)滿足。也就是要求在整個(gè)梁的范圍內(nèi)都滿足。所以只有該方程的系數(shù)項(xiàng)和自由項(xiàng)全部為零。即：∴eq\o\ac(○,2)eq\o\ac(○,3)eq\o\ac(○,4)公式中的A.、B、C…K都是待定系數(shù)。公式eq\o\ac(○,3)、eq\o\ac(○,4)中分別省掉了常數(shù)項(xiàng)和一次項(xiàng)、常數(shù)項(xiàng)。這是由于f1(y)和f2(y)分別是應(yīng)力函數(shù)中x的一次項(xiàng)和常數(shù)項(xiàng)的原因，這樣處理不會(huì)對(duì)應(yīng)力分量產(chǎn)生影響。最后求出的應(yīng)力函數(shù)為：由應(yīng)力與應(yīng)力函數(shù)的關(guān)系，可以求出各個(gè)應(yīng)力分量：由于以上求得的應(yīng)力分量滿足了平衡方程和相容方程，所以只需根據(jù)邊界確定A…K的系數(shù)，就求得了該問題的解。根據(jù)對(duì)稱性，知道為偶函數(shù)，為奇函數(shù)，所以有E=F=G=0通常梁的跨度遠(yuǎn)大于梁的深度，因此上下邊界是主要邊界，它們必須滿足。將它們代入的表達(dá)式，并且考慮E=F=G=0以上四個(gè)方程解四個(gè)未知數(shù)，求得：B=0將他們代回到應(yīng)力分量的表達(dá)式中，也就有：左右兩個(gè)邊界，由于前面已經(jīng)考慮了對(duì)稱性，所以這個(gè)僅考慮優(yōu)邊界。沒有水平力。要x=l時(shí)，σx=0，考察σx的表達(dá)式，除非q=0。而這和已知條件相違背。所以在這個(gè)邊界上我們只能要求部分滿足。運(yùn)用圣維南原理運(yùn)用等效力系代替它。（這樣產(chǎn)生的誤差只在力作用點(diǎn)附近較大）。運(yùn)用的等效力系就是合成力系為平衡力系：合力等于0合力矩等于0由第一個(gè)條件得K=0，（奇函數(shù)在對(duì)稱區(qū)間上的積分為零）；由第二個(gè)條件得可以證明剪應(yīng)力的合力為-ql。即最終求得的結(jié)果，加以整理：由于厚度為“1”，此時(shí)其慣性矩，靜矩（計(jì)算見圖）任意一點(diǎn)的彎矩剪力所以上式中的應(yīng)力分量可以改寫為：各項(xiàng)應(yīng)力的分布σx第一項(xiàng)為主應(yīng)力項(xiàng)，與材料力學(xué)中的結(jié)果完全一致。第二項(xiàng)為應(yīng)力修正項(xiàng)。當(dāng)L/h＞4時(shí)，僅占主項(xiàng)的1/60；當(dāng)L/h=2時(shí)，僅占主項(xiàng)的1/15。所以對(duì)于深梁的工程構(gòu)件不容忽視第二項(xiàng)的影響?！?.8虛功方程從前一節(jié)深梁的例子，我們可以看到，彈性力學(xué)解析求解的過程是非常復(fù)雜的。這樣的求解對(duì)實(shí)際工程情況說來，很多情況根本是不可能的。所以長(zhǎng)期以來，技術(shù)人員就一直探求數(shù)值求解的方法。有限元法是其中最成功的方法。為分析單元特性和簡(jiǎn)化分析過程，我們還需了解單元的能量關(guān)系。因?yàn)樵诹W(xué)上，很多時(shí)候從能量的角度分析，可以大大簡(jiǎn)化分析的過程。應(yīng)變能的概念有材料力學(xué)知，彈性體在受到外力作用發(fā)生變形的過程中，彈性體內(nèi)部會(huì)存儲(chǔ)變形勢(shì)能——應(yīng)變能。在單向應(yīng)力場(chǎng)中，單位體積的應(yīng)變能的計(jì)算可以表示為：對(duì)于平面問題，有三個(gè)應(yīng)力分量和與之對(duì)應(yīng)的應(yīng)變分量。由于在小變形情況下，正交力系互不影響，由力的疊加原理，所以該種情況的應(yīng)變能為：其中：整個(gè)彈性體的應(yīng)變能：上式也表示應(yīng)力在應(yīng)變上所做的總功。虛功原理和虛功方程在理論力學(xué)中，我們?cè)?jīng)學(xué)到一個(gè)虛功原理，也稱虛位移原理。其基本思想就是：假定加在物體上一個(gè)可能的、任意的、微小的位移，在平衡條件下，物體的約束反力所做虛功應(yīng)等于外力所做虛功。因?yàn)槟芰勘仨毷睾恪Ｔ谶@里所說的可能的虛位移是指位移必須滿足的約束條件。任意的是指位移的方向和類型是隨意的。把這一原理運(yùn)用到現(xiàn)在的彈性體中，衡量彈性體應(yīng)滿足的平衡能量關(guān)系就是：假定加在彈性體上一個(gè)可能的、任意的、微小的位移，在平衡條件下，彈性體內(nèi)的應(yīng)變能應(yīng)等于外力所做虛功。同樣是因?yàn)槟芰勘仨毷睾恪＿\(yùn)用這一原理，我們可以推到有限元中廣泛用到的虛功方程。假定彈性體發(fā)生、的虛位移，則由幾何方程得：現(xiàn)考察彈性體微元體和邊界處微元體上的力所做虛功：eq\o\ac(○,1)內(nèi)部微元體上的力所做虛功左面的應(yīng)力虛功右面的應(yīng)力虛功左、右兩面上的虛功之和（略去高階微量，并考慮）同理得剪應(yīng)力的虛功之和體力X的虛功同樣的考慮Y方向的以及以及體力Y的虛功，然后疊加成內(nèi)部微元體上的虛功如下：eq\o\ac(○,2)邊界上的微元體設(shè)斜邊中點(diǎn)處的虛位移、，形心處的應(yīng)力為、、那么在直角邊上的應(yīng)力和位移均有一個(gè)負(fù)的增量，如圖所示，虛功計(jì)算為：（略去了高階微量）同理的dy直角邊上的剪應(yīng)力虛功為：代入已有的幾何關(guān)系：，斜面上面力的虛功體積力的虛功同樣的方法求得另一面上正應(yīng)力與剪應(yīng)力的虛功，全部相加即得斜邊微元體上的虛功之和為：支反力處的虛位移為零，所以支反力不做功，將dw1+dw2，并對(duì)整個(gè)體積積分，可以得到整個(gè)彈性體內(nèi)的總虛功：根據(jù)平衡微分方程和靜力邊界條件，上式的第一、第二項(xiàng)都是零，所以彈性體的總虛功為：根據(jù)能量守恒，它應(yīng)與外力的虛功相等由于該等式引入了平衡方程和邊界方程，所以上式虛功方程等價(jià)于靜力平衡條件（內(nèi)部和邊界微元體）。不同之處在于它是一種能量的表示形式。為了便于有限元中方便運(yùn)用，引入廣義力和物理方程，虛功方程變形為：綜合以上推到過程，虛功方程表達(dá)的物理概念就是：“若彈性體處于平衡狀態(tài)，那么外力在虛位移上作的虛功應(yīng)等于應(yīng)力在應(yīng)變上作的虛應(yīng)變功，或者說等于虛應(yīng)變能”。§2.9平面問題單元?jiǎng)澐钟邢拊ㄔ谄矫鎲栴}進(jìn)行分析時(shí)，才采用三角形單元和四邊形單元、或者矩形單元，三角形單元的優(yōu)點(diǎn)是簡(jiǎn)單且對(duì)結(jié)構(gòu)的不規(guī)則邊界逼近好，而矩形單元卻更能反映實(shí)際彈性體內(nèi)部的應(yīng)力應(yīng)變變化。這兩點(diǎn)我們會(huì)逐漸向大家說明。所以一般說來，有限元分析，單元?jiǎng)澐值拿芏群蛦卧N類選取，對(duì)計(jì)算結(jié)果起重要作用。一般單元?jiǎng)澐衷矫芗Y(jié)果越精確。單元多也導(dǎo)致求解的線性方程組階數(shù)增高，要求計(jì)算機(jī)的內(nèi)存也更大，計(jì)算的時(shí)間也越長(zhǎng)，分析的效率就越低。解決這一矛盾的方法就是在應(yīng)力集中區(qū)域單元?jiǎng)澐置芗恍瑧?yīng)力變化梯度小的位置，劃稀疏些，這樣就能兼顧精度與效率的關(guān)系。一般的原則是：根據(jù)結(jié)構(gòu)的受力和支承特點(diǎn)，按對(duì)稱和反對(duì)稱的性質(zhì)，簡(jiǎn)化分析模型，以減少計(jì)算分析的規(guī)模。合理布局單元的密集程度，以使計(jì)算結(jié)果精度高而計(jì)算量小。在同一單元內(nèi)，單元的特性數(shù)據(jù)和材質(zhì)數(shù)據(jù)應(yīng)保持一致。集中載荷的作用點(diǎn)和載荷密度突變處應(yīng)有節(jié)點(diǎn)。在欲知道應(yīng)力狀態(tài)、內(nèi)力情況和位移值的位置應(yīng)有節(jié)點(diǎn)。單元的選取欲分析的目標(biāo)密切相關(guān)。模型的單元?jiǎng)澐趾煤?，把所有的單元和?jié)點(diǎn)按一定的規(guī)律和順序進(jìn)行編號(hào)，選擇適當(dāng)?shù)淖鴺?biāo)系（直角、柱面和球面），以方便確定各節(jié)點(diǎn)的坐標(biāo)值?！?.10節(jié)點(diǎn)位移、節(jié)點(diǎn)力和節(jié)點(diǎn)載荷彈性體在承受外力作用后，其內(nèi)力的傳遞實(shí)際是通過單元之間的邊界來實(shí)現(xiàn)的。但我們把結(jié)構(gòu)離散化后，如果單元?jiǎng)澐值米銐蛐r(shí)，可以看成為其內(nèi)力的傳遞通過單元與單元之間的節(jié)點(diǎn)進(jìn)行傳遞。對(duì)于平面問題而言，每個(gè)節(jié)點(diǎn)都有位移和力兩個(gè)未知量，這兩個(gè)量又都是x、y的函數(shù)，注意平面問題的節(jié)點(diǎn)是不能傳遞力矩的，為什么？節(jié)點(diǎn)位移對(duì)三節(jié)點(diǎn)三角形單元而言，因有三個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)的位移都有x，y兩個(gè)分量，所以一共有6個(gè)自由度。單元節(jié)點(diǎn)位移向量可表示為：節(jié)點(diǎn)力所謂節(jié)點(diǎn)力，就是單元對(duì)節(jié)點(diǎn)或節(jié)點(diǎn)對(duì)單元作用的力，它是彈性體內(nèi)部的作用力，也就是我們常說的內(nèi)力。和上面相同的道理，節(jié)點(diǎn)力向量為：節(jié)點(diǎn)載荷節(jié)點(diǎn)載荷是指作用在節(jié)點(diǎn)上的外力，包括：eq\o\ac(○,1)直接作用在節(jié)點(diǎn)上的外力eq\o\ac(○,2)經(jīng)等效處理后，移植到節(jié)點(diǎn)上的等效力?？捎肵i、Yi.表示。由力平衡條件知，節(jié)點(diǎn)要保持平衡，那么作用在節(jié)點(diǎn)上的載荷應(yīng)等于節(jié)點(diǎn)內(nèi)力的合力。即：，，所有的節(jié)點(diǎn)載荷向量表示為§2.11三節(jié)點(diǎn)三角形單元的位移模式和形函數(shù)彈性結(jié)構(gòu)受外載荷后，內(nèi)部各點(diǎn)的位移變化規(guī)律，一般都是很復(fù)雜的。很難找到一個(gè)函數(shù)來描述整個(gè)結(jié)構(gòu)內(nèi)部各點(diǎn)的位移變化規(guī)律。但當(dāng)把整個(gè)結(jié)構(gòu)離散以后，在一個(gè)相當(dāng)小的單元內(nèi)部，卻可以用簡(jiǎn)單的函數(shù)來近似描述單元內(nèi)部位移的這種變化規(guī)律。而不致造成很大的誤差。這就好比一條復(fù)雜的曲線，用一個(gè)函數(shù)很難描述，但在把這條曲線分段以后，對(duì)于每一條分段曲線，卻可以用直線或拋物線來描述。三角形單元的位移模式設(shè)單元內(nèi)任意一點(diǎn)的位移分量u、v是其位置坐標(biāo)x、y的線性函數(shù)，則：a1…a6是待定系數(shù)。改寫方程組成為矩陣的形式：?jiǎn)卧娜齻€(gè)頂點(diǎn)i、j、m的坐標(biāo)已知，分別為、和，因?yàn)樗鼈円彩菃卧系狞c(diǎn)，所以應(yīng)該滿足以上假定的位移變化規(guī)律。代入上式：解以上6個(gè)方程，求得6個(gè)待定系數(shù)。同理得順序輪換，A是三角形的面積在單元節(jié)點(diǎn)的順序號(hào)i,j,m必須是按逆時(shí)針排列，否則系數(shù)行列式是負(fù)值，而三角形的面積為負(fù)值，是不合理的。求得的6個(gè)系數(shù)可以用以下矩陣表示：形函數(shù)將所求得的6個(gè)待定系數(shù)代入位移模式表達(dá)式中：令順序輪換就有上式就是假定位移模式下導(dǎo)出的單元內(nèi)任一點(diǎn)位移表達(dá)式。該式的數(shù)學(xué)意義就是單元內(nèi)任一點(diǎn)的位移可以由單元節(jié)點(diǎn)的某種形式插值得到，其中的插值基函數(shù)就是Ni、Nj、Nm。對(duì)于我們目前假定的位移模式是線性函數(shù)，所以得出的插值基函數(shù)也是類似的線性函數(shù)。由此可以看出，插值基函數(shù)具有反映單元位移變化形態(tài)的特征，所以也稱之為位移形態(tài)函數(shù)，簡(jiǎn)稱形函數(shù)。[N]就是形函數(shù)矩陣。形函數(shù)的性質(zhì)eq\o\ac(○,1)單元內(nèi)任一點(diǎn)的三個(gè)形函數(shù)之和恒等于1。eq\o\ac(○,2)在單元的三個(gè)頂點(diǎn)處，有i節(jié)點(diǎn)處j節(jié)點(diǎn)處m節(jié)點(diǎn)處以上這些，可以通過簡(jiǎn)單的數(shù)學(xué)運(yùn)算進(jìn)行證明。位移模式收斂性的分析由于位移模式的選取是人為假定的，這種假定只能近似模擬單元內(nèi)位移的變化規(guī)律，由于單元?jiǎng)偠染仃嚨耐茖?dǎo)是以假定的位移模式展開的，那么這種假定的位移模式能否使有限元數(shù)值解收斂與精確解，在很大程度上就取決于所選的位移模式，通過數(shù)學(xué)證明，可以找出位移模式滿足收斂性的幾個(gè)條件：A）完備性位移模式必須包含單元的常應(yīng)變狀態(tài)。將u、v的表達(dá)式代入幾何方程，得：因?yàn)橄禂?shù)a2…a6都是常數(shù)，所以上面的應(yīng)變分量也是常量。這也表明所選的位移模式中包含有彈性體的常應(yīng)變狀態(tài)。在上面的表達(dá)式中不含x、y的變量，說明單元的應(yīng)變是常量，這也表明這種單元是一種常應(yīng)力單元。位移模式必須包含單元的剛體位移。彈性體位移一般包含兩個(gè)部分，即變形位移和沒有彈性變形的剛體位移。那么作為模擬單元位移狀態(tài)的位移模式，也就應(yīng)該能同時(shí)反映這兩部分的位移。在eq\o\ac(○,1)中已經(jīng)證明了彈性應(yīng)變的變形，下名說明他也包含有剛體位移的特征。改寫位移表達(dá)式如下：當(dāng)時(shí)，由eq\o\ac(○,1)知，，所以上式：*我們?cè)賮砜纯匆粋€(gè)點(diǎn)作剛體位移的運(yùn)動(dòng)方程。點(diǎn)M先轉(zhuǎn)到M1，再由M1移到M2，如左圖所示：**比較上面的*和**式，可以看出，，由此得出在三角形線性位移模式中，也包含了單元的剛體位移。B）協(xié)調(diào)性位移模式必須能夠反映位移的連續(xù)性。在彈性力學(xué)求解問題時(shí)，曾經(jīng)講到過變形協(xié)調(diào)方程，也說過它是彈性體變形后仍保持連續(xù)和不發(fā)生撕裂、侵入缺陷的條件。那么位移模式的選取，也應(yīng)該保證單元之間不出現(xiàn)撕裂和侵入的缺陷。由于上面假定的位移模式是線性函數(shù)，而兩點(diǎn)就能決定一條直線。由于相鄰單元的公共邊界上的位移由與之相連的兩個(gè)節(jié)點(diǎn)插值獲得，而相鄰的單元具有兩個(gè)公共的節(jié)點(diǎn)，所以通過這兩個(gè)節(jié)點(diǎn)所得的插值值，不可能出現(xiàn)不同。也就是不可能出現(xiàn)下圖a和圖b的情況。以上三個(gè)條件是選取位移模式必須考慮的。完備形條件是收斂的必要條件；協(xié)調(diào)條件是充分條件。在有限元中，滿足完備條件的的單元是完備單元，滿足協(xié)調(diào)條件的是協(xié)調(diào)單元?！?.12（三節(jié)點(diǎn)三角形）單元?jiǎng)偠染仃嚨耐茖?dǎo)上一節(jié)我們已經(jīng)建立了三角形單元的位移插值模式，并求得了形函數(shù)的方程，這樣就完成了單元內(nèi)任一點(diǎn)的位移由單元節(jié)點(diǎn)位移表示（插值）的工作，接下來運(yùn)用我們已學(xué)過的一系列知識(shí)，我們就可以完成單元?jiǎng)偠染仃嚨耐茖?dǎo)了。推導(dǎo)過程由位移插值函數(shù)導(dǎo)出單元應(yīng)變的單元節(jié)點(diǎn)位移表達(dá)式將代入上式，可得：如按分塊矩陣記憶，那么：其中i=i,j,m矩陣[B]稱為應(yīng)變矩陣，或稱為幾何矩陣。由以上計(jì)算公式知，它與單元的節(jié)點(diǎn)坐標(biāo)有關(guān)，但不隨點(diǎn)的坐標(biāo)變化，就是說在這一單元內(nèi)所有點(diǎn)的應(yīng)變是相同的。求得單元應(yīng)力的單元節(jié)點(diǎn)位移的表達(dá)式將幾何矩陣代入單元的物理方程，就有：彈性矩陣[D]是由材料常數(shù)組成的矩陣。令[S]=[D][B]，代入平面應(yīng)力的物理方程，就有∴也可以寫成分塊矩陣的形式i=i,j,m用虛功方程導(dǎo)出單元?jiǎng)偠染仃嚕▎蝿偩仃嚕┨摴Ψ匠碳俣▎卧暮穸葹閠，上式改寫為單元的虛功方程形式，虛應(yīng)變也可以用幾何方程表示代入上式由于虛位移元素是常量，所以可以提到積分號(hào)以外，并與左邊的消去（為什么？）。于是上式變?yōu)椋毫?，虛功方程就成為了單剛方程由于[B]、[D]都是不含有x，y的常數(shù)矩陣，所以雙重積分實(shí)際就是對(duì)面積積分了。A——是三角形面積將前面求得的應(yīng)變矩陣和彈性矩陣代入，然后作矩陣乘法。就得到我們要求得的矩陣計(jì)算公式。我們?cè)谶@里采用分塊矩陣的方法記憶。r,s=i,j,m注意以上是平面應(yīng)力狀態(tài)的單剛矩陣，如果是平面應(yīng)變問題呢？單元?jiǎng)偠染仃嚨男再|(zhì)單元?jiǎng)偠染仃嘯K]e表示單元抵抗變形的能力。它與通常的彈簧剛度系數(shù)k的物理意義本質(zhì)相同，只不過[K]e是一個(gè)6×6階的矩陣。共有36個(gè)元素。這是因?yàn)槿切蔚墓?jié)點(diǎn)力向量和節(jié)點(diǎn)位移向量都為6的緣故。物理意義分塊矩陣[Kij]表示的物理含義是：節(jié)點(diǎn)j處產(chǎn)生單位位移，而節(jié)點(diǎn)i、m被約束，此時(shí)在節(jié)點(diǎn)i處產(chǎn)生的節(jié)點(diǎn)力。我們寫出它的4個(gè)分量元素：在上面的矩陣中，元素的第一個(gè)下標(biāo)表示產(chǎn)生節(jié)點(diǎn)力的節(jié)點(diǎn)，第二個(gè)下標(biāo)表示產(chǎn)生節(jié)點(diǎn)位移的節(jié)點(diǎn)。上標(biāo)“1”表示水平分量，“2”表示豎直分量。而且上標(biāo)和下標(biāo)的關(guān)系是對(duì)應(yīng)的。也就是說第一個(gè)下標(biāo)對(duì)應(yīng)第一個(gè)上標(biāo)；第二個(gè)下標(biāo)對(duì)應(yīng)第二個(gè)上標(biāo)。如此就有：就表示節(jié)點(diǎn)j處產(chǎn)生豎直單位位移，在節(jié)點(diǎn)i處產(chǎn)生的水平方向的節(jié)點(diǎn)力。單元?jiǎng)偠染仃囀菍?duì)稱矩陣→這一點(diǎn)可以通過簡(jiǎn)單的數(shù)學(xué)證明如下：?jiǎn)卧仃嚨膶?duì)稱性，從物理學(xué)角度反映出的道理就是，“功的互等”。也就是在節(jié)點(diǎn)j處產(chǎn)生某一位移引起節(jié)點(diǎn)i處的節(jié)點(diǎn)力，應(yīng)等于在節(jié)點(diǎn)i處產(chǎn)生相同位移引起節(jié)點(diǎn)j處的節(jié)點(diǎn)力。單元?jiǎng)偠染仃囍械脑刂慌c單元的材料性質(zhì)、幾何形狀、尺寸大小有關(guān)，而與單元的位置無關(guān)。單元?jiǎng)偠染仃囍胁缓衋i、aj、am，上節(jié)對(duì)位移模式收斂性分析中，曾經(jīng)說明了a1、a4分別表示單元的平移分量，而由上式知單元的平移運(yùn)動(dòng)與ai、aj、am有關(guān)，而[K]e又與ai、aj、am無關(guān)，所以說它不隨坐標(biāo)軸的平動(dòng)而變。單元?jiǎng)偠染仃囀瞧娈惥仃?，即從力學(xué)的角度理解單元的剛度方程，當(dāng)給定位移時(shí)，可以求得力；當(dāng)給定力時(shí)，卻不能求得位移。因?yàn)閇K]e不存在逆矩陣，在單元沒有給出任何約束的情況下，除有應(yīng)變的可能性外，還同時(shí)有剛體位移的可能性。所以方程無解?！?.13載荷的節(jié)點(diǎn)移置前面對(duì)于有限元模型的分析時(shí)，曾經(jīng)說過，單元之間的力傳遞是通過節(jié)點(diǎn)進(jìn)行的。所以不在節(jié)點(diǎn)上的力，必須按靜力等效原則，把它們移置到節(jié)點(diǎn)上。靜力等效原則：原載荷在任何虛位移上所做的虛功，與移置到節(jié)點(diǎn)上的節(jié)點(diǎn)載荷所做的虛功相等。這種處理方法，和我們前面講到的圣維南原理相同。它們只會(huì)影響模型局部的應(yīng)力分布，而不會(huì)影響整個(gè)結(jié)構(gòu)的應(yīng)力。下面根據(jù)力的類型，分類說明處理的方法。集中力的移置如圖所示的受力示意圖，M處的力為，移置后的節(jié)點(diǎn)載荷為，虛功相等就是：∵，所以有∴由于我們現(xiàn)在一直局限于單元的載荷移置，所以上式中的{R}應(yīng)記為{R}e。如果將它寫成分量的形式：從上面公式可以清楚的看到，載荷移置結(jié)果與單元位移模式密切相關(guān)。體力移置體積力密度為，將單元中微元體的體力看作集中力，那么面力的移置設(shè)在單元的一個(gè)邊界上作用有分布力，面力的密度為，將微面積tds上的面力也看作是集中力，則有以上面力，體力的積分運(yùn)算較為復(fù)雜，在線性模式下，可按照理論力學(xué)中靜力學(xué)平行力分解的原理，直接求得等效的節(jié)點(diǎn)載荷：eq\o\ac(○,1)均質(zhì)等厚的三角形單元，受W的體力，則三個(gè)節(jié)點(diǎn)分別受到的節(jié)點(diǎn)載荷。eq\o\ac(○,2)i、j邊受到集度為q的均布面力載荷，則i、j邊節(jié)點(diǎn)受到的等效節(jié)點(diǎn)載荷。eq\o\ac(○,3)i、j邊受到線性分布的面力，i處集度為零，j處集度為q其合力為，那么i處的節(jié)點(diǎn)載荷為，j處的載荷為。這些簡(jiǎn)單的規(guī)則，對(duì)于今后實(shí)際中的應(yīng)用，可以提高效率?！?.14整體分析該過程是將離散分離的各個(gè)單元組集成離散的結(jié)構(gòu)物，從而建立模型的總剛方程。總剛方程和總剛矩陣的組集總剛度矩陣的組集原則整個(gè)離散結(jié)構(gòu)變形后，各個(gè)單元在節(jié)點(diǎn)處仍然協(xié)調(diào)地相互連接。即環(huán)繞某個(gè)節(jié)點(diǎn)的n個(gè)單元，在節(jié)點(diǎn)i處具有相同的位移。數(shù)學(xué)公式的描述：各個(gè)節(jié)點(diǎn)應(yīng)滿足靜力平衡條件。即每個(gè)節(jié)點(diǎn)上的節(jié)點(diǎn)力合力應(yīng)等于該節(jié)點(diǎn)的節(jié)點(diǎn)載荷。數(shù)學(xué)公式描述：此處的代表環(huán)繞節(jié)點(diǎn)i的所有單元節(jié)點(diǎn)力求和。在該原則指導(dǎo)下，實(shí)例的組集過程如圖a所示的平面問題，采用圖b的方法劃分網(wǎng)格，單元分析完成后，現(xiàn)將它們組集成原問題的有限元離散網(wǎng)格，求該問題的總剛矩陣。單元eq\o\ac(○,1)節(jié)點(diǎn)號(hào)為i，j，m=1，2，3（i節(jié)點(diǎn)從最小號(hào)開始，然后逆時(shí)針排列節(jié)點(diǎn)號(hào)），單剛方程：將它們展開完全相同的道理，得其他三個(gè)單元展開的單剛方程單元eq\o\ac(○,2)節(jié)點(diǎn)號(hào)為i，j，m=1，3，4單元eq\o\ac(○,3)節(jié)點(diǎn)號(hào)為i，j，m=3，5，4單元eq\o\ac(○,4)節(jié)點(diǎn)號(hào)為i，j，m=2，3，5實(shí)際上，剛度方程可以寫成如下的一般形式：k表示單元的節(jié)點(diǎn)組成首先我們運(yùn)用規(guī)則eq\o\ac(○,2)，即各個(gè)節(jié)點(diǎn)力的合力等于節(jié)點(diǎn)載荷所以也就有：我們將上面得到的方程按這個(gè)規(guī)律相加：對(duì)這個(gè)等式運(yùn)用規(guī)則eq\o\ac(○,1)，即并整理有：將上述的這些方程寫為矩陣的形式，就得到了該問題的總剛方程為：我們對(duì)這個(gè)總剛方程進(jìn)行分析，可以得出以下的規(guī)律總剛矩陣的規(guī)律總剛方程就是運(yùn)用所有節(jié)點(diǎn)的位移向量和總剛矩陣相乘，得到結(jié)構(gòu)的節(jié)點(diǎn)載荷，節(jié)點(diǎn)載荷向量（外界對(duì)彈性體作用的載荷）很容易求得，節(jié)點(diǎn)位移向量是未知量，關(guān)鍵就是總剛矩陣如何得出。如都按照上述的步驟，成千上萬個(gè)節(jié)點(diǎn)，該如何作呢？所以有必要總結(jié)總剛度矩陣的規(guī)律。eq\o\ac(○,1)當(dāng)r=s時(shí)，即子塊矩陣[Krs]是總剛矩陣主對(duì)角線上的子塊矩陣。由環(huán)繞節(jié)點(diǎn)r（s）各單元?jiǎng)偠染仃嚨南鄳?yīng)對(duì)角線子塊相加。eq\o\ac(○,2)當(dāng)r≠s時(shí)，但r、s是相鄰單元的公共邊上的節(jié)點(diǎn)時(shí)，子塊矩陣[Krs]等于兩相鄰單元?jiǎng)偩仃嚨南鄳?yīng)子塊矩陣相加。如上例中的[K13]、[K31]等。eq\o\ac(○,3)當(dāng)r≠s時(shí)，且r、s只是一個(gè)單元的邊上節(jié)點(diǎn)時(shí)，子塊矩陣[Krs]就等于該單元?jiǎng)偩仃嚨南鄳?yīng)子塊矩陣。如上例中的[K12]、[K21]等。eq\o\ac(○,2)當(dāng)r≠s時(shí)，且r、s不屬于一個(gè)單元的邊上的節(jié)點(diǎn)時(shí)，子塊矩陣[Krs]等于零。如上例中的[K15]、[K51]等。利用上述規(guī)律不僅可以檢驗(yàn)總剛組集的正確與否，而且可以直接組集總剛矩陣（手工組集）。但這種方法不利于計(jì)算機(jī)編程。下面介紹計(jì)算機(jī)組集總剛矩陣的方法總剛矩陣組集的步驟eq\o\ac(○,1)擴(kuò)大各單元?jiǎng)偠染仃?，使之成為與總剛矩陣相同的階數(shù)。eq\o\ac(○,2)按照總體節(jié)點(diǎn)的編號(hào)，將各單元?jiǎng)偠染仃嚨母鱾€(gè)子塊移到相應(yīng)的位置上。其余位置充零。eq\o\ac(○,3)把各個(gè)改造過的矩陣直接相加，就得到總剛矩陣。我們以單元eq\o\ac(○,4)為例講解步驟eq\o\ac(○,1)、eq\o\ac(○,2)單元4的i、j、m=2、5、3，所以單元?jiǎng)偠确匠淌菃蝿偩仃嚳倓偠染仃嚨男再|(zhì)總剛度矩陣由單元?jiǎng)偠染仃嚱M集而成，所以也具有單元?jiǎng)偠染仃嚨囊恍┬再|(zhì)，如相同的物理意義，位置無關(guān)性、對(duì)稱性和奇異性等，還具有以下性質(zhì)：稀疏性由規(guī)律4知，當(dāng)rs，且r、s不屬于同一單元的兩個(gè)節(jié)點(diǎn)是[Krs]=0，表明互不相關(guān)的節(jié)點(diǎn)數(shù)愈多，零子塊矩陣就愈多。一般說來相關(guān)節(jié)點(diǎn)數(shù)不超過9，而整個(gè)分析對(duì)象常常成百上千、上萬或幾十萬。如果整體有100個(gè)節(jié)點(diǎn)，那么百分之九是非零子塊，而百分之九十多都是零子塊。所以在總剛度矩陣中非零的子塊是很稀少的。帶狀性所謂帶狀性，就是指總剛矩陣中的非零子塊，集中分布在主對(duì)角線的兩側(cè)，呈帶狀分布。半帶寬值就是計(jì)算這帶狀寬度的數(shù)值。它是以排列元素最長(zhǎng)的一行，從第一個(gè)非零元素起至主對(duì)角線元素止，所有元素的個(gè)數(shù)。其數(shù)值可以由節(jié)點(diǎn)的總體編號(hào)算出。B=（相鄰節(jié)點(diǎn)號(hào)的最大差值+1）×（單個(gè)節(jié)點(diǎn)的自由度個(gè)數(shù)）對(duì)于平面問題說來，就是：B=（相鄰節(jié)點(diǎn)號(hào)的最大差值+1）×2下圖中的兩種編號(hào)方式，可以的分別計(jì)算其半帶寬值：圖a，B=（2+1）×2=6圖b，B=（6+1）×2=14半帶寬值影響計(jì)算機(jī)存取總剛矩陣所需的內(nèi)存大小。愈小愈好。由于半帶寬值是直接受節(jié)點(diǎn)總體編號(hào)的影響。所以我們?cè)诮⒂邢拊Ｐ蜁r(shí)，應(yīng)慎重考慮，采用優(yōu)化的方法。（現(xiàn)在通常在程序中都配有這樣的優(yōu)化程序）?？倓偩仃嚨膲嚎s存取技術(shù)由于總剛矩陣具有對(duì)稱性，所以我們只需存入主對(duì)角線上半帶或下半帶的元素，就可以完成解方程的運(yùn)算。此即為總剛矩陣的半帶寬存儲(chǔ)方法。假定取主對(duì)角線上半帶元素存儲(chǔ)，具體做法就是每行元素以主對(duì)角線上的元素開始，存儲(chǔ)每行半帶寬數(shù)值的元素個(gè)數(shù)。如此各元素的行號(hào)不變，改變的只是列號(hào)。新列號(hào)和原列號(hào)的關(guān)系式如下：新列號(hào)=原列號(hào)-行號(hào)+1前面圖示的例子，如果采用圖示表示總剛矩陣的存儲(chǔ)，就是：如下圖所示的情況。可以看出原總剛矩陣存儲(chǔ)需要24×24的2維數(shù)組，改為半帶寬存儲(chǔ)后，就成為了24×6的數(shù)組了。另有一維的壓縮存儲(chǔ)方法，比這種二維的半帶寬存儲(chǔ)方法壓縮更多。在這兒我們就不做介紹了。有興趣的同學(xué)可以看相關(guān)的參考書籍。講到這兒，實(shí)際我們已經(jīng)知道，每個(gè)元素在總剛矩陣中的位置，在節(jié)點(diǎn)編號(hào)完成后，就完全確定下來了。所以計(jì)算機(jī)對(duì)總剛矩陣的計(jì)算，是一部完成的。即“邊對(duì)號(hào)移置、邊改列號(hào)、邊累加”。總體邊界條件的處理前面介紹總剛矩陣的性質(zhì)時(shí)，說明了總剛矩陣是奇異矩陣。即。就是說總剛矩陣不存在逆矩陣。要求得節(jié)點(diǎn)位移的位移解，還必須引入邊界條件。1，邊界條件的類型eq\o\ac(○,1)節(jié)點(diǎn)固定，即eq\o\ac(○,2)給定節(jié)點(diǎn)位移值。即，處理方法1置“0、1”即首先在總剛矩陣[K]中與已知位移分量相對(duì)應(yīng)的行和列元素改為0，但主對(duì)角線上的元素改為1，然后在節(jié)點(diǎn)載荷向量的列陣中，與已知對(duì)應(yīng)元素的位移用代替，其余元素減去已知位移分別乘[K]中相應(yīng)元素。例：某結(jié)構(gòu)的總剛方程為，已知、，按照上面的方法修改如下：如果展開上式，馬上就有、第二行移項(xiàng)從中可以看出，這樣處理后不僅可以直接得到、，而且它們產(chǎn)生的效果也計(jì)入到了其他方程中。如果是固定位移情況，那么載荷向量中的變化是什么呢？（好好想想）2乘大數(shù)法在總剛矩陣[K]中，把與已知位移位移相對(duì)應(yīng)的行與列主對(duì)角線上的元素乘一個(gè)很大的數(shù)1010，然后把載荷向量中的對(duì)應(yīng)元素代以給定位移乘以相應(yīng)主對(duì)角線上的元素，再同樣乘以一個(gè)很大的數(shù)。同上例。考察第一個(gè)方程：兩邊同除以，由于＞＞k1j，所以u(píng)1=c1同理可得v2=c4，（思考該方法不能用于固定位移的情況，為什么？）必須說明的是，當(dāng)總剛矩陣采用半帶寬存儲(chǔ)的時(shí)候，以上邊界條件的引入也應(yīng)當(dāng)按照半帶寬存儲(chǔ)的格式修改，否則就是錯(cuò)誤的。實(shí)際應(yīng)用中，當(dāng)是第一種情況時(shí)，采用方法eq\o\ac(○,1)；當(dāng)?shù)诙N情況時(shí)，用方法eq\o\ac(○,2)。應(yīng)力計(jì)算總剛方程在引入足夠的邊界條件后，就可以進(jìn)行求解了。所謂足夠的邊界條件就是要使系統(tǒng)是一個(gè)幾何不變的系統(tǒng)。數(shù)值解方程的方法大多采用高斯消元法。當(dāng)求解出位移向量以后，就可以通過單元的幾何方程求解應(yīng)變，通過物理方程求解應(yīng)力了。這個(gè)過程稱之為回代過程。eq\o\ac(○,1)求解應(yīng)變eq\o\ac(○,2)求解應(yīng)力由于三角形單元是常應(yīng)力與常應(yīng)變單元，所以由上述方法求得的應(yīng)力或應(yīng)變看作是形心處的應(yīng)力或應(yīng)變。而且很明顯，這樣求得的應(yīng)力或應(yīng)變是不連續(xù)的。為了推算彈性體內(nèi)某一點(diǎn)的接近實(shí)際應(yīng)力，我們通常采用以下兩種方法，來平滑應(yīng)力的突變。eq\o\ac(○,1)繞節(jié)點(diǎn)平均法——就是將繞某一節(jié)點(diǎn)的各單元形心應(yīng)力加以平均，來表示該節(jié)點(diǎn)的應(yīng)力。如左圖所示。該方法在各單元面積相差不大，單元的內(nèi)部節(jié)點(diǎn)時(shí)較為準(zhǔn)確。而在外部邊界節(jié)點(diǎn)上，則不理想。所以對(duì)于邊界節(jié)點(diǎn)多采用三點(diǎn)插值的方法求得。也就是用內(nèi)部的節(jié)點(diǎn)來推算。eq\o\ac(○,2)兩單元平均法——把相鄰兩單元的應(yīng)力相加平均，表示兩單元公共邊中點(diǎn)的應(yīng)力值。邊界邊上的應(yīng)力值也采用插值方法得到。至此我們對(duì)于有限元求解平面問題的方法全部講解完畢了。下面通過一個(gè)例子來看看，具體如何計(jì)算的。§2.15例題計(jì)算一，有限元法求解平面問題的步驟：根據(jù)分析對(duì)象和給定的條件，按照一定的比例繪出計(jì)算簡(jiǎn)圖，該圖應(yīng)有各部位的尺寸、外力和支承情況。選擇合適的坐標(biāo)系，劃分單元準(zhǔn)備數(shù)據(jù)。數(shù)據(jù)包括節(jié)點(diǎn)數(shù)、編號(hào)和坐標(biāo)，單元數(shù)、編號(hào)和節(jié)點(diǎn)組成；材料特性常數(shù)；載荷大小及位置；邊界條件。計(jì)算單元的面積、應(yīng)變矩陣、彈性矩陣和單元的剛度矩陣。組集總剛矩陣，處理載荷移置，形成載荷向量。處理邊界條件，修改總剛矩陣，得到最后的剛度方程。求解線性方程組，得到位移向量。求解應(yīng)力值。例題設(shè)有如圖所示的正方形薄板，在對(duì)角線上作用有沿厚度均勻分布的載荷，其合力為2kN，板厚為1個(gè)單元厚度，為使計(jì)算簡(jiǎn)便，假定材料的μ=0，求該板的變形。分析：該薄板對(duì)稱于對(duì)角線，受力也是如此，所以可取其1/4來計(jì)算。由于是取1/4，所以該模型的兩個(gè)直角邊，都是對(duì)稱線。X的對(duì)稱線不能有x方向的位移，y軸的對(duì)稱線不能有y方向的位移。直角的頂點(diǎn)既是x軸對(duì)稱線上的點(diǎn)，又是y軸對(duì)稱線上的點(diǎn)，所以該點(diǎn)在x、y方向的位移都應(yīng)該約束。最后得出的計(jì)算模型如下。解：1準(zhǔn)備的數(shù)據(jù)A單元節(jié)點(diǎn)數(shù)據(jù)坐標(biāo)123456X001012Y211000B單元信息數(shù)據(jù)eq\o\ac(○,1)eq\o\ac(○,2)eq\o\ac(○,3)eq\o\ac(○,4)I1223J2455M3536所有單元的板厚是“1”C材料常數(shù)彈性模量E，泊送比μ=0D載荷數(shù)據(jù)E位移邊界2計(jì)算單元?jiǎng)偠染仃噯卧?面積形成的幾何矩陣彈性矩陣單元1的剛度矩陣同理可求得、和。需要說明的是可以利用位置無關(guān)性直接得出和，由于它們是1單元平移得到，所以單剛和一樣。3組集總剛4引入邊界條件，采用劃行劃列（由第一中方法演變而來）就是將已知邊界條件的節(jié)點(diǎn)位移所在的行和列全部劃掉。得到的就是需要求解的位移向量，最后得到的總剛方程為：解以上這個(gè)方程，可以得到以下的解。按照每個(gè)單元的節(jié)點(diǎn)組成，從總節(jié)點(diǎn)向量位移中挑選出每個(gè)單元的節(jié)點(diǎn)位移向量，運(yùn)用前面單元分析中得到的公式，回代就能計(jì)算應(yīng)變、應(yīng)力了。（是哪個(gè)公式？）§2.16四節(jié)點(diǎn)矩形單元前面我們運(yùn)用三節(jié)點(diǎn)三角形單元，獲得單元的應(yīng)變和應(yīng)力是不變的。這和實(shí)際情況有著明顯的差異（實(shí)際中應(yīng)力應(yīng)變?cè)趩卧惺沁B續(xù)變化的）。產(chǎn)生這種差異的原因，就是我們的位移模式過于簡(jiǎn)單。下面我們將就如何選擇位移模式，提高有限元計(jì)算精度，進(jìn)行說明。同時(shí)也可以就如何選擇單元，有個(gè)初步的認(rèn)識(shí)。首先講解四節(jié)點(diǎn)的矩形單元。一，矩形單元位移模式如圖所示的矩形單元，四個(gè)節(jié)點(diǎn)分別為i、j、m、p。單元節(jié)點(diǎn)位移向量對(duì)應(yīng)的會(huì)有8個(gè)節(jié)點(diǎn)力分量。由于矩形單元的節(jié)點(diǎn)位移向量是8個(gè)，那么根據(jù)前面我們學(xué)過的三角形推導(dǎo)過程知道，可以建立8個(gè)方程，求解8個(gè)待定系數(shù)。所以在位移模式中的待定系數(shù)可以是8個(gè)。即：由于我們選擇的是多項(xiàng)式插值。所以對(duì)于2元函數(shù)多項(xiàng)式項(xiàng)元的選擇，可以用帕斯卡三角形確定。如左圖。從該三角形可以看出，2元多項(xiàng)式的二次項(xiàng)有三項(xiàng)，即x2、y2、xy。如果選取前面兩項(xiàng)中的任何一項(xiàng)，都會(huì)造成位移模式偏惠與你所選擇的那個(gè)方向。從而使位移出現(xiàn)不對(duì)稱的情況。只有增加xy這一項(xiàng)，才能避免這種情況出現(xiàn)。也就滿足了我們材料各向同性的要求。選取局部坐標(biāo)系，將四個(gè)節(jié)點(diǎn)的坐標(biāo)代入該公式，并求解待定系數(shù)，最后可以得出和前面三角形單元類似的關(guān)系：其中的在上面公式中，代表節(jié)點(diǎn)坐標(biāo)的是什么？位移模式的討論。1，應(yīng)變和剛體位移由幾何方程得同三角形單元位移模式一樣，常應(yīng)變和剛體位移包括在位移模式之中。同時(shí)我們可以看到，應(yīng)變中還有隨x、y的變量應(yīng)變。所以在矩形單元不再是常應(yīng)變或常應(yīng)力單元了。2，在局部坐標(biāo)系中，當(dāng)x=±a（y=±b）時(shí)，位移模式是x（或y）的線性函數(shù)（將x=±a（y=±b）代入位移表達(dá)式中即可）。這表明在邊界上位移按線性規(guī)律變化，在公共邊界上的位移是連續(xù)的，滿足收斂性的充分條件。二，矩形單元的剛度矩陣應(yīng)力向量。分別代入平面應(yīng)力和平面應(yīng)變的彈性矩陣，就可以得出兩類問題的應(yīng)力解。剛度矩陣代入[B]和[D]矩陣，相乘后逐項(xiàng)積分即得單元?jiǎng)偠染仃嚕娊滩闹械木唧w計(jì)算公式。需要說明的是，由于[B]不再是常數(shù)矩陣（還有x、y），所以積分運(yùn)算較三角形單元復(fù)雜些。具體見如下矩陣。矩形單元的載荷向節(jié)點(diǎn)移置的方法同三角形單元總剛矩陣組集的原則和方法、邊界條件的處理等也同三角形單元。其中不同的位置僅是單元節(jié)點(diǎn)位移和節(jié)點(diǎn)載荷的數(shù)目，不再是6個(gè)而是8個(gè)。三，討論矩形單元的優(yōu)缺點(diǎn)：1）其位移模式推導(dǎo)出的應(yīng)力、應(yīng)變不再是常量，分布更接近實(shí)際物體中的分布。從而使這種單元的計(jì)算精度更高。2）對(duì)于斜邊性邊界和曲線邊界的擬合性差，且不便于在不同部位采用不同大小的單元。實(shí)際應(yīng)用中常需要三角形單元在大小不同的矩形單元之間進(jìn)行過渡。見下圖?！?.17六節(jié)點(diǎn)的三角形單元從以上對(duì)矩形單元的推導(dǎo)，可以看出，增加單元節(jié)點(diǎn)的數(shù)目可以提高單元位移模式的插值次數(shù)，進(jìn)而可以提高單元計(jì)算精度。沿用這樣的思路，我們對(duì)三角形單元也可以采用相同的辦法。不過增加的節(jié)點(diǎn)是在每條邊的中點(diǎn)。這樣就得到了一個(gè)六節(jié)點(diǎn)的三角形單元。如圖所示。一，位移模式由于這樣的單元有6個(gè)節(jié)點(diǎn)，所以節(jié)點(diǎn)位移向量的個(gè)數(shù)就是2×6=12個(gè)，那么在位移模式中可以有12個(gè)待定系數(shù)。參看帕斯卡三角形，我們選用以下的位移模式：將每個(gè)單元公共邊上的方程代入，可以求得一個(gè)二次的拋物線方程，而公共邊上有三個(gè)公共節(jié)點(diǎn)，所以可以唯一確定這條拋物線。這表明該插值模式滿足位移連續(xù)的收斂性充分條件。仍然可以按照前面三節(jié)點(diǎn)三角形單元的方法，分別將6個(gè)節(jié)點(diǎn)位移代入，然后解聯(lián)立方程組，求12個(gè)待定系數(shù)。這樣的計(jì)算，由于待定系數(shù)過多，計(jì)算過程也過于繁雜。所以實(shí)際中是采用面積坐標(biāo)的形式進(jìn)行計(jì)算。二，面積坐標(biāo)三角形三個(gè)頂點(diǎn)i、j、m，p為三角形中任意一點(diǎn)，其在三角形中的位置，可以用來確定。其中：A——三角形的面積。Ai——是三角形pjm的面積。Aj——是三角形pim的面積。Am——是三角形pij的面積。當(dāng)P在圖中虛線上任一點(diǎn)是，Li是相同的（三角形面積為底×高的一半），在該三角形的三個(gè)頂點(diǎn)，分別有：i，Li=1、Lm=Lj=0j，Lj=1、Lm=Li=0m，Lm=1、Li=Lj=0而且有，將節(jié)點(diǎn)的面積坐標(biāo)和前面我們推導(dǎo)出的形函數(shù)Ni進(jìn)行比較，可以知道，形函數(shù)實(shí)際就是面積坐標(biāo)。三，六節(jié)點(diǎn)三角形單元的位移模式推求如下：Ni在i節(jié)點(diǎn)處應(yīng)為1，在j、3節(jié)點(diǎn)處應(yīng)為零，面積坐標(biāo)Li雖在i節(jié)點(diǎn)處為1，但在2、3節(jié)點(diǎn)處卻為1/2，所以還是用Li作為Ni就不行了。為此我們構(gòu)造一個(gè)函數(shù)，考察該函數(shù)是否滿足形函數(shù)的性質(zhì)。由于Li在2、3節(jié)點(diǎn)處，所以很顯然上式等于零。在i節(jié)點(diǎn)處要為1，那么β=2。（令=1，可以計(jì)算出來）所以同理可得對(duì)于N1應(yīng)滿足j、m節(jié)點(diǎn)處為零，所以構(gòu)造函數(shù)，令=1（在節(jié)點(diǎn)1處應(yīng)為1），可以求得β=4所以，同理得綜合以上推導(dǎo)，六節(jié)點(diǎn)三角形單元的形函數(shù)如下：1、2、3、i、j、m循環(huán)四，單元?jiǎng)偠染仃囃茖?dǎo)過程類似于三角形單元的剛度矩陣，不同之處在于在對(duì)形函數(shù)求偏導(dǎo)時(shí)，要運(yùn)用復(fù)合函數(shù)求導(dǎo)的方法（因?yàn)椴逯祷鵑是面積坐標(biāo)的函數(shù)，面積坐標(biāo)才是x、y的函數(shù)）如：在求得[B][D]矩陣后，運(yùn)用就可以求得單剛矩陣。但是由于[B]矩陣的復(fù)雜性，所以求解[K]e也不再是很容易的一項(xiàng)工作了。必須運(yùn)用一些相關(guān)的積分公式。具體的計(jì)算公式可以見教材的相關(guān)內(nèi)容。五，討論在單元數(shù)目相同的情況下，六節(jié)點(diǎn)的三角形單元計(jì)算精度遠(yuǎn)比三節(jié)點(diǎn)三角形單元高，也比矩形單元高。但由于一個(gè)節(jié)點(diǎn)的相關(guān)節(jié)點(diǎn)數(shù)目，在6節(jié)點(diǎn)的三角形單元中大大增加，所以總剛度矩陣的的帶寬也較三節(jié)點(diǎn)三角形單元寬的多。所需的內(nèi)存也相應(yīng)增加。從理論上說來，運(yùn)用這種增加節(jié)點(diǎn)，改善單元計(jì)算精度的方法，可以不斷的運(yùn)用下去。如4節(jié)點(diǎn)矩形單元可以變?yōu)?節(jié)點(diǎn)矩形單元等等，但是實(shí)際運(yùn)用中，我們可以看到，節(jié)點(diǎn)數(shù)目的增加，導(dǎo)致單元計(jì)算復(fù)雜性的大大增加，有時(shí)可能求解不出單剛計(jì)算公式（在現(xiàn)有的計(jì)算分析理論上）。所以在有限元程序中，再高次數(shù)的單元運(yùn)用就比較少了。6節(jié)點(diǎn)三角形單元為二次單元，而三節(jié)點(diǎn)三角形單元為一次單元?！?.18平面問題的計(jì)算實(shí)例例1一直齒齒輪，齒數(shù)為20，模數(shù)是3，齒厚20mm，壓力角為20°，試分析其齒的受力狀況。例2變厚度圓筒的壓力容器，受有內(nèi)壓，試用有限元法計(jì)算圓筒內(nèi)外壁的應(yīng)力。

第二章內(nèi)容小結(jié)一，基礎(chǔ)理論1，外力、應(yīng)力、應(yīng)變和位移的概念2，兩類平面問題的區(qū)分3，彈性力學(xué)的3大方程4，位移相容方程的理解5，虛功方程二，有限元概念1，單元?jiǎng)澐值脑瓌t2，位移模式、形函數(shù)及其性質(zhì)、保證有限元解收斂的條件3，單元?jiǎng)偠染仃嚨耐茖?dǎo)過程4，非節(jié)點(diǎn)載荷的等效移置5，總剛矩陣的組集原則、規(guī)律和方法6，總剛矩陣的半帶寬壓縮存取方法7，邊界條件的類型及處理方法8，應(yīng)力求解及平滑處理的措施三，三種單元的不同之處1，三節(jié)點(diǎn)三角形單元2，四節(jié)點(diǎn)矩形單元3，六節(jié)點(diǎn)三角形單元

第三章軸對(duì)稱問題、空間問題和等參數(shù)元問題工程中有許多零件是對(duì)稱于一條軸線的回轉(zhuǎn)體，如飛輪、螺桿和發(fā)動(dòng)機(jī)的汽缸套等。當(dāng)它們受到的載荷和約束也是對(duì)稱于他的軸線時(shí)，彈性力學(xué)就將其歸屬于軸對(duì)稱問題求解。與平面問題類似，也是將空間問題，簡(jiǎn)化為二維問題處理。本章的主要內(nèi)容軸對(duì)稱問題的基本方程軸對(duì)稱問題的有限元方法空間問題的基本方程空間問題的有限元方法等參數(shù)元的概念§3.1軸對(duì)稱問題的基本方程根據(jù)軸對(duì)稱問題的特點(diǎn)，一般分析時(shí)采用圓柱坐標(biāo)系，而不是直角坐標(biāo)系，所以首先我們應(yīng)該清楚兩者之間的轉(zhuǎn)化關(guān)系：直角坐標(biāo)與柱面坐標(biāo)的關(guān)系：自變量子午面：就是通過roz的任一平面，它是一對(duì)稱面。軸對(duì)稱的彈性體變形后，其子午面上任意一點(diǎn)P只在該平面上發(fā)生位移，而與θ無關(guān)。即P點(diǎn)的應(yīng)力、應(yīng)變和位移都只是r、z的函數(shù)，所以軸對(duì)稱問題也是二維問題。彈性體變形后，不發(fā)生歪曲，所以任一點(diǎn)只有徑向位移u和軸向位移w。如左圖所示的微元體，其應(yīng)力分量相應(yīng)的應(yīng)變分量其中的與u、w類似于平面問題中的關(guān)系，即現(xiàn)在我們要推求的就是軸對(duì)稱問題中不同的，看左圖可以推出：綜合起來，就是軸對(duì)稱問題的幾何方程：其物理方程參照直角坐標(biāo)系中的可以寫出：轉(zhuǎn)化為應(yīng)變表示應(yīng)力的矩陣形式：虛功方程的柱面坐標(biāo)表達(dá)式：（因?yàn)棣抑胁缓鹊淖兞?，所以）?.2軸對(duì)稱問題的單元分析軸對(duì)稱問題的單元是圓環(huán)形單元，其斷面常常采用三角形或矩形，單元之間的連接不是點(diǎn)，而是園線，稱作節(jié)園。雖然其單元與平面問題中的單元形式不同，但因其可以取任一子午面進(jìn)行分析，所以在子午面上的單元界面卻構(gòu)成類似于平面問題的三角形網(wǎng)格，因此可以采用平面問題的分析方法，所不同之處是：1單元是圓環(huán)體2單元之間由節(jié)園鉸接3節(jié)點(diǎn)力與節(jié)點(diǎn)載荷是施加在節(jié)園上的均布力4單元的邊界是回轉(zhuǎn)面下面推導(dǎo)單元的剛度矩陣一，單元位移模式節(jié)點(diǎn)位移向量節(jié)點(diǎn)力向量選線性插值函數(shù)作為位移模式如果用插值基函數(shù)表示，就是：其中：i、j、m循環(huán)二，單元?jiǎng)偠染仃嚧灰坪瘮?shù)到幾何方程中其它部分類似于平面問題，只有其中的軸向應(yīng)變是其特色。所以公式中的i、j、m循環(huán)要引起注意的是，fi是r、z的函數(shù)，說明其對(duì)應(yīng)的應(yīng)變分量εθ是隨r、z坐標(biāo)變化，不是如平面問題中的那樣是常量。應(yīng)力×其中應(yīng)力矩陣的分塊矩陣為i、j、m循環(huán)上式和一起代入虛功方程令由于[B]中的fi、fj、fm是r、z的函數(shù)，所以上式應(yīng)該積分運(yùn)算。但為了消去f在r=0處的奇異性，同時(shí)簡(jiǎn)化積分運(yùn)算，我們?nèi)卧涡奶幍拇婀街械膔、z，如此以來[B]中就沒有r、z的變量了。這種近似處理，在網(wǎng)格劃分較密集時(shí)，完全能夠滿足工程計(jì)算的精度要求。是單元形心處的周長(zhǎng)，相當(dāng)于平面問題中的厚度t。仍然采用分塊矩陣的方法記憶單剛：r、s=i、j、m§3.3載荷的移置作用在環(huán)形單元上的體積力、面力和集中力都應(yīng)分別移置到單元的節(jié)園上，形成節(jié)園載荷。移置后的單元載荷表示為：移置的原則仍然是虛功相等的原則。在這里我們應(yīng)注意與平面問題的不同之處。集中力的節(jié)園載荷平面問題的是體力的節(jié)園載荷面力的節(jié)園載荷一，單元自重的移置如果自重沿z軸向下作用，那么可以推導(dǎo)出節(jié)點(diǎn)載荷計(jì)算公式：具體推導(dǎo)過程可見教材。二，單元離心力的移置離心力三，面力的移置如果作用有如圖所示的均布載荷，當(dāng)是均布載荷時(shí)那么可以推導(dǎo)出：§3.4計(jì)算實(shí)例例1運(yùn)用有限元法計(jì)算螺栓螺母之間受力的分布例2輪轂熱套在軸上緊配合，工作時(shí)由于輪緣溫升，導(dǎo)致輪轂與軸之間的壓力分布變化，要求計(jì)算這種變化?！?.5空間問題一，彈性力學(xué)空間問題的三大方程1，平衡微分方程2，幾何方程3，物理方程二，空間問題的有限元法一般采用四節(jié)點(diǎn)的四面體單元，如左圖所示其節(jié)點(diǎn)位移向量其節(jié)點(diǎn)力向量單元共12個(gè)自由度，取12個(gè)待定系數(shù)位移模式為：代入節(jié)點(diǎn)坐標(biāo)，求得的插值基函數(shù)：i、j、m、p循環(huán)其中：V——四面體體積，二，單剛矩陣推導(dǎo)過程完全同平面問題的過程，只是方程的維數(shù)增加了。得到的單剛矩陣，其分塊矩陣為其中的§3.6等參數(shù)單元一，四邊形等參數(shù)單元1，平面單元的位移模式及坐標(biāo)轉(zhuǎn)換如左圖所示的任意四邊形單元，仍采用同矩形單元相同的位移插值模式，即：現(xiàn)在考察其邊界的位移情況，由于任意四邊形的邊界不再象矩形單元的邊界，可以表示為x=±a（y=±b）的形式，而是的直線方程，代入其上位移模式方程：同理得以上方程表明在任意四邊形的邊界，其位移不再是線性方程了，而是一個(gè)拋物線方程，所以相鄰單元之間的位移連續(xù)條件不能滿足。（兩點(diǎn)不能唯一決定一條拋物線），所以這種位移模式不滿足位移收斂的充分條件。從數(shù)學(xué)的觀點(diǎn)來看待矩形和任意四邊形的關(guān)系，可以發(fā)現(xiàn)任意四邊形是矩形的畸變“映像”。就好比園與橢圓之間的關(guān)系。如果能夠找到一種映射關(guān)系，確保四邊形與矩形之間是一一對(duì)應(yīng)的關(guān)系，那么就可以在局部坐標(biāo)系中采用矩形單元分析（這樣可以保證位移連續(xù)的條件），然后用坐標(biāo)轉(zhuǎn)換的映射關(guān)系，轉(zhuǎn)換到總體坐標(biāo)系中成為任意的四邊形單元，由于這種映射是一一對(duì)應(yīng)的，所以也就保證了任意四邊形單元邊界的連續(xù)性，這樣就回避了采用更高次的位移插值模式的問題。這就是等參數(shù)元的基本思想。矩形單元只是分析問題的過度單元，稱作基本單元；任意四邊形單元是實(shí)際單元，稱作實(shí)際單元。為了描述問題的清楚方便，命局部坐標(biāo)系為，則根據(jù)前面的推導(dǎo)，得：其中：設(shè)矩形單元與任意四邊形單元的轉(zhuǎn)換關(guān)系如：將四個(gè)節(jié)點(diǎn)的局部坐標(biāo)值和總體坐標(biāo)值代入，就可以確定8個(gè)待定系數(shù)，如果我們也寫成型函數(shù)的形式，就有：*其中的插值基函數(shù)同上。這種位移以節(jié)點(diǎn)位移值插值，位置以節(jié)點(diǎn)坐標(biāo)值插值，而且插值參數(shù)相同，插值基函數(shù)也相同，具有這樣性質(zhì)的單元就稱為“等參數(shù)單元”，簡(jiǎn)稱“等參數(shù)元”。以上的轉(zhuǎn)換關(guān)系，在數(shù)學(xué)上是可以證明其正確性的。我們?cè)谶@兒只對(duì)其進(jìn)行一下說明。當(dāng)時(shí)，是矩形單元的jm邊，將該值代入式*，**該式是以節(jié)點(diǎn)j、m坐標(biāo)值為參數(shù)的直線方程，顯然它也是實(shí)際單元的jm邊。以代入式**，得、，是實(shí)際單元的m節(jié)點(diǎn)。以代入式**，得、，是實(shí)際單元的j節(jié)點(diǎn)。二等參數(shù)單元的單剛矩陣是ξ、η的函數(shù)，而ξ、η與整體坐標(biāo)x、y之間又有坐標(biāo)轉(zhuǎn)換關(guān)系式，所以它是復(fù)合函數(shù)求導(dǎo)，即：代入形函數(shù)Ni令該矩陣稱為坐標(biāo)轉(zhuǎn)換的Jacobian矩陣其中以上的推導(dǎo)過程說明了，原本形函數(shù)對(duì)整體坐標(biāo)x、y求偏導(dǎo)的運(yùn)算，轉(zhuǎn)變?yōu)榍缶植孔鴺?biāo)ξ、η的偏導(dǎo)。原本對(duì)實(shí)際單元的運(yùn)算轉(zhuǎn)為對(duì)基本單元的運(yùn)算，附加求雅可比矩陣的運(yùn)算。逆矩陣中的四個(gè)元素：由于[B]矩陣都是ξ、η的函數(shù)，所以單元?jiǎng)偠染仃囈矐?yīng)轉(zhuǎn)換到局部坐標(biāo)系下。首先看左圖。沿局部坐標(biāo)的ξ、η作微向量、，由于在ξ方向只是ξ變化，而η不變，所以同理得微向量在兩坐標(biāo)軸上的投影由它們組成的四邊形面積：而，所以如果寫為分塊矩陣的形式：以上表達(dá)式是一個(gè)十分沉繁的積分運(yùn)算，一般都只能用高斯數(shù)值積分運(yùn)算。四邊形單元的邊界逼近性較矩形單元大大提高了，但仍然避免不了折線代替曲線的情況，如果要構(gòu)造曲線邊界的單元就必須增加單元邊界上的節(jié)點(diǎn)數(shù)，使之稱為由三個(gè)節(jié)點(diǎn)插值才能求得的情況。如此就可