CN111753357B - 一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法 - Google Patents

一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法 Download PDF

Info

Publication number
CN111753357B
CN111753357B CN202010480387.0A CN202010480387A CN111753357B CN 111753357 B CN111753357 B CN 111753357B CN 202010480387 A CN202010480387 A CN 202010480387A CN 111753357 B CN111753357 B CN 111753357B
Authority
CN
China
Prior art keywords
section
shear stress
web
cross
corrugated steel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010480387.0A
Other languages
English (en)
Other versions
CN111753357A (zh
Inventor
刘超
黄钰豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN202010480387.0A priority Critical patent/CN111753357B/zh
Publication of CN111753357A publication Critical patent/CN111753357A/zh
Application granted granted Critical
Publication of CN111753357B publication Critical patent/CN111753357B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Civil Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Architecture (AREA)
  • Bridges Or Land Bridges (AREA)
  • Rod-Shaped Construction Members (AREA)

Abstract

本发明涉及一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，包括以下步骤：1)获取变截面波形钢腹板的截面平均剪应力τ；2)采用能量法计算剪力分配系数；3)获取各梁的剪应力，完成波形钢腹板剪应力的分配，并将其作为桥梁设计指标，完成波形钢腹板桥梁设计。与现有技术相比，本发明具有考虑不同腹板之间的剪应力分配效应、普适性强、提高设计安全性等优点。

Description

