胡均平，李科軍

?

考慮打樁反彈的組合錘頭沖擊過(guò)程的動(dòng)態(tài)特性分析

胡均平，李科軍

(中南大學(xué)機(jī)電工程學(xué)院，湖南長(zhǎng)沙，410083)

針對(duì)打樁過(guò)程中樁體容易損壞的現(xiàn)象，提出一種在錘體內(nèi)腔填充質(zhì)量塊的組合錘頭打樁的新型錘擊模型。運(yùn)用鍵合圖理論對(duì)錘擊系統(tǒng)建模，為模擬打樁反彈時(shí)各質(zhì)量塊與墊層接觸、分離的非線(xiàn)性特性，引入開(kāi)關(guān)類(lèi)量的鍵合圖元件功率結(jié)型結(jié)構(gòu)(SPJ)，采用集中參數(shù)搭建錘擊系統(tǒng)的鍵合圖模型，據(jù)此建立可統(tǒng)一表達(dá)系統(tǒng)所有運(yùn)動(dòng)狀態(tài)的數(shù)學(xué)模型。利用量綱一分析法推導(dǎo)錘擊力的數(shù)值解，并基于Matlab軟件仿真分析填充質(zhì)量塊和墊層參數(shù)對(duì)錘擊力的影響。通過(guò)試驗(yàn)驗(yàn)證該模型的準(zhǔn)確性。研究結(jié)果表明：采用組合錘頭時(shí)，錘擊力波形出現(xiàn)2個(gè)峰值，且峰值減小，這不僅可以保護(hù)樁體，而且可提高能量傳遞效率；各參數(shù)中，填充質(zhì)量塊和錘墊對(duì)錘擊性能的影響最大，砧坐墊的影響次之，樁墊的影響最小。該模型可用于錘擊系統(tǒng)的參數(shù)優(yōu)化研究及動(dòng)態(tài)特性分析。

組合錘頭；鍵合圖；反彈；錘擊力

液壓打樁錘以液壓能為動(dòng)力，提起錘頭后通過(guò)快速泄油或反向供油使錘頭加速下降，將樁體打入土 壤[1]。與柴油錘、振動(dòng)錘和壓樁機(jī)等樁工機(jī)械相比，液壓錘是一種振動(dòng)小、噪聲低、無(wú)油煙污染和施工適用性強(qiáng)的產(chǎn)品，而且錘擊頻率、錘擊行程可控，在樁基施工中發(fā)揮越來(lái)越重要的作用[2?3]。錘擊系統(tǒng)作為打樁錘的核心裝置，其工作性能直接關(guān)系到打樁質(zhì)量及打樁效率。為掌握樁錘沖擊過(guò)程的動(dòng)態(tài)特性對(duì)沉樁效率的影響規(guī)律，有必要對(duì)錘擊系統(tǒng)建立完整、準(zhǔn)確的數(shù)學(xué)模型。TAKE等[4?5]用集中質(zhì)量塊表示樁錘、砧座，用無(wú)質(zhì)量彈簧表示錘墊、樁墊，建立普通樁錘錘擊力的解析模型，并仿真分析了墊層參數(shù)等對(duì)錘擊力的影響。方治華等[6]基于ABAQUS軟件對(duì)混凝土預(yù)制樁的打樁過(guò)程進(jìn)行動(dòng)態(tài)模擬，得到了打樁機(jī)械參數(shù)、樁體參數(shù)對(duì)打樁過(guò)程的影響規(guī)律。羅春雷等[7]針對(duì)振動(dòng)樁錘在沉樁作業(yè)過(guò)程中的“跳機(jī)”問(wèn)題，建立了分段非線(xiàn)性的二自由度模型，對(duì)“跳機(jī)”現(xiàn)象產(chǎn)生的根本原因進(jìn)行了仿真分析和試驗(yàn)研究。為提高樁錘沖擊的貫入量、增加墊層使用壽命及使樁體不受損壞，楊永海等[8?10]研究了金屬顆粒填充樁錘的錘擊模型，發(fā)現(xiàn)金屬顆粒的加入能引發(fā)二次沖擊，生成矩形錘擊力波形，這可以降低錘擊力最大峰值及提高沖擊過(guò)程的能量轉(zhuǎn)化率，但這些研究中，有的對(duì)錘擊模型建模時(shí)忽略了沖擊過(guò)程中金屬顆粒本身的影響[8]，有的采用接觸剛度模擬金屬顆粒與錘體的相互作用，但形狀不規(guī)則的金屬顆粒給接觸剛度的計(jì)算及設(shè)計(jì)帶來(lái)了困難[9]，而且都沒(méi)有考慮打樁反彈時(shí)各質(zhì)量塊與墊層之間接觸、分離的非線(xiàn)性特性對(duì)錘擊系統(tǒng)工作性能的影響。為此，本文作者用質(zhì)量塊代替金屬顆粒填充到錘體內(nèi)腔，并在質(zhì)量塊和錘體內(nèi)腔的上、下接觸面之間放置墊層，組成一種新型錘擊模型。為模擬打樁反彈時(shí)錘擊系統(tǒng)從一個(gè)工作狀態(tài)過(guò)渡到另外一個(gè)工作狀態(tài)的非連續(xù)性工作過(guò)程,引入開(kāi)關(guān)類(lèi)量的鍵合圖元件功率結(jié)型結(jié)構(gòu)(SPJ)，采用集中參數(shù)建立可統(tǒng)一表達(dá)錘擊系統(tǒng)所有運(yùn)動(dòng)狀態(tài)的鍵合圖模型，并利用該模型仿真分析填充質(zhì)量塊和墊層參數(shù)對(duì)錘擊力的影響規(guī)律，以便為錘擊系統(tǒng)的設(shè)計(jì)和優(yōu)化提供依據(jù)。

