考虑索端刚臂的斜拉桥空间拉索非线性分析
2014-03-17邓继华邵旭东张阳彭建新
邓继华 邵旭东 张阳 彭建新
摘要：为解决斜拉桥空间拉索在塔和梁上锚固点的刚性连接计算问题，在基于微分法导出考虑垂度效应的非线性空间杆元切线刚度矩阵的基础上，根据索端刚臂在受力后只有刚体运动而本身不变形的特点，将刚臂视为空间矢量，利用空间矢量有限转动公式及微分方法，导出了结构坐标系下刚臂两端的位移之间和杆端力之间的总量及增量关系.最终获得两端带任意刚臂的斜拉索空间杆元切线刚度矩阵表达式，编制了相应的斜拉桥几何非线性有限元计算程序， 对经典桁架算例、带刚臂的三杆空间桁架和施工阶段的斜拉桥进行了空间几何非线性分析. 计算结果表明：本文方法能解决斜拉桥空间拉索在塔和梁上锚固点的刚性连接计算问题，为斜拉桥进行几何非线性精细分析提供有力工具.
关键词：空间斜拉索；刚臂；几何非线性；切线刚度矩阵；空间矢量有限转动
中图分类号：TU323.3文献标识码：A
对斜拉桥进行几何非线性分析时，斜拉索一般采用杆单元进行模拟，通过等效模量法[1]考虑拉索垂度效应.在实际斜拉桥结构中，由于拉索两端必须锚固在梁和塔上，此时拉索两端的锚固点可能并不处于梁和塔的轴线上，并且锚固点附近一般为刚度非常大的刚性区域.在有限元建模时为解决此问题，可采用多编节点和单元的方式，且需将单元的刚度取得很大，此方法易导致计算误差增大、结点和单元数增加且可能改变索长甚至非线性计算不收敛等新的问题.文献[2]和[3]对索端锚固点的刚性连接问题进行了研究，但文献[2]仅研究两端带切线刚臂的特殊情况.文献[3]在推导过程中省略了单元已有内力对切线刚度矩阵附加项的影响.文献[4]指出在非线性程度较高的时候该作法会导致计算不收敛的问题，且文献[2]和[3]的研究均局限于平面范围，对空间范围内索端锚固点的刚性连接计算问题，目前尚未见到公开发表的研究文献.鉴于此，本文将斜拉索用空间杆元来模拟，在参考文献[5]利用微分法导出能考虑垂度效应的斜拉索非线性空间杆元切线刚度矩阵的基础上，根据索端刚臂在受力后只有刚体运动而本身不变形的特点， 将刚臂视为空间矢量，利用空间矢量有限转动公式，对刚臂两端力和位移微分，在分别得到杆端力之间和位移之间全量和增量关系基础上，最终得到两端带任意刚臂的斜拉索空间杆元切线刚度矩阵显式表达式，并通过数个算例对本文方法进行了验证.
3结论
本文根据索端刚臂在受力后的运动和变形特点，基于空间矢量有限转动公式，采用微分法导出带刚臂的空间杆元的切线刚度矩阵表达式.算例表明，该方法能较好解决斜拉桥空间拉索等索结构在塔和梁上锚固点的刚性连接计算问题，且具有列式清晰、力学概念明确的优点，为斜拉桥及悬索桥进行几何非线性精细分析提供有力工具.
参考文献
[1]ERNST H J. Der Emodul von seilen unter Berueck sichtigung des durchhanges[J]. Der Bauingenieur， 1965， 40（ 2） ： 52-55.
[2]罗喜恒， 肖汝城， 项海帆. 基于精确解析解的索单元[J]. 同济大学学报：自然科学版， 2005， 33（4）： 445-450.
[3]陈常松，颜东煌，陈政清. 带刚臂的两节点精确悬链线索元的非线性分析[J]. 工程力学， 2007， 24（5）： 29-34.
[4]邓继华， 邵旭东. 基于U.L列式的带刚臂平面梁元非线性分析[J]. 湖南大学学报：自然科学版， 2012， 39（5）： 8-12.
[5]张阳，邵旭东，蔡松柏，等. 大跨桁式钢管混凝土拱桥空间非线性有限元分析[J]. 中国公路学报， 2006， 19（4）： 65-70.
[6]朱菊芬，汪海，徐胜利. 非线性有限元及其在飞机结构设计中的应用[M]. 上海：上海交通大学出版社，2011.
[7]ARGYRIS J. An excursion into large rotation[J]. Computer Methods in Applied Mechanics and Engineering， 1982， 32：85-155.
[8]CRISFIELD M A. Accelerated solution techniques and concrete cracking[J]. Computer Methods in Applied Mechanics and Engineering， 1982， 33： 585-607.
[9]RUO S R， YANG Y B. Tracing postbucking paths of structures containing multiloops[J].International Journal for Numerical Methods in Engineering，1995，38：4035-4075.
[10]陈务军， 关富铃， 袁行飞， 等. 斜拉桥施工控制分析中线性与非线性影响分析[J]. 中国公路学报， 1999，11（2）： 52-58.
[11]颜东煌， 田仲初， 李学文. 桥梁结构电算程序设计[M].长沙： 湖南大学出版社， 1999.
摘要：为解决斜拉桥空间拉索在塔和梁上锚固点的刚性连接计算问题，在基于微分法导出考虑垂度效应的非线性空间杆元切线刚度矩阵的基础上，根据索端刚臂在受力后只有刚体运动而本身不变形的特点，将刚臂视为空间矢量，利用空间矢量有限转动公式及微分方法，导出了结构坐标系下刚臂两端的位移之间和杆端力之间的总量及增量关系.最终获得两端带任意刚臂的斜拉索空间杆元切线刚度矩阵表达式，编制了相应的斜拉桥几何非线性有限元计算程序， 对经典桁架算例、带刚臂的三杆空间桁架和施工阶段的斜拉桥进行了空间几何非线性分析. 计算结果表明：本文方法能解决斜拉桥空间拉索在塔和梁上锚固点的刚性连接计算问题，为斜拉桥进行几何非线性精细分析提供有力工具.
关键词：空间斜拉索；刚臂；几何非线性；切线刚度矩阵；空间矢量有限转动
中图分类号：TU323.3文献标识码：A
对斜拉桥进行几何非线性分析时，斜拉索一般采用杆单元进行模拟，通过等效模量法[1]考虑拉索垂度效应.在实际斜拉桥结构中，由于拉索两端必须锚固在梁和塔上，此时拉索两端的锚固点可能并不处于梁和塔的轴线上，并且锚固点附近一般为刚度非常大的刚性区域.在有限元建模时为解决此问题，可采用多编节点和单元的方式，且需将单元的刚度取得很大，此方法易导致计算误差增大、结点和单元数增加且可能改变索长甚至非线性计算不收敛等新的问题.文献[2]和[3]对索端锚固点的刚性连接问题进行了研究，但文献[2]仅研究两端带切线刚臂的特殊情况.文献[3]在推导过程中省略了单元已有内力对切线刚度矩阵附加项的影响.文献[4]指出在非线性程度较高的时候该作法会导致计算不收敛的问题，且文献[2]和[3]的研究均局限于平面范围，对空间范围内索端锚固点的刚性连接计算问题，目前尚未见到公开发表的研究文献.鉴于此，本文将斜拉索用空间杆元来模拟，在参考文献[5]利用微分法导出能考虑垂度效应的斜拉索非线性空间杆元切线刚度矩阵的基础上，根据索端刚臂在受力后只有刚体运动而本身不变形的特点， 将刚臂视为空间矢量，利用空间矢量有限转动公式，对刚臂两端力和位移微分，在分别得到杆端力之间和位移之间全量和增量关系基础上，最终得到两端带任意刚臂的斜拉索空间杆元切线刚度矩阵显式表达式，并通过数个算例对本文方法进行了验证.
3结论
本文根据索端刚臂在受力后的运动和变形特点，基于空间矢量有限转动公式，采用微分法导出带刚臂的空间杆元的切线刚度矩阵表达式.算例表明，该方法能较好解决斜拉桥空间拉索等索结构在塔和梁上锚固点的刚性连接计算问题，且具有列式清晰、力学概念明确的优点，为斜拉桥及悬索桥进行几何非线性精细分析提供有力工具.
参考文献
[1]ERNST H J. Der Emodul von seilen unter Berueck sichtigung des durchhanges[J]. Der Bauingenieur， 1965， 40（ 2） ： 52-55.
[2]罗喜恒， 肖汝城， 项海帆. 基于精确解析解的索单元[J]. 同济大学学报：自然科学版， 2005， 33（4）： 445-450.
[3]陈常松，颜东煌，陈政清. 带刚臂的两节点精确悬链线索元的非线性分析[J]. 工程力学， 2007， 24（5）： 29-34.
[4]邓继华， 邵旭东. 基于U.L列式的带刚臂平面梁元非线性分析[J]. 湖南大学学报：自然科学版， 2012， 39（5）： 8-12.
[5]张阳，邵旭东，蔡松柏，等. 大跨桁式钢管混凝土拱桥空间非线性有限元分析[J]. 中国公路学报， 2006， 19（4）： 65-70.
[6]朱菊芬，汪海，徐胜利. 非线性有限元及其在飞机结构设计中的应用[M]. 上海：上海交通大学出版社，2011.
[7]ARGYRIS J. An excursion into large rotation[J]. Computer Methods in Applied Mechanics and Engineering， 1982， 32：85-155.
[8]CRISFIELD M A. Accelerated solution techniques and concrete cracking[J]. Computer Methods in Applied Mechanics and Engineering， 1982， 33： 585-607.
[9]RUO S R， YANG Y B. Tracing postbucking paths of structures containing multiloops[J].International Journal for Numerical Methods in Engineering，1995，38：4035-4075.
[10]陈务军， 关富铃， 袁行飞， 等. 斜拉桥施工控制分析中线性与非线性影响分析[J]. 中国公路学报， 1999，11（2）： 52-58.
[11]颜东煌， 田仲初， 李学文. 桥梁结构电算程序设计[M].长沙： 湖南大学出版社， 1999.
摘要：为解决斜拉桥空间拉索在塔和梁上锚固点的刚性连接计算问题，在基于微分法导出考虑垂度效应的非线性空间杆元切线刚度矩阵的基础上，根据索端刚臂在受力后只有刚体运动而本身不变形的特点，将刚臂视为空间矢量，利用空间矢量有限转动公式及微分方法，导出了结构坐标系下刚臂两端的位移之间和杆端力之间的总量及增量关系.最终获得两端带任意刚臂的斜拉索空间杆元切线刚度矩阵表达式，编制了相应的斜拉桥几何非线性有限元计算程序， 对经典桁架算例、带刚臂的三杆空间桁架和施工阶段的斜拉桥进行了空间几何非线性分析. 计算结果表明：本文方法能解决斜拉桥空间拉索在塔和梁上锚固点的刚性连接计算问题，为斜拉桥进行几何非线性精细分析提供有力工具.
关键词：空间斜拉索；刚臂；几何非线性；切线刚度矩阵；空间矢量有限转动
中图分类号：TU323.3文献标识码：A
对斜拉桥进行几何非线性分析时，斜拉索一般采用杆单元进行模拟，通过等效模量法[1]考虑拉索垂度效应.在实际斜拉桥结构中，由于拉索两端必须锚固在梁和塔上，此时拉索两端的锚固点可能并不处于梁和塔的轴线上，并且锚固点附近一般为刚度非常大的刚性区域.在有限元建模时为解决此问题，可采用多编节点和单元的方式，且需将单元的刚度取得很大，此方法易导致计算误差增大、结点和单元数增加且可能改变索长甚至非线性计算不收敛等新的问题.文献[2]和[3]对索端锚固点的刚性连接问题进行了研究，但文献[2]仅研究两端带切线刚臂的特殊情况.文献[3]在推导过程中省略了单元已有内力对切线刚度矩阵附加项的影响.文献[4]指出在非线性程度较高的时候该作法会导致计算不收敛的问题，且文献[2]和[3]的研究均局限于平面范围，对空间范围内索端锚固点的刚性连接计算问题，目前尚未见到公开发表的研究文献.鉴于此，本文将斜拉索用空间杆元来模拟，在参考文献[5]利用微分法导出能考虑垂度效应的斜拉索非线性空间杆元切线刚度矩阵的基础上，根据索端刚臂在受力后只有刚体运动而本身不变形的特点， 将刚臂视为空间矢量，利用空间矢量有限转动公式，对刚臂两端力和位移微分，在分别得到杆端力之间和位移之间全量和增量关系基础上，最终得到两端带任意刚臂的斜拉索空间杆元切线刚度矩阵显式表达式，并通过数个算例对本文方法进行了验证.
3结论
本文根据索端刚臂在受力后的运动和变形特点，基于空间矢量有限转动公式，采用微分法导出带刚臂的空间杆元的切线刚度矩阵表达式.算例表明，该方法能较好解决斜拉桥空间拉索等索结构在塔和梁上锚固点的刚性连接计算问题，且具有列式清晰、力学概念明确的优点，为斜拉桥及悬索桥进行几何非线性精细分析提供有力工具.
参考文献
[1]ERNST H J. Der Emodul von seilen unter Berueck sichtigung des durchhanges[J]. Der Bauingenieur， 1965， 40（ 2） ： 52-55.
[2]罗喜恒， 肖汝城， 项海帆. 基于精确解析解的索单元[J]. 同济大学学报：自然科学版， 2005， 33（4）： 445-450.
[3]陈常松，颜东煌，陈政清. 带刚臂的两节点精确悬链线索元的非线性分析[J]. 工程力学， 2007， 24（5）： 29-34.
[4]邓继华， 邵旭东. 基于U.L列式的带刚臂平面梁元非线性分析[J]. 湖南大学学报：自然科学版， 2012， 39（5）： 8-12.
[5]张阳，邵旭东，蔡松柏，等. 大跨桁式钢管混凝土拱桥空间非线性有限元分析[J]. 中国公路学报， 2006， 19（4）： 65-70.
[6]朱菊芬，汪海，徐胜利. 非线性有限元及其在飞机结构设计中的应用[M]. 上海：上海交通大学出版社，2011.
[7]ARGYRIS J. An excursion into large rotation[J]. Computer Methods in Applied Mechanics and Engineering， 1982， 32：85-155.
[8]CRISFIELD M A. Accelerated solution techniques and concrete cracking[J]. Computer Methods in Applied Mechanics and Engineering， 1982， 33： 585-607.
[9]RUO S R， YANG Y B. Tracing postbucking paths of structures containing multiloops[J].International Journal for Numerical Methods in Engineering，1995，38：4035-4075.
[10]陈务军， 关富铃， 袁行飞， 等. 斜拉桥施工控制分析中线性与非线性影响分析[J]. 中国公路学报， 1999，11（2）： 52-58.
[11]颜东煌， 田仲初， 李学文. 桥梁结构电算程序设计[M].长沙： 湖南大学出版社， 1999.
