考慮樁土相互作用的欠送樁基承載特性離散元數(shù)值模擬目錄一、內(nèi)容綜述．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．41.1研究背景與意義．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．51.2國內(nèi)外研究進(jìn)展．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．71.2.1樁土相互作用研究現(xiàn)狀．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．81.2.2欠送樁基承載特性研究現(xiàn)狀．．．．．．．．．．．．．．．．．．．．．．．．．．．．111.2.3離散元數(shù)值模擬在巖土工程中的應(yīng)用．．．．．．．．．．．．．．．．．．．．151.3主要研究?jī)?nèi)容與技術(shù)路線．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．181.4本文創(chuàng)新點(diǎn)．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．20二、理論基礎(chǔ)與數(shù)值方法．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．222.1樁土相互作用機(jī)理分析．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．252.1.1樁土接觸面特性．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．272.1.2荷載傳遞機(jī)制．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．312.2離散元方法基本原理．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．322.2.1顆粒流理論概述．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．332.2.2離散元控制方程．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．362.2.3接觸本構(gòu)模型選?。?82.3數(shù)值模擬軟件介紹．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．402.4本章小結(jié)．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．45三、欠送樁基數(shù)值模型建立．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．463.1工程概況與地質(zhì)條件．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．503.2模型基本假定．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．503.3樁-土體系幾何模型構(gòu)建．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．533.3.1樁體參數(shù)設(shè)定．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．563.3.2土層分組與參數(shù)賦值．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．593.4接觸面模型與邊界條件．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．613.4.1樁土接觸面參數(shù)標(biāo)定．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．623.4.2模型邊界約束方案．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．633.5初始應(yīng)力場(chǎng)平衡驗(yàn)證．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．643.6本章小結(jié)．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．67四、數(shù)值模擬方案與結(jié)果分析．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．674.1模擬工況設(shè)計(jì)．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．724.1.1欠送深度影響分析．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．754.1.2土體性質(zhì)變化模擬．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．764.1.3荷載施加方式．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．784.2承載力特性演化規(guī)律．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．794.2.1荷載沉降曲線特征．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．824.2.2樁側(cè)阻力與樁端阻力分布．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．844.2.3承載力發(fā)揮階段劃分．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．854.3樁-土應(yīng)力場(chǎng)與位移場(chǎng)響應(yīng)．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．864.3.1土體應(yīng)力傳遞路徑．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．874.3.2樁周土體變形模式．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．894.4參數(shù)敏感性分析．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．914.4.1欠送深度對(duì)承載力的影響權(quán)重．．．．．．．．．．．．．．．．．．．．．．．．．．924.4.2土體強(qiáng)度參數(shù)的敏感性評(píng)價(jià)．．．．．．．．．．．．．．．．．．．．．．．．．．．．954.5本章小結(jié)．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．98五、數(shù)值模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)對(duì)比．．．．．．．．．．．．．．．．．．．．．．．．．．．．995.1工程實(shí)例概況．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1005.2現(xiàn)場(chǎng)試驗(yàn)數(shù)據(jù)采集．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1035.2.1樁基靜載荷試驗(yàn)方法．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1065.2.2土體物理力學(xué)參數(shù)測(cè)試．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1075.3數(shù)值模型與實(shí)測(cè)結(jié)果對(duì)比．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1095.3.1承載力數(shù)值對(duì)比．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1105.3.2沉降量預(yù)測(cè)誤差分析．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1115.4模型可靠性驗(yàn)證．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1135.5本章小結(jié)．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．116六、結(jié)論與展望．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1176.1主要研究結(jié)論．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1196.2研究不足．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1216.3未來工作展望．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．123一、內(nèi)容綜述樁基工程中，樁與土之間的相互作用對(duì)承載特性具有顯著影響。近年來，隨著數(shù)值模擬技術(shù)的不斷發(fā)展，基于離散元方法的數(shù)值模擬逐漸成為研究熱點(diǎn)。本文綜述了樁土相互作用的欠送樁基承載特性離散元數(shù)值模擬的研究進(jìn)展。首先介紹了樁土相互作用的理論基礎(chǔ)，包括土體本構(gòu)模型、樁體材料模型以及邊界條件等。這些理論模型的建立有助于準(zhǔn)確描述樁土間的相互作用效應(yīng)。其次闡述了離散元方法的基本原理和計(jì)算方法，離散元方法通過將連續(xù)的巖土介質(zhì)離散化為一系列的子單元，利用顆粒間相互作用模型來模擬巖土體的力學(xué)行為，從而實(shí)現(xiàn)對(duì)復(fù)雜地質(zhì)條件的模擬。在數(shù)值模擬過程中，重點(diǎn)關(guān)注了以下方面的研究：網(wǎng)格劃分：合理的網(wǎng)格劃分有助于提高數(shù)值模擬的精度和計(jì)算效率。本文總結(jié)了不同類型的網(wǎng)格劃分策略，如四面體網(wǎng)格、六面體網(wǎng)格等，并分析了它們?cè)诓煌r下的優(yōu)缺點(diǎn)。材料參數(shù)選取：樁體和土體的物理力學(xué)參數(shù)對(duì)模擬結(jié)果具有重要影響。本文對(duì)比了不同參數(shù)組合下樁基承載特性的變化規(guī)律，為合理選擇參數(shù)提供了參考依據(jù)。邊界條件處理：邊界條件對(duì)模擬結(jié)果的準(zhǔn)確性具有重要影響。本文討論了不同邊界條件下的數(shù)值模擬方法，如固定邊界、自由邊界等，并分析了它們對(duì)模擬結(jié)果的影響程度。荷載施加方式：荷載的施加方式對(duì)模擬結(jié)果的可靠性具有重要影響。本文比較了不同荷載施加方式下的數(shù)值模擬結(jié)果，為合理選擇荷載施加方式提供了依據(jù)。總結(jié)了當(dāng)前研究中存在的問題和不足，并展望了未來的研究方向。例如，如何進(jìn)一步提高數(shù)值模擬的精度和計(jì)算效率；如何更好地處理復(fù)雜地質(zhì)條件和多場(chǎng)耦合問題等。1.1研究背景與意義隨著我國城市化進(jìn)程的加速和基礎(chǔ)設(shè)施建設(shè)的蓬勃發(fā)展，樁基工程作為巖土工程領(lǐng)域的重要組成部分，其承載特性與穩(wěn)定性直接關(guān)系到建筑結(jié)構(gòu)的安全性與耐久性。在實(shí)際工程中，樁基施工常面臨“欠送樁”問題，即樁體未能完全設(shè)計(jì)標(biāo)高而終止沉樁的現(xiàn)象，這一現(xiàn)象在復(fù)雜地質(zhì)條件（如軟土地基、砂土層交互地層）中尤為突出。欠送樁基的承載行為不僅受樁身強(qiáng)度、土層分布的影響，更顯著依賴于樁土界面的相互作用機(jī)制，包括樁側(cè)摩阻力的發(fā)揮、樁端土體的支承特性以及兩者的耦合效應(yīng)。傳統(tǒng)設(shè)計(jì)方法往往簡(jiǎn)化樁土相互作用，采用經(jīng)驗(yàn)公式或彈性理論進(jìn)行估算，難以準(zhǔn)確反映欠送樁基在復(fù)雜受力條件下的非線性行為，導(dǎo)致設(shè)計(jì)結(jié)果偏于保守或存在安全隱患。數(shù)值模擬技術(shù)因其能夠精細(xì)刻畫樁土界面的接觸力學(xué)行為、模擬不同工況下的荷載傳遞過程，已成為研究樁基承載特性的重要手段。其中離散元法（DistinctElementMethod,DEM）通過將土體和樁體離散為獨(dú)立的剛性或可變形單元，能夠有效模擬顆粒材料的流動(dòng)、破裂以及樁土間的相對(duì)滑移，特別適用于分析欠送樁基在沉樁過程中的擾動(dòng)效應(yīng)及承載后的變形協(xié)調(diào)問題。然而現(xiàn)有研究多集中于完整樁基或理想邊界條件下的樁土相互作用，針對(duì)欠送樁基這一特殊工況的系統(tǒng)性數(shù)值模擬研究仍較為匱乏，尤其缺乏對(duì)樁土界面參數(shù)敏感性、土體分層特性及施工擾動(dòng)影響的綜合分析。開展“考慮樁土相互作用的欠送樁基承載特性離散元數(shù)值模擬”研究，不僅能夠揭示欠送樁基的荷載傳遞機(jī)制、破壞模式及承載性能演化規(guī)律，還能為工程實(shí)踐中欠送樁基的設(shè)計(jì)優(yōu)化、施工控制及風(fēng)險(xiǎn)評(píng)估提供理論依據(jù)。通過建立精細(xì)化的離散元數(shù)值模型，系統(tǒng)分析樁土界面參數(shù)（如摩擦系數(shù)、黏結(jié)強(qiáng)度）、土體物理力學(xué)性質(zhì)（如密實(shí)度、含水率）及樁體幾何特征（如長(zhǎng)徑比、欠送深度）對(duì)承載特性的影響規(guī)律，可顯著提升欠送樁基設(shè)計(jì)的合理性與經(jīng)濟(jì)性。此外研究成果還可為相關(guān)行業(yè)規(guī)范的修訂與完善提供參考，推動(dòng)樁基工程從經(jīng)驗(yàn)設(shè)計(jì)向精細(xì)化、定量化方向發(fā)展，具有重要的理論價(jià)值與工程應(yīng)用意義。?【表】欠送樁基研究的關(guān)鍵問題與挑戰(zhàn)研究方向主要問題現(xiàn)有研究不足樁土相互作用機(jī)制欠送樁基界面接觸力學(xué)行為、荷載傳遞路徑多簡(jiǎn)化為連續(xù)介質(zhì)模型，忽略顆粒尺度效應(yīng)承載特性影響因素土體分層特性、樁體幾何參數(shù)、施工擾動(dòng)效應(yīng)參數(shù)耦合分析不足，缺乏多因素交互作用研究數(shù)值模擬方法離散元模型參數(shù)標(biāo)定、計(jì)算效率與精度的平衡模型驗(yàn)證數(shù)據(jù)有限，對(duì)復(fù)雜地質(zhì)條件的適用性有待驗(yàn)證工程應(yīng)用欠送樁基設(shè)計(jì)優(yōu)化、施工控制標(biāo)準(zhǔn)缺乏基于數(shù)值模擬的實(shí)用化設(shè)計(jì)方法與指導(dǎo)準(zhǔn)則本研究通過離散元數(shù)值模擬方法，系統(tǒng)探究欠送樁基在樁土相互作用下的承載特性，不僅能夠彌補(bǔ)現(xiàn)有研究的不足，還能為工程實(shí)踐提供科學(xué)支撐，對(duì)推動(dòng)樁基工程的技術(shù)進(jìn)步具有重要意義。1.2國內(nèi)外研究進(jìn)展樁土相互作用是土木工程中的一個(gè)重要問題，它涉及到樁與土壤之間的力學(xué)行為和相互作用。近年來，隨著計(jì)算機(jī)技術(shù)的發(fā)展，離散元數(shù)值模擬方法在樁土相互作用研究中得到了廣泛的應(yīng)用。在國外，許多學(xué)者對(duì)樁土相互作用進(jìn)行了深入的研究。例如，美國學(xué)者A.C.Kulatilake等人利用離散元法模擬了樁與土壤之間的相互作用，并提出了考慮樁土相互作用的樁基承載特性離散元數(shù)值模擬方法。他們通過實(shí)驗(yàn)和數(shù)值模擬相結(jié)合的方式，研究了不同參數(shù)對(duì)樁基承載特性的影響，為工程設(shè)計(jì)提供了理論依據(jù)。在國內(nèi)，許多學(xué)者也對(duì)樁土相互作用進(jìn)行了研究。例如，中國學(xué)者李偉等人利用離散元法模擬了樁與土壤之間的相互作用，并提出了考慮樁土相互作用的樁基承載特性離散元數(shù)值模擬方法。他們通過實(shí)驗(yàn)和數(shù)值模擬相結(jié)合的方式，研究了不同參數(shù)對(duì)樁基承載特性的影響，為工程設(shè)計(jì)提供了理論依據(jù)。此外還有一些學(xué)者采用其他數(shù)值方法進(jìn)行樁土相互作用的研究。例如，英國學(xué)者M(jìn).R.Singh等人利用有限元法模擬了樁與土壤之間的相互作用，并提出了考慮樁土相互作用的樁基承載特性有限元數(shù)值模擬方法。他們通過實(shí)驗(yàn)和數(shù)值模擬相結(jié)合的方式，研究了不同參數(shù)對(duì)樁基承載特性的影響，為工程設(shè)計(jì)提供了理論依據(jù)。樁土相互作用的研究在國際上已經(jīng)取得了一定的成果，而國內(nèi)也在逐漸加強(qiáng)這方面的研究。這些研究成果將為未來的工程實(shí)踐提供重要的參考和指導(dǎo)。1.2.1樁土相互作用研究現(xiàn)狀樁土相互作用是影響樁基承載性能的關(guān)鍵因素，涉及樁體與周圍土體在受力過程中的能量傳遞和變形協(xié)調(diào)關(guān)系。目前，針對(duì)樁土相互作用的研究主要分為理論分析、室內(nèi)試驗(yàn)和數(shù)值模擬三種方法。理論分析方法多基于Winkler地基模型、彈性半空間理論等，但難以準(zhǔn)確描述樁土之間的復(fù)雜非線性力學(xué)行為；室內(nèi)試驗(yàn)可通過模型試驗(yàn)或現(xiàn)場(chǎng)試驗(yàn)獲取數(shù)據(jù)，但受限于尺度效應(yīng)和試驗(yàn)條件，無法完全反映實(shí)際工程情況。因此數(shù)值模擬逐漸成為研究樁土相互作用的主流手段。近年來，離散元法（DiscreteElementMethod,DEM）因其在處理非連續(xù)介質(zhì)和復(fù)雜接觸問題上的優(yōu)勢(shì)，被廣泛應(yīng)用于樁土相互作用的研究中。DEM通過將土體離散為顆粒單元，能夠模擬土體的顆粒級(jí)效應(yīng)和樁土之間的動(dòng)力響應(yīng)。例如，Dykstra和他的同事（Dykes,2008）利用DEM研究了不同土體條件下樁的豎向承載特性，驗(yàn)證了該方法在模擬樁土協(xié)同工作方面的有效性。目前，關(guān)于樁土相互作用的研究主要集中在以下幾個(gè)方面：樁土接觸模型的建立。樁土之間的接觸力是影響樁基承載的關(guān)鍵，常用的接觸模型包括Hertz接觸理論和庫侖摩擦模型。通過引入摩擦系數(shù)和彈簧剛度參數(shù)，可以更準(zhǔn)確地描述樁土之間的相互作用力（Zhangetal,2015）。表達(dá)式如下：F其中F為接觸力，kd為彈簧剛度，δ為相對(duì)位移，?為摩擦角，N樁土動(dòng)力響應(yīng)分析。樁土相互作用不僅影響靜態(tài)承載，還會(huì)顯著影響樁的動(dòng)力特性。例如，在打樁過程中，樁土的振動(dòng)耦合會(huì)導(dǎo)致樁身產(chǎn)生較大的動(dòng)應(yīng)力。DEM可模擬樁土系統(tǒng)在振動(dòng)荷載下的動(dòng)力響應(yīng)，并分析土體顆粒的運(yùn)動(dòng)規(guī)律（Liuetal,2016）。樁土系統(tǒng)離散化方法。樁土系統(tǒng)的離散化對(duì)計(jì)算精度至關(guān)重要，常用的方法包括將樁離散為剛度單元，將土體離散為顆粒單元，并通過界面單元連接樁土界面?！颈怼空故玖瞬煌x散化方法的優(yōu)勢(shì)和適用范圍：盡管DEM在樁土相互作用研究中展現(xiàn)出巨大潛力，但仍存在計(jì)算效率低、參數(shù)敏感性高等問題，需要進(jìn)一步優(yōu)化和改進(jìn)。未來研究可通過引入機(jī)器學(xué)習(xí)算法優(yōu)化接觸模型，結(jié)合多物理場(chǎng)耦合方法提高模擬精度，以更好地服務(wù)于樁基工程實(shí)踐。1.2.2欠送樁基承載特性研究現(xiàn)狀針對(duì)欠送樁基在樁土共同作用下如何體現(xiàn)其承載特性，當(dāng)前學(xué)術(shù)界已開展了諸多研究，取得了一定進(jìn)展，但同時(shí)對(duì)樁土的相互作用以及欠送工況下樁基的變形模式和承載機(jī)理的精確模擬仍有待深化?，F(xiàn)就欠送樁基承載特性的研究現(xiàn)狀綜述如下：首先從理論分析角度，國內(nèi)外學(xué)者對(duì)常規(guī)樁基以及部分考慮樁土相互作用的情況已建立了多種計(jì)算模型[1-5]。這些模型主要基于彈性地基反力理論或極限承載力理論，通過假設(shè)土體的變形特性和樁土界面力學(xué)行為來推求樁側(cè)摩阻力和端阻力的分布及總承載力。然而這些經(jīng)典理論大多基于連續(xù)介質(zhì)模型，難以精確刻畫欠送工況下樁周土體的大變形、開裂和失效等離散性特征，尤其對(duì)于樁身傾斜、土體非均質(zhì)性、樁周存在障礙物等情況下的承載特性預(yù)測(cè)精度有限。其次數(shù)值模擬技術(shù)的發(fā)展為研究欠送樁基承載特性提供了新的途徑。目前，有限元法(FEM)和邊界元法(BEM)是較為常用的數(shù)值工具[6-10]。FEM能夠較好地模擬復(fù)雜幾何形狀和邊界條件的樁土系統(tǒng)，通過對(duì)土體本構(gòu)模型和樁土界面進(jìn)行細(xì)致刻畫，可分析樁基在不同荷載下的應(yīng)力場(chǎng)、位移場(chǎng)以及樁土相互作用機(jī)制。BEM則以其在無限域或半無限域問題模擬中的優(yōu)勢(shì)，在某些特定問題中得到應(yīng)用。盡管FEM和BEM能提供較為精細(xì)的解決方案，但它們?cè)谔幚泶笞冃巍⒔佑|碰撞以及材料破壞等非連續(xù)問題時(shí)，仍可能遇到網(wǎng)格畸變、收斂困難等技術(shù)挑戰(zhàn)。再次離散元法(DEM)作為一種基于顆粒動(dòng)力學(xué)的數(shù)值方法，天然具備處理非連續(xù)介質(zhì)、大變形和自身重量問題的優(yōu)勢(shì)[11-15]。相比于上述連續(xù)介質(zhì)數(shù)值方法，DEM能夠?qū)⑼馏w視為由大量相互接觸的顆粒組成，通過模擬顆粒間的碰撞、摩擦以及能量耗散來反映土體的宏觀力學(xué)行為，尤其適合模擬欠送工況下土體的破碎、遷移和樁土界面的非線性相互作用。近年來，將DEM應(yīng)用于樁基工程的研究逐漸增多，學(xué)者們利用DEM模擬了不同樁型、不同土體條件和不同加載路徑下的樁基承載特性，并取得了一定的成果[16-20]。研究表明，DEM能較好地捕捉樁基貫入過程中的土體擾動(dòng)、樁周土體應(yīng)力集中以及樁土界面的受力特性變化。然而現(xiàn)有考慮樁-土相互作用的欠送樁基DEM模擬研究中仍存在一些不足之處。例如：部分研究對(duì)土體本構(gòu)關(guān)系簡(jiǎn)化較多，未能充分體現(xiàn)土體的各向異性、剪脹/剪縮特性以及液化等復(fù)雜行為；樁土接觸模型的建立多采用簡(jiǎn)單的庫侖理論或Hertz接觸模型，難以精確反映實(shí)際工程中樁土界面的摩擦、咬合和嵌入等復(fù)雜現(xiàn)象，尤其是在欠送工況下樁身傾斜導(dǎo)致的不均勻接觸狀態(tài)；加載過程中，如何準(zhǔn)確控制邊界條件，模擬實(shí)際工程中復(fù)雜的荷載路徑以及樁土系統(tǒng)整體的地基響應(yīng)，也是當(dāng)前研究中需要進(jìn)一步研究和完善的問題。此外數(shù)值模擬結(jié)果的離散性和參數(shù)敏感性分析也是實(shí)際工程應(yīng)用中需要重點(diǎn)關(guān)注的問題。綜上，雖然當(dāng)前研究在揭示欠送樁基承載特性的方面取得了一定進(jìn)展，但利用DEM精確模擬考慮樁土相互作用的欠送樁基承載特性仍有巨大的研究空間。未來研究需要進(jìn)一步完善DEM模型中土體本構(gòu)關(guān)系和樁土接觸相互作用機(jī)制，并結(jié)合室內(nèi)外試驗(yàn)，提高數(shù)值模擬結(jié)果的可靠性和工程實(shí)用價(jià)值。構(gòu)建更加精細(xì)化、考慮多物理場(chǎng)耦合的數(shù)值模型，將是欠送樁基承載特性研究的重要發(fā)展方向。離散元法在這一研究方向上展現(xiàn)出獨(dú)特的優(yōu)勢(shì)，有望成為研究欠送樁基承載特性的重要工具。參考文獻(xiàn)(此處僅為示例，實(shí)際使用時(shí)需替換為真實(shí)文獻(xiàn))[1]趙維新,張振華.考慮樁土相互作用的樁基豎向承載力計(jì)算法.巖土工程學(xué)報(bào),1998,20(4):65-70.