1 功率結(jié)型結(jié)構(gòu)

現(xiàn)實(shí)中，研究的系統(tǒng)多為混雜系統(tǒng)，即系統(tǒng)中同時(shí)存在連續(xù)元件和離散元件。建立混雜系統(tǒng)鍵合圖模型的關(guān)鍵就是對(duì)系統(tǒng)中具有離散特性的元件尤其是開(kāi)關(guān)元件進(jìn)行鍵合圖描述。KARNOPP等[11]提出用轉(zhuǎn)換器MTF來(lái)描述系統(tǒng)中的離散行為，當(dāng)開(kāi)關(guān)處于閉合狀態(tài)時(shí)，MTF的調(diào)制參數(shù)取值為1；當(dāng)開(kāi)關(guān)處于斷開(kāi)狀態(tài)時(shí)，調(diào)制參數(shù)取值為0，但在開(kāi)關(guān)元件處于斷開(kāi)狀態(tài)時(shí)，結(jié)點(diǎn)和系統(tǒng)的連接很容易斷開(kāi)。BORUTZKY[12]在KARNOPP等[11]的研究基礎(chǔ)上，提出增加1個(gè)阻性元件調(diào)制轉(zhuǎn)換器MTF連接以描述系統(tǒng)離散行為。該方法能保證鍵合圖模型的因果關(guān)系不隨物理結(jié)構(gòu)的變化而變化，但的引入帶來(lái)了損耗，使得開(kāi)關(guān)元件不再呈理想特性，而且的變化范圍極大，容易導(dǎo)致系統(tǒng)模型在數(shù)值上呈剛性，增加仿真計(jì)算時(shí)間。為克服上述缺點(diǎn)，UMARIKAR等[13]提出一種表達(dá)開(kāi)關(guān)類(lèi)量的新鍵合圖元件1s結(jié)，即功率結(jié)型結(jié)構(gòu)(switch power junction, SPJ)。該元件分別在基本結(jié)型元件1結(jié)和0結(jié)的基礎(chǔ)上，增加了2個(gè)相互排斥的流通口和勢(shì)通口，通過(guò)信號(hào)鍵控制2個(gè)相互排斥的流通口或勢(shì)通口的生效。1S和0S結(jié)因果示意圖如圖1所示。從圖1(a)可見(jiàn)：1s結(jié)連接了2個(gè)流結(jié)，并且由1對(duì)布爾變量控制；變量1與鍵1相連，變量2與鍵2相連；系統(tǒng)的信號(hào)鍵決定了與布爾變量相連的鍵是否生效。當(dāng)1為1時(shí)，鍵1生效，1(流變量)為1s結(jié)的唯一流輸入；當(dāng)2為1時(shí)，鍵2生效，2(流變量)為1s結(jié)的唯一流輸入。因此，1s結(jié)上的流為11+22，通過(guò)1s結(jié)的各勢(shì)變量之和為零，組成關(guān)系式如下：

