薄壁构件刚度矩阵不同算法的影响.pdfVIP

口。褊肛羽k(shk-k)Ⅲ卢=譬。盖· L y2三i：I二．三：_至了’七2卜尚’… 即专 观察这两个矩阵，可以发现除了扭转I吹和翘曲项不同外，矩阵中其余的元素完全相同。这是冈为矩阵 位移法中的单元刚度矩阵方程是位移法中的杆端力与杆端位移间的关系式，是精确的表达式。有限元法姓 …种近似的数值解法，单元的刚度方群是按某种假定的位移状态推导出的，是近似的表达式。而二：次式的 f节移模式楚横向弯曲问题的精确何移模式．故l}j此似移模式得剑的是精确解：一：次式位移模式对扭转扛d题 足是～个近似的忙移摊式．故t-tj此伊移模式．得到的绐粜缝近似解，10精确解不同。 3 精度比较 在编制程序进行矩阵计算时我ffJ发现川有限元法推导出来的单元刚度矩阵抓转项和翘曲项都为止值， 而删位移法求出的单元刚度矩阵扭转项和翘曲项订u丁能存在负值，冈此根据出现的现象来讨论两个刚皮矩阵 ……黝口2耥…崔。鑫一商…㈧鼋 厂足一盈，z竺-0 l 2 所以只需下面三个不等式联立求k的范罔 L七一thk-0 项进行汁算． Zk)=妣一撤 躲等，(㈣‘-1_f筹卜l埘¨础)ffi-2thkfl-th (thk)”=6th2kO—th2k)一20一th二k)=一6th4七+8砌二k一2． O庇七)‘4’=一24th3k(1一th：七)+16thk(1一th2k)=24th5k一40th3k+16thk (thk)‘5’=120th4kO—th2k)一120th2七(1一th2k)+1 5 闪为th(。)=。，所以thk住零点的泰勒级数腱开为：砌七=七一三k-+215k ·I．288· kS．EPTIi等式方程组简化为 冈为sh(。)=。，ch(。)=1，所以shk在零点的泰勒级数展开为：妫足=七+苔1七3+-1-21-o-1 珥 矿 斗 ☆ 卜 ◇ 尼一2 1一∞一nu k2_10 一 辛 j + 一2旷 1“七 卜O 后_o， o尼—．41—0． 1—6 一) 2 七：．三 厂●●●√、L 一 ．．一佗一1 -． ∥一 ；r 卜0 @ ． 2 扣嗉扣 2一” 冈为k=L 1 义因为G=丽E丽，巨=而E2(1+∥) 1一∥2，所以笔一F1≥L2。只有在这个范围内进行计算才能保证由上面两(一∥) ，。 个刚度矩阵计算出来的结果相互吻合。 4 结论 (1)所谓有限元法的刚度矩阵可以视作所谓矩阵位移法的多项式逼进，所以，用有限元法的刚度矩阵 可以逼近矩阵位移法的刚度矩阵。 (2)使用有限元法的刚度矩阵时考虑扭转的薄螭杆什的单元长度应该有所限制，这一点与不考虑士H转 的薄肇杆件有明显差别