Influence of Bi-directional Horizontal Ground Motion on the Vulnerability of a Reinforced Concrete Continuous Beam Bridge
-
摘要: 以双向水平地震动作为输入,对钢筋混凝土连续箱梁高架桥开展非线性动力时程分析。建立基于双向水平地震动强度参数的桥梁结构易损性曲面,比较单向及双向水平地震动输入下桥梁结构易损性差异,分析双向水平地震动输入下横桥向地震动强度对桥梁整体易损性的影响规律。研究结果表明,双向水平地震动输入下的桥梁结构易损性明显高于单向地震动输入的情况,且随着横桥向输入地震动强度的增加,结构各破坏状态的超越概率明显增大。Abstract: The nonlinear dynamic time history analysis of reinforced concrete continuous box girder bridge is carried out with bidirectional horizontal ground motion as input. The vulnerability surface of bridge structure is established based on bidirectional horizontal ground motion intensity parameters. The difference of vulnerability of bridge structure under unidirectional and bidirectional horizontal ground motion input is compared, and the influence law of transverse seismic intensity on overall vulnerability of bridge under bidirectional horizontal ground motion input is analyzed. The results show that: the vulnerability of bridge structure under bidirectional horizontal ground motion input is significantly higher than that under unidirectional ground motion input, and with the increase of transverse ground motion input intensity, the exceeding probability of each failure state of the structure also increases significantly.
-
引言
地震前兆观测主要关注观测数据随时间的相对变化,装置系数误差不影响观测数据的相对变化,但不正确的装置系数可能导致地电阻率观测结果出现系统误差(王兰炜等,2014),因此,正确的装置系数有利于不同区域观测数据的对比和地震前兆数据的研究。自2009年起,河北大柏舍台,甘肃天水台、武都台、平凉台,陕西合阳台实施了井下地电阻率垂直观测试验,井孔深100—225m,供电极距60—120m,测量极距20—60m(刘君等,2015;王兰炜等,2015)。上述台站地电阻率垂向观测通常为1个钻孔,供电电极和测量电极均布设于1个钻孔中,部分垂向观测的供电电极A接近地表,如天水台、武都台、合阳台的供电电极A埋深仅4—5m;部分垂向观测的供电电极A埋深为40m左右,如大柏舍台。垂向地电阻率观测中的装置系数与电流的空间分布及电极位置有关,现有垂向观测装置系数计算方法依据地下点、地表点电流源产生的电场计算得出,忽略了供电电极A的埋深。江宁台深井垂向地电阻率观测装置为在2口深井进行垂向观测的试验装置,与传统垂向地电阻率观测装置不同。本文根据地下点电流源产生的电场讨论装置系数计算方法,并比较计算方法对江宁台垂向地电阻率观测的影响。
1. 台址概况
江宁台地处南京市江宁区禄口街道水荆墅村,地形开阔平坦,周围无大中型工矿企业,测区位于南京-湖熟断裂南西盘和方山-小丹阳断裂西盘的楔形地块上,东距茅山断裂带30km,西北距长江36km。测区内现有垂向地电阻率观测装置为在2口深井进行垂向观测的试验装置,井距5.17m,供电电极A、测量电极M分别布设在深275m的2号井内200m和275m处,供电电极B、测量电极N分别布设在深400m的1号井内400m和325m处(图 1)。该垂向观测系统采用ZD8BI型地电仪,根据《地震台站建设规范(地电台站第1部分)》(DB/T.18.1—2006)中关于地电阻率台站的技术要求,对新建垂向观测装置的场地进行高密度电法、电测深等测试。高密度电法探测和电测深报告中NW—SE和NS测线结果表明,观测区域电测深曲线具有K形特征,电性结构等效为3层(樊晓春等,2018),场地电性层参数见表 1。
表 1 江宁台场地电性层参数Table 1. The underground electrical structure of Jiangning geoelectric stationNW—SE测线 NS测线 层厚/m 电阻率/Ω·m 层厚/m 电阻率/Ω·m 29.40 39.29 24.56 18.71 220.94 143.06 203.42 274.52 ∞ 60.39 ∞ 54.79 2. 点电流源在非全空间产生的电场
点电流处于不完全全空间时,对点电流源位于地下和地表 2种情况进行讨论。地下点电流源产生的电场指点电流源的电流I在地下一定深度时流入地下介质中产生的电场,为不完全全空间。假设地下介质电性均匀,介质电阻率为ρ,电流I从地下A点流入(图 2),采用镜像法计算(刘昌谋等,1994;刘国兴,2005),见式(1)。地表点电流源产生的电场指地表点电流源电流I流入地下介质,从无限远处流出时在介质中产生的电场,为半无限空间。假设地下介质电性均匀,介质电阻率为ρ,电流I从地表A点流入地下,电流线的分布以A为中心向周围呈辐射状,该情况为图 2的特例(王兰炜等,2014),见式(2)。
(1)地下点电流源产生的电场(不完全全空间)
$$ {V_{{\rm{MN}}}} = \frac{{\rho I}}{{4{\rm{ \mathsf{ π} }}}}\left({\frac{1}{{\overline {AM} }} + \frac{1}{{\overline {{A_1}M} }} - \frac{1}{{\overline {{A_1}N} }} - \frac{1}{{\overline {AN} }}} \right) $$ (1) (2)地表点电流源产生的电场(半无限空间)
$$ {V_{{\rm{MN}}}} = \frac{{\rho I}}{{2{\rm{ \mathsf{ π} }}}}\left({\frac{1}{{\overline {AM} }} - \frac{1}{{\overline {AN} }}} \right) $$ (2) 3. 装置系数计算方法
装置系数是地电阻率观测中特有的参数,通常用K表示,与观测装置中电极分布情况有关,表征地电阻率是视电阻率(王兰炜等,2014)。当测区区域介质电阻率均匀分布时,地电阻率ρs与介质真实的电阻率ρ数值相同。
根据奥斯定理和镜像原理(钱家栋等,1985),假设测区为均匀介质,垂向观测中的装置系数计算方法如下:
(1)方法Ⅰ:传统垂向观测装置系数计算方法(王兰炜等,2014)
传统方法井下垂向观测装置忽略电极A的埋深,在点电源B与地面对称的位置设镜像点B1,见图 3(a)。根据式(1)和式(2),则K为:
$$ K = \frac{{4{\rm{ \mathsf{ π} }}}}{{\left({\frac{2}{{\overline {AM} }} - \frac{2}{{\overline {AN} }}} \right) - \left({\frac{1}{{\overline {BM} }} + \frac{1}{{\overline {{B_1}M} }} - \frac{1}{{\overline {{B_1}N} }} - \frac{1}{{\overline {BN} }}} \right)}} $$ (3) 不考虑江宁台垂向观测电极A埋深时,因江宁台垂向观测AM=BN,AN=BM,则:
$$ K = \frac{{4{\rm{ \mathsf{ π} }}}}{{\frac{3}{{\overline {AM} }} - \frac{3}{{\overline {AN} }} - \frac{1}{{\overline {{B_1}M} }} + \frac{1}{{\overline {{B_1}N} }}}} $$ (4) (2)方法Ⅱ:采用全空间方式的装置系数计算方法(钱家栋等,1985;王兰炜等,2014)
当电极埋深h远大于供电极长度AB时,为全空间,则K为地表观测装置系数的2倍,即:
$$ K = \frac{{4{\rm{ \mathsf{ π} }}}}{{\frac{1}{{\overline {AM} }} - \frac{1}{{\overline {AN} }} - \frac{1}{{\overline {BM} }} + \frac{1}{{\overline {BN} }}}} $$ (5) (3)方法Ⅲ:采用不完全全空间方式的装置系数计算方法
江宁台垂向观测的电极A、M、N、B分别位于埋深200m、275m、325m、400m处,应按地下点电源产生的电场模型计算(不完全全空间),如图 3(b)所示。在点电源A、B与地面对称的位置设镜像点A1、B1,忽略1号井和2号井的水平距离l,根据式(1),则供电电流I(+I和-I)在M、N间产生的电位差为:
$$ {V_{MN}} = \frac{{\rho I}}{{4{\rm{ \mathsf{ π} }}}}\left[ {\left({\frac{1}{{\overline {AM} }} + \frac{1}{{\overline {{A_1}M} }} - \frac{1}{{\overline {{A_1}N} }} - \frac{1}{{\overline {AN} }}} \right) - \left({\frac{1}{{\overline {BM} }} + \frac{1}{{\overline {{B_1}M} }} - \frac{1}{{\overline {{B_1}N} }} - \frac{1}{{\overline {BN} }}} \right)} \right] $$ (6) 则K为:
$$ K = \frac{{4{\rm{ \mathsf{ π} }}}}{{\left({\frac{1}{{\overline {AM} }} + \frac{1}{{\overline {{A_1}M} }} - \frac{1}{{\overline {AN} }} - \frac{1}{{\overline {{A_1}N} }}} \right) - \left({\frac{1}{{\overline {BM} }} + \frac{1}{{\overline {{B_1}M} }} - \frac{1}{{\overline {BN} }} - \frac{1}{{\overline {{B_1}N} }}} \right)}} $$ (7) 因江宁台垂向观测AM=BN,AN=BM,则:
$$ K = \frac{{4{\rm{ \mathsf{ π} }}}}{{\frac{2}{{\overline {AM} }} - \frac{2}{{\overline {AN} }} + \frac{1}{{\overline {{A_1}M} }} - \frac{1}{{\overline {{A_1}N} }} - \frac{1}{{\overline {{B_1}M} }} + \frac{1}{{\overline {{B_1}N} }}}} $$ (8) (4)方法Ⅳ:采用不完全全空间方式(考虑井距l)的装置系数计算方法
按照地下点电源产生电场的模型计算(不完全全空间),在点电源A、B与地面对称的位置设镜像点A1、B1,h1、h2、h3表示供电电极A、B、M的电极埋深,井距l表示2口井孔水平距离(图 3(c)),则:
$$ \overline {AM} = {h_3} - {h_1} $$ (9) $$ \overline {AN} = \sqrt {{l^2} + {{({h_2} - {h_3})}^2}} $$ (10) $$ \overline {{A_1}M} = {h_1} + {h_3} $$ (11) $$ \overline {{A_1}N} = \sqrt {{l^2} + {{({h_2} + 2{h_1} - {h_3})}^2}} $$ (12) $$ \overline {{B_1}M} = \sqrt {{l^2} + {{({h_2} + {h_3})}^2}} $$ (13) $$ \overline {{B_1}N} = 2{h_2} - {h_3} + {h_1} $$ (14) 将式(9)至式(14)代入式(8),则K变为:
$$ K = \frac{{4\pi }}{{\frac{2}{{{h_3} - {h_1}}} - \frac{2}{{\sqrt {{l^2} + {{({h_2} - {h_3})}^2}} }} + \frac{1}{{{h_1} + {h_3}}} - \frac{1}{{\sqrt {{l^2} + {{({h_2} + 2{h_1} - {h_3})}^2}} }} - \frac{1}{{\sqrt {{l^2} + {{({h_2} + {h_3})}^2}} }} + \frac{1}{{2{h_2} - {h_3} + {h_1}}}}} $$ (15) 4. 江宁台垂向地电阻率观测装置系数计算
采用方法Ⅰ—Ⅳ分别计算江宁台垂向地电阻率观测的装置系数K,结果见表 2。考虑江宁台垂向地电阻率观测的电极布设不同于传统垂向观测装置,所有电极埋深均在200m以上,本文认为应以方法Ⅳ为参考值,采用式(16)计算不同装置系数计算方法的相对误差。方法Ⅰ、Ⅱ、Ⅲ相对误差分别为-32.01%、1.37%、0.43%,可知方法Ⅰ不适用于江宁台垂向地电阻率观测,该计算方法通常仅适用于供电电极A埋深小于5m的垂向观测,如天水台、合阳台。方法Ⅱ相对误差较小,江宁台垂向观测于2018年4月25日至2018年11月22日曾采用该方法。考虑仅当电极埋深远大于供电极距时称为全空间,而江宁台垂向装置最小电极埋深仅与供电极距相当,因此,方法Ⅱ同样不适用于江宁台垂向观测。除江宁台外,大部分台站观测装置电极埋深明显小于供电极距,均不宜采用方法Ⅱ。方法Ⅲ相对误差最小,江宁台垂向地电阻率观测于2018年11月23日至2019年10月30日曾采用该方法,2018年11月出现的台阶是由调整装置系数导致的(图 4)。由于江宁台垂向地电阻率观测为多孔观测,须考虑井距才能准确计算其装置系数,因而最终采用方法Ⅳ进行计算。
$$ \sigma = \frac{{k_{方法}} - {k_{方法Ⅳ}}}{{{k_{方法Ⅳ}}}} $$ (16) 表 2 江宁台垂向观测的装置系数Table 2. The configuration coefficient of vertical geo-resistivity observation in Jiangning earthquake station计算方法 AM/m AN/m B1M/m B1N/m A1M/m A1N/m 装置系数K Ⅰ 75.000 125.000 675.000 725.000 — — 790.000 Ⅱ 75.000 125.000 — — — — 1178.000 Ⅲ 75.000 125.000 675.000 725.000 475.000 525.000 1167.000 Ⅳ 75.000 125.374 675.020 725.000 475.000 525.292 1162.000 5. 结论
本文以江宁台垂向地电阻率观测为例,提出2种以不完全全空间方式计算的新方法,并与现有垂向观测装置计算方法进行比较。研究结果表明,方法Ⅳ最符合江宁台垂向地电阻率观测装置。考虑方法Ⅳ中部分参数存在小数部分,认为保留小数点后三位能满足装置系数精度要求。方法Ⅳ除适用于2口井垂向观测装置外,同样适用于忽略井距时供电电极A埋深超过5m的单口井垂向观测装置。
致谢: 衷心感谢中国地震局地壳应力研究所王兰炜研究员对本文提出的建议和意见。 -
表 1 支座参数
Table 1. Parameter of bearings
支座类型 支座型号 支座位置 盆式固定支座 GPZ4000GX 3 盆式单向滑动支座 GPZ4000DX 1、2、4、5(上) 盆式双向滑动支座 GPZ4000SX 1、2、4、5(下) 表 2 损伤破坏准则
Table 2. Damage failure criterion
损伤状态 基本完好 轻微破坏 中等破坏 严重破坏 完全破坏 破坏准则 0<μ≤μcy1 μcy1<μ≤μcy μcy<μ≤μc4 μc4<μ≤μcmax μcmax≤μ 表 3 桥墩截面弯矩-曲率关系(单位:rad·m−1)
Table 3. Moment- curvature relationship of pier section(Unit: rad·m−1)
分析参数 фy1 фy фc4 фu 桥墩截面曲率 2.353 3.045 18.390 50.630 表 4 损伤指标与损伤等级关系
Table 4. Relationship between damage index and damage grade
损伤状态 基本完好 轻微破坏 中等破坏 严重破坏 完全破坏 破坏准则 0<μ≤1.00 1.00<μ≤1.29 1.29<μ≤3.34 3.34<μ≤6.34 6.34<μ 表 5 不同地震动强度参数的相关系数
Table 5. Determination coefficients of different ground motion intensity parameters
强度参数 Y=0 g Y=0.2 g Y=0.4 g Y=0.6 g Y=0.8 g Y=1.0 g Sa(T1) 0.837 5 0.819 4 0.763 8 0.716 5 0.683 1 0.657 5 PGA 0.485 7 0.454 0 0.419 8 0.373 6 0.336 9 0.301 7 PGV 0.813 5 0.780 1 0.714 6 0.659 5 0.611 8 0.572 2 PGD 0.813 7 0.809 1 0.769 8 0.739 1 0.714 0 0.693 9 表 6 不同地震动强度参数的均方根误差
Table 6. Root mean square error of different ground motion intensity parameters
强度参数 Y=0 g Y=0.2 g Y=0.4 g Y=0.6 g Y=0.8 g Y=1.0 g Sa(T1) 0.391 8 0.410 8 0.443 7 0.455 9 0.464 7 0.469 9 PGA 0.697 0 0.714 4 0.695 3 0.677 7 0.672 3 0.670 9 PGV 0.419 7 0.453 3 0.487 7 0.499 7 0.514 4 0.525 1 PGD 0.419 5 0.423 3 0.438 0 0.437 4 0.441 4 0.444 2 -
[1] 陈力波, 郑凯锋, 庄卫林等, 2012. 汶川地震桥梁易损性分析. 西南交通大学学报, 47(4): 558—566 doi: 10.3969/j.issn.0258-2724.2012.04.004Chen L. B. , Zheng K. F. , Zhuang W. L. , et al. , 2012. Analytical investigation of bridge seismic vulnerability in wenchuan earthquake. Journal of Southwest Jiaotong University, 47(4): 558—566. (in Chinese) doi: 10.3969/j.issn.0258-2724.2012.04.004 [2] 谷音, 李晓芳, 2019. 考虑氯离子侵蚀的近海桥梁结构地震易损性分析. 应用基础与工程科学学报, 27(5): 1019—1032Gu Y. , Li X. F. , 2019. Seismic fragility analysis of offshore bridge structure considering chloride erosion. Journal of Basic Science and Engineering, 27(5): 1019—1032. (in Chinese) [3] Hwang H, 刘晶波, 2004. 地震作用下钢筋混凝土桥梁结构易损性分析. 土木工程学报, 37(6): 47—51 doi: 10.3321/j.issn:1000-131X.2004.06.009Hwang H, 刘晶波, 2004. 地震作用下钢筋混凝土桥梁结构易损性分析. 土木工程学报, 37(6): 47—51. Hwang H. , Liu J. B. , 2004. Seismic fragility analysis of reinfo钢筋混凝土ed concrete bridges. China Civil Engineering Journal, 37(6): 47—51. (in Chinese) doi: 10.3321/j.issn:1000-131X.2004.06.009 [4] 姬淑艳, 李英民, 刘立平等, 2006. 考虑双向水平地震作用的结构设计问题研究. 地震工程与工程振动, 26(5): 68—72 doi: 10.3969/j.issn.1000-1301.2006.05.011Ji S. Y. , Li Y. M. , Liu L. P. , et al. , 2006. Study on bi-directional seismic action in structural seismic design. Earthquake Engineering and Engineering Vibration, 26(5): 68—72. (in Chinese) doi: 10.3969/j.issn.1000-1301.2006.05.011 [5] 李立峰, 吴文朋, 胡思聪等, 2016. 考虑氯离子侵蚀的高墩桥梁时变地震易损性分析. 工程力学, 33(1): 163—170 doi: 10.6052/j.issn.1000-4750.2014.06.0530Li L. F. , Wu W. P. , Hu S. C. , et al. , 2016. Time-dependent seismic fragility analysis of high pier bridge based on chloride ion induced corrosion. Engineering Mechanics, 33(1): 163—170. (in Chinese) doi: 10.6052/j.issn.1000-4750.2014.06.0530 [6] 李勇, 闫维明, 刘晶波等, 2018. 近断层高架连续梁桥地震易损性与振动台试验研究. 工程力学, 35(4): 52—60, 86 doi: 10.6052/j.issn.1000-4750.2017.03.0236Li Y. , Yan W. M. , Liu J. B. , et al. , 2018. Seismic vulnerability analysis and a shaking table test of a near-fault continuous viaduct. Engineering Mechanics, 35(4): 52—60, 86. (in Chinese) doi: 10.6052/j.issn.1000-4750.2017.03.0236 [7] 林庆利, 林均岐, 刘金龙, 2017. 汶川地震公路桥梁易损性研究. 振动与冲击, 36(4): 110—118, 126Lin Q. L. , Lin J. Q. , Liu J. L. , 2017. A study on the fragility of highway bridges in the Wenchuan earthquake. Journal of Vibration and Shock, 36(4): 110—118, 126. (in Chinese) [8] 刘洋, 吕大刚, 于晓辉, 2016. 近场地震作用下型钢-混凝土组合结构桥易损性分析. 土木工程学报, 49(S1): 56—60, 77Liu Y. , Lv D. G. , Yu X. H. , 2016. Seismic fragility analysis of steel-concrete composite bridges excited by near-fault ground motions. China Civil Engineering Journal, 49(S1): 56—60, 77. (in Chinese) [9] 宋帅, 王帅, 吴刚等, 2020. 中小跨径梁桥地震易损性研究. 振动与冲击, 39(9): 118—125Song S. , Wang S. , Wu G. , et al. , 2020. Seismic vulnerability analysis of small and medium span girder bridges. Journal of Vibration and Shock, 39(9): 118—125. (in Chinese) [10] 徐超, 温增平, 2017. 不同设防标准钢筋混凝土框架结构基于易损性分析的抗震性能评估. 震灾防御技术, 12(4): 845—857 doi: 10.11899/zzfy20170413Xu C. , Wen Z. P. , 2017. Seismic performance evaluation of 钢筋混凝土 frames with different seismic precautionary intensity based on vulnerability analysis. Technology for Earthquake Disaster Prevention, 12(4): 845—857. (in Chinese) doi: 10.11899/zzfy20170413 [11] 易方民, 高小旺, 张维嶽等, 2003. 高层建筑钢结构在多维地震动输入作用下的反应. 建筑结构学报, 24(3): 33—43 doi: 10.3321/j.issn:1000-6869.2003.03.005Yi F. M. , Gao X. W. , Zhang W. Y. , et al. , 2003. Response of tall building steel structure to multi dimensional seismic action. Journal of Building Structures, 24(3): 33—43. (in Chinese) doi: 10.3321/j.issn:1000-6869.2003.03.005 [12] 张菊辉, 2006. 基于数值模拟的规则梁桥墩柱的地震易损性分析. 上海: 同济大学.Zhang J. H., 2006. Study on seismic vulnerability analysis of normal beam bridge piers based on numerical simulation. Shanghai: Tongji University. (in Chinese) [13] 中华人民共和国交通运输部, 2015. JTG D60-2015 公路桥涵设计通用规范. 北京: 人民交通出版社.Ministry of Transport of the People's Republic of China, 2015. JTG D60-2015 General specifications for design of highway bridges and culverts. Beijing: China Communications Press. (in Chinese) [14] 中华人民共和国交通运输部, 2019. JT/T 391—2019 公路桥梁盆式支座. 北京: 人民交通出版社.Ministry of Transport of the People's Republic of China, 2019. JT/T 391—2019 Pot bearing for highway bridge. Beijing: China Communications Press. (in Chinese) [15] 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局, 2010. GB 50011—2010 建筑抗震设计规范. 北京: 中国建筑工业出版社.Ministry of Housing and Urban-Rural Development of the People's Republic of China, General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, 2010. GB 50011—2010 Code for seismic design of buildings. Beijing: China Architecture & Building Press. (in Chinese) [16] 周长东, 田苗旺, 张许等, 2017. 考虑多维地震作用的高耸钢筋混凝土烟囱结构易损性分析. 土木工程学报, 50(3): 54—64Zhou C. D. , Tian M. W. , Zhang X. , et al. , 2017. Seismic fragility analysis for high-rise 钢筋混凝土 chimney considering multi-dimensional seismic actions. China Civil Engineering Journal, 50(3): 54—64. (in Chinese) [17] Cornell C. A. , Jalayer F. , Hamburger R. O. , 2002. Probabilistic basis for 2000 SAC federal emergency management agency steel moment frame guidelines. Journal of Structural Engineering, 128(4): 526—533. doi: 10.1061/(ASCE)0733-9445(2002)128:4(526) [18] Federal Emergency Management Agency (FEMA), 2001. HAZUS 99 Estimated annualized earthquake losses for the United States. Washington: FEMA. [19] Shinozuka M. , Feng M. Q. , Lee J. , et al. , 2000. Statistical analysis of fragility curves. Journal of Engineering Mechanics, 126(12): 1224—1231. doi: 10.1061/(ASCE)0733-9399(2000)126:12(1224) [20] Shome N., 1999. Probabilistic seismic demand analysis of nonlinear structures. Stanford: Stanford University. [21] Yamazaki F., Motomura H., Hamada T., 2000. Damage assessment of expressway networks in Japan based on seismic monitoring. In: Proceeding of the 12th World Conference on Earthquake Engineering, Upper Hutt, New Zealand: WCEE, 0551. 期刊类型引用(3)
1. 李中铭. 氯离子侵蚀作用下大跨度高墩刚构桥地震易损性分析. 福建交通科技. 2024(02): 31-35 . 百度学术
2. 梁晓,姜浩然,李芳芳. 自复位预制节段拼装中空夹层钢管混凝土桥墩地震易损性分析. 震灾防御技术. 2024(03): 613-628 . 本站查看
3. 许成祥,吴永康,王理. 外包钢套加固震损双层高架桥框架式桥墩抗震性能评估与参数分析. 震灾防御技术. 2023(04): 821-832 . 本站查看
其他类型引用(2)
-