一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法
技术领域
本发明涉及波形钢腹板桥梁设计领域，尤其是涉及一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法。
背景技术
目前现有的波形腹板剪应力计算方法主要有三种：
(1)将截面总剪力除以波形钢腹板面积得到剪应力，中国和日本的规范均采用这种方法进行计算，该方法忽略了顶底板的承剪比例，认为剪力仅由波形钢腹板承担，对于钢腹板来说偏保守，但对于混凝土顶底板来说偏不安全；
目前国内进行波形钢腹板桥梁的设计时，通常依照《钢-混凝土组合桥梁设计规范：GB50917-2013》进行设计，采用的计算公式：
γ0V≤hwtwfvd (1-1)
式中：V——剪力设计值(N)；hw——钢梁腹板高度(mm)；tw——钢梁腹板厚度(mm)；fvd——钢材的抗剪强度设计值(MPa)；
贾韶丽等研究发现，国外如日本、欧洲、美国规范对于波形腹板的计算为，波形腹板承担的的剪应力等于截面受到的剪力除以腹板面积，差异在于各国取定的容许剪应力大小，中国为155MPa，日本为205MPa，欧洲为160MPa，美国为114MPa。
该方法的问题在于忽略了顶底板的承剪比例，认为剪力仅由波形钢腹板承担，对于钢腹板来说偏保守，但对于混凝土顶底板来说偏不安全。
(2)根据材料力学的方法，利用传统的计算截面静矩的方法计算剪应力，该理论仅仅考虑由剪力产生的剪应力，而忽略了弯矩和轴力产生的剪应力，对于等截面桥梁来说是合适的，但对于变截面桥梁来说准确性不足。
Figure GDA0002591985970000011
式中：Q——剪力设计值(N)；Sy0——计算点处截面静矩(mm3)；I——截面惯性矩(mm4)；b(y)——波形腹板宽度(mm)；
(3)在方法(2)的基础上，同时考虑弯矩和轴力产生的剪应力。该方法计算结果较为准确，但目前尚未推广至单箱多室桥梁。
目前国内进行波形钢腹板桥梁的设计时，忽略了顶底板的承剪比例，认为剪力仅由波形钢腹板承担，对于钢腹板来说偏保守，但对于混凝土顶底板来说偏不安全，而且传统理论对于腹板剪应力计算不准确，误差较大，所以了解波形钢腹板抗剪的机理并且给出剪应力的计算方法是提高设计安全性和准确性的现实要求。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法。
本发明的目的可以通过以下技术方案来实现：
一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，包括以下步骤：
1)获取变截面波形钢腹板的截面平均剪应力τ；
2)采用能量法计算剪力分配系数；
3)获取各梁的剪应力，完成波形钢腹板剪应力的分配，并将其作为桥梁设计指标，完成波形钢腹板桥梁设计。
所述的步骤1)中，变截面波形钢腹板的截面平均剪应力τ的计算式为：
Figure GDA0002591985970000021
其中，τ为截面平均剪应力，I为全截面惯性矩，b(y)为在高度y处截面总宽度，Sy0为坐标0～y0对应的截面静矩，x为取点沿桥轴向长度，α为轴线和水平线间夹角，A为波形钢腹板箱梁的全截面面积，Ay0为坐标0～y0对应的截面面积，N为截面上的轴力，M为截面上的弯矩。
所述的总面积的变化率
Figure GDA0002591985970000022
的表达式为：
Figure GDA0002591985970000023
其中，
Figure GDA0002591985970000024
为上顶板厚度t1沿轴向变化率，
Figure GDA0002591985970000025
为上顶板宽度b1沿轴向变化率，
Figure GDA0002591985970000026
为下顶板厚度t2沿轴向变化率，
Figure GDA0002591985970000027
为下顶板宽度b2沿轴向变化率。
所述的全截面惯性矩的变化率
Figure GDA0002591985970000028
的表达式为：
Figure GDA0002591985970000029
Figure GDA0002591985970000031
Figure GDA0002591985970000032
Figure GDA0002591985970000033
Figure GDA0002591985970000034
Figure GDA0002591985970000035
Figure GDA0002591985970000036
其中，
Figure GDA0002591985970000037
为形心位置yc沿轴向变化率，
Figure GDA0002591985970000038
为箱梁高度h沿轴向变化率，k1～k6分别为变化率系数。
当截面内所取点y0位于上翼缘板内，即0＜y0＜t1时，则有：
Figure GDA0002591985970000039
Figure GDA00025919859700000310
Figure GDA00025919859700000311
Figure GDA00025919859700000312
当截面内所取点y0处于腹板内，即t1≤y≤h-t2时，则有：
Figure GDA00025919859700000313
Figure GDA00025919859700000314
Figure GDA00025919859700000315
Figure GDA00025919859700000316
当截面内所取点y0处于腹板内，即h-t＜y0＜h时，则有：
Figure GDA00025919859700000317
Figure GDA00025919859700000318
Figure GDA00025919859700000319
Figure GDA00025919859700000320
Figure GDA00025919859700000321
所述的步骤2)中，第i号梁的剪力分配系数ki的计算式为：
Figure GDA00025919859700000322
Y(y)＝ax4+bx3+cx2+dx+e
其中，Ii为第i号梁的截面抗弯惯性矩，Y(yi)为第i号梁的横桥向结构挠度曲线，x为取点沿桥轴向长度，a、b、c、d、e为系数。
所述的系数a、b、c、d、e通过拉格朗日乘子法求偏导为0求解得到。
所述的步骤3)中，分配后各梁的剪应力τi的计算式为：
τi＝τ·ki。
与现有技术相比，本发明具有以下优点：
一、本发明解决了变截面波形钢腹板桥梁的剪应力计算问题，同时针对多箱室多腹板截面可以获取各块腹板剪力分配，并计算清楚剪力在腹板和顶底板间的剪力分配，并且同样适用于等截面波形钢腹板桥梁(等截面即为特殊的变截面)，具有普适性。
二、本发明可以准确计算变截面多箱室波形钢腹板剪应力，同时也可以考虑不同腹板之间的剪应力分配效应，采用本方法，能够大大提高相应桥梁的设计计算精度，提高设计安全性。
附图说明
图1为波形钢腹板箱梁截面，其中，图(1a)为波形钢腹板箱梁截面图，图(1b)为截面微元区段。
图2为轴向节段微元。
图3为桥梁纵向挠曲线。
图4为桥梁横向挠曲线。
图5为截面划分图，其中，图(5a)为截面划分成7段纵梁，图(5b)为边梁截面，图(5c)为中梁截面。
图6为桥梁单位横桥向抗弯惯性矩计算图示，其中，图(6a)为取出轴向单位长度节段，图(6b)为等效计算弯矩图示。
图7为波形钢腹板节段参数图。
图8为桥梁横向弯曲应变能计算图示。
图9为本发明的流程框架图。
图10为有限元模型。
图11为箱梁截面参数。
图12为波形钢腹板波形钢腹板节段。
图13为截面选取及腹板定义。
图14为有限元各腹板剪应力大小。
图15为VCM法计算腹板剪应力大小。
图16为能量法计算剪力分配系数。
图17为VCM法与EVCM法腹板剪应力对比。
图18为VCM、EVCM、FEM各腹板剪应力对比，其中，图(18a)为VCM、EVCM、FEM中腹板剪应力对比，图(18b)为VCM、EVCM、FEM次中腹板剪应力对比，图(18c)为VCM、EVCM、FEM次边腹板剪应力对比，图(18d)为VCM、EVCM、FEM边腹板剪应力对比。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
本发明提供一种基于能量法的计算变截面多箱室波形钢腹板箱梁腹板剪应力的获取方法，包括以下步骤：
1)根据桥梁实际数据建立理论模型；
2)根据能量法，计算出桥梁横向应变能与纵向应变能和外力功；
3)将三者使用拉格朗日乘子法计算各参数；
4)将计算出的参数反带回公式中，计算出每块腹板承担的剪应力比例；
5)结合波形腹板剪应力计算方法，计算出每块腹板的剪应力大小。
各步骤的具体说明如下：
一、变截面波形腹板剪应力计算
1、计算基本假定
针对变截面波形钢腹板剪应力计算，计算模型假定如下：(1)不考虑截面剪切变形，截面变形依循平截面假定；(2)钢材与混凝土材料形变均为小变形；(3)假定波形腹板和混凝土顶底板接触紧密，轴向相对滑移量忽略不计；(4)不考虑箱梁加腋造成的的影响；(5)钢和混凝土均为各向同性材料，均处于线弹性受力阶段。(6)只考虑正应力产生的剪应力，忽略畸变和扭转效应。
2、波形钢腹板变截面剪应力计算理论
图(1a)给出了一般波形钢腹板箱梁截面图，箱梁腹板数记为n，上顶板宽度记为b1，厚度记为t1，下顶板宽度记为b2，厚度记为t2，箱梁截面高度记为h，形心距上顶板距离为yc，取任意高度处截面距上顶板距离记为y0，坐标轴原点取在上顶板左上角，向右记为z轴正方向，向下记为y轴正方向。图(1b)给出了距上顶板距离为y0处截取的箱梁截面微元图，沿桥轴向长度为dx，左截面收到的轴向力合力为F，右截面受到的轴向力合力为F+dF，y0处断面存在剪应力流大小记为τ，由于微元长度很小，所以可以认为在dx长度内，剪应力流大小不变，均为τ，列出轴向平衡方程后化简有：
Figure GDA0002591985970000061
式中，b(y)表示在高度y处截面总宽度。
根据剪应力互等定理，公式(2-1)所求剪应力τ与截面剪应力相等，所以由公式(2-1)所求τ大小即可代表截面剪应力。对于图(1b)中，轴向合力F可以看成是轴向FN和由弯矩产生的弯曲轴向力FM的线性叠加，即：
F＝FN+FM(2-2)
其中，对于FN和FM，则有：
Figure GDA0002591985970000062
Figure GDA0002591985970000063
式中，N为截面上轴力的大小，M为截面上弯矩的大小，
Figure GDA0002591985970000064
为轴力在截面上产生的正应力大小，
Figure GDA0002591985970000065
为弯矩在截面上产生的正应力大小，
Figure GDA0002591985970000066
为坐标0～y0截面总面积大小，I为截面惯性矩。
对于公式(2-4)，有：
Figure GDA0002591985970000067
式中，y*为坐标0～y0截面形心坐标，
Figure GDA0002591985970000068
为坐标0～y0截面的静矩。
将公式(2-3)(2-4)(2-5)代回公式(2-2)，可得：
Figure GDA0002591985970000069
截取全截面轴向的一段小微元，示意图如图2所示，图中q为外荷载，α为轴线和水平线间夹角。
对右截面轴线处取距，由力矩平衡：
Figure GDA00025919859700000610
略去含高阶微量的(dx)2项，得：
Figure GDA00025919859700000611
将公式(2-6)代入公式(2-1)，对于公式(2-8)有：
Figure GDA0002591985970000071
通常来说，对于波形钢腹板桥梁而言，一般梁桥不存在轴力或者轴力可忽略；对于斜拉桥而言，轴力很大，但是仅在斜拉索锚固处可能发生轴力的突变，其余位置轴力均匀，不考虑斜拉索锚固截面的情况下，轴力关于轴向变化率更小，所以可以将含
Figure GDA0002591985970000072
项略去，有：
Figure GDA0002591985970000073
注意到若截面为等截面，
Figure GDA0002591985970000074
tanα、
Figure GDA0002591985970000075
项均为0，则公式(2-10)简化为：
Figure GDA0002591985970000076
公式(2-11)即为材料力学中计算截面剪应力的传统公式。
在变截面波形钢腹板箱梁截面中，剪应力由三个部分组成，分别是由弯矩造成的剪应力τM，由轴力造成的剪应力τN，由剪力造成的剪应力τQ，前两者仅存在于变截面剪应力计算中，写成公式(2-12)：
τ＝τM+τQ+τN (2-12)
对于公式(2-12)来说，目前尚未求解的变量有：全截面惯性矩的变化率
Figure GDA0002591985970000077
总面积的变化率
Figure GDA0002591985970000078
截取截面静矩的变化率
Figure GDA0002591985970000079
和截取截面面积的变化率
Figure GDA00025919859700000710
由于波形钢腹板厚度很小，通常远小于板厚，故计算全截面面积和惯性矩时可忽略波形钢腹板面积的影响。由材料力学基本公式可知，全截面面积、惯性矩、截面形心可由公式(2-13)(2-14)(2-15)表示：
A＝b1t1+b2t2 (2-13)
Figure GDA00025919859700000711
Figure GDA00025919859700000712
公式(2-13)对x求偏导，可得公式(2-16)：
Figure GDA0002591985970000081
式中，
Figure GDA0002591985970000082
为上顶板厚度沿轴向变化率，
Figure GDA0002591985970000083
为上顶板宽度沿轴向变化率，
Figure GDA0002591985970000084
为下顶板厚度沿轴向变化率，
Figure GDA0002591985970000085
为下顶板宽度沿轴向变化率。
公式(2-14)对x求偏导，可得公式(2-17)：
Figure GDA0002591985970000086
式中，
Figure GDA0002591985970000087
为形心位置沿轴向变化率，
Figure GDA0002591985970000088
为箱梁高度沿轴向变化率。各项变化率系数k1～k6分别为：
Figure GDA0002591985970000089
Figure GDA00025919859700000810
Figure GDA00025919859700000811
Figure GDA00025919859700000812
Figure GDA00025919859700000813
Figure GDA00025919859700000814
当所取的y0点处于不同位置时，
Figure GDA00025919859700000815
和
Figure GDA00025919859700000816
的计算方法不同，故分类考虑：
(a)当点y0处于上翼缘板内(0＜y0＜t1)有：
Figure GDA00025919859700000817
Figure GDA00025919859700000818
Figure GDA00025919859700000819
Figure GDA00025919859700000820
(b)当点y0处于腹板内(t1≤y≤h-t2)，则有：
Figure GDA00025919859700000821
Figure GDA00025919859700000822
Figure GDA00025919859700000823
Figure GDA00025919859700000824
(c)当点y0处于下翼缘板内(h-t＜y0＜h)
Figure GDA00025919859700000825
Figure GDA00025919859700000826
Figure GDA0002591985970000091
Figure GDA0002591985970000092
将公式(2-13)～(2-21)代入公式(2-10)即可求得截面上任意位置的剪应力大小。
二、基于能量法的变截面波形钢腹板剪力分配计算
1.挠度曲线与挠度面
以自由端受均布线荷载的悬臂梁为例，结构宽为B，悬臂长为L，f(x)为P＝1作用下悬臂梁挠度曲线，其抗弯惯距为整体截面，计算图式及理论分析如图3所示。
参考结构力学，对于悬臂梁，在自由端给定单位荷载，利用力法原理，计算桥梁纵桥向挠度曲线为：
Figure GDA0002591985970000093
式中，E为混凝土弹性模量，I为截面抗弯惯性矩，l为桥梁全长，x为取点距自由端距离。
在每一个横桥向，桥面发生横向弯曲，假设横向弯曲的曲面是四次曲线，则横桥向结构挠度曲线Y(y)为：
Y(y)＝ax4+bx3+cx2+dx+e (2-23)
式中，a、b、c、d、e为系数，图4中a1-a4分别为对应腹板的挠度
在实际结构悬臂端处横桥方向(y方向)作用均布线荷载q情况下结构的挠度面为W(x,y)；
假定挠度面W(x,y)近似可表达成：
W(x,y)＝f(x)Y(y) (2-24)
2.顺桥向弯曲应变能
设全桥抗弯惯性矩为I(x方向)，截面划分如图5所示：
则第i号截面的抗弯惯距为I(x方向)，第i号梁顺桥向弯曲应变能为：
Figure GDA0002591985970000094
若纵梁有n根，则全桥顺桥向弯曲应变能Ux计算公式如下：
Figure GDA0002591985970000095
3.横桥向弯曲应变能
横桥向结构可以看成很多空腹梁组成，而空腹梁是剪切变形为主，不能把上下弦截面看成整体刚性截面，因此必须近似换算。考虑将一室轴向单位长度的混凝土顶底板和波形钢腹板取出，如图(6a)所示，等效计算弯矩图示如图(6b)。
假定一侧刚接，另一侧施加单位力P，可利用图乘法求出施加单位力一侧腹板下挠ΔQ：
Figure GDA0002591985970000101
式中，Es为钢板的弹性模量，E为混凝土的弹性模量，
Figure GDA0002591985970000102
为单位力作用下的弯矩图，MP为外力P作用下的弯矩图。Ic为混凝土顶底板的抗弯惯性矩，Is为波形钢腹板的抗弯惯性矩，乘以1/2的原因是每块腹板的刚度被两侧箱室平分。
波形钢腹板的等效面外惯性矩Is计算公式为：
Figure GDA0002591985970000103
式中，Es为钢板的弹性模量，E为混凝土的弹性模量，t为波形钢腹板厚度，α为腹板斜板段倾角，e为波高，Ls为直板段长度，Lz为斜板段投影长度，q为波长，如图7所示。
有了ΔQ后，可以用等代法，计算出横桥向抗弯等效刚度Iy。
Figure GDA0002591985970000104
Figure GDA0002591985970000105
接下来可以计算横桥向弯曲应变能Uy，其计算图式如图8。
取dx宽的一条横梁，假定该横梁的单位长度抗弯惯距Iy已知，则该横梁的弯曲应变能：
Figure GDA0002591985970000106
所以全桥横向弯曲应变能为：
Figure GDA0002591985970000107
4.外力功计算
在实际结构自由端受均布线荷载的情况下，外力功如下式计算：
Figure GDA0002591985970000111
5.参数计算
方程中变量为假设的横桥向挠度曲线的系数a、b、c、d、e共5个，采用拉格朗日乘子法，整个系统对系数求偏导为0即可解出各个系数的值。
π＝Ux+Uy-外荷载功
Figure GDA0002591985970000112
式中，δ为Y(y)中系数，由上述方程可解出。
6.剪力分配系数
根据(2-21)式，若受到剪力大小为Q，则桥梁挠曲线方程为：
Figure GDA0002591985970000113
式中，Δ为受到剪力Q时挠曲线方程。
对于每段梁，则可写出如下表达式：
Figure GDA0002591985970000114
式中，Qi为第i号梁分配到的剪力。
将(2-35)除以(2-34)，可得到：
Figure GDA0002591985970000115
所以第i号梁的剪力分配系数为：
Figure GDA0002591985970000116
三、基于能量法的变截面波形钢腹板剪应力计算
综上所示，本发明提出基于能量法的变截面波形钢腹板剪应力计算方法。
结合式(2-37)和式(2-10)，对于第i号梁，其剪应力计算公式为：
Figure GDA0002591985970000117
式中，n表示单箱多室截面腹板数量，其余参数含义见式(2-10)。
实施例：
为了更进一步突出弯矩对腹板剪应力的影响，结合工程实际，选取了梁长50m，截面为变截面单箱六室波形钢腹板截面，结构类型为悬臂梁，软件选择为ABAQUS，有限元模型示意图如图10所示。自由端截面为小截面，固定端为大截面，中间截面参数线性过渡，箱梁截面参数如图11所示。在悬臂梁自由端施加竖向集中力1000kN，波形钢腹板示意图如图12所示。材料的选择方面，本模型只采用两种材料：钢和混凝土，均假设处于线弹性状态、各向同性材料。同时，将混凝土顶底板和波形钢腹板接触面自由度耦合，模拟无轴向滑移状态。网格划分方面，混凝土顶底板采用空间3D网格，波形钢腹板采用平面2D网格，共计划分276201个节点、207023个单元。材料属性设置见表1所示。
表5.1悬臂梁各参数设置
材料 密度/kg·m<sup>-3</sup> 弹性模量/GPa 泊松比 单元类型
混凝土 2600 35.5 0.2 C3D8(3D)
钢 7800 206.0 0.3 S4R(2D)
结果分析
(1)剪应力分析
截面选取及腹板定义如图13所示，为避免端部复杂应力影响，截取40米桥长，每隔一米取截面计算一次剪应力大小。对于该悬臂梁，对本文第二章提出的变截面计算理论VCM式(2-10)、基于能量法的变截面剪应力计算理论EVCM式(2-10)和有限元(FEM)方法的计算结果进行对比分析。
有限元模型结果如图14所示，VCM计算结果如图15所示，跨径范围内，能量法计算剪力分配系数如图16所示，VCM法与EVCM法腹板剪应力对比如图17所示。对于图17，主要可以得出三个结论：(1)EVCM法计算各腹板剪应力大小实质相当于VCM法的平移，变化趋势与VCM法相同；(2)同一位置处，EVCM法中腹板剪应力最大，次中腹板其次，次边腹板第三，边腹板最小，VCM法计算的腹板剪应力大小介于EVCM法中次边腹板和边腹板之间；(3)中腹板、次中腹板、次边腹板剪应力大小较为接近，边腹板剪应力明显较小。
VCM、EVCM、FEM各腹板剪应力对比结果如图18所示，图(18a)为中腹板结果对比，可以看出EVCM剪应力较大，EVCM与VCM结果与有限元变化趋势相同，但VCM且结果略小，可看出EVCM结果与有限元较为接近(除了在距离自由端5-10米处误差较大)。图(18b)为次中腹板结果对比，可以看出EVCM剪应力较大，EVCM与VCM结果与有限元变化趋势相同，但VCM且结果略小，可看出EVCM结果与有限元较为接近(除了在距离自由端5-10米处误差较大)。图(18c)为次边腹板结果对比，可以看出EVCM剪应力较大，EVCM与VCM结果与有限元变化趋势相同，但VCM且结果略小，可看出EVCM结果和VCM结果都与有限元较为接近，EVCM在距离自由端5-10米处误差最小。图(18d)为边腹板结果对比，可以看出EVCM剪应力较小，EVCM与VCM结果与有限元变化趋势接近，但VCM且结果略大，EVCM在距离自由端5-35米处结果小于有限元结果，在35-45米处结果大于有限元。总体来说，EVCM结果准确性高于VCM。

