第三章平面問題的有限元§3.1三角形單元§3.2矩形單元§3.3平面等參元§3.4平面單元應用比較分析第三章平面問題的有限元3.1三角形單元§3.1.1位移函數§3.1.2應變與應力矩陣§3.1.3形函數的性質§3.1.4單元剛度矩陣及其特點§3.1.5等效節點載荷§3.1.6整體剛度矩陣及其特點§3.1.7邊界條件§3.1.8算例§3.1.9收斂準則§3.1.10從能量原理推導剛度矩陣本節主要研究三節點三角形單元，首先建立以單元節點位移表示單元內各點位移的關系式。設單元e的節點編號為i、j、m，如圖所示。由彈性力學平面問題可知，每個節點在其單元平面內的位移有兩個分量，所以整個三角形單元將有六個節點位移分量，即六個自由度。用列陣可表示為式中的子向量為

（i、j、m輪換）式中：——節點i在x軸方向的位移分量；

——節點i在y軸方向的位移分量。3.1.1位移函數在有限單元法中，雖然是用離散化模型來代替原來的連續體，但每一個單元體仍是一個彈性體，所以在其內部依然是符合彈性力學基本假設的，彈性力學的基本方程在每個單元內部同樣適用。從彈性力學平面問題的解析解法中可知，如果彈性體內的位移分量函數已知，則應變分量和應力分量也就確定了。但是，如果只知道彈性體中某幾個點的位移分量的值，那么就不能直接求得應變分量和應力分量。因此，在進行有限元分析時，必須先假定一個位移模式。由于在彈性體內，各點的位移變化情況非常復雜，很難在整個彈性體內選取一個恰當的位移函數來表示位移的復雜變化，但是如果將整個區域分割成許多小單元，那么在每個單元的局部范圍內就可以采用比較簡單的函數來近似地表示單元的真實位移，將各單元的位移模式連接起來，便可近似地表示整個區域的真實位移函數。這種化繁為簡、聯合局部逼近整體的思想，正是有限單元法的絕妙之處。3.1.1位移函數基于上述思想，可以選擇一個單元位移模式，單元內各點的位移可按此位移模式由單元節點位移通過插值而獲得。假設三節點三角形單元內任意點的位移是x、y的線性函數，如下是待定常數，因三角形單元共有六個自由度，且位移函數在三個節點處的數值應該等于這些點處的位移分量的數值。假設節點i、j、m的坐標分別為（式中～、）、（）、（），代入位移函數，得3.1.1位移函數從上式中可以解得式中：——三角形i、j、m的面積。為保證求得的面積為正值，節點i、j、m的編排次序必須是逆時針方向，如圖所示。的計算公式，如下3.1.1位移函數將式所求得得系數代入位移函數，得式中：（i、j、m輪換）3.1.1位移函數同理可以求出單元內任意點的y方向位移，如下令

（i、j、m輪換）則單元內任意點的位移可以表示成如下形式式中：——單元內位移矢量3.1.1位移函數位移函數可以進一步寫成如下形式式中：

（i、j、m輪換）是坐標的函數，它們反映了單元的位移狀態，所以一般稱之為形狀函數，簡稱形函數。顯然，形函數決定了單元內的“位移模式”，反映了i節點位移對單元內任意點位移的貢獻。3.1.1位移函數有了單元的位移模式，就可以利用平面問題的幾何方程求出單元內任意點的應變，式中：——單元應變列陣；

——單元應變轉換矩陣，或稱單元幾何矩陣。3.1.2應變與應力矩陣單元應變轉換矩陣，為3×6矩陣，其子矩陣為（i、j、m輪換）可以看出，表示i節點位移對單元應變的貢獻率。當單元確定后，也就確定了，此時單元內的應變僅依賴于節點的位移。對于三節點三角形單元，由于面積和

等都是常量，所以單元應變轉換矩陣中的諸元素都是常量，因而單元中各點的應變分量也都是常量，即單元內各點的應變相同，通常稱這種單元為常應變單元。、3.1.2應變與應力矩陣求得應變之后，再將應變代入物理方程，即公式（2.19），便可推導出以節點位移表示的應力，如下式中：——單元應力矩陣。令則式中：——單元應力轉換矩陣。3.1.2應變與應力矩陣的子矩陣為單元應力轉換矩陣可以看出，表示i節點位移對單元應力的貢獻率。當單元確定后，也就確定了，此時單元內的應力僅依賴于節點的位移。對于三節點三角形單元，由于面積和

等都是常量，所以單元應力轉換矩陣中的諸元素都是常量，因而單元中各點的應力分量也都是常量，即單元內各點的應力相同，通常稱這種單元為常應力單元。、3.1.2應變與應力矩陣可見，對于三節點三角形單元，由于所選取的位移模式是線性的，因而導致單元內各點的應變、應力相同，但相鄰單元將具有不同的應力和應變，即在單元的公共邊界上應力和應變的值將會有突變，但位移卻是連續的。3.1.2應變與應力矩陣由前文可知，形函數決定了單元內的“位移模式”，反映了i節點位移對單元內任意點位移的貢獻，形函數具有如下性質：1.形函數在各單元節點上的值，具有“本點是1、它點為零”的性質即

（i、j、m輪換）證明如下：3.1.3形函數的性質在節點j、m上類似有在節點i上3.1.3形函數的性質2.在單元內任意點上，三個形函數之和等于1即證明：將單元內任意點的坐標代入上式，得，因此容易得到此結論。由此可見，三個形函數中只有二個是獨立的?？梢宰C明，，3.1.3形函數的性質例如：假設單元各節點的位移均相同并等于，則根據此性質可得單元內任意點的位移為，顯然，該性質反映了單元剛體位移，即當單元做剛體運動時，單元內任意點的位移均等于剛體位移。3.1.3形函數的性質3.三角形單元任意一條邊上的形函數，僅與該邊的兩端節點坐標有關。例如在ij邊上，有事實上，因ij

