《考虑弯曲刚度的斜拉索索力数值计算方法.pdf》由会员分享，可在线阅读，更多相关《考虑弯曲刚度的斜拉索索力数值计算方法.pdf（8页完整版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102735386 A (43)申请公布日 2012.10.17 C N 1 0 2 7 3 5 3 8 6 A *CN102735386A* (21)申请号 201210242823.6 (22)申请日 2012.07.14 G01L 5/04(2006.01) G01H 17/00(2006.01) (71)申请人福州大学 地址 350108 福建省福州市闽侯县上街镇大 学城学园路2号福州大学新区 (72)发明人吴庆雄 黄宛昆 陈宝春 (74)专利代理机构福州元创专利商标代理有限 公司 35100 代理人蔡学俊 (54) 发明名称 考虑弯曲刚度的斜拉索索力数值计算。
2、方法 (57) 摘要 本发明涉及一种考虑弯曲刚度的斜拉索索力 数值计算方法，包括以下步骤：1、对斜拉索的振 动进行测量，得到实测的振动频率；对斜拉索索 力进行预估，得到预估索力值；2、将斜拉索自身 特性值和预估索力值代入考虑弯曲刚度的斜拉索 面内运动方程，并利用有限差分的数值方法进行 求解，得到斜拉索振动的计算频率；3、计算斜拉 索振动的计算频率与实测频率之间的差值，若差 值在误差允许范围内，则判定预估索力值即是斜 拉索的实际索力值，算法结束；若差值超出误差 允许范围，则调整预估索力值，进行循环迭代，直 到计算频率与实测频率的差值落在误差允许范围 内，调整后的预估索力值即为实际索力值。该方法 。
3、有利于提高斜拉索索力测量结果的准确性。 (51)Int.Cl. 权利要求书2页 说明书3页 附图2页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 2 页 说明书 3 页 附图 2 页 1/2页 2 1.一种考虑弯曲刚度的斜拉索索力数值计算方法，其特征在于：包括以下步骤： 步骤1：对斜拉索的振动进行测量，并对测量结果进行频谱分析，得到实测的振动频 率；对斜拉索索力进行预估，得到预估索力值； 步骤2：将斜拉索弹性模量、斜拉索截面面积、斜拉索外径、斜拉索线密度、斜拉索水平 倾角、斜拉索索长、斜拉索等分节段数和预估索力值代入考虑弯曲刚度的斜拉索面内运动 方程，并利用有限差分。
4、的数值方法进行求解，得到斜拉索振动的计算频率； 步骤3：计算斜拉索振动的计算频率与实测的振动频率之间的差值，若所述差值在误 差允许范围内，则判定预估索力值即是斜拉索的实际索力值，算法结束；若所述差值超出误 差允许范围，则调整预估索力值，并返回步骤2进行循环迭代，直到斜拉索振动的计算频率 与实测的振动频率的差值落在误差允许范围内，调整后的预估索力值即为斜拉索的实际索 力值。 2.根据权利要求1所述的考虑弯曲刚度的斜拉索索力数值计算方法，其特征在于：定 义坐标系（x,z）为整体坐标系，坐标系（x*,z*）为整体割线坐标系，计算斜拉索振动的计算 频率时，只考虑方程线性项且不考虑外力，考虑弯曲刚度和垂。
5、度的斜拉索无量纲运动方程 为： (1) (2) 对方程（1）和（2）进行有限差分离散处理，可得： (3) (4) 式(1)、(2)、(3)、(4)中，s是沿斜拉索长度方向坐 标，L为斜拉索两端点间的距离，为斜拉索水平倾角，是无量纲轴向额外动张力， ，是轴向额外动张力，H是斜拉索张力的水平分量，是无量纲轴向张力， ，T是斜拉索轴向张力，u*、w*是在 x*和z*方向的动位移，是斜拉索轴向刚度和水平拉力的比值，是弯曲刚 权 利 要 求 书CN 102735386 A 2/2页 3 度与轴向刚度的比值，E是斜拉索弹性模量，A是斜拉索横截面面积，I是斜拉 索抗弯刚度，是无量纲时间，是无垂度斜弦的一阶圆。
6、频率， 是斜拉索初始张力沿x*方向的分量，m是斜拉索线密度；是斜拉索额外动张力沿x*方 向的分量，i是斜拉索离散后的节段编号，t是时间； 从而得到斜拉索运动方程离散后的矩阵形式为：KW =P，对矩阵K进行特征值和特征向 量的求解，即可计算得到斜拉索振动的计算频率和振型。 权 利 要 求 书CN 102735386 A 1/3页 4 考虑弯曲刚度的斜拉索索力数值计算方法 技术领域 0001 本发明涉及斜拉桥设计技术领域，特别是一种考虑弯曲刚度的斜拉索索力数值计 算方法。 背景技术 0002 当前普遍应用的斜拉索索力计算方法为通过对斜拉索运动方程进行近似和拟合， 获得能考虑弯曲刚度和垂度效应的索力。
7、计算公式，称之为解析公式法求解索力。解析公式 法的优点在于计算简单，不需要通过反复迭代即可获得斜拉索索力值。但是由于解析公式 在近似和拟合中存在着假设和近似，并且所应用的斜拉索运动方程未能同时考虑弯曲刚度 和垂度效应，使得其计算结果并不完全精确。随着计算技术特别是数值计算方法的不断发 展，如果需要得到更为精确的索力值，另一种方法就是采用数值方法直接求解同时考虑弯 曲刚度和垂度效应的斜拉索运动方程。 发明内容 0003 本发明的目的在于提供一种考虑弯曲刚度的斜拉索索力数值计算方法，该方法有 利于提高斜拉索索力测量结果的准确性。 0004 为实现上述目的，本发明采用的技术方案是：一种考虑弯曲刚度的。
8、斜拉索索力数 值计算方法，包括以下步骤： 步骤1：对斜拉索的振动进行测量，并对测量结果进行频谱分析，得到实测的振动频 率；对斜拉索索力进行预估，得到预估索力值； 步骤2：将斜拉索自身特性值和预估索力值代入考虑弯曲刚度的斜拉索面内运动方 程，并利用有限差分的数值方法进行求解，得到斜拉索振动的计算频率； 步骤3：计算斜拉索振动的计算频率与实测的振动频率之间的差值，若所述差值在误 差允许范围内，则判定预估索力值即是斜拉索的实际索力值，算法结束；若所述差值超出误 差允许范围，则调整预估索力值，并返回步骤2进行循环迭代，直到斜拉索振动的计算频率 与实测的振动频率的差值落在误差允许范围内，调整后的预估索力。
9、值即为斜拉索的实际索 力值。 0005 本发明的有益效果是在斜拉索索力测量过程中考虑了斜拉索自身抗弯刚度对索 力测量结果的影响，提高了斜拉索索力计算结果的准确性和精度。 附图说明 0006 图1是本发明实施例的斜拉索的几何形状示意图。 0007 图2是本发明实施例的斜拉索的离散模型示意图。 0008 图3是本发明实施例的工作流程图。 具体实施方式 说 明 书CN 102735386 A 2/3页 5 0009 本发明考虑弯曲刚度的斜拉索索力数值计算方法，如图3所示，包括以下步骤： 步骤1：对斜拉索的振动进行测量，并对测量得到的加速度时程记录进行频谱分析，得 到斜拉索振动的实测一阶频率和实测高阶。
10、频率；根据斜拉索自身特性值，结合现有技术的 弦振动张力计算公式，计算得到预估索力值。 0010 步骤2：将斜拉索自身特性值和预估索力值代入考虑弯曲刚度的斜拉索面内运动 方程，并利用有限差分的数值方法进行求解，得到斜拉索振动的计算一阶频率和计算高阶 频率。 0011 步骤3：计算斜拉索振动的计算频率与实测频率之间的差值，若所述差值在误差 允许范围内，则判定预估索力值即是斜拉索的实际索力值，算法结束；若所述差值超出误差 允许范围，则对预估索力值进行适当调整，并返回步骤2进行循环迭代计算，直到斜拉索振 动的计算频率与实测频率的差值落在误差允许范围内，调整后的预估索力值就是根据当前 频率计算得到的斜拉。
11、索的实际索力值。 0012 斜拉索的几何形状如图1所示。定义坐标系（x,z）为整体坐标系，坐标系（x*,z*） 为整体割线坐标系，即x方向为水平方向，z方向为竖直方向，x*方向为斜拉索两端点 连线方向，z*方向与x*方向垂直且与斜拉索在同一平面，L为斜拉索两端点间的距离， 为斜拉索水平倾角，在整体割线坐标系下，计算斜拉索振动的计算频率时，只考虑方程线 性项且不考虑外力，考虑弯曲刚度和垂度的斜拉索无量纲运动方程为： (1) (2) 对方程（1）和（2）进行有限差分离散处理，可得： (3) (4) 式(1)、(2)、(3)、(4)中，x*是沿x*方向坐标，z*是沿z*方 向坐标，s是沿斜拉索长度方。
12、向坐标，L为斜拉索两端点间的距离，为斜拉索水平倾角， 是无量纲轴向额外动张力，是轴向额外动张力，H是斜拉索张力的水平分量， 是无量纲轴向张力，T是斜拉索轴向张力， 说 明 书CN 102735386 A 3/3页 6 ，u*、w*是在x*和z*方向的动位移，是斜拉索轴向刚度和水平拉力的比值， ，是弯曲刚度与轴向刚度的比值，E是斜拉索弹性模量，A是斜拉 索横截面面积，I是斜拉索抗弯刚度，是无量纲时间，是无垂度斜弦的一阶圆频 率，是斜拉索初始张力沿x*方向的分量，m是斜拉索线密度；是斜拉 索额外动张力沿x*方向的分量，i是斜拉索离散后的节段编号（如图2所示），t是时间； 从而得到斜拉索运动方程离散后的矩阵形式为：KW =P，对矩阵K进行特征值和特征向 量的求解，即可计算得到斜拉索振动的计算频率和振型。 0013 以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作 用未超出本发明技术方案的范围时，均属于本发明的保护范围。 说 明 书CN 102735386 A 1/2页 7 图1 图2 说 明 书 附 图CN 102735386 A 2/2页 8 图3 说 明 书 附 图CN 102735386 A 。