Claims (9)

1.一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，其特征在于，包括以下步骤：
1)获取变截面波形钢腹板的截面平均剪应力τ；
2)采用能量法计算剪力分配系数，第i号梁的剪力分配系数ki的计算式为：
Figure FDA0003542782030000011
Y(y)＝ax4+bx3+cx2+dx+e
其中，Ii为第i号梁的截面抗弯惯性矩，Y(yi)为第i号梁的横桥向结构挠度曲线，x为取点沿桥轴向长度，a、b、c、d、e为系数；
3)获取各梁的剪应力，完成波形钢腹板剪应力的分配，并将其作为桥梁设计指标。
2.根据权利要求1所述的一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，其特征在于，所述的步骤1)中，变截面波形钢腹板的截面平均剪应力τ的计算式为：
Figure FDA0003542782030000012
其中，τ为截面平均剪应力，I为全截面惯性矩，b(y)为在高度y处截面总宽度，
Figure FDA0003542782030000013
为坐标0～y0对应的截面静矩，x为取点沿桥轴向长度，α为轴线和水平线间夹角，A为波形钢腹板箱梁的全截面面积，
Figure FDA0003542782030000014
为坐标0～y0对应的截面面积，N为截面上的轴力，M为截面上的弯矩。
3.根据权利要求2所述的一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，其特征在于，所述的总面积的变化率
Figure FDA0003542782030000015
的表达式为：
Figure FDA0003542782030000016
其中，
Figure FDA0003542782030000017
为上顶板厚度t1沿轴向变化率，
Figure FDA0003542782030000018
为上顶板宽度b1沿轴向变化率，
Figure FDA0003542782030000019
为下顶板厚度t2沿轴向变化率，
Figure FDA00035427820300000110
为下顶板宽度b2沿轴向变化率。
4.根据权利要求3所述的一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，其特征在于，所述的全截面惯性矩的变化率
Figure FDA00035427820300000111
的表达式为：
Figure FDA00035427820300000112
Figure FDA0003542782030000021
Figure FDA0003542782030000022
Figure FDA0003542782030000023
Figure FDA0003542782030000024
Figure FDA0003542782030000025
Figure FDA0003542782030000026
其中，
Figure FDA0003542782030000027
为形心位置yc沿轴向变化率，
Figure FDA0003542782030000028
为箱梁高度h沿轴向变化率，k1～k6分别为变化率系数。
5.根据权利要求4所述的一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，其特征在于，当截面内所取点y0位于上翼缘板内，即0＜y0＜t1时，则有：
Figure FDA0003542782030000029
Figure FDA00035427820300000210
Figure FDA00035427820300000211
Figure FDA00035427820300000212
6.根据权利要求4所述的一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，其特征在于，当截面内所取点y0处于腹板内，即t1≤y≤h-t2时，则有：
Figure FDA00035427820300000213
Figure FDA00035427820300000214
Figure FDA00035427820300000215
Figure FDA00035427820300000216
7.根据权利要求4所述的一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，其特征在于，当截面内所取点y0处于腹板内，即h-t＜y0＜h时，则有：
Figure FDA00035427820300000217
Figure FDA00035427820300000218
Figure FDA00035427820300000219
Figure FDA00035427820300000220
8.根据权利要求1所述的一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，其特征在于，所述的系数a、b、c、d、e通过拉格朗日乘子法求偏导为0求解得到。
9.根据权利要求1所述的一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法，其特征在于，所述的步骤3)中，分配后各梁的剪应力τi的计算式为：
τi＝τ·ki。
CN202010480387.0A 2020-05-30 2020-05-30 一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法 Active CN111753357B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010480387.0A CN111753357B (zh) 2020-05-30 2020-05-30 一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010480387.0A CN111753357B (zh) 2020-05-30 2020-05-30 一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法

