2. 北京工业大学 工程抗震与结构诊治北京市重点实验室,北京 100124
2. Beijing Key Lab. of Earthquake Eng. and Structural Retrofit,Beijing Univ. of Technol.,Beijing 100124,China
劲性支撑穹顶结构上层为索网体系,下部为刚性杆支撑,除撑杆外均为高强钢拉杆[1],该结构具有节点定位容易、施工过程相对简单、结构效率高、自支撑自平衡等优点[2]。该结构由预应力提供结构刚度使之成形,如遇不可抗力(如爆炸、冲击等),索(杆)破断将使结构发生局部破坏甚至连续倒塌,为了解结构的破坏机理及抗倒塌性能,考察局部断索(杆)后结构力学性能是十分必要的。不同位置的索(杆)发生破断,对结构的整体性能影响不同,因此对索(杆)重要性应予以关注。
劲性支撑穹顶结构是由索穹顶结构发展而来的,索穹顶结构的断索(杆)性能研究可作为其重要参考。郑君华等[3]研究了两种体系的索穹顶结构部分索松弛或破断的受力性能。黄华等[4]分析了轮辐式索膜结构局部索破坏对结构性能的影响。何键等[5]分析了索穹顶结构局部断索杆件内力的变化和结构最终到达的平衡位置,并判断结构体系是否安全以及不同位置杆件的安全等级。宗钟凌等[6]对一直径6 m的葵花型索穹顶试验模型进行了单根拉索破断时结构的响应特性进行了破断试验研究。赵冉等[7]对深圳宝安体育场屋盖索膜结构的索进行了破断分析。蔡建国等[8]针对广州站索拱结构屋盖,采用变换荷载路径法对其进行了连续倒塌分析。朱明亮等[9]考虑结构初始状态,采用等效荷载卸载法分析得到了新型环箍–穹顶全张力结构不同类型拉索破断后剩余结构的动力响应,并对其进行了动力响应分析。王功文等[10]提出了基于结构自振频率的构件重要性系数计算方法,以自振频率的下降作为构件重要性的评价指标,对索网结构、索穹顶结构等进行了分析。杨逢春[11]以概率易损性的结构响应作为性能标准,通过拆除杆件后结构体系可靠度退化程度评价杆件的重要性。
作者模拟劲性支撑穹顶结构的单根索(杆)破断的过程,通过响应变化率和动力响应系数分析结构单根索(杆)破断后的力学性能及对结构整体性能的影响。判断结构在索(杆)破断之后是否发生连续倒塌,对杆件的敏感性进行分类。通过杆件重要性系数对结构的杆件重要性进行分析,并对杆件的重要性进行排序。
1 基本理论劲性支撑穹顶结构由预应力提供刚度,索(杆)的破断是复杂突发破坏事件,属瞬态动力学问题。基于瞬态动力学分析原理,采用Newmark时间积分法进行分析,该方法没有矩阵缩减,可分析大变形等各类非线性特征问题。利用ANSYS进行模拟,用Link10单元模拟结构中的索,Link8单元模拟结构中的高强钢拉杆和撑杆,通过单元生死技术模拟杆件破断的瞬间,对于因单元杀死而可能出现“漂移”(位移过大导致计算不收敛)的节点,将其进行约束。采用自动时间步长,最小子步数取20。
为了反应劲性支撑穹顶结构断索(杆)后结构的力学性能和杆件的敏感性,定义响应变化率ρ和动力响应系数R两个指标。响应变化率是结构振动平稳后的响应量和正常工作下的响应量之差与正常工作下响应量的比值,如式(1)所示。动力响应系数是振动过程中最大的响应量与正常工作下响应量的比值,如式(2)所示。
$\rho = \frac{{{S_{{\text{稳}}}} - {S_{{\text{静}}}}}}{{{S_{{\text{静}}}}}}$ | (1) |
$R = \displaystyle\frac{{{S_{{\text{最大}}}}}}{{{S_{{\text{静}}}}}}$ | (2) |
式中:ρ为响应变化率;R为动力响应系数;S稳为振动稳定后的响应量;S静为正常工作下的响应量;S最大为振动过程中最大的响应量。
结构失效即当一根杆件断裂后,结构变成几何可变体系或其他杆件超过屈服强度而发生连续倒塌。按结构是否失效及影响范围的大小,将结构中的杆件分为3类:全局敏感性构件、局部敏感性构件、不敏感性构件。全局敏感性构件即内力的变化对整个结构的受力产生影响,结构不能继续承载,发生连续倒塌。局部敏感性构件即内力的变化只对其附近的构件造成影响,不会影响整个结构的受力。不敏感性构件即内力变化对其他构件的内力影响很小,可忽略不计[12]。
2 静力分析设2道环杆的劲性支撑穹顶结构,直径为70 m,矢高5.5 m,环向20等分,如图1所示,图1(a)为平面图,图1(b)为剖面图(图1中JS1为外脊索,JS2为中脊索,JS3为内脊索,XG1为外斜杆,XG2为中斜杆,XG3为内斜杆,HG1为外环杆,HG2为内环杆,CG1为外撑杆,CG2为中撑杆,CG3为内撑杆)。采用1 670 MPa级钢丝索,其弹性模量为1.85×105 MPa。高强钢拉杆为650 MPa级,其弹性模量为2.06×105 MPa。撑杆为305 MPa级圆钢管。周边三向固定铰支承,表面覆膜材,重度为0.01 kN/m3,屋面活荷载为0.5 kN/m2。依据快捷确定初始预应力的方法[13]确定初始预应力,见表1。
![]() |
图1 劲性支撑穹顶结构 Fig. 1 Rigid bracing dome |
表1 杆件的截面类型和初始预应力 Tab. 1 Section type and initial prestress of members |
![]() |
杆件编号的顺序为脊索、斜杆、环杆及撑杆,以1号轴对应的杆件为起点,同类杆件按先外后内逆时针顺序编制。节点编号为1号轴上的最外圈节点编为1号节点,按照逆时针顺序编写,上层和下层节点编号也如此编写,如图1(a)所示。断索(杆)均选取在1号轴上的相应杆件。
结构在静力荷载作用下的变形和杆件应力见图2。
![]() |
图2 静载作用下结构的变形图和应力图 Fig. 2 Deformation and stress of structure under static load |
3 断索(或杆)分析
劲性支撑穹顶结构的冗余度低,局部断索(或杆)可能导致结构的几何刚度大幅减小,结构发生较大的局部位移或连续倒塌。对劲性支撑穹顶结构单根构件破断后的体系进行研究,考察应力和位移的响应变化率和响应系数反应结构性能和结构抗连续倒塌能力。
1)外脊索
图3给出了1号外脊索破断后杆件应力和节点竖向位移的响应变化率图和动力响应系数图。
其中,用极坐标表述,极轴表示响应变化率或动力响应系数,角度表示与1号轴逆时针方向的夹角(即0度对应的是1号轴),杆件编号方法如前所述,图中下角标A、B、C、D、E、F分别代表外圈撑杆上节点、中圈撑杆上节点、内圈撑杆上节点、外圈撑杆下节点、中圈撑杆下节点、内圈撑杆下节点,反应索破断对整个结构构件的影响。
由图3(a)和(b)可知,21号中脊索松弛,81号中斜杆和101号内斜杆出现压应力,值分别为30.11和33.85 MPa,但未失稳,82号和100号中斜杆、102号和120号内斜杆应力响应变化率均为正值,振动稳定后应力增大,但未超过钢拉杆的抗拉强度。其他杆件稳定后的应力均减小。1号外脊索两侧外脊索和2号轴及20号轴对应的中脊索、内脊索、中斜杆、内斜杆、外撑杆及中撑杆的应力动力响应系数均大于1,但均未超过杆件的抗拉强度。由图3(c)和(d)可知,21号节点的竖向位移响应变化率为–3.364 7,竖向位移动力响应系数为–2.542,即出现向上位移,而两侧相邻节点的竖向位移响应变化率均为正值,说明产生比静载时更大的向下位移。内侧两圈节点的竖向位移响应变化率均大于零,比静载时更大的竖向位移。
![]() |
图3 1号外脊索破断 Fig. 3 No. 1 lateral upper cable is broken |
1号外脊索破断,只对断索位置附近的杆件有较大影响,其他位置的杆件影响程度均较小。结构并未此发生连续倒塌,故外脊索为局部敏感性构件。
2)中脊索
图4给出了21号中脊索破断后杆件应力和节点竖向位移的响应变化率图和动力响应系数图。
由图4(a)和(b)可知,21号中脊索破断后,41号内脊索松弛,101号内斜杆出现压应力,值为0.85 MPa,但并未发生失稳。破断索两侧的外脊索、中脊索、内脊索、内斜杆及中撑杆的应力响应变化率均为正值,且应力动力响应系数大于1,说明结构振动稳定后及振动过程中杆件内力均增大,但未超过杆件的抗拉强度。其他杆件的应力响应变化率为负,应力动力响应系数小于1,但降低的有限。由图4(c)和(d)可知,81号节点竖向位移响应变化率小于–1,说明结构振动平稳后出现向上的位移。
![]() |
图4 21号中脊索破断 Fig. 4 No. 21 middle upper cable is broken |
21号中脊索破断,只对1号轴、2号轴及20号轴上的杆件有较大影响,对其它杆件影响不大,未引起整个结构的连续倒塌。因此,中脊索为局部敏感性构件。
3)内脊索
图5给出了41号内脊索破断后杆件应力和节点竖向位移的响应变化率图和动力响应系数图。
由图5(a)和(b)可知,21号中脊索出现松弛,101号内斜杆应力响应变化率大于零,说明结构振动稳定后杆件内力增加,但未超过杆件的抗拉强度。2号轴和20号轴两侧的脊索、1号轴对应的内斜杆和中撑杆的应力动力响应系数大于1,说明振动过程中杆件的最大应力增加,但未超过杆件的抗拉强度。由图5(c)和(d)可知,21号、41号、81号、101号节点的竖向位移响应变化率大于零,竖向位移动力响应系数大于1,说明节点竖向位移比静载时大。
![]() |
图5 41号内脊索破断 Fig. 5 No. 41 internal upper cable is broken |
41号内脊索破断,只对1号、2号及20号轴上的杆件及节点有较大影响,对其他杆件及节点影响不大,并未引起整个结构的连续倒塌。因此,内脊索为局部敏感性构件。
4)外斜杆
图6为61号外环杆破断后结构位移图,结构最大位移为1.383 m,大部分脊索发生松弛,结构发生完全倒塌,因此外斜杆为全局敏感性构件。
![]() |
图6 61号外环杆破断时结构的位移图 Fig. 6 After breaking outer ring rod, displacement of the structure |
5)中斜杆
图7给出了81号中斜杆破断后杆件应力和节点竖向位移的响应变化率图和动力响应系数图。
![]() |
图7 81号中斜杆破断 Fig. 7 No.81 outer oblique rod is broken |
由图7(a)和(b)可知,41号内脊索、162号和180号外撑杆、182号和200号中撑杆的应力响应变化率大于零,应力动力响应系数大于1,说明结构振动稳定后杆件的内力增加,但未超过杆件的抗拉强度。101号内斜杆的应力响应变化率小于–1,说明内斜杆出现压应力,但并未失稳。由图7(c)和(d)可知,21号和81号节点竖向位移响应变化率分别为–3.031 7和–3.056 6,振动稳定后出现向上的竖向位移,4号轴到18号轴对应的外圈上下节点的竖向位移较静力荷载作用下有所下降。
81号中斜杆发生断裂时,结构只在断杆位置附近发生较大的振动,其他位置的节点和杆件振动幅度均较小。中斜杆断裂以后,局部发生破坏,但振动稳定后整个结构安全,因此,中斜杆为局部敏感性构件。
6)内斜杆
图8给出了101号内斜杆破断后杆件应力和节点竖向位移的响应变化率图和动力响应系数图。
![]() |
图8 101号内斜杆破断 Fig. 8 No.101 internal oblique rod is broken |
图8(a)和(b)可知,41号内脊索、102号和120号内斜杆应力响应变化率大于零,杆件拉应力增加。相应杆件的应力动力响应系数分别为1.372 6和1.371 8,但均未超过杆件的抗拉强度。182号和200号中撑杆的应力响应变化率也大于零,应力动力响应系数为1.218 1,撑杆压应力增加,但并未因此发生失稳。图8(c)和(d)可知,41号节点和101号节点的竖向位移响应变化率分别为–1.070 7和–1.067 8,均小于–1,说明振动平稳后杆件出现向上的位移。
内斜杆发生断裂时,结构只在断索位置附近发生较大的振动,其它位置的节点位移和杆件内力影响较小。振动稳定后整个结构安全,未发生连续倒塌,内斜杆为局部敏感性构件。
7)外环杆、内环杆、外撑杆
图9分别给出了外环杆、内环杆、外撑杆破断后结构的位移图。
![]() |
图9 结构的位移图 Fig. 9 Structure displacement diagram |
121号外环杆或141号内环杆或161号外撑杆发生破断,大部分脊索内力为零,出现过大位移,结构发生连续倒塌,因此外环杆、内环杆和外撑杆均为全局敏感性构件。
8)内撑杆
图10给出了181号内撑杆破断后杆件应力和节点竖向位移的响应变化率图和动力响应系数图。
图10(a)和(b)可知,41号内脊索、102号和120号内斜杆应力响应变化率分别为0.674 1和0.689 9,振动稳定后脊索内力增加,但未达到抗拉强度。101号内斜杆应力响应变化率为–1.6462,小于–1,说明内斜杆振动稳定后出现压应力,值为31.22 MPa,未发生屈服。图10(c)和(d)可知,41号节点竖向位移响应变化率为0.323 6,竖向位移动力响应系数为1.326 8,稳定后节点竖向位移较静载时大很多。101号节点竖向位移响应变化率为–0.402 3,稳定后节点上抬幅度最大。
![]() |
图10 181号内撑杆破断 Fig. 10 No.181 outer support rod is broken |
内撑杆发生破断后,结构只在断索位置附近发生较大的振动,其它位置的振动幅度均较小。内撑杆断裂振动稳定后整个结构安全,内撑杆为局部敏感性构件。
4 杆件重要性分析构件重要性反应的是结构在突发事件中,个别构件的损伤或破坏对结构整体性能影响程度的大小。结构中某些构件的重要性指标越突出,则结构整体的鲁棒性越低,结构中存在易于被袭击的致命缺陷。构件重要性的定量分析是结构鲁棒性定量研究和结构抗连续倒塌性能研究的基础。上述研究可分析破断索(或杆)后结构的力学性能和结构的敏感性,但无法全面反映杆件在整体结构中的重要性及各杆件的重要性顺序。引入杆件重要性系数
${\beta _j} = \frac{{\displaystyle\sum\limits_{i = 1}^n {(| {{\rho_{xi}}} | +| {{\rho _{yi}}} | + | {{\rho_{{\text z}i}}} |)} }}{{3n}}$ | (3) |
式中,n为节点数,j为移除的杆件,i为节点,ρxi、ρvi、ρzi分别为节点i在直角坐标系下x、y、z方向的位移响应变化率(或Rxi、Rvi、Rzi的动力响应系数)。
杆件重要性系数越大说明杆件越重要,反之则越不重要,因此可通过杆件重要性系数确定杆件的重要性并排序。
应用式(1)计算各节点3个方向的位移响应变化率(或动力响应系数),然后代入公式(3)计算杆件重要性系数,并依据杆件重要性系数的大小,对杆件的重要性进行排序,如表2所示。
表2 杆件重要性系数及排序 Tab. 2 Importance coefficients and ordering of members |
![]() |
由表2可知,环杆的重要性尤为重要,依据构件重要性系数将构件的重要性进行排序为:外环杆>内环杆>外斜杆>外撑杆>外脊索>中斜杆>中脊索>内撑杆>内斜杆>内脊索。由此可见,全局敏感性构件的重要性排在前面,其次是局部敏感性构件,最后为不敏感性构件。
5 结 论模拟了经性支撑穹顶结构的断索(杆)过程,给出了响应变化率和动力响应系数两个指标,分析了结构单根索(杆)破断后的力学性能,考察其抗连续倒塌性能,将杆件进行敏感性分类。通过杆件的重要性系数,对杆件的重要性进行分析并排序。通过研究得到以下结论:
1)提出的响应变化率和动力响应系数能充分反映断索(杆)后结构性能,并通过其判断索(杆)的敏感性。
2)外环杆、内环杆、外斜杆和外撑杆破断后结构多数索出现松弛,结构出现过大位移,发生连续倒塌,其为全局敏感性构件。
3)中斜杆、内斜杆、外脊索、中斜索、内脊索或内撑杆破断后,影响在断索(或杆)附近,其他位置的杆件影响均较小,这些杆件为局部敏感性杆件。
4)外脊索和中脊索破断后,会导致相应位置的部分斜杆出现压应力,但并未发生失稳。外脊索、中脊索和内脊索破断都会引起相应部位的索出现松弛。
5)杆件重要性系数能反映杆件对整体结构的重要性,通过其对杆件重要性进行排序。环杆的重要性尤为重要。外环杆的杆件重要性排在首位,内脊索的重要性排在末位。外侧杆件比相应的内侧杆件的杆件重要。全局敏感性杆件重要性排序较高,局部敏感性的杆件重要性居中,不敏感性构件重要性较低。通过杆件重要性分析,可为进一步连续倒塌分析提供依据。
[1] |
Gao Zhanyuan, Xue Suduo, He Yongfa. Analysis of design and construction integration of rigid bracing dome[J]. Advances in Structural Engineering, 2015, 18(11): 1947-1958. DOI:10.1260/1369-4332.18.11.1947 |
[2] |
Gao Zhanyuan, Xue Suduo, Li Xiongyan. Form-finding analysis of rigid bracing dome[J]. Journal of Sichuan university (Engineering Science Edition), 2014, 46(1): 87-94. [高占远, 薛素铎, 李雄彦. 劲性支撑穹顶结构找形分析[J]. 四川大学学报(工程科学版), 2014, 46(1): 87-94.] |
[3] |
Zheng Junhua, Yun Xingfei, Dong Shilin. Failure mode and behaveior analysis of two types of cable domes[J]. Engineering Mechanics, 2007, 24(1): 44-50. [郑君华, 袁行飞, 董石麟. 两种体系索穹顶结构的破坏形式及其受力性能研究[J]. 工程力学, 2007, 24(1): 44-50.] |
[4] |
Huang Hua, Xiao Yaoqiang, Liu Boquan. Progressive collapse property of wheel-spoke cable-membrane structure[J]. Journal of Vibration and Shock, 2015, 34(20): 27-36. [黄华, 冼耀强, 刘伯权. 轮辐式索膜结构连续倒塌性能分析[J]. 振动与冲击, 2015, 34(20): 27-36.] |
[5] |
He Jian, Yuan Xingfei, Jin Bo. Analysis of cable domes with rupture of local cable[J]. Journal of Vibration and Shock, 2010, 29(11): 13-16. [何键, 袁行飞, 金波. 索穹顶结构局部断索分析[J]. 振动与冲击, 2010, 29(11): 13-16. DOI:10.3969/j.issn.1000-3835.2010.11.003] |
[6] |
Zong Zhongling, Guo Zhengxing. Experimental research on mechanical properties and cable broken of levy cable dome[J]. Engineering Mechanics, 2013, 30(1): 271-276. [宗钟凌, 郭正兴. 葵花型索穹顶结构力学性能及拉索破断试验研究[J]. 工程力学, 2013, 30(1): 271-276.] |
[7] |
Zhao Ran,Wei Demin,Sun Wenbo,et al.Form-finding and cable-breaking analyses for cable-membrane roof structure of Shenzhen Baoan stadium[J].Engineering Mechanics,2010,27(Supp 1):266–269. [赵冉,魏德敏,孙文波,等.深圳宝安体育场屋盖索膜结构的找形和索的破断分析[J].工程力学,2010,27(增1):266–269.]
|
[8] |
Cai Jianguo, Wang Fenglan, Feng Jian. Progressive collapse analysis of cable-arch structures of the New Guangzhou Railway Station[J]. Journal of Building Structures, 2010, 31(7): 103-109. [蔡建国, 王蜂岚, 冯健. 新广州站索拱结构屋盖体系连续倒塌分析[J]. 建筑结构学报, 2010, 31(7): 103-109.] |
[9] |
Zhu Mingliang, Lu Jinyu, Guo Zhengxing. Analysis on progressive collapse resistance of novel torus-dome tensile structure due to cable rupture[J]. Journal of Southeast University (Natural Science Edition), 2016, 46(5): 1057-1062. [朱明亮, 陆金钰, 郭正兴. 新型环箍–穹顶全张力结构局部断索抗连续倒塌性能分析[J]. 东南大学学报(自然科学版), 2016, 46(5): 1057-1062. DOI:10.3969/j.issn.1001-0505.2016.05.027] |
[10] |
Wang Gongwen.A frequency-based analytical method of component importance of cable-strut tensile structures[D].Guangzhou:South China University of Technology,2015. [王功文.基于频率的索杆张力结构构件重要性分析方法[D].广州:华南理工大学,2015.]
|
[11] |
Yang Fengchun.Component importance analysis of spatial structures based on structural system reliability[D]. Hangzhou:Zhejiang University,2015. [杨逢春.基于结构体系可靠度的空间结构杆件重要性分析[D].杭州:浙江大学,2015.]
|
[12] |
郭彦林,田广宇.索结构体系设计原理与施工控制[M].北京:科学出版社,2014.
|
[13] |
Dong Shilin, Yuan Xingfei. A quick calculation method for initial prestress distribution of geiger domes[J]. Spatial Structures, 2003, 9(2): 3-8. [董石麟, 袁行飞. 肋环型索穹顶初始预应力分布的快速计算法[J]. 空间结构, 2003, 9(2): 3-8.] |