整體及總體剛度矩陣的性質(zhì)概述課件

整體分析

及總體剛度矩陣的性質(zhì)整體分析

及總體剛度矩陣的性質(zhì)1整體分析

圖示結(jié)構(gòu)的網(wǎng)格共有四個單元和六個節(jié)點。在節(jié)點1、4、6共有四個支桿支承。結(jié)構(gòu)的載荷已經(jīng)轉(zhuǎn)移為結(jié)點載荷。整體分析的四個步驟：1、建立整體剛度矩陣；2、根據(jù)支承條件修改整體剛度矩陣；3、解方程組，求節(jié)點位移；4、根據(jù)節(jié)點位移求出應(yīng)力。

單元分析得出單元剛度矩陣，下面，將各單元組合成結(jié)構(gòu)，進(jìn)行整體分析。整體分析圖示結(jié)構(gòu)的網(wǎng)格共有四個單元和六個節(jié)點。在節(jié)點12整體及總體剛度矩陣的性質(zhì)概述課件3整體及總體剛度矩陣的性質(zhì)概述課件4整體分析

1、建立整體剛度矩陣(也叫作結(jié)構(gòu)剛度矩陣)上圖中的結(jié)構(gòu)有六個節(jié)點，共有12個節(jié)點位移分量和12個節(jié)點力分量。由結(jié)構(gòu)的節(jié)點位移向量求結(jié)構(gòu)的節(jié)點力向量時，轉(zhuǎn)換關(guān)系為：分塊形式為：其中子向量和都是二階向量，子矩陣是二行二列矩陣。整體剛度矩陣[K]是12*12階矩陣。整體分析1、建立整體剛度矩陣(也叫作結(jié)構(gòu)剛度矩陣)5整體分析

2、根據(jù)支承條件修改整體剛度矩陣。

建立整體剛度矩陣時，每個節(jié)點的位移當(dāng)作未知量看待，沒有考慮具體的支承情況，因此進(jìn)行整體分析時還要針對支承條件加以處理。在上圖的結(jié)構(gòu)中，支承條件共有四個，即在節(jié)點1、4、6的四個支桿處相應(yīng)位移已知為零：建立節(jié)點平衡方程時，應(yīng)根據(jù)上述邊界條件進(jìn)行處理。3、解方程組，求出節(jié)點位移。通常采用消元法和迭代法兩種方法。4、根據(jù)節(jié)點位移求出應(yīng)力。整體分析2、根據(jù)支承條件修改整體剛度矩陣。6整體剛度矩陣的形式

整體剛度矩陣是單元剛度矩陣的集成。1、剛度集成法的物理概念：剛度矩陣中的元素，即由節(jié)點作單位位移時引起的節(jié)點力。在單元剛陣中，表示j節(jié)點單位位移，其他節(jié)點位移為零時，單元e在i節(jié)點引起的節(jié)點力；類似，在整體剛陣中，表示j節(jié)點單位位移，其他節(jié)點位移為零時，整體結(jié)構(gòu)在i節(jié)點引起的節(jié)點力（由于結(jié)構(gòu)已被離散為一系列單元，即所有與i、j節(jié)點相關(guān)的單元在i節(jié)點引起的節(jié)點力之和）。如上圖結(jié)構(gòu)，計算時，與節(jié)點2和3相關(guān)的單元有單元①和③，當(dāng)節(jié)點3發(fā)生單位位移時，相關(guān)單元①和③同時在節(jié)點2引起節(jié)點力，將相關(guān)單元在節(jié)點2的節(jié)點力相加，就得出結(jié)構(gòu)在節(jié)點2的節(jié)點力。由此看出，結(jié)構(gòu)的剛度系數(shù)是相關(guān)單元的剛度系數(shù)的集成，結(jié)構(gòu)剛度矩陣中的子塊是相關(guān)單元的對應(yīng)子塊的集成。整體剛度矩陣的形式整體剛度矩陣是單元剛度矩陣7整體剛度矩陣的形式2、剛度矩陣的集成規(guī)則：1）在整體離散結(jié)構(gòu)變形后，應(yīng)保證各單元在節(jié)點處仍然協(xié)調(diào)地相互連接，即在該節(jié)點處所有單元在該節(jié)點上有相同位移，2）整體離散結(jié)構(gòu)各節(jié)點應(yīng)滿足平衡條件。即環(huán)繞每個節(jié)點的所有單元作用其上的節(jié)點力之和應(yīng)等于作用于該節(jié)點上的節(jié)點載荷Ri，

