留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于协调翘曲场的开闭口混合薄壁截面杆件约束扭转一维有限元分析

文颖 陈泽林

文颖, 陈泽林. 基于协调翘曲场的开闭口混合薄壁截面杆件约束扭转一维有限元分析[J]. 工程力学, 2020, 37(9): 38-49. doi: 10.6052/j.issn.1000-4750.2019.11.0639
引用本文: 文颖, 陈泽林. 基于协调翘曲场的开闭口混合薄壁截面杆件约束扭转一维有限元分析[J]. 工程力学, 2020, 37(9): 38-49. doi: 10.6052/j.issn.1000-4750.2019.11.0639
Ying WEN, Ze-lin CHEN. ONE-DIMENSIONAL FINITE ELEMENT ANALYSIS OF WARPING TORSION FOR THIN-WALLED MEMBERS WITH OPEN-CLOSED CROSS SECTIONS BASED ON COMPATIBLE WARPING FIELD[J]. Engineering Mechanics, 2020, 37(9): 38-49. doi: 10.6052/j.issn.1000-4750.2019.11.0639
Citation: Ying WEN, Ze-lin CHEN. ONE-DIMENSIONAL FINITE ELEMENT ANALYSIS OF WARPING TORSION FOR THIN-WALLED MEMBERS WITH OPEN-CLOSED CROSS SECTIONS BASED ON COMPATIBLE WARPING FIELD[J]. Engineering Mechanics, 2020, 37(9): 38-49. doi: 10.6052/j.issn.1000-4750.2019.11.0639

基于协调翘曲场的开闭口混合薄壁截面杆件约束扭转一维有限元分析

doi: 10.6052/j.issn.1000-4750.2019.11.0639
基金项目: 国家自然科学基金项目(51478475,51108460);中国博士后科学基金项目(2012M511759)
详细信息
    作者简介:

    陈泽林(1995−),男,陕西人,硕士生,主要从事薄壁结构分析研究(E-mail: 13054170084@163.com)

    通讯作者: 文 颖(1981−),男,湖南人,副教授,博士,重载铁路工程结构教育部重点实验室副主任,主要从事桥梁稳定极限承载力及车桥系统振动稳定性研究(E-mail: ywen_ce@csu.edu.cn)
  • 中图分类号: TU323

