桩一土相互作用集中质量模型的土弹簧刚度计算方法

桩一土相互作用集中质量模型的土弹簧刚度计算方法

桩一土相互作用集中质量模型的

土弹簧刚度计算方法

孙利民刘东潘龙王君杰

（同济大学桥梁工程系）

[摘要]本文针对桥梁柱一土相互作用问题PushOver分析法中, 如何合理地确定土弹簧的刚度和土体的变形的课题进行分析计算。研究了不同地震强度下上弹簧刚度的变化特性，并将Penzien的方法和桥梁设计规范中的 "m法'计算结果进行比较，为桩一土相互作用问题的理论分析和参数选取提供重要的手段和依据。

关键词桩一土相互作用土弹簧刚度土体位移 Penzien模型 m法

一、引言

对于城市高架桥梁、大跨桥梁等桩承重要工程结构，除保证其上部结构的抗震安全性外，在遭受大地震作用时避免其基础受损也十分重要。近几年国外发生的大地震（如日本神户地震等）的震害表明，坐落在软弱土层上的桥梁桩基的震害十分突出，桩土相互作用这一课题又引起了人们的重视。

对于基础坐落在软弱土层上的桥梁结构来说，在地震发生时，桥梁上部结构的惯性力将通过基础反馈给地基，使地基产生局部变形。同时，地基自身也会因地震力作用而发生变形，反过来影响上部结构的反应。这即所谓地基一结构系统的相互作用。考虑地基一结构系统的相互作用的影响，不仅可以更准确地掌握桥梁上部结构的地震反应，对于正确计算土中基础的内力和变形也十分必要。

土与结构相互作用的研究已有近60～70年的历史，待别是近30年来，计算机技术的发展为其提供了有力的分析手段。桩基础是土建工程中广泛采用的基础形式之一，许多建于软土地基上的大型桥梁结构往往都采用桩基础，桩一土动力相互作用又是土一结构相互作用问题中较复杂的课题之一。至今已有不少关于桩基动力特性的研究报告，国内外研究人员[1-8]也提出了许多不同的桩一土动力相互作用计算方法。从研究成果的归类来看，理论上主要有离散理论和连续理论及两者的结合，解决的方法一般有集中质量法、有限元法、边界元法和波动场法。

60～70年代，美国学者J.penzien[9]等在解决泥沼地上大桥动力分析时提出了集中质量法，目前已在国内外得到了广泛的应用。集中质量法将桥梁上部结构多质点体系和桩一土体系的质量联合作为一个整体，来建立整体耦联的地震振动微分方程组进行求解。该模型假定桩侧土是Winkler连续介质。以半空间的Mindlin静力基本解为基础，将桩一土体系的质量按一定的厚度简化并集中为一系列质点，离散成一理想化的参数系统。并用弹簧和阻尼器模拟土介质的动力性质，形成一个包括地下部分的多质点体系。

PenZien方法的优点是可以方便地考虑成层土的非均匀性，非线性和阻尼特性等因素。其计算力学图式中，上下部结构均采用多质点有限元体系，便于直观理解。同时计算比较简便，经过适当的参数调整，该模型可以较好地反映桩的动力性能，因而在桩基桥梁抗震计算的实际工程中应用极广。

桥梁桩基础的抗震设计目前还主要采用静力的方法，土对桩基的作用通过一组等效的弹簧来表示。最近，日本等多地震国家的规范已开始建设使用pushOver的方法。该方法虽为一种非线性的静力分析

方法，但可以等效地反应结构与土相互作用的主要动力特性，而且计算简单，便于应用于工程设计。包括桩基在内的桥梁系统的PushOver计算除考虑上部结构惯性力的作用外，还要考虑地基土的水平变形对桩基的作用。已往往后者对桩基的抗震性能评价起决定性的作用。在建立计算图式时，合理地确定土弹簧的水平刚度和土的侧向变形是PushOver方法的关键。土弹簧刚度的确定，除考虑使用较为精确的有限元或边界元方法外，较为简便的方法是采用Penzien模型中提供的土弹簧计算方法或参照现行规范中土弹簧的计算方法。

我国<<公路桥规则>>(JTJ024-85)[10]用的m法在计算方法和参数选取方面比Penzien的方法要简单和方便，且为国内广大工程师所熟.m法的基本原理也是将桩作为弹性地基L的梁，按Winkler假定（梁身任一点的土抗力和该点的位移成正比）求解。但是，由于桩上相互作用的实验数据不足，土的物性取值有时亦缺乏合理性，在确定土弹簧的刚度时，仍有不少问题未能很好解决。特别是，m法中m的取值对弹簧刚度的计算结果影响很大，且不能反映地震波的频率特性和强度带来的影响。

本文针对桥梁桩土相互作用问题PushOver分析法中，如何合理地确定土弹簧的刚度和土体的变形的课题进行分析计算。研究了不同地震强度下土弹簧刚度的变化特性，并将Penzien的方法和桥梁设计规范中的"m法"计算结果进行比较，为桩土相互作用问题的理论分析和参数选取提供重要的手段和依据。

二、研究思路与计算理论

1研究思路

对同一桥梁基础实例，分别用两种方法来计算土层在水平方向的等效线弹簧刚。第一种方法为Penzien[9]的方法。用Milldlin公式确定水平土弹簧刚度，上述计算中土的剪切模量采用由SHAKE程序在不同输入地震波情况下算出的等效线性剪切模量。第二种方法直接采用规范[10]中的m法确定土的地基系数C,再由其算出土弹簧的水平刚度。 m的取值根据土的物性而定，并考虑规范所给的范围的中值和低值以比较其影响。

对于计算结果，通过比较方法1中土的等效剪切模量的变化讨论地震波种类和大小对土层刚度及变形的影响。并直接比较两种方法算出的土弹簧刚度值及讨论m法的m值取值合理性。

本文的另一主要工作是用SHAKE程序计算不同地震波输入时各土层的最大变形，为桩-土相互作用计算研究提供必要的数据。

2土层的等效刚度与水平位移计算程序SHAKE

SHAKE是一个对水平层状沉积土进行等效线性地震反应分析所应用的计算机程序，由PerSchnabel 博士和John Lysmer教授[12]于1970－1971年提出基本理论并编写。是迄今为止计算水平层状沉积土的地震反应时使用最广泛的程序．该程序计算的是半无限水平沉积上层的反应，该土层位于承受竖向传播的剪切波的均匀半空间之上。分析是在频域内完成，对于任何特性的确定该分析都是线性的。在考虑土壤的非线性行为时，SHAKE程序的计算原理中采用了等效线性的概念，即使用了剪切模量G和阻尼比λ与十的应交关系曲线。这两类关系曲线可从土质试验取得，G随着剪应变的增大而减小，λ随着剪应变的增大而增大。地震发生时土的剪切模量和阻尼在体系运动过程中都不是常数，从而存在土的非线性问题。对上述非线性特性作进一步的处理，认为体系在整个地震作用期的反应则以用一个平均意义下的不变的剪切模量和阻尼比来计算，采用迭代法可以得到此两项数值针对某一地震波的收敛值，即土的等