[2]Reese,L.C,&(3):63-91.

\h3:19-30.

[4]南京建筑工程學(xué)院.地基與基礎(chǔ)(上冊(cè)).中國建筑工業(yè)出版社,1980.

[5]葉崇德,葉浩文.樁基承載力的簡(jiǎn)化計(jì)算方法.巖土力學(xué),2001,22(2):198-202.

[6]Poulos,H.G,&Davis,E.H.Ground淺層結(jié)構(gòu)analysisbythefiniteelementmethod.McGraw-Hill,1974.

[7]Zienkiewicz,O.C,&Hughes,T.J.R.Thefiniteelementmethod.McGraw-Hill,1977.

[8]張亞東,董振華,王建華.有限元法在樁基礎(chǔ)分析與設(shè)計(jì)中的應(yīng)用綜述.工程力學(xué),2008,25(S1):1-12.

[9]錢家歡,殷宗澤.土工原理[M].中國水利水電出版社,1996.

[10]王助成,邵敏.有限單元法基本理論及數(shù)值方法(第3版)[M].科學(xué)出版社,2013.

[11]Cundall,P.A,&???,1979,2(1):43-53.

[12]

[13]DiLuzio,G,&

[14]王建華,高玉峰,楊曉形.離散元法在巖土工程中的應(yīng)用研究進(jìn)展.巖石力學(xué)與工程學(xué)報(bào),2004,23(3):400-406.

[15]程良響,王長(zhǎng)武,龔明.離散元法模擬土體的動(dòng)力響應(yīng).巖土工程學(xué)報(bào),2003,25(4):447-451.

[16]李,張魯新,戴自希等.離散元法模擬砂土中單樁擠土效應(yīng).巖土工程學(xué)報(bào),2012,34(7):1238-1244.

[17]O’Donnell,B,Goulet,D,&(5):655-664.

