用MIDAS模拟桩土相互作用

1、用MIDAS模拟桩-土相互作用（“m法”确定土弹簧刚度）北京迈达斯技术有限公司2009年05月1、引言土与结构相互作用的研究已有近6070年的历史，待别是近30年来，计算机技术的发展为其提供了有力的分析手段。桩基础是土建工程中广泛采用的基础形式之一，许多建于软土地基上的大型桥梁结构往往都采用桩基础，桩-土动力相互作用又是土-结构相互作用问题中较复杂的课题之一。至今已有不少关于桩基动力特性的研究报告，国内外研究人员也提出了许多不同的桩-土动力相互作用计算方法。从研究成果的归类来看，理论上主要有离散理论和连续理论及两者的结合，解决的方法一般有集中质量法、有限元法、边界元法和波动场法。6070年代，

2、美国学者J.penzien等在解决泥沼地上大桥动力分析时提出了集中质量法，目前已在国内外得到了广泛的应用。集中质量法将桥梁上部结构多质点体系和桩一土体系的质量联合作为一个整体，来建立整体耦联的地震振动微分方程组进行求解。该模型假定桩侧土是Winkler连续介质。以半空间的Mindlin静力基本解为基础，将桩-土体系的质量按一定的厚度简化并集中为一系列质点，离散成一理想化的参数系统。并用弹簧和阻尼器模拟土介质的动力性质，形成一个包括地下部分的多质点体系。土弹簧刚度的确定，除考虑使用较为精确的有限元或边界元方法外，较为简便的方法是采用Penzien模型中提供的土弹簧计算方法或参照现行规范中土弹簧的

3、计算方法。我国公路桥涵地基与基础设计规范(JTG D63-2007)用的“m法”计算方法和参数选取方面比Penzien的方法要简单和方便，且为国内广大工程师所熟. “m法”的基本原理是将桩作为弹性地基梁，按Winkler假定（梁身任一点的土抗力和该点的位移成正比）求解。但是，由于桩-土相互作用的实验数据不足，土的物性取值有时亦缺乏合理性，在确定土弹簧的刚度时，仍有不少问题未能很好解决。特别是，“m法”中m的取值对弹簧刚度的计算结果影响很大，且不能反映地震波的频率特性和强度带来的影响。本次介绍的土弹簧的模拟是采用规范中的“m法”确定土的地基系数C（ m的取值根据土的物性而定）,再由其算出土弹簧的

4、水平刚度。 2、采用“m”法，确定土弹簧刚度桩在横轴向荷载作用下桩身的内力和位移计算，国内外学者提出了许多方法。目前较为普遍的是桩侧土采用文克尔假定，通过求解挠曲微分方程，再结合力的平衡条件，求出桩各部位的内力和位移，该方法称为弹性地基梁法。以文克尔假定为基础的弹性地基梁法从土力学观点看是不够严密的，但其基本概念明确，方法简单，所得结果一般较安全，在国内外工程界得到广泛应用。我国公路、铁路在桩基础的设计中常用的“m”法、就属此种方法，本文将主要根据“m”法来确定土弹簧刚度。2.1地基系数的概念及确定方法 （1）概念地基系数C表示单位面积土在弹性限度内产生单位变形时所需施加的力，单位为

5、kN/m3或MN/m3。（2）确定方法地基系数大小与地基土的类别、物理力学性质有关。地基系数C值是通过对试桩在不同类别土质及不同深度进行实测及后反算得到。大量的试验表明，地基系数C值不仅与土的类别及其性质有关，而且也随着深度而变化。由于实测的客观条件和分析方法不尽相同等原因，所采用的C值随深度的分布规律也各有不同。常采用的地基系数分布规律有图下所示的几种形式，因此也就产生了与之相应的基桩内力和位移的计算方法。图1 地基系数变化规律现将桩的几种有代表性的弹性地基梁计算方法概括在表下中。 桩的几种典型的弹性地基梁法 表4-1计算方法图 号地基系数随深度分布地基系数C表达式说 明m法4-1a）与深度

6、成正比C=mZm为地基土比例系数K法4-1b）桩身第一挠曲零点以上抛物线变化，以下不随深度变化C=KK为常数C值法4-1c）与深度呈抛物线变化C=cZ0.5c为地基土比例系数张有龄法4-1d）沿深度均匀分布C=K0K0为常数上述的四种方法各自假定的地基系数随深度分布规律不同，其计算结果是有差异的。实验资料分析表明，宜根据土质特性来选择恰当的计算方法。2.2桩的计算宽度1定义计算桩的内力与位移时不直接采用桩的设计宽度（直径），而是换算成实际工作条件下相当于矩形截面桩的宽度b1，b1称为桩的计算宽度。2采用计算宽度的原因为了将空间受力简化为平面受力，并综合考虑桩的截面形状及多排桩桩间的相互遮蔽作用