(2)

(3)

類(lèi)似地，從圖1(b)可見(jiàn)0s結(jié)上有2個(gè)勢(shì)輸入，0s結(jié)上的勢(shì)為11+22。通過(guò)0s結(jié)的各流變量之和為0，組成關(guān)系式如下：

(5)

(7)

(8)

根據(jù)需要，運(yùn)用布爾算法，可以得1，2，…，u等個(gè)變量推導(dǎo)出1s結(jié)和0s結(jié)上的組成關(guān)系式。

(a) 1s結(jié)；(b) 2s結(jié)

2 錘擊系統(tǒng)建模

假設(shè)樁體為無(wú)限長(zhǎng)桿，據(jù)一維波動(dòng)力學(xué)理 論[14?16]，樁頂?shù)腻N擊力與波速成正比，即=(式中：為樁的波阻，與樁的彈性模量p、樁的橫截面積p及波在樁內(nèi)的軸向傳播速度p有關(guān)，=pp/p)。因此，建模時(shí)樁體可以用1個(gè)阻尼器等效替代。

組合錘頭的錘擊模型如圖2所示。分別用不同的集中質(zhì)量塊表示填充質(zhì)量、錘體、砧座，用并聯(lián)的彈簧、阻尼器表示砧坐墊，用彈簧表示錘墊、樁墊，用阻尼器表示樁體。

圖2 組合錘頭錘擊模型

2.1 數(shù)學(xué)模型

結(jié)合前面的功率結(jié)型結(jié)構(gòu)及組合錘頭錘擊模型，建立圖3所示的錘擊系統(tǒng)鍵合圖模型。在錘頭沖擊過(guò)程中，重力相對(duì)于錘擊力要小得多，因此，建模時(shí)不予考慮。為避免因果矛盾，在1()結(jié)上，代表樁體等效阻尼器的元件上的勢(shì)變量為輸入、流變量為輸出。圖3中：r1為填充質(zhì)量塊的質(zhì)量；r2為錘體的質(zhì)量；a為砧座質(zhì)量；1，2和3分別為錘墊、砧坐墊、樁墊的剛度系數(shù)；2和3分別為砧坐墊阻尼系數(shù)和樁體阻尼系數(shù)。各慣性、容性及阻性元件都具有線(xiàn)性特征，利用各結(jié)約束條件及標(biāo)定的因果關(guān)系，可以推導(dǎo)出系統(tǒng)的狀態(tài)方程：

圖3 錘擊系統(tǒng)鍵合圖模型

式中：1，2和3分別為填充質(zhì)量塊、錘體、砧座的動(dòng)量；，，和分別為填充質(zhì)量塊、錘體、砧座阻尼器的運(yùn)動(dòng)速度。

由于填充質(zhì)量塊上下接觸面都放置了墊層，可以認(rèn)為打樁過(guò)程中填充質(zhì)量塊始終與彈簧1保持接觸，仿真計(jì)算時(shí)1取1，2取0。式(9)中其他2對(duì)布爾變量的取值和轉(zhuǎn)換條件如圖4所示，可知錘擊系統(tǒng)有22即4種運(yùn)動(dòng)狀態(tài)。需要指出的是：并不是每種運(yùn)動(dòng)狀態(tài)都會(huì)發(fā)生，但該式統(tǒng)一表達(dá)了錘擊過(guò)程中可能出現(xiàn)的所有運(yùn)動(dòng)狀態(tài)，且鍵合圖模型的因果關(guān)系在系統(tǒng)處于任何運(yùn)動(dòng)狀態(tài)時(shí)均保持不變。

圖4 錘擊系統(tǒng)運(yùn)動(dòng)狀態(tài)轉(zhuǎn)換