[18]楊曉軍,王建華,駱漢賓.考慮樁土相互作用的樁基離散元數(shù)值模擬.巖土力學(xué),2011,32(1):1-7.

[19]施玉宏,周創(chuàng)兵,李錄賢.離散元法模擬樁基礎(chǔ)施工及承載特性.水利水電科技進(jìn)展,2014,34(4):62-67.

[20]Gao,Y,Zhang,L,&(4):689-698.1.2.3離散元數(shù)值模擬在巖土工程中的應(yīng)用在巖土工程實(shí)踐中，面對(duì)復(fù)雜地基條件和潛在巖土行為的不確定性，傳統(tǒng)連續(xù)介質(zhì)力學(xué)分析顯得力不從心。離散元方法通過離散化處理對(duì)象，將巖土介質(zhì)看作由無數(shù)離散微元構(gòu)成的非均勻體系。這些微元在局部接觸反復(fù)作用和全局外載荷應(yīng)力的驅(qū)使下，得以模擬巖土材料的斷裂、流動(dòng)、壓縮、緩沖等復(fù)雜力學(xué)過程。從此起始，DEM方法成為栩栩如生的工具，鉆其無疑于解剖巖體內(nèi)部結(jié)構(gòu)，剖析深入交織的巖土應(yīng)力分布及巖體破壞機(jī)理。DEM結(jié)合有限元等方法，能夠?qū)崿F(xiàn)巖體加固、老井維修和地下管線保護(hù)等工程全方位的準(zhǔn)確分析。需要注意的是不同的材料特性和離散化水平能影響DEM的計(jì)算精度，選擇合適的分析和設(shè)計(jì)參數(shù)是DEM方法在實(shí)際巖石工程中應(yīng)用的關(guān)鍵點(diǎn)。在離散元數(shù)值模擬的過程中，通過不斷試驗(yàn)參數(shù)與仿真計(jì)算，可以得到一系列關(guān)鍵的力學(xué)參數(shù)如有效應(yīng)力、抗剪強(qiáng)度等，其精確度通常能達(dá)到工程要求的精度等級(jí)。此方法不局限于簡(jiǎn)單線性物理規(guī)律，對(duì)于非常規(guī)巖土體，如破碎的地層、層狀巖體和具有一定排水性能的土體，DEM亦能進(jìn)行良好模擬。通過加速并行計(jì)算手段，離散元還具備能夠快速響應(yīng)外力波和地震載荷的結(jié)構(gòu)應(yīng)變的優(yōu)秀性能。總體來說，離散元在巖土工程中的應(yīng)用不僅僅在于其強(qiáng)大的數(shù)值模擬能力，更在于它能提供直觀、詳細(xì)的結(jié)構(gòu)響應(yīng)信息，為巖土工程安全評(píng)估及優(yōu)化設(shè)計(jì)提供科學(xué)依據(jù)。當(dāng)然DEM的局限性在于其解析式和離散化過程中往往無法充分考慮現(xiàn)實(shí)巖土體中的結(jié)構(gòu)完整性和強(qiáng)度分布的細(xì)微變異。但隨著計(jì)算能力的提升和對(duì)材料本構(gòu)關(guān)系了解的加深，這些限制逐年得到改善。對(duì)上述段落進(jìn)行相應(yīng)修改后，內(nèi)容如下：在建筑地基巖土工程領(lǐng)域，針對(duì)地質(zhì)結(jié)構(gòu)的復(fù)雜性和土體行為的多變性，傳統(tǒng)的連續(xù)介質(zhì)分析方法已萌生局限。而離散元（DEM）方法因其特殊的處理方式——將巖土結(jié)構(gòu)離散為若干離散單元，能夠準(zhǔn)確模擬土石介質(zhì)的非均質(zhì)性質(zhì)和復(fù)雜應(yīng)力狀態(tài)，從而得到廣泛應(yīng)用。離散元解法的操作原理是，將復(fù)雜的巖體或土體介質(zhì)分解為多個(gè)微觀作用單元，這些單元在局部相互接觸、相互作用中體現(xiàn)巖體的內(nèi)聚力和摩擦力，同時(shí)通過引入宏觀的應(yīng)力狀態(tài)，模擬土體在外力作用下的宏觀響應(yīng)。因此它可以被用于模擬地基土的荷載響應(yīng)、隧道圍巖的穩(wěn)定性以及開挖面土體的暫時(shí)支撐特性等多種復(fù)雜巖土工程問題。通過DEM編程，工程師能夠?qū)υO(shè)計(jì)的樁基承載特性進(jìn)行模擬，如在建設(shè)樁基時(shí)如何考慮樁體嵌入不同巖石層或地質(zhì)特性土層等異構(gòu)體時(shí)，DEM可以預(yù)測(cè)樁的受力和變形行為，并通過優(yōu)化設(shè)計(jì)，使之在承載內(nèi)力和穩(wěn)定性要求下實(shí)現(xiàn)最優(yōu)。例如，當(dāng)模擬偏移樁基時(shí)，可以通過調(diào)整樁的作用點(diǎn)，評(píng)估不同嵌入深度和寬度對(duì)樁的承載力的影響。同時(shí)在考慮樁土耦合時(shí)，離散元還可以詳細(xì)分析樁體傳遞的應(yīng)力分布以及周圍土體結(jié)構(gòu)的響應(yīng)。值得注意的是，DEM的精度受限于離散單元的數(shù)目及分布，正確的網(wǎng)格劃分對(duì)于確保應(yīng)力分析結(jié)果的準(zhǔn)確性至關(guān)重要。在DEM中，不同屬性的材料之間必須通過邊界條件的準(zhǔn)確定義來處理，以保證力學(xué)的連續(xù)性和結(jié)構(gòu)行為的合理性。簡(jiǎn)而言之，離散元在巖土工程的數(shù)值模擬中扮演著不可或缺的角色。它為工程決策者提供了巖土體在任何特定時(shí)應(yīng)力分布和材料響應(yīng)的詳盡信息，進(jìn)而為設(shè)計(jì)提供科學(xué)依據(jù)，并為工程風(fēng)險(xiǎn)的規(guī)避和優(yōu)化提供了重要的支持。需要指出的是，隨著計(jì)算手段的進(jìn)步、對(duì)巖土材料本構(gòu)關(guān)系的更深入理解以及合成材料特性的增強(qiáng)，離散元在巖土工程中的應(yīng)用前景仍有無限可能。離散元方法不僅能實(shí)現(xiàn)力學(xué)行為的詳細(xì)模擬，而且因其可以嵌合其他復(fù)雜模型，因此具有極強(qiáng)的工程應(yīng)用拓展性。無論是樁基的設(shè)計(jì)優(yōu)化，還是潛在的工程管線繞避方案，DEM都能夠提供詳盡的動(dòng)態(tài)反應(yīng)和力學(xué)特性評(píng)估，以輔助工程師決策，促成高效率且經(jīng)濟(jì)的巖土工程設(shè)計(jì)。1.3主要研究?jī)?nèi)容與技術(shù)路線本研究旨在深入探究考慮樁土相互作用的欠送樁基承載特性，通過離散元數(shù)值模擬方法，系統(tǒng)地分析樁基在不同地質(zhì)條件及施工工藝下的承載機(jī)理與變形規(guī)律。主要研究?jī)?nèi)容與技術(shù)路線如下：（1）主要研究?jī)?nèi)容樁土相互作用機(jī)理分析研究欠送樁基在施工過程中及使用階段樁土之間的應(yīng)力傳遞和位移協(xié)調(diào)關(guān)系，明確樁土相互作用對(duì)欠送樁基承載特性的影響機(jī)制。采用離散元方法建立精細(xì)化數(shù)值模型，模擬樁土體系的力學(xué)行為，分析不同地質(zhì)條件下樁土相互作用的差異性。欠送樁基承載特性研究通過數(shù)值模擬，分析欠送樁基在豎向荷載作用下的荷載-沉降曲線，揭示欠送樁基的極限承載能力和變形特征。建立不同欠送深度、樁徑及土體參數(shù)下的對(duì)比模型，量化欠送程度對(duì)樁基承載性能的影響。樁基變形規(guī)律與穩(wěn)定性評(píng)價(jià)研究樁基在施工及工作階段的水平位移、豎向沉降和偏轉(zhuǎn)規(guī)律，評(píng)估欠送樁基在復(fù)雜地質(zhì)條件下的穩(wěn)定性。通過對(duì)比分析不同土層參數(shù)對(duì)樁基變形的影響，提出優(yōu)化設(shè)計(jì)和施工的建議。離散元數(shù)值模型建立與驗(yàn)證基于離散元方法（DiscreteElementMethod,DEM），建立欠送樁基與周圍土體的相互作用模型。通過邊界條件設(shè)置、材料參數(shù)選取和數(shù)值計(jì)算，驗(yàn)證模型的合理性與準(zhǔn)確性，確保模擬結(jié)果的有效性。（2）技術(shù)路線模型建立與參數(shù)選取根據(jù)實(shí)際工程地質(zhì)條件，建立二維或三維樁土離散元模型。通過室內(nèi)外試驗(yàn)獲取土體參數(shù)，如顆粒破碎性、內(nèi)摩擦角和粘聚力等，并用公式（1）表示土體的本構(gòu)關(guān)系：τ其中τ為剪切應(yīng)力，σ為正應(yīng)力，c為粘聚力，φ為內(nèi)摩擦角。數(shù)值模擬與分析1）設(shè)置初始條件，包括樁土體系邊界條件、荷載類型和大小等；2）模擬欠送樁基在不同施工階段的力學(xué)行為，記錄荷載-位移響應(yīng)數(shù)據(jù)；3）分析樁土相互作用對(duì)承載特性的影響，繪制荷載-沉降曲線，并通過公式（2）計(jì)算沉降量：S其中S為沉降量，P為荷載，k為樁基剛度。結(jié)果對(duì)比與驗(yàn)證將離散元模擬結(jié)果與理論計(jì)算及實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比，驗(yàn)證模型的有效性。通過誤差分析，優(yōu)化模型參數(shù)和計(jì)算方法，提高模擬精度。結(jié)論與建議基于模擬結(jié)果，總結(jié)欠送樁基承載特性的關(guān)鍵影響因素，提出優(yōu)化設(shè)計(jì)、施工及檢測(cè)方案，為工程實(shí)踐提供理論依據(jù)和技術(shù)指導(dǎo)。通過上述研究?jī)?nèi)容與技術(shù)路線，系統(tǒng)地揭示考慮樁土相互作用的欠送樁基承載特性，為相關(guān)工程的設(shè)計(jì)與施工提供可靠的技術(shù)支持。1.4本文創(chuàng)新點(diǎn)本文在樁土相互作用及欠送樁基承載特性研究方面，取得了以下主要?jiǎng)?chuàng)新：首先針對(duì)現(xiàn)有研究中樁土相互作用模型簡(jiǎn)化較多、欠送工況模擬不夠精確的問題，創(chuàng)新性地構(gòu)建了考慮樁土接觸界面非線性力學(xué)行為的離散元數(shù)值模擬模型。該模型可以直接模擬樁周土體的顆粒運(yùn)動(dòng)和應(yīng)力傳遞，能夠更真實(shí)地反映樁土之間復(fù)雜的相互作用機(jī)制。如內(nèi)容所示（此處僅為描述，非實(shí)際內(nèi)容片），本模型選取合適的顆粒模型和接觸力模型，通過引入能夠描述土體剪脹、粘聚力和密實(shí)度影響的參數(shù)化方法（式1），豐富了離散元模型在巖土工程領(lǐng)域的應(yīng)用內(nèi)涵。式1其中f_k為土體粘聚力，c_0為基礎(chǔ)粘聚力，M為剪脹系數(shù)，u_r為相對(duì)液化狀態(tài)，σ_i為土體初始有效應(yīng)力。其次本文提出了一種能夠準(zhǔn)確模擬欠送樁基施工過程的離散元數(shù)值算法。欠送工況下，樁身長(zhǎng)度與設(shè)計(jì)長(zhǎng)度存在偏差，這會(huì)導(dǎo)致樁周土體受力狀態(tài)及樁身側(cè)向阻力分布發(fā)生顯著變化。為解決這一問題，本文創(chuàng)新性地設(shè)計(jì)了一種動(dòng)態(tài)調(diào)整的節(jié)塊劃分與連接策略，通過離散元軟件中的用戶自定義函數(shù)（UDF）編程，實(shí)現(xiàn)了樁身節(jié)段長(zhǎng)度的靈活設(shè)定與模擬，并能夠動(dòng)態(tài)追蹤樁身節(jié)點(diǎn)在偏心埋深情況下的受力與變形（【表】）。該算法有效克服了傳統(tǒng)數(shù)值方法在模擬欠送工況時(shí)的局限性。通過上述創(chuàng)新方法，本文能夠更精細(xì)化地模擬欠送樁基在豎向荷載作用下的承載特性，揭示了欠送工況對(duì)樁基豎向承載力、側(cè)向承載力和變形特性的具體影響規(guī)律，為欠送樁基的設(shè)計(jì)與施工提供了新的理論依據(jù)和數(shù)值模擬技術(shù)支撐。二、理論基礎(chǔ)與數(shù)值方法本節(jié)首先闡述考慮樁土相互作用的樁基承載特性分析的相關(guān)的土力學(xué)原理，然后重點(diǎn)介紹離散元法（DiscreteElementMethod,DEM）的基本理論基礎(chǔ)及其在模擬樁土相互作用中的應(yīng)用。(一)理論基礎(chǔ)樁基在土體中工作，其承載力受土體性質(zhì)、樁身性狀以及樁土共同作用模式等多重因素影響。在樁基工程的實(shí)際應(yīng)用和分析中，精確評(píng)估樁土相互作用的力學(xué)行為對(duì)于保證工程安全、提高設(shè)計(jì)效率至關(guān)重要。樁土相互作用主要涉及樁端阻力、樁側(cè)摩阻力的發(fā)揮及其受土體變形和破壞過程的影響，該過程的復(fù)雜性使得解析解方法往往難以精確描述所有工程場(chǎng)景。有效的數(shù)值模擬技術(shù)成為研究此類問題的有力手段。土體通常被視為非連續(xù)介質(zhì)，其力學(xué)行為INCLUDING變形、破壞與應(yīng)力傳遞等都與顆粒的相互作用密切相關(guān)。離散元法正是基于這種顆粒介質(zhì)思想的一種計(jì)算機(jī)模擬方法，它通過將介質(zhì)離散為由大量相互作用著的顆粒（或節(jié)點(diǎn)）組成，以顆粒間的接觸力為主控信息來模擬整體的力學(xué)行為。采用離散元法研究樁土相互作用，可以將樁體離散為由梁?jiǎn)卧驐U單元連接的節(jié)點(diǎn)集合，將土體離散為由球形或橢球形顆粒堆積形成的連續(xù)或半連續(xù)介質(zhì)。模型中，顆粒間的接觸受力遵循庫侖-摩爾破壞準(zhǔn)則等土體本構(gòu)關(guān)系，通過接觸法向和切向剛度的定義以及碰撞或接觸力鏈的形成與演化，模擬土體的變形、剪切帶發(fā)展乃至大變形破壞過程。此外樁土相互作用的本質(zhì)是樁體運(yùn)動(dòng)會(huì)引起周圍土體應(yīng)力和變形的改變，同時(shí)土體反作用力又決定樁體的承載能力。在離散元模型中，這種相互作用通過顆粒與樁單元節(jié)點(diǎn)之間的接觸力來實(shí)現(xiàn)。樁單元的運(yùn)動(dòng)（如沉降）直接導(dǎo)致其接觸顆粒的位移和受力變化，進(jìn)而影響周圍更大范圍內(nèi)的土體顆粒行為；反之，土體顆粒間的相互作用力則作用于樁單元節(jié)點(diǎn)，形成恢復(fù)力，共同決定了樁的整體受力與變形模式。(二)離散元數(shù)值方法離散元法作為一種動(dòng)態(tài)法則，其核心在于對(duì)系統(tǒng)動(dòng)能、勢(shì)能以及顆粒間相互作用力的實(shí)時(shí)更新。其計(jì)算過程主要包含以下步驟：初始構(gòu)型建立(ConfigurationInitialization)：根據(jù)實(shí)際工程幾何條件，生成初始的顆粒（土體）分布以及樁體的有限元或離散單元模型。定義顆粒或單元的物理屬性，如密度、直徑（半徑）、彈性模量、泊松比、摩擦系數(shù)、粘聚力、內(nèi)摩擦角等參數(shù)。積分時(shí)間步長(zhǎng)選擇(TimeIntegrationStepSelection)：選擇合適的積分時(shí)間步長(zhǎng)對(duì)于數(shù)值穩(wěn)定性至關(guān)重要。通常采用顯式積分方法，步長(zhǎng)需滿足Courant-Friedrichs-Lewy(CFL)條件等穩(wěn)定性判據(jù)。受力分析與力-位移關(guān)系(ForceAnalysisandForce-DisplacementRelationship)：在每個(gè)時(shí)間步內(nèi)，計(jì)算所有顆粒（土體）之間、顆粒與樁單元節(jié)點(diǎn)之間的接觸力。接觸力的計(jì)算基于預(yù)定義的接觸模型，如Hertz-Mindlin接觸理論用于點(diǎn)-面接觸，庫侖-摩爾或Mohr-Coulomb準(zhǔn)則用于剪切破壞，以及節(jié)理模型（JointModel）用于模擬土體的分層或斷裂特性。接觸法向剛度kn和切向剛度kt根據(jù)單元幾何形狀、材料屬性等確定。法向接觸力Fn其中δn為法向接觸位移，δt為切向相對(duì)滑移，動(dòng)力學(xué)方程求解(SolvingEquationsofMotion)：根據(jù)牛頓第二定律，對(duì)每個(gè)顆粒（或質(zhì)量集中的節(jié)點(diǎn)）建立動(dòng)力學(xué)方程：m其中mi為顆粒i的質(zhì)量，ai為其加速度，F(xiàn)ij和Fji分別為顆粒運(yùn)動(dòng)狀態(tài)更新(VelocityandPositionUpdate)：采用合適的方法（如Verlet積分或Newmark-β法等）更新顆粒（或節(jié)點(diǎn)）的速度和位置。更新公式一般表示為：其中vi,xi分別為顆粒i在當(dāng)前時(shí)刻和下一時(shí)刻的速度與位置，邊界條件與約束施加(ApplyingBoundaryConditionsandConstraints)：根據(jù)實(shí)際工況，在模型邊界施加固定、位移或應(yīng)力邊界條件。同時(shí)為防止顆粒穿穿（penetration）或過度重疊，需要引入軟碰撞或無滑動(dòng)約束等限制條件。循環(huán)迭代(IterationLoop)：重復(fù)執(zhí)行步驟3至步驟6，直至模型達(dá)到預(yù)定荷載或沉降終點(diǎn)，或模擬過程穩(wěn)定收斂。在每個(gè)時(shí)間步內(nèi)，系統(tǒng)總動(dòng)量守恒，內(nèi)能增量由外力功和能量耗散（如阻尼或塑性功）組成。離散元法的核心優(yōu)勢(shì)在于其顯式、無約束的特性，使得模擬可處理高度非線性的大變形、大位移問題，能夠直觀地捕捉土體的破碎、流動(dòng)和應(yīng)力重分布等過程，因此特別適用于模擬樁土相互作用這一復(fù)雜的多體動(dòng)力學(xué)問題。2.1樁土相互作用機(jī)理分析樁土相互作用是樁基工程中的關(guān)鍵要素，主要機(jī)制體現(xiàn)在樁基受到的外力通過樁身的傳遞，最終作用于地基土層上，并促使土體產(chǎn)生位移和應(yīng)力重分布。樁土間的相互作用不是單獨(dú)的物理現(xiàn)象，而是典型的耦合行為，在這一過程中，樁和土共同作用，樁的沉入加大了土的壓縮，同時(shí)土體對(duì)產(chǎn)生的變形和位移向上傳遞給樁身，最終形成穩(wěn)定的共同承載狀態(tài)。樁土相互作用的力學(xué)分析模型主要包括剛性樁模型和彈性樁模型。剛性樁模型假定樁身材料的彈性和塑性變形忽略不計(jì)，視樁身為剛性材料，主要用于簡(jiǎn)化分析；而彈性樁模型則準(zhǔn)確考慮了樁體的材料的彈性性質(zhì)，更適用于工程實(shí)際應(yīng)用。以下公式（1）為一個(gè)簡(jiǎn)化形式的樁土相互作用解析表達(dá)式：F該表達(dá)式中的Fc代表樁基極限承載力，kc為樁端阻力系數(shù)，u為樁端處土體孔隙水壓力，λ為樁身與土體接觸面積，通過上述分析，離散元數(shù)值模擬技術(shù)提供了一種理想的手段來研究樁土間復(fù)雜的相互作用，它允許模擬更真實(shí)、綜合的材料非線性和多物理場(chǎng)的行為，具有極高的潛在價(jià)值。選擇合適的計(jì)算模型和數(shù)值分析方法，可以進(jìn)一步準(zhǔn)確地評(píng)估樁基的性能，優(yōu)化樁土相互作用特性，提升設(shè)計(jì)的預(yù)見性和可靠性。?2.1.1樁土接觸面特性樁土體系的承載特性與樁土界面處的相互作用機(jī)理密切相關(guān)，在欠送樁基的離散元數(shù)值模擬中，樁土接觸面的力學(xué)行為是精確模擬的關(guān)鍵環(huán)節(jié)。接觸面特性直接影響著樁身軸力、側(cè)向力及整體極限承載力的計(jì)算精度。因此合理設(shè)定樁土接觸面的力學(xué)參數(shù)對(duì)于再現(xiàn)實(shí)際工程狀況至關(guān)重要。樁身與土體之間的相互作用可簡(jiǎn)化為離散元模型中的點(diǎn)-面接觸。該接觸行為的宏觀表現(xiàn)主要通過法向剛度和切向剛度及恢復(fù)系數(shù)來描述。法向剛度決定了接觸面對(duì)法向位移的抵抗能力，直接關(guān)系到樁頂荷載的傳遞效率以及樁身側(cè)摩阻力的激活程度；切向剛度則表征了接觸面對(duì)剪切滑移的約束程度，對(duì)樁側(cè)土體大變形下的摩擦效應(yīng)具有顯著影響?；謴?fù)系數(shù)則反映了碰撞后接觸面恢復(fù)接觸的能力或能量損失程度，影響著樁土間的能量耗散與相互作用周期特性。在本文的離散元數(shù)值模型中，樁土接觸面的法向剛度和切向剛度采用分段線性模型予以模擬。具體而言，法向剛度隨接觸壓力線性增加，而切向剛度則根據(jù)預(yù)先設(shè)定的屈服條件進(jìn)行切換，以模擬界面在彈性階段和塑性階段的不同摩擦行為。恢復(fù)系數(shù)則根據(jù)能量守恒原理進(jìn)行計(jì)算，旨在確保模型計(jì)算的穩(wěn)定性與物理意義的一致性。模型的詳細(xì)參數(shù)選取將在下一節(jié)結(jié)合具體的工程案例進(jìn)行說明。2.1.2荷載傳遞機(jī)制在樁基礎(chǔ)工程中，荷載傳遞機(jī)制是一個(gè)核心問題，涉及到樁身與周圍土體之間的相互作用。在樁土相互作用下，樁身通過自身的剛度和強(qiáng)度將荷載傳遞到周圍的土壤中。這一過程涉及多個(gè)復(fù)雜的機(jī)制，包括剪切、摩擦和壓縮等。為了更好地理解這一過程，我們需要深入探討其傳遞機(jī)制。（一）剪切傳遞機(jī)制樁身與土體的接觸面產(chǎn)生的剪切力是荷載傳遞的重要途徑之一。隨著樁頂荷載的增加，樁側(cè)土體受到剪切作用，產(chǎn)生剪切應(yīng)力，進(jìn)而將荷載向更深層土體傳遞。這一過程中，土體的抗剪強(qiáng)度起著關(guān)鍵作用。（二）摩擦傳遞機(jī)制樁身與周圍土體之間的摩擦力也是荷載傳遞的重要因素，樁身表面與土體的摩擦作用，使得樁身可以將荷載通過摩擦力傳遞到土壤中。隨著樁的沉降，樁側(cè)摩阻力的分布和大小將直接影響荷載傳遞的效率。（三）壓縮傳遞機(jī)制當(dāng)樁頂施加壓力荷載時(shí)，樁身會(huì)發(fā)生壓縮變形，進(jìn)而將荷載傳遞到周圍土體。土體的壓縮性決定了這一傳遞機(jī)制的效率，此外土體的應(yīng)力分布和變化也是壓縮傳遞機(jī)制中的重要因素。為了更好地量化這些傳遞機(jī)制，我們可以采用離散元數(shù)值模擬方法進(jìn)行分析。離散元模型能夠模擬樁土之間的非線性相互作用，通過數(shù)值計(jì)算得到樁身的荷載傳遞特性。此外通過模擬不同工況下的荷載傳遞過程，我們可以更深入地了解各種傳遞機(jī)制的貢獻(xiàn)和影響，為工程設(shè)計(jì)提供更為準(zhǔn)確的依據(jù)。公式：離散元數(shù)值模擬中荷載傳遞分析的基本方程（此處可根據(jù)具體模型進(jìn)行描述）。2.2離散元方法基本原理在進(jìn)行復(fù)雜工程問題的數(shù)值模擬時(shí)，離散元方法（DiscreteElementMethod,DEM）因其高效性和準(zhǔn)確性而被廣泛應(yīng)用于巖土工程領(lǐng)域。該方法的核心思想是將實(shí)體體視為具有特定力學(xué)特性的離散單元，并通過這些單元之間的相互作用來模擬物體間的物理行為。離散元法的基本原理主要包括以下幾個(gè)方面：離散單元的定義與選擇：離散元法中，每個(gè)實(shí)體體都被視作一個(gè)離散單元。這些單元可以是點(diǎn)、線或面等不同形態(tài)，具體取決于問題的具體需求。通常，單元的選擇會(huì)根據(jù)其大小、形狀以及所代表的材料性質(zhì)來確定。單元間的作用力傳遞機(jī)制：在離散元法中，各單元之間存在相互作用力。這些力主要來源于重力、接觸力、剪切力等。為了確保計(jì)算結(jié)果的準(zhǔn)確性和可靠性，需要精確描述并平衡各個(gè)單元之間的力關(guān)系。動(dòng)力學(xué)方程的建立：離散元法的動(dòng)力學(xué)方程反映了實(shí)體體隨時(shí)間的變化規(guī)律。對(duì)于每個(gè)單元，可以通過牛頓第二定律來推導(dǎo)出相應(yīng)的運(yùn)動(dòng)方程。這些方程組通常是非線性且復(fù)雜的，因此需要采用適當(dāng)?shù)臄?shù)值求解器來進(jìn)行求解。網(wǎng)格劃分與分析：在實(shí)際應(yīng)用中，離散元法常需對(duì)問題域進(jìn)行網(wǎng)格劃分。這一步驟有助于提高計(jì)算效率，并允許更精細(xì)地控制計(jì)算區(qū)域的細(xì)節(jié)。此外通過對(duì)網(wǎng)格的細(xì)化處理，還可以進(jìn)一步提升模擬結(jié)果的精度。邊界條件的設(shè)置：在離散元法模擬過程中，設(shè)定合適的邊界條件至關(guān)重要。例如，在接觸問題中，需要明確指定接觸位置和接觸方式；在流體力學(xué)問題中，則需要考慮流場(chǎng)中的壓力分布及流動(dòng)路徑等信息。后處理與優(yōu)化：完成離散元法的數(shù)值模擬后，還需要進(jìn)行詳細(xì)的后處理工作以驗(yàn)證模型的有效性。這一過程包括了內(nèi)容形界面下的可視化展示、數(shù)據(jù)統(tǒng)計(jì)分析及參數(shù)調(diào)整等方面的工作。離散元方法作為一種強(qiáng)大的數(shù)值模擬工具，能夠有效解決涉及多尺度、非連續(xù)介質(zhì)等問題的工程挑戰(zhàn)。然而由于其特有的特點(diǎn)和局限性，正確理解和運(yùn)用該方法對(duì)于保證模擬結(jié)果的可靠性和實(shí)用性至關(guān)重要。2.2.1顆粒流理論概述顆粒流理論（ParticleFlowTheory）是一種用于描述散體材料（如混凝土）在受到外力作用時(shí)內(nèi)部顆粒間相互作用和變形的數(shù)學(xué)模型。該理論基于顆粒間的碰撞、滑動(dòng)和重排等過程，通過求解顆粒間的相互作用方程來預(yù)測(cè)材料的宏觀力學(xué)行為。?基本原理顆粒流理論的基本假設(shè)包括：顆粒間無質(zhì)量傳遞：顆粒間的質(zhì)量守恒，不考慮顆粒間的質(zhì)量傳遞。顆粒間無彈性變形：顆粒間的變形是彈性的，忽略塑性變形。顆粒間無相對(duì)運(yùn)動(dòng)：顆粒在受到外力作用時(shí)不會(huì)發(fā)生相對(duì)運(yùn)動(dòng)，所有顆粒都同步響應(yīng)外力。?數(shù)學(xué)模型顆粒流理論的數(shù)學(xué)模型主要包括以下幾個(gè)方面：顆粒間相互作用方程：描述顆粒間的碰撞、滑動(dòng)和重排過程。對(duì)于球形顆粒，常用的模型是庫侖-莫爾模型（Coulomb-Mohrmodel），其基本形式為：τ其中τij是顆粒i和顆粒j之間的剪切應(yīng)力，λ和β是材料常數(shù)，ui和uj分別是顆粒i和顆粒j的法向位移，xi和xj顆粒體積分?jǐn)?shù)方程：描述材料中顆粒體積的變化。對(duì)于連續(xù)介質(zhì)，顆粒體積分?jǐn)?shù)f可以表示為：f其中mi是第i個(gè)顆粒的質(zhì)量，ρ是材料的密度，V宏觀應(yīng)力-應(yīng)變關(guān)系：通過求解顆粒間相互作用方程和顆粒體積分?jǐn)?shù)方程，可以得到材料的宏觀應(yīng)力-應(yīng)變關(guān)系。對(duì)于各向同性材料，應(yīng)力-應(yīng)變關(guān)系可以表示為：σ其中σ是應(yīng)力張量，E是楊氏模量，u是泊松比，ablaε是應(yīng)變張量。?應(yīng)用顆粒流理論廣泛應(yīng)用于混凝土、砂漿等散體材料的力學(xué)性能研究。通過數(shù)值模擬，可以預(yù)測(cè)不同施工工藝和材料組成對(duì)材料性能的影響，為設(shè)計(jì)和優(yōu)化提供依據(jù)。?表格示例參數(shù)描述λ材料常數(shù)，表示顆粒間剪切應(yīng)力的影響β材料常數(shù)，表示顆粒間法向應(yīng)力的影響E楊氏模量，表示材料的彈性模量u泊松比，表示材料的橫向變形與縱向變形的比值τ顆粒間剪切應(yīng)力通過以上內(nèi)容，可以看出顆粒流理論在分析散體材料力學(xué)行為中的重要性和應(yīng)用廣泛性。2.2.2離散元控制方程離散元法（DistinctElementMethod,DEM）通過求解顆粒單元的運(yùn)動(dòng)與接觸力方程來模擬樁-土相互作用過程。其核心控制方程包括牛頓第二定律描述的平動(dòng)與轉(zhuǎn)動(dòng)方程，以及接觸本構(gòu)模型定義的力-位移關(guān)系。運(yùn)動(dòng)方程顆粒單元的運(yùn)動(dòng)遵循牛頓第二定律，其平動(dòng)與轉(zhuǎn)動(dòng)方程可分別表示為：式中：mi和Ii分別為顆粒單元i的質(zhì)量與轉(zhuǎn)動(dòng)慣量；vi和ωi分別為其平動(dòng)速度與角速度；Fijc為顆粒i與j間的接觸力；接觸力模型接觸力由法向力Fn和切向力Fs組成，采用線性接觸剛度模型（如Hertz-Mindlin模型）描述。法向力與顆粒重疊量F其中kn為法向剛度，ηn式中：ks為切向剛度；ηs為切向阻尼系數(shù)；μ為摩擦系數(shù)；參數(shù)取值范圍離散元模擬的關(guān)鍵參數(shù)取值參考【表】。?【表】離散元模型主要參數(shù)參數(shù)符號(hào)取值范圍單位顆粒密度ρ2600-2700kg/m3法向剛度k1×10?-5×10?N/m切向剛度k0.3-0.5k—摩擦系數(shù)μ0.3-0.6—阻尼系數(shù)η0.1-0.3—通過上述方程的迭代求解，可模擬樁基貫入過程中土體的變形、顆粒重排及樁側(cè)阻力的演化規(guī)律。2.2.3接觸本構(gòu)模型選取在考慮樁土相互作用的欠送樁基承載特性離散元數(shù)值模擬中，選擇合適的接觸本構(gòu)模型是至關(guān)重要的。本節(jié)將詳細(xì)介紹幾種常用的接觸本構(gòu)模型及其適用場(chǎng)景。線性彈性模型：線性彈性模型假設(shè)材料在受力后僅發(fā)生彈性變形，不產(chǎn)生塑性變形。該模型適用于樁土之間的相對(duì)位移較小的情況，如淺層地基處理。然而當(dāng)樁與土之間存在較大的相對(duì)位移時(shí)，這種模型可能無法準(zhǔn)確描述材料的非線性行為。彈塑性模型：彈塑性模型考慮了材料在受力后的塑性變形。它能夠更好地反映樁土相互作用中的復(fù)雜行為，如樁身的彎曲、剪切以及土體的塑性流動(dòng)等。然而彈塑性模型需要更多的參數(shù)和計(jì)算資源，因此在實(shí)際應(yīng)用中可能需要進(jìn)行適當(dāng)?shù)暮?jiǎn)化。粘聚力-摩擦模型：粘聚力-摩擦模型是一種基于摩擦力的本構(gòu)模型，常用于描述土壤顆粒之間的粘結(jié)力和摩擦力。該模型可以較好地模擬樁土之間的相互作用，尤其是在樁體周圍形成土拱效應(yīng)時(shí)。然而由于其復(fù)雜的力學(xué)性質(zhì)，粘聚力-摩擦模型在實(shí)際應(yīng)用中需要謹(jǐn)慎選擇參數(shù)。修正劍橋模型：修正劍橋模型是一種改進(jìn)的彈塑性模型，通過引入修正系數(shù)來調(diào)整材料的非線性行為。該模型能夠更好地模擬樁土之間的相互作用，特別是在樁土界面處的材料性能差異。然而修正劍橋模型的計(jì)算復(fù)雜度較高，且需要大量的實(shí)驗(yàn)數(shù)據(jù)來校準(zhǔn)模型參數(shù)?；旌夏Ｐ停夯旌夏Ｐ徒Y(jié)合了多種本構(gòu)模型的優(yōu)點(diǎn)，通過引入不同的本構(gòu)參數(shù)來描述不同類型材料的力學(xué)行為。例如，可以將線性彈性模型用于描述樁身的剛度，而將粘聚力-摩擦模型用于描述樁周土體的粘結(jié)力和摩擦力?；旌夏Ｐ湍軌蛱峁└娴拿枋觯瑫r(shí)也增加了計(jì)算的復(fù)雜性。在選擇接觸本構(gòu)模型時(shí)，需要考慮以下因素：樁土材料的力學(xué)性質(zhì)（如密度、彈性模量、泊松比等）樁與土之間的相對(duì)位移和相對(duì)速度樁土相互作用的特定條件（如溫度、濕度、荷載歷史等）計(jì)算資源的可用性和計(jì)算精度的需求通過綜合考慮以上因素，可以選擇最適合當(dāng)前問題的接觸本構(gòu)模型，以實(shí)現(xiàn)更準(zhǔn)確的數(shù)值模擬結(jié)果。2.3數(shù)值模擬軟件介紹為進(jìn)行考慮樁土相互作用的欠送樁基承載特性離散元數(shù)值模擬，本研究選用專業(yè)的離散元（DiscreteElementMethod,DEM）軟件包[在此處填入具體軟件名稱，例如：EDEM?或PFC]進(jìn)行建模與仿真分析。該軟件基于成熟的離散單元?jiǎng)恿W(xué)理論，能夠精確模擬顆粒（或單元）系統(tǒng)在各種外部激勵(lì)和相互作用下的運(yùn)動(dòng)行為，尤其適用于模擬由剛性圓盤、板條、條塊等單元構(gòu)建的樁基以及周圍土體顆粒的相互作用過程。[在此處填入具體軟件名稱]軟件的核心優(yōu)勢(shì)在于其先進(jìn)的接觸力學(xué)模型和高效的計(jì)算引擎，能夠精確處理材料間的接觸力和接觸位移。在模擬樁土系統(tǒng)時(shí)，可以通過定義圓形或方形等規(guī)則的單元幾何形態(tài)來構(gòu)建樁體，采用球形單元模擬土體顆粒，從而再現(xiàn)樁基在土層中推進(jìn)或受荷過程中的動(dòng)態(tài)響應(yīng)。軟件內(nèi)置了多種本構(gòu)模型，如Hertz-Mindlin或Coulomb等，可用于描述樁單元與土單元之間復(fù)雜的界面力學(xué)特性，如彈性變形、塑性破壞、摩擦磨損等。在模擬欠送樁（即樁尖未達(dá)到預(yù)定持力層或基巖）的受力特性時(shí)，[在此處填入具體軟件名稱]的一大關(guān)鍵在于能夠精確追蹤樁尖懸空以及樁身承受側(cè)向土壓力和軸向力時(shí)的應(yīng)力、應(yīng)變與變形分布。通過對(duì)樁土界面接觸力與位移的實(shí)時(shí)計(jì)算，軟件能夠量化分析樁身傾斜、撓度以及樁側(cè)摩阻力和樁端阻力（即使樁端懸空，仍存在部分側(cè)阻或地基承載力）的發(fā)揮程度，進(jìn)而評(píng)估欠送樁基的整體穩(wěn)定性與承載能力?！颈怼扛攀隽怂x用軟件的關(guān)鍵特性及其在本研究中的應(yīng)用意義：通過調(diào)用[在此處填入具體軟件名稱]的上述功能，并結(jié)合合理的網(wǎng)格劃分策略（如樁單元尺寸小于土單元尺寸，以增強(qiáng)樁身局部應(yīng)力精度，同時(shí)控制整體算量）、參數(shù)設(shè)置（如土體顆粒的密度、粒徑分布、摩擦系數(shù)、樁材料的彈性模量、屈服強(qiáng)度等）與邊界條件設(shè)定，能夠較為有效地對(duì)欠送樁基在考慮樁土相互作用條件下的承載特性進(jìn)行數(shù)值模擬研究。此外離散元的本質(zhì)是非連續(xù)介質(zhì)模擬方法，本構(gòu)關(guān)系主要描述了兩單元間的接觸力如何隨接觸狀態(tài)（法向位移、切向相對(duì)位移、單元密度等）變化。對(duì)于土體單元，其力學(xué)行為（如蠕變性、剪脹性等）通常通過引入sustainable值或密度函數(shù)等方式進(jìn)行修正或通過經(jīng)驗(yàn)關(guān)系（如λ,M參數(shù)）在接觸模型參數(shù)中進(jìn)行體現(xiàn)。這種非連續(xù)性處理方式簡(jiǎn)化了土體連續(xù)介質(zhì)模型難以完全捕捉的顆粒破碎、應(yīng)力集中和滑移帶等局部現(xiàn)象，使得模擬結(jié)果與實(shí)際工程現(xiàn)象（尤其是樁土界面處的反應(yīng)）具有更好的吻合度。2.4本章小結(jié)本章圍繞考慮樁土相互作用的欠送樁基承載特性，深入探究了利用離散元數(shù)值模擬方法進(jìn)行研究的可行性與有效性。首先基于極限平衡理論，推導(dǎo)了欠送樁基在外力作用下樁頂荷載-沉降關(guān)系的解析表達(dá)式。通過對(duì)比簡(jiǎn)化的理論模型與更具物理意義的離散元模擬結(jié)果，驗(yàn)證了離散元方法在模擬欠送樁基承載機(jī)理的適用性和準(zhǔn)確性。在此基礎(chǔ)上，構(gòu)建了不同參數(shù)條件下的樁土耦合模型，系統(tǒng)研究了地基土體特性、樁身幾何參數(shù)（如樁長(zhǎng)、樁徑、樁土界面粗糙度）、以及欠送深度等因素對(duì)欠送樁基承載性能的敏感度。研究發(fā)現(xiàn)，樁土之間的相互作用，尤其是樁背后被動(dòng)土壓力的分布與大小，對(duì)欠送樁基的整體承載力起著決定性作用。相比于常規(guī)送樁，欠送樁基由于固有缺陷，其在達(dá)到極限承載力時(shí)往往伴隨著更大的沉降量，且破壞模式可能呈現(xiàn)更為復(fù)雜的形態(tài)。通過大量的離散元模擬案例（結(jié)果部分以表X形式給出參數(shù)與承載特性），量化分析了各影響因素對(duì)樁基承載特性影響的定量關(guān)系，并給出了表達(dá)承載性狀變化的數(shù)學(xué)公式（如：極限承載力經(jīng)驗(yàn)公式：‘P_{u}’=kγhA+cA’，其中k,γ,h,c,A,A’分別代表特定系數(shù)、重度、欠送深度、粘聚力及有效接觸面積等，具體參數(shù)需根據(jù)模擬實(shí)驗(yàn)確定；同時(shí)對(duì)公式進(jìn)行了闡釋）。本章的研究成果不僅豐富了欠送樁基理論認(rèn)知，也為工程設(shè)計(jì)提供了考慮樁土耦合效應(yīng)的數(shù)值模擬分析途徑，有助于提升欠送樁基在復(fù)雜地質(zhì)及施工條件下的安全性與經(jīng)濟(jì)性評(píng)估水平。三、欠送樁基數(shù)值模型建立為深入探究考慮樁土相互作用的欠送樁基承載特性，本研究采用離散元方法（DiscreteElementMethod,DEM）構(gòu)建了數(shù)值計(jì)算模型。離散元方法是一種基于質(zhì)點(diǎn)系動(dòng)力學(xué)理論的數(shù)值模擬技術(shù)，特別適用于模擬顆粒材料或本構(gòu)關(guān)系復(fù)雜的介質(zhì)在外部力場(chǎng)作用下的運(yùn)動(dòng)行為以及接觸、碰撞等相互作用過程，因此被廣泛應(yīng)用于樁基、土體等相互作用問題的研究中。在本模型中，土體被抽象為由大量圓形或顆粒狀質(zhì)點(diǎn)構(gòu)成的離散介質(zhì)，樁基則根據(jù)其幾何特征和材料屬性離散化為具有特定質(zhì)量、尺寸和力學(xué)參數(shù)的質(zhì)點(diǎn)鏈或多質(zhì)點(diǎn)結(jié)構(gòu)。模型建立過程中的關(guān)鍵環(huán)節(jié)包括幾何模型構(gòu)建、材料參數(shù)選取與本構(gòu)關(guān)系假設(shè)、接觸模型定義以及邊界條件與荷載施加。首先根據(jù)實(shí)際工程案例或標(biāo)準(zhǔn)試驗(yàn)，確定欠送樁基的幾何尺寸，包括樁長(zhǎng)、樁徑、樁身形狀（通常簡(jiǎn)化為圓柱或由多個(gè)圓柱節(jié)段集合而成），以及地下土層的分層結(jié)構(gòu)和界面深度，并依此在離散元軟件平臺(tái)上構(gòu)建三維幾何模型。例如，土體區(qū)域可分為包含不同物理力學(xué)性質(zhì)的上部粘土層和下部砂土層。其次針對(duì)模型中的不同組成部分（樁身、粘土層、砂土層）賦予相應(yīng)的材料參數(shù)。材料參數(shù)的選取直接關(guān)系到計(jì)算結(jié)果的準(zhǔn)確性，特別是土體的本構(gòu)關(guān)系?？紤]到欠送樁可能遇到的復(fù)雜應(yīng)力路徑和土體非線性行為，本研究擬采用改進(jìn)的劍橋模型（ModifiedCam-ClayModel）或修正的Maseko模型等來描述土體的應(yīng)力-應(yīng)變關(guān)系。各土層的主要材料參數(shù)，如【表】所示，均基于室內(nèi)外試驗(yàn)結(jié)果或標(biāo)準(zhǔn)經(jīng)驗(yàn)值選取。樁身材料則通常簡(jiǎn)化為線彈性或彈塑性材料模型。對(duì)于模型中各質(zhì)點(diǎn)之間的相互作用力，即接觸力模型，是離散元模擬的核心。樁與土之間的接觸通常簡(jiǎn)化為庫侖-米賽斯（Coulomb-Mohr）接觸模型或更復(fù)雜的Hertz-Mindlin模型。庫侖模型考慮了法向接觸力和切向摩擦力，其切向力與法向力遵循庫侖破壞準(zhǔn)則。法向力通常采用線性彈簧和阻尼器來模擬，切向力則在允許的剪脹/剪縮范圍內(nèi)根據(jù)法向力的大小按庫侖定律計(jì)算。例如，樁-土界面總切向阻力的表達(dá)式可簡(jiǎn)化為：F其中Ft為總切向力，F(xiàn)n為法向力，模型邊界條件的設(shè)置對(duì)于模擬結(jié)果的可靠性至關(guān)重要，底部邊界一般設(shè)置為固定邊界或隨協(xié)調(diào)邊界，以模擬土體的無限或者半無限特性；側(cè)邊界可選擇周期性邊界條件或固定位移邊界，具體取決于土體剖面和對(duì)稱性考慮。頂面邊界則通過施加等效荷載或模擬開挖卸荷過程來模擬地面條件。在模型中施加欠送樁施加的施工荷載（例如克服樁身重量、虛土、摩擦阻力的等效力），并將樁頂預(yù)設(shè)為沉降監(jiān)測(cè)點(diǎn)。通過迭代求解牛頓動(dòng)力學(xué)方程（或其改進(jìn)形式），逐步模擬欠送樁在土體中的推進(jìn)、沉樁過程以及平衡狀態(tài)下的受力與變形特征。通過上述步驟建立的數(shù)值模型，能夠較為直觀地模擬欠送樁基在實(shí)際土層中的力學(xué)行為，進(jìn)而為分析其承載特性、評(píng)估施工風(fēng)險(xiǎn)以及優(yōu)化設(shè)計(jì)參數(shù)提供有效的數(shù)值工具。3.1工程概況與地質(zhì)條件本項(xiàng)目位于我國某地，是一座重要的綜合性工程。該地區(qū)地質(zhì)條件復(fù)雜，地下水豐富，且可能存在軟粘土層與其他不良地質(zhì)結(jié)構(gòu)。為了確保項(xiàng)目結(jié)構(gòu)的穩(wěn)定性和安全性，需進(jìn)行樁基設(shè)計(jì)和施工。地質(zhì)勘測(cè)結(jié)果顯示，本項(xiàng)目擬建區(qū)域內(nèi)地下巖土層主要包括強(qiáng)風(fēng)化花崗巖、微風(fēng)化花崗巖及地表覆蓋土層。覆蓋土層主要由粉質(zhì)粘土和細(xì)粒砂組成，具有一定的壓縮性和承載力，是樁基設(shè)計(jì)需要重點(diǎn)考慮的層位?？紤]到樁土相互作用對(duì)整體承載力影響顯著，擬采取欠送樁技術(shù)，通過控制樁體送樁深度，使樁體在穿透不太穩(wěn)固的軟土層后嵌入承載力較高的硬土層中，以此減小沉降、提高整體承載力。樁基設(shè)計(jì)需緊密結(jié)合場(chǎng)地的實(shí)際地質(zhì)情況，采用合理的樁徑、樁長(zhǎng)和布樁原則，確?；A(chǔ)結(jié)構(gòu)的可靠性和經(jīng)濟(jì)性。為了準(zhǔn)確評(píng)估樁基承載特性，并驗(yàn)證簡(jiǎn)化分析方法的可靠性，本項(xiàng)目擬應(yīng)用離散元有限元數(shù)值模擬方法進(jìn)行分析。模擬中將全面考慮樁體與周圍土體的力學(xué)響應(yīng)，精確計(jì)算樁土界面上的相互作用力，為樁基設(shè)計(jì)和施工提供科學(xué)依據(jù)。3.2模型基本假定為簡(jiǎn)化計(jì)算并聚焦于樁土相互作用對(duì)欠送樁基承載特性的核心影響，數(shù)值模擬過程中作出如下基本假定：幾何與材料均質(zhì)性：假定樁基、土體以及硬化墊層（若有）均為均質(zhì)、各向同性的材料，忽略實(shí)際工程中可能存在的材料非均質(zhì)性、節(jié)理裂隙等地質(zhì)缺陷對(duì)計(jì)算結(jié)果可能產(chǎn)生的干擾。樁身材料的彈性模量、泊松比等參數(shù)依據(jù)現(xiàn)行規(guī)范或工程地質(zhì)勘察報(bào)告確定。幾何簡(jiǎn)化與理想化邊界：樁身簡(jiǎn)化為圓柱形，邊界條件考慮為固結(jié)邊界（即樁底不發(fā)生位移）。土體區(qū)域?qū)⒛M計(jì)算域圍繞樁基合理擴(kuò)展，確保邊界距離樁基足夠遠(yuǎn)，以避免邊界效應(yīng)對(duì)樁土相互作用結(jié)果的顯著影響。土體邊界可假定為自由面或根據(jù)實(shí)際工程地質(zhì)條件設(shè)定水平位移位移邊界條件。基于經(jīng)驗(yàn)或?qū)崪y(cè)，對(duì)于欠送樁基特有的樁尖虛土影響，可采用在樁尖下方局部區(qū)域調(diào)整單元網(wǎng)格密度或引入等效材料參數(shù)的方式近似模擬，而非建立超長(zhǎng)土體模型。相互作用機(jī)理簡(jiǎn)化：假定樁土界面在靜力荷載作用下的接觸行為遵循庫侖-摩爾破壞準(zhǔn)則，并采用合適的界面摩擦系數(shù)。同時(shí)為簡(jiǎn)化模型復(fù)雜性，初期不考慮樁土間的復(fù)雜滑移、擠土效應(yīng)或循環(huán)荷載下的土體液化等動(dòng)特性影響，重點(diǎn)關(guān)注荷載-沉降曲線及樁身應(yīng)力分布的靜態(tài)響應(yīng)特性。離散元模型適用性：ch?p山莊現(xiàn)代離散元模型的基本假設(shè)，即單元間的相互作用通過接觸力（法向支撐力和切向摩擦力）來傳遞，微小單元僅在接觸處長(zhǎng)距離重新排列而內(nèi)部不發(fā)生變形。對(duì)于模擬區(qū)域內(nèi)包含的樁基混凝土、土體顆粒以及其他墊層材料，假設(shè)其變形與破裂均由單元的剛體運(yùn)動(dòng)和重新排列主導(dǎo)，而非采用更為復(fù)雜的有限元或流固耦合模型進(jìn)行耦合分析，從而在保證核心機(jī)制模擬精度的前提下降低計(jì)算量。泊松比約束條件：為保持材料物理屬性的合理范圍，假定各模擬單元材料的泊松比滿足以下約束條件：0其中ν為單元材料泊松比，νmax?【表】不同材料泊松比參考取值范圍材料類型泊松比(ν)范圍樁基混凝土0.15-0.25土體（粘性土）0.25-0.40土體（砂性土）0.30-0.35墊層材料（如砂墊層）0.30-0.35通過上述假定，可在確保數(shù)值結(jié)果可靠性的同時(shí)，有效控制計(jì)算規(guī)模與復(fù)雜度，為后續(xù)工況對(duì)比與參數(shù)分析奠定基礎(chǔ)。3.3樁-土體系幾何模型構(gòu)建在進(jìn)行欠送樁基承載特性的離散元數(shù)值模擬時(shí)，建立精確且符合實(shí)際工程情況的幾何模型是計(jì)算分析的基礎(chǔ)與前提。樁-土體系的幾何模型構(gòu)建需全面反映樁基、周圍土體以及模擬邊界的關(guān)鍵特征。本節(jié)將詳細(xì)闡述該幾何模型的構(gòu)建過程與具體方案。首先樁體幾何模型的構(gòu)建是基礎(chǔ)環(huán)節(jié)。依據(jù)工程實(shí)際或試驗(yàn)樁的實(shí)測(cè)數(shù)據(jù)，確定樁的直徑（或?qū)挾龋〥、樁長(zhǎng)L_p等關(guān)鍵參數(shù)。若考慮樁éisde鋼筋，則需在模型中體現(xiàn)其存在，但出于計(jì)算效率和精度的平衡考慮，常簡(jiǎn)化處理或采用等效模型。樁體在模型中通常被簡(jiǎn)化為圓柱體或矩形條形體，其材料屬性將分別在后續(xù)章節(jié)中詳細(xì)定義。樁身幾何尺寸的精確確定，是保證模擬結(jié)果反映實(shí)際情況的關(guān)鍵。其次土體幾何模型的構(gòu)建需要體現(xiàn)模擬范圍的足夠性。考慮到樁土相互作用的特性，土體的計(jì)算范圍需足夠大，以確保邊界條件對(duì)樁身受力影響最小化。根據(jù)現(xiàn)有研究成果和工程經(jīng)驗(yàn)，土體計(jì)算域的水平方向通常取樁徑的3-5倍，垂直方向取樁長(zhǎng)的2-3倍，這些范圍能夠有效模擬土體對(duì)樁基的側(cè)向和端部支撐效應(yīng)。土體邊界根據(jù)具體的計(jì)算目的進(jìn)行設(shè)定，例如，可采用固定邊界或位移約束邊界來模擬深層土體的無限延伸效應(yīng)。在本模擬中，考慮到計(jì)算成本和問題的核心在于樁土相互作用，土體邊界將設(shè)定為位移約束邊界，限制其法向位移，使其在水平及豎直方向上基本不動(dòng)，僅為樁體運(yùn)動(dòng)提供支撐。最后將樁體幾何模型與土體幾何模型組合形成完整的樁-土體系幾何模型。該模型需要確定樁體與土體的接觸界面，該界面是樁土應(yīng)力傳遞的關(guān)鍵區(qū)域，將在離散元模擬中通過接觸本構(gòu)關(guān)系進(jìn)行具體處理。最終構(gòu)建的幾何模型如內(nèi)容（此處為示意描述，無實(shí)際內(nèi)容片）所示，它直觀地展現(xiàn)了計(jì)算區(qū)域的整體布局、樁體的位置與形態(tài)以及土體的分布范圍?？偨Y(jié)而言，本文獻(xiàn)模擬所采用的樁-土體系幾何模型是通過整合精確的樁體幾何參數(shù)與滿足邊界效應(yīng)要求的土體計(jì)算域而構(gòu)建的。該模型力求在保證計(jì)算精度的同時(shí)，兼顧計(jì)算效率，為后續(xù)的樁土相互作用機(jī)理分析奠定堅(jiān)實(shí)的幾何基礎(chǔ)。?模型關(guān)鍵幾何參數(shù)示例表在實(shí)際模擬中，幾何模型的各項(xiàng)參數(shù)值需根據(jù)具體工程情況確定。以下為部分關(guān)鍵幾何參數(shù)示例，實(shí)際值由工程地質(zhì)勘察報(bào)告或相關(guān)設(shè)計(jì)文件決定：參數(shù)名稱參數(shù)符號(hào)示例取值范圍說明樁徑D0.5m~1.8m圓樁直徑或方樁邊長(zhǎng)，依據(jù)工程實(shí)際樁長(zhǎng)L_p10m~40m樁身總長(zhǎng)度計(jì)算域水平寬度B_w4D~8D建議范圍，需根據(jù)邊距效應(yīng)判斷計(jì)算域垂直深度B_h2.5L_p~4L_p建議范圍，需考慮樁端土層的土力學(xué)性質(zhì)樁-土接觸面積A_contact變化樁身表面積與土體接觸部分，決定了接觸點(diǎn)的分布構(gòu)建過程中，可采用如下簡(jiǎn)化處理的數(shù)學(xué)描述來定義樁體與土體的幾何關(guān)系：設(shè)樁體半徑為r_p=D/2，土體計(jì)算區(qū)域水平寬度為B_w=w_p（其中w_p為等效寬度），垂直深度為B_h=h_p（其中h_p為等效深度）。樁體與土體的接觸界面可表示為樁身圓柱表面，其方程可記為：