邊的直線方程方程為3.1.3形函數的性質代入及，得則3.1.3形函數的性質另外，由可以求得利用形函數的這一性質可以證明，相鄰單元的位移分別進行線性插值之后，在其公共邊上將是連續的。例如，對圖3.2所示的單元具有ij公共邊，可知，在ij邊上有這樣，不論按哪個單元來計算，公共邊ij上的位移均由下式表示3.1.3形函數的性質由此可見，在公共邊上的位移u、v將完全由公共邊上的兩個節點i、j的位移所確定，因而相鄰單元的位移是連續的。例3.1求如圖3.3所示等邊三角形的形函數。三角形的面積將以上數據代入式（3.10），得3.1.3形函數的性質①將各節點坐標代入以上三式，有即滿足“本點是1、它點為零”的性質。3.1.3形函數的性質即滿足“在單元內任意點上，三個形函數之和等于1”的性質。②將三個形函數相加，得，即在ij邊上任意位置的位移，將完全由節點i、j的位移所確定。取其它邊同樣可以得出其結論，即滿足“三角形單元任意一條邊上的形函數，僅與該邊的兩端節點坐標有關”的性質。③以ij邊為例，在ij邊上，則3.1.3形函數的性質單元在節點處受力，單元會發生變形，也就是說單元在節點處受到的力與單元節點位移之間有必然的聯系。單元間正是通過與節點間的相互作用力連接起來成為整體，而一個單元仍是一個彈性體，如果將每一個單元的受力－位移關系找到，則整體的受力－位移關系也容易清楚。為了推導單元的節點力和節點位移之間的關系，可應用虛位移原理對單元進行分析。單元在節點處受到的力稱為單元節點力，單元在節點力的作用下處于平衡，節點力可以表示為

式中：——節點i沿x方向的節點力分量；

——節點i沿y方向的節點力分量。3.1.4單元剛度矩陣及其特點單元節點的虛位移為相應的單元內的虛位移場為

單元內的虛應變為

設單元僅在節點上受力，則單元節點力在節點虛位移上的虛功為

3.1.4單元剛度矩陣及其特點單元內的應力在虛應變上所做的功為式中：——在體積上積分。令單元剛度矩陣為由虛功原理可知，式（3.27）和（3.28）相等，再由虛位移的任意性，可得3.1.4單元剛度矩陣及其特點上式為單元剛度矩陣的定義式，其推導過程與形函數的具體表達形式、節點個數均無關，這說明該表達式具有普遍意義，則上式是單元節點力與單元節點位移之間的關系式，即單元節點平衡方程，或稱單元剛度方程，具有普遍意義。對于三節點三角形單元，由于、中的元素都是常量，這里假設單元厚度均勻，則式中：——單元面積；

——單元厚度。3.1.4單元剛度矩陣及其特點上式進一步可以表示為將幾何矩陣式（3.16）和彈性矩陣式（2.19）代入，可得每個子塊的矩陣表達式，如下3.1.4單元剛度矩陣及其特點單元剛度矩陣具有如下特點：1.中的每一個元素都是一個剛度系數假設單元的節點位移為，由，得到節點力為

3.1.4單元剛度矩陣及其特點由此可得：表示i節點在水平方向產生單位位移時，在節點i的水平方向上需要施加的節點力；表示i節點在水平方向產生單位位移時，在節點i的垂直方向上需要施加的節點力。選擇不同的單元節點位移，可以得到單元剛度矩陣中每個元素的物理含義，如下：表示s節點在水平方向產生單位位移時，在節點r的水平方向上需要施加的節點。表示s節點在水平方向產生單位位移時，在節點r的垂直方向上需要施加的節點力；表示s節點在垂直方向產生單位位移時，在節點r的水平方向上需要施加的節點力；表示s節點在水平方向產生單位位移時，在節點r的垂直方向上需要施加的節點力；因此，中的每一個元素都是一個剛度系數，表示單位節點位移分量所引起的節點力分量，中的每一個子塊表示單元某個節點位移矢量對單元某個節點力矢量的貢獻率。3.1.4單元剛度矩陣及其特點2.是對稱矩陣由式（3.32）可知，是對稱矩陣。根據對稱性，可以減少計算量和存儲量。

3.是奇異矩陣，即根據這一性質，當已知單元節點位移時，可以從式（3.31）中求出單元節點力。反之，由于單元剛陣奇異，不存在逆矩陣，因此，當已知單元節點力時不能求出單元上的節點位移。從物體變形的實際情況來說，單元剛陣奇異是必須的，因為單元除了產生變形外，還會產生剛體位移，即僅僅依靠節點力是無法唯一地確定剛體位移的。事實上，當單元的節點力為零時，單元仍可做剛體運動。3.1.4單元剛度矩陣及其特點例如，假定單元產生了x方向的單位剛體位即并假設此時對應的單元節點力為零，則由

得可以得到，在單元剛度矩陣中1，3，5列中對應行的系數相加為零，由行列式的性質可知，。3.1.4單元剛度矩陣及其特點

4.單元的剛度不隨單元或坐標軸的平行移動而改變單元的剛度取決于單元的形狀、大小、方向和彈性系數，而與單元的位置無關，即不隨單元或坐標軸的平行移動而改變。因此，只要單元的形狀、大小、方向和彈性系數相同，無論單元出現在任何位置均有相同的單元剛度矩陣，根據對稱性，可以減少計算量。，厚度為t。例3.2求圖3.4所示等腰直角三角形單元的剛度矩陣，設解根據1（0，a）、2（0，0）、3（a，0）的坐標和式（3.8）可得3.1.4單元剛度矩陣及其特點單元面積，將以上各數據代入幾何矩陣式（3.16）為由于，則平面應力與平面應變問題的彈性矩陣相同，為3.1.4單元剛度矩陣及其特點由式（3.18）可得應力轉換矩陣，如下則單元剛度矩陣為3.1.4單元剛度矩陣及其特點從以上求得的結果可以看出，單元剛度矩陣是對稱的。1、3、5列相加等于零，2、4、6列相加也等于零，單元剛度矩陣是奇異的。另外在幾何矩陣中所涉及到的元素均為“坐標差量”，而與坐標的具體值無關，從而也說明了單元剛度矩陣不隨單元或坐標軸的平行移動而改變。1）獲取單元節點信息；2）計算單元面積及參數、等；3）計算幾何矩陣中各元素；4）根據問題類型（平面應力、平面應變）計算彈性矩陣；5）計算應力轉換矩陣；6）計算單元剛度矩陣。從單元剛度矩陣的計算公式可以看出，在求解單元剛度矩陣時，可以不必求得應力轉換矩陣，但在后期求解單元應力時要用到該矩陣，所以在此處一并求出，也可以根據實際情況具體考慮。3.1.4單元剛度矩陣及其特點實際上作用在物體上的外力可以直接作用在節點上，也可以不作用在節點上，然而在有限元法中要求作用在物體上的各種外力必須用作用在節點上的力表示。這一點體現了有限法中“離散化”這一概念，即將連續體離散成只有在節點處相連的單元體，單元與單元之間的聯系只能通過節點。如前文所述，單元內任意點的位移、應力、應變等變量最終都用單元節點位移來表示。同樣，作用在物體上的各種外力也必須用作用在節點上的力表示，這一過程稱為外力的靜力等效移置，所得到的節點力稱為等效節點力。

