支座激励下小垂度悬索动力响应试验及数值模拟

王国威 张鹏 卢文胜

王国威, 张鹏, 卢文胜. 支座激励下小垂度悬索动力响应试验及数值模拟 [J]. 工程科学与技术, 2022, 54(6): 166-175. doi: 10.15961/j.jsuese.202101004
引用本文: 王国威, 张鹏, 卢文胜. 支座激励下小垂度悬索动力响应试验及数值模拟 [J]. 工程科学与技术, 2022, 54(6): 166-175. doi: 10.15961/j.jsuese.202101004
WANG Guowei, ZHANG Peng, LU Wensheng. Experiment and Numerical Simulation on the Dynamic Responses for Small-sagged Cables Under Support Excitations [J]. Advanced Engineering Sciences, 2022, 54(6): 166-175. doi: 10.15961/j.jsuese.202101004
Citation: WANG Guowei, ZHANG Peng, LU Wensheng. Experiment and Numerical Simulation on the Dynamic Responses for Small-sagged Cables Under Support Excitations [J]. Advanced Engineering Sciences, 2022, 54(6): 166-175. doi: 10.15961/j.jsuese.202101004

支座激励下小垂度悬索动力响应试验及数值模拟

基金项目: 国家自然科学基金项目(52178501;51638012);山东省高校协同创新重点课题(XTZ201903);江西理工大学校级纵向科研项目(205200100550)
详细信息
    • 收稿日期:  2021-10-08
    • 网络出版时间:  2022-11-09 11:25:00
  • 作者简介:

    王国威(1995—),男,博士生. 研究方向:玻璃幕墙柔性支撑体系. E-mail:guowei@tongji.edu.cn

    通信作者:

    卢文胜, 教授,E-mail: wally@tongji.edu.cn

  • 中图分类号: TU323.4

