变截面波形钢腹板悬臂梁抗剪性能研究

王方旭 1，刘世忠 1?，杨霞林 1，秦翱翱 2，刘欣益 1; WANG Fangxu1，LIU Shizhong1?，YANG Xialin1，QIN Aoao2，LIU Xinyi1

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

变截面波形钢腹板悬臂梁抗剪性能研究 PDF

- ORCID：
王方旭 ¹
- ORCID：
刘世忠 ¹
✉
- ORCID：
杨霞林 ¹
- ORCID：
秦翱翱 ²
- ORCID：
刘欣益 ¹

1. 兰州交通大学土木工程学院，甘肃兰州 730070； 2. 枣庄学院城市与建筑工程学院，山东枣庄 277160

中图分类号： U443.3

最近更新：2025-01-21

DOI： 10.16339/j.cnki.hdxbzkb.2025010

摘要

为了在变截面波形钢腹板悬臂梁的剪应力计算时考虑翼板、腹板以及组合箱梁截面转角差异，首先从翼板、腹板以及组合箱梁截面转角出发建立位移函数，采用能量变分法将翼板与腹板承剪进行剥离，利用弯矩等效分离顶、底板承剪；其次结合变截面波形钢腹板组合箱梁的梁段分析单元刚度矩阵及节点荷载列阵建立了剪应力求解程序；最后分析了悬臂梁在不同荷载工况下顶、底板及腹板剪应力与承剪比.结果表明：相较于有限元计算结果，考虑转角差异后顶板、腹板和底板承剪比计算结果较已有变截面剪应力计算结果精度最大可提高3.48%、3.43%和6.91%；悬臂梁各组件承剪比取决于荷载形式，梁端集中荷载作用下，顶、底板各自承剪比达到最大，分别为自由端的12.82%和固定端的60.81%；均布荷载作用下，腹板承剪比达到最大，为自由端的78.11%.

关键词

剪应力; 能量变分法; 有限梁段法; 承剪比

波形钢腹板组合箱梁因其有效地解决了混凝土腹板开裂问题而被广泛使用^［

1-3］.但波形钢腹板手风琴效应及剪切模量折减易引起剪切变形、剪切屈曲等问题^{［参考文献 4-6}4-6］，而变截面波形钢腹板组合箱梁的该类问题更为突出且应用范围更为广泛^{［参考文献 7-10}7-10］.因此，研究变截面波形钢腹板组合箱梁抗剪性能对该类桥梁的应用推广具有十分重要的力学意义.

在波形钢腹板组合箱梁的发展过程中，国内外众多学者对波形钢腹板组合箱梁的抗剪性能展开了一系列的研究.抗剪强度方面Luo等^［

11］研究了变截面波形钢腹板梁极限强度；Yi等^{［参考文献 12

百度学术}12］研究了波形钢腹板通用经典板屈曲；Papangelis等^{［参考文献 13

百度学术}13］研究了变截面波形钢腹板梁的抗剪承载力设计；Hassanein等^{［参考文献 14-16}14-16］研究了波形钢腹板梁强度设计.剪应力计算方面Moon等^{［参考文献 17

百度学术}17］通过公式代换将波形钢腹板整体屈曲荷载及合成剪切屈曲应力计算公式发展成了一般的经典板屈曲计算公式；Nie等^{［参考文献 18

百度学术}18］基于线弹性屈曲分析进行了广泛的参数研究，得出用于计算考虑三种不同剪切屈曲模式的波形钢腹板弹性剪切屈曲强度的简化公式；Zhou等^{［参考文献 7-10}7-10］在不考虑波形钢腹板翼板界面滑移的前提下考虑附加剪力，求得变截面波形钢腹板组合箱梁剪应力；Ye等^{［参考文献 19

百度学术}19］对超高性能混凝土-波形钢腹板预应力变截面组合箱梁的剪应力计算公式进行推导；与此同时，张玉元等^{［参考文献 20

百度学术}20］推导了利用各内力分量引起的剪应力分量叠加求解变截面薄壁箱梁剪应力的计算公式.承剪比方面Kadotani等^{［参考文献 21

百度学术}21］通过试验研究得出简支梁65%的剪力由波形钢腹板承担；Shitou等^{［参考文献 22

百度学术}22］指出翼板开裂导致波形钢腹板承剪比增大，但屈服后这一比例减小；Kato等^{［参考文献 23

百度学术}23］考虑了波形钢腹板与翼板之间的剪力分布，提出了波形钢腹板组合箱梁的设计方法.虽然已有文献对剪应力计算开展了比较深入的研究，但经典等截面剪应力计算对于腹板剪应力计算过于保守，变截面剪应力计算未对顶、底板承剪进行剥离，使得顶底板承剪无法量化分析；同时已有文献推导变截面薄壁箱梁及变截面波形钢腹板组合箱梁剪应力计算公式均通过求解各内力分量引起的剪应力分量并叠加而得，且波形钢腹板组合箱梁手风琴效应会使翼板、腹板以及组合箱梁截面转角产生差异，因此亟须一种能体现波形钢腹板组合箱梁特有手风琴效应的剪应力计算方法.