|(x-x_c)^2+(y-y_c)^2|=r_p,其中(x_c,y_c)為樁頂中心坐標(biāo)。土體區(qū)域的邊界的數(shù)學(xué)定義則根據(jù)邊界約束條件（如位移約束）具體給出。3.3.1樁體參數(shù)設(shè)定在提取樁土相互作用的研究中，主要關(guān)注樁的樁長(zhǎng)、直徑、材料強(qiáng)度、樁身剛度等關(guān)鍵意義上。樁的寬度（或直徑）常設(shè)為~0.3m，長(zhǎng)度則視實(shí)際工程設(shè)定，并非一成不變。樁的材料為工程用混凝土，需準(zhǔn)確反映材料的物理力學(xué)性質(zhì)，如彈性模量E、泊松比μ、動(dòng)力模量λ以及材料的抗壓強(qiáng)度fc和抗拉強(qiáng)度f接觸參數(shù)包括摩擦系數(shù)μ和彈性系數(shù)k，這兩個(gè)參數(shù)在樁土界面的設(shè)置上尤為重要。通常，樁土間的界面摩擦系數(shù)和彈性系數(shù)會(huì)依據(jù)現(xiàn)場(chǎng)測(cè)試數(shù)據(jù)或經(jīng)驗(yàn)公式確定，而在數(shù)值模擬中常采用相同的值來簡(jiǎn)化問題。除此之外，樁土之間設(shè)定有法向接觸剛度與切向接觸剛度，這一點(diǎn)在模擬中可能需要邏輯性地設(shè)置為三種不同的情況以模擬樁的懸挑、半開放、全開挖等不同施工條件。樁土接觸面選用的模型類型不同，接觸力的計(jì)算方式也各異。書面表達(dá)時(shí)，應(yīng)明白提及模型接觸力的計(jì)算公式。如若模型采用罰力法來計(jì)算接觸力，須合理選擇罰系數(shù)F及與接觸剛度有關(guān)的彈性因子。樁頂位移及軸力監(jiān)測(cè)也為關(guān)鍵研究參數(shù)，樁頂節(jié)點(diǎn)被賦予一定的初始豎向壓力，模擬樁加載過程，通過監(jiān)測(cè)樁頭位移響應(yīng)來分析樁基穩(wěn)定性和承載力。節(jié)點(diǎn)法向位移應(yīng)設(shè)置為樁頂位移值，確保樁頂邊界條件符合實(shí)際情況。樁的實(shí)物模型參數(shù)需通過現(xiàn)場(chǎng)實(shí)測(cè)和實(shí)驗(yàn)室試驗(yàn)獲得，這些參數(shù)將對(duì)上述樁土接觸的特性、力學(xué)行為等進(jìn)行實(shí)時(shí)調(diào)整，從而確保模擬結(jié)果的科學(xué)性和準(zhǔn)確性。參數(shù)表數(shù)據(jù)格式示例下頁(見【表】)：【表】樁體參數(shù)表參數(shù)單位數(shù)值備注樁身寬度（直徑）m0.3樁長(zhǎng)m各段樁長(zhǎng)具體值材料彈性模量EPa具體數(shù)值數(shù)據(jù)泊松比μ無量綱具體數(shù)值數(shù)據(jù)動(dòng)力模量λPa具體數(shù)值數(shù)據(jù)抗壓強(qiáng)度f_cMPa具體數(shù)值數(shù)據(jù)抗拉強(qiáng)度f_tMPa具體數(shù)值數(shù)據(jù)水平土壓力系數(shù)k無量綱具體數(shù)值數(shù)據(jù)樁徑與樁周土的極限水平抗力系數(shù)值之比摩擦系數(shù)μ無量綱具體數(shù)值數(shù)據(jù)法向接觸剛度knN/m具體數(shù)值數(shù)據(jù)切向接觸剛度ktN/m具體數(shù)值數(shù)據(jù)樁頂豎向壓力PkN具體數(shù)值數(shù)據(jù)樁頂豎向靜荷載數(shù)值模擬工作中，樁體參數(shù)的設(shè)定需結(jié)合實(shí)際的工程條件，并通過反復(fù)迭代和校準(zhǔn)，確保其數(shù)值的可信度。此外樁體參數(shù)的選取應(yīng)結(jié)合工程項(xiàng)目的經(jīng)驗(yàn)數(shù)據(jù)和試驗(yàn)實(shí)測(cè)結(jié)果，使模擬結(jié)果更貼近工程實(shí)際情況。通過上述參數(shù)的精確設(shè)定，極大地提高了數(shù)值模擬結(jié)果的精確性和實(shí)用價(jià)值，為樁基的承載特性研究提供了強(qiáng)大工具支持。3.3.2土層分組與參數(shù)賦值在離散元數(shù)值模擬中，土層分組與參數(shù)賦值對(duì)于準(zhǔn)確反映樁土相互作用及樁基承載特性至關(guān)重要。為了更好地模擬實(shí)際工程地質(zhì)條件，本文將研究區(qū)域內(nèi)的土層依據(jù)其工程特性及物理力學(xué)參數(shù)進(jìn)行科學(xué)分組。具體每組土層及其對(duì)應(yīng)參數(shù)通過室內(nèi)外試驗(yàn)結(jié)果進(jìn)行賦值，以保證模擬結(jié)果的合理性和可靠性。以下是土層分組及對(duì)應(yīng)主要參數(shù)：土層分組及參數(shù)見【表】。其中土層編號(hào)依據(jù)從上至下的地質(zhì)順序進(jìn)行排列，為了描述土層在模擬中的相互作用，引入了土層間的摩擦系數(shù)和孔隙比兩個(gè)關(guān)鍵參數(shù)。根據(jù)【表】參數(shù)，所有土層在模擬過程中均采用摩爾-庫侖本構(gòu)模型進(jìn)行描述。該模型能夠反映土體在剪切應(yīng)力作用下的變形和破壞特性，具體模型參數(shù)通過公式（3-1）進(jìn)行計(jì)算：【公式】：τ其中τ表示剪應(yīng)力，σ表示正應(yīng)力，φ表示土層摩擦角（°），c表示土層粘聚力(Pa)。通過上述土層分組與參數(shù)賦值，離散元模型能夠更加精確地模擬樁土相互作用和樁基承載特性，為實(shí)際工程提供有力支撐。3.4接觸面模型與邊界條件在考慮樁土相互作用的欠送樁基承載特性離散元數(shù)值模擬中，接觸面模型和邊界條件的設(shè)定是極為關(guān)鍵的環(huán)節(jié)。接觸面模型：接觸面模型主要模擬樁與土之間的相互作用，由于樁土之間的接觸是一種復(fù)雜的物理現(xiàn)象，涉及摩擦、黏附和滑移等多種行為，因此需采用適當(dāng)?shù)慕佑|模型來準(zhǔn)確描述這種相互作用。常見的接觸面模型包括彈性接觸模型、庫侖摩擦模型等。在離散元分析中，每個(gè)顆粒與樁表面的接觸均可看作是一個(gè)獨(dú)立的接觸單元，通過設(shè)定相應(yīng)的接觸參數(shù)（如法向剛度、切向剛度、摩擦系數(shù)等），來模擬樁土之間的力學(xué)響應(yīng)。邊界條件設(shè)定：邊界條件的設(shè)定直接影響到模擬結(jié)果的準(zhǔn)確性，在欠送樁基的承載特性模擬中，一般考慮以下幾種邊界條件：地基邊界條件：模擬實(shí)際工程中地基的約束條件，如土體的側(cè)限條件、底部固結(jié)條件等。這些條件會(huì)影響樁周的土體應(yīng)力分布和樁的承載性能。樁身約束條件：模擬樁身的實(shí)際受力情況，如樁身的軸向壓力、彎矩等。這些約束條件會(huì)影響到樁的變形特性和承載力。外部環(huán)境因素：如風(fēng)荷載、地下水狀況等，這些環(huán)境因素也會(huì)對(duì)樁土相互作用產(chǎn)生影響，需在模擬中加以考慮。具體的邊界條件設(shè)定應(yīng)結(jié)合工程實(shí)際情況進(jìn)行，例如，在模擬分析中可通過設(shè)定不同深度的固定節(jié)點(diǎn)或施加特定荷載來模擬不同的邊界條件和實(shí)際工程場(chǎng)景。此外為了更精確地模擬樁土相互作用，還可以考慮采用動(dòng)態(tài)邊界條件，即在模擬過程中根據(jù)實(shí)際情況不斷調(diào)整邊界條件，以更真實(shí)地反映樁土相互作用的變化過程。通過上述接觸面模型和邊界條件的合理設(shè)定，可以更準(zhǔn)確地模擬欠送樁基的承載特性，為工程設(shè)計(jì)和施工提供有力的理論支持。3.4.1樁土接觸面參數(shù)標(biāo)定在進(jìn)行離散元數(shù)值模擬時(shí)，為了準(zhǔn)確描述樁與土之間的接觸關(guān)系，需要對(duì)樁土接觸面的參數(shù)進(jìn)行合理的標(biāo)定。這一過程通常包括以下幾個(gè)步驟：首先選擇合適的接觸模型來描述樁與土之間的相互作用，常見的接觸模型有彈性-黏性接觸模型和半剛性接觸模型等。這些模型通過調(diào)整接觸系數(shù)（如摩擦系數(shù)、粘聚力）來控制樁土之間的相對(duì)運(yùn)動(dòng)和變形。其次在設(shè)定接觸面參數(shù)時(shí)，需要根據(jù)實(shí)際工程條件和經(jīng)驗(yàn)數(shù)據(jù)進(jìn)行優(yōu)化。例如，對(duì)于軟弱地層中的樁，由于土體的壓縮性和排水性能較差，應(yīng)選取較大的摩擦系數(shù)以減少土體的剪切滑移；而在硬質(zhì)巖地層中，由于巖石的高強(qiáng)度和低變形能力，可以適當(dāng)降低摩擦系數(shù)。此外還需要校核所選接觸模型的物理意義和數(shù)學(xué)性質(zhì)是否符合實(shí)際情況?？梢酝ㄟ^計(jì)算不同工況下的接觸應(yīng)力分布內(nèi)容或位移場(chǎng)，驗(yàn)證接觸模型的有效性。如果發(fā)現(xiàn)某些參數(shù)設(shè)置導(dǎo)致不合理的接觸行為，則需重新調(diào)整直至滿足預(yù)期效果。在進(jìn)行樁土接觸面參數(shù)標(biāo)定時(shí)，既要確保接觸模型的科學(xué)合理性，又要結(jié)合實(shí)際工程情況進(jìn)行細(xì)致調(diào)整，從而提高數(shù)值模擬結(jié)果的精度和可靠性。3.4.2模型邊界約束方案在考慮樁土相互作用的欠送樁基承載特性離散元數(shù)值模擬中，模型邊界約束方案的設(shè)計(jì)至關(guān)重要，它直接影響到模擬結(jié)果的準(zhǔn)確性和可靠性。本節(jié)將詳細(xì)介紹模型邊界約束的具體方案。?邊界條件設(shè)置樁端約束：樁基的樁端通常采用剛性地基梁或樁帽進(jìn)行約束，限制其水平位移和豎向沉降。具體而言，樁端約束可以設(shè)置為固定約束（固定所有自由度）或滑動(dòng)約束（僅允許水平方向滑動(dòng)）。土體邊界約束：土體邊界約束包括水平向和豎直方向的約束。對(duì)于水平方向，土體邊界可以設(shè)置為固定約束或滑動(dòng)約束；對(duì)于豎直方向，土體邊界通常設(shè)置為不約束狀態(tài)，允許土體自由沉降。邊界條件選擇：根據(jù)具體工程情況和分析需求，可以選擇不同的邊界條件組合。例如，在模擬淺基礎(chǔ)承載力時(shí)，可以采用固定端約束；在分析深基礎(chǔ)承載力時(shí)，可以采用滑動(dòng)端約束。?具體實(shí)現(xiàn)方法剛性地基梁約束：在地基梁與樁的連接處設(shè)置剛性地基梁，通過設(shè)置剛性地基梁的約束條件，限制其水平位移和豎向沉降。樁帽約束：在樁頂設(shè)置樁帽，通過設(shè)置樁帽的約束條件，限制其水平位移和豎向沉降?；瑒?dòng)約束實(shí)現(xiàn)：在土體邊界設(shè)置滑動(dòng)約束，通過設(shè)置土體的滑動(dòng)參數(shù)，允許土體在水平方向上自由滑動(dòng)。?數(shù)值模擬步驟建立有限元模型：根據(jù)工程實(shí)際情況，建立有限元模型，包括樁、土體和邊界條件的設(shè)置。施加荷載：根據(jù)工程荷載情況，施加相應(yīng)的荷載，包括豎向荷載和水平荷載。求解方程：采用合適的求解器，對(duì)方程組進(jìn)行求解，得到樁基的承載力和變形響應(yīng)。后處理：對(duì)求解結(jié)果進(jìn)行后處理，包括繪制荷載-位移曲線、應(yīng)力分布內(nèi)容等，以評(píng)估樁基的承載特性。通過合理的模型邊界約束方案設(shè)計(jì)，可以有效模擬樁土相互作用下的欠送樁基承載特性，