Experiment and Numerical Simulation on the Dynamic Responses for Small-sagged Cables Under Support Excitations

  • 摘要: 索网幕墙采光顶通常支承于两相邻结构之间,但主体结构的相对位移振动将作为支座激励导致悬索构件发生过度张拉或松弛现象。为研究悬索在端部支座简谐及随机激励下的动力响应,首先,以综合考虑垂跨比、索力及轴向刚度对悬索振动影响的无量纲变形特征参数 $ \alpha {b^2} $ 为设计变量,对3组悬索试件开展支座简谐激励下的动力响应试验,通过改变激励幅值和激励频率设置两组加载工况,分析在不同的 $ \alpha {b^2} $ 下,悬索位移和索力响应的变化规律;其次,采用ANSYS软件进行数值模拟,计算结果与试验结果吻合较好,验证了数值模拟方法的正确性;最后,以理想高斯白噪声模拟悬索支座随机环境激励,基于已验证的有限元模型计算悬索动力响应,并讨论 $ \alpha {b^2} $ 、阻尼比及激励强度对位移、索力响应谱及均方差的影响。结果表明:1)在支座简谐及随机激励下,悬索位移和索力响应均随 $ \alpha {b^2} $ 、阻尼比的增加而降低,随激励强度的增加而增大;且对于 $ \alpha {b^2} $ 较大的悬索,其位移和索力响应对支座的随机激励更敏感。2)当阻尼比小于0.1时,随着阻尼比的增加,悬索在随机激励下的位移和索力响应分别降低60.9%和67.5%,此阶段下悬索的动力响应对阻尼比变化的敏感性较高,且增加阻尼带来的抑制效果显著。3)当支座简谐激励和索力相同时,垂跨比的变化对悬索的动位移影响较低,但对索力的影响较大,其中垂跨比较小的悬索索力变化更剧烈。

     

    Abstract: The cable-net supported glass roof is commonly located between adjacent structures, however, the relative displacement vibration of the main structure will be considered as the support excitation leading to over-tensioning or relaxation of the small-sagged cables. In order to investigate dynamic responses of small-sagged cables under support harmonic and stochastic excitations. Dynamic responses of three sets of specimens under the support harmonic excitation are tested with the dimensionless deformation parameter $ \alpha {b^2} $ , as the design variable, which integrates the influence of the sag-span ratio, tensions, and linear stiffness on the vibration of small-sagged cables. The two loading cases are set up by adjusting the loading amplitude and frequency for analyzing displacements and cable tensions under different $ \alpha {b^2} $ . Then, by using ANSYS, the calculated results agree well with the experimental results, thus verifying the correctness of the numerical simulation method. Finally, taking the ideal Gaussian white noise as the stochastic excitation, the paper calculates the dynamic response based on the validated numerical model and discusses the effects of $ \alpha {b^2} $ , the damping ratio and the excitation intensity on response spectrums and mean square deviations. Results show that: 1) Displacements and tension responses under support harmonic and stochastic excitation decrease with increasing $ \alpha {b^2} $ and damping ratio, and increase with increasing excitation intensity. Especially for cables with large $ \alpha {b^2} $ , their displacement and tension responses are more sensitive to excitations; 2) When the damping ratio is less than 0.1, with the increased damping ratio, the displacement and tension responses of the cable under stochastic excitation decrease by 60.9% and 67.5%, respectively, indicating that the dynamic responses under this stage are more sensitive to the damping ratio, and the effect brought by increasing damping is significant; 3) When harmonic excitations and cable force are the same, the sag-span ratio has less effect on the dynamic displacement of small-sagged cables but has a significant effect on cable tensions, in which the change of tensions are more drastic for the cable with a smaller sag-span ratio.

     

  • 索网幕墙采光顶是近年来兴起的一种建筑外围护体系,因其美观轻盈、通透性能优越,被广泛应用于多功能商业中心、体育场馆和机场航站楼等大型公共设施,如北京金成大厦[1]、联想融科资讯中心[2]和德国海尔布隆车站[3]等。其中,悬索作为此类结构的基本承力单元,具有刚度低、非线性强的特点,在风、地震及环境振动等诸多不利因素作用下,有可能发生大变形和索力突变,尤其是连接于两主体结构的采光顶悬索,在强震作用下,主体结构若发生相对位移振动,悬索会发生过度张拉或松弛现象,甚至会出现断丝破坏。因此,悬索的静动力特性和动力响应研究受到广泛关注。

    国内外学者对悬索的静动力特性开展了深入的研究。Irvine等[4]根据线性理论,建立小垂跨比(小于1∶8)悬索的自由振动方程,得到综合反映悬索几何形状及材性的Irvine系数。Rega等[5]进一步考虑悬索的强非线性效应,建立了弹性悬索的大振幅自由振动偏微分方程,研究其垂跨比在较大范围内的非线性现象,并指出悬索的软、硬化行为取决于自身特性和振幅。沈世钊等[6]系统地阐述了悬索的静力解析和有限元计算方法,并通过部分算例验证方法的有效性。为考虑不同边界条件的影响,刘海涛等[7]采用多尺度法研究索–梁耦合结构的非线性行为,分析得出主共振下两者刚度比较大时,索从刚度硬化转变为软化,而梁的非线性行为消失;Zhang等[89]对比分析端部固定和增设弹簧的悬索结构的静动力特性,计算表明,边界约束降低会导致悬索变形增加,索力降低。

    在索网幕墙采光顶等实际工程中,支承悬索的主结构质量远大于索结构本身,故主结构的运动也被视为支座激励施加于悬索端部。Perkins[10]通过对悬索平面内、外耦合模态的近似描述,检验了所提出的两自由度悬索非线性振动方程的正确性。Benedettini等[11]基于连续化理论,建立支座激励下悬索的四自由度模型,并讨论了几何性质等参数对系统动力响应的影响。Guo等[12]考虑几何非线性和索–支承耦合效应建立悬索振动简化模型,并通过数值方法进行验证。此外,在索结构振动测试方面,胡俊等[13]实测了某悬索桥的风振响应,监测数据表明,日常风环境下吊索为高频涡激振动。贺文宇等[14]以两端固定的斜索为例,通过有限差分法计算索的模态参数,探究了斜拉索倾角对索力识别精度的影响。桂成中等[15]提出一种索力增量校验系数的方法,用于评定拉索和吊杆的技术状况。悬索在随机荷载下的响应研究则较少,Georgakis等[16]基于有限元法分析了斜拉索在竖向激励下的面内、外振动位移。任淑琰等[17]采用解析方法对斜拉索在轴向随机激励下的参数振动进行分析。王波等[18]研究了斜拉索在横桥向激励下面内的耦合振动特性,结果表明,在一定的横向激励下,拉索由于耦合作用其面内振动将被激发。

    悬索作为采光顶的主要受力构件因支撑玻璃面板具有一定的垂度,主体结构的相对位移振动将作为支座激励导致悬索构件产生不利的变形。因此,本文主要研究支座激励下具有一定垂度悬索的动力特性和动力响应。综上,现有研究大多集中在理论分析层面,从试验角度对悬索在支座激励下的动力响应研究不足,且多数文献尚未考虑支座随机激励及垂度对悬索动力特性和响应的影响。为进一步开展支座简谐激励下悬索动力特性、响应的试验和数值分析,对比了两者结果验证数值模型的正确性,并对支座随机激励下悬索的动力响应进行了参数研究,对相邻结构间索网幕墙采光顶的设计与应用具有一定的工程参考价值。

    为获取悬索在支座激励下的动力特性和动力响应,在同济大学土木工程防灾国家重点实验室设计并搭建试验平台,包括试验模型、加载装置及力与位移传感器等。试件采用直径为10 mm的1×19镀锌钢绞线,跨度为7.3 m,最小破断力为69.5 kN,其中单根钢丝强度为1320 MPa。将悬索试件一端通过耳板铰接于反力墙上,另一端同样与加载装置铰接连接,且在索上附加若干数量的配重块,由MTS作动器直接通过加载装置对悬索施加水平支座激励。

    首先,建立小垂度悬索面内的无阻尼自由振动频率方程[19],其归一化后的表达式为:

    $$ \tan \left( {{{{\omega _{n}}} / 2}} \right) = {{{\omega _{n}}} / 2} - \frac{{\omega _{n}^3}}{{128\alpha {b^2}}},{n = 2,4,6,\cdots} $$ (1)

    式中: ${\omega _{n}}$ 为悬索的正对称模态频率;n为频率阶数; $ \alpha \text{=}E\cdot A/H $ 为悬索轴向刚度与索力比, $ H $ 为初始索力, $ E $ $ A $ 分别为弹性模量和截面面积; $b = {w / l}$ 为悬索跨中垂度( $w$ )与跨度(l)之比。

    由表达式(1)可知, $ \alpha {b^2} $ 可以综合考虑垂跨比、索力及材料常数对悬索振动特性的影响。因此,下文将其作为关键参数研究悬索在支座激励下的动力特性和动力响应。

    为研究变形特征参数 $ \alpha {b^2} $ 的影响,设计3组悬索试件参数( $ \alpha {b^2} $ = 0.04、0.08、0.19),保持跨度不变,通过调整索力及配重数改变垂跨比。悬索Ⅰ、Ⅱ等间隔布置8个配重,悬索Ⅲ等间隔布置17个配重,单个配重为24 kg,具体悬索试件参数见表1

    表  1  悬索试件参数
    Table  1  Parameters of the small-sagged cable
    悬索
    编号
    $ \alpha {b^2} $ 索力
    H/kN
    跨度
    l/m
    静载下跨中
    垂度w/mm
    截面面积
    A/mm2
    弹性模量
    E/MPa
    0.04 24 7.3 80 58.5 1.25×105
    0.08 20 7.3 109 58.5 1.25×105
    0.19 24 7.3 180 58.5 1.25×105

    支座激励方式为简谐位移荷载,荷载工况分变幅值和变频率加载,考虑MTS作动器的加载能力,两组加载工况列于表2。试验的整个加载过程为:首先,通过悬索端部张紧装置施加1 kN的初始拉力,固定限位构件以防止安装配重过程中悬索发生松弛;其次,释放限位约束,检查各传感器及采集仪工作状态,并通过作动器控制索力在预定值附近;最后,采用正弦波位移激励: ${X_{\rm{d}}}\left( t \right) = D\sin \left( {2{\text{π}}ft} \right)$ ,使MTS作动器连续往复加载30周,两组工况的总加载时间t分别为25和60 s,加载完毕后让悬索自由振动,直至振动停止,测量悬索的频率和阻尼比。

    表  2  试验加载工况
    Table  2  Cases of the testing
    工况编号 激励幅值D/mm 激励频率f/Hz
    1 5 1.2
    2 12 0.5

    试验采用力传感器和位移计2种传感器。在悬索两端分别安装一个量程为5 t的力传感器,实时量测索力大小;在悬索跨度的1/3、1/2和2/3位置处各布置一个量程为750 mm的位移计,用来测量悬索面内的竖向动位移;力传感器和位移计的采样频率均设定为50 Hz,并通过数据采集仪记录各传感器数据。在悬索下方搭设防护架以防止配重块坠落损坏试验设备,试验测点及装置示意图和加载装置现场照片如图12所示。

    图  1  试验测点及装置示意图
    Fig.  1  Schematic diagram of test measuring points and device
    下载: 全尺寸图片
    图  2  悬索试验装置
    Fig.  2  Testing setup of the small-sagged cable specimen
    下载: 全尺寸图片
    1.2.1   悬索动力特性

    图3为各悬索在自由振动阶段由位移计实测数据计算所得的基频,从图3中可以看出:随着 $ \alpha {b^2} $ 的增加,悬索基频 $ {\omega _1} $ 从2.07 Hz降低至1.90 Hz。在边界条件相同的情况下,悬索的刚度随垂跨比的增加而降低。在支座简谐激励完毕后,以悬索Ⅰ ( $ \alpha {b^2} $ = 0.04)和悬索Ⅲ( $ \alpha {b^2} $ = 0.19)为例,提取跨中自由振动曲线并计算悬索等效黏滞阻尼比 $ \xi $ [20],自由衰减法计算公式如下:

    图  3  不同 ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$ 下悬索的基频 ${{\boldsymbol{\omega}} _{\bf{1}}}$
    Fig.  3  Base frequencies of the small-sagged cable with different ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$
    下载: 全尺寸图片
    $$ \ln \left( {\frac{{{v_r}}}{{{v_{r + s}}}}} \right) = \frac{{2s{\text{π}}\xi }}{{\sqrt {1 - {\xi ^2}} }} $$ (2)

    式中, $ {v_r} $ $ {v_{r + s}} $ 分别为rr + s时刻出现的相邻正波峰。

    图4为悬索的自由振动曲线,随着变形特征参数 $ \alpha {b^2} $ 的增加,自由振动曲线的指数分布规律特征逐渐明显,通过指数拟合得到悬索Ⅰ与Ⅲ的阻尼比分别为1.2%、0.8%。

    图  4  悬索Ⅰ、Ⅲ的自由振动曲线
    Fig.  4  Free vibration of the small-sagged Cable Ⅰ and Ⅲ
    下载: 全尺寸图片
    1.2.2   悬索动力响应

    以悬索跨度1/2位置处为例,选取悬索Ⅰ和Ⅲ为研究对象,分析 $ \alpha {b^2} $ 对支座激励下悬索动力响应及动力特性的影响,获得动力响应时程曲线,如图5所示。为便于对照,工况1和2均截取前25 s的时程曲线。由图5(a)可知:在工况1支座激励下,悬索面内动位移幅值 $ \left| {\Delta {w_{\max }}} \right| $ $ \alpha {b^2} $ 的增加而增大, $ \left| {\Delta {w_{\max }}} \right| $ 增加27.7 %。相反地,动索力幅值 $ \left| {\Delta {H_{\max }}} \right| $ $ \alpha {b^2} $ 的增加而明显减小, $ \left| {\Delta {H_{\max }}} \right| $ 降低69.1 %;当 $ \alpha {b^2} $ = 0.19时,除0 ~ 5 s时间段内的索力变化较大外,5 ~ 25 s时间段内, $ \left| {\Delta {H_{\max }}} \right| $ 均在1.5 kN附近波动,分析主要原因是:悬索的初始垂度较大,施加相同的支座激励,索力变化较小。由图5(b)可知,在工况2支座激励下,面内 $ \left| {\Delta {w_{\max }}} \right| $ $ \left| {\Delta {H_{\max }}} \right| $ 的变化也表现出相同的规律,动位移幅值略有增大,但 $ \left| {\Delta {H_{\max }}} \right| $ 降低29.4 %。

    图  5  不同 ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$ 下悬索动力响应时程曲线
    Fig.  5  Dynamic response curves of the small-sagged cable with different ${\boldsymbol{\alpha }}{{\boldsymbol{b}}^{\bf{2}}}$
    下载: 全尺寸图片

    采用ANSYS软件建立悬索数值模型,模型尺寸(l×A)、配重块数目及弹性模量E均与试验构件相同,泊松比μ=0.3,密度取钢材材料密度ρ=7 800 kg/m3。悬索截面较小,忽略自身抗弯刚度对面内振动的影响,采用Link 10仅受拉单元模拟,通过初位移 $ \Delta L $ 施加初始索力,配重块采用Mass 21结构质量单元模拟,不考虑转动惯量效应。悬索模型一端(A端)为铰接约束,另一端(B端)释放X方向约束,由于主要研究悬索面内(XY向)的振动,故同时约束索体Z向平动位移。采用绝对位移加载,在支座位置处施加与试验相同的正弦波位移激励 $ {X_{\text{d}}}\left( t \right) $ ,计算模型如图6所示。

    图  6  支座激励下悬索振动计算模型
    Fig.  6  Computing model of a small-sagged cable under support excitations
    下载: 全尺寸图片

    具体建模过程如下:

    1)根据试验构件建立几何模型,定义并赋予各部分材料属性,将模型划分为50个单元,并施加初始索力及设置模型边界条件。

    2)进行悬索静力计算,打开“NLGEOM”命令,考虑重力和几何非线性影响,对悬索进行找形分析,获得静载下悬索 $ \alpha {b^2} $ 的试验和有限元对比值。

    3)修正模型坐标,定义计算和扩展模态数,计算悬索的1阶频率,并与试验实测结果对比。

    4)采用完全法瞬态动力分析,假定各阶模态阻尼比相同,将实测阻尼比 $ \xi $ 转化为Rayleigh阻尼,并通过命令“ALPHAD”和“BETAD”定义阻尼系数;时间步长 $ \Delta t $ =0.05 s,采用水平支座简谐激励进行模拟试验加载;最后,提取悬索响应结果曲线并与试验结果对比分析。

    表3给出了3组悬索试件静载下索力H、挠度b $ \alpha {b^2} $ 及1阶频率的对比结果。由表3可知,有限元计算结果与试验结果非常接近,对于悬索Ⅰ、Ⅱ和Ⅲ,各参数误差均在5%~12%以内。

    表  3  静力及频率计算对比结果
    Table  3  Comparison results of static and frequency calculation
    悬索编号 H/kN b/mm $ \alpha {b^2} $ 1阶频率/Hz
    试验 有限元 试验 有限元 试验 有限元 试验 有限元
    24.0 24.0 80 81.4 0.04 0.04 2.07 2.12
    20.0 20.0 109 97.9 0.08 0.07 2.04 2.05
    24.0 23.9 180 162.4 0.19 0.16 1.90 1.84

    以悬索Ⅱ跨度1/2位置为例,图7给出了工况1和2下动力响应试验结果和数值模拟结果的对比,由图7可知:1)悬索Ⅱ在工况1下的动位移幅值 $ \left| {\Delta {w_{\max }}} \right| $ 和索力幅值 $ \left| {\Delta {H_{\max }}} \right| $ 变化较工况2更为明显, $ \left| {\Delta {w_{\max }}} \right| $ $ \left| {\Delta {H_{\max }}} \right| $ 分别为44.37 mm和4.02 kN,数值模拟和试验结果在部分峰值处略有差异,但平均索力和平均位移变化幅值误差均在12 %以内,主要原因是:在高频的支座激励下,悬索上的配重块刚心与质心不完全重合,配重发生少量偏转,对刚度有一定的影响,从而导致计算与试验存在误差,但两者结果的幅值变化趋势和相位均保持同步。2)相比于工况1,工况2的位移和索力变化相对平稳,且数值模拟和试验结果吻合较好,数值模拟和试验的 $ \left| {\Delta {w_{\max }}} \right| $ $ \left| {\Delta {H_{\max }}} \right| $ 的误差分别为5.4%、6.8%。综合上述结果可知,悬索动力特性及响应试验结果与数值模拟结果均吻合较好,验证了有限元模型的准确性。

    图  7  悬索Ⅱ( ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$ = 0.08)试验结果与数值模拟结果对比
    Fig.  7  Comparison of test results and numerical simulations results for Cable Ⅱ( ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$ = 0.08)
    下载: 全尺寸图片

    在前文验证有限元建模方法正确的前提下,以悬索Ⅰ( $ \alpha {b^2} $ =0.04)的有限元模型为例,令支座简谐激励频率 $f$ =2.07 Hz,研究不同阻尼比 $ \xi $ 下悬索的1阶共振响应,总加载时间为15 s,获得共振响应时程曲线如图8所示。从图8可以看出:当阻尼比 $ \xi $ = 0.01时,动位移和动索力呈周期性的“震荡”变化,但其幅值 $ \left| {\Delta {w_{\max }}} \right| $ $ \left| {\Delta {H_{\max }}} \right| $ 总体呈下降的趋势;当阻尼比 $ \xi $ = 0.20时,位移和索力变化均呈现平稳变化, $ \left| {\Delta {w_{\max }}} \right| $ $ \left| {\Delta {H_{\max }}} \right| $ 分别稳定在20 mm、4 kN附近。值得注意的是,随着阻尼比的增加,悬索1阶共振下位移和索力的响应均受到不同程度的控制,位移控制较索力下降明显, $ \left| {\Delta {w_{\max }}} \right| $ 降低64.0 %, $ \left| {\Delta {H_{\max }}} \right| $ 降低33.2 %。从侧面反映出,若在工程中的索网采光顶悬索端部或连接件部位布置黏滞、黏弹或摩擦阻尼器等装置,将会有效地降低悬索的位移和索力响应。

    图  8  不同阻尼比 $ {\boldsymbol{\xi}} $ 下悬索1阶共振响应时程曲线( $ {\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)
    Fig.  8  First-order resonance responses of the small-sagged cable with different $ {\boldsymbol{\xi}} $ ( ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)
    下载: 全尺寸图片

    风、地震及环境振动等属于非确定性作用,均具有一定的随机性,因此,研究悬索在随机激励下的动力响应具有重要的现实意义。为讨论悬索在支座随机激励下的动力响应,在前文工作的基础上,以理想高斯白噪声模拟随机环境激励,分析悬索控制参数 $ \alpha {b^2} $ 、阻尼比 $ \xi $ 和激励强度S0对悬索随机动力响应的影响。

    基于有限元法将悬索离散,对于阻尼较小的悬索结构,由振型j产生的响应与振型k的响应近似满足统计独立[6],则位移响应谱密度 $ {S _{{\text{u}}i}}\left( \omega \right) $ 可表示为:

    $$ \begin{aligned}[b]& {S _{{\text{u}}i}}\left( \omega \right) \approx \sum\limits_{j = 1}^J {{S _{{\text{u}}ij}}\left( \omega \right)} = \sum\limits_{j = 1}^J {\phi _{ij}^2} {\left| {{H_j}\left( {i\omega } \right)} \right|^2}{S _{{F_j}}}\left( \omega \right) , \\& {\left| {{H_j}\left( {i\omega } \right)} \right|^2}{\text{ = }}\frac{1}{{\omega _j^4\left\{ {{{\left[ {1 - {{\left( {{\omega / {{\omega _j}}}} \right)}^2}} \right]}^2}{\text{ + }}{{\left[ {2\xi {\omega / {{\omega _j}}}} \right]}^2}} \right\}}} \end{aligned} $$ (3)

    式中, $ {S _{{F_j}}}\left( \omega \right) $ 为激励的谱密度函数, $ \omega $ 为激励圆频率, $ {\omega _j} $ 为悬索第j振型圆频率, $ {\phi _{ij}} $ 为第i自由度的第j振型值,J为振型总数, $\left| {{H_j}\left( {i\omega } \right)} \right|^2 $ 为传递函数的模。对式(3)进行积分求和,可得位移响应均方差:

    $$ {\sigma _{{\text{u}}i}} = {\left[ {\sum\limits_{j = 1}^n {\int_{ - \infty }^\infty {S _{{\text{u}}ij}\left( \omega \right){\text{d}}\omega } } } \right]^{1/2}} $$ (4)

    随机激励下悬索第e个单元的索力响应 $ \Delta {T_e} $ 为:

    $$ \Delta {T_e} = {{EA( {l_e^1 - l_e^0} )} / {l_e^0}} $$ (5)

    式中, $ l_e^0 $ $ l_e^1 $ 分别是静载和随机激励下的索长,可通过e单元两端pq节点的位置坐标计算。考虑到悬索面内振动位移起控制作用,当忽略位移响应的2阶矩对索力响应 $ \Delta {T_e} $ 的影响时,其索力响应均方差可由pq节点面内位移 $ {w_p} $ $ {w_q} $ 及其均方差 $ {\sigma _{wp}} $ $ {\sigma _{wq}} $ 表示:

    $$ {\sigma _{\Delta {T_e}}} = {{EA\left( {{w_p} - {w_q}} \right)\left( {{\sigma _{wp}} - {\sigma _{wq}}} \right)} / {l_e^2}} $$ (6)

    采用均值为0、激励强度为S0的理想高斯白噪声 $ D\left( t \right) $ 作为支座水平激励,其自相关函数可表示为: $E\left[ {D\left( t \right)D\left( {t + \tau } \right)} \right] = 2{\text{π }} {S _0} \delta \left( \tau \right)$ ,其中, ${\delta}\left( \tau \right)$ 为Dirac–δ函数, $ \tau $ 为时间间隔, $ E[\cdot] $ $ [\cdot] $ 的期望。以试验中悬索Ⅰ的几何和材料参数为例,在ANSYS中建立支座随机激励下悬索的数值模型,并进行随机振动计算(PSD功率谱密度分析)。主要建模过程如下:

    1)悬索几何建模、单元类型、单元划分及边界约束与简谐激励相同,将初始索力H直接施加在释放X方向自由度的索端支座处。

    2)静力计算完毕后,定义计算模态和扩展模态,同时令“ELCALC = YES”计算应力结果。

    3)采用位移谱作为水平支座激励,由于软件输入的要求,将双边谱转换为单边谱进行加载,定义计算频率范围 $ \bar f $ 为0.1~10.0 Hz,激励强度S0为0.01~0.20 ${{{{\text{m}}^{\text{2}}}} \mathord{\left/ {\vphantom {{{{\text{m}}^{\text{2}}}} {{\text{Hz}}}}} \right. } {{{\rm{(rad·s^{-1})}}}}}$

    4)在PSD分析过程中,支座激励的缩放系数取为1,采用自动时间步,阻尼比 $ \xi $ 通过命令“DMPRAT”施加,采用“PSDCOM”命令计算悬索位移及索力响应的均方差(1σ),即位于–1σ~1σ区间发生概率为68.3 %的结果量。

    首先,研究参数 $ \alpha {b^2} $ 对支座随机激励下悬索动力响应的影响,计算不同 $ \alpha {b^2} $ 下的悬索位移和索力响应,结果如图9所示。

    图  9  不同 $ {\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ 下悬索的位移、索力响应(S0=0.01, ${\boldsymbol{ \xi}} $ =0.05)
    Fig.  9  Responses of the small-sagged cable with different ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ (S0=0.01, ${\boldsymbol{ \xi}} $ =0.05)
    下载: 全尺寸图片

    图9(a)(b)为3组 $ \alpha {b^2} $ 的悬索跨中1/2位置处的位移、索力响应谱,可以发现,两者结果符合正态分布。由于 $ \alpha {b^2} $ 与结构特性直接相关,随着 $ \alpha {b^2} $ 的增加,悬索动力响应谱幅值出现左移,说明结构频率逐渐降低,这也与图3中试验的实测结果相一致。在相同的激励强度S0和阻尼比 $ \xi $ 下,悬索动力响应随 $ \alpha {b^2} $ 的增加而显著加大,位移响应谱幅值 $\left| {\Delta {{\overline w}_{\max }}} \right|$ 增加了3倍,索力响应谱幅值 $\left| {\Delta {{\overline H}_{\max }}} \right|$ 甚至增加了20倍。图9(c)为悬索响应的均方差与 $ \alpha {b^2} $ 的关系,可知: $ \alpha {b^2} $ 从0.04增加至1.0的过程中,位移和索力的均方差呈上升趋势,但增加的幅度随 $ \alpha {b^2} $ 的增加而逐渐降低;当 $ \alpha {b^2} $ = 1.0时,悬索的位移均方差为0.339 m,增加了3.3倍,索力均方差为2.1 kN,增加了3.17倍。该规律反映出 $ \alpha {b^2} $ 较大的悬索,其位移和索力响应对支座的随机激励更加敏感。

    取不同的悬索阻尼比 $ \xi $ ,分析其对悬索位移及索力响应的影响。不同 $ \xi $ 下悬索的动力响应谱如图10(a)(b)所示,可知:当 $ \xi $ =0.01时,位移和索力响应谱幅值最大,分别为0.58 m2/Hz和11.55 kN2/Hz;当 $ \xi $ =0.05时,位移和索力响应谱幅值已降低至0.02 m2/Hz和0.47 kN2/Hz。对比两组结果,可以发现,悬索位移和索力响应已被显著抑制。从图10(c)悬索响应的均方差曲线也能看出,位移和索力均方差随着阻尼比 $ \xi $ 的增加而降低,当 $ \xi $ = 0.20时,两者分别减小68.1 %和75.3 %;另外,两者曲线下降趋势逐渐变缓,表明当 $ \xi $ < 0.10时,位移和索力响应对阻尼比的变化敏感性较高,增加 $ \xi $ 带来的抑制效果显著。上述规律与支座简谐激励下悬索1阶共振响应随阻尼比 $ \xi $ 的变化相同,也进一步说明增加悬索系统阻尼将明显降低位移和索力响应。

    图  10  不同 ${\boldsymbol{ \xi }}$ 下悬索的位移、索力响应(S0=0.01, ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)
    Fig.  10  Responses of the small-sagged cable with different ${\boldsymbol{ \xi}} $ (S0=0.01, ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)
    下载: 全尺寸图片

    若支座所受的随机环境激励强度S0发生改变,系统的动力响应也将随之变化。计算不同激励强度S0下的悬索动力响应,如图11所示。

    图  11  不同S0下悬索的位移、索力响应( $ {\boldsymbol{\xi }}$ =0.05, ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)
    Fig.  11  Responses of the small-sagged cable with different S0 ( $ {\boldsymbol{\xi }} $ =0.05, ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)
    下载: 全尺寸图片

    图11表明,在相同的 $ \alpha {b^2} $ 和阻尼比 $ \xi $ 下, $ \left| {\Delta {{\overline w}_{\max }}} \right| $ $ \left| {\Delta {{\overline H}_{\max }}} \right| $ 及位移、索力均方差均随S0的增加而增大;对应的 $ \left| {\Delta {{\overline w}_{\max }}} \right| $ 增加了7.1倍, $ \left| {\Delta {{\overline H}_{\max }}} \right| $ 增加了6.9倍,响应均方差曲线斜率逐渐减小,但总体呈上升趋势。

    本文首先进行了悬索在支座简谐激励下动力特性及响应的试验和数值模拟研究,对比两者结果验证模型的正确性;其次,通过数值方法计算悬索在支座随机荷载激励下的动力响应,并分析响应谱及均方差随 $ \alpha {b^2} $ 、阻尼比 $ \xi $ 和激励强度S0关键参数的变化规律,得到的主要结论如下:

    1)设计并开展悬索在支座简谐激励下的动力响应试验,分析不同 $ \alpha {b^2} $ 下悬索动力响应的变化:在两组工况下,悬索动位移随 $ \alpha {b^2} $ 的增加而小幅增加;相反地,动索力随 $ \alpha {b^2} $ 的增加而明显增大。结果表明,相同的支座简谐激励下,相比于位移响应,索力响应对悬索的关键参数 $ \alpha {b^2} $ 更加敏感。

    2)建立悬索试件有限元模型,模型结果与悬索自由振动试验实测频率、静载下实测垂跨比及各工况下悬索试件动力响应试验结果吻合较好,从而验证了有限元模型的正确性。

    3)对支座简谐激励下悬索的1阶共振响应进行有限元计算,结果表明,悬索动位移受阻尼比 $ \xi $ 的影响更大,且时程响应变化平稳。

    4)在支座随机激励下,通过数值方法讨论了 $ \alpha {b^2} $ $ \xi $ S0对悬索位移及索力响应的影响。计算表明,两种响应谱值和均方差均随参数 $ \alpha {b^2} $ 、激励强度S0的增加而增大,随阻尼比 $ \xi $ 的增加而减小,且变化幅度逐渐降低,说明3种参数对悬索动力响应的影响会逐步减弱。

  • 图  1   试验测点及装置示意图

    Fig.  1   Schematic diagram of test measuring points and device

    下载: 全尺寸图片

    图  2   悬索试验装置

    Fig.  2   Testing setup of the small-sagged cable specimen

    下载: 全尺寸图片

    图  3   不同 ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$ 下悬索的基频 ${{\boldsymbol{\omega}} _{\bf{1}}}$

    Fig.  3   Base frequencies of the small-sagged cable with different ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$

    下载: 全尺寸图片

    图  4   悬索Ⅰ、Ⅲ的自由振动曲线

    Fig.  4   Free vibration of the small-sagged Cable Ⅰ and Ⅲ

    下载: 全尺寸图片

    图  5   不同 ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$ 下悬索动力响应时程曲线

    Fig.  5   Dynamic response curves of the small-sagged cable with different ${\boldsymbol{\alpha }}{{\boldsymbol{b}}^{\bf{2}}}$

    下载: 全尺寸图片

    图  6   支座激励下悬索振动计算模型

    Fig.  6   Computing model of a small-sagged cable under support excitations

    下载: 全尺寸图片

    图  7   悬索Ⅱ( ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$ = 0.08)试验结果与数值模拟结果对比

    Fig.  7   Comparison of test results and numerical simulations results for Cable Ⅱ( ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}}$ = 0.08)

    下载: 全尺寸图片

    图  8   不同阻尼比 $ {\boldsymbol{\xi}} $ 下悬索1阶共振响应时程曲线( $ {\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)

    Fig.  8   First-order resonance responses of the small-sagged cable with different $ {\boldsymbol{\xi}} $ ( ${\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)

    下载: 全尺寸图片

    图  9   不同 $ {\boldsymbol{\alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ 下悬索的位移、索力响应(S0=0.01, ${\boldsymbol{ \xi}} $ =0.05)

    Fig.  9   Responses of the small-sagged cable with different ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ (S0=0.01, ${\boldsymbol{ \xi}} $ =0.05)

    下载: 全尺寸图片

    图  10   不同 ${\boldsymbol{ \xi }}$ 下悬索的位移、索力响应(S0=0.01, ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)

    Fig.  10   Responses of the small-sagged cable with different ${\boldsymbol{ \xi}} $ (S0=0.01, ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)

    下载: 全尺寸图片

    图  11   不同S0下悬索的位移、索力响应( $ {\boldsymbol{\xi }}$ =0.05, ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)

    Fig.  11   Responses of the small-sagged cable with different S0 ( $ {\boldsymbol{\xi }} $ =0.05, ${\boldsymbol{ \alpha}} {{\boldsymbol{b}}^{\bf{2}}} $ =0.04)

    下载: 全尺寸图片

    表  1   悬索试件参数

    Table  1   Parameters of the small-sagged cable

    悬索
    编号
    $ \alpha {b^2} $ 索力
    H/kN
    跨度
    l/m
    静载下跨中
    垂度w/mm
    截面面积
    A/mm2
    弹性模量
    E/MPa
    0.04 24 7.3 80 58.5 1.25×105
    0.08 20 7.3 109 58.5 1.25×105
    0.19 24 7.3 180 58.5 1.25×105

    表  2   试验加载工况

    Table  2   Cases of the testing

    工况编号 激励幅值D/mm 激励频率f/Hz
    1 5 1.2
    2 12 0.5

    表  3   静力及频率计算对比结果

    Table  3   Comparison results of static and frequency calculation

    悬索编号 H/kN b/mm $ \alpha {b^2} $ 1阶频率/Hz
    试验 有限元 试验 有限元 试验 有限元 试验 有限元
    24.0 24.0 80 81.4 0.04 0.04 2.07 2.12
    20.0 20.0 109 97.9 0.08 0.07 2.04 2.05
    24.0 23.9 180 162.4 0.19 0.16 1.90 1.84
  • [1] 张付奎,张自刚.金成大厦采光顶幕墙工程设计及施工[J].建筑结构,2009,39(7):94–96.

    Zhang Fukui,Zhang Zigang.Design and construction in the glass roof of Jincheng Mansion[J].Building Structure,2009,39(7):94–96
    [2] 陈彬磊,刘笛,钱基宏,等.融科C座悬索玻璃大堂结构设计[J].建筑结构,2006,36(6):69–72.

    Chen Binlei,Liu Di,Qian Jihong,et al.Structural design of suspended curtain hall in raycom C project[J].Building Structure,2006,36(6):69–72
    [3] Schober H,Schneider J.Developments in structural glass and glass structures[J].Structural Engineering International,2004,14(2):84–87. doi: 10.2749/101686604777964044
    [4] Irvine H M,Caughey T K.The linear theory of free vibrations of a suspended cable[J].Proceedings of the Royal Society A,1974,341(1626):299–315.
    [5] Rega G,Vestroni F,Benedettini F.Parametric analysis of large amplitude free vibrations of a suspended cable[J].International Journal of Solids and Structures,1984,20(2):95–105. doi: 10.1016/0020-7683(84)90001-5
    [6] 沈世钊,徐崇宝,赵臣.悬索结构设计[M].2版.北京:中国建筑工业出版社,2006.
    [7] 刘海涛,魏明海,肖仪清,等.索–梁耦合结构非线性分析[J].振动与冲击,2015,34(14):147–152. doi: 10.13465/j.cnki.jvs.2015.14.025

    Liu Haitao,Wei Minghai,Xiao Yiqing,et al.Nonlinear response analysis of a cable–beam coupled system[J].Journal of Vibration and Shock,2015,34(14):147–152 doi: 10.13465/j.cnki.jvs.2015.14.025
    [8] Zhang Peng,Zhu Huanian,Lu Wensheng,et al.Vibration analysis of shallow cable with horizontal spring and dashpot at one end[J].Engineering Structures,2020,211:110452. doi: 10.1016/j.engstruct.2020.110452
    [9] 张鹏,卢文胜,吕西林.附加弹簧保护装置的悬索主共振响应分析[J].振动工程学报,2020,33(3):550–558. doi: 10.16385/j.cnki.issn.1004-4523.2020.03.014

    Zhang Peng,Lu Wensheng,Lü Xilin.Primary resonance responses study of suspended cables with supplemental spring-protection devices[J].Journal of Vibration Engineering,2020,33(3):550–558 doi: 10.16385/j.cnki.issn.1004-4523.2020.03.014
    [10] Perkins N C.Modal interactions in the non-linear response of elastic cables under parametric/external excitation[J].International Journal of Non-Linear Mechanics,1992,27(2):233–250. doi: 10.1016/0020-7462(92)90083-J
    [11] Benedettini F,Rega G,Alaggio R.Non-linear oscillations of a four-degree-of-freedom model of a suspended cable under multiple internal resonance conditions[J].Journal of Sound and Vibration,1995,182(5):775–798. doi: 10.1006/jsvi.1995.0232
    [12] Guo Tieding,Kang Houjun,Wang Lianhua,et al.Cable dynamics under non-ideal support excitations:Nonlinear dynamic interactions and asymptotic modelling[J].Journal of Sound and Vibration,2016,384:253–272. doi: 10.1016/j.jsv.2016.08.020
    [13] 胡俊,欧进萍.某大跨悬索桥风振响应现场实测与理论对比分析[J].振动与冲击,2013,32(11):149–154. doi: 10.3969/j.issn.1000-3835.2013.11.030

    Hu Jun,Ou Jinping.Comparison between field measurement and theoretical analysis results for a long-span suspension bridge’s wind-induced vibration responses[J].Journal of Vibration and Shock,2013,32(11):149–154 doi: 10.3969/j.issn.1000-3835.2013.11.030
    [14] 贺文宇,孟凡成,任伟新.斜拉索倾角对振动法测索力的精度影响[J].振动、测试与诊断,2021,41(4):792–796.

    He Wenyu,Meng Fancheng,Ren Weixin.Influence of inclination angle on accuracy of vibration based cable force identification[J].Journal of Vibration,Measurement & Diagnosis,2021,41(4):792–796
    [15] 桂成中,宋国华,柯在田,等.拉索索力增量校验系数检定方法研究[J].工程科学与技术,2019,51(4):75–83.

    Gui Chengzhong,Song Guohua,Ke Zaitian,et al.Investigation on rating method of cable force increment adjustment factor[J].Advanced Engineering Sciences,2019,51(4):75–83
    [16] Georgakis C T,Taylor C A.Nonlinear dynamics of cable stays.Part 2:Stochastic cable support excitation[J].Journal of Sound and Vibration,2005,281(3/4/5):565–591.
    [17] 任淑琰,顾明.斜拉索在轴向白噪声随机激励下的响应[J].振动工程学报,2008,21(6):535–541. doi: 10.3969/j.issn.1004-4523.2008.06.001

    Ren Shuyan,Gu Ming.Parametric vibration of stay cables induced by axial white noise excitation[J].Journal of Vibration Engineering,2008,21(6):535–541 doi: 10.3969/j.issn.1004-4523.2008.06.001
    [18] 王波,张海龙,徐丰,等.随机横桥向激励下斜拉索面内耦合振动特性研究[J].振动与冲击,2008,27(8):59–63. doi: 10.3969/j.issn.1000-3835.2008.08.014

    Wang Bo,Zhang Hailong,Xu Feng,et al.In-plane coupled vibration of inclined cables under random transverse excitation[J].Journal of Vibration and Shock,2008,27(8):59–63 doi: 10.3969/j.issn.1000-3835.2008.08.014
    [19] Nayfeh A,Arafat H,Chin C M,et al.Multimode interactions in suspended cables[C]//Proceedings of the 19th AIAA Applied Aerodynamics Conference.Anaheim:AIAA,2001:1660.
    [20] Kawashima K,Unjoh S,Tunomoto M.Estimation of damping ratio of cable-stayed bridges for seismic design[J].Journal of Structural Engineering,1993,119(4):1015–1031. doi: 10.1061/(ASCE)0733-9445(1993)119:4(1015)
图(11)  /  表(3)

本文结构

    /

    返回文章
    返回