计算复杂弹性地基板的半解析、半数值方法.pdfVIP

第六届全国岩土力学数值分析与解析方法讨论会，中国，广州 1998年8月4～6日 计算复杂弹性地基板的半解析、半数值方法 王元汉谢华 张佑启 (华中理工大学武授430074) (普蒲大学 香港) 捶薹本文将不同弹性地基上作用集中力的位移沉降公式作为基本解，用解析方法研究了弹性地基与板的相互作用．葡 时用等参有聚元公式分析了地基与板界面接触匣力与爷点力的关系。因此，仅需将弹性地基上的板划分为有限单元．选种 半解析，半数值方法什算量小．精度较高，是分析地基与基础相互作用的有效方法。 燕-谰弹性地基板‘半解析、半数值方法；基本解{等参元法 ，戋￡。； ．}+口《寻 1前 言 弹性地基板的应用十分普遍。通常采用有限元法将板和地基都划分为单元进行计算。由于地蕤为弹 性半空间，计算时往往人为地取一条远处边界。这不仅需大羞计算工作量，计算效果也不理想。 早在60年代，张佑启教授等“3将薄板用矩形单元离散，板单元的每个节点与地基用链杆相联接．用 Boussinesq公式建立了单元节点力与地基反力之间的近似关系。这样仅需将板划分若干矩形单元，按照有 限元法进行计算。 本文进一步推广这种方法于复杂弹性地基板的计算。主要改进之处有：一是采用集中力作用于不同弹 性地基上的解作为基本解，用解析方法计算复杂地基与板的相互作用I二是采用八节点等参元法，用解析 方法研究地基反力与节点力之间的关系，便于不同形状扳的计算{三是采用Mindlin扳理论，同时适用予 薄板与厚板。 2板的弯曲理论 经典薄板的一个基本假设是：平板弯曲后其中面法线仍垂直于中面。Mindlin板理论放弃这一假设， 中面法线不再垂直于弯曲后的平面。这样既方便构造协调单元，又考虑扳的横向剪切变形。同时适合予薄 板与厚板“]。 根据Mindlin板理论，某点位移分量包括#方向挠度w和确定中面法线位置的转角以和日，。倒如，单 元节点位移列阵为 毒一{w，以，口，)! (1) 与i点对应的广义节点力刊阵为 R。一{F，M，眠汀 (2) 弹性基础上的板通常同时作用有上部载荷，和地基反力Q。因此，用板的总刚度矩阵置和总的节点 位移列阵5所表示的平衡方程可表示为叫： Ka=F—Q (3) 如果Q与5的关系已知，例如 Q=K万 (4) 则式(3)可写为 (K+Kr)5一F (5) 80 计算复杂弹性地基板的半解析、半数值方法 王元汉谢华张佑启 解此用位移表示的线性方程组，可得节点位移5，进丽可由式(4)计算板与基础的相互作用，以及板中的成 变、应力和内力。 设单元上荷载p(x。，)与单元节点处荷载P的关系为 P一ⅣP (6) 单元节点力g公式为 g一|1N7pdxdy (7) 将式(6)代凡式(7)，考虑到P可以提到积分式外，则 ∥=CP (8) 其中 C一{JN7Ndxdy (9) 如果地基采用均匀、各向同性弹性半空间模型，根据Boussinesq公式，地基对板单元上作用载荷户引 起节点i处位移为：