Publications (2)

Publication Number Publication Date
CN111753357A CN111753357A (zh) 2020-10-09
CN111753357B true CN111753357B (zh) 2022-08-09

Family

ID=72674451

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010480387.0A Active CN111753357B (zh) 2020-05-30 2020-05-30 一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法

Country Status (1)

Country Link
CN (1) CN111753357B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112784333B (zh) * 2020-12-22 2023-12-05 中冶南方城市建设工程技术有限公司 一种现浇混凝土箱梁横梁的设计方法
CN112699448B (zh) * 2021-01-13 2023-07-18 中铁大桥勘测设计院集团有限公司 一种全焊钢桁梁整体节点撕破计算方法
CN117094063B (zh) * 2023-10-17 2024-01-26 山东大学 波形钢腹板模拟方法及系统、组合梁桥建模方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000017613A (ja) * 1998-06-29 2000-01-18 Ps Corp 波形鋼板ウエブ桁の接続方法
CN103743630A (zh) * 2013-12-25 2014-04-23 广西科技大学 波形钢腹板组合梁的剪切性能测试方法
CN104695339A (zh) * 2015-02-09 2015-06-10 河南省交通科学技术研究院有限公司 采用波形钢腹板进行混凝土箱梁抗剪加固的施工方法
CN107247856A (zh) * 2017-08-01 2017-10-13 西安电子科技大学 一种单滚柱包络环面蜗杆副时变啮合刚度解析方法
CN108519433A (zh) * 2018-04-10 2018-09-11 长沙理工大学 一种波形钢腹板之间横向连接结构的健康监测装置及方法
CN109858071A (zh) * 2018-12-06 2019-06-07 山西大学 一种考虑剪力滞后作用的薄壁箱梁结构动力特性分析方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000017613A (ja) * 1998-06-29 2000-01-18 Ps Corp 波形鋼板ウエブ桁の接続方法
CN103743630A (zh) * 2013-12-25 2014-04-23 广西科技大学 波形钢腹板组合梁的剪切性能测试方法
CN104695339A (zh) * 2015-02-09 2015-06-10 河南省交通科学技术研究院有限公司 采用波形钢腹板进行混凝土箱梁抗剪加固的施工方法
CN107247856A (zh) * 2017-08-01 2017-10-13 西安电子科技大学 一种单滚柱包络环面蜗杆副时变啮合刚度解析方法
CN108519433A (zh) * 2018-04-10 2018-09-11 长沙理工大学 一种波形钢腹板之间横向连接结构的健康监测装置及方法
CN109858071A (zh) * 2018-12-06 2019-06-07 山西大学 一种考虑剪力滞后作用的薄壁箱梁结构动力特性分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
变截面波形钢腹板箱梁剪应力计算理论;刘超 等;《同济大学学报（自然科学版）》;20190430;第47卷(第4期);第475-481页 *