ONE-DIMENSIONAL FINITE ELEMENT ANALYSIS OF WARPING TORSION FOR THIN-WALLED MEMBERS WITH OPEN-CLOSED CROSS SECTIONS BASED ON COMPATIBLE WARPING FIELD

  • 摘要: 开口及闭口薄壁杆件约束扭转问题已由经典Timoshenko和Benscoter理论解决。然而,开闭口混合薄壁截面杆件约束扭转分析必须考虑开、闭口部分翘曲能力的差异,翘曲剪流形成机理有待进一步研究。该文假定开、闭口截面翘曲分别满足Vlasov和Umanskii假定,考虑开、闭口截面公共节点翘曲连续性要求,建立含有待定翘曲参数的协调翘曲模型。由截面受力平衡,确定翘曲参数显式列式,提出开闭口混合薄壁截面杆件约束扭转分析的一维有限元模型。算例及参数分析结果表明,基于Umanskii第二理论的I类方法在悬臂板及闭口周边引入附加剪流,影响翘曲剪应力精度。基于Umanskii第二理论的II类方法只能计算截面板件平均剪应力,无法反映真实翘曲剪流分布。基于Vlasov约束扭转假定的Beam-189单元忽略闭口周边约束效应产生的附加翘曲及剪流,影响翘曲正应力和剪应力精度。该文方法与Shell-63单元能得到基本吻合的变形与应力结果,说明一维梁元能正确反映开闭口混合薄壁截面杆件约束扭转及翘曲刚度。
  • 图  1  开闭口混合截面布置

    Figure  1.  Cross sectional layout for open-closed profile

    图  2  开闭口混合截面主扇性坐标

    Figure  2.  Principal sectorial coordinates for open-closed profile

    图  3  开闭口混合薄壁截面杆件单元

    Figure  3.  Thin-walled element with open-closed cross section

    图  4  简支箱梁跨中截面翘曲正应力随悬臂板宽度变化曲线

    Figure  4.  Correlations between the warping normal stresses at mid-span and the width of the cantilever slab of a simply- supported box girder

    图  5  简支箱梁端截面(x=30)剪应力随悬臂板宽度变化曲线

    Figure  5.  Correlations between the shear stresses at the right end and the width of the cantilever slab of a simply-supported box girder

    图  6  简支箱梁跨中截面翘曲正应力随梁高变化曲线

    Figure  6.  Correlations between the warping normal stresses at mid-span and the height of a simply-supported box girder

    图  7  简支箱梁端截面(x=30)剪应力随梁高变化曲线

    Figure  7.  Correlations between the shear stresses at the right end and the height of a simply-supported box girder

    表  1  悬臂箱梁约束扭转位移和正应力

    Table  1.   Displacements and normal stresses for warping torsion of a cantilever box girder

    计算模型自由端扭转角/(×10−5 rad)自由端翘曲/(×10−5 m)/初等梁理论翘曲/(×10−5 m)固定端正应力/MPa/初等梁理论应力/MPa
    本文一维梁元−4.95−1.838/−1.84−10.72/−10.71−14.21/−14.24−0.091/−0.092−0.545/−0.536−0.687/−0.712
    Shell-63单元−5.05−1.7/−1.84−11.0/−10.71−14.0/−14.24−0.092/−0.092−0.54/−0.536−0.68/−0.712
    文献[7]−0.091/−0.092−0.598/−0.536−0.55/−0.712
    注:翘曲位移和应力的正负号与Ox坐标轴指向保持一致。
    下载: 导出CSV

    表  2  均布扭矩作用下简支箱梁约束扭转位移和应力

    Table  2.   Displacements and stresses for warping torsion of a simply-supported box girder under the action of uniformly distributed torques

    计算模型跨中扭转角/
    (×10−5 rad)
    跨中正应力/kPa梁端剪应力/kPa
    本文一维梁元2.2511.15−9.425.14−3.73−74.15
    Beam-189模型2.1510.86−9.245.01−16.69−51.79
    Shell-63模型2.2210.91−9.365.22−2.82−74.14
    Umanskii第二理论[8]2.2511.41−9.605.32−4.79−71.18
    注:与文献[8]不同,⑦点剪应力包含自由扭转剪应力的贡献。
    下载: 导出CSV

    表  3  跨中集中偏载作用下简支箱梁正应力和扭转角

    Table  3.   Normal stresses and angle of twist of a simply-supported box girder subjected to an eccentric transverse load applied at mid-span

    计算模型跨中扭转角/
    (×10−2rad)
    跨中正应力/(×108Pa)/
    弯曲正应力/(×108Pa)
    本文一维梁元2.802−2.69/−2.93−3.14/−2.933.69/3.57
    Beam-189模型2.678−2.70/−2.93−3.12/−2.933.67/3.57
    Shell-63模型2.783−2.82/−2.93−3.32/−2.934.01/3.57
    Umanskii第二理论2.812−2.79/−2.93−3.05/−2.933.63/3.57
    下载: 导出CSV
  • [1] Vlasov V Z. Thin-walled elastic beams, second edition [M]. Jerusalem: Israel Program for Scientific Translations, 1961: 3 − 51, 222 − 255.
    [2] Timoshenko S P, Gere J M. Theory of elastic stability [M]. 2nd ed. New York: McGraw-Hill Book Company, 1961: 212 − 224.
    [3] Murray N W. Introduction to the theory of thin-walled structures [M]. New York: Oxford University Press, 1984: 54 − 134.
    [4] 徐勋, 叶华文, 强士中. 带悬臂板薄壁箱梁的扭转和畸变分析[J]. 铁道学报, 2015, 37(10): 83 − 91. doi:  10.3969/j.issn.1001-8360.2015.10.012

    Xu Xun, Ye Huawen, Qiang Shizhong. Torsion and distorsion analysis of thin-walled box girder with cantilever flanges [J]. Journal of the China Railway Society, 2015, 37(10): 83 − 91. (in Chinese) doi:  10.3969/j.issn.1001-8360.2015.10.012
    [5] 郭金琼, 房贞政, 郑振. 箱形梁设计理论 [M]. 第2版. 北京: 人民交通出版社, 2008: 51 − 92.

    Guo Jinqiong, Fang Zhenzheng, Zheng Zhen. Design theory of box girder [M]. 2nd ed. Beijing: China Communication Press, 2008: 51 − 92. (in Chinese)
    [6] 周履. 关于薄壁箱形梁的悬臂翼缘板在翘曲扭转中作用的探讨[J]. 世界桥梁, 2003(3): 28 − 33. doi:  10.3969/j.issn.1671-7767.2003.03.009

    Zhou Lü. Exploration of contribution of cantilever flanges of thin-walled box girder to its behavior in warping torsion [J]. World Bridges, 2003(3): 28 − 33. (in Chinese) doi:  10.3969/j.issn.1671-7767.2003.03.009
    [7] 张元海, 林丽霞. 偏心轴向荷载作用下薄壁箱梁的约束扭转翘曲应力研究[J]. 工程力学, 2013, 30(3): 32 − 36.

    Zhang Yuanhai, Lin Lixia. Research on warping stress of restrained torsion for thin-walled box girder under eccentric axial load [J]. Engineering Mechanics, 2013, 30(3): 32 − 36. (in Chinese)
    [8] 马俊军, 蔺鹏臻. 考虑悬臂板贡献的箱梁约束扭转效应研究[J]. 铁道建筑, 2018, 58(3): 10 − 14. doi:  10.3969/j.issn.1003-1995.2018.03.03

    Ma Junjun, Lin Pengzhen. Study on restrained torsion effect for box girder considering contribution of cantilever slab [J]. Railway Engineering, 2018, 58(3): 10 − 14. (in Chinese) doi:  10.3969/j.issn.1003-1995.2018.03.03
    [9] 聂国隽, 钱若军. 考虑约束扭转的薄壁梁单元刚度矩阵[J]. 计算力学学报, 2002, 19(3): 344 − 348. doi:  10.3969/j.issn.1007-4708.2002.03.018

    Nie Guojuan, Qian Ruojun. Element stiffness matrix of thin-walled beam considering restrained torsion [J]. Chinese Journal of Computational Mechanics, 2002, 19(3): 344 − 348. (in Chinese) doi:  10.3969/j.issn.1007-4708.2002.03.018
    [10] Dikaros I C, Sapountzakis E J, Argyridi A K. Generalized warping effect in the dynamic analysis of beams of arbitrary cross section [J]. Journal of Sound and Vibration, 2016, 369: 119 − 146. doi:  10.1016/j.jsv.2016.01.022
    [11] Qiao P, Di J, Qin F J. Warping torsional and distorsional stress of composited box girder with corrugated steel webs [J]. Mathematical Problems in Engineering, 2018, 2018(6): 1 − 13.
    [12] Cambronero-Barrientos F, Díaz-del-Valle J, Martínez-Martínez J A. Beam element for thin-walled beams with torsion, distortion, and shear [J]. Engineering Structures, 2017, 143: 571 − 588. doi:  10.1016/j.engstruct.2017.04.020
    [13] 胡启平, 涂佳黄, 梁经群. 组合断面薄壁杆件弯扭耦合分析[J]. 工程力学, 2010, 27(7): 52 − 55.

    Hu Qiping, Tu Jiahuang, Liang Jingqun. Bending-torsion coupling analysis of the composite section thin-walled bar structure [J]. Engineering Mechanics, 2010, 27(7): 52 − 55. (in Chinese)
    [14] 王晓峰, 杨庆山. 基于Timoshenko梁理论的薄壁梁弯扭耦合分析[J]. 工程力学, 2008, 25(5): 12 − 16, 21.

    Wang Xiaofeng, Yang Qingshan. Coupled bend and torsion analysis of the spatial thin-walled beam using Timoshenko theory [J]. Engineering Mechanics, 2008, 25(5): 12 − 16, 21. (in Chinese)
    [15] Sapountzakis E J, Mokos V G. Warping shear stresses in non-uniform torsion by BEM [J]. Computational Mechanics, 2003, 30: 131 − 142. doi:  10.1007/s00466-002-0373-4
    [16] Günay M G, Timarci T. Static analysis of thin-walled laminated composite closed-section beams with variable stiffness [J]. Composite Structures, 2017, 182: 67 − 78. doi:  10.1016/j.compstruct.2017.08.092
    [17] Wan H X, Mahendran M. Behaviour and strength of hollow flange channel sections under torsion and bending [J]. Thin-Walled Structures, 2015, 94: 612 − 623. doi:  10.1016/j.tws.2015.05.013
    [18] 陈伯真, 胡毓仁. 薄壁结构力学 [M]. 上海: 上海交通大学出版社, 1998: 65 − 67.

    Chen Bozhen, Hu Yuren. Mechanics of thin-walled structures [M]. Shanghai: Shanghai Jiaotong University Press, 1998: 65 − 67. (in Chinese)
    [19] Kindmann R, Kraus M. Steel structures-design using FEM [M]. Berlin: Ernst & Sohn, 2011: 261 − 263.
  • [1] 魏亚, 梁思明, 吴泽弘.  徐变对水泥路面板湿度翘曲及应力影响研究 . 工程力学, 2017, 34(10): 106-115. doi: 10.6052/j.issn.1000-4750.2016.05.0397
    [2] 王丽娟, 胡昌斌, 曾宇鑫.  水泥混凝土路面板早龄期翘曲行为数值分析研究 . 工程力学, 2017, 34(7): 146-155. doi: 10.6052/j.issn.1000-4750.2016.02.0121
    [3] 张清东, 卢兴福, 张晓峰.  具有初始翘曲缺陷冷轧薄带钢板形瓢曲变形行为研究 . 工程力学, 2014, 31(8): 243-249,249. doi: 10.6052/j.issn.1000-4750.2013.04.0340
    [4] 樊春雷, 郝际平, 田炜烽.  H型截面钢梁约束扭转的简化计算 . 工程力学, 2013, 30(11): 180-184,220. doi: 10.6052/j.issn.1000-4750.2012.07.0529
    [5] 张元海, 林丽霞.  偏心轴向荷载作用下薄壁箱梁的约束扭转翘曲应力研究 . 工程力学, 2013, 30(3): 32-36. doi: 10.6052/j.issn.1000-4750.2011.10.0719
    [6] 张元海, 林丽霞.  一种考虑剪滞效应的斜支承连续箱梁挠曲扭转分析方法 . 工程力学, 2012, 29(2): 94-100.
    [7] 童根树, 王继涛, 张 磊.  箱形柱与工字梁连接节点翘曲变形的传递分析 . 工程力学, 2011, 28(6): 96-102,.
    [8] 刘计良, 张学森, 王正中, 韩彦宝.  支座约束及荷载分布对薄壁深梁应力的影响 . 工程力学, 2011, 28(8): 23-029.
    [9] 蔺鹏臻, 周世军, 刘凤奎.  抛物线型剪滞翘曲位移函数引起的附加轴力分析 . 工程力学, 2010, 27(8): 90-093,.
    [10] 刘计良, 王正中, 韩彦宝, 张学森.  均布荷载作用下工字形截面单跨超静定深梁应力计算方法研究 . 工程力学, 2010, 27(03): 174-179.
    [11] 夏 飞, 辛克贵.  薄壁曲梁振动分析的样条有限杆元法 . 工程力学, 2009, 26(3): 1-005.
    [12] 张元海, 李 乔.  曲线箱梁考虑剪滞效应时的弯扭分析 . 工程力学, 2009, 26(10): 123-129.
    [13] 徐 翔, 郝际平.  从Levinson高阶梁理论的一致变分到高次翘曲梁理论 . 工程力学, 2008, 25(2): 0-061.
    [14] 甘亚南, 周广春.  薄壁箱梁纵向剪滞翘曲函数精度选择的研究 . 工程力学, 2008, 25(6): 0-106.
    [15] 安 新, 张效松, 苏 波, 张锦辉.  薄壁杆件翘曲剪应力的边界元精确积分解法 . 工程力学, 2008, 25(1): 0-096.
    [16] 王正中, 朱军祚, 谌 磊, 郭佳陇, 谭东岳, 米文静.  集中力作用下深梁弯剪耦合变形应力计算方法 . 工程力学, 2008, 25(4): 0-120.
    [17] 谈梅兰, 王鑫伟, 吴光.  三维曲井内钻柱的接触非线性有限元分析 . 工程力学, 2006, 23(6): 162-166.
    [18] 胡少伟, 聂建国, 朱林森.  钢-混凝土组合梁复合弯扭作用下非线性分析 . 工程力学, 2005, 22(2): 1-5,26.
    [19] 段海娟, 张其林.  考虑翘曲效应的薄壁曲梁几何非线性分析 . 工程力学, 2004, 21(5): 157-160,.
    [20] 鲍永方, 黄文彬.  矩形箱梁约束扭转理论的分析与比较 . 工程力学, 1997, 14(3): 132-137.
  • 加载中
图(7) / 表 (3)
计量
  • 文章访问数:  66
  • HTML全文浏览量:  47
  • PDF下载量:  9
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-11-02
  • 修回日期:  2020-01-14
  • 网络出版日期:  2020-09-07
  • 刊出日期:  2020-09-25

基于协调翘曲场的开闭口混合薄壁截面杆件约束扭转一维有限元分析

