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.
-
0. 引言
非饱和状态是影响土体力学特性的重要因素。基质吸力的增加直接降低了土体的含水率,提高了土体的抗剪强度和模量,进而对岩土工程的稳定性产生显著影响。因此,土水特征研究一直是非饱和土研究的热点。
过去,人们对细粒土(如黏土、粉质黏土、砂土等)的土水特征进行了大量研究,但对大颗粒粗粒堆积土的土水特征研究较少,主要是受限于现有测试设备的尺寸,它们只能对小尺寸的土样进行测试。然而,粗粒堆积土含较多10~60 mm的颗粒,因此小尺寸试验设备不能满足大颗粒土体试验的基本要求(试样土样的直径需大于或等于土样最大粒径的5倍)。因此,目前鲜有报道粗粒堆积土土-水特征的研究成果。
而且,人们普遍认为,粗粒堆积土的持水能力较弱,对土体的变形和强度几乎没有影响,然而,目前的研究结果表明,土体的SWCC受许多因素的影响,如密度、干湿循环、应力历史、细粒土含量等[1-2]。其中,细粒土的含量对土体SWCC有更显著的影响。粗粒堆积土中不仅含有大颗粒,而且含有大量细粒土。随着细粒土含量的增加,土体的SWCC不可避免地会发生变化,土体的力学特性将逐渐受到影响。因此,在细粒土含量增加的条件下,对粗粒堆积土土水特性的研究具有重要的理论意义和实际应用价值。
土水特性试验是一种耗时的试验,因为在每个阶段的气体压力下,土体中的气液平衡时间较长。因此,一般来说,研究土体SWCC的快速方法是建立SWCC的经验公式[3],例如Van Genuishen公式和Fredlund公式属于这种预测方法。这些模型只需要使用少量的试验数据来拟合经验公式中的关键参数,就可以得到完整的SWCC。然而,使用经验公式获得SWCC也是基于一些试验数据。为了进一步简化SWC的应用,一些学者在总结了大量土体的SWCC后,建立了与土颗粒分布相关的SWCC预测模型。例如,Arya等[4]、Tyler等[5]、Sweijen等[6]都建立了许多基于土颗粒级配的SWCC预测模型。这些模型基于细粒土的SWCC预测,其对粗粒堆积土的适用性仍需进一步研究。
针对粗粒堆积土SWCC研究的不足,本文首先介绍了一种新开发的非饱和粗粒堆积土三轴试验系统,该系统可以对粗粒堆积土进行大尺寸SWCC试验。在此基础上,总结了粗粒堆积土的SWCC试验方法和吸力平衡标准,研究了细粒土含量增加条件下粗粒堆积土SWCC的变化规律。最后,利用经典的AP模型,分析并得出了适用于粗粒堆积土SWCC的预测方法。
1. 非饱和粗粒堆积土三轴仪
大型非饱和粗粒堆积土动静态三轴试验系统由压力控制系统、轴压加载系统、数据采集系统以及软件控制与分析系统组成。其中,涉及粗粒堆积土SWCC试验的主要系统和部件包括:陶土板底座、气压控制装置和排水测量系统,见图 1。大型非饱和三轴仪试样尺寸为150 mm×300 mm,该试验系统允许的最大颗粒粒径最大可达30 mm,这就大大消除了试验时替代大颗粒土所带来的偏差,试得试验结果更加符合实际。陶土板底座由3个具有高进气值的小陶土板构成,在有效隔绝土样气体排放的同时,避免了装样过程中陶土板受冲击破裂。
2. 非饱和粗粒堆积土SWCC试验方案
2.1 非饱和粗粒堆积土试样制备
试验土样取自某高铁站粗粒堆积土高填方边坡。为了研究不同细粒料含量(颗粒粒径小于0.075 mm)对粗粒堆积土土水特征曲线的影响,设计了3组细粒料含量不同的试样,细粒料质量含量分别为4.525%,9.05%,13.575%,试样的土颗粒级配见图 2。3组试样在粒径为5~30 mm的占比几乎相同,颗粒级配差异主要体现在5 mm粒径以下部分,并重点关注0.075 mm以下部分的含量。3组试样的干密度均控制为1.8 g/cm3。将配置好的试样搅拌均匀,分5次平均加入制样桶中。每次加样后应压实试样,控制单次加样高度为60 mm。
2.2 试验步骤与吸力平衡标准
粗粒土SWCC试验步骤参考细粒土试验步骤,不同的是,由于粗粒堆积土样品尺寸大、饱和含水率高,试验期间的气液平衡标准不同于细粒土试验。
小颗粒土壤水分特征试验的基质吸力平衡准则大致可分为两类:第一类是根据排水量与样本量的比例进行控制。例如,在褚进晶等[7]对黏土水的特性试验中,平衡标准是每2 h排水量小于样品体积的0.05%。第二类是控制排水量。例如,Pham[8]采用的平衡标准是,24 h内土壤样品的排水质量小于0.1 g。
由于粗粒堆积土样本的含水率相对较高,第二类控制标准对于大样本粗粒堆积土来说过于严格,因此根据第一类平衡标准来判断粗粒堆积土吸力平衡更合适。经过多次试验,发现以1 mL/h的出水速率作为吸力平衡的控制标准更为合理,即每1 h的排量小于样品体积的0.04%。
3. 细颗粒含量对粗粒堆积土SWCC影响
3.1 不同细颗粒含量试样试验结果
图 3给出了不同细颗粒含量试样的SWCC试验结果。由于试验数据点略少,采用经典VG模型对测得的试验数据进行拟合,以研究土壤水分特征参数,实测数据及拟合曲线见图 3(a)~(c)。
从图 3可以看出,与粉质黏土等细粒土相比,粗粒堆积土的持水能力相对较弱,土体中的水可以在较小的气压下排出,土体中的基质吸力相对较小。然而,随着细粒土含量的增加,粗粒堆积土的持水能力不断增加,土壤的进气值也相应增加,即Sa3 > Sa2 > Sa1,Sa1趋于0,未在图中展出。此外,随着细土含量的增加,土体的残余基质吸力也在增加,即ψ3 > ψ2 > ψ1,这也表明土壤的持水能力在增加,土体持水能力的提高必然会影响土体的变形和强度特性。
3.2 试验设备及方法可靠性验证
目前,关于粗粒堆积土土水特征的研究很少,仅通过前人极其有限的研究结果[9]作为对比数据,验证本文试验数据的合理性。对比试验数据见图 4。
通过对比可以发现,尽管土颗粒的分布不同,但粗粒堆积土的土水分布特征曲线具有明显的相似性,基质吸力范围基本相同。随着粗粒堆积土不均匀系数的增加,相同基质吸力下的饱和度逐渐升高,表明土体持水能力不断提高。上述对比结果说明了实验数据的合理性,证明了利用该三轴仪器进行试验的可靠性。
4. 基于AP模型的粗粒堆积土SWCC曲线预测方法
4.1 Arya-Paris模型
Arya-Paris模型[4]从土样土水特征曲线与颗分曲线的相似性出发,将土样视为由圆球状颗粒以及圆柱状孔隙组成的多孔介质,按照试样颗粒组成划定不同的粒组。将不同的颗粒分组视为若干等效粒径为Ri的相互独立的独立体,假设各独立体孔隙比与试样整体的孔隙比相同。各个独立体内的孔隙均为开口孔隙,并简化为毛细管通道。经此划分和简化后,结合一定合理化假设,即可求解土水特征曲线。具体过程可归纳为两条关键步骤:①通过求各组孔隙体积Vvi与总孔隙体积Vp的比值,推算各组孔隙充满水后的体积含水率θvi;②通过各组独立体颗粒个数ni和独立体孔隙半径ri,推算各级孔隙充满水后的基质吸力ψi。
综上所述,根据AP理论,在确定好试样颗粒级配后,就可以预测试样的土水特征曲线进行。具体过程可参见文献[17]。
4.2 Arya-Paris模型关键参数与预测效果分析
基于AP模型预测SWCC的过程中,考虑到土样理想孔隙和实际孔隙结构之间的差异,引入了一个经验参数,将每个独立体中的理想球形颗粒数ni转化为更真实的。这一参数是决定AP模型预测土样土水特征曲线实际效果的关键,现对其进行探讨。
关于值,Arya[4]进行了大量测试后得出结论,值在1.31~1.43,并建议取值为1.38,适用于大多数土体。随着对AP模型越来越深入的研究发现,取恒定值并不合理。一些研究人员认为,的值应符合土颗粒粒径越大,值越大的规则。为了更加准确地预测土-水分布特征曲线,研究者提出了3种确定经验参数值的方法[10-11],即线性拟合方法、非线性拟合方法及函数法。
图 5比较了基于非线性拟合方法、线性拟合方法以及函数法所确定的试样的SWCC曲线,对比结果总结如下:①通过函数法确定经验参数的方式,最终预测得到的试样土-水特征曲线与实测值偏差最大。②对于常数,当细粒土含量较低时,预测效果相对较好,但当细粒土含量增加时,常数取值的预测曲线与试验数据偏差较大,细粒土含量越高,偏差越大。③对于线性拟合方法得到的,当细粒土含量较低时,预测效果并不理想。但当细粒土含量增加时,预测效果相对理想,细土含量越高,预测曲线越好。④对于非线性拟合方法得到的,无论细粒土含量多少,预测曲线的方法都与试验数据吻合良好。因此,建议使用该方法预测粗粒堆积土的SWCC。
5. 结论
针对目前粗粒堆积土SWCC测试设备的缺乏和SWCC研究的不足,本文首先介绍了新开发的大型非饱和粗粒堆积土力学特性测试系统,并总结了粗粒堆积土SWCC的测试方法。然后分析了不同细粒土含量的粗粒堆积土SWCC的变化规律。最后,分析了不同α测定方法的AP模型对粗粒堆积土SWCC预测的影响,得到以下3点结论。
(1)对于粗粒堆积土,建议以1 mL/h的出水速率作为更合理的吸力平衡控制标准,即每1 h排水量小于样品体积的0.04%。
(2)对于粗粒堆积土,细粒土含量越高,迫使试样排水所需的空气压力越高(进气值越高),土体的残余基质吸力越高。这一规律意味着细粒土含量越高,粗粒堆积土的保水性能越好。
(3)AP模型中关键参数的确定方法,由非线性值得到的粗粒堆积土的SWCC曲线与试验数据最为吻合。建议使用该方法,预测粗粒堆积土的SWCC曲线。
-
表 1 土层物理力学参数
Table 1 Physical and mechanical parameters of soil layer
土类 密度/(kg·m-3) 泊松比 黏聚力/kPa 内摩擦角/(°) 弹性模量/MPa 粉土 1950 0.29 16 24 217.76 表 2 文峰塔两种模型前二阶自振频率及差值
Table 2 First two natural frequencies and differences of three models for Wenfeng Pagoda
模态阶数 非倾斜模型 倾斜模型 实测频率/Hz 频率/Hz 误差e1-0/% 频率/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 注:。 -
[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.