0 引言

大跨弯连续刚构桥，其结构属于墩梁永久固结体系，该桥型具有较强的抗弯能力和抗扭能力，在我国西部山岭地区应用广泛^[1]。由于跨越能力较大，桥梁在外荷载作用下，动力响应较为突出。尤其是在车载作用下，桥梁典型截面时程曲线和冲击系数等指标反应较为敏感。因此，开展大跨弯连续刚构桥车桥耦合分析显得尤为重要。近些年，各国学者对车桥耦合进行广泛研究，以实现桥梁的正常运行^[2-4]。Y. B. Yang and B. H. Lin等^[5]通过有限元建模分析车辆与简支梁桥的相互作用接触力，并通过有限元插值函数得到桥梁单元内部任一点处的荷载，相关研究成果为本研究理论分析提供思路。凌敏^[6]基于西北某高速连续刚构桥梁为依托工程，分析车速、桥面不平整度等因素对桥梁结构典型截面动挠度影响；沈锐利^[7]采用ANSYS分别建立平面3轴车型和桥梁模型，考虑桥面不平整度的影响，并通过数理统计的方法得到桥梁冲击系数的概率分布。邓露^[8]通过建立13座不同截面的简支梁桥，研究桥梁基频、桥面不平整度、车速等因素对冲击系数的影响。综上可知：车桥耦合针对常规桥型的理论分析较为丰富，但是针对大跨弯连续刚构桥的的耦合分析还相对较少。因此，针对此类特殊桥梁的特点，开展多敏感元参数下大跨弯连续刚构桥冲击系数研究尤为重要，本研究采用主成分-逐步回归法确定冲击系数的计算公式，相关研究成果可为此类桥梁动力分析和冲击系数的计算提供参考。

1 冲击系数计算方法 1.1 理论法

车辆荷载作用桥梁结构产生动态响应^[9]，为简化计算，通常采用冲击系数加以分析：

(1)

式中，Y_dmax为车辆过桥时，桥梁动态响应的最大值；Y_jmax为桥梁对应位置的静力效应最大。

1.2 规范法

中国公路桥涵通用规范^[10] (JTG D60—2015)：冲击系数与结构基频f的函数表达式为：

(2)

式中f为桥梁结构的基频。

2 车桥耦合振动基本理论 2.1 车-桥耦合基本假定^[11]

(1) 桥梁的结构刚度和质量特性应均匀分配桥跨长度范围内。

(2) 忽略梁横截面的变形。

(3) 桥梁的阻尼采用瑞利阻尼。

2.2 车辆振动方程建立

本研究采用3轴5自由度平面车辆模型，其计算图示如图 2所示。

图中，M_c为车体质量; z为车辆竖向位移; α为车辆横轴旋转的自由度; m₁，m₂，m₃为车辆悬架和车轮等效质量；z₁，z₂，z₃为车辆悬架和车轮等效质量的竖向位移; z′₁，z′₂，z′₃为支承车体点的竖向位移。

车辆的待求未知位移向量为：

(2)

车辆模型中惯性力、悬置力(F_si，i=1, 2，3)、轮胎力(F_ti，i=1，2，3)。由牛顿第二定律可知：

(3)

由几何关系可得：

(4)

由广义虚功原理得：

(5)

将公式(3)、(4)代入公式(5)，整理后，如式(6)所示：

(6)

式中，[M_v]为车体质量矩阵；[C_v]为车体阻尼矩阵；[K_v]为车体刚度矩阵；{F_bv}为车桥耦合荷载向量；{G_v}为重力荷载向量。

(7)

2.3 桥梁动力平衡方程建立

同理，可建立桥梁结构动力平衡方程：

(8)

式中，M_b为桥梁质量矩阵；C_b为桥梁阻尼矩阵；K_b为桥梁刚度矩阵；，，Y为桥梁节点处加速度、速度和位移向量；F为桥梁节点力列向量^[11]。

2.4 桥面不平整度

桥面不平整度，又称桥面粗糙度，表征桥面恶化程度的指标。在实际工程中，由于桥面不平整度的影响，使得车桥振动较为明显，过大的振动会降低乘车人舒适性。因此，在进行车桥耦合分析时，应考虑该因素的影响。本研究根据三角级数叠加法^[12]进行模拟，其公式如(9)所示：

(9)

式中，R(x)为桥面不平整度样本点；N为采样频段数；θ_k为服从[0, 2π]均匀分布的随机相位角；S(Ω_k)为功率谱密度函数；Ω_k为车辆空间频率；ΔΩ-Ω_k为间隔带宽^[14]；