12i

3412i

Ri34整體剛度矩陣的形式2、剛度矩陣的集成規(guī)則：12i348整體剛度矩陣的形式2、整體剛度矩陣的集成方法具體集成方法是：先對每個單元求出單元剛度矩陣，然后將其中的每個子塊送到結(jié)構(gòu)剛度矩陣中的對應(yīng)位置上去，進(jìn)行迭加之后即得出結(jié)構(gòu)剛度矩陣[K]的子塊，從而得出結(jié)構(gòu)剛度矩陣[K]。關(guān)鍵是如何找出中的子塊在[K]中的對應(yīng)位置。這需要了解單元中的節(jié)點編碼與結(jié)構(gòu)中的節(jié)點編碼之間的對應(yīng)關(guān)系。整體剛度矩陣的形式2、整體剛度矩陣的集成方法9整體剛度矩陣的形式

結(jié)構(gòu)中的節(jié)點編碼稱為節(jié)點的總碼，各個單元的三個節(jié)點又按逆時針方向編為i,j,m,稱為節(jié)點的局部碼。單元剛度矩陣中的子塊是按節(jié)點的局部碼排列的，而結(jié)構(gòu)剛度矩陣中的子塊是按節(jié)點的總碼排列的。因此，在單元剛度矩陣中，把節(jié)點的局部碼換成總碼，并把其中的子塊按照總碼次序重新排列。整體剛度矩陣的形式結(jié)構(gòu)中的節(jié)點編碼稱為節(jié)點的總碼，各10整體剛度矩陣的形式

以單元②為例，局部碼i,j,m對應(yīng)于總碼5,2,4，因此子塊按照總碼重新排列后，得出擴(kuò)大矩陣為：而相應(yīng)的單元剛度方程為（或節(jié)點力表達(dá)式）：整體剛度矩陣的形式以單元②為例，局部碼i,j,m對應(yīng)11整體剛度矩陣的形式用同樣的方法可得出其他單元的擴(kuò)大的單元剛度方程:據(jù)節(jié)點力平衡，各個單元相應(yīng)節(jié)點力疊加：整理可得，整體平衡方程：整體剛度矩陣的形式用同樣的方法可得出其他單元的擴(kuò)大的單元剛度12整體剛度矩陣的形式整體平衡方程：

1）其中[K]為將各單元的擴(kuò)大矩陣迭加所得出的結(jié)構(gòu)剛度矩陣：集成包含搬家和迭加兩個環(huán)節(jié)：A、將單元剛度矩陣中的子塊搬家，得出單元的擴(kuò)大剛度矩陣。B、將各單元的擴(kuò)大剛度矩陣迭加，得出結(jié)構(gòu)剛度矩陣[K]。2）為節(jié)點載荷向量，為節(jié)點位移向量。整體剛度矩陣的形式整體平衡方程：13整體及總體剛度矩陣的性質(zhì)概述課件142.整體剛度矩陣的特點

在有限元法中，整體剛度矩陣的階數(shù)通常是很高的，在解算時常遇到矩陣階數(shù)高和存貯容量有限的矛盾。找到整體剛度矩陣的特性達(dá)到節(jié)省存貯容量的途徑。1、對稱性。只存貯矩陣的上三角部分，節(jié)省近一半的存貯容量。2、稀疏性。矩陣的絕大多數(shù)元素都是零，非零元素只占一小部分。2.整體剛度矩陣的特點15整體剛度矩陣的特點