3.1.5等效節點載荷靜力等效原則：

對于剛體來說，所謂靜力等效原則就是單元上原有的外力系和將外力系向各節點移置所得的等效節點力，二者向同一點簡化應具有相同的主失和主矩；對于彈性體來說，所謂靜力等效原則就是指單元上的外力系和將該力系向各節點移置后的等效節點力，二者在虛位移上的虛功相等，也即外力作用在單元上所引起的變形能和移置后的等效節點力在單元上引起的變形能相等，在一定的位移模式下這種移置是唯一的。3.1.5等效節點載荷1.彈性體靜力等效原則——虛功原理虛功等效的思路是：就一個單元來說，把作用在單元上的外力移置到節點上后，應當與原來的實際外力所作虛功等效。其計算方法是：對任意允許的微小虛位移，原外力所作虛功等于移置后的等效節點力所作虛功。即式中：——節點虛位移；——單元虛位移場；——集中力作用處產生的虛位移。3.1.5等效節點載荷等式左邊表示等效節點力在節點虛位移上所作虛功；右邊第一項表示作用在單元任意點c處的集中力所作虛功、第二項表示表面力所作的虛功、第三項表示體積力所作的虛功。將代入上式，并考慮的任意性，得式中：——由集中力等效的單元等效節點力；——由表面力等效的單元等效節點力；——由體積力等效的單元等效節點力；——形函數在集中力作用處的值。3.1.5等效節點載荷單元等效節點力、和的表達式分別為

以上三式不但適用三角形單元，同樣適用其它單元，具有普遍意義。下面根據以上三式，分別進行討論。3.1.5等效節點載荷可以看出，當外力作用在節點上時，如作用在i節點上，此時，，則計算所得到的等效節點力與原外力在數值上是相同的。3.1.5等效節點載荷（1）作用在單元任意點的集中力3.1.5等效節點載荷（2）作用在單元上的表面力現舉例說明，如圖所示的單元e，在ij邊上作用有表面力，假設ij邊的長度為l，單元厚度t為常數，其上任一點P距節點i的距離為s。有，，，則例3.3如圖所示，設31邊的長度為l，單元厚度為t，其上作用沿x方向的均布荷載，求其等效節點荷載。解：在31邊上，，

，則等效節點荷載為從計算結果來看，等效后的節點力是將原來作用在31邊上的均布面力平分到1、3節點上。3.1.5等效節點載荷（3）作用在單元上的體積力設單元所受體積力為，單元厚度t為常數，如所示，則3.1.5等效節點載荷當單元體是均質、等厚、比重為時，則、，有即當單元有均勻自重時，在三節點三角形單元中每個節點沿y方向的等效節點力等于單元總重量的1/3。3.1.5等效節點載荷2.剛體靜力等效原則——力系簡化（1）作用在單元任意點的集中力若集中力作用在節點上，則可直接將集中力認為是等效節點力即可。若集中力作用于單元的某一邊界上，如圖所示，則先將分解為、，然后按線段的比例把分別移置到i，j兩點上。即3.1.5等效節點載荷（2）作用在單元上的表面力如圖所示，已知在ij邊受有均布面力，單元厚度t為常數，則移置到i、j節點上的等效節點力為如圖所示，當某一邊上有三角形分布的面力時，則等效節點力為

3.1.5等效節點載荷（3）作用在單元上的體積力如圖所示，如果三角形單元ijm的重心c上受有體積力，單元厚度t為常數，則按剛體靜力等效原理可把體積力直接移置到i，j，m三個節點上，即根據圣維南原理可知，按剛體靜力等效原則移置后的等效節點力所引起的應力誤差只是局部的不會影響到整體，即物體內部距實際外力作用處相當遠的各點處應力，與其外力的具體分布情況關系很小。在線性位移模式下，彈性體和剛體按各自的靜力等效原則移置外力的結果是一致的。3.1.5等效節點載荷

整體節點載荷列陣：由各節點所受等效節點力按節點號碼以從小到大的順序排列組成的列陣。等效節點力是由集中力、表面力和體積力共同移置構成的，其中集中力包括直接作用在彈性體上的外力和邊界約束力，如支座反力。前文對單元體進行了分折，得到了單元剛度方程，但要解決問題，還必須進一步建立整個計算模型的整體剛度方程。完成這一步的關鍵，在于怎樣將單元的剛度矩陣和節點荷載列陣，分別“組裝”成整體剛度矩陣和整體節點荷載列陣。這里通過研究任意節點的平衡來建立整體剛度矩陣，該方法不但比較直觀、易懂，而且對怎樣編寫計算機程序是很有幫助的。為了研究整體剛度矩陣的組裝過程，先引入兩個概念。

整體節點位移列陣：由各節點位移按節點號碼以從小到大的順序排列組成的列陣。3.1.6整體剛度矩陣及其特點式中：，

3.1.6整體剛度矩陣及其特點不失一般性，僅考慮計算模型中有4個單元，如圖所示。四個單元的整體節點位移列陣為對每個單元都可以寫出相應的單元剛度方程，即單元節點平衡方程。例如，對①號單元，有式中：——①號單元中第i（i=1,2,3）節點所受力。為了便于組裝整體剛度矩陣，將上式以整體節點位移表示，即

