[发明专利]考虑弯曲刚度的斜拉索索力数值计算方法

说明书

技术领域

本发明涉及斜拉桥设计技术领域，特别是一种考虑弯曲刚度的斜拉索索力数值计算方法。

背景技术

当前普遍应用的斜拉索索力计算方法为通过对斜拉索运动方程进行近似和拟合，获得能考虑弯曲刚度和垂度效应的索力计算公式，称之为解析公式法求解索力。解析公式法的优点在于计算简单，不需要通过反复迭代即可获得斜拉索索力值。但是由于解析公式在近似和拟合中存在着假设和近似，并且所应用的斜拉索运动方程未能同时考虑弯曲刚度和垂度效应，使得其计算结果并不完全精确。随着计算技术特别是数值计算方法的不断发展，如果需要得到更为精确的索力值，另一种方法就是采用数值方法直接求解同时考虑弯曲刚度和垂度效应的斜拉索运动方程。

发明内容

本发明的目的在于提供一种考虑弯曲刚度的斜拉索索力数值计算方法，该方法有利于提高斜拉索索力测量结果的准确性。

为实现上述目的，本发明采用的技术方案是：一种考虑弯曲刚度的斜拉索索力数值计算方法，包括以下步骤：

步骤1：对斜拉索的振动进行测量，并对测量结果进行频谱分析，得到实测的振动频率；对斜拉索索力进行预估，得到预估索力值；

步骤2：将斜拉索自身特性值和预估索力值代入考虑弯曲刚度的斜拉索面内运动方程，并利用有限差分的数值方法进行求解，得到斜拉索振动的计算频率；

步骤3：计算斜拉索振动的计算频率与实测的振动频率之间的差值，若所述差值在误差允许范围内，则判定预估索力值即是斜拉索的实际索力值，算法结束；若所述差值超出误差允许范围，则调整预估索力值，并返回步骤2进行循环迭代，直到斜拉索振动的计算频率与实测的振动频率的差值落在误差允许范围内，调整后的预估索力值即为斜拉索的实际索力值。

本发明的有益效果是在斜拉索索力测量过程中考虑了斜拉索自身抗弯刚度对索力测量结果的影响，提高了斜拉索索力计算结果的准确性和精度。

附图说明

图1是本发明实施例的斜拉索的几何形状示意图。

图2是本发明实施例的斜拉索的离散模型示意图。

图3是本发明实施例的工作流程图。

具体实施方式

本发明考虑弯曲刚度的斜拉索索力数值计算方法，如图3所示，包括以下步骤：

步骤1：对斜拉索的振动进行测量，并对测量得到的加速度时程记录进行频谱分析，得到斜拉索振动的实测一阶频率和实测高阶频率；根据斜拉索自身特性值，结合现有技术的弦振动张力计算公式，计算得到预估索力值。

步骤2：将斜拉索自身特性值和预估索力值代入考虑弯曲刚度的斜拉索面内运动方程，并利用有限差分的数值方法进行求解，得到斜拉索振动的计算一阶频率和计算高阶频率。

步骤3：计算斜拉索振动的计算频率与实测频率之间的差值，若所述差值在误差允许范围内，则判定预估索力值即是斜拉索的实际索力值，算法结束；若所述差值超出误差允许范围，则对预估索力值进行适当调整，并返回步骤2进行循环迭代计算，直到斜拉索振动的计算频率与实测频率的差值落在误差允许范围内，调整后的预估索力值就是根据当前频率计算得到的斜拉索的实际索力值。

斜拉索的几何形状如图1所示。定义坐标系（x,z）为整体坐标系，坐标系（x*,z*）为整体割线坐标系，即x方向为水平方向，z方向为竖直方向，x*方向为斜拉索两端点连线方向，z*方向与x*方向垂直且与斜拉索在同一平面，L为斜拉索两端点间的距离，                                               为斜拉索水平倾角，在整体割线坐标系下，计算斜拉索振动的计算频率时，只考虑方程线性项且不考虑外力，考虑弯曲刚度和垂度的斜拉索无量纲运动方程为：

    (1)

    (2)

对方程（1）和（2）进行有限差分离散处理，可得：

    (3)

    (4)