Optimization of anchorage location of three dimensional slopes based on improved differential evolution algorithm
-
摘要: 基于改进的差分进化算法,在三维锚固力计算基础上考虑锚固力的三维空间分布状态,对三维边坡锚固位置进行合理优化,从而有效节约锚杆数量。利用集中力作用下的半无限弹性空间内一点应力解答,将锚固力等效叠加到边坡滑面正应力和剪应力上,由此推导了关于锚固位置坐标的锚固力函数,并以此作为锚固位置优化的目标函数。对于标准差分进化算法(DE)的改进,将自适应差分进化算法与一般约束乘子算法(PHR)结合,使算法能够解决锚固边坡约束优化问题。算例和工程实例表明,在边坡加固安全系数相同且单根锚杆受力保持不变的情况下,通过本文方法优化锚固位置能节约20%~25%的锚固点数量,该方法可为工程应用提供参考。Abstract: Based on the improved differential evolution algorithm, considering the 3D spatial distribution of anchoring force through calculation of the 3D anchoring force, the anchorage position of 3D slope is optimized reasonably to save the number of bolts. By using the stress solution of one point in the semi-infinite elastic space under the action of concentrated force, the anchoring force is superimposed on the normal stress and shear stress of the slope surface. Then the function for anchoring force about the coordinates of the anchorage position is derived, and taken as the objective function of the optimization of the anchorage position. For the improvement of the standard differential evolution algorithm (DE), the adaptive differential evolution algorithm is combined with the general constraint multiplier algorithm (PHR) to solve the problem of constraint optimization of anchoring slopes. Numerical and engineering examples show that with the same safety factor of slope reinforcement and unchanged single anchoring force, through the optimization of anchorage position, the number of anchorage points can be saved by 20% ~ 25%, and the proposed method can provide reference for engineering application.
-
Keywords:
- slope /
- limit equilibrium /
- 3D stability /
- safety factor
-
0. 引言
砾性土是对砾质土、砾砂、粉砾土、砂砾土和砾类土的统称[1],是含砾量从大于10%到100%的宽级配粗粒土,是工程建设的重要材料,地理分布广泛。砾土在自然和人工场地都曾发生地震液化或失稳破坏等现象,造成极大经济损失,如1975年唐山地震、2008年汶川地震与2016年新西兰凯库拉地震[2-3]。其力学特性与本构研究受人重点关注,但研究需基于单元试验。但橡皮膜作为单元试验关键影响因素,相关研究不甚充分,目前研究无法考虑并消除其影响,结果可靠性不足,严重制约粗粒土动力研究发展。因此,橡皮膜影响消除是粗粒土力学特性研究中的关键课题。
三轴不排水试验中,橡皮膜主要存在顺变效应[4],即不排水试验中,有效应力变化导致橡皮膜嵌入体积顺势改变而影响孔隙水压力真实表现的现象。砾性土液化试验中,随孔压增长,有效应力降低,橡皮膜从初始嵌入位置逐渐回弹,部分孔隙水运移填补橡皮膜回弹空间,导致试样局部排水,而影响土孔压发展和动强度。已有研究表明,橡皮膜对试验结果的影响是全方位的:最大可能导致85%体积测量误差[5-6];可能导致低估最多70%动孔隙水压力[7-8]以及最高达50%的抗液化强度误差[4,8-11];最高达50%静强度误差[12];对土的阻尼比和泊松比的测量也有明显影响[13]。
液化试验橡皮膜顺变性校正工作,分两个阶段:①第一阶段,橡皮膜嵌入体积的确定,是后续工作阶段的基础与前提,提供关键应用数据。②第二阶段,橡皮膜顺变性的消除,需基于前阶段结果,采取科学合理的方法实现对橡皮膜顺变性的校正。
其中,橡皮膜顺变性消除方法主要有物理缓解、仪器补偿和计算修正3类。物理缓解主要为处理试样外表面凹陷,如外表面涂抹硅橡胶[6]或黏土[14]、细砂充填孔隙[4,9]或冷冻试样表面光滑切割[15]等。这类方法对试样带来的附加影响难以估计,且如涂胶法无法100%抵消橡皮膜顺变[6]。仪器补偿需根据提前确定的橡皮膜嵌入体积,向试样内[11, 16-17]或特制的双层膜间[18]进行及时补水,补偿因橡皮膜顺变导致试样局部排水影响,被认为可100%消除橡皮膜影响[19],但大多需要特殊开发新设备,适用性较低,难作为常规手段推广。计算修正基于三轴试验,从孔压发展角度对结果修正[8, 20],得到相对真实的孔压发展规律,方法需根据大量复杂的工作建立合理孔压模型并获取可靠修正系数,在土动力学早期难以广泛开展,至今仍鲜有应用。综上,因应用效果、适用范围或理论发展等因素,橡皮膜顺变性消除方法在过去近20 a应用甚少,未能取得更深入发展,阻碍着液化理论深化,导致当前的土壤液化研究,尤其是大颗粒粗粒土的液化单元试验中,对其常避而不谈或选择性忽略。
本文基于橡皮膜顺变效应影响机理,以仪器补偿作为消除橡皮膜顺变性研究的主要手段,分析现有补偿手段存在的理论缺陷,提出一种对饱和试样进行补偿的新原理和新方法,结合砾性土的橡皮膜嵌入体积测量的前序成果,利用大型动三轴仪,以汶川地震中液化砾性土为基础料,开展新型补偿试验,通过与已有补偿方法、计算修正后无橡皮膜顺变性的理论结果等对比,验证新方法的技术可行性,理论可靠性和方法适用性。研究成果为开展粗粒动力试验研究提供新途径,为深化土动力学理论发展提供新动力。
1. 橡皮膜顺变性现有补偿理论方法
1.1 修正理论与计算方法
Martin等[20]提出了基于孔压产生及发展规律的针对橡皮膜顺变效应的理论和计算方法。理想不排水条件下,饱和土受动荷载所产生的孔隙水压力满足
Δu0=Er⋅Δεvd, (1) 式中,
Δu0 为完全不排水条件下的残余孔压增量,Er 为某动力初始阶段土的体积回弹模量,Δεvd 为排水条件下土在同等应力条件下受到相同动荷载作用后所产生的体应变增量。当进一步考虑了橡皮膜顺变性所导致土局部排水后,由式(1)得到
Δum=Δεvd1Er+1Em=Er⋅Δεvd(1+Cr), (2) 式中,
Δum 为存在橡皮膜顺变效应影响的残余孔压增量,Em 为某动力初始阶段时对应的橡皮膜的变形模量,Cr 为橡皮膜顺变性修正系数,Cr=Er/Em ,可由橡皮膜嵌入体积测量试验确定。结合式(1)可得Δu0=Δum⋅(1+Cr)。 (3) 式(3)即为考虑橡皮膜顺变性的孔隙水压力修正计算基础公式。其理论科学合理,被视为橡皮膜顺变性校正的理论结果,曾得到短暂发展和应用。但近20 a在粗粒土试验中鲜有深入研究,一方面,粗粒土孔压模型发展有限,如Finn模型的代表性参数由二氧化硅砂所得[21],能否推广至更宽级配粗粒土尚未可知。另外该理论对孔压全量模型并不适用,曾有学者误对累计残余孔压时程修正橡皮膜顺变效应影响[22],与孔压增长规律相悖。同时,确定修正系数
Cr 依赖于能够准确确定橡皮膜嵌入体积[4,23-24],这在以往难度较大。因此计算修正方法原理科学可靠但实际应用较少。1.2 现有仪器补偿方法及缺陷
Ramana等[17]提出一种利用压力控制器的补偿方法,简称Ramana方法(下同),成为后续学者研发更高级补水仪器和系统的核心理论和技术源头[11, 16],被认为可达100%补偿效果[19]。其无需进行常规液化试验,可直接由补偿液化试验得到土的真实抗液化强度。基本原理是将液化试验过程中,某动力阶段用于回弹橡皮膜的水再补回试样中,来保证土试样本身的不排水状态,可表达为
Vcp=Vd, (4) 式中,
Vcp 为向试样补偿水的体积,Vd 为因产生有效应力变化Δσ′3 而支持橡皮膜回弹的水的体积,由橡皮膜嵌入体积与有效应力关系确定,如图1中黑色实线所示。操作方法为,获取膜嵌入体积与有效应力关系,如图1黑色实线,随后由某动力阶段有效应力变量
Δσ′3 ,读取其导致的橡皮膜回弹体积ΔVm1 ,将该体积的水及时补充到试样内,可得无橡皮膜顺变影响的试验结果。考虑注水引起新的孔压增长导致橡皮膜顺变,如图1箭头与虚线所示需多回合注水,反复逼近使总注水体积与有效应力增量重新在曲线上重合,一次补偿工作完成,并可以开始下一阶段动力试验。为方便讨论,提出孔隙水体积与回弹水体积两概念。孔隙水体积为在存在于土颗粒孔隙未与橡皮膜接触水的体积,根据孔压增长机理,土受动荷载过程中孔隙水对孔压发展起正向贡献作用。回弹水定义为有效应力变化后运移至土颗粒与橡皮膜之间填补膜回弹空间的水的体积,等于橡皮膜回弹体积,这部分水体积的增加将对孔隙水压力的发展起负向作用。
Ramana方法符合对橡皮膜顺变效应机理认识,但对补水体积认识存在误区:向试样内补水不仅需要补偿动载阶段回弹水的体积。补偿的目的是对试样内饱和土的孔隙水进行补偿。如图2,随补水的进行,橡皮膜二次回弹对所注入水量起到即时分流作用,即向试样补偿的水不能完全用于补偿土孔隙中的水分流失。从总注水体积角度而言,因时刻存在的分流作用,Ramana方法将动力阶段橡皮膜回弹水体积作为补水总量,必然且始终有一部分顺应膜再回弹成为回弹水,而实际孔压贡献的孔隙水体积仍然得不到充分补偿,故无法100%消除橡皮膜顺变效应。且该方法未提出明确的注水截止标准,实际操作易出现“越补偿待补体积越大”的死循环,试验效率低,不易操作。
由前文可见,无论计算修正或仪器补偿方法,都需要开展橡皮膜嵌入体积测量作为基础,为确定修正系数或补偿体积提供关键数据,对此笔者已在前续文章展开研究[23],本文不再赘述研究细节。
2. 橡皮膜顺变性补偿新原理方法
2.1 仪器补偿新原理
本文提出一种考虑橡皮膜实时分流的试样补偿新原理。新原理以孔隙水得到100%补偿为目标。补偿过程中,向试样内部的总注水体积
V试样-cp 包含对孔隙水补偿体积V孔隙水-cp 以及补水过程中支持橡皮膜进一步回弹的回弹水体积V回弹水-cp 。其中对土的孔隙水补偿体积V孔隙水-cp 等于动力阶段产生有效应力变化值Δσ′3 所对应的回弹水的体积V回弹水-Δσ′3 ,具体可表达为V试样-cp=V孔隙水-cp+V回弹水-cp, (5) V孔隙水-cp=V回弹水-Δσ′3。 (6) 即若考虑橡皮膜的分流作用,向试样内补偿水的体积应大于动力阶段橡皮膜的回弹体积,才能够对试样土骨架内真正的完全补水。
补偿水体积可根据橡皮膜嵌入体积测量试验结果进一步确定。不排水条件下,有效应力、孔隙水体积、回弹水体积三者协调变化。以汶川地震液化砾性土为例,由双尺寸法可确定三轴试验中有效应力、各类型水体积随总注水体积变化曲线[23],如图3所示。基于这种协调关系与式(5),(6),可确定如下补水关系:
V试样-cp=k⋅V孔隙水-cp ,V孔隙水-cp=V回弹水-Δσ′3 ,k=V孔隙水-Δσ3'+V回弹水-Δσ′3V孔隙水-Δσ′3 ,} (7) 式中,k为体积调整系数,可根据预先的补水试验结果确定,
V回弹水-Δσ′3 与V孔隙水-Δσ′3 需由某动力阶段初始有效应力σ′3 、动力阶段有效应力变化值Δσ′3 与前序橡皮膜嵌入体积测量试验结果联合确定,见图3。2.2 仪器补偿新方法
根据仪器补偿新原理,新补偿方法将在饱和试样受动载试验中,开展多阶段单回合注水补偿,如图4。
单阶段补偿具体操作为:当饱和试样在某阶段动荷载作用下产生一定孔压增量后,暂停动力加载;随后,根据此动力阶段的初始有效应力
σ′3 、孔隙水压力增量Δu ,由预先在橡皮膜嵌入体积测量中得到的有效应力、孔隙水、回弹水三者关系[23],结合式(7),确定待补水体积V试样-cp ;最后一次性缓慢将这部分水注入试样,得到稳定后的孔压值即为该动力阶段无橡皮膜顺变性影响的孔隙水压力,并可继续动加载。补偿节点:试验中可根据施加的动应力的大小、实际孔隙水压力的增长程度、振动次数预判等综合决定补偿节点。建议两种节点选取方法:①孔压增量阈值,建议以单次或多次振动后,以残余孔隙水压力增长达到或大于5 kPa作为一次补偿节点,或根据研究需求适当放宽至10 kPa作为阈值;②振动次数预设,现有三轴仪器均可满足预设振动次数功能,因此试验中可根据动应力大小与孔压增长水平,依实际工况选取如5,10或15等固定振次作为补偿节点。以上两种节点选取方法可根据实际研究需求单独或混合使用。
补水速率:包含本文在内的补偿方法皆以待补水总体积作为某补偿阶段控制目标,但以往鲜对补水速率作出要求。对本文将应用的300 mm直径砾性土试样,参考反压饱和过程中压力控制器注水速率(0~150 mm3/s),以其中后期平稳注水阶段的5~50 mm3/s作为补偿注水速率,该阶段主要特点为孔压/反压接近预设值,注水体积缓慢增长至逐渐稳定。其他尺寸试样,可参考自身此阶段注水速率。具体试验中可根据含砾量调整注水速率(渗透系数会影响孔隙水压力场的均匀程度和平衡速率),具体原则:补水速率应随含砾量降低而逐步减小。同时本文建议利用试样双向同步补水[16]、控制补水速率与增加透水石厚度的手段相结合以缓解注水导致孔压场不均匀分布或局部孔压瞬时较大的可能。
新方法由于需要中断动荷载加载,可能会改变动应力历史。但对于均等固结条件下受正弦荷载的单元试验而言,土的密实度、排水条件、累计孔压增量、动应力(变)幅值、土的回弹模量等因素的影响相对动应力加载历史对结果影响更显加显著。如Finn模型[25]、Byrne模型[26]中以土的一维回弹模量、塑性体应变及增量、相对密度作为主要计算参数,Ishihara等[27]也认为应力作用速率在一定范围内影响很小。相比之下,橡皮膜顺变性对液化结果的影响则更加显著,如含砾量对液化势影响研究中,因橡皮膜顺变性的存在,现有研究成果[4, 28-30]截然不同甚至自相矛盾。据此,本文认为相比橡皮膜顺变性对土抗液化强度影响的显著性,可适当放松对动应力加载历史要求,后文将结合试验结果做进一步论证。
综上,本文仪器补偿新方法,原理上更科学合理;操作上,暂停式补偿与Ramana等[17]和Tokimatsu等[11]补偿方法研究一致,但单次补偿效率高,补偿截止标准明确,可实施性显著增强;成本上,相比Seed等[16]特制研发的低频振动补水系统,新方法由现有通用三轴仪即可完成补偿,开展成本大大降低。
3. 补偿试验设计
利用GDS大型动三轴仪,基于常规无补偿、Ramana、本文新补偿3种试验方法,开展砾性土液化试验,验证本文方法的理论可靠性和技术可行性。仪器装备ADVDPC型号高精度压力/控制器,为本文研究提供重要技术保障。
以汶川地震中发生实际液化的松柏村砾性土为基础料,制备3种含砾量的砾性土,级配曲线如图5,土的基本物性指标如表1。含砾量为5 mm以上粒组占土总质量百分比。试样干法装填,分十层击实或振密,相对密度50%,注意尖锐棱角朝向以防刺破橡皮膜,采用二氧化碳与反压饱和法,保证试样B值达到并维持在0.97以上,试样两端为4 cm厚粗烧陶瓷透水石。
表 1 试验砾性土基本参数及工况Table 1. Basic parameters of soils and tests名称 d10 /mmd20 /mmd50 /mm含砾量/% 试样直径/mm 橡皮膜厚度/mm 试验CSR 无补偿 Ramana 本文 SBC-G65(真液化土) 0.43 0.95 13.6 65 300 2 0.075 ~ 0.2 0.1/ 0.15 0.1/ 0.15 SBC-G80 0.8 5 25 80 300 2 0.125 ~ 0.25 0.1/ 0.15 0.1/ 0.15 SBC-G100 10 13 32 100 300 2 0.15 ~ 0.275 0.1/ 0.15 0.1/ 0.15 试验基本工况设计如表1所示,首先根据振动三轴试验方法[31]开展常规无补偿条件下砾性土的液化试验,以孔压等于初始有效应力100 kPa作为液化标准,另选择循环剪应力比CSR(
σd/2σ′3 )为0.1,0.15的工况分别开展两种橡皮膜补偿方法的液化-补偿试验,探究本文补偿方法在不同液化振次条件下的适用及可行情况。本文以孔隙水压力增长大于或等于5 kPa作为补偿节点,孔压稳定增长阶段辅助以预设振次5或10次提升工作效率,而孔压增长初期和接近液化阶段则以“少振次,多补偿”原则增加补偿节点,以记录残余孔压变化中的细节规律。振动频率1 Hz。4. 新方法验证
本节将从新方法的补偿过程与补偿结果两方面进行论证。补偿结果方面,通过确定砾性土孔压增量模型参数及橡皮膜修正系数
Cr ,进行橡皮膜顺变效应计算修正,作为补偿试验结果的理论参照;通过对比本文与Ramana两种补偿方法液化试验结果,探究新补偿方法的适用性与结果可靠性。4.1 新补偿方法对试样的影响分析
砾性土SBC-G65某补偿阶段,总补水体积、孔隙水压力、轴应变时程,如图6所示。补水速率为25 mm3/s。孔隙水压力随注水进程逐渐升高,未出现瞬时激增现象。注水停止后,孔压基本保持稳定,未出现显著下降。表明补偿过程中试样内部孔压场均匀变化,未因端部补水而出现局部孔压显著升高现象。可见本文增厚透水石、控制补水速率、双向同步补水措施可有效实现孔隙水压力均匀、同步变化。同时本文3种试验砾性土,含砾量较高,自身渗透系数较大,对注水补偿能够快速实现内部孔压自平衡。说明本文补偿方法未对试样内部孔压分布造成显著影响。
试验随注水补偿回弹,轴应变有所降低。根据笔者前续研究中的补水试验研究[23],这部分变形为有效应力下降导致的弹性回弹,基本不含塑性变形。当注水停止后,试样后续未出现显著变形。以图6为例,注水时轴应变由0.6073%回弹至0.6049%,稳定阶段轴应变变化区间为0.6048%~0.6050%,都属于小应变范畴。说明本文新方法在补偿过程中,试样不会出现显著变形甚至影响试样自身力学性质。
综上,新补偿方法在补偿过程中能够保证不对试样造成显著附加影响,从而确保补偿结果可靠性。
4.2 砾性土橡皮膜计算修正的关键参数
基于丰万玲等[32]和孙锐等[33]给出的均等固结饱和土孔压增量模型,可以建立砾性土孔压增量计算式:
¯UN=ΔuNσ′3-uN-1=C1,0NC2,0(τNσ′3-uN-1)A4,0, (8) 式中,
¯UN 为第N次动力荷载的孔压增量比,ΔuN 为第N次荷载的孔压增量,σ′3 为初始有效固结应力,uN-1 为第N-1动循环后的累计残余孔压,τN 为第N次动应力幅值,C1,0 ,C2,0 ,A4,0 为待定系数,与土性有关。该理论模型成熟,以此为基础发展考虑橡皮膜顺变性的计算修正模型,可作为补偿试验的理论参考。砾性土由常规无补偿试验确定的孔压增量模型基本参数如表2所示,参数确定细节及讨论将在后续文章中开展,本文不做赘述。但如图7,存在顺变效应影响的残余孔压实测值与基于表2的孔压计算值非常接近,变化规律一致,说明本文孔压增量模型参数可以有效地描述均等固结下不同含砾量砾性土的孔压变化情况。
表 2 砾性土橡皮膜顺变性计算修正的关键参数Table 2. Key parameters for calculation correction of membrane compliance of gravelly soils含砾量/% 孔压增量模型参数 橡皮膜顺变性修正系数 Cr C1,0 C2,0 A4,0 65 10.80 0.710 2.67 1.688 80 6.97 0.777 2.73 2.930 100 5.985 0.950 2.95 5.820 根据式(3),还需确定橡皮膜顺变性修正系数
Cr 。由笔者前续研究,已利用双尺寸法[23]获得几种砾性土在不同有效应力下单位面积橡皮膜嵌入体积Δvm 与土骨架体应变两εv 两个关键指标[23],如图8所示。Δvm 和εv 同有效应力σ′3 在半对数坐标系内呈线性关系,Δvm 关于lgσ′3 的斜率用Sm 表示,εv 随lgσ′3 变化斜率用Sr 表示,其本质为土骨架回弹系数,与回弹模量互为倒数。并结合各参数物理意义,可推导Cr 的计算式:Cr=ErEm=1/Sr1/(Sm⋅A/V)=Sm⋅ASr⋅V=4SmSr⋅D, (9) 式中,
Sr ,Sm 可由图8拟合得到,A为试样表面积,V为试样体积,D为试样直径。3种砾性土橡皮膜顺变性修正系数Cr ,如表2。由式(3),(8)与表2,可确定考虑橡皮膜顺变性的孔压结果。4.3 新方法的结果检验
不同动荷载条件下,砾性土由两种补偿方法、常规无补偿的三轴液化试验孔压结果对比,如图9,10所示。其中黑色实心圆点为存在橡皮膜顺变性的常规液化残余孔压结果,空心圆点为新补偿方法所得残余孔压,空心三角为Ramana方法所得残余孔压。
联立式(3),(8)与表2,可得对橡皮膜计算修正的残余孔压结果,如图10黑色实心方点所示。其中式(3),即Martin等[20]的基础模型,符合孔压发展机理,被作为橡皮膜校正的孔压理论解与参考值[10-11,14]。而式(8)为成熟的孔压增量模型,满足式(3)中
Δum 定义,且兼顾了周景星等[8]提出孔压增量模型需能够考虑已有的累积孔压水平的要求,能够反映动应力历史。由此,图10中理论计算结果可作为橡皮膜顺变性校正的理论参考值,并检验各补偿方法的试验结果。由图9孔压时程,本文新方法可成功使试样达到初始液化状态,但相比常规三轴试验,达到液化状态所需振动次数大大减少,验证了本文新方法的可行性,也说明橡皮膜顺变性对粗粒土液化特性影响显著。
图10对比新、老两补偿方法结果。本文新方法相比Ramana方法在相同试验条件下达到液化所需振动次数始终更少,土更迅速达到初始液化状态,此现象与两方法的总补水体积的大小关系相符。Ramana方法把回弹水作为全部补偿水,而新方法仅将回弹水用于补偿试样中土的孔隙水,总补水体积更高,也因此获得了数值更高、发展更快的补偿后孔隙水压力。
整理各方法达到初始液化所需振次于表3。结合图9,10,新方法达到液化所需平均振次仅为计算修正结果1.3倍,最大不超过1.7倍,而Ramana方法则平均需要2倍的振动次数,最高时振次接近3倍。可见新方法所得残余孔压规律更接近理论参考值,且不受含砾量影响,能够得到更贴近真实的动力试验结果。因此结合前文关于Ramana方法在补水体积方面的讨论,该方法并非以往认知中的100%补偿方法[19]。相较下,本文新方法能够更好消除橡皮顺变性对饱和土动力反应的影响,是更接近100%的补偿方法。
表 3 不同补偿方法达到初始液化所需要的振动次数Table 3. Vibration times required by different compensation methods to achieve initial liquefaction方法名称 SBC-G65 SBC-G80 SBC-G100 CSR=0.1 CSR =0.15 CSR =0.1 CSR =0.15 CSR =0.1 CSR =0.15 计算修正法 29.5 5 64.5 6.5 383.5 7.5 本文新补偿法 39.5 7 66.5 11.5 342.0 10.5 Ramana方法 58.0 14 100.0 17.5 635.0 11.5 注: 单次振动周期内试样达液化,按0.5次计。图10中,新补偿方法结果在中、高振次条件更接近理论计算值,低振次结果略有偏差,但与理论计算结果的一致性仍优于Ramana方法。新补偿方法同理论计算结果间的差异,本文认为与仪器施加动应力并非完全恒定有关。低振次工况,孔压增长迅速,土体刚度变化较快,三轴仪难以稳定施加大幅值动应力,即便配置应力自适应功能,但临近液化阶段动应力施加值与预设值差异会明显增加,达到液化所需振次也有所增加。中、高振次对应动应力较低,孔压发展较慢,动应力幅值较为稳定,所得结果与计算参考值更加一致。因此本文建议新补偿方法在较低振次的液化工况时,可适当降低振动频率,使动应力施加更充分,有利于获取更接近理论值的液化结果。另外膜嵌入体积测量结果可靠性对橡皮膜校正结果具有较大影响。
对比新方法、理论修正与常规试验结果发现,新方法补偿所得残余孔压结果仍呈现反S型趋势,土的基本动力特性未因新方法的暂停补偿而受明显影响,说明暂停式补偿新方法技术可行,所得结果有效合理。
5. 结论
橡皮膜顺变效应是影响砾性土三轴液化试验结果可靠性的关键,直接影响粗粒土动力特性研究的发展。
(1)提出橡皮膜顺变性补偿新原理。将土孔隙水得到100%补偿作为目标,以动力阶段橡皮膜回弹水体积与补偿过程补充土孔隙水体积二者相等为准则,认为总补水体积为补偿土孔隙水体积与补偿过程中伴随橡皮膜回弹水体积之和,合理考虑了橡皮膜对补偿水的分流作用,克服了现有补偿原理仅以动力导致的回弹水作为总补水体积而补偿不足的缺点。
(2)提出橡皮膜顺变性补偿新方法。液化试验中当孔压达到预设补偿节点后暂停动荷载,基于新原理由橡皮膜嵌入体积、阶段性初始有效应力与孔压增量确定待补水总体积,通过控制补给速率对试样进行单回合补水作为一次完整补偿环节;与现有补偿方法相比,新方法单补偿阶段仅需注水一次,无需多回合逼近,补偿、截止标准明确,可操作性更强,试验效率更高,开展成本更低。
(3)基于新仪器补偿方法成功开展液化试验,补偿过程对试样附加影响小;所得结果与成熟的计算修正理论结果进行对比,二者接近,均呈现一致的反S型;与现有方法比较,新方法得到的残余孔压发展更加接近理论修正结果,表明新方法是一种更接近试样100%补偿的试验方法。
-
表 1 新坐标轴的方向余弦
Table 1 Direction cosines of new axis
坐标 x y z x' l1 m1 n1 y' l2 m2 n2 z' l3 m3 n3 表 2 滑面及滑体方程提取
Table 2 Extraction of equations for sliding surface and sliding body
计算断面 滑面方程提取 K102+240 (x+1.71)2+(y−16.74)2=16.962 坐标点提取:(13.05,10.28),(17.10,12.40),(20.57,15.05),(23.21,19.05),(24.91,24.10) K102+260 (x+18.86)2+(y−44.65)2=45.552 坐标点提取:(12.12,10.24),(19.68,14.12),(25.24,18.50),(31.02,25.34),(34.70,31.96) K102+280 (x+8.41)2+(y−46.80)2=47.172 坐标点提取:(13.04,10.44),(20.16,14.04)、(26.24,18.56),(31.80,25.06),(35.54,32.04) 滑体方程 x2900+y21089+(z−21.72)2471.76=1 表 3 相关方程及初始参数
Table 3 Equations and initial parameters
方程 初始参数 滑面方程 z=21.72−(21.722−0.773x2−0.524y2)0.5 坡面方程 z=0.879x(0 m<x<8.986 m) z=7.92(8.986 m<x<11 m) z=1.31x-5.28(11 m<x<20 m) z=0.4x+11.72(x>20 m) 土层方程 土层1:(0<z<0.2143x+9.66) 土层2:(0.2143x+9.66<z<0.8793x) 参数初始化 种群规模NP=30,初始变异系数0.5, crl =0.1,cru =0.6;迭代次数G=8000锚固优化参数取值 锚固角度15°;vd≥2.5 m,hd≥2.5 m;2.5 m≤x≤16.5 m,-26 m≤y≤23 m -
[1] 郑颖人, 杨明成. 边坡稳定安全系数求解格式的分类统一[J]. 岩石力学与工程学报, 2004, 23(16): 2836-2841. doi: 10.3321/j.issn:1000-6915.2004.16.031 ZHENG Ying-ren, YANG Ming-cheng. Classification unification of calculation formulation for safety factor of slope stability[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(16): 2836-2841. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.16.031
[2] 康景文, 朱大勇. 一种边坡锚固力新的计算方法[J]. 四川建筑科学研究, 2005, 31(5): 72-74. https://www.cnki.com.cn/Article/CJFDTOTAL-ACZJ200505021.htm KANG Jing-wen, ZHU Da-yong. A new method for calculating the anchoring force of slope[J]. Sichuan Building Science, 2005, 31(5): 72-74. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ACZJ200505021.htm
[3] 赵洪领, 台佳佳, 朱大勇. 边坡稳定性计算中锚固力计入方式讨论[J]. 中国水利水电科学研究院学报, 2006, 4(4): 271-276. doi: 10.3969/j.issn.1672-3031.2006.04.006 ZHAO Hong-ling, TAI Jia-jia, ZHU Da-yong. Discussion on calculation method of anchorage force in slope stability calculation[J]. Journal of China Institute of Water Resources and Hydropower Research, 2006, 4(4): 271-276. (in Chinese) doi: 10.3969/j.issn.1672-3031.2006.04.006
[4] WANG Wen-jie, SONG Qian-qiang. Mechanical behavior of fully grouted GFRP rock bolts under the joint action of pre-tension load and blast dynamic load[J]. Tunneling and Underground Space Technology, 2018, 38(3): 82-91.
[5] ZHANG R, ZHAO J, WANG G X. Stability analysis of anchored soil slope based on finite element limit equilibrium method[J]. Mathematical Problems in Engineering, 2016(20): 1-8.
[6] 刘华丽, 李宏伟, 朱大勇. 锚固力在边坡加固中的工程应用[J]. 岩土力学, 2015, 36(2): 509-514. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201502032.htm LIU Hua-li, Li Hong-wei, ZHU Da-yong. A method for calculation of anchor force in slope reinforcement[J]. Rock and Soil Mechanics, 2015, 36(2): 509-514. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201502032.htm
[7] HARIYADI S M, YOSHIMI S. Experimental analysis of anchor bolt in concrete under the pull-out loading[J]. Procedia Engineering, 2017, 171: 926-933. doi: 10.1016/j.proeng.2017.01.391
[8] LAURA B M, TIJANI M, HADJ-HASSEN F. A new analytical solution to the mechanical behaviour of fully grouted rock bolts subjected to pull-out tests[J]. Construction & Building Materials, 2011, 25(2): 749-755.
[9] 雷军, 肖世国. 考虑滑体中锚拉力扩散效应的锚固边坡稳定性分析方法[J]. 岩土工程学报, 2019, 41(9): 1724-1730. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201909019.htm LEI Jun, XIAO Shi-guo. Stability analysis method for anchored slopes considering stress dispersion effect in sliding mass due to anchor tensile forces[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(9): 1724-1730. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201909019.htm
[10] 卢坤林, 王运敏, 朱大勇. 三维滑体锚固力计算方法及工程应用[J]. 岩土力学, 2017, 38(2): 501-506, 516. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201702026.htm LU Kun-lin, WANG Yun-min, ZHU Da-yong. A method for calculating reinforcing forces of 3D sliding mass and its engineering applications[J]. Rock and Soil Mechanics, 2017, 38(2): 501-506, 516. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201702026.htm
[11] LI Z, WEI J, YANG J. Stability calculation method of slope reinforced by prestressed anchor in process of excavation[J]. The Scientific World Journal, 2014(10): 1-7.
[12] 曹国金, 姜弘道, 熊红梅. 一种确定拉力型锚杆支护长度的方法[J]. 岩石力学与工程学报, 2003, 22(7): 1141-1145. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200307018.htm CAO Guo-jin, JIANG Hong-dao, XIONG Hong-mei. Calculation method of support length for stretched bolts[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(7): 1141-1145. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200307018.htm
[13] YANG G, ZHONG Z, ZHANG Y, et al. Optimal design of anchor cables for slope reinforcement based on stress and displacement fields[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2015, 7(4): 411-420.
[14] LI C, WU J, WANG J, et al. Layout and length optimization of anchor cables for reinforcing rock wedges[J]. Bulletin of Engineering Geology and the Environment, 2016, 75(4): 1399-1412.
[15] 石露, 李小春, 任伟, 等. 蚁群算法与遗传算法融合及其在边坡临界滑动面搜索中的应用[J]. 岩土力学, 2009, 30(11): 3486-3492. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200911057.htm SHI Lu, LI Xiao-chun, REN Wei, et al. Hybrid of ant colony algorithm and genetic algorithm and its application to searching critical slope slip surface[J]. Rock and Soil Mechanics, 2009, 30(11): 3486-3492. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200911057.htm
[16] HUANG Q J, ZHANG K, SONG J C, et al. Adaptive differential evolution with a Lagrange interpolation argument algorithm[J]. Information Sciences, 2018, 9(4): 180-202.
[17] 徐东方, 郭战伟, 侯春娟. 约束尺度和算子自适应变化的差分进化算法[J]. 数学的实践与认识, 2017(16): 151-158. https://www.cnki.com.cn/Article/CJFDTOTAL-SSJS201716018.htm XU Dong-fang, GUO Zhan-wei, HOU Chun-juan. Differential evolutionary algorithm of constraint scale and operator adaptive change[J]. Mathematics in Practice and Theory, 2017(16): 151-158. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SSJS201716018.htm
[18] 王凯. 层状弹性体系的力学分析与计算[M]. 北京: 科学出版社, 2009. WANG Kai. Mechanical Analysis and Calculation of Layered Elastic System[M]. Beijing: Science Press, 2009. (in Chinese)
[19] 卢坤林, 朱大勇, 甘文宁, 等. 一种边坡稳定性分析的三维极限平衡法及应用[J]. 岩土工程学报, 2013, 35(12): 2276-2282. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201312023.htm LU Kun-lin, ZHU Da-yong, GAN Wen-ning, et al. 3D limit equilibrium method for slope stability analysis and its application[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(12): 2276-2282. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201312023.htm
[20] 岩土锚杆与喷射混凝土支护工程技术规范:GB 50086—2015[S]. 2015. Technical code for engineering of ground anchoranges and shotcrete support: GB 50086—2015[S]. 2015. (in Chinese)
[21] 基于分片光滑技术三维边坡安全系数计算软件V2.0(登记号:2019SR0071874)[CP]. 3D slope safety factor calculation software based on piecewise smoothing technique V2.0 registration. No. 2019SR0071874[CP]. (in Chinese)
-
期刊类型引用(15)
1. 高睿,关鹏. 复杂地质条件下城市基坑开挖支护与变形研究. 建筑技术. 2025(01): 69-73 . 百度学术
2. 王浩. 改进型1+1复合冲孔灌注桩工艺在超深入岩排桩支护工程中的应用. 工程技术研究. 2025(02): 68-70 . 百度学术
3. 王旭东. 非平衡荷载下基坑开挖对围护结构的影响分析. 交通科技与管理. 2025(06): 152-155 . 百度学术
4. 叶征远. 不同偏压条件下基坑开挖对支护排桩变形特性的影响研究. 东北水利水电. 2024(07): 43-46+62+72 . 百度学术
5. 白锦烽. 基坑边坡土钉墙支护效果研究. 河南科技. 2024(17): 61-64 . 百度学术
6. 王金伦,刘慧敏,樊晨光. 某输变电工程深基坑降水止水措施及开挖变形研究. 河南科技. 2024(18): 64-68 . 百度学术
7. 邓盛双. 邻近公路深基坑围护结构设计与应用分析. 安徽建筑. 2023(01): 165-167 . 百度学术
8. 崔巍. 拉森钢板桩在逆作法深基坑中的应用研究. 辽宁省交通高等专科学校学报. 2023(01): 21-25 . 百度学术
9. 张秦军. 富水砂层条件下偏压深基坑开挖变形特性研究. 科技与创新. 2023(20): 1-5 . 百度学术
10. 赵平. 基于非对称荷载作用下的基坑变形特性研究. 菏泽学院学报. 2022(02): 65-70 . 百度学术
11. 刘涛,胡柏春,杨迅. 偏压条件下基坑排桩围护结构工作性状研究. 建筑结构. 2022(S1): 2595-2601 . 百度学术
12. 程馨玉,付成华,刘健,蒋林杰. 基坑开挖变形及其对周围建筑的影响分析. 人民珠江. 2022(09): 128-136 . 百度学术
13. 周国强,吴明明,王晶,顾文超,徐启良,方敏杰. 偏压作用下坑底土体加固对围护墙变形特性的影响研究. 岩土工程学报. 2022(S1): 201-206 . 本站查看
14. 张铎,惠海峰,王飞,李龙,谢忠. 偏压荷载下临近地铁深基坑变形特性数值模拟研究. 建筑结构. 2022(S2): 2262-2266 . 百度学术
15. 黄汉祥,张海亮,魏宏亮,许野,郑迪,边栋杰,李中浩. 偏压荷载对深基坑变形特性实测分析. 建筑结构. 2022(S2): 2284-2288 . 百度学术
其他类型引用(8)