2.2 參數(shù)量綱一化分析

為便于研究各參數(shù)變化對(duì)錘擊系統(tǒng)的反彈和錘擊力影響，對(duì)式(9)的參數(shù)進(jìn)行量綱一化處理。選擇質(zhì)量r(r=r1+r2)、樁體等效阻尼3及錘擊初速度0為基本量，可得量綱一時(shí)間為*=3/r；量綱一質(zhì)量為，，；量綱一彈簧剛度為，，；量綱一力為，，；量綱一動(dòng)量為，，；量綱一阻尼為，；量綱一位移為，，，。將以上量綱一參數(shù)代入式(9)可得：

在樁錘撞擊樁體瞬間，錘體及腔內(nèi)填充質(zhì)量塊的初始速度均為0，對(duì)應(yīng)的量綱一速度為：，。砧座、樁體的初始速度、加速度均為0，因此，錘擊系統(tǒng)的初始條件為

(11)

3 打樁反彈影響

1—考慮打樁反彈；2—考慮打樁反彈。

(a) u3；(b) u5

4 參數(shù)對(duì)錘擊系統(tǒng)動(dòng)態(tài)特性的影響

利用建立的數(shù)學(xué)模型可研究各參數(shù)組合對(duì)錘擊系統(tǒng)動(dòng)態(tài)特性的影響。由狀態(tài)方程式(10)可知：樁頂?shù)腻N擊力主要由7個(gè)量綱一參數(shù)確定，即填充質(zhì)量塊質(zhì)量、錘體質(zhì)量、砧板質(zhì)量、錘墊剛度系數(shù)、砧坐墊剛度系數(shù)、阻尼系數(shù)和樁墊的剛度系數(shù)。

4.1 填充質(zhì)量塊參數(shù)

Mr1*：1—0；2—0.3；3—0.5；4—0.6；5—0.7。

4.2 墊層參數(shù)

K1*：1—2；2—4；3—8；4—16；5—32。

K2*：1—2；2—4；3—8；4—16；5—32。

R2*：1—0.1；2—0.2；3—0.4；4—0.8。

K3*：1—2；2—4；3—8；4—16；5—32。

5 試驗(yàn)

為驗(yàn)證組合錘頭打樁效果及數(shù)學(xué)模型的準(zhǔn)確性，采用湖南長(zhǎng)河機(jī)械有限公司生產(chǎn)的ZCYL20-6液壓打樁機(jī)對(duì)錘擊力進(jìn)行測(cè)量。填充質(zhì)量塊、錘體、砧座的質(zhì)量分別為3.0，3.0和1.5 t；錘墊、砧坐墊和樁墊剛度分別為1×109，4×109和2×109N/m；樁錘初速度為5.6 m/s。樁采用預(yù)應(yīng)力高強(qiáng)混凝土管(PHC)樁，AB型，外直徑為500 mm，壁厚為125 mm，長(zhǎng)為13 m，彈性模量約為36.0 GPa；應(yīng)力波在樁身中傳播的速度為3.8~4.0 km/s；樁體最大承載力為5.6 MN左右。普通錘頭和組合錘頭產(chǎn)生的錘擊力見(jiàn)圖12。從圖12可見(jiàn)：普通錘頭錘擊力超過(guò)了樁體承載力，容易損壞樁體，而組合錘頭錘擊力峰值要小，始終沒(méi)有穿越樁體最大承載力所在的參照線(xiàn)，保護(hù)了樁體，且錘擊力波形接近矩形，增加了有效工作持續(xù)時(shí)間；組合錘頭錘擊力的計(jì)算值和實(shí)測(cè)值具有相同的變化趨勢(shì)，2次沖擊的錘擊力峰值也較接近，表明建立的數(shù)學(xué)模型是正確、可行的。

1—普通錘頭計(jì)算值；2—樁體最大承載力；3—組合錘頭計(jì)算值；4—組合錘頭實(shí)測(cè)值。

6 結(jié)論