7、。3计算方法桩的计算宽度b1请按照现行桥涵地基与基础设计规范(JTG D63-2007)P83P84页要求确定。以上的计算方法比较复杂，理论和实践的根据也是不够的，因此国内有些规范建议简化计算。圆形桩：当d1m时，b1=0.9（1.5d+0.5）；当d1m时，b1=0.9（d+1）。方形桩：当边宽b1m时，b1=1.5b+0.5；当边宽1m时，b1=b+1。而国外有些规范更为简单：柱桩及桩身尺寸直径0.8m以下的灌注桩，b1=d+1（m）；其余类型及截面尺寸的桩，b1=1.5d+0.5（m）。2.3刚性桩与弹性桩为计算方便起见，按照桩与土的相对刚度，将桩分为刚性桩和弹性桩。1弹性桩当桩的入土深

8、度时，这时桩的相对刚度小，必须考虑桩的实际刚度，按弹性桩来计算。其中称为桩的变形系数，2刚性桩当桩的入土深度时，则桩的相对刚度较大，计算时认为属刚性桩，2.43土的地基比例系数“m”地基土水平抗力系数的比例系数m值宜通过桩的水平静载试验确定。但由于试验费用、时间等原因，某些建筑物不一定进行桩的水平静载试验，可采用规范提供的经验值如下表所示。 在应用上表时应注意以下事项图4-5 比例系数m的换算1由于桩的水平荷载与位移关系是非线性的，即m值随荷载与位移增大而有所减小，因此，m值的确定要与桩的实际荷载相适应。一般结构在地面处最大位移不超过10mm，对位移敏感的结构、桥梁工程为6mm。位移较大时，应

9、适当降低表列m值。2当基桩侧面由几种土层组成时，从地面或局部冲刷线起，应求得主要影响深度hm=2（d+1）米范围内的平均m值作为整个深度内的m值（见图4-5），对于刚性桩，hm采用整个深度h。当基础侧面地面或局部冲刷线以下当hm（hm>2.5m时）深度内存在两层不同土时： (4-5)当hm深度内存在三层不同土时： （4-6）3承台侧面地基土水平抗力系数Cn Cn=m·hn （4-7）式中：m承台埋深范围内地基土的水平抗力系数，MN/m4； hn承台埋深，m。4地基土竖向抗力系数C0、Cb和地基土竖向抗力系数的比例系数m0（1）桩底面地基土竖向抗力系数C0 C0=m0h （4-8

10、）式中：m0桩底面地基土竖向抗力系数的比例系数，kN/m4，近似取m0=m； h桩的入土深度(m)，当h小于10m时，按10m计算。（2）承台底地基土竖向抗力系数Cb Cb=m0hn （4-9）式中：hn承台埋深(m)，当hn小于1m时，按1m计算。 岩石地基竖向抗力系数C0 表 3-17单轴极限抗压强度标准值RC（MPa）C0（MN/m3）12530015000注：当RC为表列数值的中间值时，C0采用插入法确定。3、单桩计算1）采用正方形桩宽为1m的，桩长为15m。桩顶假设与地面齐平。根据规范，各种参数计算如下表。地基土的比例系数m(N/m4)基础宽度b（m）形状换算系数k空间工作系数k0=

11、1+1/b基础的计算宽度b1基础弹模E(N/m2)2000000011223E+10基础抗弯I(m4)EI(KNm2)土中基础变形系数(1/m),=(mb1)/(EI).2基础入土深度h（m）h<2.5时按照刚性基础计算，否则弹性基础0.0833333333E+060.437156.5601722）采用正方形桩宽为1.5m的，桩长为15m。桩顶假设与地面齐平。根据规范，各种参数计算如下表。地基土的比例系数m(N/m4)基础宽度b（m）形状换算系数k空间工作系数k0=1+1/b基础的计算宽度b1基础弹模E(N/m2)200000001.511.66672.53E+10基础抗弯I(m4)EI(KNm2)土中基础变形系数(1/m),=(mb1)/(EI).2基础入土深度h（m）h<2.5时按照刚性基础计算，否则弹性基础0.4218751E+070.331154.9593443）在上述参数条件下，桩基土弹簧刚度计算见表“M值法计算地基弹性刚度”。根据表中数据在civil中建立桩基模型，进行计算，所得结果与桥博比较，得出下列表格。单桩,桩顶在地面，桩顶受到100KN水平力桩顶位移（mm）桩基最大弯矩（KN/m)midas Civil