3.1.6整體剛度矩陣及其特點——①號單元的擴大剛度矩陣或稱為單元貢獻矩陣。同理，對于②單元，有式中：——②號單元中第i（i=1,3,4）節點所受力；——②號單元的擴大剛度矩陣。3.1.6整體剛度矩陣及其特點對于④單元，有式中：——④號單元中第i（i=3,4,5）節點所受力；——④號單元的擴大剛度矩陣。3.1.6整體剛度矩陣及其特點對于任意一個節點，可能承受兩種力的作用，一種是其它單元給予該節點的反作用力；另一種是作用在節點上的等效節點力。對整體而言，前者屬于內力，后者屬于外力，每個節點在兩種力的作用下處于平衡。將各單元剛度方程左邊相加，即將各節點所受力相加，由于對于整體而言，單元給予節點的反作用力屬于內力，在相加過程中相互抵消，所以各節點所受力相加的結果只有外力，即等效節點力，從而得到整體節點荷載列陣，如下3.1.6整體剛度矩陣及其特點將各單元剛度方程右邊相加，從而得到整體剛度矩陣，如下3.1.6整體剛度矩陣及其特點通過以上分析得，整體節點載荷與整體節點位移之間的關系式，即結構整體有限元方程，如下式中：——整體剛度矩陣。3.1.6整體剛度矩陣及其特點整體剛度矩陣組裝的基本步驟：1）將單元剛度矩陣中的每個子塊放到在整體剛度矩陣中的對應位置上，得到單元的擴大剛度矩陣。注意對于單元剛度矩陣是按照局部編碼排列的，即對應單元剛度矩陣中的i、j、m；對于整體剛度矩陣是按照整體編碼排列的，即按節點號碼以從小到大的順序排列。在組裝過程中，必須知道單元節點的局部編碼與該節點在整體結構中的整體編碼之間的關系，才能得到單元剛度矩陣中的每個子塊在整體剛度矩陣中的位置。將單元剛度矩陣中的每個子塊按總體編碼順序重新排列后，可以得到單元的擴大矩陣。例如在圖中，單元②的局部編碼為i、j、m，對應整體編碼為1、3、4，然后將單元②剛度矩陣中的每個子塊按總體編碼順序重新排列后，可以得到單元的擴大矩陣。注意有些書籍中將局部編碼表示為1、2、3或1，2，3等；2）將全部單元的擴大矩陣相加得到整體剛度矩陣。3.1.6整體剛度矩陣及其特點通過以上組裝過程可以得到組裝整體剛度矩陣的一般規則：1）結構中的等效節點力是相關單元結點力的疊加，整體剛度矩陣的子矩陣是相關單元的單元剛度矩陣子矩陣的集成；2）當整體剛度矩陣中的子矩陣中r=s時，該節點（節點r或s）被哪幾個單元所共有，則就是這幾個單元的剛度矩陣中的子矩陣的相加。如應該是單元①－④中對應子矩陣的集成，即3.1.6整體剛度矩陣及其特點3.1.6整體剛度矩陣及其特點3）當中時，若rs邊是組合體的內邊，則就是共用該邊的兩相鄰單元剛度矩陣中的子矩陣的相加。如13邊為單元①和②的共用邊，則4）當中r和s不同屬于任何單元時，則=0。如節點r＝1和s＝5不同屬于任何單元，此時=0。上述組裝基本步驟和規則具有普遍意義，對于不同類型、不同形式的單元，只是相應節點的子矩陣的階數（節點自由度×節點自由度）可能不同，至于組裝整體剛度矩陣的規律仍是相同的。正是應為有了這種組裝規律，使得有限元法能夠很方便地應用電子計算機進行計算。例3.4如圖所示有限元模型，彈性模量為，厚度為，為簡化計算取，求整體剛度矩陣。3.1.6整體剛度矩陣及其特點單元編號①②③④整體編碼1、2、32、4、55、3、23、5、6局部編碼i、j、mi、j、mi、j、mi、j、m以整體編碼表示的單元剛度矩陣子塊解：該模型中共有6個節點，4個單元，各單元的信息如表3.1所示。表3.1各單元信息3.1.6整體剛度矩陣及其特點同例3.2類似的分析，得根據單元剛度矩陣的性質，可知，由于幾何矩陣，所以有，整體剛度矩陣中的各子塊是對所有單元相應的子塊求和得到的（實際只是對相關單元求和），其中各子塊矩陣均為2行×2列，整體剛度矩陣用子塊矩陣可以表示為3.1.6整體剛度矩陣及其特點上式中任意一子塊矩陣均為2行×2列，在計算過程中，無需將每個單元剛度矩陣進行擴大，只需判斷整體剛度矩陣子塊的下標，然后利用組裝整體剛度矩陣的一般規則進行計算，如，由圖形可知，節點2由單元①、②和③所共有，則3.1.6整體剛度矩陣及其特點，由圖形可知，25邊為單元②和③的共用邊，則，由圖形可知，節點1、5不同屬于任何單元，則采用同樣的方法進行計算，得到整體剛度矩陣為3.1.6整體剛度矩陣及其特點3.1.6整體剛度矩陣及其特點4.是奇異矩陣，在排除剛體位移后，它是正定陣3.是稀疏矩陣，非零元素呈帶狀分布用有限元方法分析復雜工程問題時，節點的數目比較多，整體剛度矩陣的階數通常也是很高的。那么在進行計算時，如果存儲整體剛度矩陣的全部元素，將會浪費較大的資源、降低計算效率。如果根據整體剛度矩陣的特點進行編寫程序，可以大大節省資源、并提高計算效率。因此有必要了解和掌握整體剛度矩陣的特點，整體剛度矩陣具有以下幾個顯著的特點：1.是對稱矩陣2.中主對角元素總是正的3.1.6整體剛度矩陣及其特點3.1.6整體剛度矩陣及其特點1.是對稱矩陣由單元剛度矩陣的對稱性和整體剛度矩陣的組裝過程，可知整體剛度矩陣必為對稱矩陣。利用對稱性，在計算機編寫程序時，只存儲整體剛度矩陣上三角或下三角部分即可。2.中主對角元素總是正的例如，剛度矩陣中的元素表示節點2在x方向產生單位位移，而其它位移均為零時，在節點2的x方向上必須施加的力；表示節點2在y方向產生單位位移，而其它位移均為零時，在節點2的y方向上必須施加的力。很顯然在此情況下力的方向應該與位移方向一致，故應為正號。3.1.6整體剛度矩陣及其特點3.是稀疏矩陣，非零元素呈帶狀分布如果遵守一定的節點編號規則，就可使矩陣的非零元素都集中在主對角線附近呈帶狀。整體剛度矩陣中的子矩陣只有當下標s等于r或者s與r同屬于一個單元時才不為零，這就說明，在第r雙行中非零子矩陣的塊數，應該等于節點r周圍直接相鄰的節點數目加1。可見，中元素一般都不是填滿的，是稀疏矩陣，且非零元素呈帶狀分布。以下圖所示的單元網格為例，其整體剛度矩陣中的非零子塊（每個子塊為2行2列）的分布情況如下圖所示，圖中陰影部分表示該子塊不為零，其它子塊部位均為零。3.1.6整體剛度矩陣及其特點3.1.6整體剛度矩陣及其特點顯然，帶狀剛度矩陣的帶寬取決于單元網格中相鄰節點號碼的最大差值D。把半個斜帶形區域中各行所具有的非零元素的最大個數叫做剛度矩陣的半帶寬（包括主對角元），用B表示，如下