Also Published As

Publication number Publication date
CN111753357A (zh) 2020-10-09

Similar Documents

Publication Publication Date Title
CN111753357B (zh) 一种变截面多箱室波形钢腹板箱梁腹板剪应力的分配方法
CN110032829B (zh) 钢-混凝土组合梁的应力计算方法
CN109241604B (zh) 考虑滑移作用的钢-混组合梁桥横向分布系数有限元计算方法
JP6699639B2 (ja) 床構造に用いる鉄骨梁の設計方法、床構造
CN108487542A (zh) 钢管混凝土束组合构件稳定承载力计算方法
CN106894328B (zh) 一种π形结合梁剪力滞的处理方法
CN109446715B (zh) 一种大跨度桥梁地震响应时程分析方法
CN114638046B (zh) 一种铁路桥墩数字孪生变截面仿真计算方法
JP2018131883A (ja) 床構造
CN114218658A (zh) 一种适用于锚索框架结构的内力变形解析计算方法
Bozdogan et al. An approximate method for lateral stability analysis of wall-frame buildings including shear deformations of walls
CN106940736A (zh) 一种考虑侧向载荷作用的纵骨多跨失稳的极限载荷确定方法
CN114218655A (zh) 一种变截面波形钢腹板组合梁桥的剪应力实用计算方法
Wang et al. Non-linear buffeting response analysis of long-span suspension bridges with central buckle
Han et al. Analysis of Steel Bridge by means of Specially Orthotropic Plate Theory
CN113656861A (zh) 一种由铁路桥梁变形导致的无砟轨道空间映射变形的计算方法
CN105046027A (zh) 一种多肋式t形梁桥断面的优化设计方法
JP7348509B2 (ja) 連続梁の評価方法及び連続梁の評価プログラム
CN111723428A (zh) 基于位移影响线分配的连续刚构桥成桥预拱度设置方法
CN111611642A (zh) 一种计算外包波纹侧板混凝土组合梁的抗弯刚度方法
CN115630458A (zh) 一种基于弹性地基梁理论的现浇连续梁的方法及其应用
CN111241727B (zh) 一种利用单跨有限元模型进行豪华邮轮极限强度计算的方法
Cheng et al. Comparison of numerical techniques for 3D flutter analysis of cable-stayed bridges
CN109235274A (zh) 一种新型顶推施工临时垫梁系统及其施工方法
CN107700336A (zh) 一种混凝土斜拉桥主梁施工阶段剪力滞的确定方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant