鄭　聰， 趙　威， 張海豐， 馬保松

(1. 中國(guó)地質(zhì)大學(xué)(武漢)工程學(xué)院， 湖北 武漢　430074；

2. 中交第二公路勘察設(shè)計(jì)研究院有限公司， 湖北 武漢　430056)

?

布置參數(shù)對(duì)TBM邊緣滾刀受力的影響

鄭聰1， 趙威2， 張海豐1， 馬保松1

(1. 中國(guó)地質(zhì)大學(xué)(武漢)工程學(xué)院， 湖北 武漢430074；

2. 中交第二公路勘察設(shè)計(jì)研究院有限公司， 湖北 武漢430056)

摘要：邊緣滾刀安裝在全斷面隧道掘進(jìn)機(jī)(TBM)刀盤的最外緣，在破巖時(shí)邊緣滾刀的受力情況十分復(fù)雜，工作條件惡劣，容易損壞。為研究布置參數(shù)對(duì)邊緣滾刀的影響，采用有限元分析軟件ABAQUS建立邊緣滾刀破巖的非線性動(dòng)力學(xué)模型，分析邊緣滾刀的受力情況。結(jié)果表明： 邊緣滾刀在破巖時(shí)所受的作用力會(huì)劇烈振蕩變化； 安裝傾角和刀盤過渡圓弧半徑對(duì)邊緣滾刀受力的影響顯著； 邊緣滾刀受力隨著安裝傾角和刀盤過渡圓弧半徑的增加而呈現(xiàn)近似線性的增長(zhǎng)。

關(guān)鍵詞：TBM； 邊緣滾刀； 安裝傾角； 過渡圓弧半徑； 有限元

0引言

盤形滾刀是全斷面隧道掘進(jìn)機(jī)(TBM)刀盤上的主要切削刀具。根據(jù)滾刀在刀盤上的安裝位置，將滾刀分為中心滾刀、正滾刀和邊緣滾刀。邊緣滾刀安裝在刀盤外緣，具有一定的安裝傾角。由于其位置特殊、受力情況十分復(fù)雜、線速度大和工作條件惡劣，故邊緣滾刀易損壞。對(duì)邊緣滾刀破巖過程中的受力情況進(jìn)行分析，可以幫助設(shè)計(jì)人員優(yōu)化邊緣滾刀的布置參數(shù)。

國(guó)內(nèi)學(xué)者對(duì)盤形滾刀進(jìn)行了大量的研究。夏毅敏等[1]采用遺傳算法分析了復(fù)合式土壓平衡盾構(gòu)盤形滾刀布置的基本規(guī)律； 沈斌[2]從破巖機(jī)理的角度研究了滾刀破巖的關(guān)鍵影響因素，提出了邊緣滾刀布置的優(yōu)化方法； 宋克志等[3]分析了刀盤上滾刀推力分布規(guī)律，得出刀盤上不同位置的滾刀推力差異較大； 張厚美[4]總結(jié)滾刀受力的計(jì)算方法，根據(jù)各種計(jì)算方法，編寫了掘進(jìn)性能預(yù)測(cè)程序，用于TBM刀具布置優(yōu)化； 吳元等[5]從破巖效率的角度出發(fā)，分析了邊緣滾刀刀刃角等關(guān)鍵參數(shù)的影響。但目前研究邊緣滾刀破巖過程受力變化情況的文獻(xiàn)較少。

本文采用ABAQUS軟件建立不同安裝傾角和刀盤過渡圓弧半徑條件下的邊緣滾刀破巖的非線性動(dòng)力學(xué)模型，研究不同條件下邊緣滾刀所受垂直力、滾動(dòng)力和扭矩的變化規(guī)律。

1邊緣滾刀的位置參數(shù)及巖石模型

1.1TBM刀具類型及位置