B=2(D+1)通常的有限元程序，一般都利用剛度矩陣的對稱和稀疏帶狀的特點，在計算求解中，只存儲上半帶的元素，即所謂的半帶存儲。因此，在劃分完有限元網格進行節點編號時，要采用合理的編碼方式，使同一單元中相鄰兩節點的號碼差盡可能小，以便節省存儲空間、提高計算效率。3.1.6整體剛度矩陣及其特點對于同樣的有限元單元網格，按照圖（a）的結點編碼，最大的半帶寬為B=2×(10－4+1)=14；按照圖（b）的結點編碼，最大的半帶寬為B=2×(10－2+1)＝18；按照圖（c）的結點編碼，最大的半帶寬為B=2×(10－6+1)＝10。(a)(b)(c)3.1.6整體剛度矩陣及其特點4.是奇異矩陣，在排除剛體位移后，它是正定陣無約束的彈性體（或結構物）的整體剛度矩陣是奇異的，不存在逆矩陣，即關于位移的解不唯一。這是因為彈性體在外力的作用下處于平衡，外力的分量應該滿足三個靜力平衡方程。這反映在整體剛度矩陣中就意味著存在三個線性相關的行或列，所以是個奇異陣，不存在逆矩陣。例如：設彈性體在外力的作用下處于平衡，這時相應的解為，然后在給予彈性體以剛體位移而相應的節點位移，這時，仍是問題的解，因為剛體位移不會破壞平衡。注：當排除剛體位移后，整體剛度矩陣是正定矩陣。前文已經提到在排除剛體位移后，整體剛度矩陣是正定的，方程才可求得唯一解。排除剛體位移可以通過引入邊界約束條件來實現，這里介紹兩種比較簡單的引入已知位移的方法。1．代入法2．乘大數法3.1.7邊界條件3.1.7邊界條件1．代入法該方法保持方程組仍為2n×2n系統，僅對整體剛度矩陣和整體載荷列陣進行修正。下面以一個只有四個方程的簡單例子加以說明，方程如下假定系統中節點位移、，則當引入這些節點的已知位移之后，方程就變成3.1.7邊界條件若，則然后，用這組維數不變的方程來求解所有的節點位移。顯然，其解答仍為原方程的解答。在手算時，可直接將零位移約束所對應的整體剛度矩陣中的行和列直接劃去，使得整體剛陣的維數變小，更便于手算。3.1.7邊界條件2．乘大數法將中與指定的節點位移相對應的主對角元素乘上一個大數，同時將中的對應元素換成指定的節點位移值、該大數與節點位移相對應的主對角元素三者的乘積。若把此方法用于上面的例子，則方程就變成該方程組的第一個方程為解得，這種方法就是使中相應行的修正項遠大于非修正項。3.1.7邊界條件1．代入法2．乘大數法在以上的兩種方法中，代入法接近人工解法，雖然該方法比較直觀，但該方法對剛度矩陣改變較多，程序效率不高。乘大數法對剛度矩陣改變較少，工作量較小，但相乘的“大數”若取得過大，求解時會發生“溢出”、若取得太小則會引起較大的誤差。根據前面的討論，現以三角形常應變單元為例來說明應用有限元法求解彈性力學平面問題的具體步驟。1）力學模型的確定。根據工程實際情況確定問題的力學模型，并按一定比例繪制結構圖、注明尺寸、載荷和約束情況等；2）結構進行離散化。將彈性體劃分為許多單元，并對節點進行編號，確定全部節點的坐標值；對單元進行編號，并列出各單元三個節點的節點號；3）計算等效節點載荷，形成整體載荷列陣；4）計算單元剛度矩陣；5）組裝整體剛度矩陣；6）處理約束，引入邊界條件；7）求解線性方程組，得到節點位移和節點力；8）計算單元應力；9）整理計算結果（后處理部分）。3.1.8算例例3.5如圖所示為一厚度t=1cm的均質正方形薄板，邊長為2m，上下受均勻拉力=106N/m2，材料彈性模量為E，泊松比，不記自重，試用有限元法求其應力分量。

3.1.8算例解：1）力學模型的確定。由于此結構長、寬遠大于厚度，而載荷作用于板平面內，且沿板厚均勻分布，故可按平面應力問題處理，考慮到結構和載荷的對稱性，可取結構的1/4來研究。3.1.8算例表3.2節點坐標值節點坐標1234x0110y00112）結構進行離散化。為了計算簡便并能說明整個計算過程，將該1/4結構離散為兩個三角形單元。節點編號、單元劃分如圖所示，各節點的坐標值如表3.2所示，各單元編碼如表3.3所示。3.1.8算例表3.3單元編碼單元號①②整體編碼1、2、33、4、1局部編碼i、j、mi、j、m3）計算等效節點載荷，形成整體載荷列陣。按照彈性體靜力等效原則，應用公式（3.41），并參考例3.3可得單元②由均布載荷等效的節點力為3.1.8算例若按照剛體靜力等效原則，參考公式（3.45），同樣可得到以上結果?？紤]到支座反力，整體載荷列陣為4）計算單元剛度矩陣。對于單元①有，局部編碼i、j、m對應整體編碼1、2、3，利用公式（3.34）并參考例3.2，可得單元①的剛度矩陣為3.1.8算例對于單元②有，局部編碼i、j、m對應整體編碼3、4、1，由單元剛度矩陣的性質可知，單元②的剛度矩陣與單元①的剛度矩陣相同。

5）組裝整體剛度矩陣。在組裝整體剛度矩陣時，要注意單元局部編碼與整體編碼之間的關系，按照組裝整體剛度矩陣的一般規則和基本步驟，可得整體剛度矩陣為3.1.8算例

6）處理約束，引入邊界條件。由以上分析可知，整體有限元方程為3.1.8算例根據邊界約束條件：，采用代入法引入邊界條件，劃去整體剛度矩陣中1、2、4、7的行和列，得3.1.8算例

