板桁组合结构分析中的板梁单元
陈玉骥, 叶梅新
(中南大学 土木建筑学院,湖南 长沙,410075)
摘要: 为了能合理和方便地分析板桁组合结构的受力性能,将纵、横梁视为混凝土板的加劲梁,建立考虑板梁共同作用和相对滑移板梁单元、无共同作用和过渡区的板梁单元。其中,将矩形板视为单元的基本部分,梁视为附属部分,即以矩形板4个结点的结点位移以及每条板梁结合面(有钉区)中线上的3个相对滑移自由度作为板梁单元的基本自由度,以构造板和梁的位移模式,由此得到板梁单元的单元刚度矩阵。
关键词: 板桁组合结构; 板梁单元; 位移模式; 刚度矩阵
中图分类号:U448.38 文献标识码:A 文章编号: 1672-7207(2005)04-0715-06
Plate-beam element in calculating plate-truss composite structure
CHEN Yu-ji, YE Mei-xin
(School of Civil and Architectural Engineering, Central South University, Changsha 410075, China)
Abstract: For analyzing the properties of plate-truss composite structures expediently, the plate-beam elements were set up in which the action and slip between the plates and beams were considered. Composite elements in which plates have no action to beams and excessive elements were also established. In these elements, plates were the basic parts, and beams were the adjunctive part. That is to say, the displacement pattern of plates and beams was constructed with the basic free degrees of the displacement of 4 nodes on rectangular plate and the relative slips of 3 nodes at combining surface between beams and plates. Thereby, stiffness matrixes of plate-beam elements were obtained.
Key words: plate-truss composite structure; plate-beam element; displacement pattern; stiffness matrix
板桁组合结构是桥面系纵、横梁与混凝土板结合的一种桥梁结构形式,它具有建筑高度较低,行车时燥声和震动较小,刚度较大等优点,是高速铁路桥梁中比较理想的结构形式之一。对于这种结构,可以用工程软件进行有限元分析[1-8],主要采用桁单元、梁单元、平板壳单元和块体单元。但采用这些单元时,板梁结合面除结点外钢梁与混凝土板的位移很难协调,且难以考虑钢梁与混凝土板在结合面上的相对滑移[9-14]。可见,采用常规单元离散桥面系不甚理想。为了能较合理和方便地分析板桁组合结构的受力性能,作者将纵、横梁视为混凝土板的加劲梁,建立了包含混凝土板与纵、横梁的组合单元(以下简称为板梁单元)。该单元以矩形板4个结点的结点位移和每条板梁结合面中线上3个相对滑移自由度作为板梁单元的独立自由度来构造板与梁的位移模式。在无钉区,若不考虑钢梁与混凝土板之间的纵向摩擦,则钢梁与混凝土板在结合面不存在[CM(22] 纵向剪力,故板与梁仅竖向位移和法向转角相互协调,而纵向、横向位移则相互独立,由此即可构造这种无钉区单元的位移模式。此外,为了使有钉区单元和无钉区单元能有效过渡,还必须建立过渡区板梁单元的单元刚度矩阵。根据板桁组合结构桥面系纵梁通过剪力连接件与混凝土板完全结合,横梁与混凝土板部分结合的特点,有钉区板梁单元应有2种形式:一种为双梁板梁单元,另一种为单梁板梁单元。在此,作者分别讨论有钉区、无钉区和过渡区板梁单元的位移模式及其单元刚度矩阵。
1 板梁单元的单元分析
1.1 考虑板梁共同作用的板梁单元
图1所示为双梁板梁单元,结点1~4分别为板中面的4个角点,s,t,m,n和o分别为板梁结合面的中线的端点和中点。因混凝土桥面板除发生弯曲变形外,还产生面内变形,故该板可用平板壳单元进行模拟,其应力状态由平面应力问题和薄板小挠度弯曲问题迭加而成,而梁的位移模式可根据板与梁的变形协调关系建立。
图 1 双梁板梁单元
Fig. 1 Plate-beam element with two beams
1.1.1 板的位移模式及其对单元刚度矩阵的贡献
令 ,将坐标进行无量纲处理。为简化分析,引入如下假定:
a. 板与梁不产生竖向分离;
b. 板与梁在结合面的剪力连接件为连续分布。
如前所述,板采用平板壳单元,而常规4个结点矩形平板壳单元是以结点挠度及其一阶导数(转角)作为单元的自由度,它是一种非协调元,即单元边界上的挠度及其沿边界方向的偏导数是连续的,但在单元交界线之间,挠度沿垂直于边界方向的偏导数是不连续的,即单元边界上的法向导数是不连续的。为了使相邻板梁单元的梁在相连接的点(端点)位移协调,应将矩形平板壳单元位移(挠度)模式修正为:
其中: 为修正项:
而Nx(ζ,η),Nxi(ζ,η)和Nyi(ζ,η)与常规的4结点矩形板弯曲单元形函数相同;kx和ky为常数,当单元中含有双梁(纵、横梁)时,kx=ky=1;当板中只有与x轴(或y轴)平行的单梁时,kx=1,ky=0(或ky=1,kx=0)。
板中面的面内位移为:
式中:Mi (ζ,η)为平面应力问题4个结点矩形单元的形函数;[AKu-D]i(i为s,o,t),[AKv-D]i(i为m,o,n)为板与梁在结合面中线上结点i的相对滑移;
分别为考虑板与梁之间相对滑移不引起纵向力而引入的系数;ELx,ELy,ALx和ALy分别为与x轴平行的梁和与y轴平行的梁的弹性模量和截面积;EP,APx和APy分别为板的弹性模量和板垂直于x轴和y轴截面的面积;而Lxi(ζ,η)和Lyi(ζ,η)为板与梁在结合面相对滑移对应的形函数,因结合面处的正应变沿梁轴为线性函数,与其相对应的滑移应变亦以采用线性函数为宜,故组合单元中每一个板梁结合面的中线设置3个结点(2个端点和正中点),以这3点的滑移构造滑移对单元位移的影响函数,有:
不难看出,以上形函数 , 在单元边界上的值均为零,故上述位移模式不会改变常规的4个结点矩形板单元在结点处的位移连续条件。
由几何方程可得板的应变为:
式中:B=[Bu Bw Bd],为应变矩阵;δ=[δu δw δd]T,为板梁单元的结点位移阵;δu=[u1 v1 u2 v2 u3 v3 u4 v4]T,为平面应力问题的结点位移列阵;δw=[w1 θx1 θy1 w2 θx2 θy2 w3 θx3 θy3 w4 θx4 θy4]T,为薄板小挠度弯曲问题的结点位移列阵; ,为钢与混凝土相对滑移列阵。
根据最小势能原理可得板对单元刚度矩阵的贡献:
1.1.2 梁的位移模式及其对单元刚度矩阵的贡献
首先分析与x轴平行的梁。根据变形协调条件,可将梁的位移用板的结点位移和相对滑移表示为:
式中: ,为与系数fx对应的系数。
梁的应变为:
故梁对单元刚度矩阵的贡献为:
若梁与y轴平行,推导过程及结果与上类似。
1.1.3 剪力连接件对单元刚度矩阵的贡献
根据板与梁的相对滑移,可写出剪力连接件的变形能为:
式中:kdx和kdy分别为剪力连接件单位长度的剪切刚度;
故对剪力连接件对板梁单元的单元刚度矩阵的贡献为:
若板梁单元中只有1条 (平行于x轴或平行于y轴)梁(见图2),其单元刚度矩阵可按以上双梁板梁单元的结果推出。即若单元仅有与x轴平行的梁,则其单元刚度矩阵可将以上单元刚度矩阵中的与滑移δdy对应行列的全部元素置零,并令kx=1,ky=0即可;若仅有与y轴平行的梁,则其单元刚度矩阵可将以上单元刚度矩阵中与滑移δdx对应行列的全部元素置零,并令kx=0,ky=1即可。
图 2 单梁板梁单元
Fig. 2 Plate-beam elements with one beam
1.2 板梁无共同作用的板梁单元
在无钉区,梁与板无共同作用,但仍假设两者不能分离,即板梁的竖向挠度相同,但梁两端点的面外位移为独立位移。这样,板的位移(挠度)模式与上述考虑板梁共同作用的单梁板梁组合单元的相同,中面面内位移与传统平面应力问题的位移模式相同。下面分析梁对单元刚度矩阵的贡献。
考虑与x轴平行的梁(见图2),其挠度与板η=0处的挠度相同,为:
梁轴线的面外(y方向)位移可根据梁轴线两端点的位移及转角用Hermit插值写出:
梁轴线的轴向位移为:
式中: ;h0为板中面至梁轴线的距离。
梁的应变为:
式中:
若梁平行于y轴,同理可按上述相同方法推导。
1.3 过渡区板梁单元
为了使上述有钉区(板梁有共同作用)和无钉区(板梁无共同作用)的板梁单元能有效过渡,还需建立过渡区板梁单元。图2所示也可视为一个过渡区板梁组合单元,其中s端为有钉端,该端纵梁与板存在共同作用; t端为无钉端,该端梁与板仅竖向挠度及腹板面内转角相同,但纵向、横向位移相互独立。故板挠度与前述单梁(平行于x轴)板梁组合单元的w(ζ,η) 表达式相同,而梁t端的水平位移为独立位移,s端的位移可由板的位移按照变形协调条件用板的结点位移和滑移表示。为了使位移模式在单元内不产生突变,设板与梁的共同作用由s端到t端逐渐减弱至无共同作用,故板的位移模式为:
式中: ,为反映板与梁共同作用程度的函数。
梁的挠度与板中线(η=0)处的相同,梁轴线沿y方向的位移和轴向位移为:
梁的应变为:
式中:δ为过渡区板梁单元的结点位移列阵 ,为梁t 端的结点位移列阵。
2 实验验证
取文献[1]中的4个节间简支下承式板桁组合桥实验模型(见图3),用上述公式进行有限元分析。其中,纵梁为I20a工字钢,横梁为I25b工字钢,主桁上、下弦杆和斜杆16,26,59和49为H150×150宽翼H型钢,中间斜杆27,37,38和48为I14工字钢,上、下平联和制动撑为L63×63×5等边角钢;混凝土等级为C50,板厚10 cm;荷载P=100 kN,钢的弹性模量为ES=2.1×105 MPa,混凝土的弹性模量为EC=0.355×105 MPa;纵梁上布满剪力钉,横梁除靠主桁两端各0.335 m范围外,其他位置均布满剪力钉。单元划分情况为:主桁各杆采用梁单元,上、下平联和制动撑采用桁元,在桥面系纵横梁处采用作者所提出的板梁单元,其他地方的混凝土板采用常规矩形板壳元。主要计算结果和试验结果见表1 和图4~7。
图 3 下承式钢杵结合梁模型
Fig. 3 Model of through truss composite beam
图 4 主桁竖向位移
Fig. 4 Vertical displacements of main truss
图 5 横梁2下翼缘内侧端点正应力
Fig. 5 Stress at bottom fringe of traverse beam 2
图 6 纵梁下翼缘内侧端点正应力
Fig. 6 Stress at bottom fringe of longitudinal beam
图 7 节间2-3跨中混凝土板截面正应力
Fig. 7 Stress of concrete slab at middle of space 2-3
表 1 主桁各杆的轴应力
Table 1 Axial forces in each beam of main truss
3 结 论
a. 根据板桁组合结构中桥面系的变形特征,给出了考虑板梁共同作用(有钉区)、无共同作用(无钉区)和过渡区板梁单元的位移模式。
b. 所提出的板梁单元既可用于计算板桁组合结构,也可用于计算板梁组合结构和带肋加劲板,是一种适用于具有板梁组合性质的结构分析的较好单元。
参考文献:
[1]叶梅新,陈玉骥,张晔芝,等.下承式桁梁结合梁模型试验和有限元分析[R].长沙:中南大学土木建筑学院,2003.
YE Mei-xin, CHEN Yu-ji, ZHANG Ye-zhi, et al. Test of concrete-steel through truss composite girder models and analysis of finite element method[R]. Changsha: School of Civil and Architectural Engineering, Central South University, 2003.
[2]周胜利,林亚超.高速铁路钢桁梁桥桥面结构设计及减小噪音的结构措施[J].国外桥梁, 1996(3): 22-32.
ZHOU Sheng-li, LIN Ya-chao. High-speed railway steel truss bridge deck structural design and structural measure of reducing the noise[J]. Foreign Bridges, 1996(3): 22-32.
[3]周胜利,林亚超.日本北陆新干线犀川桥[J].国外桥梁,1996(3): 1-11.
ZHOU Sheng-li, LIN Ya-chao. Xichuan Bridge on the northland new railway of Japan[J]. Foreign Bridges, 1996(3): 1-11.
[4]彭月桑. 下弦杆和钢桥面结合的低高度桥面系铁路桁梁桥[J]. 国外桥梁, 1995(2): 105-107.
PENG Yue-sang. Railway steel truss bridge with low height floor system of chord member combining with steel deck[J]. Foreign Bridges, 1995(2): 105-107.
[5]叶梅新,江锋. 芜湖桥板桁组合结构的研究[J]. 铁道学报, 2001, 23(5):65-69.
YE Mei-xin, JIANG Feng. Study on the concrete slab-steel truss composite structures for Wuhu River Bridge[J]. Journal of the China Railway Society, 2001, 23(5):65-69.
[6]谭莹,田启贤.板桁组合结构的受力特性及空间分析方法[J] .铁道建筑, 2001(8):2-5.
TAN Ying, TIAN Qi-xian. Features of force bearing and method of spatial analysis of plate and truss composite structure[J]. Railway Construction, 2001(8): 2-5.
[7]何畏,李乔.板桁组合结构体系受力特性及计算方法研究[J] .中国铁道科学, 2001, 22(5): 65-72.
HE Wei, LI Qiao. Study on features of force bearing of plate and truss composite structure system and the calculation method[J]. China Railway Science, 2001, 22(5): 65-72.
[8]周惟德,陈辉求. 钢—混凝土组合桁架有限元计算模型[J]. 工业建筑, 1996, 26(7): 14-19.
ZHOU Wei-de, CHEN Hui-qiu. The finite element calculating model of composite truss[J]. Industrial Architecture, 1996,26(7):14-19.
[9]张凤鸣,熊健民,周金枝,等. 板桁组合结构分析的超级有限元法[J].湖北工学院学报, 2000, 15(1): 21-23.
ZHANG Feng-ming, XIONG Jian-min, ZHOU Jin-zhi, et al. Static analysis of plate-truss composite structure by super finite element[J]. Journal of Hubei Polytechnic University, 2000, 15(1): 21-23.
[10]朱伯芳.有限单元法原理及应用(第2版)[M].北京:中国水利水电出版社,1998.
ZHU Bo-fang. The theory and application of finite element method [M]. Beijing: China Water and Electricity Press, 1998.
[11]刘世恩. 加劲板的计算方法研究[J]. 桥梁建设,1998, 15(1): 21-23.
LIU Shi-en. Research on the computer method of stiffened plates[J]. Bridge Construction, 1998,15(1):21-23.
[12]赵鸿铁. 钢与混凝土组合结构[M].北京:科学出版社,2001.
ZHAO Hong-tie. Steel and concrete composite structure[M]. Beijing: Science Press, 2001.
[13]陈玉骥,叶梅新,伍洪泽.钢桁梁桥纵梁和横梁的受力分析[J]. 中国有色金属学报, 1996, 6(2):151-155.
CHEN Yu-ji, YE Mei-xin, WU Hong-ze.The stress analysis of longitudinal and transverse beams of railway steel truss bridges[J]. The Chinese Journal of Nonferrous Metals, 1996, 6(2):151-155.
[14]陈玉骥,叶梅新. 钢—混凝土连续结合梁的温度效应[J].中南大学学报(自然科学版), 2004, 35(1): 142-146.
CHEN Yu-ji, YE Mei-xin. Temperature response of continuous composite beam[J]. Journal of Central South University(Science and Technology), 2004, 35(1): 142-146.
收稿日期:2004-10-15
基金项目:铁道部科技研究开发计划项目(99G33)
作者简介:陈玉骥(1962-),男,海南文昌人,教授,博士,从事钢—混凝土组合结构和弹塑性力学及其应用研究
论文联系人: 陈玉骥,教授,博士;电话:0731-2655757(H),13874999878(手机);E-mail:Chenyj@mail.csu.edu.cn