弹性地基杆系有限元法单元刚度的若干讨论研究.pdfVIP

第六届全国岩土力学数值分析与解析方法讨论会，中国，广州 1998年8月4^6F 弹性地基杆系有限元法单元刚度的若干讨论 杨 敏 艾智勇 张永谋 {同济大学地下建筑与工程系 上海 200092) 握要本文对弹性地基杆系有限单元法中的几种杆单元刚度矩阵作了讨论，得出一些有益的结论 关.词 弹性地基，杆系有限元;单元刚度矩阵 ‘、 /r, 1 引 言 弹性地基杆系有限单元法目前已成为深基坑工程中柔性支挡结构计算内力和变形的主要方法川，同 时也是受水平荷载桩和支撑在地基上的梁等结构物的主要方法。然而人们运用该法进行计算时所取的单 元刚度矩阵往往有所不同。这样自然会使最终计算结果有所差异，本文试图对这一问题进行讨论，以期得 到一些有益的结论。 Z 单元刚度矩阵的推导 2.1 甚床系数K为常数时的单刚 为了得到弹性地基梁的杆单元刚度矩阵，常常有以下几种处理方法:一种方法是将基床系数为K的 分布弹赞简化为梁单元的结点弹簧，应用结构力学的方法得到单刚;另一种方法是假定满足边界条件的形 函数，直接运用虚功原理推导而得;第三种方法是利用弹性地基梁挠曲线方程的解析解而得。 Y,Y Y,Y 图 1 弹性地蚤梁 图2 简化为结点弹簧 2.1.1 分布基床系致简化为结点弹获(荟床系教)时的单元刚度矩阵 如图1，设单元长为l,梁宽为B，基床系数为K，则图2中的简化结点弹赞刚度为K, K，一音·“·‘ (1) 单元节点力为 {F`}二 {Y;从 Y，从}r 而单元节点位移为 万焚丫丁 弹性地基杆系有限元法单元刚度的若干讨论 杨 饭 艾搜舅 张永谋 {子}二{Y 坎 y 即T 节点力与节点位移之间的关系为 tF)=[K]{6-) “(2) 按照结构力学方法获得单刚及总刚的方法可参见一般的结构力学教材或土力学及基础工程学教材 下面运用虚功原理来推导单元刚度矩阵 K〔〕 下面运用虚功原理来推导[K] 单元内任一点的挠度Y可表示为;z7 f 飞 y i . 1 . ! 伙 . I Y二 N「,NoN3N 今 J 飞 才 、 1 广 t 1 ~ ~ N)][{ (3) 1 为 1 ‘ 毛 . 色 J J J J 其中 N 片 1一 Nz=x一州 一 1 + 一 尸 尸 h3x+2Lx} N 3x 2x