随着科学的发展，数值方法在工程上的应用越来越广泛，其中有限元方法在工程上的应用最广。但有限元在处理大变形的非线性力学问题中存在网格畸变或缠结等不足，同时有限元在求解液体振动、裂纹扩展、材料相变和成形等不定边界或可动边界问题时，需要重新划分网格，而新旧网格之间物理量转换将产生新的误差。因此，无网格方法应运而生。无网格方法不需要借助任何单元，直接利用节点构造插值函数，避免了有限元方法的以上缺陷。在多种无网格方法中自然单元法以其独特的优势获得了快速的发展。当前，超长桩大量应用于超高层建筑、大跨桥梁和深水港口工程中。然而，超长桩现阶段的研究相对滞后，工程中仍按普通桩理论进行设计，现有理论及分析模型不能充分考虑分层土特性或横纵荷载共同作用的影响，也无法很好地反映超长桩的承载性状。特别在港口工程中，超长桩通常承受很大的水平荷载，并且在横纵荷载共同作用下，往往发生大变形，而传统有限元方法在模拟计算大变形的超长桩和土体的相互作用时存在网格畸变或缠结等问题。针对以上问题本文在自然单元法理论基础上进一步改进了插值函数，并将其应用到成层地基中超长桩水平承载力问题上。主要的研究工作及创新点如下：1、改进了自然单元法的插值函数，编制三维无网格自然单元法计算程序。本文采用的是Voronoi图中二阶结构的边元素作为插值变量，大幅度提高了计算效率。无网格程序不需要借助任何单元，直接利用节点构造插值函数，避免了有限元中的网格畸变等问题，弥补了有限元在处理大变形、非线性力学等问题中存在的明显缺陷，具有较高的工程实用价值和广阔的应用前景。将程序计算结果与解析解及有限元计算结果相比较，验证了本程序的正确性和合理性。2、基于变分原理计入横纵向荷载共同作用下超长桩的P-Δ效应，推导出单元刚度修正矩阵。针对横纵向荷载共同作用下超长桩的大变形问题，建立考虑非线性大变形的无网格自然单元法计算方法。通过对实例的计算分析，验证了该方法的正确性和合理性，得到P-Δ效应特征。结果表明：由于P-Δ效应，桩的位移响应与荷载为非线性关系。当地基土质较差，桩自由长度较大时，P-Δ效应对桩身位移和内力的影响不可忽略。本方法为计算时考虑P-Δ效应提供了一种简便易行的方法，在工程设计和施工中具有一定的价值。3、基于无网格自然单元法计算并分析了荷载大小、加载顺序、长径比、桩顶约束条件、水平荷载作用位置、桩土相对刚度比和荷载分布形式等诸多参数对超长桩水平承载力的影响，并得出相关结论。4、当层状地基中采用p-y曲线法模拟桩周土非线性特性时，桩周土水平位移的计算误差随荷载增大而增大。针对该问题本文采用层状弹性体系理论考虑土体纵向连续性，并利用层状各向同性体的研究方法，建立了适用于桩周层状地基的水平位移系数传递矩阵解法。根据有限单元等效载荷的计算原理，推导出水平位移系数矩阵，并给出外荷载较大时，p-y曲线法桩周土体水平位移的修正式。根据所建立的模型编制程序，对某三层地基土水平位移系数矩阵进行了计算和分析，验证了该方法的正确性和合理性。由于计入了土体的纵向连续性，本文方法所得的水平位移影响系数曲线在荷载作用点的及其邻近区域较Mindlin解的曲线更平滑，位移小于Mindlin解的结果。这表明当土层间性质差异较大时，本文方法能更好地体现层状土体实际分布差异的影响和临近土层间的相互作用。