2、稀疏性。矩陣的絕大多數(shù)元素都是零，非零元素只占一小部分。節(jié)點5只與周圍的六個節(jié)點(2、3、4、6、8、9)用三角形單元相連，它們是5的相關(guān)節(jié)點。只有當(dāng)這七個相關(guān)節(jié)點產(chǎn)生位移時，才使該節(jié)點產(chǎn)生節(jié)點力，其余節(jié)點發(fā)生位移時并不在該節(jié)點處引起節(jié)點力。因此，在矩陣[K]中，第5行的非零子塊只有七個(即與相關(guān)節(jié)點對應(yīng)的七個子塊)。整體剛度矩陣的特點2、稀疏性。節(jié)點5只與周圍的16整體剛度矩陣的特點

2、稀疏性。

一般，一個節(jié)點的相關(guān)結(jié)點不會超過九個，如果網(wǎng)格中有200個節(jié)點，則一行中非零子塊的個數(shù)與該行的子塊總數(shù)相比不大于9/200，即在5%以下，如果網(wǎng)格的節(jié)點個數(shù)越多，則剛度矩陣的稀疏性就越突出。利用矩陣[K]的稀疏性，可設(shè)法只存貯非零元素，從而可大量地節(jié)省存貯容量。整體剛度矩陣的特點2、稀疏性。一般，一個節(jié)點17整體剛度矩陣的特點

3、帶形分布規(guī)律。上圖中，矩陣[K]的非零元素分布在以對角線為中心的帶形區(qū)域內(nèi)，稱為帶形矩陣。在半個帶形區(qū)域中(包括對角線元素在內(nèi))，每行具有的元素個數(shù)叫做半帶寬，用d表示。半帶寬的一般計算公式是：半帶寬d=(相鄰結(jié)點碼的最大差值+1)*2上圖中相鄰節(jié)點碼的最大差值為4，故d=(4+1)*2=10利用帶形矩陣的特點并利用對稱性，可只存貯上半帶的元素，叫半帶存貯。

整體剛度矩陣的特點3、帶形分布規(guī)律。18整體剛度矩陣的特點

圖(a)中的矩陣[K]為n行n列矩陣，半帶寬為d。半帶存貯時從[K]中取出上半帶元素，按圖(b)中的矩陣的排列方式進(jìn)行存貯，即將上半部斜帶換成豎帶。存貯量n*d，存貯量與[K]中元素總數(shù)之比為d/n，d值越小，則存貯量約省。矩陣[K]矩陣對角線第1列r行r行r列45度斜線r行s列r行s-r+1列元素元素整體剛度矩陣的特點圖(a)中的矩陣[K]為n行n列矩19整體剛度矩陣的特點

同一網(wǎng)格中，如果采用不同的節(jié)點編碼，則相應(yīng)的半帶寬d也可能不同。如圖，是同一網(wǎng)格的三種節(jié)點編碼，相鄰節(jié)點碼的最大差值分別為4、6、8，半帶寬分別為10、14、18。因此，應(yīng)當(dāng)采用合理的節(jié)點編碼方式，以便得到最小的半帶寬，從而節(jié)省存貯容量。整體剛度矩陣的特點同一網(wǎng)格中，如果采用不同的節(jié)點編碼20整體分析

及總體剛度矩陣的性質(zhì)整體分析

及總體剛度矩陣的性質(zhì)21整體分析

圖示結(jié)構(gòu)的網(wǎng)格共有四個單元和六個節(jié)點。在節(jié)點1、4、6共有四個支桿支承。結(jié)構(gòu)的載荷已經(jīng)轉(zhuǎn)移為結(jié)點載荷。整體分析的四個步驟：1、建立整體剛度矩陣；2、根據(jù)支承條件修改整體剛度矩陣；3、解方程組，求節(jié)點位移；4、根據(jù)節(jié)點位移求出應(yīng)力。

單元分析得出單元剛度矩陣，下面，將各單元組合成結(jié)構(gòu)，進(jìn)行整體分析。整體分析圖示結(jié)構(gòu)的網(wǎng)格共有四個單元和六個節(jié)點。在節(jié)點122整體及總體剛度矩陣的性質(zhì)概述課件23整體及總體剛度矩陣的性質(zhì)概述課件24整體分析

