常应变三角形单元课件.pptVIP

不是對於局部座標的，而是對於整體座標的，需要轉換。是沿著一個空間曲面進行積分，也需要轉換成局部座標。常應變三角形單元單元剛度矩陣形成單元本身是平衡的[k]δ:位移引起的結點力,其中[k]為相應的剛度。主元：結點本身位移引起的本結點的結點力輔元：結點2的位移引起結點1的結點力有了剛度值，就可以由位移求出力（非應力）有了彈性模量值，就可以由應變求出應力σεpδ斜率為彈性模量斜率為剛度彈性地基粱為例一般彈性體為例整體各個單元剛度矩陣分析124356①②③④總體剛度矩陣的疊加124356①②③④表示：整體結構中，整體結點編號為1的結點位移引起的本身的結點力，相應的剛度表示：整體結構中，整體結點編號為2的結點位移引起結點1的結點力，相應的剛度表示：整體結構中，整體結點編號為5的結點位移引起結點3的結點力，相應的剛度整體剛度矩陣的性質整剛表示整體位移與力的關係總體剛度矩陣的行數和列數＝總結點數×單元結點自由度數的方陣，其排列的順序按照結點標號依次從小到大排列，每個結點的按照自由度的順序排列。上例中總剛度矩陣排列見後，這樣排列是為了與單元結點位移向量對應。結點的多少決定著整剛的大小。約束處理就是將整剛中相應的約束結點的約束自由度上的行、列劃去，將荷載向量中相應項劃去即可。其實質是：強制使得u＝0，或v＝0；或u＝0，v=0對於給定位移約束條件，即：u＝c1或v＝c2 是將整體剛度矩陣中相應的主元乘以一個大值A，等效結點向量相應的值設為A×c1或A×c2，利用“大數吃小數”的原理，得到u＝c1或v＝c2。在上例中，若整體結點號為2的水準位移u2＝c2，則有：相對於A很小§4.5等參元以上三角形單元分析中可以看出，從位移模式到形函數的構造上面，需要解三元一次方程，若位移模式中取得的項數進一步增加到4～20項，其形函數的構造非常麻煩，而且不能夠統一編程，因此需要找更方便的辦法在三角形單元中，對於等效結點荷載向量的求解等就很複雜了，對於其他複雜單元其求解將更加困難。這就引出來另外一種思路。即：先將實際單元通過座標轉換函數轉化成一個母元（即：形狀規則、簡單的單元），所有位移模式、形函數、單元分析，等效結點荷載向量建立等均在該母元上進行，使得求解簡單化，並便於變成。然後將求出的量再轉化成實際單元的量，求出最後的結果。這種思路概念清楚、便於理解，並且便於標準化。這就是等參元的思路。 其中：座標轉換函數＝母元的形函數時稱為等參元。實際單元母單元母單元實際單元2×22×2×2等參元的形函數母元對於平面問題空間8結點等參元母元空間20結點等參元母元母元位移模式與形函數的構造位移模式的待定參數的個數應該與結點個數一致，才能保證待定參數可以解出來。對於4邊形單元，位移模式取為（在局部座標中分析）：將母元各個結點的局部座標值代入位移模式中，得到：以上公式就是平面4結點等參元的形函數對於空間8結點得到其形函數為： 有些專著中還提出8－20結點單元的通用形函數公式。座標轉換函數的建立座標轉換函數的作用在於將幾何形狀不規則的真實單元轉換成為標準單元（數學上稱為映射），等參元是將座標轉換函數與形函數一致，即：形函數N座標轉換函數M等參元中：Ni＝Mi我們分析一下形函數的性質：即：i結點的形函數Ni在j結點處為0在單元邊界上位移是局部座標的一次函數，線性的。座標轉換函數形式：正是由於形函數有以上的性質才能夠使得其可以將不規則的4邊形轉換成為一個2×2的正方形。 我們不妨驗證一下：對於母元而言，其3－4邊的表示方程為η=1。 將η=1（3-4邊）代入到座標轉換方程： 得到：化解得到：將聯立方程中ζ消去得到：直線方程已知兩點座標，求連接兩點直線的方程這是連接（y3,z3）和（y4,z4）兩點的直線。——這說明母元3－4邊的方程η=1，經過座標轉換函數轉換以後，得到在實際座標下連接3－4點的直線方程。同理也可以證明實際座標的直線邊，可以轉化成母元的邊。從而驗證了座標轉換函數可以實現將實際不規則單元轉換為母元的功能。雅克比矩陣(Jacobi)[J]就是真實座標對母元局部座標的導數矩陣對於平面問題：對於空間問題：單元類型不同，其Jacobi矩陣由於座標轉換函數不同而不同。Jacobi矩陣地意義得到：[J]表示座標轉換函數對局部座標的導數與座標轉換函數對整體座標的導數的關係矩陣那麼： Jacobi矩陣的逆陣就表示：座標轉換函數對整體座標的導數與座標轉換函數對局部座標的導數的關係矩陣。 在計算[B]