• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

考虑土-结构相互作用倾斜状态文峰塔地震响应分析

安一劲, 韩鹏举, 秦健栋, 白向灵, 何斌, 王效渊

安一劲, 韩鹏举, 秦健栋, 白向灵, 何斌, 王效渊. 考虑土-结构相互作用倾斜状态文峰塔地震响应分析[J]. 岩土工程学报, 2023, 45(S2): 201-207. DOI: 10.11779/CJGE2023S20028
引用本文: 安一劲, 韩鹏举, 秦健栋, 白向灵, 何斌, 王效渊. 考虑土-结构相互作用倾斜状态文峰塔地震响应分析[J]. 岩土工程学报, 2023, 45(S2): 201-207. DOI: 10.11779/CJGE2023S20028
AN Yijing, HAN Pengju, QIN Jiandong, BAI Xiangling, HE Bin, WANG Xiaoyuan. Seismic response analysis of leaning Wenfeng Pagoda considering soil-structure interaction[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(S2): 201-207. DOI: 10.11779/CJGE2023S20028
Citation: AN Yijing, HAN Pengju, QIN Jiandong, BAI Xiangling, HE Bin, WANG Xiaoyuan. Seismic response analysis of leaning Wenfeng Pagoda considering soil-structure interaction[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(S2): 201-207. DOI: 10.11779/CJGE2023S20028

考虑土-结构相互作用倾斜状态文峰塔地震响应分析  English Version

基金项目: 

山西省文物局科技项目 RH2100001662

详细信息
    作者简介:

    安一劲(1999—),男,博士,主要从事建筑遗产健康评估等方面的科研工作。E-mail: anyijing0020@163.com

    通讯作者:

    韩鹏举, E-mail: 13834569544@163.com

  • 中图分类号: TU311.3

Seismic response analysis of leaning Wenfeng Pagoda considering soil-structure interaction

  • 摘要: 为研究在土-结构相互作用下倾斜状态古塔地震响应,利用ABAQUS有限元软件,采用等效线性化地基土及非线性化砖石材料,考虑土体与结构分离滑移的几何非线性,借助python语言添加黏弹性人工边界,建立了倾斜及未倾斜状态永祚寺文峰塔两种有限元模型,通过时程分析方法对比研究倾斜因素对古塔抗震性能影响。结果表明:多遇地震作用下,文峰塔历史最大倾斜量对其抗震性能影响较小,峰值位移、层间位移角及加速度放大系数较未倾斜状态增加较小;设防地震作用下,倾斜状态文峰塔残余位移大幅增加,下层层间位移角普遍放大,加速度放大系数有所提高,塔身在倾斜侧损伤面积及损伤程度急剧增大,同时基础损伤也不容忽视。分析结果可为同类倾斜状态高层密檐砖塔抗震保护提供借鉴作用。
    Abstract: In order to study the seismic response of the leaning pagoda under soil-structure interaction, two kinds of finite element models for the leaning and un-leaning Wenfeng Pagoda of Yongzuo Temple are established by using the ABAQUS finite element software, adopting the equivalent linearization of foundation soils and non-linearization of masonry materials, considering the geometric non-linearities of the soil-structure separation and slip, and adding the visco-elastic artificial boundaries with the help of python language. By comparing the two models through the time-range analysis method, the effects of inclination factor on the seismic performance of the pagoda are investigated. The results show that under the action of small earthquakes, the historical maximum tilt has less impact on its seismic performance, and the peak displacement, section displacement angle and acceleration amplification coefficient have small increase. Under the action of middle earthquakes, the residual displacement of the leaning pagoda increases greatly, the lower layer of the pagoda section displacement angle is generally enlarged, the acceleration amplification coefficient increases, the damage area and degree of the pagoda in the leaning side increase dramatically, and the foundation damage should not be neglected. The analysis results can provide a reference for the seismic protection of similar leaning high-rise dense eaves brick pagodas.
  • 砾性土是对砾质土、砾砂、粉砾土、砂砾土和砾类土的统称[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应用甚少,未能取得更深入发展,阻碍着液化理论深化,导致当前的土壤液化研究,尤其是大颗粒粗粒土的液化单元试验中,对其常避而不谈或选择性忽略。

    本文基于橡皮膜顺变效应影响机理,以仪器补偿作为消除橡皮膜顺变性研究的主要手段,分析现有补偿手段存在的理论缺陷,提出一种对饱和试样进行补偿的新原理和新方法,结合砾性土的橡皮膜嵌入体积测量的前序成果,利用大型动三轴仪,以汶川地震中液化砾性土为基础料,开展新型补偿试验,通过与已有补偿方法、计算修正后无橡皮膜顺变性的理论结果等对比,验证新方法的技术可行性,理论可靠性和方法适用性。研究成果为开展粗粒动力试验研究提供新途径,为深化土动力学理论发展提供新动力。

    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],这在以往难度较大。因此计算修正方法原理科学可靠但实际应用较少。

    Ramana等[17]提出一种利用压力控制器的补偿方法,简称Ramana方法(下同),成为后续学者研发更高级补水仪器和系统的核心理论和技术源头[11, 16],被认为可达100%补偿效果[19]。其无需进行常规液化试验,可直接由补偿液化试验得到土的真实抗液化强度。基本原理是将液化试验过程中,某动力阶段用于回弹橡皮膜的水再补回试样中,来保证土试样本身的不排水状态,可表达为

    Vcp=Vd, (4)

    式中,Vcp为向试样补偿水的体积,Vd为因产生有效应力变化Δσ3而支持橡皮膜回弹的水的体积,由橡皮膜嵌入体积与有效应力关系确定,如图1中黑色实线所示。

    图  1  传统橡皮膜补偿的工作示意图
    Figure  1.  Working diagram of traditional membrane compensation

    操作方法为,获取膜嵌入体积与有效应力关系,如图1黑色实线,随后由某动力阶段有效应力变量Δσ3,读取其导致的橡皮膜回弹体积ΔVm1,将该体积的水及时补充到试样内,可得无橡皮膜顺变影响的试验结果。考虑注水引起新的孔压增长导致橡皮膜顺变,如图1箭头与虚线所示需多回合注水,反复逼近使总注水体积与有效应力增量重新在曲线上重合,一次补偿工作完成,并可以开始下一阶段动力试验。

    为方便讨论,提出孔隙水体积与回弹水体积两概念。孔隙水体积为在存在于土颗粒孔隙未与橡皮膜接触水的体积,根据孔压增长机理,土受动荷载过程中孔隙水对孔压发展起正向贡献作用。回弹水定义为有效应力变化后运移至土颗粒与橡皮膜之间填补膜回弹空间的水的体积,等于橡皮膜回弹体积,这部分水体积的增加将对孔隙水压力的发展起负向作用。

    Ramana方法符合对橡皮膜顺变效应机理认识,但对补水体积认识存在误区:向试样内补水不仅需要补偿动载阶段回弹水的体积。补偿的目的是对试样内饱和土的孔隙水进行补偿。如图2,随补水的进行,橡皮膜二次回弹对所注入水量起到即时分流作用,即向试样补偿的水不能完全用于补偿土孔隙中的水分流失。从总注水体积角度而言,因时刻存在的分流作用,Ramana方法将动力阶段橡皮膜回弹水体积作为补水总量,必然且始终有一部分顺应膜再回弹成为回弹水,而实际孔压贡献的孔隙水体积仍然得不到充分补偿,故无法100%消除橡皮膜顺变效应。且该方法未提出明确的注水截止标准,实际操作易出现“越补偿待补体积越大”的死循环,试验效率低,不易操作。

    图  2  补水过程中橡皮膜回弹的分流作用
    Figure  2.  Effects of distributary of membrane rebound during compensation

    由前文可见,无论计算修正或仪器补偿方法,都需要开展橡皮膜嵌入体积测量作为基础,为确定修正系数或补偿体积提供关键数据,对此笔者已在前续文章展开研究[23],本文不再赘述研究细节。

    本文提出一种考虑橡皮膜实时分流的试样补偿新原理。新原理以孔隙水得到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),可确定如下补水关系:

    图  3  试样中各类型水体积、有效应力随注水体积变化示意图
    Figure  3.  Diagram of volume change of different water and effective stresses with water injection
    V-cp=kV-cp  ,V-cp=V-Δσ3  ,k=V-Δσ3'+V-Δσ3V-Δσ3 ,} (7)

    式中,k为体积调整系数,可根据预先的补水试验结果确定,V-Δσ3V-Δσ3需由某动力阶段初始有效应力σ3、动力阶段有效应力变化值Δσ3与前序橡皮膜嵌入体积测量试验结果联合确定,见图3

    根据仪器补偿新原理,新补偿方法将在饱和试样受动载试验中,开展多阶段单回合注水补偿,如图4

    图  4  新补偿方法操作示意图
    Figure  4.  Operation diagram of new compensation method

    单阶段补偿具体操作为:当饱和试样在某阶段动荷载作用下产生一定孔压增量后,暂停动力加载;随后,根据此动力阶段的初始有效应力σ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]特制研发的低频振动补水系统,新方法由现有通用三轴仪即可完成补偿,开展成本大大降低。

    利用GDS大型动三轴仪,基于常规无补偿、Ramana、本文新补偿3种试验方法,开展砾性土液化试验,验证本文方法的理论可靠性和技术可行性。仪器装备ADVDPC型号高精度压力/控制器,为本文研究提供重要技术保障。

    以汶川地震中发生实际液化的松柏村砾性土为基础料,制备3种含砾量的砾性土,级配曲线如图5,土的基本物性指标如表1。含砾量为5 mm以上粒组占土总质量百分比。试样干法装填,分十层击实或振密,相对密度50%,注意尖锐棱角朝向以防刺破橡皮膜,采用二氧化碳与反压饱和法,保证试样B值达到并维持在0.97以上,试样两端为4 cm厚粗烧陶瓷透水石。

    图  5  砾性土试样级配曲线
    Figure  5.  Gradation curves of gravelly soils
    表  1  试验砾性土基本参数及工况
    Table  1.  Basic parameters of soils and tests
    名称d10/mmd20/mmd50/mm含砾量/%试样直径/mm橡皮膜厚度/mm试验CSR
    无补偿Ramana本文
    SBC-G65(真液化土)0.430.9513.66530020.075 ~ 0.20.1/ 0.150.1/ 0.15
    SBC-G800.85258030020.125 ~ 0.250.1/ 0.150.1/ 0.15
    SBC-G10010133210030020.15 ~ 0.2750.1/ 0.150.1/ 0.15
    下载: 导出CSV 
    | 显示表格

    试验基本工况设计如表1所示,首先根据振动三轴试验方法[31]开展常规无补偿条件下砾性土的液化试验,以孔压等于初始有效应力100 kPa作为液化标准,另选择循环剪应力比CSR(σd/2σ3)为0.1,0.15的工况分别开展两种橡皮膜补偿方法的液化-补偿试验,探究本文补偿方法在不同液化振次条件下的适用及可行情况。本文以孔隙水压力增长大于或等于5 kPa作为补偿节点,孔压稳定增长阶段辅助以预设振次5或10次提升工作效率,而孔压增长初期和接近液化阶段则以“少振次,多补偿”原则增加补偿节点,以记录残余孔压变化中的细节规律。振动频率1 Hz。

    本节将从新方法的补偿过程与补偿结果两方面进行论证。补偿结果方面,通过确定砾性土孔压增量模型参数及橡皮膜修正系数Cr,进行橡皮膜顺变效应计算修正,作为补偿试验结果的理论参照;通过对比本文与Ramana两种补偿方法液化试验结果,探究新补偿方法的适用性与结果可靠性。

    砾性土SBC-G65某补偿阶段,总补水体积、孔隙水压力、轴应变时程,如图6所示。补水速率为25 mm3/s。孔隙水压力随注水进程逐渐升高,未出现瞬时激增现象。注水停止后,孔压基本保持稳定,未出现显著下降。表明补偿过程中试样内部孔压场均匀变化,未因端部补水而出现局部孔压显著升高现象。可见本文增厚透水石、控制补水速率、双向同步补水措施可有效实现孔隙水压力均匀、同步变化。同时本文3种试验砾性土,含砾量较高,自身渗透系数较大,对注水补偿能够快速实现内部孔压自平衡。说明本文补偿方法未对试样内部孔压分布造成显著影响。

    图  6  某补偿阶段孔压、轴应变随注水体积变化时程
    Figure  6.  Time histories of pore pressure and axial strain with water injection during a compensation stage

    试验随注水补偿回弹,轴应变有所降低。根据笔者前续研究中的补水试验研究[23],这部分变形为有效应力下降导致的弹性回弹,基本不含塑性变形。当注水停止后,试样后续未出现显著变形。以图6为例,注水时轴应变由0.6073%回弹至0.6049%,稳定阶段轴应变变化区间为0.6048%~0.6050%,都属于小应变范畴。说明本文新方法在补偿过程中,试样不会出现显著变形甚至影响试样自身力学性质。

    综上,新补偿方法在补偿过程中能够保证不对试样造成显著附加影响,从而确保补偿结果可靠性。

    基于丰万玲等[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的孔压计算值非常接近,变化规律一致,说明本文孔压增量模型参数可以有效地描述均等固结下不同含砾量砾性土的孔压变化情况。

    图  7  不同含砾量土的残余孔压实测值与计算值代表性对比(σ3=100 kPa)
    Figure  7.  Representative comparison between measured and calculated residual pore pressures of different gravelly soils
    表  2  砾性土橡皮膜顺变性计算修正的关键参数
    Table  2.  Key parameters for calculation correction of membrane compliance of gravelly soils
    含砾量/%孔压增量模型参数橡皮膜顺变性修正系数Cr
    C1,0C2,0A4,0
    6510.800.7102.671.688
    806.970.7772.732.930
    1005.9850.9502.955.820
    下载: 导出CSV 
    | 显示表格

    根据式(3),还需确定橡皮膜顺变性修正系数Cr。由笔者前续研究,已利用双尺寸法[23]获得几种砾性土在不同有效应力下单位面积橡皮膜嵌入体积Δvm与土骨架体应变两εv两个关键指标[23],如图8所示。Δvmεv同有效应力σ3在半对数坐标系内呈线性关系,Δvm关于lgσ3的斜率用Sm表示,εvlgσ3变化斜率用Sr表示,其本质为土骨架回弹系数,与回弹模量互为倒数。并结合各参数物理意义,可推导Cr的计算式:

    图  8  Δvmεv随有效应力σ3变化的关系图
    Figure  8.  Relationship betweenΔvm, εvand effective stress σ3
    Cr=ErEm=1/Sr1/(SmA/V)=SmASrV=4SmSrD, (9)

    式中,Sr,Sm可由图8拟合得到,A为试样表面积,V为试样体积,D为试样直径。3种砾性土橡皮膜顺变性修正系数Cr,如表2。由式(3),(8)与表2,可确定考虑橡皮膜顺变性的孔压结果。

    不同动荷载条件下,砾性土由两种补偿方法、常规无补偿的三轴液化试验孔压结果对比,如图9,10所示。其中黑色实心圆点为存在橡皮膜顺变性的常规液化残余孔压结果,空心圆点为新补偿方法所得残余孔压,空心三角为Ramana方法所得残余孔压。

    图  9  新补偿方法孔压时程与常规三轴试验结果对比
    Figure  9.  Comparison of time histories of pore pressure from new compensation method and traditional traixial test
    图  10  两种补偿方法所得孔压时程与理论计算结果对比(σ3=100 kPa)
    Figure  10.  Comparison of time histories of pore pressure from 2 compensation tests and theoretical results

    联立式(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-G65SBC-G80SBC-G100
    CSR=0.1CSR =0.15CSR =0.1CSR =0.15CSR =0.1CSR =0.15
    计算修正法29.5564.56.5383.57.5
    本文新补偿法39.5766.511.5342.010.5
    Ramana方法58.014100.017.5635.011.5
    注:单次振动周期内试样达液化,按0.5次计。
    下载: 导出CSV 
    | 显示表格

    图10中,新补偿方法结果在中、高振次条件更接近理论计算值,低振次结果略有偏差,但与理论计算结果的一致性仍优于Ramana方法。新补偿方法同理论计算结果间的差异,本文认为与仪器施加动应力并非完全恒定有关。低振次工况,孔压增长迅速,土体刚度变化较快,三轴仪难以稳定施加大幅值动应力,即便配置应力自适应功能,但临近液化阶段动应力施加值与预设值差异会明显增加,达到液化所需振次也有所增加。中、高振次对应动应力较低,孔压发展较慢,动应力幅值较为稳定,所得结果与计算参考值更加一致。因此本文建议新补偿方法在较低振次的液化工况时,可适当降低振动频率,使动应力施加更充分,有利于获取更接近理论值的液化结果。另外膜嵌入体积测量结果可靠性对橡皮膜校正结果具有较大影响。

    对比新方法、理论修正与常规试验结果发现,新方法补偿所得残余孔压结果仍呈现反S型趋势,土的基本动力特性未因新方法的暂停补偿而受明显影响,说明暂停式补偿新方法技术可行,所得结果有效合理。

    橡皮膜顺变效应是影响砾性土三轴液化试验结果可靠性的关键,直接影响粗粒土动力特性研究的发展。

    (1)提出橡皮膜顺变性补偿新原理。将土孔隙水得到100%补偿作为目标,以动力阶段橡皮膜回弹水体积与补偿过程补充土孔隙水体积二者相等为准则,认为总补水体积为补偿土孔隙水体积与补偿过程中伴随橡皮膜回弹水体积之和,合理考虑了橡皮膜对补偿水的分流作用,克服了现有补偿原理仅以动力导致的回弹水作为总补水体积而补偿不足的缺点。

    (2)提出橡皮膜顺变性补偿新方法。液化试验中当孔压达到预设补偿节点后暂停动荷载,基于新原理由橡皮膜嵌入体积、阶段性初始有效应力与孔压增量确定待补水总体积,通过控制补给速率对试样进行单回合补水作为一次完整补偿环节;与现有补偿方法相比,新方法单补偿阶段仅需注水一次,无需多回合逼近,补偿、截止标准明确,可操作性更强,试验效率更高,开展成本更低。

    (3)基于新仪器补偿方法成功开展液化试验,补偿过程对试样附加影响小;所得结果与成熟的计算修正理论结果进行对比,二者接近,均呈现一致的反S型;与现有方法比较,新方法得到的残余孔压发展更加接近理论修正结果,表明新方法是一种更接近试样100%补偿的试验方法。

  • 图  1   黏弹性边界示意图

    Figure  1.   Schematic diagram of viscoelastic boundary

    图  2   文峰塔土-结构相互作用模型

    Figure  2.   Model for soil-structure interaction of Wenfeng Pagoda

    图  3   地震波加速度时程曲线

    Figure  3.   Seismic acceleration-time curves

    图  4   多遇烈度El-Centro波作用下文峰塔x向位移时程曲线

    Figure  4.   Time-history curves of x-directional displacement under action of El-Centro wave with small intensities

    图  5   设防烈度El-Centro波作用下文峰塔x向位移时程曲线

    Figure  5.   Time-history curves of x -directional displacement under action of El-Centro wave with middle intensities

    图  6   层间位移角

    Figure  6.   Section displacement angles

    图  7   加速度放大系数

    Figure  7.   Acceleration amplification factors

    图  8   El-Centro波作用下未倾斜状态古塔损伤云图

    Figure  8.   Damage cloud of Wenfeng Pagoda in untilted state under action of El-Centro wave

    图  9   El-Centro波作用下倾斜状态古塔损伤云图

    Figure  9.   Damage cloud of Wenfeng Pagoda in inclined state under.action of El-Centro wave

    图  10   古塔基础受拉损伤因子时程

    Figure  10.   Time histories of tensile damage factors of foundation

    表  1   土层物理力学参数

    Table  1   Physical and mechanical parameters of soil layer

    土类 密度/(kg·m-3) 泊松比 黏聚力/kPa 内摩擦角/(°) 弹性模量/MPa
    粉土 1950 0.29 16 24 217.76
    下载: 导出CSV

    表  2   文峰塔两种模型前二阶自振频率及差值

    Table  2   First two natural frequencies and differences of three models for Wenfeng Pagoda

    模态阶数 非倾斜模型 倾斜模型
    实测频率f0/Hz 频率f1/Hz 误差e1-0/% 频率f2/Hz 差值e1-0/%
    一阶X 0.732 0.718 -1.91 0.662 -7.80
    一阶Z 0.732 0.719 -1.78 0.700 -2.64
    二阶X 3.125 3.393 8.58 3.391 -0.06
    二阶Z 3.125 3.410 9.12 3.407 -0.09
    注:eij=(fifj)/(fifj)fjfj
    下载: 导出CSV
  • [1] 李犇. 倾斜砖石古塔的有限元模拟及危险性分析[D]. 南京: 东南大学, 2022.

    LI Ben. Finite Element Simulation and Risk Analysis of Inclined Masonry Ancient Tower[D]. Nanjing: Southeast University, 2022. (in Chinese)

    [2] 崔晓玲, 李凯, 卢俊龙, 等. 倾斜砖石古塔随机振动测试分析[J]. 建筑结构, 2017, 47(2): 70-74. https://www.cnki.com.cn/Article/CJFDTOTAL-JCJG201702013.htm

    CUI Xiaoling, LI Kai, LU Junlong, et al. Analysis of random vibration test of an inclined ancient masonry pagoda[J]. Building Structure, 2017, 47(2): 70-74. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JCJG201702013.htm

    [3] 严绍军, 方云, 唐朝晖. 云南大理崇圣寺南北塔倾斜成因研究[J]. 岩土力学, 2012, 33(11): 3367-3373. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201211031.htm

    YAN Shaojun, FANG Yun, TANG Zhaohui. Research on inclination reason of north and south pagodas of Chongsheng Temple, Dali, Yunnan Province[J]. Rock and Soil Mechanics, 2012, 33(11): 3367-3373. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201211031.htm

    [4] 李宏男, 徐静. 考虑桩–土–结构动力相互作用的输电塔线体系简化抗震计算模型[J]. 岩土工程学报, 2009, 31(11): 1763-1767. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200911024.htm

    LI Hongnan, XU Jing. Simplified aseismic calculation model for transmission tower-line system considering pile-soil dynamic interaction[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(11): 1763-1767. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200911024.htm

    [5]

    CASOLO S, DIANA V, UVA G. Influence of soil deformability on the seismic response of a masonry tower[J]. Bulletin of Earthquake Engineering, 2017, 15(5): 1991-2014. doi: 10.1007/s10518-016-0061-y

    [6] 刘羽. 汶川地震中龙护舍利塔的损伤特征研究[D]. 扬州: 扬州大学, 2012.

    LIU Yu. Research on Seismic Damage Feature of Longhu Pagoda in Wenchuan Earthquake[D]. Yangzhou: Yangzhou University, 2012. (in Chinese)

    [7]

    BAO X, LIU J B, LI S T, et al. A new viscoelastic artificial boundary with improved numerical stability in explicit calculation of wave propagation problems in infinite domains[J]. Computers and Geotechnics, 2022, 145: 104698. doi: 10.1016/j.compgeo.2022.104698

    [8] 马笙杰, 迟明杰, 陈红娟, 等. 黏弹性人工边界在ABAQUS中的实现及地震动输入方法的比较研究[J]. 岩石力学与工程学报, 2020, 39(7): 1445-1457. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202007013.htm

    MA Shengjie, CHI Mingjie, CHEN Hongjuan, et al. Implementation of viscous-spring boundary in ABAQUS and comparative study on seismic motion input methods[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(7): 1445-1457. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202007013.htm

    [9] 沈继美. 砌体匀质化过程的数值模拟方法与应用研究[D]. 长沙: 长沙理工大学, 2012.

    SHEN Jimei. Study on Homogenization Process of Masonry and its Applications using Numerical Simulation[D]. Changsha: Changsha University of Science & Technology, 2012. (in Chinese)

    [10] 杨卫忠. 砌体受压本构关系模型[J]. 建筑结构, 2008, 38(10): 80-82. https://www.cnki.com.cn/Article/CJFDTOTAL-JCJG200810026.htm

    YANG Weizhong. Constitutive relationship model for masonry materials in compression[J]. Building Structure, 2008, 38(10): 80-82. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JCJG200810026.htm

    [11] 郑妮娜, 李英民, 潘毅. 芯柱式构造柱约束的低层砌体结构抗震性能[J]. 西南交通大学学报, 2011, 46(1): 24-29, 55. https://www.cnki.com.cn/Article/CJFDTOTAL-XNJT201101006.htm

    ZHENG Nina, LI Yingmin, PAN Yi. Seismic behavior of low masonry structure with core-tie-columns[J]. Journal of Southwest Jiaotong University, 2011, 46(1): 24-29, 55. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XNJT201101006.htm

    [12] 李世温, 李庆玲. 太原永祚寺东塔维修后的抗震能力初探[C]//建筑物鉴定与加固改造第五届全国学术讨论会. 南京, 2000.

    LI Shiwen, LI Qinglin. A preliminary investigation of the seismic resistance of the repaired Yongzuo Temple East Pagoda in Taiyuan[C]// Proceedings of the Fifth National Symposium on Building Identification and Strengthening and Renovation. Nanjing, 2000. (in Chinese)

    [13] 施楚贤. 砌体结构理论与设计[M]. 2版. 北京: 中国建筑工业出版社, 2003: 78-80.

    SHI Chuxian. Theory and Design of Masonry Structure[M]. 2nd ed. Beijing: China Architecture & Building Press, 2003: 78-80. (in Chinese)

    [14] 杜修力, 许紫刚, 许成顺, 等. 基于等效线性化的土–地下结构整体动力时程分析方法研究[J]. 岩土工程学报, 2018, 40(12): 2155-2163. doi: 10.11779/CJGE201812001

    DU Xiuli, XU Zigang, XU Chengshun, et al. Time-history analysis method for soil-underground structure system based on equivalent linear method[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(12): 2155-2163. (in Chinese) doi: 10.11779/CJGE201812001

    [15]

    BARDET J, ICHII K, LIN C. EERA-A Computer Program for Equivalent-linear Earthquake Site Response Analyses of Layered Soil Deposits[R]. Los Angeles: University of Southern California, 2000.

    [16] 魏剑伟, 李铁英. 古砖塔动力特性测试与分析研究[J]. 山西建筑, 2002, 28(1): 18-19. https://www.cnki.com.cn/Article/CJFDTOTAL-JZSX200201010.htm

    WEI Jianwei, LI Tieying. Research on the test and analysis of ancient brick tower dynamic characteristics[J]. Shanxi Agchitecture, 2002, 28(1): 18-19. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JZSX200201010.htm

    [17]

    FLORA A, CHIARADONNA A, DE SANCTIS L, et al. Understanding the Damages Caused by the 1999 Kocaeli Earthquake on One of the Towers of the Theodosian Walls of Constantinople[J]. International Journal of Architectural Heritage, 2022, 16(7): 1076-1100. doi: 10.1080/15583058.2020.1864512

    [18]

    BI Z H, CUI B H, ZHAI Y F. Research on seismic input method of layered ground foundation[J]. IEEE Access, 2021, 9: 54820-54830. doi: 10.1109/ACCESS.2021.3070479

    [19]

    JIANLI Y, SHENGCAI L. Study of the seismic damage regularity of ancient masonry pagodas in the 2008 Wenchuan Earthquake[C]// WIT Transactions on The Built Environment", "Earthquake Resistant Engineering Structures IX. A Coruña, 2013.

  • 期刊类型引用(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)

图(10)  /  表(2)
计量
  • 文章访问数:  161
  • HTML全文浏览量:  23
  • PDF下载量:  14
  • 被引次数: 23
出版历程
  • 收稿日期:  2023-11-29
  • 网络出版日期:  2024-04-19
  • 刊出日期:  2023-11-30

目录

/

返回文章
返回