7）求解線性方程組，得到節點位移。求解上面方程組可得出相應得節點位移為3.1.8算例所以，整體節點為列陣為代入整體有限元方程，便可得整體載荷列陣為若要求得每個單元所受的節點力，則可以從整體節點位移列陣中提取單元所對應的節點位移，然后在利用，便可求得單元節點所受的力。如單元①所對應的節點位移為3.1.8算例求得單元①所對應的節點力為單元②所對應的節點位移為求得單元②所對應的節點力為注意，這里求得的單元節點位移、和單元所對應的節點力、均對應局部編碼。很顯然，每個單元均滿足力的平衡。8）計算單元應力。先求出各單元的應力矩陣、，然后再求得各單元的應力分量，如下3.1.8算例9）整理計算結果（后處理部分）。對于三節點三角形單元，單元內各點的應力值相等，算出的應力一般作為單元形心處的應力。由于單元應力為常數，整個結構的應力場呈階梯狀，在單元之間不連續。而工程上往往更加關心邊界和節點上的受力情況，因此，必須對所得到的應力再次進行處理，得到更加合理的應力場，并得到所需點上的應力值。這里介紹兩種簡單的方法，一種方法稱為節點平均法，即把環繞某一節點的各單元的應力加以平均作為該節點的應力值。例如圖中節點3的應力為3.1.8算例為了使通過這樣平均得來的應力比較接近實際情況，要求環繞節點的各單元尺寸不應相差太大。這種做法，對內節點比較好，對邊界點則可能很差。因此，邊界節點處的應力不宜直接由單元應力平均來獲得，而要根據內節點的應力構造插值函數推算出來。例如圖中邊界點1的應力，可以先用節點平均法求得節點2、3、4處的應力，在構造相應的插值函數推算邊界點1的應力，如常用的拋物線插值公式如下3.1.8算例另一種方法稱為單元平均法，即把兩相鄰單元的應力加以平均，用以表示公共邊界中點的應力。為了由這樣平均所得到的應力具有較好的精度，兩相鄰單元的面積不應相差太大。如圖中單元②和③邊界的中點處的應力為3.1.8算例在不同的有限元軟件中均具有各自的后處理方法，但無論后期怎樣處理，應力場來源于應力的計算，應力的精度主要依賴于單元的尺寸、單元的類型（位移模式）。對于有限元這種數值計算方法，一般總是希望隨著網格的逐步細分所得到的解能夠收斂于問題的精確解。根據前面的分析，可知在有限元分析中，一旦確定了單元的形狀之后，位移模式的選擇將是非常關鍵的。由于載荷的移置、應力矩陣和剛度矩陣的建立等等，都依賴于單元的位移模式，所以，如果所選擇的位移模式與真實的位移分布有很大的差別，那么就很難獲得良好的數值解。為了保證解答的收斂性，要求位移模式必須滿足以下三個條件，即1）位移模式必須包含單元的剛體位移。2）位移模式必須包含單元的常應變。3）位移模式在單元內要連續、且在相鄰單元之間的位移必須協調。3.1.9收斂準則3.1.9收斂準則

1）位移模式必須包含單元的剛體位移。也就是說，當節點位移是由某個剛體位移所引起時，彈性體內將不會產生應變。所以，位移模式不但要具有描述單元本身形變的能力，而且還要具有描述由于其它單元形變而通過節點位移引起單元剛體位移的能力。2）位移模式必須包含單元的常應變。每個單元的應變一般都是包含著兩個部分：一部分是與該單元中各點的坐標位置有關的應變（即所謂各點的變應變）；另一部分是與位置坐標無關的應變（即所謂的常應變）。從物理意義上看，當單元尺寸無限縮小時，每個單元中的應變應該趨于常量。因此，在位移模式中必須包含有這些常應變，否則就不可能使數值解收斂于正確解。3.1.9收斂準則3）位移模式在單元內要連續、且在相鄰單元之間的位移必須協調。當選擇多項式來構成位移模式時，單元內的連續性要求總是得到滿足的，單元間的位移協調性，就是要求單元之間既不會出現開裂也不會出現重疊的現象。通常，當單元交界面上的位移取決于該交界面上節點的位移時，就可以保證位移的協調性。在有限單元法中，把能夠滿足條件1和2的單元，稱為完備單元，完備單元是收斂的必要條件。滿足條件3的單元，叫做協調單元或保續單元。前面討論過的三節點三角形單元，均能同時滿足上述三個條件，因此屬于完備的協調單元，完備的協調單元是收斂的充分條件。在某些梁、板及殼體分析中，要使單元滿足條件3比較困難，所以實踐中有時也出現一些只滿足條件1和2的單元，其收斂性往往也能夠令人滿意。特別是放松條件3的單元，即完備而不協調的單元，這時要保證其收斂必須通過分片試驗的測試，關于分片試驗的內容請參見相關書籍。目前，完備而不協調的單元，已獲得了很多成功應用。3.1.9收斂準則在有限元分析中，將實際連續體分成許多單元體的組合后，根據線性或非線性位移的假定，人為地選擇一個位移場，通過這些措施所得到的模型比實際連續體的剛性要高。因而，近似模型的剛度是實際連續體剛度的上界。若選擇不協調單元，那么這種模型可能由于單元分離、疊加或單元之間形成鉸而降低剛性，這種影響就有可能使得不協調元比應用協調元所得的結果要好。不過，應用不協調單元事先不能肯定所得的剛度是真實剛度的上界。換句話說，不協調元不一定象協調單元那樣剛硬，可能比較柔軟，因此有可能會比協調單元收斂得快。3.1.9收斂準則經驗證明，根據巴斯卡（Pascal）三角形來選擇二維多項式的各項。在二維多項式中，如果包含有對稱軸一邊的某一項，那么就必須同時包含有另一邊的對稱項。選擇多項式位移模式時，還應該要考慮的一個因素是，多項式中的項數必須等于或稍大于單元邊界上的外節點的自由度數。通常是取項數與單元的外節點的自由度數相等，取過多的項數是不恰當的。從節點平衡得到結構整體有限元方程式的方法物理概念清晰，但這種方法在數學上不夠嚴謹。本節從勢能原理導出有限元方程，該方法是一種更加嚴謹、適用性更加廣泛的方法。由第2章可知，系統的總勢能為如果有集中力作用等式右邊還應加一項集中力所作的功，如下3.1.10從能量原理推導剛度矩陣當把連續體離散為有限個單元，并在各單元內假設了位移模式，在單元內應力與應變用節點位移（或單元中）來描述，、，在連續體上的整體積分等于在每個單元上分別積分再求和，如下由勢能原理可知，對于處于平衡狀態的彈性體，位移真實解應總使總勢能取極小值，反之依然，則由，得3.1.10從能量原理推導剛度矩陣上式進一步可寫為如果有集中力作用，有上式進一步可寫為上式同式（3.55），即為整體結構有限元方程。

