发明名称 考虑弯曲刚度的斜拉索索力数值计算方法
摘要 本发明涉及一种考虑弯曲刚度的斜拉索索力数值计算方法,包括以下步骤:1、对斜拉索的振动进行测量,得到实测的振动频率;对斜拉索索力进行预估,得到预估索力值;2、将斜拉索自身特性值和预估索力值代入考虑弯曲刚度的斜拉索面内运动方程,并利用有限差分的数值方法进行求解,得到斜拉索振动的计算频率;3、计算斜拉索振动的计算频率与实测频率之间的差值,若差值在误差允许范围内,则判定预估索力值即是斜拉索的实际索力值,算法结束;若差值超出误差允许范围,则调整预估索力值,进行循环迭代,直到计算频率与实测频率的差值落在误差允许范围内,调整后的预估索力值即为实际索力值。该方法有利于提高斜拉索索力测量结果的准确性。
申请公布号 CN102735386B 申请公布日期 2014.05.14
申请号 CN201210242823.6 申请日期 2012.07.14
申请人 福州大学 发明人 吴庆雄;黄宛昆;陈宝春
分类号 G01L5/04(2006.01)I;G01H17/00(2006.01)I 主分类号 G01L5/04(2006.01)I
代理机构 福州元创专利商标代理有限公司 35100 代理人 蔡学俊
主权项 1.一种考虑弯曲刚度的斜拉索索力数值计算方法,其特征在于:包括以下步骤:步骤1:对斜拉索的振动进行测量,并对测量结果进行频谱分析,得到实测的振动频率;对斜拉索索力进行预估,得到预估索力值;步骤2:将斜拉索弹性模量、斜拉索截面面积、斜拉索外径、斜拉索线密度、斜拉索水平倾角、斜拉索索长、斜拉索等分节段数和预估索力值代入考虑弯曲刚度的斜拉索面内运动方程,并利用有限差分的数值方法进行求解,得到斜拉索振动的计算频率;步骤3:计算斜拉索振动的计算频率与实测的振动频率之间的差值,若所述差值在误差允许范围内,则判定预估索力值即是斜拉索的实际索力值,算法结束;若所述差值超出误差允许范围,则调整预估索力值,并返回步骤2进行循环迭代,直到斜拉索振动的计算频率与实测的振动频率的差值落在误差允许范围内,调整后的预估索力值即为斜拉索的实际索力值;定义坐标系(x,z)为整体坐标系,坐标系(x*,z*)为整体割线坐标系,计算斜拉索振动的计算频率时,只考虑方程线性项且不考虑外力,考虑弯曲刚度和垂度的斜拉索无量纲运动方程为:<img file="82651DEST_PATH_IMAGE001.GIF" wi="480" he="52" />(1)<img file="453108DEST_PATH_IMAGE002.GIF" wi="481" he="52" />(2)对方程(1)和(2)进行有限差分离散处理,可得:<img file="694734DEST_PATH_IMAGE003.GIF" wi="438" he="46" />(3)<img file="907540DEST_PATH_IMAGE004.GIF" wi="569" he="97" />(4)式(1)、(2)、(3)、(4)中,<img file="551011DEST_PATH_IMAGE005.GIF" wi="58" he="42" />,<img file="823861DEST_PATH_IMAGE006.GIF" wi="57" he="42" />,<img file="552782DEST_PATH_IMAGE007.GIF" wi="44" he="42" />,<i>s</i>是沿斜拉索长度方向坐标,<i>L</i>为斜拉索两端点间的距离,<img file="569280DEST_PATH_IMAGE008.GIF" wi="14" he="20" />为斜拉索水平倾角,<img file="2012102428236100001DEST_PATH_IMAGE009.GIF" wi="16" he="18" />是无量纲轴向额外动张力,<img file="67257DEST_PATH_IMAGE010.GIF" wi="81" he="42" />,<img file="2012102428236100001DEST_PATH_IMAGE011.GIF" wi="14" he="16" />是轴向额外动张力,<i>H</i>是斜拉索张力的水平分量,<img file="776587DEST_PATH_IMAGE012.GIF" wi="17" he="21" />是无量纲轴向张力,<img file="727226DEST_PATH_IMAGE013.GIF" wi="238" he="64" />,T是斜拉索轴向张力,<img file="547414DEST_PATH_IMAGE014.GIF" wi="60" he="42" />,<img file="899898DEST_PATH_IMAGE015.GIF" wi="64" he="42" />,u*、w*是在<i>x*</i>和<i>z*</i>方向的动位移,<img file="842446DEST_PATH_IMAGE016.GIF" wi="21" he="22" />是斜拉索轴向刚度和水平拉力的比值,<img file="218064DEST_PATH_IMAGE017.GIF" wi="90" he="42" />,<img file="638681DEST_PATH_IMAGE018.GIF" wi="13" he="20" />是弯曲刚度与轴向刚度的比值,<img file="111251DEST_PATH_IMAGE019.GIF" wi="78" he="46" />,E是斜拉索弹性模量,A是斜拉索横截面面积,I是斜拉索抗弯刚度,<img file="224700DEST_PATH_IMAGE020.GIF" wi="48" he="25" />,<img file="87614DEST_PATH_IMAGE022.GIF" wi="9" he="15" />是无量纲时间,<img file="2012102428236100001DEST_PATH_IMAGE023.GIF" wi="21" he="25" />是无垂度斜弦的一阶圆频率,<img file="311922DEST_PATH_IMAGE024.GIF" wi="92" he="48" />,<img file="308172DEST_PATH_IMAGE025.GIF" wi="26" he="18" />是斜拉索初始张力沿x*方向的分量,m是斜拉索线密度;<img file="592523DEST_PATH_IMAGE026.GIF" wi="24" he="18" />是斜拉索额外动张力沿x*方向的分量,i是斜拉索离散后的节段编号,t是时间;从而得到斜拉索运动方程离散后的矩阵形式为:KW =P,对矩阵K进行特征值和特征向量的求解,即可计算得到斜拉索振动的计算频率和振型。 
地址 350108 福建省福州市闽侯县上街镇大学城学园路2号福州大学新区
您可能感兴趣的专利