弹性地基梁单元刚度矩阵及固端内力的求解方法研究.pdfVIP

第九届全国建筑工程计算机应用学术会议论文集 黑龙江 ·哈尔滨 1998 弹性地基梁单元刚度矩阵及固端内力的求解方法 舒 宣 武 (华南理工大学建筑设计研究院) 【摘要】 本文引进弹性地基梁单元特性矩阵和积分常数矩阵的概念，提出了一种计算弹性地荃梁单元 刚度矩阵的方法以及求其上作用均布荷载的固端内力的方法，并给出了相应的FORTRAN语言程序。本 文给出的计算机程序，可直接用于有限元的结构分析程序中。 一、 引言 弹性地基梁的解析解法和实用的计算表格在很多文献中都有详细的介绍}1]。不过，这 种解法只适用于手算，且计算工作量很大，解题范围有限。若要对任意弹性地基梁系进行 结构分析，特别是考虑地基基础同上部结构共同作用时，用查表的方法进行计算是根本行 不通的。如果能象一般有限元法一样，用简便的方法得到弹性地基梁单元的单元刚度矩阵， 问题就迎刃而解了。 本文给出了文克尔(Winkler)弹性地基梁单元刚度矩阵以及其上作用均布荷载的固端 内力的实用算法，并用FORTRAN语言编写了程序。上机计算表明，效果很好。 二、弹性地基梁的单元刚度矩阵 图1弹性地基梁单元 弹性地基梁单元如图1所示，其刚度矩阵[K]定义为: I l ‘ 戈 e s l l 仁 气 尤 . 尤 . 4 1 l 厂 . . . . △ 1 、 ! 一 十 w e l l 从 w e l 气 暇 k k 4 2 e s e s e s △ ， ‘ 了 、 一 一 e s l ‘ ! 1 丛 o t l 凡 札 k k 4 3 e s e s [ △ 、 ， e s l 1 凡 J 气 札 k k 4 4 J l △ 4 或 [A]=[幻{A} 文克尔地基上的弹性地基梁挠曲微分方程为1【] 第九届全国建筑工程计算机应用学术会议论文集 黑龙江 ·哈尔滨 1998 dy(x) — 十 4刀4Y(x)=0 dc 上式中 、一n4E1一a4kEB1 为弹性地基梁的柔度特征值，其中，k为地基土的基床系数，B为基础底板宽，K=kB为 集中基床系数。 式2的解可写成如 厂形式: Y(x)=艺a,H,(ftc) 上式中， 仁 a 拭(口=ofsink J l b 从(动=eecosh f 气 C 丛(动=esink ︺