3.1.10從能量原理推導剛度矩陣3.2矩形單元雖然三角形單元具有很好的“適應性”，幾乎任何復雜邊界的彈性體總可以劃分為三角形，并且三角形單元計算公式簡單，但精度較低。三角形單元間雖然能夠保證位移連續，但應力的精度較差，不能很好的反映彈性體內應力的準確分布規律。為了提高計算精度，準確反映彈性體內的應力狀態，可以采用一些較精密的單元類型。本節將介紹常用的矩形單元，它采用了比常應變三角形單元次數更高的位移模式，因而可以更好地反映彈性體中的位移狀態和應力狀態。另外，對一些邊界比較規則且呈直線的平面結構的分析，采用矩形單元較合適。這時單元總數可以減少，相應的原始數據準備工作和單元特征計算工作均可節省。3.2.1位移函數如圖所示的矩形單元，不失一般性，令矩形單元的長、寬分別為2a、2b。矩形單元有4個節點，共8個自由度，即共有8個節點位移，采用類似三角形單元的分析方法，同樣可以完成對矩形單元的力學特性分析。這里引入一個局部坐標系、，這樣可以推出比較簡潔的結果。如圖所示，取矩形單元的形心o為局部坐標系的原點，和軸分別與整體坐標軸x和y平行，兩坐標系存在有以下的坐標變換關系式中：、——矩形形心處坐標。矩形形心處坐標以及矩形長、寬可由下式計算3.2.1位移函數在局部坐標系中，節點i的坐標是，其值分別為±1。如節點1在局部坐標系下的坐標為（－1，－1）。由于矩形有4個節點，共8個自由度，可以選擇有8個待定參數的位移模式，如下采用類似三角形的分析方法，將各節點的局部坐標代入上式，便可以求得相應的待定參數，然后將帶回式位移模式，可得到用節點位移表示的位移模式，如下3.2.1位移函數式中：——矩形單元的形函數，i＝1，2，3，4；——形函數矩陣；——單元節點位移列陣，，i＝1，2，3，4。3.2.1位移函數（i＝1，2，3，4）形函數的表達式為3.2.1位移函數引入符號，，i＝1，2，3，4，則上式可以統一寫為可以看出，矩形單元的形函數具有和三角形單元形函數同樣的性質，即：形函數在各單元節點上的值，具有“本點是1、它點為零”的性質；在單元內任意點上，四個形函數之和等于1；單元任意一條邊上的形函數，僅與該邊的兩端節點坐標有關。3.2.1位移函數有了單元的位移模式，就可以利用平面問題的幾何方程求出單元內任意點的應變，將位移代入幾何方程，得式中的應變轉換矩陣的子塊（i＝1，2，3，4）為3.2.2應變與應力矩陣求得應變之后，再將應變代入物理方程，便可推導出以節點位移表示的應力，如下式中，應力矩陣為其子塊（i＝1，2，3，4）為3.2.2應變與應力矩陣由前面的討論可以發現，四邊形單元的位移模式，式（3.77）比常應變三角形單元所采用的線性位移模式增添了項（即相當于xy項），把這種位移模式稱為雙線性模式。在這種模式下，單元內的應變分量將不再是常量，這一點可以從的表達式中看出。另外四邊形單元的位移模式中的與三角形單元相同，它反映了剛體位移和常應變，而且在單元的邊界上(=±1或

=±1)，位移是按線性變化的，顯然在兩個相鄰單元的公共邊界上，其位移是連續的。3.2.2應變與應力矩陣由單元的應力矩陣表達式還可以看出，矩形單元中的應力分量也都不是常量。正應力、和剪應力均沿、兩個方向線性變化，即沿x、y兩個方向線性變化。正因為如此，若在彈性體中采用相同數目的節點時，矩形單元的精度要比常應變三角形單元的精度高。但是，矩形單元也有一些明顯的缺點，矩形單元不能適應斜交的邊界和曲線邊界，不便于對不同部位采用不同大小的單元，以便提高有限元分析計算的效率和精度。3.2.2應變與應力矩陣矩形單元剛度矩陣的推導過程與三節點三角形單元類似，同樣可以表示為式（3.30）的形式，即，由前文可知的推導過程與形函數的具體表達形式、節點個數均無關，該表達式具有普遍意義。若單元厚度t

為常量，則可以進一步表示為將單元剛度矩陣寫成子塊的形式，如下3.2.3單元剛度矩陣上式中每一個子塊矩陣均為2行×2列，單元剛度矩陣中的子塊矩陣的表達式為

