一种fcec施工方法引起土体水平位移的计算方法
【专利摘要】本发明公开了一种FCEC施工方法引起土体水平位移的计算方法，采用FCEC全回转清障机，包括以下步骤：S1：根据土体的土质参数，计算得到所述FCEC全回转清障机的钢套筒周边土体的塑性区半径，所述土质参数包括：不排水抗剪强度，弹性模量，泊松比，容重；S2：根据步骤S1中的所述塑性区半径和土质参数，计算得到弹性区和塑性区交界面处的径向应力和径向位移；S3：根据步骤S1中的所述塑性区半径、步骤S2中的所述径向应力和径向位移，计算得到塑性区域的位移场和弹性区域的位移场。本发明综合考虑了在钢套筒进行旋转切割时对其周围土体造成环向剪切的影响，以此通过进一步计算来确定FCEC施工方法引起土体的水平位移。
【专利说明】
-种FCEC施工方法引起±体水平位移的计算方法
技术领域
[0001] 本发明设及粧基工程领域，特别是设及一种FCEC施工方法引起±体水平位移的计 算方法。
【背景技术】
[0002] FCECXF'ast Qit Evacuate Construction,全回转快速切割清障）施工方法是目前 工程领域中较先进的拔粧技术之一，该拔粧技术采用的设备是FCEC全回转清障机。该设备 通过电动机头驱动薄壁钢套筒对粧侧±体进行旋转切割，使粧与周围±体分离，当钢套筒 钻进深度超过粧长后，通过特殊装置将钢丝绳送到粧底，与粧身锁扣牢固后将旧粧进行拔 出。考虑整个施工过程，钢套筒对±体的旋转切割是FCEC施工方法在施工时引起±体水平 位移的主要来源，而FCEC施工方法引起的±体水平位移极易对钢套筒周围的已有地下结构 物产生破坏。如果在拔粧施工前，对±体的水平位移产生的扰动范围不进行预测，可能会对 已有地下结构物产生破坏，发生工程事故。在钢套筒旋转切割时，对其周围±体会造成环向 剪切的影响，因此，需要考虑FCEC施工方法在施工过程中的环向剪切作用，来确定FCEC施工 方法引起±体的水平位移，W获得对水平位移更加准确的计算。
[0003] 在现有技术中，虽然该FCEC施工方法在工程中已有较多的应用，但对其施工引起 ±体水平位移的研究还非常不足，仅处于初步阶段，理论分析与现场试验方面的资料极其 缺乏。经过对现有国内外技术文献的检索，尚未发现用来解决FCEC施工方法在施工中对钢 套筒旋转切割引起±体水平位移计算方法的文献。因此，业界需要一种对FCEC施工方法引 起±体水平位移进行计算的方法，应用在敏感环境下提供工程理论指导，避免工程事故的 发生。
【发明内容】
[0004] 本发明所要解决的技术问题是提供一种对FCEC施工方法引起±体水平位移进行 计算的方法。
[0005] 为了实现上述目的，本发明采用如下技术方案予W实现：
[0006] 一种FCEC施工方法引起±体水平位移的计算方法，采用FCEC全回转清障机，包括 W下步骤：
[0007] S1:根据±体的±质参数，计算得到所述FCEC全回转清障机的钢套筒周边±体的 塑性区半径，所述±质参数包括:不排水抗剪强度，弹性模量，泊松比，容重；
[000引S2:根据步骤S1中的所述塑性区半径和±质参数，计算得到弹性区和塑性区交界 面处的径向应力和径向位移；
[0009] S3:根据步骤S1中的所述塑性区半径、步骤S2中的所述径向应力和径向位移，计算 得到塑性区域的位移场和弹性区域的位移场。
[0010] 优选的，在步骤S1中，通过公式1计算得到所述塑性区半径：
[00川公式1:
[0012]
[001引式中：町为上体的塑性区半径，单位为m;
[0014] Ru为钢套筒的外径，单位为m;
[0015] Δ为塑性区平均体积应变；
[0016] B功中间变量，所述中间变量Bi通过公式2计算：
[0017] 公式 2:
[001 引
[0019] 式中:6为±体的剪切模量，单位为kPa;
[0020] 3。为±体的不排水抗剪强度，单位为kPa;
[0021 ] Δ为塑性区平均体积应变；
[0022] Ru为钢套筒的外径，单位为m;
[0023] Ro为钢套筒的内径，单位为m。
[0024] 优选的，通过公式3计算得到所述±体的剪切模量：
[002引公式3:
[0026] G = E/(2+2v)
[0027] 式中，6为±体的弹性模量，单位为kPa;
[0028] v为±体的泊松比。
[0029] 优选的，在步骤S2中，通过公式4计算得到所述弹性区和塑性区交界面处的径向应 力：
[0030] 公式 4:
[0031]
[0032] 式中：py为弹性区和塑性区交界面处的径向应力，单位为kPa;
[0033] Ru为钢套筒的外径，单位为m;
[0034] 3。为±体的不排水抗剪强度，单位为kPa;
[0035] 。为±体的塑性区半径，单位为m;
[0036] 口〇为±体中初始各向同性应力，单位为kPa。
[0037] 优选的，通过公式5计算得到±体中初始各向同性应力，
[003引公式5:
[0039] p〇= γ h
[0040] 式中，h为±体的深度，单位为m;
[0041] 丫为±体的容重，单位为N/m3。
[0042] 优选的，通过公式6计算得到所述弹性区和塑性区交界面处的径向位移：
[00创公式6:
[0044]
[0045] 式中，Up为弹性区和塑性区交界面处的径向位移，单位为m;
[0046] py为弹性区和塑性区交界面处的径向应力，单位为kPa;
[0047] 。为±体的塑性区半径，单位为m;
[004引6为±体的剪切模量，单位为kPa;
[0049] 口〇为±体中初始各向同性应力，单位为kPa。
[0050] 优选的，在步骤S3中，通过公式7计算得到所述弹性区域和塑性区域的位移场： [0化1] 公式7:
[0化2]
[0053] 式中，化为弹性区域和塑性区域的水平位移，单位为m;
[0054] Py为弹性区和塑性区交界面处的径向应力，单位为kPa;
[0055] Up为弹性区和塑性区交界面处的径向位移，单位为m;
[0056] 。为±体的塑性区半径，单位为m;
[0057] 6为±体的剪切模量，单位为kPa;
[005引口0为±体中初始各向同性应力，单位为kPa。
[0059] 优选的，在步骤S1中，通过原位十字板剪切试验，测定所述±质的不排水抗剪强 度。
[0060] 优选的，在步骤S1中，通过静力触探试验，测定所述±质的不排水抗剪强度。
[0061 ]优选的，通过公式8计算得到所述±质的不排水抗剪强度：
[0062] 公式 8:
[0063] su = 0.05Ps
[0064] 式中，Ps为静力触探试验得到的比贯入阻力。
[0065] 与现有技术相比，本发明的技术方案:根据±体的±质参数，计算得到所述FCEC全 回转清障机的钢套筒周边±体的塑性区半径;进一步计算得到弹性区和塑性区交界面处的 径向应力和径向位移;再通过所述塑性区半径、所述径向应力和所述径向位移，计算得到塑 性区域的位移场和弹性区域的位移场。本技术方案综合考虑了在钢套筒进行旋转切割时对 其周围±体造成环向剪切的影响，W此通过进一步计算来确定FCEC施工方法引起±体的水 平位移，提供了一种适用于FCEC施工方法所引起±体水平位移的计算方法，且该计算方法 更加准确、便捷。
【附图说明】
[0066] 图1是本发明所述FCEC施工方法的模型示意图；
[0067] 图2是本发明一实施例所述FCEC施工方法引起±体水平位移的曲线示意图。
[006引图中所示:ro、钢套筒内径;钢套筒外径;Pu、钢套筒外壁与±体接触界面上的压 力;町、塑性区半径;S、塑性区；T、弹性区。
【具体实施方式】
[0069] 下面结合附图对本发明作详细描述：
[0070] 参照图1所示，本发明的FCEC施工方法引起±体水平位移的计算方法，采用FCEC全 回转清障机，包括W下步骤：
[0071] S1:收集工作区已有地质资料，掌握场地±质情况，确定±体的±质参数及所述 FCEC全回转清障机的钢套筒内径r日和外径。，所述±质参数包括:不排水抗剪强度su，弹性 模量E，泊松比V，容重丫。所述FCEC全回转清障机的钢套筒内径ro表示初始的±体圆柱孔半 径、钢套筒外径。表示扩张后的±体圆柱孔半径，钢套筒外壁与±体接触界面上的压力为 Pu。在弹塑性分析中，钢套筒周围分为弹性区T和塑性区S，扩张后±体的塑性区S半径为 弹性区T和塑性区S边界面处的径向位移和径向应力分别为up、py，扩孔过程发生的塑性区S 平均体积应变为Δ。求解时，规定应力、应变W压为正。根据所述±质参数，计算得到所述 FCEC全回转清障机的钢套筒周边±体的塑性区S半径。
[0072] 作为一种实施方法，通过原位十字板剪切试验，测定所述±质的不排水抗剪强度 Su 〇
[0073] 采用上述技术，由于十字板剪切试验不需要采取±样，避免了±样扰动及天然应 力状态的改变，是一种有效的现场测试方法。
[0074] 作为另一种实施方法，通过静力触探试验，测定所述±质的不排水抗剪强度Su。在 所述静力触探试验中，通过公式：Su = 0.05Ps，计算得到所述±质的不排水抗剪强度Su，其 中，Ps为静力触探试验得到的比贯入阻力。
[0075] 采用上述技术，可W在现场直接得到各±层的贯入阻力指标，从而能够了解±层 在原始状态下的有关物理力学参数，具有连续、快速和准确的优点。
[0076] 进一步的，通过公式1计算得到所述塑性区S半径：
[0077] 公式 1:
[007引
[0079] 式中：町为±体的塑性区S半径，单位为m;
[0080] Ru为钢套筒的外径，单位为m;
[0081 ] Δ为塑性区S平均体积应变，可取为0.015;
[0082] Bi为中间变量；
[0083] 进一步的，所述中间变量Bi通过公式2计算：
[0084] 公式 2:
[0085]
[0086] 式中:6为±体的剪切模量，单位为kPa;
[0087] 3。为±体的不排水抗剪强度，单位为kPa;
[0088] Δ为塑性区S平均体积应变；
[0089] Ru为钢套筒的外径，单位为m;
[0090] Ro为钢套筒的内径，单位为m。
[0091 ]进一步的，通过公式3计算得到所述±体的剪切模量：
[0092] 公式 3:
[0093] G = E/(2+2v)
[0094] 其中，6为±体的弹性模量，单位为kPa;
[00巧]v为±体的泊松比，可取值为0.3。
[0096] 采用上述技术，通过±体的弹性模量和剪切模量的转换，通过采集±体的±质弹 性模量来获得其剪切模量，具有更易测量和实现的优点，提高了本发明的实用性。
[0097] S2:根据步骤S1中的所述塑性区S半径。和±质参数，计算得到弹性区T和塑性区S 交界面处的径向应力Py和径向位移Up。
[0098] 进一步的，通过公式4计算得到所述弹性区T和塑性区S交界面处的径向应力：
[0099] 公式 4:
[0100]
[0101] 式中：Py为弹性区T和塑性区S交界面处的径向应力，单位为kPa;
[0102] Ru为钢套筒的外径，单位为m;
[0103] 3。为±体的不排水抗剪强度，单位为kPa;
[0104]。为±体的塑性区S半径，单位为m;
[0105] 口〇为±体中初始各向同性应力，单位为kPa。
[0106] 进一步的，通过公式5计算得到±体中初始各向同性应力，
[0107] 公式 5:
[010引 po二丫 h
[0109] 式中，]1为±体的深度，单位为m;
[0110] 丫为±体的容重，单位为N/m3。
[0111] 采用上述技术，可W根据参数h的不同取值，来确定不同±体深度对应的初始各向 同性应力P0,进一步计算该深度h对应的弹性区T和塑性区S的位移场，能够实现空间计算的 效果。
[0112] 进一步的，通过公式6计算得到所述弹性区T和塑性区S交界面处的径向位移：
[0113] 公式 6:
[0114]
[0115] 式中，Up为弹性区T和塑性区S交界面处的径向位移，单位为m;
[0116] Py为弹性区T和塑性区S交界面处的径向应力，单位为kPa;
[0117] 。为±体的塑性区S半径，单位为m;
[0118] 6为±体的剪切模量，单位为kPa;
[0119] 口〇为±体中初始各向同性应力，单位为kPa。
[0120] S3:根据步骤S1中的所述塑性区S半径步骤S2中的所述径向应力Py和径向位移 Up,计算得到塑性区S的位移场和弹性区T的位移场。
[0121] 进一步的，通过公式7计算得到所述弹性区Τ和塑性区S的位移场：
[0122] 公式 7:
[0123]
[0124] 式中，化为弹性区T和塑性区S的水平位移，单位为m;
[0125] py为弹性区T和塑性区S交界面处的径向应力，单位为kPa;
[0126] Up为弹性区T和塑性区S交界面处的径向位移，单位为m;
[0127] 。为±体的塑性区S半径，单位为m;
[0128] 6为±体的剪切模量，单位为kPa;
[0129] 口〇为±体中初始各向同性应力，单位为kPa。
[0130] 本发明的技术方案根据±体的±质参数，计算得到所述FCEC全回转清障机的钢套 筒周边±体的塑性区S半径进一步计算得到弹性区T和塑性区S交界面处的径向应力Py和 径向位移Up;再通过所述塑性区S半径所述径向应力Py和所述径向位移Up,计算得到塑性 区S的位移场和弹性区T的位移场。
[0131] 在实施例中，假设各向同性的均质饱和粘±受到初始均匀同性应力P0的作用，p〇 = lOOkPa，初始剪应力为0，±体的弹性模量E = 25MPa，泊松比v = 0.3，不排水抗剪强度su = 25MPa，假设塑性区S±体的平均体积应变Δ=0.005。钢套筒内径为1.5m，外径为1.54m。由 公式1~7,可最后计算得到塑性区S的位移场和弹性区T的位移场，其计算结果曲线图通过 Matlab软件显示如图2所示，由图2可W判断出钢套筒旋入±体中的位移影响范围。
[0132] 本技术方案综合考虑了在钢套筒进行旋转切割时对其周围±体造成环向剪切的 影响，W此通过进一步计算来确定FCEC施工方法引起±体的水平位移，本技术方案提供了 一种适用于FCEC施工方法所引起±体水平位移的计算方法，且该计算方法更加准确、便捷。 本技术方案为FCEC拔粧施工技术在敏感环境下提供工程理论指导，避免工程事故的发生。
【主权项】
1. 一种FCEC施工方法引起土体水平位移的计算方法，采用FCEC全回转清障机，其特征 在于，包括以下步骤： S1:根据土体的土质参数，计算得到所述FCEC全回转清障机的钢套筒周边土体的塑性 区半径，所述土质参数包括:不排水抗剪强度，弹性模量，泊松比，容重； S2:根据步骤S1中的所述塑性区半径和土质参数，计算得到弹性区和塑性区交界面处 的径向应力和径向位移； S3:根据步骤S1中的所述塑性区半径、步骤S2中的所述径向应力和径向位移，计算得到 塑性区域的位移场和弹性区域的位移场。2. 根据权利要求1所述的计算方法，其特征在于，在步骤S1中，通过公式1计算得到所述 塑性区半径： 公式1:式中:rP为土体的塑性区半径，单位为m; Ru为钢套筒的外径，单位为m; Δ为塑性区平均体积应变； Βι为中间变量，所述中间变量Βι通过公式2计算： 公式2:式中:G为土体的剪切模量，单位为kPa; su为土体的不排水抗剪强度，单位为kPa; Δ为塑性区平均体积应变； Ru为钢套筒的外径，单位为m; Ro为钢套筒的内径，单位为m。3. 根据权利要求1所述的计算方法，其特征在于，通过公式3计算得到所述土体的剪切 模量： 公式3: G = E/(2+2v) 式中，E为土体的弹性模量，单位为kPa; v为土体的泊松比。4. 根据权利要求1所述的计算方法，其特征在于，在步骤S2中，通过公式4计算得到所述 弹性区和塑性区交界面处的径向应力： 公式4:式中:py为弹性区和塑性区交界面处的径向应力，单位为kPa; Ru为钢套筒的外径，单位为m; su为土体的不排水抗剪强度，单位为kPa; rP为土体的塑性区半径，单位为m; po为土体中初始各向同性应力，单位为kPa。5. 根据权利要求1所述的计算方法，其特征在于，通过公式5计算得到土体中初始各向 同性应力， 公式5: ρ〇= γ h 式中，h为土体的深度，单位为m; γ为土体的容重，单位为N/m3。6. 根据权利要求1所述的计算方法，其特征在于，在步骤S2中，通过公式6计算得到所述 弹性区和塑性区交界面处的径向位移： 公式6:式中，uP为弹性区和塑性区交界面处的径向位移，单位为m; py为弹性区和塑性区交界面处的径向应力，单位为kPa; rP为土体的塑性区半径，单位为m; G为土体的剪切模量，单位为kPa; po为土体中初始各向同性应力，单位为kPa。7. 根据权利要求1所述的计算方法，其特征在于，在步骤S3中，通过公式7计算得到所述 弹性区域和塑性区域的位移场： 公式7:式中，ur为弹性区域和塑性区域的水平位移，单位为m; py为弹性区和塑性区交界面处的径向应力，单位为kPa; uP为弹性区和塑性区交界面处的径向位移，单位为m; rP为土体的塑性区半径，单位为m; G为土体的剪切模量，单位为kPa; po为土体中初始各向同性应力，单位为kPa。8. 根据权利要求1所述的计算方法，其特征在于，在步骤S1中，通过原位十字板剪切试 验，测定所述土质的不排水抗剪强度。9. 根据权利要求1所述的计算方法，其特征在于，在步骤S1中，通过静力触探试验，测定 所述土质的不排水抗剪强度。10. 根据权利要求9所述的计算方法，其特征在于，通过公式8计算得到所述土质的不排 水抗剪强度： 公式8: su = 0.05Ps 式中，PS为静力触探试验得到的比贯入阻力。
【文档编号】G06F17/50GK105825012SQ201610150863
【公开日】2016年8月3日
【申请日】2016年3月16日
【发明人】孟振, 陈锦剑, 王建华
【申请人】上海交通大学