doi: 10.6052/j.issn.1000-4750.2019.11.0639
    基金项目:  国家自然科学基金项目(51478475,51108460);中国博士后科学基金项目(2012M511759)
    作者简介:

    陈泽林(1995−),男,陕西人,硕士生,主要从事薄壁结构分析研究(E-mail: 13054170084@163.com)

    通讯作者: 文 颖(1981−),男,湖南人,副教授,博士,重载铁路工程结构教育部重点实验室副主任,主要从事桥梁稳定极限承载力及车桥系统振动稳定性研究(E-mail: ywen_ce@csu.edu.cn)
  • 中图分类号: TU323

摘要: 开口及闭口薄壁杆件约束扭转问题已由经典Timoshenko和Benscoter理论解决。然而,开闭口混合薄壁截面杆件约束扭转分析必须考虑开、闭口部分翘曲能力的差异,翘曲剪流形成机理有待进一步研究。该文假定开、闭口截面翘曲分别满足Vlasov和Umanskii假定,考虑开、闭口截面公共节点翘曲连续性要求,建立含有待定翘曲参数的协调翘曲模型。由截面受力平衡,确定翘曲参数显式列式,提出开闭口混合薄壁截面杆件约束扭转分析的一维有限元模型。算例及参数分析结果表明,基于Umanskii第二理论的I类方法在悬臂板及闭口周边引入附加剪流,影响翘曲剪应力精度。基于Umanskii第二理论的II类方法只能计算截面板件平均剪应力,无法反映真实翘曲剪流分布。基于Vlasov约束扭转假定的Beam-189单元忽略闭口周边约束效应产生的附加翘曲及剪流,影响翘曲正应力和剪应力精度。该文方法与Shell-63单元能得到基本吻合的变形与应力结果,说明一维梁元能正确反映开闭口混合薄壁截面杆件约束扭转及翘曲刚度。

English Abstract