（r、s＝1，2，3，4）將應變轉換矩陣子塊和彈性矩陣，代入上式，得3.2.3單元剛度矩陣式中：（r、s＝1，2，3，4）3.2.3單元剛度矩陣例3.6如圖3.22所示，該模型中有兩個四邊形單元，彈性模量為＝210GPa，厚度為＝0.025m，泊松比＝0.3，＝1kN，求單元所受應力。3.2.4算例解：單元①所對應的節點為2、4、3、1，單元②所對應的節點為4、6、5、3。在有限元分析過程中，首先求解單元①和②的剛度矩陣；然后組裝整體剛度矩陣，組裝的過程同三角形單元，此處不再給出；最后引入邊界條件（，＝0）并結合受力情況，求得整體節點位移列陣＝10-5×{0，0，0，0，0.1162，-0.1674，-0.1149，-0.1628，0.1514，-0.4707，-0.1568，-0.4978}T。為了求解單元應力，需要先求得每個單元所對應的節點位移列陣，結合單元的節點編號，可從整體位移列陣中提取單元①和單元②的位移列陣，如下＝10-5×{0，0，-0.1149，-0.1628，0.1162，-0.1674，0，0}T＝10-5×{-0.1149，-0.1628，-0.1568，-0.4978，0.1514，-0.4707，0.1162，-0.1674}T3.2.4算例注意：整體節點位移列陣是按照節點編號由小到大排列的，而單元位移列陣是按照單元節點編號排列的，如單元①所對應的節點為2、4、3、1，則單元①的位移列陣中的前兩個數則表示節點2的x和y方向位移。將單元①和單元②的位移列陣代入應力矩陣，可求得單元①和單元②的應力，如下3.2.4算例通過以上兩式便可以求得單元①和單元②內任意點的應力，以單元①為例，若取，，則表示單元①的13邊中點處的應力；若取，，則表示24邊中點處的應力。若計算單元形心處的應力，則取，為通過分析結果可知，單元內任意點的應力是坐標的函數，若在彈性體中采用相同數目的節點時，矩形單元的精度顯然要高于常應變三角形單元的精度。3.2.4算例3.3平面等參元前一節討論的矩形單元具有精度較高、形狀規則、便于計算自動化等優點，對于結構復雜的曲邊外形，只能通過減小單元尺寸，增加單元數量進行逐漸逼近。這樣，自由度的數目隨之增加，計算時間長，工作量大。另外，這些單元的位移模式是線性模式，是實際位移模式的最低級逼近形式，問題的求解精度受到限制。為了克服以上缺點，人們試圖找出這樣一種單元：一方面，單元能很好地適應曲線邊界和曲面邊界，準確地模擬結構形狀；另一方面，這種單元要具有較高次的位移模式，能更好地反映結構的復雜應力分布情況，即使單元網格劃分比較稀疏，也可以得到比較好的計算精度。等參數單元（等參元）就具備了以上兩條優點，等參元為任意四邊形單元，其網格劃分不受邊界形狀的限制，單元大小可以不相等，是一種精度高而且應用廣泛的單元。3.3平面等參元然而直接對任意四邊形進行單元分析是困難的，這是由于它的幾何形狀不規則，沒有統一的形狀，對各個單元逐個按不同的公式計算，因其工作量大而難以進行。為了解決這一問題，人們采用坐標變換的方法，如圖所示，把坐標系ξη（局部坐標系）內形狀規則單元（母單元）變換為另一坐標系xy（整體坐標系）中的任意四邊形單元（子單元）。然后導出關于局部坐標系形狀規整的單元（母單元）的高階位移模式的形函數，然后利用形函數進行坐標變換，得到關于整體坐標系的復雜形狀的單元（子單元）。a）母單元b）子單元3.3平面等參元從圖形變換的角度看，坐標系ξη和坐標系xy可以分別看成是母單元和子單元這兩個不同單元的坐標系，它們都是直角坐標系。而從另一角度看，坐標系ξη和坐標系xy又可以看成是同一單元（子單元）的兩種不同的坐標系。坐標系xy是子單元的直角坐標系，而坐標系ξη可看成是子單元的曲線坐標系。可以看出坐標系xy始終扮演同一角色，即子單元的直角坐標；而坐標系ξη則扮演兩種角色，它既是母單元的直角坐標，又是子單元的曲線坐標。在有限元分析中，兩者的作用是不同的。直角坐標系xy在整個結構的所有子單元中共同采用，所以稱為整體坐標系。而曲線坐標系ξη則只適用于單個獨立的子單元，所以稱為局部坐標。整體坐標在整體分析中采用，局部坐標則在單元分析中采用。平面問題中常用的等參元有四節點四邊形單元、八節點曲邊四邊形單元和6~8可變節點曲邊四邊形單元等，本節以八節點曲邊四邊形等參元為例介紹平面等參元的分析過程。八節點曲邊四邊形等參元的母單元是二維二次單元。八個節點分別為正方形的四個角點和四個邊中點，母單元（邊長為2的正方形）采用直角坐標系ξη。由于母單元為16個自由度，因此單元的位移模式為3.3.1等參元剛度矩陣式中：式中的16個常數用8個節點的位移（，）表示后，則上式為3.3.1等參元剛度矩陣采用坐標變換可使母單元的八個節點坐標與等參元的八個節點坐標建立一一對應的關系，整體坐標和局部坐標的變換式為式中：——母單元的形函數?？芍竼卧先我庵边?，例如2－3邊，有ξ＝1，則有關的形函數N2、N3和N6均是η的二次函數，其余的形函數均為零。這樣變換就成為二次非線性變換，它可將母單元的直邊2－3映射成子單元的曲邊2－3。根據等參元的思想，利用上述的形函數，可得等參元的位移模式為3.3.1等參元剛度矩陣有了單元的位移模式，就可以利用平面問題的幾何方程求出單元內任意點的應變，將上式代入式（2.12），得等參元的應變矩陣為式中：——i節點位移，（i=1，2，…，8）；3.3.1等參元剛度矩陣——單元幾何矩陣，其中為（i=1，2，…，8）在上式中，是ξ、η的函數，因此必須用坐標變換式來轉換導數關系，根據復合函數的求導法則，有（i=1，2，…，8）3.3.1等參元剛度矩陣式中的J稱為雅可比（Jacobi）矩陣，為則3.3.1等參元剛度矩陣雅可比（Jacobi）矩陣的逆矩陣為將單元應變代入平面問題的物理方程式，就得到平面八節點等參元的應力列陣，為（i=1，2，…，8）3.3.1等參元剛度矩陣等參元剛度矩陣仍可用式（3.30）計算，即在上式中，應變矩陣為局部坐標系ξη下的顯式，因此在計算上式時應將微面積dxdy變換成dξ和dη所包圍的面積dA，由于所以，有3.3.1等參元剛度矩陣則等參元剛度矩陣為應該指出，上式是對ξ和η的重積分，盡管其積分區域十分簡單，但其被積函數卻比較復雜，一般該積分沒有顯式，需要采用數值積分法求解。至于等參元的等效節點力的計算一般也要采用數值積分的方法求解，這里不再討論，可參考相關書籍。由以上分析可知，在劃分單元時，只需確定單元節點的整體坐標值，而不必畫出等參元的具體形狀，因為在計算中實際使用的只有單元八個節點在整體坐標系下的位移值。等參元變換的條件為，因此在有限元網格劃分時，要特別注意這一點。3.3.1等參元剛度矩陣如圖所示，在圖b中3、4點退化為一個點，在該點，同理在圖c中2、3點退化為一個點，在該點，在圖d在點1、2、3處>0，在點4處<0，又因為在單元內連續變換，所以單元內肯定存在＝0，這是由于單元過分歪曲造成的。a）正常b）不正常c）不正常