本文从翼板、腹板以及组合箱梁截面转角出发建立位移函数，基于能量变分法及有限梁段法结合弯矩等效，将顶板、底板及腹板承剪进行剥离并求解，在此基础上，分析悬臂梁在不同荷载工况下顶板、底板、钢腹板剪应力及承剪比.

1 传统剪应力计算公式

1.1 等截面箱梁剪应力计算

假设剪力全部由腹板承担且剪应力沿高度方向均匀分布时：

τ = \frac{Q}{A_{w}}

（1）

式中：τ为腹板竖向剪应力；Q为全截面竖向剪力；A_w为腹板承剪面积.

1.2 变截面箱梁剪应力计算

参考文献［

7，20］，取微段dx为研究对象，如图1所示，利用微元体平衡、剪应力互等定理可得变截面组合箱梁任意点的剪应力计算公式为：

τ = \frac{Q S_{a}}{I b_{3}} + \frac{N}{b_{3}} (\frac{S_{a} t a n ϕ}{I} + \frac{1}{A} \frac{d A_{a}}{d x} - \frac{A_{a}}{A^{2}} \frac{d A}{d x}) +

\frac{M}{I b_{3}} (\frac{d S_{a}}{d x} - \frac{S_{a}}{I} \frac{d I}{d x}) = τ_{Q} + τ_{N} + τ_{M}

（2）

式中：Q、N、M分别表示全截面剪力、轴力和弯矩；ϕ为截面形心连线的水平倾角；b₃为剪应力计算点的截面宽度；A_a为剪应力计算点至梁顶的截面面积；A、I分别为横截面面积及对形心轴的惯性矩；S_a为A_a对全截面形心轴的静矩；τ_Q、τ_N、τ_M分别为由剪力、轴力、弯矩引起的剪应力分量.

图1 微段平衡示意图

Fig.1 Micro segment equilibrium diagram

式（2）存在 $\frac{d A_{a}}{d x}$ ， $\frac{d A}{d x}$ ， $\frac{d S_{a}}{d x}$ ， $\frac{d I}{d x}$ 四项微分项，且随着计算位置的不同，A_a与S_a表达式也不相同，则可计算不同位置处剪应力.

2 顶板、底板、腹板剪力分离及计算

2.1 基本假定

1）结构处于弹性工作阶段，材料受力符合胡克定律；

2）弹性工作阶段拟平截面假定成立；

3）考虑波形钢腹板的手风琴效应，忽略波形钢腹板上的法向应变，并认为混凝土顶、底板主要贡献抗弯强度，波形钢腹板仅提供抗剪强度；

4）忽略钢腹板混凝土翼板连接界面处的剪切滑移；

5）剪应力沿高度方向均匀分布，不考虑顶、底板的有效分布宽度及剪切变形.

2.2 组合箱梁位移函数

波形钢腹板组合箱梁截面示意如图2所示.图中坐标原点o位于截面形心处，y、z为横截面的形心主惯性轴，b₁、b₂分别表示顶板、底板的1/2宽度，h_s、h_x分别表示截面形心到顶板底面、底板顶面的距离，h_w表示顶、底板间净距，z_s、z_x分别表示截面形心到顶、底板中心的距离；t₁、t₂分别为顶板、底板厚度.

图2 组合箱梁示意图

Fig .2 Composite box beam diagram