文颖, 陈泽林. 基于协调翘曲场的开闭口混合薄壁截面杆件约束扭转一维有限元分析[J]. 工程力学, 2020, 37(9): 38-49. doi: 10.6052/j.issn.1000-4750.2019.11.0639
引用本文: 文颖, 陈泽林. 基于协调翘曲场的开闭口混合薄壁截面杆件约束扭转一维有限元分析[J]. 工程力学, 2020, 37(9): 38-49. doi: 10.6052/j.issn.1000-4750.2019.11.0639
Ying WEN, Ze-lin CHEN. ONE-DIMENSIONAL FINITE ELEMENT ANALYSIS OF WARPING TORSION FOR THIN-WALLED MEMBERS WITH OPEN-CLOSED CROSS SECTIONS BASED ON COMPATIBLE WARPING FIELD[J]. Engineering Mechanics, 2020, 37(9): 38-49. doi: 10.6052/j.issn.1000-4750.2019.11.0639
Citation: Ying WEN, Ze-lin CHEN. ONE-DIMENSIONAL FINITE ELEMENT ANALYSIS OF WARPING TORSION FOR THIN-WALLED MEMBERS WITH OPEN-CLOSED CROSS SECTIONS BASED ON COMPATIBLE WARPING FIELD[J]. Engineering Mechanics, 2020, 37(9): 38-49. doi: 10.6052/j.issn.1000-4750.2019.11.0639
  • 桥梁工程中广泛应用的箱形梁属于开闭口混合薄壁截面杆件,封闭箱体能提供显著双向弯曲和扭转刚度,适应纵、横向偏载作用,开口悬臂板则为桥面结构提供良好支撑,满足桥梁功能要求。箱梁设计必须解决好空间弯扭分析问题。开口及全封闭薄壁截面杆件约束扭转分析分别由基于Vlasov假定的Timoshenko理论和基于Umanskii假定的Benscoter理论解决[1-3]。研究发现开、闭口断面翘曲刚度差异明显。因此,箱梁约束扭转分析必须考虑开口悬臂板和闭口周边翘曲能力的区别[4]

    具有任意截面形状的薄壁杆件扭转翘曲位移一般可写成截面翘曲函数和翘曲参数的乘积形式。现有文献认为截面扇性坐标和广义扇性坐标分别适用于开口和闭口周边的翘曲函数,反映二者翘曲能力的差异,而描述翘曲沿杆件纵向变化的翘曲参数却是一致的。徐勋等[4]、郭金琼等[5]、周履[6]、张元海和林丽霞[7]、马俊军和蔺鹏臻[8]、聂国隽和钱若军[9]、Dikaros等[10]和Qiao等[11]采用Umanskii第二理论,即开、闭口部分采用相同的待定翘曲参数,进行薄壁构件约束扭转分析。但是,关于翘曲剪流沿开、闭口部分如何分布仍存在分歧,例如文献[5, 7-8, 12]认为开口部分参与闭口周边翘曲剪流重分配,导致悬臂板产生附加剪流,而文献[4, 6]则从悬臂板自由端剪应力为零的条件出发,认为闭口周边将产生平衡悬臂板翘曲正应力的附加翘曲剪流。胡启平等[13]和王晓峰和杨庆山[14]认为开口截面翘曲剪流将引起附加翘曲。产生分歧的原因在于采用Umanskii假定后,悬臂板翘曲剪流究竟由剪应变还是根据与翘曲正应力间的平衡条件来决定。考虑到闭口截面自由扭转刚度远大于约束扭转刚度,Sapountzakis和Mokos[15]、Günay和Timarci[16]及Cambronero-Barrientos等[12]采用Umanskii第一理论,假定开、闭口部分约束扭转翘曲参数取为截面扭转角变化率。然而,它却无法反映翘曲剪应力引起附加翘曲,将产生不可允许的误差[5]。综上所述,现有文献在开、闭口部分约束扭转翘曲参数如何确定的问题上并未达成一致,势必影响约束扭转分析可靠性。

    考虑到薄壁箱梁悬臂板宽度有限,且在闭口周边面内弯翘连续性条件影响下,板件中面内形成附加翘曲剪流,本文假定截面开口部分翘曲参数等于扭转角变化率(考虑悬臂板面内弯翘,忽略剪翘),而闭口部分翘曲参数则满足Umanskii假定。考虑开、闭口截面公共节点翘曲连续性条件,建立包含待定翘曲参数的协调翘曲模型。基于截面内外力平衡条件,建立闭口部分待定翘曲参数与刚性扭转角导数间的显式关系。依据势能不变值原理,建立开闭口混合薄壁截面杆件约束扭转分析的一维有限元列式。通过数值算例分析并与Umanskii第二理论、Shell-63单元、基于Vlasov假定的Beam-189单元结果进行比较,论证本文方法计算翘曲位移和应力具有良好的精度。通过变化悬臂板宽及梁高,揭示不同约束扭转分析模型计算开、闭口截面翘曲刚度变化条件下杆件翘曲应力变化规律。

    • 开闭口混合薄壁截面杆件位移场的建立基于如下假定:

      1) 截面发生整体轴向位移和双向弯曲转动后仍然保持为平面;

      2) 截面面内位移满足刚性周边假定,即忽略畸变和横向弯曲的影响;

      3) 截面闭口周边约束扭转翘曲函数采用以剪心为极点的广义主扇性坐标,翘曲参数取为待定函数η(x);

      4) 截面开口部分约束扭转翘曲函数采用以剪心为极点的主扇性坐标,翘曲参数取为扭转角φ的变化率dφ/dx

      首先建立描述截面空间位移的笛卡儿坐标系Oxyz(如图1所示),它以截面形心O为坐标原点。由于箱梁截面一般关于Oy轴左右对称,则OyOz为截面惯性主轴。取S点(αy, 0)为截面剪心。由假定2可知,截面上任意点P(y, z)的面内位移(vp, wp)可用S点沿yz方向的位移分量(vs, ws)以及截面绕S点刚性转角φ来表示,故有:

      $$\begin{split}& {v_{\rm {p}}} = {v_{\rm {s}}} - {\textit z}\varphi ,\\& {w_{\rm {p}}} = {w_{\rm {s}}} + ( {y - {\alpha _{\rm {y}}}} )\varphi \end{split} $$ (1)

      图  1  开闭口混合截面布置

      Figure 1.  Cross sectional layout for open-closed profile

      假设P点沿截面中线切线方向(称为s轴)位移为$v_{\rm{p}}^{\rm{s}}$,并令s轴对z轴的倾角为α(以自z轴按右手法则转到s轴为正),则有:

      $$v_{\rm {p}}^{\rm{s}} = - {v_{\rm {p}}}\sin \alpha + {w_{\rm {p}}}\cos \alpha $$ (2)

      将式(1)代入式(2),可得:

      $$\left\{ {\begin{aligned} & {v_{\rm {p}}^{\rm{s}} = - {v_{\rm {s}}}\sin \alpha + {w_{\rm {s}}}\cos \alpha + {\rho _{\rm {s}}}\varphi}\\& {{\rho _{\rm {s}}} = {\textit{z}}\sin \alpha + ( {y - {\alpha _{\rm {y}}}} )\cos \alpha} \end{aligned}} \right. $$ (3)

      式中,ρs表示剪心S到过Ps轴的距离。

      首先考虑开闭口混合截面发生自由扭转,此时开口部分的中面剪流为零,而闭口部分中面存在自由扭转Bredt剪流qf。因此,由剪应变定义可得:

      对于闭口部分:

      $$\frac{{{q_{\rm {f}}}}}{{Gt}} = \frac{{\partial {u_{{\rm {p}}0}}}}{{\partial s}} + \frac{{\partial v_{{\rm {p}}0}^{\rm{s}}}}{{\partial x}} $$ (4)

      其中:

      $$\left\{ {\begin{aligned} & {v_{{\rm {p}}0}^{\rm{s}} = - {v_{\rm {s}}}\sin \alpha + {w_{\rm {s}}}\cos \alpha + {\rho _{{\rm {s}}0}}\varphi}\\& {{\rho _{{\rm {s}}0}} = {{\textit{z}}_{\rm {s}}}\sin \alpha + ( {{y_{\rm {s}}} - {\alpha _{\rm {y}}}} )\cos \alpha} \end{aligned}} \right. $$ (5)

      式中:(up0, $v_{{\rm{p}0}} ^{\rm{s}}$)表示与P点对应的截面中面(ys, zs)位移;ρs0表示剪心S到位于截面中面s轴的距离。取Oy轴与顶板中面交点⑦(y0, z0)为s轴原点,则有:

      $$\begin{split}& {u_{\rm {p}}}\! =\! {u_{{\rm {p}}0}} \!= \!{u_0} - \frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}}\int_0^s {\frac{{{\rm d}y}}{{{\rm d}s}}{\rm d}s} - \frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}}\int_0^s {\frac{{{\rm d}{\textit{z}}}}{{{\rm d}s}}{\rm d}s} - \frac{{{\rm d}\varphi }}{{{\rm d}x}}\varpi = \\& \quad \, {u_0} - \frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}}\left( {{y_{\rm {s}}} - {y_0}} \right) - \frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}}\left( {{{\textit{z}}_{\rm {s}}} - {{\textit{z}}_0}} \right) - \frac{{{\rm d}\varphi }}{{{\rm d}x}}\varpi \qquad\qquad(6) \end{split}$$

      式中,u0表示s轴原点的轴向翘曲。式(6)假定翘曲沿壁厚均布,即忽略截面次翘曲。整理式(6)得:

      $${u_{\rm {p}}} = {\overline u _0} - {y_{\rm {s}}}\frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}} - {{\textit{z}}_{\rm {s}}}\frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}} - \frac{{{\rm d}\varphi }}{{{\rm d}x}}\varpi $$ (7)

      其中:

      $$\left\{ {\begin{aligned} & {{\overline u _0} = {u_0} + \frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}}{y_0} + \frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}}{{\textit{z}}_0}}\\& {\varpi = \int_0^s {\left( {{\rho _{{\rm {s}}0}} - \frac{\Omega }{{t\displaystyle\oint {\dfrac{{{\rm d}s}}{t}} }}} \right){\rm d}s}} \end{aligned}} \right.$$ (8)

      式中,$\varpi $是以剪心为极点广义主扇性坐标[5]

      对于开口部分:

      $$\frac{{\partial {u_{{\rm {p}}0}}}}{{\partial s}} + \frac{{\partial v_{{\rm {p}}0}^{\rm{s}}}}{{\partial x}} = 0$$ (9)

      取悬壁板与顶板中面交点②(y2, z2)或⑤(y5, z5)为积分起点,则有:

      $${u_{\rm {p}}} = {u_{2(5)}} - \frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}}\left( {{y_{\rm {s}}} - {y_{2(5)}}} \right) - \frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}}\left( {{{\textit{z}}_{\rm {s}}} - {{\textit{z}}_{2(5)}}} \right) - \frac{{{\rm d}\varphi }}{{{\rm d}x}}\omega $$ (10)

      其中:

      $$\omega = \int_0^s {{\rho _{{\rm {s}}0}}{\rm d}s} $$ (11)

      式中,u2(5)表示积分起点的轴向翘曲。尽管式(10)忽略次翘曲,其对自由扭转剪应力贡献将在下文考虑。由式(7)可得:

      $${u_{2(5)}} = {\overline u _0} - {y_{2(5)}}\frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}} - {{\textit{z}}_{2(5)}}\frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}} - \frac{{{\rm d}\varphi }}{{{\rm d}x}}{\varpi _{2(5)}}$$ (12)

      将式(12)代入式(10),得到:

      $${u_{\rm {p}}} = {\overline u _0} - \frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}}{y_{\rm {s}}} - \frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}}{{\textit{z}}_{\rm {s}}} - \frac{{{\rm d}\varphi }}{{{\rm d}x}}\left( {{\varpi _{2(5)}} + \omega } \right)$$ (13)

      不难发现,式(7)和式(13)的区别在于开、闭口部分主扇性坐标不同。图1所示的开闭口混合截面主扇性坐标分布如图2所示。

      图  2  开闭口混合截面主扇性坐标

      Figure 2.  Principal sectorial coordinates for open-closed profile

      当开闭口混合截面发生约束扭转时,由假定3可知,闭口周边任意点翘曲位移为:

      $${u_{\rm {p}}} = {\overline u _0} - {y_{\rm {s}}}\frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}} - {{\textit{z}}_{\rm {s}}}\frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}} - \eta \varpi $$ (14)

      则悬壁板与顶板交点②(y2, z2)或⑤(y5, z5)的翘曲为:

      $${u_{2(5)}} = {\overline u _0} - {y_{2(5)}}\frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}} - {{\textit{z}}_{2(5)}}\frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}} - \eta {\varpi _{2(5)}}$$ (15)

      基于假定4,并将式(15)代入式(10),得到开口截面任意点翘曲位移为:

      $${u_{\rm {p}}} = {\overline u _0} - \frac{{{\rm d}{v_{\rm {s}}}}}{{{\rm d}x}}{y_{\rm {s}}} - \frac{{{\rm d}{w_{\rm {s}}}}}{{{\rm d}x}}{{\textit{z}}_{\rm {s}}} - \eta {\varpi _{2(5)}} - \frac{{{\rm d}\varphi }}{{{\rm d}x}}\omega $$ (16)

      不难发现,由式(14)和式(16)表示的开闭口混合截面翘曲场是连续的。

    • 由式(14)和式(16)可知,开闭口混合截面任意点翘曲应力的表达式为:

      对于闭口部分:

      $${\sigma _{\rm {p}}} = E\left( {{{\overline u }_0^{\prime}} - {y_{\rm {s}}}{v_{\rm {s}}^{\prime \prime }} - {{\textit{z}}_{\rm {s}}}{w_{\rm {s}}^{\prime \prime }} - \eta '\varpi } \right)$$ (17)

      对于开口部分:

      $${\sigma _{\rm {p}}} = E\left( {{{\overline u }_0^\prime} - {y_{\rm {s}}}{v_{\rm {s}}^{\prime \prime }} - {{\textit{z}}_{\rm {s}}}{w_{\rm {s}}^{\prime \prime }} - {\varpi _{2(5)}}\eta ' - \omega \varphi ''} \right)$$ (18)

      式中,撇号(′)代表对坐标x求导。由图2及式(17)和式(18)可知,${\overline u _0'}$表示截面平均翘曲应变。假设开闭口混合截面杆件因承受外扭矩作用发生约束扭转,则杆件剪心轴不挠曲($v_{\rm{s}}'' $=0, $w_{\rm{s}}'' $=0),由于无轴力作用,截面平均翘曲应变为零(${\overline u _0'}$=0)。因此,式(17)和式(18)变为:

      对于闭口部分:

      $${\sigma _{\rm {p}}} = - E\eta '\varpi $$ (19)

      对于开口部分:

      $${\sigma _{\rm {p}}} = - E\left( {{\varpi _{2(5)}}\eta ' + \omega \varphi ''} \right)$$ (20)

      式(19)和式(20)中,只有翘曲参数η缺乏明确物理概念。为了建立ηφ的关系,考虑杆件截面受力平衡关系如下:

      $$\left\{ {\begin{aligned} & {N = \int_A {{\sigma _{\rm {p}}}t{\rm{d}}s} = 0}\\& {{M_{\rm{z}}} = - \int_A {{\sigma _{\rm {p}}}{y_{\rm {s}}}t{\rm{d}}s} = 0}\\& {{M_{\rm {y}}} = \int_A {{\sigma _{\rm {p}}}{{\textit{z}} _{\rm {s}}}t{\rm{d}}s} = 0} \end{aligned}} \right. $$ (21)

      式中,A表示全截面积。将式(19)和式(20)代入式(21),并考虑图2所示的主扇性坐标分布,可知式(21)前两项积分恒为零。由第3项积分可得:

      $$\eta '\left( {\oint_{{A_{\rm c}}} {\varpi {{\textit{z}}_{\rm {s}}}t{\rm{d}}s} + \int_{{A_{\rm o}}} {{\varpi _{2(5)}}{{\textit{z}}_{\rm {s}}}t{\rm{d}}s} } \right) + \varphi ''\int_{{A_{\rm{o}}}} {\omega {{\textit{z}}_{\rm {s}}}t{\rm{d}}s} = 0$$ (22)

      考虑自由扭转下开闭口混合截面受力平衡,故有:

      $$\int_A {\varpi {{\textit{z}}_{\rm {s}}}{\rm{d}}A} = \oint_{{A_{\rm c}}} {\varpi {{\textit{z}}_{\rm {s}}}t{\rm{d}}s} + \int_{{A_{\rm o}}} {\left( {{\varpi _{2(5)}} + \omega } \right){{\textit{z}}_{\rm {s}}}t{\rm{d}}s} = 0$$ (23)

      式中,AcAo分别表示闭口和开口部分截面积。将式(23)代入式(22),得:

      $$\eta ' = \varphi ''$$ (24)

      由于约束扭转条件下,截面闭口部分剪应力包含沿壁厚线性分布的圣维南剪应力τs、沿壁厚均匀分布的Bredt剪应力τB及沿壁厚均匀分布的翘曲剪应力τωτs相比τB对截面扭矩的贡献很小,忽略τs对截面内外扭矩平衡的影响。沿壁厚均匀分布的τBτω可根据剪应变及材料弹性本构关系得到,故有:

      $${\tau _{\rm {B}}} + {\tau _{\text{ω}} } = G\gamma _{\rm yz}^{\rm c},\quad \gamma _{\rm yz}^{\rm c} = \frac{{\partial {u_{{\rm {p}}0}}}}{{\partial s}} + \frac{{\partial v_{{\rm {p}}0}^{\rm{s}}}}{{\partial x}}$$ (25)

      将式(5)和式(14)代入式(25),可得:

      $${\tau _{\rm {B}}} + {\tau _{\text{ω}} } = G\left[ { - \eta \left( {{\rho _{{\rm {s}}0}} - \frac{\Omega }{{t\displaystyle \oint {\dfrac{{{\rm d}s}}{t}} }}} \right) + {\rho _{{\rm {s}}0}}\varphi '} \right]$$ (26)

      因此,闭口周边中面剪应力对扭心的合力矩为:

      $$\begin{split}& {M^{\rm c}} = \oint {\left( {{\tau _{\rm {B}}} + {\tau _{\text{ω}} }} \right){\rho _{{\rm {s}}0}}t{\rm{d}}s} = \\[-5pt]& \quad \;\; G\oint {{\rho _{{\rm {s}}0}}^2t{\rm{d}}s} \left( {\varphi ' - \eta } \right) + G\eta \frac{{{\Omega ^2}}}{{\displaystyle \oint {\dfrac{{{\rm{d}}s}}{t}} }} = \\[-5pt]& \quad \;\;G{I_{{\rm {pc}}}}\left( {\varphi ' - \eta } \right) + G{J_{\rm {B}}}\eta \end{split}$$ (27)

      其中:

      $${I_{{\rm {pc}}}} = \oint {{\rho _{{\rm {s}}0}}^2t {\rm{d}} s} ,\quad {J_{\rm {B}}} = \frac{{{\Omega ^2}}}{{\displaystyle \oint {\dfrac{{ {\rm{d}} s}}{t}} }}$$ (28)

      式中:Ipc表示闭口部分极常数[3]JB表示闭口部分抗扭惯性矩。

      形成闭口部分扭矩Mc的剪应力(τB+τω)还可由与翘曲正应力σp的平衡关系而得到[5],故有:

      $${\tau _{\rm {B}}} + {\tau _{\text{ω}} } = \frac{1}{t}\left( {{\tau _0}{t_0}} \right) - \frac{1}{t}\int_0^P {\frac{{\partial {\sigma _{\rm {p}}}}}{{\partial x}}t {\rm{d}} s} $$ (29)

      式中,τ0t0表示位于s轴起点的剪力流。将式(29)代入式(27),可得:

      $$\begin{split} {M^{\rm{c}}} = & \oint {\left( {{\tau _{\rm {B}}} + {\tau _{\text{ω}} }} \right){\rho _{{\rm {s}}0}}t {\rm{d}} s} = \\& {\tau _0}{t_0}\Omega + E\eta ''\left( {\oint {{S_\varpi }{\rho _{{\rm {s}}0}} {\rm{d}} s} + \underline {{\varpi _5}{A_{5 - 6}} \cdot {\rm{d}} {\omega _{5 - 2}}} } \right) + \\& \underline {E\varphi ''' {\rm{d}} {S_{{\text{ω}} 5 - 6}} \cdot {\rm{d}} {\omega _{5 - 2}}} \\[-12pt] \end{split} $$ (30)

      其中:

      $$\left\{ {\begin{aligned} & {{A_{5 - 6}} = \int_5^6 {t {\rm{d}} s} ,\; {\rm{d}} {\omega _{5 - 2}} = \int_5^2 {{\rho _{{\rm {s}}0}} {\rm{d}} s}}\\& {{S_\varpi } = \int_0^P {\varpi t {\rm{d}} s} ,\;{\rm{d}}{S_{\omega 5 - 6}} = \int_5^6 {\omega t {\rm{d}} s}} \end{aligned}} \right. $$ (31)

      式中,积分上、下限数字表示如图1所示截面上的关键点。带下划线项表示悬臂板翘曲正应力引起闭口周边附加剪流形成的翘曲扭矩。联立式(25)、式(29)及式(30),并考虑截面板件面内弯翘连续性条件:

      $$\oint {\frac{{\partial {u_{{\rm {p}}0}}}}{{\partial s}} {\rm{d}} s} = 0$$ (32)

      可得:

      $$\begin{split} {M^{\rm{c}}} = &G{J_{\rm {B}}}\varphi ' - E{I_{{\text{ϖ}} {\rm c}}}\eta '' + \\& \underline {E\left( {{\varpi _2} - {\varpi _5}} \right)\left( {\eta ''{\varpi _5}{A_{5 - 6}} + \varphi ''' {\rm{d}} {S_{{\text{ω}} 5 - 6}}} \right)} \end{split}$$ (33)

      其中:

      $${I_{{\text{ϖ}} {\rm{c}}}} = \oint_{{A_{\rm{c}}}} {{\varpi ^2}t {\rm{d}} s} $$ (34)

      需要注意的是,倘若忽略悬臂板对闭口部分翘曲剪流的影响,即略去式(33)中带下划线项,则得到纯闭口截面内外扭矩平衡条件。

      由式(27)和式(33)的恒等条件,并考虑式(24),得:

      $$\eta = \varphi ' - \frac{{E\left( {{I_{\text{ϖc}}} + {I_{\text{ωϖ}}}} \right)}}{{G\left( {{J_{\rm {B}}} - {I_{{\rm {pc}}}}} \right)}}\varphi '''$$ (35)

      其中:

      $${I_{\text{ϖω}}} = \int_{{A_{\rm{o}}}} {{\varpi _{2\left( 5 \right)}} \cdot( {{\varpi _{2\left( 5 \right)}} + \omega } )t{\rm{d}}s} $$ (36)

      式(35)给出了开闭口混合截面闭口部分翘曲参数与截面刚性扭转角导数间的关系。当悬臂板宽为零,即截面恢复纯闭口形态时,由式(22)可知,式(24)不再成立。联立式(27)和式(33)后,式(35)应修正为:

      $$\eta = \varphi ' - \frac{{E{I_{\text{ϖc}}}}}{{G( {{J_{\rm {B}}} - {I_{{\rm {pc}}}}} )}}\eta ''$$ (37)

      由式(27)可得:

      $$m = G{I_{{\rm {pc}}}}( {\eta ' - \varphi ''} ) - G{J_{\rm {B}}}\eta '$$ (38)

      式中,m(=−dMc/dx)表示单位长度外扭矩。将式(37)代入式(38),可得:

      $$E{I_{\text{ϖc}}}\eta ''' - vG{J_{\rm {B}}}\eta ' = vm$$ (39)

      式中,v=1−JB/Ipc,即为文献[5, 17]中的截面约束系数。式(39)表示经典Umanskii第二理论给出的闭口截面杆件约束扭转翘曲微分方程。

    • 由前两节讨论可知,式(1)、式(14)和式(16)所描述的开闭口混合截面任意点空间位移依赖于截面平均翘曲位移${\overline u _0}$、剪心线位移(vs,ws)和绕剪心扭转角φ。考虑如图3所示开闭口混合薄壁截面杆件单元ij,假定杆件截面平均轴向位移沿杆件长度方向线性变化,杆件剪心轴线挠曲及扭转角沿杆长方向按Hermitian三次多项式变化。因此,杆件任意截面独立位移参数u与节点自由度d关系如下:

      $${{u}} = {{Cd}}$$ (40)

      其中:

      $$\!\left\{ {\begin{aligned} & {{{u}} = {\{ {\begin{array}{*{20}{c}} {{{\overline u }_0}}&{{v_{\rm {s}}}}&{{w_{\rm {s}}}}&\varphi \end{array}} \}^{\rm{T}}}}\\& {{{d}} = {\{ {\begin{array}{*{20}{c}} {{{{d}}_i}}&{{{{d}}_j}} \end{array}} \}^{\rm{T}}}}\\& {{{{d}}_i} = {\{ {\begin{array}{*{20}{c}} {{{\overline u }_{0i}}}&{{v_{{\rm {s}}i}}}&{{w_{{\rm {s}}i}}}&{{\varphi _i}}&{{\varphi _{{\rm y}i}}}&{{\varphi _{{\rm z}i}}}&{{{\varphi _i'}}} \end{array}}\}^{\rm{T}}}}\\& {{{{d}}_j} = {\{ {\begin{array}{*{20}{c}} {{{\overline u }_{0j}}}&{{v_{{\rm {s}}j}}}&{{w_{{\rm {s}}j}}}&{{\varphi _j}}&{{\varphi _{{\rm y}j}}}&{{\varphi _{{\rm z}j}}}&{{{\varphi _j'}}} \end{array}} \}^{\rm{T}}}}\\& {{{C}} \!\!=\! \!\left[\!\!\! \!{\begin{array}{*{20}{c}} {{N_1}{{{I}}_1}}&{{0}}&{{0}}&{{N_4}{{{I}}_1}}&{{0}}&{{0}} \\ {{0}}&{{N_2}{{{I}}_3}}&{{N_3}{{J}}}&{{0}}&{{N_5}{{{I}}_3}}&{{N_6}{{J}}} \end{array}} \!\!\!\!\right]}\\& {{{J}} \!\!=\!\! \left[\!\!\! {\begin{array}{*{20}{c}} 0&1&0 \\ { - 1}&0&0 \\ 0&0&1 \end{array}} \!\!\!\right],\;{N_1} \!= \!1 \!-\! \xi ,\;{N_2} \!= \!1\! -\! 3{\xi ^2} \!+\! 2{\xi ^3}}\\& {{N_3} = x( {1 - 2\xi + {\xi ^2}} ),\;{N_4} = \xi ,\;{N_5} = 3{\xi ^2} - 2{\xi ^3}}\\& {{N_6} = x( { - \xi + {\xi ^2}} ),\;\xi = \frac{x}{l}} \end{aligned}} \right.$$ (41)

      式中:Ii表示第i阶单位矩阵;l表示单元长度。由式(3)、式(5)、式(14)、式(16)、式(24)和式(35),得截面应力场。

      对于闭口部分:

      $${{\sigma }} = {{DBu}}$$ (42)

      其中:

      $$\left\{ {\begin{aligned} & {{{\sigma }} = {\{ {\begin{array}{*{20}{c}} {{\sigma _{\rm {p}}}}&{{\tau _{\rm {p}}}} \end{array}} \}^{\rm{T}}},\;{{D}} = \left[ {\begin{array}{*{20}{c}} E&0 \\ 0&G \end{array}} \right]}\\[-2pt]& {{{B}} = {{{B}}_{{\rm c}0}}{{{B}}_{{\rm c}1}}}\\[-2pt]& {{{{B}}_{{\rm c}0}} = \left[ {\begin{array}{*{20}{c}} {{{B}}_0^1}&{{0}} \\ {{0}}&{{{B}}_0^2} \end{array}} \right]}\\[-2pt]& {{{B}}_0^1 = \{ {\begin{array}{*{20}{c}} 1&{ - {y_{\rm {s}}}}&{ - {{\textit{z}} _{\rm {s}}}}&{ - \varpi } \end{array}} \}}\\[-2pt]& {{{B}}_0^2 = \left\{ {\begin{array}{*{20}{c}} {\dfrac{\Omega }{{t\oint {\dfrac{{{\rm d}s}}{t}} }}}&{\dfrac{{E( {{I_{\text{ϖc}}} + {I_{\text{ωϖ}}}} )}}{{G( {{J_{\rm {B}}} - {I_{{\rm {pc}}}}} )}} \cdot \dfrac{{{\rm d}\varpi }}{{{\rm d}s}}} \end{array}} \right\}}\\[-2pt]& {{{{B}}_{{\rm c}1}} = {\left[ {\begin{array}{*{20}{c}} {\dfrac{{\rm d}}{{{\rm d}x}}}&{{0}}&0&0 \\ {{0}}&{\dfrac{{{{\rm d}^2}}}{{{\rm d}{x^2}}}{{{I}}_3}}&{\dfrac{{\rm d}}{{{\rm d}x}}{{{e}}^{\rm{T}}}_3}&{\dfrac{{{{\rm d}^3}}}{{{\rm d}{x^3}}}{{{e}}^{\rm{T}}}_3} \end{array}} \right]^{\rm{T}}}}\\[-2pt]& {{{{e}}_3} = \{ {\begin{array}{*{20}{c}} 0&0&1 \end{array}} \}} \end{aligned}} \right.$$ (43)

      图  3  开闭口混合薄壁截面杆件单元

      Figure 3.  Thin-walled element with open-closed cross section

      对于开口部分:

      $${{\sigma }} = {{DBu}}$$ (44)

      其中:

      $$\left\{ {\begin{aligned} & {{{B}} = {{{B}}_{{\rm o}0}}{{{B}}_{{\rm o}1}}}\\[-2pt]& {{{{B}}_{{\rm o}0}} = \left[ \!\!\!{\begin{array}{*{20}{c}} 1&{ - {y_{\rm {s}}}}&{ - {{\textit{z}}_{\rm {s}}}}&{ - {\varpi _{2(5)}} - \omega }&0 \\ 0&0&0&0&{{\rho _{{\rm {s}}0}} - {\rho _{\rm {s}}}} \end{array}}\!\! \right]}\\[-2pt]& {{{{B}}_{{\rm o}1}} = {\left[ {\begin{array}{*{20}{c}} {\dfrac{{\rm d}}{{{\rm d}x}}}&{{0}}&0 \\ {{0}}&{\dfrac{{{{\rm d}^2}}}{{{\rm d}{x^2}}}{{{I}}_3}}&{\dfrac{d}{{{\rm d}x}}{{{e}}^{\rm{T}}}_3} \end{array}} \right]^{\rm{T}}}} \end{aligned}} \right. $$ (45)

      由虚功原理,可得开闭口混合薄壁截面杆件单元节点平衡条件如下:

      $$\delta U + \delta V = 0$$ (46)

      其中:

      $$\left\{ {\begin{aligned} & {\delta U = \delta {U_{\rm c}} + \delta {U_{\rm o}} }\\[-2pt]& {\delta {U_{\rm c}} = \delta {{{d}}^{\rm{T}}}\int_0^l {{{\left( {{{{B}}_{{\rm c}1}}{{C}}} \right)}^{\rm{T}}}\left( {\oint_{{A_{\rm{c}}}} {{{{B}}_{{\rm c}0}^{\rm{T}}}{{D}}{{{B}}_{{\rm{c}}0}} {\rm{d}} A} } \right){{{B}}_{{\rm{c}}1}}{{C}} {\rm{d}} x} {{d}}}\\[-2pt]& {\delta {U_{\rm{o}}} = \delta {{{d}}^{\rm{T}}}\int_0^l {{{\left( {{{{B}}_{{\rm{o}}1}}{{C}}} \right)}^{\rm{T}}}\left( {\int_{{A_{\rm{o}}}} {{{{B}}_{{\rm{o}}0}^{\rm{T}}}{{D}}{{{B}}_{{\rm{o}}0}} {\rm{d}} A} } \right){{{B}}_{{\rm{o}}1}}{{C}} {\rm{d}} x} {{d}}}\\[-2pt]& \delta V = - \sum\limits_i {F_{\rm{x}}^iu_{{\rm {p}}0}^i + F_{\rm {y}}^iv_{\rm {p}}^i + F_{\rm{z}}^iw_{\rm {p}}^i + M_{\rm{x}}^i{\varphi ^i}} - \\[-5pt]& \quad \quad \int_0^l {( {{q_{\rm{x}}}{u_{q0}} + {q_{\rm {y}}}{v_q} + {q_{\rm{z}}}{w_q} + {m_{\rm{x}}}\varphi } ) {\rm{d}} x} = - \delta {{{d}}^{\rm{T}}}{{F}} \end{aligned}} \right. $$ (47)

      式中:UcUo分别表示单元截面闭、开口部分应变能;V表示外荷载势能;F表示与单元节点位移对应的广义力。因此,单元平衡方程为:

      $$\delta {{{d}}^{\rm{T}}}{{kd}} = \delta {{{d}}^{\rm{T}}}{{F}}$$ (48)

      其中:

      $$\left\{ {\begin{aligned} & {{{k}} = {{{k}}_{\rm{c}}} + {{{k}}_{\rm{o}}} }\\& {{{{k}}_{\rm{c}}} = \int_0^l {{{\left( {{{{B}}_{{\rm{c}}1}}{{C}}} \right)}^{\rm{T}}}\left( {\oint_{{A_{\rm{c}}}} {{{{B}}_{{\rm{c}}0}^{\rm{T}}}{{D}}{{{B}}_{{\rm{c}}0}}{\rm{d}}A} } \right){{{B}}_{{\rm{c}}1}}{{C}}{\rm{d}}x}}\\& {{{{k}}_{\rm{o}}} = \int_0^l {{{\left( {{{{B}}_{{\rm{o}}1}}{{C}}} \right)}^{\rm{T}}}\left( {\int_{{A_{\rm{o}}}} {{{{B}}_{{\rm{o}}0}^{\rm{T}}}{{D}}{{{B}}_{{\rm{o}}0}}{\rm{d}}A} } \right){{{B}}_{{\rm{o}}1}}{{C}}{\rm{d}}x}} \end{aligned}} \right.$$ (49)

      式中,kcko分别表示闭、开口截面对单元刚度贡献。单元刚度矩阵k的显式表达式见附录。

    • 为了验证本文提出的开闭口混合薄壁截面杆件约束扭转分析方法和一维有限元刚度矩阵正确性,本节通过典型算例分析,并与ANSYS Beam-189单元、Shell-63单元及文献结果(基于Umanskii第二理论[5, 18])进行对比。要确保Shell-63单元结果为准确解,通过施加刚性扭转荷载及沿箱梁纵向增设横隔板,排除箱梁截面畸变、横向弯曲对应力分布的影响。

    • 考虑悬臂箱梁承受偏心轴向线荷载,该轴向均布荷载作用于箱梁截面②号点并指向ox轴负方向。箱梁基本几何及材料参数详见文献[7]。截面形心和剪心几何位置参数为yo=82.35 mm,αy=−6.72 mm。截面常数为A=6.80×103 mm2Iy=6.80×107 mm4Iz=4.46×107 mm4Iω=2.10×1010 mm6Ipc=5.63×107 mm4Iωc=8.58×109 mm6Iωω=−5.85×109 mm6J≈JB=5.49×107 mm4。 截面关键点主扇性坐标为ω1=−6386.22 mm2ω2=2521.18 mm2ω3=−335.97 mm2。本文取60个薄壁梁单元离散悬臂梁,运用805个Shell-63单元离散箱梁板件。为最大限度减小畸变对约束扭转翘曲及应力的影响,Shell-63模型沿纵向均匀增设8道横隔板(t=6 mm)。本文方法、Shell-63单元及文献[7]结果列于表1。不难发现,本文方法和板壳元结果吻合良好,二者均说明箱梁约束扭转对截面翘曲位移及正应力的影响不超过5%,而文献[7]认为箱梁约束扭转使得顶板与腹板交点正应力增大12%,悬臂板自由端正应力则减小30%。事实上,本例箱梁截面翘曲系数[19]$ \varepsilon_{\rm T}=l\sqrt{GJ_{\rm{B}}/EI_{\rm{w}}}$=40.1>10,箱梁约束扭转翘曲位移和正应力相比初等梁理论计算结果可以忽略。

      表 1  悬臂箱梁约束扭转位移和正应力

      Table 1.  Displacements and normal stresses for warping torsion of a cantilever box girder

      计算模型自由端扭转角/(×10−5 rad)自由端翘曲/(×10−5 m)/初等梁理论翘曲/(×10−5 m)固定端正应力/MPa/初等梁理论应力/MPa
      本文一维梁元−4.95−1.838/−1.84−10.72/−10.71−14.21/−14.24−0.091/−0.092−0.545/−0.536−0.687/−0.712
      Shell-63单元−5.05−1.7/−1.84−11.0/−10.71−14.0/−14.24−0.092/−0.092−0.54/−0.536−0.68/−0.712
      文献[7]−0.091/−0.092−0.598/−0.536−0.55/−0.712
      注:翘曲位移和应力的正负号与Ox坐标轴指向保持一致。
    • 文献[8]分析了均布扭矩(mT=23.5 kN·m/m)作用下30 m简支箱梁的约束扭转问题。箱梁基本几何、材料参数及截面主扇性坐标见文献[5]。截面形心和剪心位置参数为yo=0.955 m,αy=−0.296 m。与抗扭特性相关的截面常数包括:Iω=2.37 m6Ipc=9.85 m4Iωc=1.56 m6Iωω=−0.18 m6J=8.07 m4JB=8.05 m4。本文方法采用2个薄壁梁单元离散箱梁,分别采用60个Beam-189单元和11340个Shell-63单元建立简支箱梁有限元模型。通过施加刚性扭转荷载并沿梁长均匀布置11道横隔板(t=0.25 m),最大限度地减小箱梁截面畸变对约束扭转正应力和剪应力的影响。本文方法、Beam-189单元模型、Shell-63单元模型及文献[8]的结果列于表2。Beam-189单元模型因忽略闭口周边剪切变形,导致扭转角、约束扭转正应力及翘曲剪应力偏小,而增大截面开口部分剪应力。本文方法和Shell-63单元模型与文献[8]的正应力结果较接近,但剪应力差异较大。例如,Umanskii第二理论计算出位于箱梁闭口周边⑦号点剪应力的相对误差为4%,而开口部分⑧号点剪应力的相对误差则达到70%。原因是悬臂板产生的附加剪流引起悬臂板自由端出现剪应力(由Umanskii第二理论算得⑥号点剪应力为−3.87 kPa),这显然与物理事实不符,势必改变悬臂板剪应力分布。同时悬臂板根部剪力流变化将改变闭口周边剪力流分布。

      表 2  均布扭矩作用下简支箱梁约束扭转位移和应力

      Table 2.  Displacements and stresses for warping torsion of a simply-supported box girder under the action of uniformly distributed torques

      计算模型跨中扭转角/
      (×10−5 rad)
      跨中正应力/kPa梁端剪应力/kPa
      本文一维梁元2.2511.15−9.425.14−3.73−74.15
      Beam-189模型2.1510.86−9.245.01−16.69−51.79
      Shell-63模型2.2210.91−9.365.22−2.82−74.14
      Umanskii第二理论[8]2.2511.41−9.605.32−4.79−71.18
      注:与文献[8]不同,⑦点剪应力包含自由扭转剪应力的贡献。
    • 考虑在单侧腹板顶端(图1截面⑤号点)集中力(P=150 kN)作用下30 m简支箱梁的约束扭转问题。箱梁基本几何及材料参数如下:b1=0.94 m,b2=0.48 m,h=0.42 m,tu=0.04 m,tl=0.07 m,tw=0.06 m,E=3.4×104 N/mm2G=1.45×104 N/mm2。与算例4.1和4.2相比,截面布置更为紧凑(具有更小的宽跨比和高跨比),目的是极大程度地减小板件局部变形(如剪力滞和截面畸变效应)。截面形心和剪心几何位置参数为yo=0.19 m,αy=−0.06 m。截面关键点主扇性坐标为ω6=0.065 m2ω5=−0.055 m2ω4=0.03 m2。与抗扭特性相关的截面常数包括:Iω=1.51×10−4 m6Ipc=0.018 m4Iωc=0.996×10−4 m6Iωω=−0.12×10−4 m6JJB=1.29×10−2 m4。本文方法、Beam-189单元模型、Shell-63单元模型及Umanskii第二理论关于截面扭转角和翘曲正应力计算结果列于表3。再次发现Beam-189模型得到的扭转角偏小,导致截面扭转翘曲峰值正应力相比Shell-63单元结果减小10%,而Umanskii第二理论因在悬臂板形成附加剪流,在本例中减小截面翘曲,导致翘曲正应力偏小。本例箱梁截面翘曲系数$\varepsilon_{\rm T}=l\sqrt{GJ_{\rm{B}}/EI_{\rm{w}}} $=180.66>10,箱梁约束扭转翘曲正应力相比初等梁理论计算结果都控制在10%以内。

      表 3  跨中集中偏载作用下简支箱梁正应力和扭转角

      Table 3.  Normal stresses and angle of twist of a simply-supported box girder subjected to an eccentric transverse load applied at mid-span

      计算模型跨中扭转角/
      (×10−2rad)
      跨中正应力/(×108Pa)/
      弯曲正应力/(×108Pa)
      本文一维梁元2.802−2.69/−2.93−3.14/−2.933.69/3.57
      Beam-189模型2.678−2.70/−2.93−3.12/−2.933.67/3.57
      Shell-63模型2.783−2.82/−2.93−3.32/−2.934.01/3.57
      Umanskii第二理论2.812−2.79/−2.93−3.05/−2.933.63/3.57
    • 第4节所给数值算例仅考虑特定截面下不同计算模型/方法的分析精度,本节将揭示箱梁悬臂板宽及梁高变化对采用不同计算模型/方法得到的箱梁开、闭口截面应力影响规律。以算例4.2的简支箱梁为研究对象,假定箱梁仍然承受均布扭矩(mT=23.5 kN·m/m)的作用。

      首先考虑箱梁悬臂板宽变化对应力的影响。当箱梁顶板宽度保持不变时,为改善箱梁顶板及悬臂板横向受力,悬臂板与顶板宽度之比通常不超过0.5。因此,假定箱梁悬臂板宽从0 m变化到3.5 m。图4给出简支箱梁跨中截面④、⑤和⑥点(如图1所示)正应力随悬臂板宽的变化规律。

      图  4  简支箱梁跨中截面翘曲正应力随悬臂板宽度变化曲线

      Figure 4.  Correlations between the warping normal stresses at mid-span and the width of the cantilever slab of a simply- supported box girder

      图5给出简支箱梁端部截面(x=30)顶板正中⑦号点和悬臂板中部⑧号点剪应力随悬臂板宽的变化曲线。

      图  5  简支箱梁端截面(x=30)剪应力随悬臂板宽度变化曲线

      Figure 5.  Correlations between the shear stresses at the right end and the width of the cantilever slab of a simply-supported box girder

      图4可知,无论箱梁悬臂板宽度如何变化,本文一维梁元和Umanskii第二理论都能准确计算跨中截面翘曲正应力,得到和Shell-63单元模型一致的结果,这是因为它们均能适用于描述箱梁板件翘曲变形。然而,当箱梁悬臂板宽较小时,Beam-189单元模型无法描述箱梁闭口周边对板件翘曲的约束效应,如图4(a)图4(b)所示,闭口周边④和⑤号点正应力的计算误差较大。随着悬臂板宽度增大,截面闭口周边约束效应逐渐减小,Beam-189单元结果趋向于准确解(当悬臂板宽b2=3.5 m时,与Shell-63单元结果完全重合)。此外,Beam-189单元模型能较好地模拟悬臂板翘曲,如图4(c)所示,悬臂板自由端⑥号点正应力结果具有较高精度,说明闭口约束效应对开口部分翘曲的影响可以忽略,这也进一步证明本文针对箱梁截面开、闭口部分采取不同翘曲参数是可行的。

      Umanskii第二理论既能运用平衡关系(简称为U-I类方法),也可由剪应变及弹性本构关系(简称为U-II类方法)计算剪应力。如图5(a)所示采用U-I类方法计算箱梁顶板⑦号点剪应力是依据文献[5, 8, 18]给出的公式。从图5(a)中可以看出,U-I类方法由于忽略箱梁悬臂板翘曲而产生的闭口周边附加剪流,将影响闭口周边翘曲剪流精度,导致⑦号点剪应力随着悬臂板宽的增加偏离准确解。U-II类方法的精度有所改善,原因是并未考虑悬臂板剪流对闭口剪流分布的影响,其实质是计算箱梁板件平均剪应力。一维梁元及Beam-189单元模型计算⑦号点剪应力则考虑闭口周边受力平衡和开、闭口公共节点(⑤号点)剪流连续性条件。Beam-189单元模型因忽略闭口约束效应产生的附加剪流,算得的剪应力均减小30%。一维梁元与Shell-63单元模型结果吻合良好。此外,闭口周边⑦号点剪应力对箱梁悬臂板宽度变化不敏感,这是因为决定截面抗扭刚度的闭口周边扭转与翘曲特征未发生变化。

      由于本文提出的一维梁元及Beam-189单元模型均基于悬臂板中面剪应变为零的假定计算翘曲变形,故悬臂板约束扭转剪应力只能通过平衡关系而求得。如图5(b)所示,一维梁元与Shell-63单元模型结果基本吻合,U-II类方法由于受悬臂板附加剪流影响(算得的悬臂板自由端剪应力不为零),剪应力结果已偏离准确解,只有当悬臂板宽较大时(附加剪流影响减小),才能靠近准确解。采用U-I类方法,剪应力精度有了明显提升,但当悬臂板宽增大时,附加剪流引起的附加翘曲导致剪应力结果偏离准确解。Beam-189单元模型算得的⑧号点剪应力在不同悬臂板宽条件下均大于Shell-63单元模型结果,这是因为其忽略闭口周边附加剪切变形,降低截面自由扭转力矩,从而增大约束扭转力矩,导致约束扭转剪应力结果偏大。

      接下来考虑箱梁高度变化对应力的影响。假定箱梁悬臂板及顶板宽度与算例4.2一致。为了避免闭口周边宽高比过大或过小带来畸变效应,影响Shell-63单元模型约束扭转分析精度‚在板壳元模型内部等间距设置横隔板。图6给出箱梁跨中截面典型节点(④、⑤和⑥点)翘曲正应力随梁高的变化曲线。可以看出各模型/方法计算的截面正应力在不同梁高条件下均吻合良好,说明本文一维梁元和Umanskii第二理论能准确地描述不同梁高条件下箱梁板件翘曲特征。前文已指出,当箱梁悬臂板较宽(b2>2.4 m)时,闭口周边约束效应已减弱,即便梁高增加后(扭转和抗翘刚度增加),由于箱梁腹板难以产生弯翘(板件面内刚度增加),闭口周边约束效应无法显现,Beam-189单元模型结果充分接近准确解。同时,这也解释了梁高增加后,闭口周边翘曲位移及正应力可以忽略(如图6(a)图6(b)所示)。

      图  6  简支箱梁跨中截面翘曲正应力随梁高变化曲线

      Figure 6.  Correlations between the warping normal stresses at mid-span and the height of a simply-supported box girder

      箱梁顶板⑦号点剪应力随梁高变化规律如图7(a)所示。与变化箱梁悬臂板宽度的情形类似,Beam-189单元模型结果因忽略闭口约束效应产生的附加剪流,算得的剪应力偏小。当梁高增加后,由于闭口周边约束效应减弱,剪应力计算精度明显提高。对于U-I类方法,当梁高较小时,附加剪流将引起剪应力计算误差。随着梁高增加,因封闭周边抗翘能力显著增强,剪应力结果趋近准确解。U-II类方法因未考虑悬臂板剪流对闭口周边剪流分布的影响,相比U-I类方法精度有所改善。

      图  7  简支箱梁端截面(x=30)剪应力随梁高变化曲线

      Figure 7.  Correlations between the shear stresses at the right end and the height of a simply-supported box girder

      箱梁悬臂板⑧号点剪应力随梁高变化曲线如图7(b)所示。本文一维梁元与Shell-63单元模型结果在梁高h<2.5 m时吻合良好。随着梁高增加,因闭口周边翘曲迅速减小,悬臂板受力趋向平面应力状态,一维梁元精度有所降低。U-I类方法因受悬臂板附加剪流的影响,在梁高较小时偏离准确解,随着梁高增大,附加剪流的影响不断减弱,又趋向准确解。U-II类方法因无法反映剪应力沿悬臂板宽变化,与准确解存在明显差别。Beam-189单元模型计算的翘曲剪应力显著偏大,因其忽略闭口周边附加剪切变形,从而增大约束扭转力矩。当梁高增大后,由于闭口周边抗扭和抗翘刚度显著增加,附加剪切变形减小,剪应力与准确解的差距随之减小。

    • 本文针对开闭口混合薄壁截面杆件约束扭转分析问题,假定开、闭口截面采用不同翘曲参数。考虑开、闭口部分翘曲连续性条件和截面受力平衡,提出协调翘曲模型。依据势能不变值原理,建立开闭口混合薄壁截面杆件约束扭转分析的一维有限元列式。通过算例分析,论证本文提出的一维梁元具有较高精度。随后开展了不同悬臂板宽和梁高条件下薄壁箱梁截面应力精度分析,对一维梁元、Umanskii第二理论、Beam-189单元模型和Shell-63单元模型性能进行了比较。本文提出的基于协调翘曲场一维梁元具有如下特点:

      (1) 与Beam-189单元模型相比,可以考虑闭口周边对截面翘曲的约束效应,提高箱梁截面正应力与剪应力计算精度。

      (2) 与Umanskii第二理论相比,可以满足开口部分自由端剪应力为零的条件,消除悬臂板和闭口周边“虚拟”附加剪流的影响,提高截面剪应力计算精度。

参考文献 (19)

目录

    /

    返回文章
    返回