1、建立整體剛度矩陣(也叫作結(jié)構(gòu)剛度矩陣)上圖中的結(jié)構(gòu)有六個節(jié)點，共有12個節(jié)點位移分量和12個節(jié)點力分量。由結(jié)構(gòu)的節(jié)點位移向量求結(jié)構(gòu)的節(jié)點力向量時，轉(zhuǎn)換關(guān)系為：分塊形式為：其中子向量和都是二階向量，子矩陣是二行二列矩陣。整體剛度矩陣[K]是12*12階矩陣。整體分析1、建立整體剛度矩陣(也叫作結(jié)構(gòu)剛度矩陣)25整體分析

2、根據(jù)支承條件修改整體剛度矩陣。

建立整體剛度矩陣時，每個節(jié)點的位移當(dāng)作未知量看待，沒有考慮具體的支承情況，因此進(jìn)行整體分析時還要針對支承條件加以處理。在上圖的結(jié)構(gòu)中，支承條件共有四個，即在節(jié)點1、4、6的四個支桿處相應(yīng)位移已知為零：建立節(jié)點平衡方程時，應(yīng)根據(jù)上述邊界條件進(jìn)行處理。3、解方程組，求出節(jié)點位移。通常采用消元法和迭代法兩種方法。4、根據(jù)節(jié)點位移求出應(yīng)力。整體分析2、根據(jù)支承條件修改整體剛度矩陣。26整體剛度矩陣的形式

整體剛度矩陣是單元剛度矩陣的集成。1、剛度集成法的物理概念：剛度矩陣中的元素，即由節(jié)點作單位位移時引起的節(jié)點力。在單元剛陣中，表示j節(jié)點單位位移，其他節(jié)點位移為零時，單元e在i節(jié)點引起的節(jié)點力；類似，在整體剛陣中，表示j節(jié)點單位位移，其他節(jié)點位移為零時，整體結(jié)構(gòu)在i節(jié)點引起的節(jié)點力（由于結(jié)構(gòu)已被離散為一系列單元，即所有與i、j節(jié)點相關(guān)的單元在i節(jié)點引起的節(jié)點力之和）。如上圖結(jié)構(gòu)，計算時，與節(jié)點2和3相關(guān)的單元有單元①和③，當(dāng)節(jié)點3發(fā)生單位位移時，相關(guān)單元①和③同時在節(jié)點2引起節(jié)點力，將相關(guān)單元在節(jié)點2的節(jié)點力相加，就得出結(jié)構(gòu)在節(jié)點2的節(jié)點力。由此看出，結(jié)構(gòu)的剛度系數(shù)是相關(guān)單元的剛度系數(shù)的集成，結(jié)構(gòu)剛度矩陣中的子塊是相關(guān)單元的對應(yīng)子塊的集成。整體剛度矩陣的形式整體剛度矩陣是單元剛度矩陣27整體剛度矩陣的形式2、剛度矩陣的集成規(guī)則：1）在整體離散結(jié)構(gòu)變形后，應(yīng)保證各單元在節(jié)點處仍然協(xié)調(diào)地相互連接，即在該節(jié)點處所有單元在該節(jié)點上有相同位移，2）整體離散結(jié)構(gòu)各節(jié)點應(yīng)滿足平衡條件。即環(huán)繞每個節(jié)點的所有單元作用其上的節(jié)點力之和應(yīng)等于作用于該節(jié)點上的節(jié)點載荷Ri，

12i

3412i

Ri34整體剛度矩陣的形式2、剛度矩陣的集成規(guī)則：12i3428整體剛度矩陣的形式2、整體剛度矩陣的集成方法具體集成方法是：先對每個單元求出單元剛度矩陣，然后將其中的每個子塊送到結(jié)構(gòu)剛度矩陣中的對應(yīng)位置上去，進(jìn)行迭加之后即得出結(jié)構(gòu)剛度矩陣[K]的子塊，從而得出結(jié)構(gòu)剛度矩陣[K]。關(guān)鍵是如何找出中的子塊在[K]中的對應(yīng)位置。這需要了解單元中的節(jié)點編碼與結(jié)構(gòu)中的節(jié)點編碼之間的對應(yīng)關(guān)系。整體剛度矩陣的形式2、整體剛度矩陣的集成方法29整體剛度矩陣的形式

