螺旋钻头掘削土壤的阻力矩研究.pdf
螺旋钻头掘削土壤的阻力矩研究 ① 胡均平, 张政华, 郭 勇, 李亮红 (中南大学 机电工程学院, 湖南 长沙 410083) 摘 要 考虑进给和回转运动共同对切削厚度和切削角的影响,引用进给回转速比 m,结合 McKyes⁃Ali 土壤三维失效模型和极限 力学平衡方程,建立了新的螺旋切削阻力矩解析公式。 理论和实验研究表明,螺旋切削阻力矩与进给回转速比 m 成非线性关系;m 的存在使实际切削角、实际切削厚度及实际失效角都发生改变,导致钻头切削阻力矩理论值与不考虑 m 影响时不同;两理论值偏差 随 m 没有确定的变化规律,偏差的大小与土壤力学参数、钻具结构参数和上层土预压力等有关;实验值与考虑进给回转速比影响的 理论值有更高的吻合度,证明了理论的正确性。 关键词 螺旋钻头; 进给回转速比; 切削阻力矩 中图分类号 U415.5文献标识码 Adoi10.3969/ j.issn.0253-6099.2014.02.004 文章编号 0253-6099(2014)02-0015-05 Research on Resistance Torque of a Working Auger Head HU Jun⁃ping, ZHANG Zheng⁃hua, GUO Yong, LI Liang⁃hong (School of Mechanical and Electrical Engineering, Central South University, Changsha 410083, Hunan, China) Abstract In view of the impacts of feeding and rotary motion on real cutting thickness and cutting angle, the parameter m, that is the feeding⁃rotational speed ratio, was introduced in the study. In combination with McKyes⁃Ali′s three⁃ dimensional soil failure pattern and the corresponding limited equilibrium mechanics equations, a new analytical formula for cutting resistance torque of a working auger head was deduced. Theoretical and experimental researches show that, the relationship between feeding⁃rotational speed ratio and cutting resistance torque is nonlinear. The introduction of the parameter m leads to changes in the real cutting angle, real cutting thickness and real failure angle, resulting in definitely different cutting resistance torque. The difference between two theoretical values doesn′t follow any regular pattern as m value changes, but it is connected with mechanical parameters of soil, structural parameters of drill and pre⁃loaded pressure of topsoil. The experimental data being in a good agreement with theoretical value, has further testified the feasibility of such theory. Key words auger head; feeding⁃rotational speed ratio; cutting resistance moment 螺旋掘削作为基础设施建设、采矿工程和农林业 机械等领域的成熟技术具有广泛的运用基础,其基本 原理即采用电动或液压顶驱式动力头为钻具提供扭矩 切削岩土。 螺旋掘削力学行为属于机具岩土交互作用 研究内容,其基础理论可追溯到 20 世纪中期的平面切 削阻力理论。 1965 年 Reece 首次提出了二维土壤失 效方式的阻抗力预测公式,但其中众多未知的 N 影响 系数,使公式推广受阻[1]。 1984 年 McKyes 在 Godwin⁃ Spoor 窄齿切削土壤三维失效模型基础上对中心失效 域做出楔形块简化假设,运用连续介质力学得到不计 土壤与刀齿粘滞力的阻抗力公式[2]。 1995 年 Ji Zhang 和 Kushwaha 重新考虑该粘滞力并简化 McKyes 公 式[3],被后续研究者广泛运用到挖掘机挖斗力学建 模[4]和农机耕具设计[5-6]等领域。 2006 年 Wheeler 和 Godwin 经实验证实了切削速度引起惯性力对阻抗力 的影响,并指出惯性力产生影响的临界值[7]。 国内对 该领域内的研究较晚, 1998 年开始赵伟民对土的三维 掘削形态做了大量研究,并提出了引入切削比阻力的 回转斗掘削能耗半经验公式[8-9]。 2004 年开始胡长胜 基于这些公式和 SPT 图研究了钻机钻进参数智能匹 ①收稿日期 2013-10-30 基金项目 湖南省科技计划项目(2011WK3041) 作者简介 胡均平(1965-),男,湖南邵阳人,教授,博士,主要研究方向为液压机械与控制。 作者简介 张政华(1987-),男,湖南益阳人,硕士,主要研究方向为液压机械。 第 34 卷第 2 期 2014 年 04 月 矿矿 冶冶 工工 程程 MINING AND METALLURGICAL ENGINEERING Vol.34 №2 April 2014 配系统设计[10-11]。 2007 年,郭峰结合试验首次提出旋 挖切削土阻力的解析公式[12]。 2009 年暨智勇运用同 样的模型对盾构掘进机切刀切削软岩和土壤的受力解 析解进行了理论和试验研究[13]。 这些研究对确定掘 进设备的受力有重要意义,但目前对螺旋掘削阻力矩 的解析表达式均假设为刀齿在指定厚度和切削角下的 二维回转切削情况,忽略了螺旋掘削是三维掘削的事 实,即进给和回转运动共同影响实际切削厚度和切 削角。 本文依托 McKyes⁃Ali 土壤三维失效模型,重新推 导了土壤三维螺旋掘削的稳态切削阻力矩公式,并分 析了影响螺旋掘削阻力矩的变量。 1 实际切削参数 螺旋掘削过程中钻具进给速度 v(m/ s)与回转速 度 n(r/ min)之比 m(简称进给回转速比)不仅影响钻 进设备的施工工艺,并且是影响刀齿实际切削角和切 削厚度的关键变量。 例如地锚钻机和螺纹桩机施工时 需严格控制钻进速度和回转速度使 60m=p(钻头叶片 螺距),保证钻头每回转一周就进给一个螺距[14]。 长 螺旋钻机、套管钻机和旋挖钻机等施工时 60m<p,刀具 对土壤起切削作用,土屑堆积在螺旋叶片或因孔壁摩 擦和下层土屑挤压作用沿叶片上升[15]。 理论设计时 不存在 60m>p 的情况。 考虑钻头以特定的 m 切削土壤,取一个回转周期 内钻头刀齿的工作轨迹并展开,如图 1 所示,刀齿可简 化为工作于与水平面成角度 θ 的斜面,且沿斜面方向 实际切削土壤厚度为 d,其中 θ = arctan 60v nπD = arctan 60m πD d = 60v in cos(θ) = 1 i 1 πD 2 + 1 60m 2 { (1) 式中 v 为钻头进给速度,m/ s;n 为钻头回转速度, r/ min;i 为钻头头数,通常为单头或双头;D 为刀齿回 转直径,m。 由式(1)可知 θ 和切削厚度 d 均为 m 的单 调递增函数。 图 1 刀齿螺旋运动轨迹及展开图 2 力学模型与理论分析 参考 McKyes⁃Ali 土壤三维失效模型,将单个刀齿 切削土壤的失效域分为中心失效域及两侧对称的侧向 失效域,如图 2 所示。 图 2 土壤的三维失效模型 假设① 刀齿螺旋切削简化为斜面切削; ② 中心 失效域简化为一宽度与刀齿等宽的楔形块,失效面 BCGH 与面 GCEF 构成失效角 β;③ 侧向失效域轮廓 线 AB 简化为半径为 r 的圆弧,对应圆心角为 ρ′,A 到 中心失效域距离为 s。 考虑的力学因素包括土壤重力、上层覆盖土的预 压力、土壤内部阻抗力、土壤与刀齿的外部阻抗力。 因 岩土掘削设备施工时刀齿线速度多未达到文献[7]所 述的临界速度,故不考虑土壤惯性力。 将中心失效域的各作用力等效于偏转 θ 后的水平 切削情况,如图 3 所示,α′为刀齿安装切削角,此时的 实际切削角为 α = α′ - θ = α′ - arctan 60m πD (2) 图 3 中心失效域受力分析 61矿 冶 工 程第 34 卷 由中心失效域极限力学平衡方程 ∑Fx = 0 ∑Fy = 0 { (3) 即 P1sin(α + δ) + cawdcotα + qwrsinθ + γwdrsinθ 2 = R 1sin(β + φ) + cwdcotβ P1cos(α + δ) + cawd + R1cos(β + φ) = qwrcosθ + γgwdr 2 cosθ + cwd (4) 得刀齿对中心失效域的作用力为 P1= w γgdr 2 + qr N1+ cdN2 - c adN3 N4 (5) 其中 N1= cosθ - sinθcot(β + φ)(6) N2= 1 + cotβcot(β + φ)(7) N3= 1 - cotαcot(β + φ)(8) N4= cos(α + δ) + sin(α + δ)cot(β + φ)(9) 对侧向失效域取图 2(c)所示角度 ρ 处的微元,作 与中心失效域同样的处理,由极限力学平衡方程 ∑Fx = 0 ∑Fy = 0 { (10) 即 dp2sin(α + δ) + γgdr2dρsinθ 6 + qr2dρsinθ 2 = dR2sin(β + φ) + crddρcotβ 2 dp2cos(α + δ) + dR2cos(β + φ) = γgdr2dρcosθ 6 + crddρ 2 + qr2dρcosθ 2 (11) 得 dP2= 1 6 γgdr2+ 1 2 qr2dρ N1+ crddρ 2 N2 N4 (12) 刀齿对单个侧向失效域对的作用力为 P2=∫ ρ′ 0 dP2cosρ = sinρ′ 1 6 γgdr2+ 1 2 qr2 N1+ crd 2 N2 N4 (13) 由失效模型的几何关系得 s = dcot2β + 2cotαcotβ(14) r = d(cotα + cotβ)(15) sinρ′ = s r (16) 总阻抗力为 P = P1+ 2P2 = (γgd2Nγ+ qdNq+ cdNc + c adNca)w (17) 其中 Nγ= 1 2 N1N5(1 + 2d 3wN7) (18) Nq = N 1N5(1 + d w N7)(19) Nc = N 6(1 + d w N7)(20) Nca= cos(α + β + φ) sinαsin(α + δ + β + φ) (21) N5= (cotα + cotβ)sin(β + φ) sin(α + δ + β + φ) (22) N6= cosφ sinβsin(α + δ + β + φ) (23) N7=cot2β + 2cotαcotβ(24) 上述各式中 α 为实际切削角,();β 为失效角,(); δ 为土与刀齿摩擦角,();φ 为土壤内部摩擦角,();γ 为土壤重度,kg/ m3;c 为土内粘聚力,kPa;ca为土与刀 齿粘聚力,kPa;w 为刀齿宽度,m;g 为重力加速度,m/ s2。 上述各式求解的关键为失效角 β,根据被动土压 理论的最小阻力原则,β 为 Nγ取最小值时对应的角 度,即求解∂Nγ/ ∂β=0,结果可由 MATLAB 求根函数求 取。 β 求出后可得单个刀齿螺切削的水平阻抗力为 PH= Psin(α + θ + δ) + cawdcos(α + θ) / sinα (25) 一般螺旋掘削钻具如 AGS 型双头螺旋钻头、BSS 型旋挖钻斗[15]和短螺旋钻头等都安装有多个刀齿,其 切削扭矩为多个刀齿螺旋切削阻力矩之和 T = ∑ ni i = 1 PHi(Ri+ wi 2 )(26) 式中 ni为刀齿数量;Ri表示任意一个刀齿内侧到轴心 的距离,m;wi为刀齿宽,m。 由上述各式可知受刀齿螺旋运动的影响,土重系 数 Nγ和预压力系数 Nq相对 McKyes⁃Ali 平面切削阻 力公式多出[cosθ-sinθcot(β+φ)]项,平面二维切削实 际为取 θ=0 的特殊情况。 新项的引入同实际切削厚 度和实际切削角势必引起失效角 β 的变化。 表 1 是不 同 m 下某刀齿的实际切削厚度 d,实际切削角 α 和实 际失效角 β 与不考虑进给回转速比引起偏斜角 θ 影响 (即 θ=0)的对应理论比较值。 设刀齿在双头螺旋钻 头直径 595 mm 处,安装切削角 60。 71第 2 期胡均平等 螺旋钻头掘削土壤的阻力矩研究 表 1 理论比较值 m 值 α / () α′ / () d / mm d′ / mm β / () β′ / () 0.7558.66022.522.530.129.7 1.257.860363631.831 1.6557.06049.449.533.031.8 2.156.16062.96333.932.4 2.5555.36076.276.534.732.8 3.054.56089.69035.433.2 由表 1 可知实际切削角 α 较 θ = 0 时偏小,且 α 是 m 的单调递减函数;实际切削厚度受 θ 影响小,可 认为 d=d′,且 d 是 m 的单调递增函数;实际失效角 β 较 θ=0 时偏大,β 随 m 增大而增大。 3 螺旋切土试验 试验在湖南长河机械有限公司液压螺旋钻进试验 台上完成。 系统功率 3 kW,采用变量泵⁃定量马达的 容积变量调速方式,调节输出转速为 5~20 r/ min,用 E2E 型接近开关感应安装在钻具连接法兰处的磁性物 体,通过计算脉冲个数换算输出转速;采用工作压力 16 MPa、精度0.5%FS 的压力变送器测量动力头马达进油 腔压力并换算为扭矩。 由进给油缸控制进给运动以保 证速度,比例阀调节进油量,该实验控制钻具进给速度 15 mm/ s,以得到不同进给回转速比 m,采用量程 0~45 L/ min、精度 0.5%FS 的流量计检测流量并转换为进给 速度。 为保证模型的有效性,应控制 m 使 θ <α′-ε (ε 为刀齿刃角),否则刀齿背面承受土壤阻抗力,与建 模不符,为此在本试验中限定钻具转数至少为 5 r/ min。 试验钻头为双头结构,如图 4 所示,无输土螺 旋叶片,消除叶片输土对扭矩的影响。 6 个刀齿中心 对称布置,安装切削角可调,分 30、45、60三种,保证 不同直径处刀齿有相同安装切削角,单个刀齿宽 50 mm,钻孔直径 700 mm(钻具设计直径一般比钻孔直径 小 10~50 mm[15])。 图 4 钻头结构示意 筛选粒度直径 5 mm 的土样,均匀搅拌,分层倒入 实验槽内喷水并夯实平整。 测量给定 m 下钻进土样深 度约为100 mm 时的进油腔压力值,选取代表性土样,通 过实验测得该深度范围附近土样力学参数如表 2 所示。 表 2 土样力学参数 含水量 / % 内摩擦角 / () 外摩擦角 / () 内粘聚力 / kPa 黏着力 / kPa 土重度 / (kNm -3 ) 45.237.327.38.53.215.3 针对每组试验选三次较理想的数据取平均值,测 得切削扭矩如图 5 所示,图中实线为考虑进给回转速 比影响的理论值 1,虚线为不考虑该影响的理论值 2, 黑点为测量值。 由图可知,安装切削角相同时,随着回 转速度的减小(即 m 的增大),两理论切削阻力均增大, 本实验中在同样的进给速度下,回转速度大于 10 r/ min 时会有较小的切削扭矩,表明较大的回转速度有利于 图 5 不同进给回转速比 m 下的切削扭矩 (a) α′=30; (b) α′=45; (c) α′=60 81矿 冶 工 程第 34 卷 切削;随着安装切削角的增大,在同样回转速度(即同 样 m)下,两理论切削扭矩均增大,这与文献[12]的结 论一致;两理论值偏差与 m 没有确定的变化规律,当 α′= 30时,理论值 1 较理论值 2 偏大,最大偏差为 13 3%;当 α′=45和 α′= 60时,理论值 1 较理论值 2 均偏小,最大偏差分别为-11.1%和-23.5%。 由实验 值分布可知,考虑进尺回转速比 m 引起倾角 θ 的理论 值与实验数据有更高的吻合度。 4 结 语 以 McKyes⁃Ali 土壤失效模型为原型研究了进给 速度和回转速度共同对螺旋掘削时实际切削厚度、实 际切削角和实际失效角的影响,推导出了单齿螺旋切 削阻力公式,进而得到了钻头在给定进给回转速比时 的稳态切削扭矩,通过理论及试验研究发现 1) 螺旋切削阻力矩与进给回转速比 m 成非线性 关系,随着 m 增大,切削阻力矩增大。 2) m 的存在使实际切削角、实际切削厚度及实际 失效角较不考虑 m 影响时不同,导致切削阻力理论值 有偏差。 3) 两理论值偏差随 m 没有确定的变化关系,偏 差的大小与土壤力学参数、钻具结构参数及上层土预 压力等相关。 参考文献 [1] Hettiaratchi D R P,Witney B D,Reece A R. The calculation of pas⁃ sive pressure in two⁃dimensional soil failure[J]. Geotechnique,1966, 11(2)89-107. [2] McKyes E,Desir F L. Prediction and field measurements of tillage tool draft forces and efficiency in cohesive soils[J]. Soil Tillage Research, 1984,4(5)459-470. [3] Zhang J,Kushwaha R L. A modified model to predict soil cutting re⁃ sistance[J]. Soil and Tillage Research,1995,34157-168. [4] Bhaveshkumar P Patel. Soil⁃tool interaction as a review for digging op⁃ eration of mini hydraulic excavator[J]. International Journal of Engi⁃ neering Science and Technology,2011,3(2)894-901. [5] Liu J,Lobb D A,Chen Y,et al. Steady⁃state models for the movement of soil and straw during tillage with a single sweep[J]. Transactions of the ASABE,2008,51(3)781-789. [6] Manuwa S I. Performance evaluation of tillage tines operating under different depths in a sandy clay loam soil[J]. Soil and Tillage Re⁃ search,2009,103(2)309-405. [7] Godwin R J. A review of the effect of implement geometry on soil fail⁃ ure and implement forces[J]. Soil and Tillage Research,2007,97 (2)331-340. [8] 赵伟民,顾迪民. 螺旋钻具上土的动力分析[J]. 哈尔滨建筑大学 学报,1998(3)90-93. [9] 赵伟民,陆念力,施引蕃,等. 回转斗掘削岩土的功率消耗计算 [J]. 建筑机械,2000(10)41-44. [10] 胡长胜,赵伟民,李瑰贤,等. 基于数据融合的钻孔机保护与钻具 动力智能匹配系统设计[J]. 中国工程机械学报,2004,2(1)75 -79. [11] 胡长胜,赵伟民,李瑰贤,等. SPT 图在钻孔机扭矩⁃转速控制系 统中的应用[J]. 哈尔滨工业大学学报,2005,37(11)1564- 1566. [12] 郭 峰,赵伟民,李瑰贤. 旋挖切削阻力的理论分析与实验研究 [J]. 岩土力学,2007,28(12)2666-2670. [13] 暨智勇. 盾构掘进机切刀切削软岩和土攘受力模型研究及实验 验证[D]. 长沙中南大学机电工程学院,2009. [14] 李波扬,吴 敏. 一种新型的全螺旋灌注桩⁃螺纹桩[J]. 建筑结 构,2004,34(8)55-57. [15] 刘家荣. 复杂地层桩孔钻进工艺及机具研究[D]. 北京中国地 质大学工程技术学院,2010. [16] 刘三意. 多工艺旋挖钻进技术研究[D]. 北京中国地质大学工 程技术学院,2008. 中国知网版权声明 本刊已许可中国学术期刊(光盘版)电子杂志社在中国知网及其系列数据库产品中以数字化方式 复制、汇编、发行、信息网络传播本刊全文。 该社著作权使用费与本刊稿酬一并支付。 作者向本刊提交 文章发表的行为即视为同意本刊上述声明。 矿冶工程编辑部 2014 年 4 月 91第 2 期胡均平等 螺旋钻头掘削土壤的阻力矩研究