用M I D A S模拟桩土相
互作用
-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN
用MIDAS模拟桩-土相互作用（“m法”确定土弹簧刚度）
北京迈达斯技术有限公司
2009年05月
2
1、引言
土与结构相互作用的研究已有近60～70年的历史，待别是近30年来，计算机技术的发展为其提供了有力的分析手段。桩基础是土建工程中广泛采用的基础形式之一，许多建于软土地基上的大型桥梁结构往往都采用桩基础，桩-土动力相互作用又是土-结构相互作用问题中较复杂的课题之一。至今已有不少关于桩基动力特性的研究报告，国内外研究人员也提出了许多不同的桩-土动力相互作用计算方法。从研究成果的归类来看，理论上主要有离散理论和连续理论及两者的结合，解决的方法一般有集中质量法、有限元法、边界元法和波动场法。
60～70年代，美国学者J.penzien等在解决泥沼地上大桥动力分析时提出了集中质量法，目前已在国内外得到了广泛的应用。集中质量法将桥梁上部结构多质点体系和桩一土体系的质量联合作为一个整体，来建立整体耦联的地震振动微分方程组进行求解。该模型假定桩侧土是Winkler连续介质。以半空间的Mindlin静力基本解为基础，将桩-土体系的质量按一定的厚度简化并集中为一系列质点，离散成一理想化的参数系统。并用弹簧和阻尼器模拟土介质的动力性质，形成一个包括地下部分的多质点体系。
土弹簧刚度的确定，除考虑使用较为精确的有限元或边界元方法外，较为简便的方法是采用Penzien模型中提供的土弹簧计算方法或参照现行规范中土弹簧的计算方法。我国公路桥涵地基与基础设计规范(JTG D63-2007)用的“m 法”计算方法和参数选取方面比Penzien的方法要简单和方便，且为国内广大
3
工程师所熟.“m法”的基本原理是将桩作为弹性地基梁，按Winkler假定（梁身任一点的土抗力和该点的位移成正比）求解。但是，由于桩-土相互作用的实验数据不足，土的物性取值有时亦缺乏合理性，在确定土弹簧的刚度时，仍有不少问题未能很好解决。特别是，“m法”中m的取值对弹簧刚度的计算结果影响很大，且不能反映地震波的频率特性和强度带来的影响。
本次介绍的土弹簧的模拟是采用规范中的“m法”确定土的地基系数C（m 的取值根据土的物性而定）,再由其算出土弹簧的水平刚度。?
2、采用“m”法，确定土弹簧刚度
桩在横轴向荷载作用下桩身的内力和位移计算，国内外学者提出了许多方法。目前较为普遍的是桩侧土采用文克尔假定，通过求解挠曲微分方程，再结合力的平衡条件，求出桩各部位的内力和位移，该方法称为弹性地基梁法。
以文克尔假定为基础的弹性地基梁法从土力学观点看是不够严密的，但其基本概念明确，方法简单，所得结果一般较安全，在国内外工程界得到广泛应用。我国公路、铁路在桩基础的设计中常用的“m”法、就属此种方法，本文将主要根据“m”法来确定土弹簧刚度。
2.1地基系数的概念及确定方法
（1）概念
4
5
地基系数C 表示单位面积土在弹性限度内产生单位变形时所需施加的力，单位为kN/m 3或MN/m 3。
（2）确定方法
地基系数大小与地基土的类别、物理力学性质有关。
地基系数C 值是通过对试桩在不同类别土质及不同深度进行实测z x 及zx 后反算得到。大量的试验表明，地基系数C 值不仅与土的类别及其性质有关，而且也随着深度而变化。由于实测的客观条件和分析方法不尽相同等原因，所采用的C 值随深度的分布规律也各有不同。常采用的地基系数分布规律有图下所示的几种形式，因此也就产生了与之相应的基桩内力和位移的计算方法。
图1 地基系数变化规律
现将桩的几种有代表性的弹性地基梁计算方法概括在表下中。
桩的几种典型的弹性地基梁法 表4-1
计算方法
图 号 地基系数随深度分布 地基系数C 表达式 说 明 m 法 4-1a ） 与深度成正比 C =mZ m 为地基土比例系数
上述的四种方法各自假定的地基系数随深度分布规律不同，其计算结果是有差异的。实验资料分析表明，宜根据土质特性来选择恰当的计算方法。
2.2桩的计算宽度
1．定义
计算桩的内力与位移时不直接采用桩的设计宽度（直径），而是换算成实际工作条件下相当于矩形截面桩的宽度b1，b1称为桩的计算宽度。
2．采用计算宽度的原因
为了将空间受力简化为平面受力，并综合考虑桩的截面形状及多排桩桩间的相互遮蔽作用。
3．计算方法
桩的计算宽度b1请按照现行桥涵地基与基础设计规范(JTG D63-
2007)P83~P84页要求确定。
以上的计算方法比较复杂，理论和实践的根据也是不够的，因此国内有些规范建议简化计算。圆形桩：当d≤1m时，b1=0.9（1.5d+0.5）；当d＞1m 时，b1=0.9（d+1）。方形桩：当边宽b≤1m时，b1=1.5b+0.5；当边宽＞1m 时，b1=b+1。而国外有些规范更为简单：柱桩及桩身尺寸直径0.8m以下的灌注桩，b1=d+1（m）；其余类型及截面尺寸的桩，b1=1.5d+0.5（m）。
2.3刚性桩与弹性桩
6