通过调研，我国桥面不平整度大部分处于A~C级，因此，本研究模拟生成3种桥面不平整度曲线, 如图 3所示。

2.5 结构动力响应求解

基于分离迭代法原理，采用Matlab编译车桥耦合分析模块，实现车辆过桥振动响应分析。

2.5.1 位移联系方程建立

假设车轮在运动时始终与桥面保持接触，则车轮位移可通过桥梁的相对位移表示^[13]，如图 4所示：

(10)

式中，y_i为轮胎由平衡位置对应竖向位移；r_i为车轮作用点处的不平整度竖向坐标；w_i为轮胎作用下桥梁瞬时竖向变位；Δy_i为车辆轴悬挂弹簧与轮胎竖向变形的相对竖向位移。

2.5.2 车桥相互作用

将车辆和桥梁视作两个分离体系，二者之间耦合作用通过轮胎与桥面间的相互作用联系。其相互作用力为：

(11)

式中，k_ti为第i个轮胎的刚度系数；c_tj为第i个轮胎的阻尼系数。

3 主成分-逐步回归基本原理 3.1 主成分分析

在大跨连续刚构桥动力特性的敏感参数分析时，由于车桥耦合作用受多种因素的共同作用，在计算分析时，如果采用较多参数进行分析，会使得计算过程较为繁琐，且变量间的多重相关性会影响估计结果，造成模型的误差增大。

主成分分析将多指标参数进行降维处理^[14]，采用几个变量来替代多个变量的统计方法。

(1) 主成分标准化

在对建立模型进行主成分分析时，考虑多个变量间存在不同数量级和单位，在进行相互比较时，为消除量纲对计算模型的影响，需要提前对相关变量进行标准化变换。

(12)

式中，i=1, 2, …, n；j=1, 2, …, k；；

经标准化后的样本矩阵均为无量纲矩阵，且每个列向量的均值为0，标准差为1。

(2) 求解特征值λ_i和特征向量e_j

特征方程通常采用雅可比法(Jacobi)求解特征值λ_i，并根据特征值的大小进行排列，其中当λ₁≥λ₂≥…≥λ_p>0时，相应正则化特征向量为：

(13)

(3) 主成分方差贡献率的确定

在实际工程中，主要通过方差贡献率的累积量确定主成分的个数。其中λ_i的方差贡献率可通过式(14)求解。

(14)

式中，α_i为代表第i个因素的方差贡献率；α(m)为前m个因素的方差累计贡献率。

由式(14)可知：主成分的方差贡献率越大，表明该主成分保留因变量的信息能力越强，如果前m(m≤p)个主成分的累计值达到一定数值(本研究取方差累计贡献率T_i≥90%)，表明这m个主成分涵盖了样本的绝大部分信息，则这m个主成分可作为因变量分析的主要敏感参数。

3.2 逐步回归分析

当冲击系数受m个因素影响时，可通过多元线性方程来描述，其关系式为：

(15)

逐步回归^[15]是一种双向筛选法，其本质属于前进法，即变量从少到多引入时，采用最小二乘法建立多元回规方程，基于偏回归平方和F检验, 判断方程是否显著，剔除其中不显著的影响因素，其主要步骤如下：

(1) 对主成分分析中的相关自变量进行回归分析，分别计算因变量与自变量间样本点的拟合度R²。

(2) 对各个自变量所对应的R²进行排序。

(3) 选取最大值R²变量A_i建立回归模型，计算R₁，F₁。

(4) 选取排序第二R²的变量A_j的回归模型, 计算R₂，F₂。

若R₁ < R₂，F₁ < F₂, 且A_i仍然显著，则A_j的引入使得各个变量指标更好，则引入A_j，否则删掉A_j。

(5) 重复步骤(4)，进行下一个变量筛选，直到所有主成分变量都筛选后，结束计算。

4 敏感参数下大跨弯连续刚构桥冲击系数分析 4.1 工程简介

本研究以某高墩大跨连续弯刚构桥为工程背景，桥梁全长239 m，跨径组合(65+108+65) m，上部结构采用单箱单室截面，曲率半径R=960 m，下部结构采用薄壁空心墩，墩高60 m，主梁采用C50混凝土，桥墩采用C40混凝土^[11]，其相关尺寸和仿真模型如图 6~7所示。

4.2 敏感参数主成分分析

根据相关学者的研究^[16-17]，本研究选取了桥面不平整度、曲率半径、桥墩高度、行车速度、车辆重量作为大跨弯连续刚构桥车桥耦合分析的敏感参数。为简化计算，本研究对5个敏感参数进行主成份分析，确定其相对重要的敏感参数，见表 1。