結(jié)構(gòu)中的節(jié)點編碼稱為節(jié)點的總碼，各個單元的三個節(jié)點又按逆時針方向編為i,j,m,稱為節(jié)點的局部碼。單元剛度矩陣中的子塊是按節(jié)點的局部碼排列的，而結(jié)構(gòu)剛度矩陣中的子塊是按節(jié)點的總碼排列的。因此，在單元剛度矩陣中，把節(jié)點的局部碼換成總碼，并把其中的子塊按照總碼次序重新排列。整體剛度矩陣的形式結(jié)構(gòu)中的節(jié)點編碼稱為節(jié)點的總碼，各30整體剛度矩陣的形式

以單元②為例，局部碼i,j,m對應(yīng)于總碼5,2,4，因此子塊按照總碼重新排列后，得出擴(kuò)大矩陣為：而相應(yīng)的單元剛度方程為（或節(jié)點力表達(dá)式）：整體剛度矩陣的形式以單元②為例，局部碼i,j,m對應(yīng)31整體剛度矩陣的形式用同樣的方法可得出其他單元的擴(kuò)大的單元剛度方程:據(jù)節(jié)點力平衡，各個單元相應(yīng)節(jié)點力疊加：整理可得，整體平衡方程：整體剛度矩陣的形式用同樣的方法可得出其他單元的擴(kuò)大的單元剛度32整體剛度矩陣的形式整體平衡方程：

1）其中[K]為將各單元的擴(kuò)大矩陣迭加所得出的結(jié)構(gòu)剛度矩陣：集成包含搬家和迭加兩個環(huán)節(jié)：A、將單元剛度矩陣中的子塊搬家，得出單元的擴(kuò)大剛度矩陣。B、將各單元的擴(kuò)大剛度矩陣迭加，得出結(jié)構(gòu)剛度矩陣[K]。2）為節(jié)點載荷向量，為節(jié)點位移向量。整體剛度矩陣的形式整體平衡方程：33整體及總體剛度矩陣的性質(zhì)概述課件342.整體剛度矩陣的特點

在有限元法中，整體剛度矩陣的階數(shù)通常是很高的，在解算時常遇到矩陣階數(shù)高和存貯容量有限的矛盾。找到整體剛度矩陣的特性達(dá)到節(jié)省存貯容量的途徑。1、對稱性。只存貯矩陣的上三角部分，節(jié)省近一半的存貯容量。2、稀疏性。矩陣的絕大多數(shù)元素都是零，非零元素只占一小部分。2.整體剛度矩陣的特點35整體剛度矩陣的特點

2、稀疏性。矩陣的絕大多數(shù)元素都是零，非零元素只占一小部分。節(jié)點5只與周圍的六個節(jié)點(2、3、4、6、8、9)用三角形單元相連，它們是5的相關(guān)節(jié)點。只有當(dāng)這七個相關(guān)節(jié)點產(chǎn)生位移時，才使該節(jié)點產(chǎn)生節(jié)點力，其余節(jié)點發(fā)生位移時并不在該節(jié)點處引起節(jié)點力。因此，在矩陣[K]中，第5行的非零子塊只有七個(即與相關(guān)節(jié)點對應(yīng)的七個子塊)。整體剛度矩陣的特點2、稀疏性。節(jié)點5只與周圍的36整體剛度矩陣的特點

2、稀疏性。

一般，一個節(jié)點的相關(guān)結(jié)點不會超過九個，如果網(wǎng)格中有200個節(jié)點，則一行中非零子塊的個數(shù)與該行的子塊總數(shù)相比不大于9/200，即在5%以下，如果網(wǎng)格的節(jié)點個數(shù)越多，則剛度矩陣的稀疏性就越突出。利用矩陣[K]的稀疏性，可設(shè)法只存貯非零元素，從而可大量地節(jié)省存貯容量。整體剛度矩陣的特點