Seismic Damage Evolution Simulation of Long-Span Deck Steel Truss Arch Bridge
-
摘要: 为评估大跨度钢拱桥在多维地震作用下的损伤状态,以一座铁路上承式钢桁架拱桥为研究对象,首先借助OpenSEES平台建立桥梁弹塑性纤维单元模型;接着选取“双参数”损伤指标和应变指标进行损伤评价;最后通过增量动力分析,从构件和全桥两个层面进行损伤评估,并对比两种指标对桥梁主要构件(钢管、钢管混凝土、钢筋混凝土)的适用性. 研究结果表明:基于变形和能量的双参数损伤指标比应变指标的评价结果会严重1~2个等级;相比于地震波纵向作用,大跨度钢桁架拱桥在横向地震作用时损伤更小,在峰值加速度(PGA)为1.5g时全桥仍处于轻微损伤状态;3种类型构件中,钢管构件损伤数量较多但是损伤程度轻,钢管混凝土构件(拱脚)损伤较早但是损伤程度轻,钢筋混凝土构件(交界墩)损伤最为严重,在抗震设计时,应对相应位置予以重视.Abstract: To evaluate the damage state of a long-span steel arch bridge under multi-dimensional earthquake loadings, a deck steel truss arch bridge was studied. A bridge elastoplastic fiber element model was established by OpenSEES code. Two damage index, i.e., strain and two-parameter damage index were selected for damage evaluation. The damage assessment of both members and the whole bridge were obtained by incremental dynamic analysis (IDA). Moreover, the applicability of two damage index on the main members of the bridge (steel pipe, concrete-filled steel tube, and reinforced concrete member) was compared. The results show that the assessment result by the two-parameter damage index based on deformation and energy may lead to 1−2 damage assessment levels worse than by strain damage index. Less damage is caused when the identical earthquake loading is input in transverse direction than longitudinal direction. The full bridge is still in a slightly damaged state when the peak ground acceleration (PGA) is 1.5g. Among three types of members, the damage quantity of steel tube members is large but damage state is slight. The damage of concrete-filled steel tubular members (arch feet) comes early but damage state is slight. The damage state of reinforced concrete members (joint balance pier) is the most serious. Attention should be paid to the corresponding structural locations in seismic design.
-
Key words:
- seismic response /
- arch bridges /
- damage index /
- incremental dynamic analysis /
- strain
-
吸力锚系统主要由吸力锚、泵机组、应急解脱系统和控制系统组成,是集机械、电气、液压、控制于一体的复合系统[1-3]. 作为一种新型的海洋工程基础,吸力锚系统的可靠性直接关系到发射平台的可靠性、稳定性以及发射任务能否顺利进行等[4]. 然而,吸力锚系统在实际工程中往往受到循环载荷作用,针对循环载荷产生的累积效应,传统的基于应力-强度干涉模型的静态可靠性方法无法准确地给予描述. 因此,考虑吸力锚系统在循环载荷作用下的累积效应,开展时变可靠性研究十分必要.
目前,国内外学者对时变可靠性分析理论方法进行了很多研究并取得了丰硕的研究成果;Rice[5]在分析随机噪声时创造性地提出了首次跨越概率计算方法,为时变可靠性分析奠定了基础;Crandall等[6]利用数值分析方法求解了跨越率;Andrieu-Renaud等[7]针对跨越率的求解提出了PHI2方法,对时变可靠性的求解提供了一种高效的计算方法;Singh等[8]针对非单调功能函数的累积失效概率计算提出了一种改进的PHI2算法;Hu等[9-10]提出了一种精度更高的联合跨越率的求解方法;Jiang等[11]研究了基于随机过程离散的时变可靠性分析方法;Li等[12]基于Gamma过程研究了数控机床刀具的时变可靠性问题;Cazuguel等[13]研究了在随机载荷作用下的平板时变可靠性模型;Kopustinskas等[14]基于随机过程的微分方程研究了动态可靠性建模方法;张德权等[15]引入区间理论处理模型中参数不确定性问题,进而提出了I-PHI2方法近似求解时变可靠性模型. 更多的时变可靠性分析理论方法参见文献[16-19].
时变可靠性目前已得到了广泛研究,但大都在理论层面,实际工程应用方面还相当欠缺. 基于此,本文以LS-FA-211001吸力锚为例,考虑循环载荷作用的累积效应,开展吸力锚的时变可靠性研究. 首先,基于广义应力-强度干涉理论,建立循环载荷作用下的时变可靠性模型;其次,结合LS-FA-211001吸力锚抗拔力和锚链拉力的不确定性量化数据,对其开展时变可靠性分析;最后,利用蒙特卡罗模拟(MCS)方法对本文结果进行对比验证.
1. 循环载荷作用下的时变可靠性模型
传统的结构可靠性模型通常采用广义应力-强度干涉模型,即在外部载荷作用下产生的应力超过强度时,结构失效,通常将应力和强度描述为随机变量,从而基于概率分布得到应力小于强度的概率,即为可靠度. 当用fs(s)、Fs(s)分别表示应力随机变量s的概率密度函数和累积分布函数,fσ(σ)、Fσ(σ)分别表示强度随机变量σ的概率密度函数和累积分布函数时,结构的可靠度表示为
R=P(s<σ)=∫+∞−∞fσ(σ)∫σ−∞fs(s)dsdσ=∫+∞−∞fσ(σ)Fs(σ)dσ. (1) 值得指出的是,式(1)是结构在单一失效模式下的可靠性模型,而且该模型描述的是结构在某一个具体时刻的可靠度,并没有表达出结构在外部载荷循环作用下的累积效应,是一个静态的可靠性模型.
当考虑外部循环载荷的累积效应时,结构的可靠性模型(即时变可靠性模型)可以表述为
R(w)=P(si<σi, ∀i=1,2,⋯,w), (2) 式中:w为外部循环载荷次数,si为外部循环载荷在第i次作用时的应力随机变量,σi为外部循环载荷在第i次作用时的强度随机变量.
此时,应力和强度表现为随作用次数变化的离散随机过程. 显然,式(2)考虑了结构在外部循环载荷作用下的累积效应,即结构在外部循环载荷作用w次的可靠性模型不再是简单地计算第w次应力小于强度的概率,而是要保证外部循环载荷在w次作用中任何一次作用应力都要小于强度.
当外部循环载荷来源于同一个随机变量,且在不考虑结构强度退化的情况下,有:
R(w)=P(si<σi, ∀i=1,2,⋯,w)=P(s1<σ,⋯,sw<σ)=P(smax<σ), (3) 式中:smax=max{s1,s2,⋯,sw},为最大顺序统计量.
由于s1,s2,⋯,sw是来源于同一个随机变量循环载荷作用下的广义应力,因此,可以把s1,s2,⋯,sw看作是来源于同一母体的样本观测值. 由概率论[20]可知,smax概率分布函数fsmax(smax)和累积分布函数Fsmax(smax)分别为
fsmax(smax)=w(Fs(s))w−1fs(s), (4) Fsmax(smax)=(Fs(s))w. (5) 此时,结合传统的广义应力-强度干涉模型可得外部循环载荷作用下的结构在单一失效模式且强度不退化下的时变可靠性模型,如式(6).
R(w)=P(smax<σ)=∫+∞−∞fσ(σ)Fsmax(σ)dσ=∫+∞−∞fσ(σ)(Fs(σ))wdσ. (6) 2. 吸力锚的时变可靠性分析
在实际工程中,LS-FA-211001吸力锚及其系统结构如图1所示,其承受着复杂的载荷,包括自身重力、土体地基的承载力、锚体侧壁与土体的摩擦力、锚链的拉力以及风、浪、流等的综合作用等. 由于条件的限制,前往锚泊现场进行水文环境监测和调研不具备可行性. 在此情形下对吸力锚展开精确的载荷分析存在困难. 因此,本文将吸力锚所受载荷进行简化,如图2所示.
在图2中,吸力锚主要受力包括:水对顶盖的压力Fw、自身重力Wp、土体的竖向承载力Fn、锚链拉力Ft以及筒壁侧摩阻力即抗拔力Tk. 在安装的沉贯阶段,吸力锚尚未连接平台,此时吸力锚在自身重力、土体的竖向承载力和水的压力作用下逐渐达到静力平衡状态. 若只考虑竖直方向载荷,则吸力锚与平台建立连接后受力可简化为锚链拉力的竖直分量Fty和抗拔力Tk.
2.1 抗拔力的不确定性量化
本文收集并计算了5个不同勘探点(点A、C、G、M、O)处的10组抗拔力数据,具体计算结果如表1所示. 此外,不同勘探点的平面图以及抗拔力计算流程分别如图3、4所示. 其中:l1、l2、l3和l4依次为某勘探点处吸力锚完全吸附时与淤泥质土、淤泥混砂、粉土-粉质黏土及粉细砂4个土层的接触厚度,L1、L2和L3分别为淤泥质土、淤泥混砂及粉土-粉质黏土3个土层的实际厚度,h为给定值,本文取值13.5.
表 1 各勘探点附近的抗拔力值Table 1. Pulling resistance values near each exploratory pointkN 组号 点A 点C 点G 点M 点O 1 1330.9777 2007.6882 2247.3961 1826.8866 1835.3776 2 1339.7254 1998.8052 2261.1515 1820.1196 1840.0185 3 1336.7282 1982.7419 2274.1373 1814.4576 1843.6222 4 1330.8861 1964.7617 2286.1875 1809.4007 1846.4546 5 1323.8135 1947.6050 2297.1344 1804.6813 1848.9262 6 1319.7456 1954.0500 2191.5069 1800.0756 1851.5450 7 1319.9079 1970.4304 2187.9786 1795.7828 1854.8203 8 1322.8437 1987.0357 2182.8988 1792.0530 1859.1009 9 1328.1747 2003.6953 2176.1155 1791.3837 1864.3240 10 1338.2259 2020.2961 2167.4843 1791.0338 1869.6512 实际工程中吸力锚的抗拔力与其载荷环境息息相关,其计算过程涉及到土层厚度等多种地勘数据,相对复杂. 为了节约篇幅,本文直接给出了5个不同勘探点处的抗拔力计算结果,其详细计算流程及载荷环境的相关原始地勘数据可以参考文献[21]. 进一步利用K-S检验,可以判断该类型吸力锚的抗拔力服从正态分布,其均值和方差具体计算为
u=1nn∑j=1xj, (7) σ2=1nn∑j=1(xj−u)2, (8) 式中:xj为用于计算某勘探点处抗拔力的样本,n为样本总量.
针对计算的不同勘探点处的原始抗拔力数据,可以估算出其对应的不确定性参数,如表2所示.
表 2 各勘探点附近抗拔力分布参数Table 2. Parameters of pulling resistance distribution near each exploratory point勘探点 u σ A 1329.1029 7.4612 C 1983.7110 24.2016 G 2227.1991 48.0464 M 1804.5875 12.7894 O 1851.3841 10.7934 基于上述估算的不确定性参数可以获得不同勘探点处抗拔力的正态分布概率密度函数和累积分布函数,分别如式(9)、(10).
fTk(x)=1√2πσe−12(x−uσ)2, (9) FTk(x)=∫x−∞1√2πσe−12(x−uσ)2dx. (10) 不同探勘点处量化的抗拔力正态分布概率密度函数如图5所示. 从图中可知,点A处的抗拔力最小,说明该勘探点处的吸力锚是最容易失效的.
2.2 吸力锚时变可靠性计算结果
本文依据LS-FA-211001吸力锚的设计指标,并结合工程经验,假设锚链拉力的竖直分量Fty服从正态分布,且其均值uS=1000,均方差σS=100,其对应的概率密度函数如图6所示.
从而,在吸力锚作用w次后,该型吸力锚的时变可靠度可以表示为
R(w)=P(Fty(w)<Tk(w))=∫+∞−∞fTk(x)(FFty(x))wdx, (11) 式中:FFty(x)为锚链拉力竖直分量Fty的累积分布函数.
最终,利用给出的时变可靠性模型,针对LS-FA-211001吸力锚在不同勘探点(点A、C、G、M、O)处的抗拔力量化数据,进行时变可靠性分析,具体计算结果如图7所示.
由图7可知:在勘探点A处,吸力锚的可靠性最差,随着使用次数的增加,其可靠度下降,再次说明了在勘探点A处的吸力锚最容易失效. 然而,即使在恶劣的勘探点A处,在可靠度需达到95%的要求下,LS-FA-211001吸力锚仍然能够使用100次左右,说明了其具有较好的环境适应性. 而在其他勘探点处的良好环境下,该型吸力锚的寿命接近于无限使用次数,说明了LS-FA-211001吸力锚具有较高的可靠性. 本文的计算结果是基于一种简化的吸力锚受力分析展开的,与实际工程中的真实情况存在一定偏差. 同时,图7中的结果仅为实际工程中决策者对吸力锚的可靠性了解及制定相应的运行维护周期提供参考.
3. 模型验证
由LS-FA-211001吸力锚的时变可靠性分析结果可知,勘探点A处的吸力锚最容易失效,时变可靠度最低. 因此,本节主要对勘探点A处的时变可靠度结果开展验证分析,以说明本文所得吸力锚时变可靠度结果的有效性. MCS方法是一种经典的数理统计方法,由于其具有较高的估计精度,因而往往用于新模型的验证. MCS方法中用于新模型可靠性结果的验证误差为
ε=|Pf−PfMCS|PfMCS×100%, (12) 式中:Pf为利用新模型评估出来的失效概率,PfMCS为采用MCS方法评估出来的失效概率.
采用MCS方法计算吸力锚在勘探点A处的时变可靠度,并与本文方法所得结果进行对比,具体结果如图8所示.
由图8可知:利用本文给出的时变可靠性模型计算的勘探点A处可靠度结果与MCS方法计算出来的结果基本一致,验证了本文方法的有效性;MCS方法与本文给出的吸力锚时变可靠性模型的最大时变可靠度计算误差为2.15%,进一步说明了本文方法的准确性,可为工程实际中对吸力锚的可靠性研究提供一定指导.
4. 结 论
本文建立了循环载荷作用下吸力锚的时变可靠性模型,利用MATLAB软件编制了对应的程序代码. 同时,收集计算了LS-FA-211001吸力锚在不同勘探点的抗拔力,并进行了不确定性量化分析. 进一步建立了不同勘探点处吸力锚应力-强度干涉模型,进行了循环载荷作用下吸力锚的时变可靠性分析,得到主要结论如下:
1) 吸力锚在勘探点A的抗拔力相对较小,在其他勘探点处的抗拔力相对差异不大,说明勘探点A处的吸力锚是最容易失效的.
2) 即使在恶劣的勘探点A,可靠度95%以上的条件下,LS-FA-211001吸力锚仍然能够使用100次左右,而在其他勘探点处基本是无限使用次数,说明了该型吸力锚可靠性高.
3) 利用MCS方法对勘探点A处的时变可靠度结果进行验证,本文方法的最大验证误差为2.15%,满足精度要求,说明了本文方法的有效性.
本文所得结果是基于一种吸力锚简化的受力分析展开的,与实际工程存在一定偏差,仅为工程实际中吸力锚的前期设计提供指导,为工程决策者提供参考.
-
表 1 材料参数
Table 1. Material parameters
项目 fpc/MPa epsc0 fpcu/MPa epsU ft/MPa Ets/MPa 下拱脚 −61.58 −0.0036 −31.24 −0.004 6.16 3079 上拱脚 −50.32 −0.0035 −15.42 −0.004 5.03 2516 交界墩横梁 −23.85 −0.0020 −4.77 −0.004 3.01 2340 交界墩冒梁 −23.40 −0.0020 −4.68 −0.004 3.01 2340 交界墩立柱 −24.28 −0.0021 −4.86 −0.004 3.01 2340 表 2 自振周期和振型
Table 2. Fundamental periods and modes
振型 频率/Hz 周期/s 振型描述 T1 0.28 3.61 拱梁一阶正对称横弯 T2 0.44 2.26 拱梁一阶反对称横弯 T5 0.91 1.10 左侧交界墩纵向弯曲 T14 1.25 0.80 拱上墩纵向弯曲 T30 2.16 0.46 左侧交界墩扭转 表 3 地震损伤等级划分
Table 3. Classification of seismic damages
损伤等级 钢管构件应变指标 钢筋混凝土构件 钢管混凝土构件钢管应变指标 钢管混凝土与钢筋混凝土构件双参数指标 混凝土应变指标 钢筋应变指标 基本完好 ε ≤ εy ε ≤ 0.0020 ε ≤ 0.00194 ε ≤ 0.00189 (0, 0.2] 轻度损伤 εy < ε ≤ 2.0εy (0.0020, 0.0035] (0.00194, 0.01000] (0.00189, 0.01000] (0.2, 0.4] 中等损伤 2.0εy < ε ≤ 8.4εy (0.0035, 0.0050] (0.01000, 0.03000] (0.01000, 0.03000] (0.4, 0.8] 严重损伤 ε > 8.4εy (0.0050, 0.0080] (0.03000, 0.05000] (0.03000, 0.05000] (0.8, 1.0] 倒塌 — ε > 0.0080 ε > 0.05000 ε > 0.05000 D > 1.0 表 4 全桥损伤指数
Table 4. Damage index of bridge
PGA/(×g) 纵向 横向 双向 损伤
指数损伤
状态损伤
指数损伤
状态损伤
指数损伤
状态0.5 0.02 − 0.07 − 0.16 − 0.6 0.11 − 0.07 − 0.25 轻微 0.7 0.13 − 0.12 − 0.34 轻微 0.8 0.21 − 0.17 − 0.34 轻微 0.9 0.28 轻微 0.17 − 0.48 中等 1.0 0.33 轻微 0.17 − 0.56 中等 1.1 0.33 轻微 0.17 − 0.58 中等 1.2 0.39 轻微 0.23 轻微 0.71 中等 1.3 0.4 轻微 0.23 轻微 0.71 中等 1.4 0.72 中等 0.37 轻微 0.73 中等 1.5 0.85 严重 0.37 轻微 0.83 严重 -
谢旭,唐站站,胡欣科,等. 纤维模型在钢拱桥抗震设计中的适用性研究[J]. 中国公路学报,2015,28(2): 33-42. doi: 10.3969/j.issn.1001-7372.2015.02.005XIE Xu, TANG Zhanzhan, HU Xinke, et al. Study on applicability of fiber model in seismic design for steel arch bridge[J]. China Journal of Highway and Transport, 2015, 28(2): 33-42. doi: 10.3969/j.issn.1001-7372.2015.02.005 LIANG C Y, CHEN A. A method for examining the seismic performance of steel arch deck bridges[J]. Frontiers of Architecture and Civil Engineering in China, 2010, 4(3): 311-320. doi: 10.1007/s11709-010-0080-8 USAMI T, GE H. A performance-based seismic design methodology for steel bridge systems[J]. Journal of Earthquake and Tsunami, 2009, 3(3): 175-193. doi: 10.1142/S179343110900055X TANG Z, XIE X, WANG T, et al. Study on FE models in elasto-plastic seismic performance evaluation of steel arch bridge[J]. Journal of Constructional Steel Research, 2015, 113: 209-220. doi: 10.1016/j.jcsr.2015.06.009 谢开仲,吕文高,覃乐勤,等. 钢管混凝土拱桥地震破坏评估研究[J]. 中国公路学报,2012,25(2): 53-59. doi: 10.3969/j.issn.1001-7372.2012.02.009XIE Kaizhong, LU Wengao, QIN Leqin, et al. Research on seismic damage evaluation of CFST arch bridges[J]. China Journal of Highway and Transport, 2012, 25(2): 53-59. doi: 10.3969/j.issn.1001-7372.2012.02.009 中国钢铁工业协会. 桥梁用结构钢: GB/T 714—2015[S]. 北京: 中国国家标准化管理委员会, 2015. 赵金钢,杜斌,占玉林,等. OpenSees中混凝土本构模型用于模拟结构滞回性能的对比[J]. 桂林理工大学学报,2017,37(1): 59-67. doi: 10.3969/j.issn.1674-9057.2017.01.008ZHAO Jingang, DU Bin, ZHAN Yulin, et al. Comparison of constitutive concrete models in OpenSees for hysteretic behavior of structures[J]. Journal of Guilin University of Technology, 2017, 37(1): 59-67. doi: 10.3969/j.issn.1674-9057.2017.01.008 陈建伟,边瑾靓,苏幼坡,等. 应用OpenSEES模拟方钢管混凝土柱的抗震性能[J]. 世界地震工程,2015,31(3): 71-77.CHEN Jianwei, BIAN Jinliang, SU Youpo, et al. Application of OpenSEES on calculating lateral force-displacement hysteretic curves of concrete-filled rectangular steel tubular columns[J]. World Earthquake Engineering, 2015, 31(3): 71-77. PARK-ANG Y J, ANG A H S. Mechanistic seismic damage model for reinforced concrete[J]. Journal of Structural Engineering, 1985, 111(4): 722-739. doi: 10.1061/(ASCE)0733-9445(1985)111:4(722) 过镇海. 钢筋混凝土原理[M]. 3版. 北京: 清华大学出版社, 2013: 312-313. 郭蓉,王铁成,赵少伟,等. 方钢管混凝土柱的地震损伤模型[J]. 河北农业大学学报,2007,30(3): 109-112. doi: 10.3969/j.issn.1000-1573.2007.03.027GUO Rong, WANG Tiecheng, ZHAO Shaowei, et al. Seismic damage model of concrete- filled rectangular tubular column[J]. Journal of Agricultural University of Hebei, 2007, 30(3): 109-112. doi: 10.3969/j.issn.1000-1573.2007.03.027 周长东,田苗旺,张许,等. 考虑多维地震作用的高耸钢筋混凝土烟囱结构易损性分析[J]. 土木工程学报,2017,50(3): 54-64.ZHOU Changdong, TIAN Miaowang, ZHANG Xu, et al. Seismic fragility analysis for high-rise RC chimney considering multi-dimensional seismic actions[J]. China Civil Engineering Journal, 2017, 50(3): 54-64. 谷音,卓卫东. 基于IDA和纤维模型的高墩大跨连续刚构桥梁地震反应分析[J]. 土木工程与管理学报,2011,28(3): 235-240. doi: 10.3969/j.issn.2095-0985.2011.03.053GU Yin, ZHOU Weidong. Study on seismic vulnerability of long-span continuous rigid frame bridge with high pier[J]. Journal of Civil Engineering and Management, 2011, 28(3): 235-240. doi: 10.3969/j.issn.2095-0985.2011.03.053 刘震. 中承式钢管混凝土拱桥地震易损性分析[J]. 公路交通科技,2015,32(8): 72-79,88. doi: 10.3969/j.issn.1002-0268.2015.08.013LIU Zhen. Analysis on seismic fragility of half-through CFST arch bridge[J]. Journal of Highway and Transportation Research and Development, 2015, 32(8): 72-79,88. doi: 10.3969/j.issn.1002-0268.2015.08.013 中国地震局. 中国地震烈度表: GB/T 17742—1999[S]. 北京: 中国国家标准化管理委员会, 2008. -