层状地基轴向受荷单桩的边界单元法奇异处理.docxVIP

层状地基轴向受荷单桩的边界单元法奇异处理 1 桩-土作用研究 桩基面是土木工程中常见的基本形式。目前，它被应用于线性理论的分析和计算。Poulos提出了以Mindlin解为基础的比较完整的弹性理论法来分析桩基础,随后Poulos和Davis又对该方法进行了系统的归纳和总结,给出了一系列图表,供设计使用。Muki和Sternberg提出了虚拟桩的概念来分析桩-土相互作用问题,并把该问题最终归结为第2类Fredholm积分方程的求解。Butterfield和Banerjee运用边界单元法来分析桩-土作用问题,通过引入桩表面的未知虚拟力分布函数(包括径向和竖向),根据桩-土变形协调条件建立桩-土共同作用的奇异积分方程,通过对奇异积分方程的求解,得出桩的轴力和位移。对于层状地基中的桩-土作用问题,国内外不少学者也进行了研究。在国外,Chin等和Lee等分别研究了2层地基中的桩基计算问题,Lee和Small运用有限层法来分析了层状地基中的桩基;在国内,金波等、艾智勇、王建华等也分别对层状地基中的桩-土相互作用问题进行了分析。 运用边界元法对桩-土相互作用问题进行研究有其独特的优点:该方法只在边界上剖分单元,并通过基本解把域内未知量化为边界未知量来求解,这就使得自由度数目大大减少,而且由于基本解本身的奇异性特点,使得用边界元法解决奇异问题时的精度较高,但合理反映真实地基性状的基本解的选取和对其奇异性的处理是运用该方法分析问题成功的关键。天然地基往往呈层状分布,故本文以层状地基内部作用一竖向集中力时的广义Mindlin解作为边界单元法的基本解,运用边界元法来分析层状地基中的轴向受荷单桩,克服了以往边界元法只能分析均质地基中的桩基础的缺陷;另外,本文采用精度控制来自动划分奇异积分区间的方法来处理基本解的奇异性。在此基础上,编制了相应的计算程序,并进行了数值分析和计算,研究结果表明,该方法具有较快的计算速度和良好的计算精度。 2 桩-土位移接触及虚拟竖向位移 如图1所示,层状地基中的一根轴向受荷圆形单桩,长度为L,半径为a,桩周围第i层土的弹性模量、泊松比、厚度分别为Esi、vsi、ΔHi。根据间接边界元法的原理,假定一个虚应力系统ψ,当这个虚应力系统作用在层状地基中时,它会产生一个相应的边界位移和应力,而且这个边界位移和应力应该与实际桩的边界位移和应力条件完全一样。 假设桩杆的竖向虚应力为ψs,并且它在层状地基深度为h处沿桩杆边界均匀分布。由桩杆上的一个微量高度dh处的竖向虚应力引起点D(r,z)的竖向位移为dw1(r,z),它可通过由点荷载(ψsadθdh)引起竖向位移的积分获得 式中:w(h,r1,z)为r=r1时w(h,r,z)的值;w(h,r,z)为竖向集中力作用下层状地基竖向位移的基本解,具体的求解过程和表达式可参考文献。 整个桩杆的竖向虚应力在D(r,z)点处引起的竖向位移为 假设在桩端底面上的点o′(λ,θ,L)处作用竖向虚应力ψb,同理可得由于桩端虚应力ψb作用而引起点D(r,z)的竖向位移: 因此,由于桩杆和桩端的虚应力引起的点D(r,z)处的竖向位移为 式(4)可以通过下列边界条件求解:在D(r,0)处,σz=τrz=0;在D(a,z)处,且0≤z≤L,W(a,z)=f1(z);在D(r,L)处,且0≤r≤a,W(r,L)=f1(L)=常数。 上述边界条件中f1(z)是一个预先假设的桩表面的虚拟竖向位移,且f1(L)为当z=L时f1(z)的值。 通过预先假定的桩-土位移边界条件,可以求出桩表面处的虚应力ψs和ψb,具体过程详见下一节。 将相应的桩-土位移边界条件代入式(4),可得 3 单桩地震差分求解 式(5)为层状地基中桩基础的精确积分方程,然而,由于基本解的复杂性,要想获得其解析解是不可能的,只能获得其数值解。把桩杆划分为n个相等的单元,桩端划分为m个圆环。对于桩杆的第i个单元,其竖向位移可以表示为 式中:i=1,2,3,…,n。 对于桩端,同理有: 式中:i=1,2,3,…,m。 式(6)、(7)中不同的Kij值可以通过简单的积分求得,其表达式如下: 未知的虚应力为 对于刚性单桩来说,在给定桩顶位移的情况下,其桩侧及桩端所有的竖向位移都是相等的,所以,通过式(13)可以求得假设的虚应力,进而可以获得该问题的解答,但对于可压缩性单桩来说,由于桩的压缩性的影响,各个单元的竖向位移是不相等的。对于一维桩杆,有: 式中:Ep为桩的弹性模量;Pz为深度为z处的桩的轴向荷载,它可以由下式获得: 由于W是有关z的微分函数,因此,可以把它展开成关于z的泰勒级数形式(忽略3阶以上的无穷小量),其形式如下: 对于单元2到n-1的表达式,可以将式(16)、(17)相减得到: 式中:当h→0时,则O(h3)将以h的3次幂的形式趋向于0,可