1) 提出一種采用組合錘頭打樁的新型錘擊模型,為模擬沖擊過(guò)程中打樁反彈的運(yùn)動(dòng)狀態(tài)，引入開(kāi)關(guān)類(lèi)量的鍵合圖元件功率結(jié)型結(jié)構(gòu)(SPJ)，采用集中參數(shù)建立了全面反映系統(tǒng)運(yùn)動(dòng)狀態(tài)的鍵合圖模型；比較分析了考慮打樁反彈、不考慮打樁反彈2種數(shù)學(xué)模型的錘擊力仿真曲線(xiàn)。所提出的模型能更好地模擬打樁的錘擊機(jī)理。

2) 與普通錘頭相比，組合錘頭能產(chǎn)生2次沖擊，錘擊力峰值減小，錘擊力波形更接近矩形，有效沖擊時(shí)間增加，這不僅可以保護(hù)樁體，而且能提高能量傳遞效率；各參數(shù)中，填充質(zhì)量塊和錘墊對(duì)錘擊性能的影響最大，砧坐墊的影響次之，樁墊的影響最小。

3)所建立的數(shù)學(xué)模型能較準(zhǔn)確地反映各參數(shù)對(duì)錘擊力的影響規(guī)律，對(duì)完善錘擊系統(tǒng)設(shè)計(jì)理論和提高打樁效率有參考價(jià)值。

[1] 嚴(yán)東兵, 胡均平, 梅麗仿, 等. 氮爆式液壓打樁錘結(jié)構(gòu)原理與分析[J]. 工程機(jī)械, 2009, 40(4): 35?38. YAN Dongbing, HU Junping, MEI Lifang, et al. Structural principle and analysis of nitrogen-exploding type hydraulic pile hammer[J]. Construction Machinery and Equipment, 2009, 40(4): 35?38.

[2] GUO Y, HU J P, ZHANG L Y. Finite-element analysis of multi-body contacts for pile driving using a hydraulic pile hammer[J]. Proceedings of the Institution of Mechanical Engineers: Part C, 2011, 225(5): 1153?1161.

[3] 劉成沛. 液壓打樁錘壓力反饋配流系統(tǒng)及回油特性研究[D]. 長(zhǎng)沙: 中南大學(xué)機(jī)電工程學(xué)院, 2012: 1?10. LIU Chengpei. Research on pressure feedback flow distribution system and oil return characteristics of hydraulic pile hammer[D]. Changsha: Central South University. School of Mechanical and Electrical Engineering, 2012: 1?10.

[4] TAKE W A, VALSANGKAR A J, RANDOLPH M F. Analytical solution for pile hammer impact[J]. Computers and Geotechnics, 1999, 25(2): 57?74.