4.2.1 敏感参数降维分析

在敏感参数降维分析时，采用德尔菲法对敏感参数进行评估，具体步骤如下:

(1) 确定调研主旨和目的，设计敏感参数调查问卷。

(2) 选取在车桥耦合方面有较为丰富科研或工程经验的专家。

(3) 发放调查问卷，并与专家建立函询工作，搜集相关结果。

(4) 对结果进行统计汇总，并将计算结果进行主成分降维分析。

(5) 搜集原始资料

按照敏感因素对车桥耦合冲击系数的影响程度划分10个等级，使用0~10进行评分，见表 2所示。

向上述专家进行问卷调查共计40人，回收筛选出有效样本36份，进一步整理数据，得到待评价矩阵{A_ij}_n×k，其中k为评价因素，即k=5, n为有效样本。

4.2.2 敏感参数主成分分析

采用SPSS22.0进行数据处理，得到相关系数矩阵R特性值，方差贡献率和累计方程贡献率，如表 3所示。

通过表 3可知，特征值大于1的主成分个数g=3，此时方差累计贡献率为95.852%，超过了主成分分析的90%，即涵盖了大部分信息，这表明前3个指标能代表5个指标来分析大跨弯连续刚构桥冲击系数变化情况，通过文献[18]可知，车速对于桥梁位移、弯矩和扭矩冲击系数的影响幅度不大，验证了降维方法的正确性，因此，本研究后续对前3个指标进行敏感性分析。

4.3 敏感参数下桥梁动力时程曲线

大跨弯连续刚构桥冲击系数受多个敏感参数共同作用。本研究通过拟定，曲率半径取值为250，500，750和960 m；墩高取值为60，80和100 m；路面不平整度等级取值为A，B和C级。求解桥梁关键截面的挠度、弯矩、扭矩的时程曲线，如图 8所示，其中车辆的基本参数如表 4所示。

4.4 大跨弯连续刚构桥冲击系数回归分析 4.4.1 关键截面冲击系数最大值确定

通过公式(1)计算各个关键截面处挠度、弯矩和扭矩冲击系数，从中找出挠度，弯矩和扭矩冲击系数最大值，如图 9所示。

4.4.2 桥梁冲击系数回归分析

通过采用SPSS22.0对图 9中相关数据进行线性逐步回归分析，得出各个敏感参数对冲击系数的贡献性大小，结果如表 5所示。

由表 5可知：

(1) 检验回归系数全部为0，小于sig显著性分析界限值0.05，表明自变量参数与因变量参数的线性相关性较好，模型1挠度冲击系数分析时，主要取边跨跨中截面和中跨跨中截面处挠度冲击系数的最大值作为因变量，其因变量受桥面不平整度、墩高的影响性较为显著，而曲率半径对其影响相对较小，可忽略该因素的影响；模型2弯矩冲击系数主要考虑跨中截面正弯矩和支点负弯矩的最大值作为因变量，其值受桥面不平整度的影响较为显著，墩高和曲率半径对其影响相对较小，可忽略不计；模型3扭矩冲击系数分析时，其因变量受桥面不平整度和曲率半径的影响较大，而墩高影响较小，可忽略不计。

(2) 挠度、弯矩和扭矩冲击系数的回归公式相关系数R均在0.94以上，表明该回归公式相关性较好；且各回归公式均与桥面不平整度正相关，其标准化相关系数Beta均在0.9以上，表明该因素对大跨弯连续刚构桥冲击系数的影响很大，而现有规范未考虑该因素的影响，建议后续规范修订中考虑该因素的影响。

5 结论

(1) 采用ANSYS软件建立桥梁模型，通过Matlab软件编译车桥耦合分析模块，基于位移协调方程构建车桥耦合分析模型。

(2) 采用主成分-逐步回归分析方法去确定桥面不平整度、曲率半径和墩高作为桥梁冲击系数分析的主要敏感参数，其方差累计贡献率达到95%，并通过SPSS22.0进行逐步回归，确定了桥梁挠度、弯矩和扭矩冲击系数的计算公式。

(3) 冲击系数回归公式相关系数R均在0.94以上，表明该回归公式相关性较好，且各个公式中桥面不平整度标准化相关系数Beta均在0.9以上，表明该因素对大跨弯连续刚构桥冲击系数的影响很大，而现有规范未考虑该因素的影响，建议后续规范修订中考虑该因素的影响。