圖1為刀盤上滾刀刀具類型與安裝位置示意圖。邊緣滾刀安裝在刀盤最外緣的過渡圓弧區(qū)域，且有一定的安裝傾角。邊緣滾刀距離刀盤中心較遠(yuǎn)，掘進(jìn)時(shí)繞刀盤中心點(diǎn)公轉(zhuǎn)的半徑較大，在破碎巖石時(shí)線速度較大，工作條件差，受力情況復(fù)雜，也更容易損壞[6-7]。邊緣滾刀與正滾刀在刀盤上的分布位置及具體參數(shù)如圖2所示(其中，R為過渡圓弧中心到刀盤中心軸線的距離，r為過渡圓弧半徑，θ為滾刀安裝傾角)。

圖1　滾刀刀具類型與位置示意圖

圖2　邊緣滾刀位置參數(shù)

1.2滾刀模型

滾刀模型采用使用最廣泛的17 in盤形滾刀(半徑r0=216 mm)為原型，并進(jìn)行一定的簡(jiǎn)化。圖3為滾刀模型剖面圖。滾刀一般采用合金鋼材料，其強(qiáng)度和硬度遠(yuǎn)大于巖石。在仿真過程中，不考慮滾刀磨損，將滾刀設(shè)定為剛體，采用各向同性的線彈性材料來模擬邊緣滾刀材料。滾刀與巖石材料的各項(xiàng)物理力學(xué)參數(shù)見表1。

圖3　滾刀模型剖面圖(單位： mm)

1.3巖石模型

本文主要考慮邊緣滾刀與巖石的相互作用，而邊緣滾刀安裝于刀盤過渡圓弧段，所以選取與刀盤過渡圓弧段相適應(yīng)的具有一定厚度的弧形巖石模型。巖石材料的力學(xué)性質(zhì)十分復(fù)雜，其應(yīng)力應(yīng)變分為3個(gè)階段： 彈性階段、塑性階段和斷裂損傷階段[8]。巖石的應(yīng)力-應(yīng)變曲線如圖 4所示：AB段材料處于彈塑性變形階段，材料的應(yīng)變不斷增大，應(yīng)力也不斷增大； 當(dāng)應(yīng)力增加到B點(diǎn)時(shí)，即D=0時(shí)(D為損傷因子)，材料開始進(jìn)入損傷階段，材料的應(yīng)變繼續(xù)增大而應(yīng)力逐漸減??； 直至C點(diǎn)，即D=1時(shí)，材料破壞

表1　滾刀與大理石材料參數(shù)

圖4　巖石的應(yīng)力-應(yīng)變曲線

本次仿真分析中，為了較好地模擬巖石材料的塑性本構(gòu)關(guān)系，依據(jù)德魯克-普拉格塑性屈服準(zhǔn)則，采用包含單元?jiǎng)h除功能的損傷失效準(zhǔn)則模擬巖石材料的屈服和破壞。擴(kuò)展線性德魯克-普拉格模型在π平面上的屈服軌跡如圖 5所示，屈服函數(shù)為

F=t-ptanβ-d=0。

(1)

式中： t為偏應(yīng)力； p為等效壓應(yīng)力； β為屈服面在p-t應(yīng)力空間上的傾角，等于材料的摩擦角； d為屈服面在應(yīng)力空間t軸上的截距，等于材料的黏聚力。

以大理石作為本次模擬中巖石的材料原型，大理石的物理力學(xué)參數(shù)見表1[9]。

2邊緣滾刀切削巖石的動(dòng)力學(xué)模型

圖6為邊緣滾刀破巖的動(dòng)力學(xué)模型。在考慮實(shí)際工況及數(shù)值模擬效率的情況下，設(shè)定切削深度h=5mm。參考點(diǎn)N(如圖6(b)所示)設(shè)在滾刀的質(zhì)心處，建立其與整個(gè)滾刀的剛體約束關(guān)系，對(duì)N施加約束及邊界條件來控制整個(gè)滾刀的運(yùn)動(dòng)。滾刀隨切削刀盤轉(zhuǎn)動(dòng)時(shí)的運(yùn)動(dòng)軌跡是一小段圓弧，由于邊緣滾刀距刀盤中心點(diǎn)較遠(yuǎn)，在仿真中以一小段直線代替這段圓弧，即假定滾刀的運(yùn)動(dòng)軌跡是直線。對(duì)N施加沿z軸負(fù)方向的平動(dòng)速度v和繞自身中軸線的轉(zhuǎn)動(dòng)角速度ω，將其他方向的自由度全部固定。固定巖石模型中除開挖面之外的所有平面的自由度，放開開挖面的自由度。滾刀與巖石之間采用非光滑的接觸來模擬，法向性質(zhì)為允許分開的硬接觸，切向性質(zhì)為采用罰函數(shù)的摩擦算法，摩擦因數(shù)為0.325[10]。