波形钢腹板组合箱梁承受任意荷载，考虑腹板手风琴效应，各组件变形如图3所示.图中：α（x）为顶、底板中心连线与z轴的夹角， $w^{'}$ （x）为顶、底板竖向挠曲转角，θ（x）为波形钢腹板挠曲转角，u、w分别表示组合箱梁的纵向和竖向位移.

图3 波形钢腹板变形模式示意图

Fig.3 Deformation mode diagram of corrugated steel web

由图2和图3可得，上下翼板的纵向位移为：

u_{s} (z, x) = z_{s} α (x) - (z + z_{s}) w^{'} (x)

（3）

u_{x} (z, x) = - z_{x} α (x) - (z - z_{x}) w^{'} (x)

（4）

由式（3）、式（4）可得：

θ (z) = \frac{u_{s} (- h_{s}, x) - u_{x} (h_{x}, x)}{h_{s} + h_{x}} =

λ α (x) + (1 - λ) w^{'} (x)

（5）

式中： $λ = (z_{s} + z_{x}) / (h_{s} + h_{x}) = h_{u} / h_{w}$ .

由式（3）~式（5）可得腹板任意点的纵向位移为：

u_{w} (z, x) = u_{s} (- h_{s}, x) - (z + h_{s}) θ (x) =

[z_{s} - λ (z + h_{s})] α (x) -

[(z_{s} - h_{s}) - χ (z + h_{s})] w^{'} (x)

（6）

式中： $χ = λ - 1$ .

2.3 控制微分方程

由式（3）~式（6）可得：

顶板正应变为：

ε_{x}^{t} = z_{s} α^{'} (x) - (z + z_{s}) w^{″} (x)

（7）

底板正应变为：

ε_{x}^{b} = - z_{x} α^{'} (x) - (z - z_{s}) w^{″} (x)

（8）

腹板剪应变为：

γ_{w} = λ [w^{'} (x) - α (x)]

（9）

应变能为：

\begin{array}{l} U = \frac{1}{2} \int_{l} \int_{A} {E_{c} [(ε_{x}^{t})^{2} + (ε_{x}^{b})^{2}] + G_{w} (γ_{w})^{2}] d A d x = \\ \frac{1}{2} E_{c} \int_{l} (I_{y 0} {θ^{'}}^{2} + I_{y c} {w^{″}}^{2}) d x + \end{array}

\frac{1}{2} G_{w} \int_{l} λ^{2} A_{w} (w^{'} {- α)}^{2} d x

（10）

其中：

I_{y 0} = z_{s}^{2} A_{s} + z_{x}^{2} A_{x}

（11）

I_{y c} = \frac{1}{12} t_{1}^{2} A_{s} + \frac{1}{12} t_{2}^{2} A_{x}

（12）

G_{w} = \frac{a + b}{a + c} G_{s}, A_{w} = h_{w} t_{w}

（13）

式中：A_s、A_x和A_w分别为组合箱梁顶板、底板和腹板的面积；E_c为翼板混凝土弹性模量；G_s为钢材的剪切模量；G_w表示组合箱梁钢腹板的换算剪切模量；a、b、c为波形参数，如图4所示.

图4 波形参数示意图

Fig.4 Schematic diagram of waveform parameters

总势能为：

Π = \frac{1}{2} E_{c} \int_{l} (I_{y 0} {α^{'}}^{2} + I_{y c} {w^{″}}^{2}) d x +

\frac{1}{2} G_{w} \int_{l} λ^{2} A_{w} (w^{'} {- α)}^{2} d x - \int_{0}^{l} q (x) w (x) d x

（14）

由最小势能原理δΠ=0，可得：

\begin{array}{l} δ Π = - \int_{0}^{l} [E_{c} I_{y 0} α^{″} + λ^{2} G_{w} A_{w} (w^{'} - α)] δ α d x + \\ \int_{0}^{l} [E_{c} I_{y c} {w^{″}}^{″} - λ^{2} G_{w} A_{w} (w^{″} - α^{'}) - q] δ w d x + \\ E_{c} I_{y 0} α^{'} δ α |_{0}^{l} + E_{c} I_{y c} w^{″} δ w^{'} |_{0}^{l} + \\ [λ^{2} G_{w} A_{w} (w^{'} - α) - E_{c} I_{y c} w^{‴}] δ w |_{0}^{l} \end{array}

（15）

控制微分方程为：

E_{c} I_{y c} {w^{″}}^{″} - λ^{2} G_{w} A_{w} (w^{″} - α^{'}) = q

（16）

E_{c} I_{y 0} α^{″} + λ^{2} G_{w} A_{w} (w^{'} - α) = 0

（17）

边界条件为：

E_{c} I_{y c} w^{″} δ w^{'} |_{0}^{l} = 0

（18）

E_{c} I_{y 0} α^{'} δ α |_{0}^{l} = 0

（19）

[λ^{2} G_{w} A_{w} (w^{'} - α) - E_{c} I_{y c} w^{‴}] δ w |_{0}^{l} = 0

（20）

由边界条件可知，顶板和底板变形产生的弯矩为：

M_{α} = E_{c} I_{y 0} α^{'}

（21）

顶板和底板绕自身中心的弯矩为：

M_{w} = E_{c} I_{y c} w^{″}

（22）

截面总剪力为：

Q = λ^{2} G_{w} A_{w} (w^{'} - α) - E_{c} I_{y c} w^{‴}

（23）

其中，波形钢腹板承担的剪力为：

Q_{w} = \int_{A} G_{w} γ_{w} d A = λ G_{w} A_{w} (w^{'} - α)

（24）

顶板和底板承担的剪力为：

Q_{c} = λ χ G_{w} A_{w} (w^{'} - α) - E_{c} I_{y c} w^{‴}

（25）

波形钢腹板剪应力为：

τ_{w} = \frac{Q_{w}}{A_{w}}

（26）

求出w（x）、α（x）、θ（x），便可求解各剪力.

2.4 顶、底板承剪分离

如图5所示，Q为全截面剪力，Q_cs1、Q_w1、Q_cx1分别为顶板、腹板、底板本身所受剪力，Q_cs2、Q_w2、Q_cx2分别为顶板、腹板、底板附加剪力，且Q_cs2+Q_w2+Q_cx2=0.

图5 箱梁受剪示意图

Fig.5 Shear diagram of box beam

Q_cs1、Q_w1、Q_cx1计算方法如下：

Q_{w 1} = \frac{Q S}{I b_{0}} \cdot A_{w}

（27）

Q_{c s 1} = (Q - Q_{w 1}) \cdot \frac{A_{s}}{A_{s} + A_{x}}

（28）

Q_{c x 1} = (Q - Q_{w 1}) \cdot \frac{A_{x}}{A_{s} + A_{x}}

（29）

式中：S为横截面上该点以外部分的面积对中性轴的静距；I为横截面对中性轴的惯性矩；b₀为截面有效宽度.

参考文献［

7］，组合箱梁弯矩M可以等效为一对作用在顶、底板截面形心上的平衡力系｛F_t，F_c｝，如图6所示.

图6 等效弯矩示意图

Fig.6 Diagram of equivalent bending moment

则：

F_{t} = \frac{M}{0.5 t_{1} + h_{w} + 0.5 t_{2}} = - F_{c}

（30）

Q_{c x 2} = F_{c} t a n φ

（31）

Q_{c x} = Q_{c x 1} + Q_{c x 2}

（32）

Q_{c s} = Q - Q_{w} - Q_{c x}

（33）

式中：φ为底板形心轴与水平线的夹角；Q_cx为底板所受剪力；Q_cs为顶板所受剪力.

顶板剪应力为：

τ_{s} = \frac{Q_{c s}}{A_{s}}

（34）

底板剪应力为：

τ_{x} = \frac{Q_{c x}}{A_{x}}

（35）

2.5 承剪比

承剪比为：

ρ = \frac{Q_{i}}{Q}

（36）

式中：Q_i为各组件所受剪力.

3 转角位移求解

为求解w（x）、α（x）、θ（x），建立2节点6自由度的梁段单元，如图7所示.

图7 有限梁段模型

Fig.7 Finite beam segment model

位移向量δ^e及其所对应的力向量F^e为：

δ^{e} = [w_{i} w_{i}^{'} α_{i} w_{j} w_{j}^{'} α_{j}]^{T}

（37）

F^{e} = [Q_{i} M_{w i} M_{α i} Q_{j} M_{w j} M_{α j}]^{T}

（38）

式中：w_i和w_j分别表示梁段单元i端和j端的竖向位移； $w_{i}^{'}$ 和 $w_{j}^{'}$ 分别表示单元i端和j端顶、底板截面形心的转角；α_i和α_j分别表示单元i端和j端顶、底板绕各板中心的转角；Q_i和Q_j分别表示单元i端和j端截面的全部剪力；M_wi和M_wj分别表示单元i端和j端顶、底板绕截面形心转动产生的弯矩；M_αi和M_αj分别表示单元i端和j端顶、底板绕各自板中心转动产生的弯矩.

由式（16）和式（17）得：

w^{(6)} - k^{2} w^{(4)} = \frac{q}{k_{0}} + \frac{q^{″}}{E_{c} I_{y c}}

（39）

其中：

k^{2} = \frac{λ^{2} G_{w} A_{w} I_{y}}{E_{c} I_{y c} I_{y 0}}

（40）

k_{0} = \frac{E_{c}^{2} I_{y c} I_{y 0}}{λ^{2} G_{w} A_{w}}

（41）

微分方程（39）的解为：

w (x) = C_{1} + C_{2} x + C_{3} x^{2} + C_{4} x^{3} +

C_{5} c o s h k x + C_{6} s i n h k x + C_{0} x^{4}

（42）

式中：C₀x⁴为均布荷载作用时的特解，C₁~C₆为积分常数.

由式（18）、式（19）得：

\frac{E_{c} I_{y c} I_{y 0}}{λ^{2} G_{w} A_{w}} w^{(4)} - E_{c} I_{y} w^{″} - \frac{q (x) E_{c} I_{y 0}}{λ^{2} G_{w} A_{w}} = M (x)

（43）

对无载梁段：

M (x) = - 2 E_{c} I_{y} C_{3} - 6 E_{c} I_{y} C_{4} x

（44）

由式（20）可得：

α (x) = w^{'} - \frac{E_{c} I_{y} w^{‴}}{λ^{2} G_{w} A_{w}} - \frac{Q}{λ^{2} G_{w} A_{w}}

（45）

求解α（x）可得：

α (x) = C_{2} + 2 C_{3} x + C_{4} (6 ζ + 3 x^{2}) +

C_{5} ς s i n h k x + C_{6} ς c o s h k x

（46）

其中：

ζ = \frac{E_{c} I_{y 0}}{λ^{2} G_{w} A_{w}}

（47）

ς = - k \frac{I_{y c}}{I_{y c}}

（48）

由式（39）~式（41）可得：

Q = M^{'} (x) = 6 E_{c} I_{y} C_{4}

（49）

M_{w i} = - E_{c} I_{y c} (2 C_{3} + 6 C_{4} x +

k^{2} C_{5} c o s h k x + k^{2} C_{6} s i n h k x)

（50）

M_{α i} = - E_{c} I_{y 0} (2 C_{3} + 6 C_{4} x +

k ς C_{5} c o s h k x + k ς C_{6} s i n h k x)

（51）

将式（37）改写为：

δ^{e} = A_{E} C

（52）

向量C及A_E为：

C = [C_{1} C_{2} C_{3} C_{4} C_{5} C_{6}]^{T}

（53）

A_{E} = [\begin{matrix} 1 & 0 & 0 & 0 & 1 & 0 \\ 0 & 1 & 0 & 0 & 0 & k \\ 0 & 1 & 0 & 6 ζ & 0 & ς \\ 1 & l & l^{2} & l^{3} & c o s h k l & s i n h k l \\ 0 & 1 & 2 l & 3 l^{2} & k s i n h k l & k c o s h k l \\ 0 & 1 & 2 l & 6 ζ + 3 l^{2} & ς s i n h k l & ς c o s h k l \end{matrix}]

（54）

将式（38）改写为：

F^{e} = B_{E} C

（55）

B_{E} = [\begin{matrix} 0 & 0 & 0 & 6 E_{c} I_{y} & 0 & 0 \\ 0 & 0 & - 2 E_{c} I_{y c} & 0 & - k^{2} E_{c} I_{y c} & 0 \\ 0 & 0 & - 2 E_{c} I_{y 0} & 0 & - k ς E_{c} I_{y 0} & 0 \\ 0 & 0 & 0 & - 6 E_{c} I_{y} & 0 & 0 \\ \begin{array}{l} 0 \end{array} & \begin{array}{l} 0 \end{array} & \begin{array}{l} 2 E_{c} I_{y c} \end{array} & \begin{array}{l} 6 E_{c} I_{y c} l \end{array} & \begin{array}{l} k^{2} E_{c} I_{y c} \\ s i n h k l \end{array} & \begin{array}{l} k^{2} E_{c} I_{y c} \\ c o s h k l \end{array} \\ \begin{array}{l} 0 \end{array} & \begin{array}{l} 0 \end{array} & \begin{array}{l} 2 E_{c} I_{y 0} \end{array} & \begin{array}{l} 6 E_{c} I_{y 0} l \end{array} & \begin{array}{l} k ς E_{c} I_{y 0} \\ s i n h k l \end{array} & \begin{array}{l} k ς E_{c} I_{y 0} \\ c o s h k l \end{array} \end{matrix}]

（56）

将系数矩阵代入式（52）可得：

F^{e} = B_{E} {A_{E}}^{- 1} δ^{e} = K^{e} δ^{e}

（57）

单元任一点处的位移为：

u^{δ} = K_{u} {A_{E}}^{- 1} δ^{e} = N δ^{e} = [w w^{'} {α]}^{T}

（58）

位移形函数Ν 为：

Ν = [\begin{array}{l} N_{11} N_{12} N_{13} N_{14} N_{15} N_{16} \\ N_{21} N_{22} N_{23} N_{24} N_{25} N_{26} \\ N_{31} N_{32} N_{33} N_{34} N_{35} N_{36} \end{array}]

（59）

K_{u} = [\begin{matrix} 1 & x & x^{2} & x^{3} & c o s h k x & s i n h k x \\ 0 & 1 & 2 x & 3 x^{2} & k s i n h k x & k c o s h k x \\ 0 & 1 & 2 x & 6 ζ + 3 x^{2} & ς s i n h k x & ς c o s h k x \end{matrix}]

（60）

根据虚功原理，可求得均布荷载作用时的等效节点荷载为：

F^{e q} = \int_{0}^{l} q N^{T} d x =

- q \int_{0}^{l} [N_{11} N_{12} N_{13} N_{14} N_{15} N_{16}]^{T} d x

（61）

F_{q} = [F_{1}^{q} F_{2}^{q} F_{3}^{q} F_{4}^{q} F_{5}^{q} F_{6}^{q}]^{T}

（62）

式中：

\{\begin{array}{l} F_{1}^{q} = F_{4}^{q} = \frac{q l}{2} \\ F_{2}^{q} = - F_{5}^{q} = \frac{q l^{2} υ}{12} \\ F_{3}^{q} = - F_{6}^{q} = \frac{q l^{2} (1 - υ)}{12} \end{array}

（63）

υ = \frac{6 k l + ς k l^{2} + 12}{(ς - k) k l^{2}} + \frac{12 e^{k l}}{(1 - e^{k l}) (ς - k) l}

（64）

转角位移的求解，采用计算软件编写程序进行.

4 算例分析

4.1 几何参数

参考文献［

7］中计算模型，其中悬臂梁的跨度30 m，顶板宽4 m，底板宽2 m，梁高由固定端4 m线性变化至自由端2 m，顶板厚度0.3 m，底板厚由固定端0.5 m线性变化至自由端0.3 m，波形钢腹板为1000型（CSW1000），钢腹板厚度0.01 m.如图8所示.

（a）纵向示意图

（b）截面示意图

图8 变截面梁几何尺寸（单位：mm）

Fig.8 Geometric dimension of variable section beam （unit：mm）

4.2 有限元模型

根据模型尺寸，采用ANSYS软件建立三维有限元分析模型，如图9所示.其中顶板、底板采用实体单元SOLID65模拟，波形钢腹板采用壳单元SHELL63模拟，顶板、底板混凝土材料弹性模量E=3.45×10⁴ MPa，泊松比μ=0.2；腹板钢材弹性模量E=2.1×10⁵ MPa，泊松比μ=0.3.顶板、底板实体单元与腹板壳单元采用共节点连接，固定端（x=0）约束平动和转动自由度.

图9 有限元模型示意图

Fig.9 Diagram of finite element model

4.3 转角位移求解

采用MATLAB计算软件编写对应的梁段分析程序，进行转角位移的求解.模型梁划分30个单元，单元尺寸1 m，计算流程如图10所示.

图10 计算流程图

Fig.10 Calculation flowchart diagram

4.4 剪应力分析

4.4.1 均布荷载

在悬臂梁上施加均布荷载q=35 kN/m，采用对称加载，如图11所示.为减少边界效应带来的误差，将固定端l号截面和自由端2号截面向内侧移动1 m.同时，对式（2）、式（26）和有限元模型的计算结果进行对比，验证模型的可靠性.

图11 均布荷载加载示意图

Fig.11 Diagram of uniformly load

将本文推导的计算结果与式（1）、式（2）、有限元计算结果作对比，同时以有限元计算结果作为参考进行误差分析，分析结果如图12~图14所示.

图12 均布荷载作用下的顶板剪应力

Fig.12 Shear stress of top flange under uniform load

图13 均布荷载作用下的腹板剪应力

Fig.13 Shear stress of web under uniform load

图14 均布荷载作用下的底板剪应力

Fig.14 Shear stress of bottom flange under uniform load

由图12可知，本文式（34）计算结果介于式（1）与式（2）与有限元结果结果之间，且与有限元结果吻合良好，误差在6%以内；式（1）、式（2）结果与有限元结果差值较大，式（2）最大差值达到2倍以上甚至更大.

由图13可知，本文式（26）计算结果介于公式（2）与有限元结果之间，且与有限元结果的误差在5%以内.式（1）忽略翼板对竖向剪力的贡献，使得截面1处腹板剪应力达到了有限元结果的1.59倍，这与实际抗力组成存在较大差异，随着跨度变化，其差值由截面1的58.95%减小到截面2的28.42%.截面1处式（2）与有限元结果差值为7.29%，式（26）为4.17%.计算结果中式（1）最大，其次为式（2）、式（26）、有限元.

由图14可知，本文式（35）计算结果与有限元以及式（2）结果变化趋势大致相同，且与有限元值吻合良好，误差在6%以内.

4.4.2 集中荷载

在悬臂梁端部作用一竖向集中力P=500 kN，加载示意图如图15所示，计算结果如图16~图18所示.

图15 集中荷载加载示意图

Fig.15 Diagram of concentrated load

图16 集中荷载作用下的顶板剪应力

Fig.16 Shear stress of top flange under concentrated load

图17 集中荷载作用下的腹板剪应力

Fig.17 Shear stress of web under concentrated load

图18 集中荷载作用下的底板剪应力

Fig.18 Shear stress of bottom flange under concentrated load

由图16可知，本文式（34）计算结果介于式（1）与式（2）结果之间，且与有限元结果吻合良好，式（2）与有限元差值较大.

由图17可知，本文式（26）计算结果介于有限元与式（2）结果之间，三者吻合良好，且与有限元结果误差在5%以内.计算结果中式（1）最大，其次为式（2）、式（26）、有限元.随着跨度的变化，式（1）与有限元差值与有限元计算结果比值由截面1处的2.44减小到截面2处的1.31.

由图18可知，本文式（35）计算结果与有限元结果以及式（2）结果变化趋势大致相同，且与有限元结果吻合良好，误差在6%以内.

4.5 承剪比分析

为了更直观地表达各组件的承剪情况，用承剪比来进一步说明各组件对抗剪的贡献，选取4种荷载工况（工况Ⅰ q=35 kN/m，工况Ⅱ P=500 kN，工况Ⅲ q=17 kN/m，工况Ⅳ P=300 kN），其中工况Ⅰ与工况Ⅱ在截面1处弯矩近似相等；工况Ⅱ与工况Ⅲ在截面1处剪力近似相等；工况Ⅰ与工况Ⅲ荷载形式相同，大小不同；工况Ⅱ与工况Ⅳ荷载形式相同，大小不同.利用式（36）求得顶、底板、波形钢腹板的承剪比.不同工况下顶、底板和波形钢腹板承剪比如图19~图21及表1所示.

图19 顶板承剪比

Fig.19 Shear capacity ratio of top flange

图20 腹板承剪比

Fig.20 Shear capacity ratio of web

图21 底板承剪比

Fig.21 Shear capacity ratio of bottom flange

表1 工况Ⅰ~Ⅳ 承剪比

Tab.1 Condition Ⅰ~Ⅳ shear capacity ratio

	工况	承剪比/%
		式（1）		式（2）		ANSYS		式（33）/（24）/（32）
		截面1	截面2	截面1	截面2	截面1	截面2	截面1	截面2
顶板	Ⅰ	0	0	5.16	14.63	1.41	10.95	1.90	11.15
	Ⅱ	0	0	1.39	14.25	0.16	10.19	0.12	12.82
	Ⅲ	0	0	5.16	14.63	1.41	10.95	1.90	11.15
	Ⅳ	0	0	1.39	14.25	0.16	10.19	0.12	12.82
腹板	Ⅰ	100	100	65.49	81.54	60.53	76.12	62.53	78.11
	Ⅱ	100	100	41.94	80.10	39.06	74.12	39.07	76.71
	Ⅲ	100	100	65.49	81.54	60.53	76.12	62.53	78.11
	Ⅳ	100	100	41.94	80.10	39.06	74.12	39.07	76.71
底板	Ⅰ	0	0	29.35	3.83	38.06	12.93	35.57	10.74
	Ⅱ	0	0	56.67	5.65	60.78	15.69	60.81	10.47
	Ⅲ	0	0	29.35	3.83	38.06	12.93	35.57	10.74
	Ⅳ	0	0	56.67	5.65	60.78	15.69	60.81	10.47

由图19~图21、表1及本文计算公式可知：

1）工况Ⅰ、Ⅲ各截面承剪比相同，工况Ⅱ、Ⅳ各截面承剪比相同，但工况Ⅰ、Ⅲ与Ⅱ、Ⅳ不同，由此说明组合悬臂梁承剪比取决于荷载形式.

2）随着跨度增加，顶板和腹板的承剪比逐渐增大而底板的承剪比逐渐减小.梁端集中荷载作用下，顶、底板各自承剪比达到最大，分别为自由端的12.82%和固定端的60.81%；均布荷载作用下，腹板承剪比达到最大，为自由端的78.11%.

3）荷载作用下，顶板承剪比与式（2）结果最大差值为工况Ⅰ、Ⅲ下截面2处的3.48%，腹板最大差值为工况Ⅰ、Ⅲ下截面2处3.43%，底板最大差值为工况Ⅰ、Ⅲ下截面2处6.91%.结合有限元模型，本文计算方法可使顶板计算精度最大提高3.48%，腹板计算精度最大提高3.43%，底板计算精度最大提高6.91%.

5 结论

本文从翼板、腹板以及组合箱梁截面转角差异出发，求解变截面波形钢腹板悬臂梁顶板、底板、腹板剪应力及承剪比，结合算例分析得出如下结论：

1）采用等截面剪应力计算公式高估腹板承剪的最大误差达到2.44倍.

2）变截面波形钢腹板悬臂梁由固定端到自由端，顶板、腹板承剪比逐渐增大，底板承剪比逐步减小.

3）变截面波形钢腹板悬臂梁顶板、底板及腹板承剪比取决于荷载形式；梁端集中荷载作用下，顶、底板各自承剪比达到最大，分别为自由端的12.82%和固定端的60.81%；均布荷载作用下，腹板承剪比达到最大，为自由端的78.11%.

4）从截面转角出发可使均布荷载下顶板、腹板和底板承剪比的计算精度最多提高3.48%、3.43%和6.91%.

参考文献

李丽园. 考虑剪切变形影响的波形钢腹板组合箱梁挠曲力学特性及其试验研究［D］. 兰州：兰州交通大学， 2019. [百度学术]

LI L Y. Vertical flexural mechanical characteristics and experimental study of composite box girder with corrugated steel web considering shear deformation［D］. Lanzhou： Lanzhou Jiaotong University， 2019. （in Chinese） [百度学术]

JIANG R J，WU Q M，XIAO Y F，et al．The shear lag effect of composite box girder bridges with corrugated steel webs［J］．Structures，2023，48：1746-1760． [百度学术]

INAAM Q， UPADHYAY A．Accordion effect in bridge girders with corrugated webs［J］. Journal of Constructional Steel Research， 2022， 188： 107040． [百度学术]

OH J Y，LEE D H，KIM K S．Accordion effect of prestressed steel beams with corrugated webs［J］. Thin-Walled Structures， 2012，57：49-61． [百度学术]

SU X L，ZHOU M．Shear force transfer efficiency of tapered bridges with trapezoidal corrugated steel webs［J］．Journal of Constructional Steel Research， 2023， 211： 108135． [百度学术]

LEBLOUBA M，JUNAID M T，BARAKAT S，et al．Shear buckling and stress distribution in trapezoidal web corrugated steel beams［J］．Thin-Walled Structures， 2017， 113： 13-26． [百度学术]

ZHOU M，ZHANG J D，ZHONG J T，et al．Shear stress calculation and distribution in variable cross sections of box girders with corrugated steel webs［J］. Journal of Structural Engineering， 2016， 142（6）： 04016022． [百度学术]

ZHOU M，YANG D Y，ZHANG J D，et al．Stress analysis of linear elastic non-prismatic beams with corrugated steel webs［J］．Thin-Walled Structures， 2017， 119： 653-661． [百度学术]

ZHOU M，LIAO J C ， ZHONG J T， et al. Unified calculation formula for predicting the shear stresses in prismatic and non-prismatic beams with corrugated steel webs［J］. Structures， 2021，29： 507-518． [百度学术]

SU X L，ZHOU M，TANWARI K A．Shear behavior of prismatic and nonprismatic beams with trapezoidal steel web corrugations considering secondary internal force induced by external prestressing［J］．Structures， 2023， 47： 638-647． [百度学术]

LUO R， EDLUND B. Ultimate strength of girders with trapezoidally corrugated webs under patch loading［J］. Thin-Walled Structures，1996， 24（2）： 135-156． [百度学术]

YI J，GIL H，YOUM K，et al．Interactive shear buckling behavior of trapezoidally corrugated steel webs［J］. Engineering Structures，2008， 30（6）： 1659-1666． [百度学术]

PAPANGELIS J，TRAHAIR N，HANCOCK G．Direct strength method for shear capacity of beams with corrugated webs［J］．Journal of Constructional Steel Research， 2017， 137： 152-160． [百度学术]

HASSANEIN M F，KHAROOB O F．Shear buckling behavior of tapered bridge girders with steel corrugated webs［J］. Engineering Structures， 2014， 74： 157-169． [百度学术]

HASSANEIN M F， KHAROOB O F. Linearly tapered bridge girder panels with steel corrugated webs near intermediate supports of continuous bridges［J］. Thin-Walled Structures， 2015， 88： 119-128． [百度学术]

ZEVALLOS E，HASSANEIN M F，REAL E，et al. Shear evaluation of tapered bridge girder panels with steel corrugated webs near the supports of continuous bridges［J］．Engineering Structures， 2016， 113： 149-159． [百度学术]

MOON J， YI J， CHOI B H， et al. Shear strength and design of trapezoidally corrugated steel webs［J］. Journal of Constructional Steel Research， 2009， 65（5）： 1198-1205． [百度学术]

NIE J G， ZHU L， TAO M X， et al. Shear strength of trapezoidal corrugated steel webs［J］. Journal of Constructional Steel Research， 2013， 85： 105-115． [百度学术]

YE M，LI L F，YOO D Y，et al．Shear performance of prestressed composite box beams with ultra-high-performance concrete and corrugated steel webs under different loading conditions［J］．Thin-Walled Structures， 2023， 186： 110675． [百度学术]

张玉元，张元海，张慧．基于换算剪力的变截面箱梁弯曲剪应力计算方法［J］．东南大学学报（自然科学版）， 2021， 51（1）：16-22． [百度学术]

ZHANG Y Y， ZHANG Y H， ZHANG H．Calculation method of bending shear stress of variable section box girder based on converted shear force［J］. Journal of Southeast University （Natural Science Edition）， 2021， 51（1）： 16-22．（in Chinese） [百度学术]

KADOTANI T， AOKI K， ASHIZUKA K， et al. Shear buckling behavior of prestressed concrete girders with corrugated steel webs［C］//Proc. 1st Fib Congress Sessions： Composite Structures. Osaka， Japan， 2002，： 269-276. [百度学术]

SHITOU K，NAKAZONO A，SUZUKI N，et al．Experimental research on shear behavior of corrugated steel web bridge［J］．Doboku Gakkai Ronbunshuu A， 2008， 64（2）： 223-234． [百度学术]

KATO H， NISHIMURA N. Practical analysis methods for continuous girder and cable stayed bridges composed of beams with corrugated steel webs［J］. Structural Engineering， 2004， 21（2）： 207S-222S． [百度学术]

作者稿件一经被我刊录用，如无特别声明，即视作同意授予我刊论文整体的全部复制传播的权利，包括但不限于复制权、发行权、信息网络传播权、广播权、表演权、翻译权、汇编权、改编权等著作使用权转让给我刊，我刊有权根据工作需要，允许合作的数据库、新媒体平台及其他数字平台进行数字传播和国际传播等。特此声明。

关闭

首页

期刊简介

编委会

作者中心

下载中心

学术道德

常见问题

版权声明

联系我们

English

变截面波形钢腹板悬臂梁抗剪性能研究 PDF

摘要

关键词

1 传统剪应力计算公式

1.1 等截面箱梁剪应力计算

1.2 变截面箱梁剪应力计算

2 顶板、底板、腹板剪力分离及计算

2.1 基本假定

2.2 组合箱梁位移函数

2.3 控制微分方程

2.4 顶、底板承剪分离

2.5 承剪比

3 转角位移求解

4 算例分析

4.1 几何参数

4.2 有限元模型

4.3 转角位移求解

4.4 剪应力分析

4.5 承剪比分析

5 结论

参考文献

变截面波形钢腹板悬臂梁抗剪性能研究 PDF

摘要

关键词

1 传统剪应力计算公式

1.1 等截面箱梁剪应力计算

1.2 变截面箱梁剪应力计算

2 顶板、底板、腹板剪力分离及计算

2.1 基本假定

2.2 组合箱梁位移函数

2.3 控制微分方程

2.4 顶、底板承剪分离

2.5 承剪比

3 转角位移求解

4 算例分析

4.1 几何参数

4.2 有限元模型

4.3 转角位移求解

4.4 剪应力分析

4.5 承剪比分析

5 结 论

参考文献

5 结论