[5] 謝永建, 朱合華, 王懷忠, 等. 打樁錘擊模型的解析算法[J]. 巖石力學(xué)與工程學(xué)報(bào), 2005, 24(1): 171?176. XIE Yongjian, ZHU Hehua, WANG Huaizhong, et al. Analytical solution for model of pile hammer impact[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(1): 171?176.

[6] 方治華, 程建華. 打樁過(guò)程的非線(xiàn)性有限元分析[J]. 內(nèi)蒙古科技大學(xué)學(xué)報(bào), 2009, 128(2): 184?188. FANG Zhihua, CHENG Jianhua. The non-linear finite element analysis in the process of piling[J]. Journal of Inner Mongolia University of Science and Technology, 2009, 128(2): 184?188.

[7] 羅春雷, 吳偉傳, 丁吉, 等. 振動(dòng)沉樁過(guò)程“跳機(jī)”現(xiàn)象動(dòng)力學(xué)特性[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 43(11): 4274?4279. LUO Chunlei, WU Weichuan, DING Ji, et al. Dynamic characteristics of vibratory hammer “dancing” in piling[J]. Journal of Central South University (Science and Technology), 2012, 43(11): 4274?4279.

[8] 楊永海, 呂景忠, 隋振, 等. 加速下落沖擊式液壓樁錘設(shè)計(jì)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2005, 36(5): 150?152. YANG Yonghai, Lü Jingzhong, SUI Zhen, et al. Design of Accelerated Impact Pile Hammer[J]. Transactions of the Chinese Society for Agricultural Machinery, 2005, 36(5): 150?152.

[9] 胡均平, 郭勇, 宋光偉, 等. 金屬顆粒填充樁錘錘擊模型的數(shù)值分析[J]. 中國(guó)公路學(xué)報(bào), 2010, 23(2): 122?126. HU Junping, GUO Yong, SONG Guangwei, et al. Numerical analysis for impact model of pile hammer filled with metal particles[J]. China Journal of Highway and Transport, 2010, 23(2): 122?126.

[10] 劉德順, 李夕兵, 朱萍玉. 沖擊機(jī)械動(dòng)力學(xué)與反演設(shè)計(jì)[M]. 北京: 科學(xué)出版社, 2007: 38?75.LIU Deshun, LI Xibing, ZHU Pingyu. Dynamics and inverse design of impact machines[M]. Beijing: Science Press, 2007: 38?75.

[11] KARNOPP D C, MARGOLIS D L, ROSENBERG R C. System dynamics: modeling, simulation, and control of mechatronic systems[M]. New York: John Wiley & Sons, 2012: 17?33.

[12] BORUTZKY W. Bond graph methodology: development and analysis of multidisciplinary dynamic system models[M]. Berlin: Springer, 2010: 305?349.

[13] UMARIKAR A C, UMMLAND L. Modeling of switching systems in bond graphs using the concept of switched power junctions[J]. Journal of the Franklin Institute, 2005, 342(2): 131?147.

[14] 苗潤(rùn)田. 基于波動(dòng)力學(xué)的鑿巖機(jī)效率研究[D]. 長(zhǎng)沙: 中南大學(xué)機(jī)電工程學(xué)院, 2006: 4?22. MIAO Shuntian. The research of rock drill efficiency based on the wave mechanics[D]. Changsha: Central South University. School of Mechanical and Electrical Engineering, 2006: 4?22.

[15] 雷林源. 樁基動(dòng)力學(xué)[M]. 北京: 冶金工業(yè)出版社, 2000: 35?41. LEI Linyuan. Pile foundation dynamics[M]. Beijing: Metallurgical Industry Press, 2000: 35?41.

[16] 劉德順, 楊志高, 楊書(shū)儀. 電錘沖擊系統(tǒng)的波動(dòng)力學(xué)研究[J]. 振動(dòng)工程學(xué)報(bào), 2013, 26(6): 864?870. LIU Deshun, YANG Zhigao, YANG Shuyi. Wave mechanics for an electric hammer impact system[J]. Journal of Vibration Engineering, 2013, 26(6): 864?870.

(編輯 陳燦華)

Dynamic characteristics analysis of impact model of composite hammer considering rebound in piling

HU Junping, LI Kejun

(School of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China)

In view of the phenomenon of pile body damage in the process of pile driving, the impact model of composite hammer filled with mass in the inner cavity was proposed. The lumped parameter method was used to establish the bond graph model of impact model, which introduced the concept of switched power junctions (SPJ) to represent the nonlinear characteristics of contact and separation between masses and cushions, and then the system equations were derived from the model in a unified way. The method of dimensionless analysis was used to analyze the effect of various parameters on the impact force within the Matlab environment. The model was verified experimentally. The results show that the composite hammer causes a lower crest and a longer effective impact time as it has more than one force wave, which also protects the pile body and improves the efficiency of energy transfer. The filling mass and the hammer cushion have the greatest effects on the performance of the hammer, the anvil cushion is the second, and the pile cushion is the smallest. The method can be used to the optimum design and dynamic analysis of the hammer impact system.

composite hammer; bond graph; rebound; hammer impact force

10.11817/j.issn.1672?7207.2017.02.012

TH113.1；U415.5

A

1672?7207(2017)02?0348?07

2016?02?10；

2016?04?22

國(guó)家自然科學(xué)基金資助項(xiàng)目(51175518)；湖南省科技型中小企業(yè)技術(shù)創(chuàng)新基金資助項(xiàng)目(12C26214305029) (Project(51175518) supported by the National Natural Science Foundation of China; Project(12C26214305029) supported by Small and Mediu-sized Enterprise Technology Innovation Fund of Hunan Province)

李科軍，博士研究生，從事液壓機(jī)械與控制研究；E-mail：likejuncsu@126.com