圖5　擴(kuò)展線性德魯克-普拉格模型p-t關(guān)系曲線

(a) 平面圖

(b) 三維圖

本次模擬中設(shè)定過渡圓弧中心到刀盤中心軸線的距離R=2 150mm，轉(zhuǎn)速ω0=1.026 3rad/s[11]。為了研究安裝傾角和刀盤過渡圓弧半徑對(duì)邊緣滾刀切削巖石的影響，設(shè)置4種不同的安裝傾角θ=15°、30°、45°和60°及4種不同的刀盤過渡圓弧半徑r=200、400、600、800mm。假定滾刀與巖石之間無相對(duì)滑移，即滾刀運(yùn)動(dòng)為純滾動(dòng)，由此得滾刀的平動(dòng)速度v和轉(zhuǎn)動(dòng)角速度ω的計(jì)算表達(dá)式為：

l=R+r·sinθ;

(2)

v=ω0·l;

(3)

ω=v/r0。

(4)

式中 l為邊緣滾刀中心點(diǎn)到刀盤軸線的距離，mm。

邊緣滾刀線速度及角速度的計(jì)算結(jié)果見表 2。

表2　邊緣滾刀線速度及角速度

3模擬結(jié)果分析

3.1邊緣滾刀受力動(dòng)態(tài)變化

如圖7所示，滾刀在沿著切削進(jìn)程方向的滾動(dòng)力Fr作用下向前滾動(dòng)破碎巖石，在受到垂直于巖石表面方向的垂直力Fv的作用下壓入巖石，并保持一定的切削深度；F1為巖石對(duì)滾刀刀刃頂部的作用力；F2和F3為巖石對(duì)滾刀刀刃側(cè)面的作用力；M為繞滾刀軸線的扭矩。

圖8是過渡圓弧半徑r=800 mm、邊緣滾刀安裝傾角θ=60°時(shí)，邊緣滾刀的垂直力、滾動(dòng)力和扭矩的變化情況。整個(gè)切削過程中，邊緣滾刀受力是不斷振蕩變化的。垂直力最大值為2 595 N； 滾動(dòng)力最大值為298 N，約為垂直力的1/10； 扭矩最大值為49 012 N·m。滾刀受力出現(xiàn)高頻率的振蕩，反映了切削巖石這一過程是不平穩(wěn)的，這是由巖石材料的物理力學(xué)性能決定的。巖石材料與土體材料不同，其具有一定的彈脆性，塑性變形的能力較小，往往在塑性變形不大的情況下突然發(fā)生破壞。滾刀壓入巖石后，巖石在外力作用下產(chǎn)生彈性變形，隨著滾刀壓入深度的不斷加大，產(chǎn)生的彈性變形越來越大，直至巖石材料屈服進(jìn)入塑性變形階段，此階段滾刀的受力是不斷增加的。巖石材料可產(chǎn)生的塑性變形量很小，所以在進(jìn)入塑性變形階段之后很快會(huì)破壞，滾刀受力大幅度減小。滾刀向前移動(dòng)不斷切削巖石，巖石不斷產(chǎn)生彈性變形，繼而塑性變形，直至破壞，使得滾刀受力曲線在整體上出現(xiàn)高頻率的陡增和回落。紀(jì)昌明等[12]認(rèn)為滾刀破巖時(shí)，巖石的破壞具有不斷階躍破碎的特點(diǎn)，本文的數(shù)據(jù)結(jié)果也驗(yàn)證了這一點(diǎn)。

(a)(b)

圖7滾刀受力示意圖

Fig. 7Stressing of disc cutter

研究邊緣滾刀破巖時(shí)的受力情況，除了應(yīng)分析滾刀整體受力的振蕩變化規(guī)律之外，也應(yīng)當(dāng)重點(diǎn)分析各項(xiàng)受力的最大值。為了分析在切削過程中的各向受力最大值的變化情況，作垂直力峰值、滾動(dòng)力峰值及扭矩峰值包絡(luò)線，如圖8所示。垂直力峰值為1 600～2 700 N，滾動(dòng)力峰值為180～330 N，扭矩的峰值為25 000～50 000 N·m。

3.2邊緣滾刀受力單因素分析

在邊緣滾刀的布置參數(shù)中，安裝傾角和刀盤過渡圓弧半徑對(duì)邊緣滾刀受力的影響較大。邊緣滾刀的受力情況是不斷變化的，分析中應(yīng)考慮整個(gè)過程中的受力平均值及受力峰值平均值。將整個(gè)受力過程中的受力取平均作為受力平均值； 取受力最大的20個(gè)峰值，將其取平均設(shè)定為受力峰值平均值。

3.2.1邊緣滾刀受力與安裝傾角的關(guān)系

圖9為邊緣滾刀不同安裝傾角下的垂直力、滾動(dòng)力和扭矩的變化情況。從圖9中可以看出，在相同的安裝傾角和刀盤過渡圓弧半徑下，邊緣滾刀受力的峰值平均值約為其平均值的2～3倍，這與李輝等[13]的研究結(jié)果是一致的。

(a) 垂直力

(b) 滾動(dòng)力

(c) 扭矩

圖8邊緣滾刀垂直力、滾動(dòng)力和扭矩與時(shí)間的關(guān)系曲線(r=800 mm，θ=60°)

Fig. 8Vertical stresses, rolling stresses and torques of gage disc cutter vs. time (r=800 mm，θ=60°)

在一定過渡圓弧半徑下，安裝傾角θ=15°時(shí)邊緣滾刀受力最小，安裝傾角θ=60°時(shí)邊緣滾刀受力最大。垂直力、滾動(dòng)力和扭矩的平均值及峰值平均值都隨著邊緣滾刀安裝傾角的增大而呈現(xiàn)近似線性的增長(zhǎng)。

3.2.2邊緣滾刀受力與刀盤過度段圓弧半徑的關(guān)系

圖10為邊緣滾刀在不同刀盤過渡圓弧半徑下的垂直力、滾動(dòng)力和扭矩的變化情況。在一定安裝傾角的情況下，過渡圓弧半徑r=200 mm時(shí)邊緣滾刀受力最小，r=800 mm時(shí)邊緣滾刀受力最大。垂直力、滾動(dòng)力和扭矩的平均值及峰值平均值都隨著刀盤過渡圓弧半徑的增大而呈現(xiàn)近似線性的增長(zhǎng)。

(a) 垂直力

(b) 滾動(dòng)力

(c) 扭矩

Fig. 9Vertical stresses, rolling stresses and torques of gage disc cutter under different installation angles

(a) 垂直力

(b) 滾動(dòng)力

(c) 扭矩

Fig. 10Vertical stresses, rolling stresses and torques of gage disc cutter under different transition arch radii

3.3綜合考慮安裝傾角及過渡圓弧半徑下的邊緣滾刀受力分析

從圖9可以看出，不同過渡圓弧半徑時(shí)的邊緣滾刀受力隨傾角的變化是不同的。r=800 mm時(shí)邊緣滾刀受力隨傾角增長(zhǎng)的速率大于r=200 mm時(shí)。過渡圓弧半徑越大，滾刀受力回歸直線的斜率越大，邊緣滾刀受力隨安裝傾角的增長(zhǎng)越迅速。從圖10可以看出，不同安裝傾角時(shí)的邊緣滾刀受力與過渡圓弧半徑的關(guān)系也是不同的。θ=60°時(shí)滾刀受力隨傾角增長(zhǎng)的速率要大于θ=15°時(shí)。邊緣滾刀安裝傾角越大，滾刀受力回歸直線的斜率越大，邊緣滾刀受力隨過渡圓弧半徑的增長(zhǎng)越迅速。在分析邊緣滾刀的受力情況時(shí)，除了應(yīng)分析單因素作用下的邊緣滾刀受力規(guī)律之外，還應(yīng)將安裝傾角和過渡圓弧半徑2種因素結(jié)合在一起綜合考慮。采用多元線性回歸分析可得：

(5)

(6)

(7)

(8)

(9)

(10)

4結(jié)論與討論

1)巖石具有一定的彈脆性，在邊緣滾刀切削巖石過程中，會(huì)出現(xiàn)階躍破碎的特點(diǎn)，滾刀受力曲線表現(xiàn)為高頻率振蕩，不斷出現(xiàn)陡增和回落。

2)安裝傾角和刀盤過渡圓弧半徑是影響邊緣滾刀受力的重要因素。邊緣滾刀受力平均值及峰值平均值都隨著安裝傾角和刀盤過渡圓弧半徑的增加而呈現(xiàn)近似線性的增長(zhǎng)。

3)邊緣滾刀的受力情況將直接影響邊緣滾刀的使用壽命。目前關(guān)于邊緣滾刀的布置參數(shù)對(duì)其受力情況影響的研究較少，這方面有待繼續(xù)深入研究。

參考文獻(xiàn)(References)：

[1]夏毅敏, 林賚貺, 羅德志, 等. 復(fù)合式土壓平衡盾構(gòu)盤形滾刀布置規(guī)律[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013,44(9): 3652-3657.(XIAYimin,LINLaikuang,LUODezhi,etal.DisccutterlayoutlawforcompositeEPBshield[J].JournalofCentralSouthUniversity(ScienceandTechnology), 2013,44(9): 3652-3657.(inChinese))

[2]沈斌.TBM邊緣滾刀組合優(yōu)化布置研究[D]. 長(zhǎng)沙： 中南大學(xué), 2014.(SHENBin.ResearchonoptimallayoutofTBMgagedisccutterscombination[D].Changsha:CentralSouthUniversity, 2014. (inChinese))

[3]宋克志, 李傳明, 袁大軍, 等. 盾構(gòu)隧道施工盤形滾刀推力分布規(guī)律研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2008,27(增刊2): 3875-3881. (SONGKezhi,LIChuanming,YUANDajun,etal.Studyondistributionregularitiesofdisccutterthrustforceforshieldtunneling[J].ChineseJournalofRockMechanicsandEngineering, 2008,27(S2): 3875-3881. (inChinese))

[4]張厚美.TBM的掘進(jìn)性能數(shù)值仿真研究[J]. 隧道建設(shè), 2006, 26(增刊2): 1-7. (ZHANGHoumei.Studyonnumericalsimulationsofperformanceoftunnelboringmachines(TBM)[J].TunnelConstruction, 2006, 26(S2): 1-7.(inChinese))

[5]吳元, 夏毅敏, 郭金成, 等.TBM邊緣滾刀關(guān)鍵參數(shù)對(duì)破巖效率影響規(guī)律研究[J]. 現(xiàn)代隧道技術(shù), 2015, 52(1): 119-126.(WUYuan,XIAYimim,GUOJincheng,etal.InfluenceofkeyparametersofaTBMgagedisccutteronrock-breakingefficiency[J].ModernTunnellingTechnology, 2015, 52(1): 119-126.(inChinese))

[6]王旭, 趙羽, 張寶剛, 等.TBM滾刀刀圈磨損機(jī)理研究[J]. 現(xiàn)代隧道技術(shù), 2010, 47(5): 15-19.(WANGXu,ZHAOYu,ZHANGBaogang,etal.ResearchontheringwearmechanismofTBMdisccutter[J].ModernTunnellingTechnology, 2010, 47(5): 15-19. (inChinese))

[7]楊宏欣. 淺談隧道施工中TBM滾刀失效及刀具管理[J]. 建筑機(jī)械化, 2008, 28(12): 60-62.(YANGHongxin.BrieflytalkingaboutTBMcuttermanagementinconstructingtunnel[J].ConstructionMechanization, 2008, 28(12): 60-62. (inChinese))

[8]高瑋．巖石力學(xué)[M]. 北京: 北京大學(xué)出版社，2010: 27-30.(GAOWei.Rockmechanics[M].Beijing:PekingUniversityPress, 2010: 27-30.(inChinese))

[9]朱震寰.TBM盤形滾刀磨損演變數(shù)值模擬研究[D]. 長(zhǎng)沙： 中南大學(xué), 2013.(ZHUZhenhuan.Numericalsimulationofdisccutter’swearevolutionforTBM[D].Changsha:CentralSouthUniversity, 2013. (inChinese))

[10]鄭聰, 趙威, 張海豐, 等. 不同貫入度下盤形滾刀受力數(shù)值研究[J]. 礦山機(jī)械, 2015,43(8): 118-123.(ZHENG Cong, ZHAO Wei, ZHANG Haifeng, et al. Numerical analysis of stress acting on disc cutter at various penetration depths[J]. Mining & Processing Equipment, 2015,43(8): 118-123.(in Chinese))

[11]宋克志, 袁大軍, 王夢(mèng)恕. 隧道掘進(jìn)機(jī) (TBM) 刀盤轉(zhuǎn)速的討論[J]. 建筑機(jī)械, 2005(8): 63-65.(SONG Kezhi, YUAN Dajun, WANG Mengshu. Discussion on cutterhead rotation speed of tunnel boring machine (TBM)[J]. Construction Machinery, 2005(8): 63-65.(in Chinese))

[12]紀(jì)昌明, 張照煌, 葉定海. 盤形滾刀刀間距對(duì)巖石躍進(jìn)破碎參數(shù)的影響[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào), 2008, 16(2): 255-263.(JI Changming, ZHANG Zhaohuang, YE Dinghai. The influence of the disk cutter space on rock’s jump break coefficients[J]. Journal of Basic Science and Engineering, 2008, 16(2): 255-263.(in Chinese))

[13]李輝, 王樹林, 汪加科. TBM 盤形滾刀受力分析[J]. 現(xiàn)代隧道技術(shù), 2012, 49(3): 193-197.(LI Hui, WANG Shulin, WANG Jiake. Force analysis of the disc cutter on a TBM[J]. Modern Tunnelling Technology, 2012, 49(3): 193-197.(in Chinese))

Influence of Layout Parameters on Stress of TBM Gage Disc Cutter

ZHENG Cong1， ZHAO Wei2， ZHANG Haifeng1， MA Baosong1

(1.FacultyofEngineering,ChinaUniversityofGeosciences(Wuhan),Wuhan430074,Hubei,China;2.CCCCSecondHighwayConsultantCo.,Ltd.,Wuhan430056,Hubei,China)

Abstract:Gage disc cutters installed on the edge of cutterhead of full-face tunnel boring machine (TBM) are subjected to complexity of stress due to the bad working conditions. Nonlinear dynamic models are established by using ABAQUS software so as to study the influence of layout parameters on stress of gage disc cutter; and the stress on the TBM gage disc cutter is analyzed. The results show that: 1) The stress on the gage disc cutter during rock breaking varies violently. 2) The installation angle and transition arc radius have a significantly influence on the stress of gage disc cutter. 3) The stress on the gage disc cutter increases approximately linearly with the installation angle and transition arc radius increase.

Keywords:TBM; gage disc cutter; installation angle; transition arc radius; finite elements

中圖分類號(hào)：U 455.3

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1672-741X(2016)03-0349-06

DOI:10.3973/j.issn.1672-741X.2016.03.017

作者簡(jiǎn)介：第一 鄭聰(1993—)，男，湖北荊州人，中國(guó)地質(zhì)大學(xué)(武漢)地質(zhì)工程專業(yè)在讀碩士，主要研究方向?yàn)榈刭|(zhì)工程。E-mail: zhengcongcug@163.com。

收稿日期：2015-09-14; 修回日期： 2015-11-08