高强钢端板连接节点弯矩⁃转角曲线数学模型

摘要

为预测高强钢端板连接节点在常温、火灾下及火灾后的弯矩⁃转角曲线，以四参数指数模型为基础，基于组件法和等效T型连接，提出了初始刚度、屈服后刚度和抗弯极限承载力的计算方法，将计算结果代入四参数指数模型得到预测高强钢端板连接节点弯矩⁃转角曲线的方法；与足尺试验中Q690和Q960高强钢端板连接节点在常温、火灾下及火灾后的共10组试验结果进行对比。结果表明，抗弯极限承载力的计算结果与试验值较吻合，其中Q960的相对误差在5%以内，初始刚度的计算结果比欧洲钢结构设计规范（EN 1993—1‒8）的结果更接近试验值，预测的Q960高强钢端板连接节点弯矩⁃转角曲线与试验结果吻合，而Q690预测结果偏于安全。

关键词

高强钢; 端板连接节点; 弯矩⁃转角曲线; 数学模型

1 研究背景

钢结构梁柱端板连接节点具有较高的承载力、一定的延性和较好的抗震耗能能力，是多高层钢结构与钢结构厂房中常用的装配式节点形式之一^[1]。它属于半刚性节点，常采用组件法^[2]研究。我国和欧洲现行钢结构设计规范都建议采用节点的弯矩⁃转角（M⁃θ）曲线作为半刚性框架结构设计分析的依据^[3,4]。端板连接节点的典型M⁃θ曲线如图1所示。节点力学性能的主要指标包括：①初始刚度k_ini；②抗弯承载力M_y；③抗弯承载力对应的转角θ_y；④屈服后刚度k_p；⑤抗弯极限承载力M_max；⑥抗弯极限承载力对应的转角θ_max；⑦极限转角Ф_c。

图1 端板连接节点的M⁃θ曲线

Fig.1 M⁃θ curve of endplate connection

我国钢结构设计标准^[4]规定：在内力分析时，应假定半刚性连接节点的M⁃θ曲线，并在节点设计时，保证其构造与假定的M⁃θ曲线符合。因此，探究半刚性节点M⁃θ曲线的数学表达，对于内力分析和节点设计非常必要。最简单的数学表达莫过于线性模型，其包括单线性、双线性和多线性3种形式。单线性模型采用节点的初始刚度k_ini表征节点全部加载范围的力学特性；双线性模型在某一弯矩处，采用更小的刚度代替连接的初始刚度^[5,6,7,8,9]；多线性模型则采用一组直线来逼近节点非线性的M⁃θ曲线。双线性和多线性模型更接近节点的实际力学行为，但在转折点处刚度存在突变，容易造成计算不收敛。总的来说，线性模型过于简单，很难准确描述M⁃θ曲线的非线性段，即图1中的BE段，而OA为线性段。

为寻求节点M⁃θ曲线更准确的数学表达，Frye等^[10]提出了奇次方多项式数学模型；Jones等^[11,12]采用B样条曲线对节点的M⁃θ曲线进行拟合，保证了曲线在交点处一、二阶导数连续；Kishi等^[13]和Attiogbe等^[14]分别提出了三参数和四参数幂函数模型；郭兵等^[15,16]提出了两参数指数模型；Wu等^[17]提出了三参数指数模型，来描述顶、底角钢连接或带双腹板角钢的顶、底角钢连接的M⁃θ曲线；Yee等^[18]针对端板无加劲肋的外伸式端板连接提出了四参数指数模型，参数分别为节点的初始刚度k_ini、屈服后刚度k_p、抗弯极限承载力M_max和形状参数c，如下式：

$$M=M_{\mathrm{m}\mathrm{a}\mathrm{x}}\left\{\mathrm{1}-\mathrm{e}\mathrm{x}\mathrm{p}\left[\frac{-\left(k_{\mathrm{i}\mathrm{n}\mathrm{i}}-k_{\mathrm{p}}+c\theta \right)\theta }{M_{\mathrm{m}\mathrm{a}\mathrm{x}}}\right]\right\}+k_{\mathrm{p}}\theta$$

(1)

因此，现有研究针对节点M⁃θ曲线提出的数学模型主要有：线性模型、多项式模型、B样条模型、幂函数模型和指数模型。

目前，对于高强钢的定义，欧洲钢结构设计规范Eurocode 3（EC3）^[19]定义高强钢的名义屈服强度不小于460 MPa；澳大利亚规范AS 4100^[20]定义高强钢的名义屈服强度大于450 MPa；香港规范^[21]定义高强钢的名义屈服强度在460~690 MPa之间；我国国家标准《钢分类》^[22]定义屈服强度大于等于420 MPa的钢材为高强钢。

高强钢已经在国内外土木工程中得到较广泛应用^[23]，但因其力学性能与普通钢材存在差异，高强钢端板连接节点的力学性能也与普通钢端板连接节点不同。另一方面，高温会使钢材性能发生退化，火灾对钢结构建筑的影响一直是钢结构研究领域的一个关键问题^[23]。然而，上述已有研究中节点M⁃θ曲线的数学模型仅针对普通钢常温时的情况。Al⁃Jabri等^[24]对普通钢端板连接节点高温试验研究成果进行了总结，并采用Ramberg⁃Osgood模型拟合了不同温度下节点M⁃θ曲线；王卫永等^[25]在Yee等^[18]研究的基础上，通过引入钢材高温下力学性能折减系数，以式(1)为基础，提出了普通钢端板连接节点火灾下M⁃θ曲线的数学模型。然而，目前高强钢端板连接节点在常温、火灾下及火灾后M⁃θ曲线的数学模型尚未见报道。

本文提出了高强钢端板连接节点在常温、火灾下及火灾后的初始刚度k_ini、屈服后刚度k_p和抗弯极限承载力M_max的计算方法，与作者前期的试验结果^[26,27,28,29]进行对比，验证了所提k_ini和M_max计算方法的准确性；将计算得到的k_ini、k_p、M_max和c代入式(1)，得到节点M⁃θ曲线的预测结果，将其与试验结果作对比，验证了k_p和c计算方法的准确性，以及采用式(1)预测高强钢端板连接节点M⁃θ曲线的合理性。

2 初始刚度计算方法

2.1　各组件初始刚度计算方法

根据组件法，将节点拆分为多个基本组件，每个组件采用线性或非线性弹簧模拟，通过弹簧的串、并联组合计算，获得节点整体的力学行为^[2]，如图2所示。EN 1993—1‒8^[3]采用等效T型连接模拟端板连接节点受拉区的力学行为，根据EN 1993—1‒8^[3]，T型连接的失效模式有3种：翼缘在螺栓处和腹板处屈服（失效模式1）、翼缘在腹板处屈服的同时螺栓拉断（失效模式2）和螺栓拉断（失效模式3），如图3所示。图中，T_y为外荷载作用，Q为撬力，M_p为端板的全截面塑性抗弯承载力，B为螺栓受到的拉力，B_y为螺栓的屈服承载力，n为螺栓轴线至翼缘板边缘的距离，m为螺栓轴线至翼缘板根部的距离。为计算端板连接节点的初始转动刚度k_ini，需求得以下组件的初始刚度：①柱腹板抗压刚度k_cw,C；②柱腹板抗剪刚度k_cw,v；③柱腹板抗拉刚度k_cw,T；④柱翼缘抗弯刚度k_cf；⑤端板抗弯刚度k_ep；⑥螺栓抗拉刚度k_b。

图2 初始刚度的计算模型

Fig.2 Calculation model for initial stiffness

图3 T型连接的失效模式

Fig.3 Failure mode of T-stubs

2.1.1　柱腹板抗压刚度k_cw,C

对于无横向加劲肋的柱腹板，EN 1993—1‒8^[3]给出其抗压刚度k_cw,C的计算方法为

$$k_{\mathrm{c}\mathrm{w},\mathrm{C}}=\frac{\mathrm{0.7}E{b}_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{e}\mathrm{f}\mathrm{f}}{t}_{\mathrm{c}\mathrm{w}}}{h_{\mathrm{c}\mathrm{w}}}$$

(2)

式中：E为材料的弹性模量；h_cw为柱腹板的高度，h_cw=h_c-2t_cf；h_c为柱截面高度；t_cf为柱翼缘厚度；t_cw为柱腹板厚度；b_cw,C,eff为柱腹板受压区的有效高度，若为热轧型钢柱，其计算方法为

$$b_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{e}\mathrm{f}\mathrm{f}}=t_{\mathrm{b}\mathrm{f}}+\mathrm{2}{t}_{\mathrm{e}\mathrm{p}}+\mathrm{2}\sqrt[]{\mathrm{2}}{h}_{\mathrm{f},\mathrm{e}\mathrm{p}}+\mathrm{5}\left(t_{\mathrm{c}\mathrm{f}}+r_{\mathrm{c}}\right)$$

(3)

若为焊接型钢柱，其计算方法为

$$b_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{e}\mathrm{f}\mathrm{f}}=t_{\mathrm{b}\mathrm{f}}+\mathrm{2}{t}_{\mathrm{e}\mathrm{p}}+\mathrm{2}\sqrt[]{\mathrm{2}}{h}_{\mathrm{f},\mathrm{e}\mathrm{p}}+\mathrm{5}\left(t_{\mathrm{c}\mathrm{f}}+h_{\mathrm{f},\mathrm{c}\mathrm{f}}\right)$$

(4)

式（3）、（4）中：t_bf为梁翼缘厚度；t_ep为端板厚度；h_f,ep为钢梁与端板连接焊缝的有效高度；h_f,cf为柱翼缘与腹板连接焊缝的有效高度；r_c为柱翼缘根部半径。

对于有横向加劲肋的柱腹板，文献[30]给出其抗压刚度计算方法为

$$k_{\mathrm{c}\mathrm{w},\mathrm{C}}=\frac{E\left(\mathrm{0.7}{b}_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{e}\mathrm{f}\mathrm{f}}{t}_{\mathrm{c}\mathrm{w}}+t_{\mathrm{s}}{b}_{\mathrm{s}}\right)}{h_{\mathrm{c}\mathrm{w}}}$$

(5)

式中：t_s为横向加劲肋的厚度；b_s为横向加劲肋的宽度。

2.1.2　柱腹板抗拉刚度k_cw,T

与柱腹板抗压刚度的计算方法相似，对于无横向加劲肋的柱腹板，其抗拉刚度k_cw,T的计算方法为

$$k_{\mathrm{c}\mathrm{w},\mathrm{T}}=\frac{\mathrm{0.7}E{b}_{\mathrm{c}\mathrm{w},\mathrm{T},\mathrm{e}\mathrm{f}\mathrm{f}}{t}_{\mathrm{c}\mathrm{w}}}{h_{\mathrm{c}\mathrm{w}}}$$

(6)

式中：b_cw,T,eff为柱腹板受拉区的有效高度，b_cw,T,eff=h_cw-b_cw,C,eff。若为热轧型钢柱，其计算方法见式(3)；若为焊接型钢柱，其计算方法见式(4)。

对于有横向加劲肋的柱腹板，其抗拉刚度计算方法为

$$k_{\mathrm{c}\mathrm{w},\mathrm{T}}=\frac{E\left(\mathrm{0.7}{b}_{\mathrm{c}\mathrm{w},\mathrm{T},\mathrm{e}\mathrm{f}\mathrm{f}}{t}_{\mathrm{c}\mathrm{w}}+t_{\mathrm{s}}{b}_{\mathrm{s}}\right)}{h_{\mathrm{c}\mathrm{w}}}$$

(7)

2.1.3　柱腹板抗剪刚度k_cw,v

EN 1993—1‒8^[3]给出柱腹板抗剪刚度k_cw,v的计算方法为

$$k_{\mathrm{c}\mathrm{w},\mathrm{v}}=\frac{\mathrm{0.38}E{h}_{\mathrm{c}\mathrm{w}}{t}_{\mathrm{c}\mathrm{w}}}{h_{\mathrm{b}\mathrm{f}}}$$

(8)

式中：h_bf为梁上、下翼缘中心间的距离，h_bf=h_b-t_bf； h_b为梁截面高度。

2.1.4　端板抗弯刚度k_ep

采用作者前期基于连续梁模型推导的考虑螺栓抗弯刚度的T型连接初始刚度计算方法^[31,32,33]计算端板的抗弯刚度，对于螺栓无预紧力的端板连接节点，其端板抗弯刚度k_ep的计算方法为

$$k_{\mathrm{e}\mathrm{p}}=\frac{\mathrm{48}}{Z_{\mathrm{e}\mathrm{p}}\left[\mathrm{1}-{\rho }_{\mathrm{2}}\left(\mathrm{3}{\alpha }_{\mathrm{e}\mathrm{p}}-\mathrm{4}{\alpha }_{\mathrm{e}\mathrm{p}}^{\mathrm{3}}\right)-{\rho }_{\mathrm{1}}{\rho }_{\mathrm{2}}\left(\mathrm{3}-\mathrm{12}{\alpha }_{\mathrm{e}\mathrm{p}}^{\mathrm{2}}\right)\right]}$$

(9)

$$Z_{\mathrm{e}\mathrm{p}}=\frac{l_{\mathrm{e}\mathrm{p}}^{\mathrm{3}}}{E{I}_{\mathrm{e}\mathrm{p}}}=\frac{\mathrm{12}{l}_{\mathrm{e}\mathrm{p}}^{\mathrm{3}}}{E{l}_{\mathrm{e}\mathrm{f}\mathrm{f},\mathrm{e}\mathrm{p}}{t}_{\mathrm{e}\mathrm{p}}^{\mathrm{3}}}$$，$$l_{\mathrm{e}\mathrm{p}}=\mathrm{2}\left(m_{\mathrm{e}\mathrm{p}}+n_{\mathrm{e}\mathrm{p}}\right)$$，

$${\alpha }_{\mathrm{e}\mathrm{p}}=\frac{n_{\mathrm{e}\mathrm{p}}}{l_{\mathrm{e}\mathrm{p}}}$$，$${\alpha }_{\mathrm{e}\mathrm{p}\mathrm{1}}=\frac{{\alpha }_{\mathrm{e}\mathrm{p}}}{\mathrm{8}}-\frac{{\alpha }_{\mathrm{e}\mathrm{p}}^{\mathrm{3}}}{\mathrm{6}}$$，$${\alpha }_{\mathrm{e}\mathrm{p}\mathrm{2}}=\frac{{\alpha }_{\mathrm{e}\mathrm{p}}^{\mathrm{2}}}{\mathrm{2}}-\frac{\mathrm{2}{\alpha }_{\mathrm{e}\mathrm{p}}^{\mathrm{3}}}{\mathrm{3}}$$，

$${\alpha }_{\mathrm{e}\mathrm{p}\mathrm{3}}=\frac{{\alpha }_{\mathrm{e}\mathrm{p}}}{\mathrm{2}}-{\alpha }_{\mathrm{e}\mathrm{p}}^{\mathrm{2}}$$，$${\alpha }_{\mathrm{e}\mathrm{p}\mathrm{4}}=\frac{\mathrm{1}}{\mathrm{8}}-\frac{{\alpha }_{\mathrm{e}\mathrm{p}}^{\mathrm{2}}}{\mathrm{2}}$$，$${\alpha }_{\mathrm{e}\mathrm{p}\mathrm{5}}=\frac{\mathrm{1}}{\mathrm{2}}-{\alpha }_{\mathrm{e}\mathrm{p}}$$，

$${\rho }_{\mathrm{1}}=\frac{Z_{\mathrm{e}\mathrm{p}}{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{2}}+\frac{\mathrm{1}}{k_{\mathrm{b}}}-Z_{\mathrm{e}\mathrm{p}}\frac{{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{1}}{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{3}}}{{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{4}}}}{Z_{\mathrm{e}\mathrm{p}}\frac{{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{1}}{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{5}}}{{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{4}}}+\frac{l_{\mathrm{e}\mathrm{p}}^{\mathrm{2}}{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{1}}}{k_{\mathrm{b}\mathrm{b}}{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{4}}}-Z_{\mathrm{e}\mathrm{p}}{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{3}}}$$，

$${\rho }_{\mathrm{2}}=\frac{Z_{\mathrm{e}\mathrm{p}}{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{1}}}{Z_{\mathrm{e}\mathrm{p}}{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{2}}+Z_{\mathrm{e}\mathrm{p}}{\rho }_{\mathrm{1}}{\alpha }_{\mathrm{e}\mathrm{p}\mathrm{3}}+\frac{\mathrm{1}}{k_{\mathrm{b}}}}$$

式中：l_ep为与端板等效的T型连接的连续梁模型的梁长；I_ep为梁模型的截面惯性矩；l_eff,ep为梁模型塑性铰线的计算长度；m_ep为与端板等效的T型连接的螺栓轴线至翼缘根部的距离；n_ep为梁模型中用于模拟螺栓的弹簧至支座约束的距离；k_b和k_bb分别为螺栓的轴向刚度和抗弯刚度。m_ep及l_eff,ep的确定方法见EN 1993—1‒8^[3]；n_ep根据文献[32]给出的n的确定方法计算得到。

对于螺栓有预紧力的端板连接节点，需依据文献[33]求得栓板单元的割线刚度k_bs，将k_bs替代k_b代入式(9)，即可求得螺栓有预紧力的端板连接节点的端板抗弯刚度。

2.1.5　柱翼缘抗弯刚度k_cf

与端板抗弯刚度的计算方法相似，对于螺栓无预紧力的端板连接节点，其柱翼缘抗弯刚度k_cf的计算方法为

$$k_{\mathrm{c}\mathrm{f}}=\frac{\mathrm{48}}{Z_{\mathrm{c}\mathrm{f}}\left[\mathrm{1}-{\rho }_{\mathrm{2}}\left(\mathrm{3}{\alpha }_{\mathrm{c}\mathrm{f}}-\mathrm{4}{\alpha }_{\mathrm{c}\mathrm{f}}^{\mathrm{3}}\right)-{\rho }_{\mathrm{1}}{\rho }_{\mathrm{2}}\left(\mathrm{3}-\mathrm{12}{\alpha }_{\mathrm{c}\mathrm{f}}^{\mathrm{2}}\right)\right]}$$

(10)

$$Z_{\mathrm{c}\mathrm{f}}=\frac{l_{\mathrm{c}\mathrm{f}}^{\mathrm{3}}}{E{I}_{\mathrm{c}\mathrm{f}}}=\frac{\mathrm{12}{l}_{\mathrm{c}\mathrm{f}}^{\mathrm{3}}}{E{l}_{\mathrm{e}\mathrm{f}\mathrm{f},\mathrm{c}\mathrm{f}}{t}_{\mathrm{c}\mathrm{f}}^{\mathrm{3}}}$$，$$l_{\mathrm{c}\mathrm{f}}=\mathrm{2}\left(m_{\mathrm{c}\mathrm{f}}+n_{\mathrm{c}\mathrm{f}}\right)$$，

$${\alpha }_{\mathrm{c}\mathrm{f}}=\frac{n_{\mathrm{c}\mathrm{f}}}{l_{\mathrm{c}\mathrm{f}}}$$，$${\alpha }_{\mathrm{c}\mathrm{f}\mathrm{1}}=\frac{{\alpha }_{\mathrm{c}\mathrm{f}}}{\mathrm{8}}-\frac{{\alpha }_{\mathrm{c}\mathrm{f}}^{\mathrm{3}}}{\mathrm{6}}$$，$${\alpha }_{\mathrm{c}\mathrm{f}\mathrm{2}}=\frac{{\alpha }_{\mathrm{c}\mathrm{f}}^{\mathrm{2}}}{\mathrm{2}}-\frac{\mathrm{2}{\alpha }_{\mathrm{c}\mathrm{f}}^{\mathrm{3}}}{\mathrm{3}}$$，

$${\alpha }_{\mathrm{c}\mathrm{f}\mathrm{3}}=\frac{{\alpha }_{\mathrm{c}\mathrm{f}}}{\mathrm{2}}-{\alpha }_{\mathrm{c}\mathrm{f}}^{\mathrm{2}}$$，$${\alpha }_{\mathrm{c}\mathrm{f}\mathrm{4}}=\frac{\mathrm{1}}{\mathrm{8}}-\frac{{\alpha }_{\mathrm{c}\mathrm{f}}^{\mathrm{2}}}{\mathrm{2}}$$，$${\alpha }_{\mathrm{c}\mathrm{f}\mathrm{5}}=\frac{\mathrm{1}}{\mathrm{2}}-{\alpha }_{\mathrm{c}\mathrm{f}}$$，

$${\rho }_{\mathrm{1}}=\frac{Z_{\mathrm{c}\mathrm{f}}{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{2}}+\frac{\mathrm{1}}{k_{\mathrm{b}}}-Z_{\mathrm{c}\mathrm{f}}\frac{{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{1}}{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{3}}}{{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{4}}}}{Z_{\mathrm{c}\mathrm{f}}\frac{{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{1}}{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{5}}}{{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{4}}}+\frac{l_{\mathrm{c}\mathrm{f}}^{\mathrm{2}}{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{1}}}{k_{\mathrm{b}\mathrm{b}}{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{4}}}-Z_{\mathrm{c}\mathrm{f}}{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{3}}}$$，

$${\rho }_{\mathrm{2}}=\frac{Z_{\mathrm{c}\mathrm{f}}{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{1}}}{Z_{\mathrm{c}\mathrm{f}}{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{2}}+Z_{\mathrm{c}\mathrm{f}}{\rho }_{\mathrm{1}}{\alpha }_{\mathrm{c}\mathrm{f}\mathrm{3}}+\frac{\mathrm{1}}{k_{\mathrm{b}}}}$$

式中：l_cf为与柱翼缘等效的T型连接的连续梁模型的梁长；I_cf为梁模型的截面惯性矩；l_eff,cf为梁模型塑性铰线的计算长度；m_cf为与柱翼缘等效的T型连接的螺栓轴线至翼缘根部的距离；n_cf为梁模型中用于模拟螺栓的弹簧至支座约束的距离。m_cf及l_eff,cf的确定方法见EN 1993—1‒8^[3]；n_cf根据文献[33]给出的n的确定方法计算得到。

对于螺栓有预紧力的端板连接节点，可依据文献[33]求得栓板单元的割线刚度k_bs，将k_bs替代k_b并代入式(10)，即可求得螺栓有预紧力的端板连接节点的柱翼缘抗弯刚度。

2.1.6　螺栓抗拉刚度k_b

由于采用 基于连续梁模型推导的T型连接初始刚度计算方法^[31,32,33]计算端板和柱翼缘的抗弯刚度，计算得到的端板和柱翼缘抗弯刚度中已包含螺栓的抗拉刚度，因此无需单独计算螺栓的抗拉刚度。

2.2　节点初始刚度的计算方法

对于常用的梁受拉翼缘上、下各一排螺栓的端板连接节点，其各排螺栓处的等效抗拉刚度k_TS,i为

$$k_{\mathrm{T}\mathrm{S},i}=\frac{\mathrm{1}}{\frac{\mathrm{1}}{k_{\mathrm{e}\mathrm{p},i}}+\frac{\mathrm{1}}{k_{\mathrm{c}\mathrm{f},i}}+\frac{\mathrm{1}}{k_{\mathrm{c}\mathrm{w},\mathrm{T},i}}}$$

(11)

式中：k_ep,i为第i排螺栓处的端板抗弯刚度；k_cf,i为第i排螺栓处的柱翼缘抗弯刚度；k_cw,T,i为第i排螺栓处的柱腹板抗拉刚度。

在弯矩M作用下，可建立如下平衡方程：

$$k_{\mathrm{e}\mathrm{q}}{h}_{\mathrm{e}\mathrm{q}}\theta =\stackrel{\mathrm{2}}{\sum _{i=\mathrm{1}}}{k}_{\mathrm{T}\mathrm{S},i}{h}_i\theta$$

(12)

$$k_{\mathrm{e}\mathrm{q}}{h}_{\mathrm{e}\mathrm{q}}^{\mathrm{2}}\theta =\stackrel{\mathrm{2}}{\sum _{i=\mathrm{1}}}{k}_{\mathrm{T}\mathrm{S},i}{h}_i^{\mathrm{2}}\theta$$

(13)

式中：h_i为第i排螺栓的螺栓轴线至梁下翼缘中面的距离；k_eq为等效抗拉弹簧的刚度；h_eq为等效抗拉弹簧距梁下翼缘中面的距离。

联立式(12)和(13)可得

$$h_{\mathrm{e}\mathrm{q}}=\frac{\stackrel{\mathrm{2}}{\sum _{i=\mathrm{1}}}{k}_{\mathrm{T}\mathrm{S},i}{h}_i^{\mathrm{2}}}{\stackrel{\mathrm{2}}{\sum _{i=\mathrm{1}}}{k}_{\mathrm{T}\mathrm{S},i}{h}_i}$$

(14)

$$k_{\mathrm{e}\mathrm{q}}=\frac{\stackrel{\mathrm{2}}{\sum _{i=\mathrm{1}}}{k}_{\mathrm{T}\mathrm{S},i}{h}_i}{h_{\mathrm{e}\mathrm{q}}}=\frac{{\left(\stackrel{\mathrm{2}}{\sum _{i=\mathrm{1}}}{k}_{\mathrm{T}\mathrm{S},i}{h}_i\right)}^{\mathrm{2}}}{\stackrel{\mathrm{2}}{\sum _{i=\mathrm{1}}}{k}_{\mathrm{T}\mathrm{S},i}{h}_i^{\mathrm{2}}}$$

(15)

在h_eq处，由于柱腹板受拉、柱翼缘受弯、端板受弯和螺栓受拉而产生的变形δ_T为

$${\delta }_{\mathrm{T}}=\frac{M}{h_{\mathrm{e}\mathrm{q}}{k}_{\mathrm{e}\mathrm{q}}}$$

(16)

柱腹板受压产生变形δ_C为

$${\delta }_{\mathrm{C}}=\frac{M}{h_{\mathrm{e}\mathrm{q}}{k}_{\mathrm{c}\mathrm{w},\mathrm{C}}}$$

(17)

柱腹板由于剪切变形而在h_eq高度范围内产生的变形量δ_v为

$${\delta }_{\mathrm{v}}=\frac{M}{h_{\mathrm{e}\mathrm{q}}{k}_{\mathrm{c}\mathrm{w},\mathrm{v}}}$$

(18)

由式(16)~(18)可得节点的转角θ为

$$\begin{array}{c}\theta =\frac{{\delta }_{\mathrm{T}}+{\delta }_{\mathrm{C}}+{\delta }_{\mathrm{v}}}{h_{\mathrm{e}\mathrm{q}}}=\\ \frac{\mathrm{1}}{h_{\mathrm{e}\mathrm{q}}}\left(\frac{M}{h_{\mathrm{e}\mathrm{q}}{k}_{\mathrm{e}\mathrm{q}}}+\frac{M}{h_{\mathrm{e}\mathrm{q}}{k}_{\mathrm{c}\mathrm{w},\mathrm{C}}}+\frac{M}{h_{\mathrm{e}\mathrm{q}}{k}_{\mathrm{c}\mathrm{w},\mathrm{v}}}\right)\end{array}$$

(19)

综上可得，节点初始刚度k_ini为

$$k_{\mathrm{i}\mathrm{n}\mathrm{i}}=\frac{M}{\theta }=\frac{h_{\mathrm{e}\mathrm{q}}^{\mathrm{2}}}{\frac{\mathrm{1}}{k_{\mathrm{e}\mathrm{q}}}+\frac{\mathrm{1}}{k_{\mathrm{c}\mathrm{w},\mathrm{C}}}+\frac{\mathrm{1}}{k_{\mathrm{c}\mathrm{w},\mathrm{v}}}}$$

(20)

3 抗弯极限承载力计算方法

3.1　各组件抗弯极限承载力计算方法

端板连接节点的抗弯极限承载力取决于节点中最薄弱组件的极限承载力，故为计算节点的抗弯极限承载力，需计算：①柱腹板抗拉极限承载力F_cw,T,u；②柱腹板抗压极限承载力F_cw,C,u；③柱腹板抗剪极限承载力F_cw,v,u；④柱翼缘抗弯极限承载力T_cf,u；⑤端板抗弯极限承载力T_ep,u；⑥螺栓抗拉极限承载力B_u。

3.1.1　柱腹板抗压极限承载力F_cw,C,u

对于无横向加劲肋的柱腹板，其抗压极限承载力F_cw,C,u的计算方法为

$$F_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{u}}=f_{\mathrm{u}}{b}_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{e}\mathrm{f}\mathrm{f}}{t}_{\mathrm{c}\mathrm{w}}$$

(21)

对于有横向加劲肋的柱腹板，其抗压极限承载力F_cw,C,u的计算方法为

$$F_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{u}}=f_{\mathrm{u}}{b}_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{e}\mathrm{f}\mathrm{f}}{t}_{\mathrm{c}\mathrm{w}}+f_{\mathrm{u}\mathrm{s}}{t}_{\mathrm{s}}{b}_{\mathrm{s}}$$

(22)

式中：f_u为节点所用钢材的抗拉强度；f_us为加劲肋所用钢材的抗拉强度。

3.1.2　柱腹板抗拉极限承载力F_cw,T,u

对于无横向加劲肋的柱腹板，其抗拉极限承载力F_cw,T,u的计算方法为

$$F_{\mathrm{c}\mathrm{w},\mathrm{T},\mathrm{u}}=f_{\mathrm{u}}{b}_{\mathrm{c}\mathrm{w},\mathrm{T},\mathrm{e}\mathrm{f}\mathrm{f}}{t}_{\mathrm{c}\mathrm{w}}$$

(23)

对于有横向加劲肋的柱腹板，其抗拉极限承载力F_cw,T,u的计算方法为

$$F_{\mathrm{c}\mathrm{w},\mathrm{T},\mathrm{u}}=f_{\mathrm{u}}{b}_{\mathrm{c}\mathrm{w},\mathrm{T},\mathrm{e}\mathrm{f}\mathrm{f}}{t}_{\mathrm{c}\mathrm{w}}+f_{\mathrm{u}\mathrm{s}}{t}_{\mathrm{s}}{b}_{\mathrm{s}}$$

(24)

3.1.3　柱腹板抗剪极限承载力F_cw,v,u

柱腹板抗剪极限承载力F_cw,v,u的计算方法为

$$F_{\mathrm{c}\mathrm{w},\mathrm{v},\mathrm{u}}=\frac{f_{\mathrm{u}}}{\sqrt[]{\mathrm{3}}}{h}_{\mathrm{c}\mathrm{w}}{t}_{\mathrm{c}\mathrm{w}}$$

(25)

式中：h_cw为柱腹板的高度，h_cw=h_c-2t_cf；h_c为柱截面高度；t_cf为柱翼缘厚度；t_cw为柱腹板厚度。

3.1.4　端板抗弯极限承载力T_ep,u

根据EN 1993—1‒8^[3]，节点受拉区端板对应失效模式1和模式2（如图3a和3b所示）的抗弯承载力T_{ep,y,f 1}和T_{ep,y,f 2}的计算方法为

$$T_{\mathrm{e}\mathrm{p},\mathrm{y},\mathrm{f}\mathrm{1}}=\frac{\mathrm{4}{M}_{\mathrm{p},\mathrm{e}\mathrm{p}}}{m_{\mathrm{e}\mathrm{p}}}$$

(26)

$$T_{\mathrm{e}\mathrm{p},\mathrm{y},\mathrm{f}\mathrm{2}}=\frac{\mathrm{2}{M}_{\mathrm{p},\mathrm{e}\mathrm{p}}+n_{\mathrm{e}\mathrm{p},\mathrm{r}}\sum B_{\mathrm{y}}}{m_{\mathrm{e}\mathrm{p}}+n_{\mathrm{e}\mathrm{p},\mathrm{r}}}$$

(27)

$$M_{\mathrm{p},\mathrm{e}\mathrm{p}}=\frac{l_{\mathrm{e}\mathrm{f}\mathrm{f},\mathrm{e}\mathrm{p}}{t}_{\mathrm{e}\mathrm{p}}^{\mathrm{2}}{f}_{\mathrm{y}}}{\mathrm{4}}$$

(28)

$$B_{\mathrm{y}}=f_{\mathrm{b}\mathrm{y}}{A}_{\mathrm{s}}$$

(29)

式中：M_p,ep为端板的全截面塑性抗弯承载力；f_by为螺栓钢材屈服强度；A_s为螺栓的净截面面积；l_eff,ep为塑性铰线的计算长度；m_ep和n_ep,r分别为与端板等效的T型连接的螺栓轴线至翼缘根部的距离和螺栓轴线至翼缘边缘的距离。m_ep、n_ep,r及l_eff,ep的确定方法见EN 1993—1‒8^[3]。

其抗弯承载力T_ep,y取T_{ep,y,f 1}和T_{ep,y,f 2}中的较小值，即

$$T_{\mathrm{e}\mathrm{p},\mathrm{y}}=\mathrm{m}\mathrm{i}\mathrm{n}\left(T_{\mathrm{e}\mathrm{p},\mathrm{y},\mathrm{f}\mathrm{1}},T_{\mathrm{e}\mathrm{p},\mathrm{y},\mathrm{f}\mathrm{2}}\right)$$

(30)

若考虑螺栓头与端板接触面之间的均布面应力对节点承载能力的影响，则节点受拉区端板对应失效模式1的抗弯承载力为

$$T_{\mathrm{e}\mathrm{p},\mathrm{y},\mathrm{f}\mathrm{1}}=\frac{\mathrm{32}{n}_{\mathrm{e}\mathrm{p},\mathrm{r}}-\mathrm{2}{d}_{\mathrm{w}}}{\mathrm{8}{m}_{\mathrm{e}\mathrm{p}}{n}_{\mathrm{e}\mathrm{p},\mathrm{r}}-d_{\mathrm{w}}\left(m_{\mathrm{e}\mathrm{p}}+n_{\mathrm{e}\mathrm{p},\mathrm{r}}\right)}{M}_{\mathrm{p},\mathrm{e}\mathrm{p}}$$

(31)

式中：d_w为螺栓垫片的直径。

Jaspart^[34]和Piluso等^[35]基于材料的真实应力⁃应变本构关系（如图4所示）和节点受拉区端板抗弯承载力的计算方法，推导了端板对应失效模式1和失效模式2的抗弯极限承载力T_{ep,u,f 1}和T_{ep,u,f 2}的计算公式为

$$T_{\mathrm{e}\mathrm{p},\mathrm{u},\mathrm{f}\mathrm{1}}=\frac{\mathrm{4}{M}_{\mathrm{u},\mathrm{e}\mathrm{p}}}{m_{\mathrm{e}\mathrm{p}}}$$

图4 材料的真实应力⁃应变曲线

Fig.4 True stress-strain curve of material

或

$$T_{\mathrm{e}\mathrm{p},\mathrm{u},\mathrm{f}\mathrm{1}}=\frac{\mathrm{32}{n}_{\mathrm{e}\mathrm{p},\mathrm{r}}-\mathrm{2}{d}_{\mathrm{w}}}{\mathrm{8}{m}_{\mathrm{e}\mathrm{p}}{n}_{\mathrm{e}\mathrm{p},\mathrm{r}}-d_{\mathrm{w}}\left(m_{\mathrm{e}\mathrm{p}}+n_{\mathrm{e}\mathrm{p},\mathrm{r}}\right)}{M}_{\mathrm{u},\mathrm{e}\mathrm{p}}$$

(32)

$$T_{\mathrm{e}\mathrm{p},\mathrm{u},\mathrm{f}\mathrm{2}}=\frac{\mathrm{2}{M}_{\mathrm{u},\mathrm{e}\mathrm{p}}+n_{\mathrm{e}\mathrm{p},\mathrm{r}}\sum B_{\mathrm{u}}}{m_{\mathrm{e}\mathrm{p}}+n_{\mathrm{e}\mathrm{p},\mathrm{r}}}$$

(33)

$$M_{\mathrm{u},\mathrm{e}\mathrm{p}}=\frac{\xi l_{\mathrm{e}\mathrm{f}\mathrm{f},\mathrm{e}\mathrm{p}}{t}_{\mathrm{e}\mathrm{p}}^{\mathrm{2}}{f}_{\mathrm{y}}}{\mathrm{6}}$$

(34)

$$B_{\mathrm{u}}=\frac{B_{\mathrm{y}}}{\mathrm{0.9}}=\frac{f_{\mathrm{b}\mathrm{y}}{A}_{\mathrm{s}}}{\mathrm{0.9}}$$

(35)

$$\begin{array}{c}\xi =\frac{\mathrm{1}}{\mathrm{2}}\left[\mathrm{3}-{\left(\frac{{\epsilon }_{\mathrm{p}}}{{\epsilon }_{\mathrm{t}}}\right)}^{\mathrm{2}}\right]+\\ \frac{E_{\mathrm{s}}}{\mathrm{2}E}\left(\frac{{\epsilon }_{\mathrm{t}}-{\epsilon }_{\mathrm{y}}}{{\epsilon }_{\mathrm{p}}}\right)\left(\mathrm{1}-\frac{{\epsilon }_{\mathrm{y}}}{{\epsilon }_{\mathrm{t}}}\right)\left(\mathrm{2}+\frac{{\epsilon }_{\mathrm{y}}}{{\epsilon }_{\mathrm{t}}}\right)-\\ \frac{E_{\mathrm{s}}-E_{\mathrm{t}}}{\mathrm{2}E}\left(\frac{{\epsilon }_{\mathrm{t}}-{\epsilon }_{\mathrm{s}}}{{\epsilon }_{\mathrm{p}}}\right)\left(\mathrm{1}-\frac{{\epsilon }_{\mathrm{s}}}{{\epsilon }_{\mathrm{t}}}\right)\left(\mathrm{2}+\frac{{\epsilon }_{\mathrm{s}}}{{\epsilon }_{\mathrm{t}}}\right)\end{array}$$

(36)

其抗弯极限承载力T_ep,u取T_{ep,u,f 1}和T_{ep,u,f 2}中的较小值，即

$$T_{\mathrm{e}\mathrm{p},\mathrm{u}}=\mathrm{m}\mathrm{i}\mathrm{n}\left(T_{\mathrm{e}\mathrm{p},\mathrm{u},\mathrm{f}\mathrm{1}},T_{\mathrm{e}\mathrm{p},\mathrm{u},\mathrm{f}\mathrm{2}}\right)$$

(37)

3.1.5　柱翼缘抗弯极限承载力T_cf,u

与端板抗弯极限承载力的计算方法相似，其抗弯极限承载力的计算方法为

$$T_{\mathrm{c}\mathrm{f},\mathrm{u},\mathrm{f}\mathrm{1}}=\frac{\mathrm{4}{M}_{\mathrm{u},\mathrm{c}\mathrm{f}}}{m_{\mathrm{c}\mathrm{f}}}$$

或

$$T_{\mathrm{c}\mathrm{f},\mathrm{u},\mathrm{f}\mathrm{1}}=\frac{\mathrm{32}{n}_{\mathrm{c}\mathrm{f},\mathrm{r}}-\mathrm{2}{d}_{\mathrm{w}}}{\mathrm{8}{m}_{\mathrm{c}\mathrm{f}}{n}_{\mathrm{c}\mathrm{f},\mathrm{r}}-d_{\mathrm{w}}\left(m_{\mathrm{c}\mathrm{f}}+n_{\mathrm{c}\mathrm{f},\mathrm{r}}\right)}{M}_{\mathrm{u},\mathrm{c}\mathrm{f}}$$

(38)

$$T_{\mathrm{c}\mathrm{f},\mathrm{u},\mathrm{f}\mathrm{2}}=\frac{\mathrm{2}{M}_{\mathrm{u},\mathrm{c}\mathrm{f}}+n_{\mathrm{c}\mathrm{f},\mathrm{r}}\sum B_{\mathrm{u}}}{m_{\mathrm{c}\mathrm{f}}+n_{\mathrm{c}\mathrm{f},\mathrm{r}}}$$

(39)

$$M_{\mathrm{u},\mathrm{c}\mathrm{f}}=\frac{{\xi }_{\mathrm{e}\mathrm{f}\mathrm{f},\mathrm{c}\mathrm{f}}{t}_{\mathrm{c}\mathrm{f}}^{\mathrm{2}}{f}_{\mathrm{y}}}{\mathrm{6}}$$

(40)

$$T_{\mathrm{c}\mathrm{f},\mathrm{u}}=\mathrm{m}\mathrm{i}\mathrm{n}\left(T_{\mathrm{c}\mathrm{f},\mathrm{u},\mathrm{f}\mathrm{1}},T_{\mathrm{c}\mathrm{f},\mathrm{u},\mathrm{f}\mathrm{2}}\right)$$

(41)

式中：l_eff,cf为塑性铰线的计算长度；m_cf和n_cf,r分别为与柱翼缘等效的T型连接的螺栓轴线至翼缘根部的距离和螺栓轴线至翼缘边缘的距离。m_cf、n_cf,r及l_eff,cf的确定方法详见EN 1993—1‒8^[3]。

3.1.6　螺栓抗拉极限承载力B_u

螺栓的抗拉极限承载力计算方法见式(35)。

3.2　节点抗弯极限承载力的计算方法

对于常用的梁受拉翼缘上、下各一排螺栓的端板连接节点，各排螺栓处的抗拉极限承载力F_u,T,i为

$$F_{\mathrm{u},\mathrm{T},i}=\mathrm{m}\mathrm{i}\mathrm{n}\left(T_{\mathrm{e}\mathrm{p},\mathrm{u},i},T_{\mathrm{c}\mathrm{f},\mathrm{u},i},\mathrm{2}{B}_{\mathrm{u},i}\right)$$

(42)

式中：T_ep,u,i、T_cf,u,i和2B_u,i分别为第i排螺栓处的端板抗弯极限承载力、柱翼缘抗弯极限承载力和螺栓抗拉极限承载力。

则节点受拉区的极限承载力F_u,T为

$$F_{\mathrm{u},\mathrm{T}}=\mathrm{m}\mathrm{i}\mathrm{n}\left(F_{\mathrm{c}\mathrm{w},\mathrm{T},\mathrm{u}},\stackrel{\mathrm{2}}{\sum _{i=\mathrm{1}}}{F}_{\mathrm{u},\mathrm{T},i}\right)$$

(43)

节点受压区的极限承载力F_u,C为

$$F_{\mathrm{u},\mathrm{C}}=F_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{u}}$$

(44)

节点受剪区的极限承载力F_u,v为

$$F_{\mathrm{u},\mathrm{v}}=F_{\mathrm{c}\mathrm{w},\mathrm{v},\mathrm{u}}$$

(45)

若节点的失效模式为节点受拉区破坏，则节点的抗弯极限承载力M_max,T可由下式计算：

$$M_{\mathrm{m}\mathrm{a}\mathrm{x},\mathrm{T}}=\mathrm{m}\mathrm{i}\mathrm{n}\left(F_{\mathrm{c}\mathrm{w},\mathrm{T},\mathrm{u}}{h}_{\mathrm{b}\mathrm{f}},\stackrel{\mathrm{2}}{\sum _{i=\mathrm{1}}}{F}_{\mathrm{u},\mathrm{T},i}{h}_i\right)$$

(46)

若节点的失效模式为节点受压区破坏，则节点的抗弯极限承载力M_max,C可由下式计算：

$$M_{\mathrm{m}\mathrm{a}\mathrm{x},\mathrm{C}}=F_{\mathrm{c}\mathrm{w},\mathrm{C},\mathrm{u}}{h}_{\mathrm{b}\mathrm{f}}$$

(47)

若节点的失效模式为节点受剪区破坏，则节点的抗弯极限承载力M_max,v可由下式计算：

$$M_{\mathrm{m}\mathrm{a}\mathrm{x},\mathrm{v}}=F_{\mathrm{c}\mathrm{w},\mathrm{v},\mathrm{u}}{h}_{\mathrm{b}\mathrm{f}}$$

(48)

节点抗弯极限承载力M_max取M_max,T、M_max,C和M_max,v中的最小值，即

$$M_{\mathrm{m}\mathrm{a}\mathrm{x}}=\mathrm{m}\mathrm{i}\mathrm{n}\left(M_{\mathrm{m}\mathrm{a}\mathrm{x},\mathrm{T}},M_{\mathrm{m}\mathrm{a}\mathrm{x},\mathrm{C}},M_{\mathrm{m}\mathrm{a}\mathrm{x},\mathrm{v}}\right)$$

(49)

4 屈服后刚度的计算方法

对于端板连接节点屈服后刚度k_p，王卫永等^[25]假定节点在受力变形过程中不发生应变强化，即假定k_p=0。然而，从节点的M⁃θ曲线（图1）可看出，曲线存在明显的应变强化段；Yee等^[18]也明确指出，当节点的失效模式为柱腹板受压区破坏时，节点往往具有明显的后屈曲强度，即节点会发生应变强化，故假定k_p=0不甚合理。Yee等^[18]给出了节点屈服后刚度k_p的计算方法，但该计算方法较为复杂。井泉^[36]采用Yee等^[18]提出的四参数指数模型（如式(1)所示）对普通钢端板连接节点的M⁃θ曲线进行拟合，结果表明，k_p取初始刚度k_ini的1/50较为合理。因此，本文采用文献[36]的建议，取k_p=0.02k_ini。

由文献[21,32]可知，形状参数c对式(1)计算结果的影响较小，可忽略，故取c=0。

5 弯矩⁃转角曲线截断方法

式(1)描述的是节点承受的弯矩M与节点转角θ间的关系，θ无限大时，式(1)计算得到的M也会无限大，从而使M⁃θ曲线的预测结果无限长。然而，节点的承载能力和转动能力都是有限的。为准确模拟节点的实际力学行为，有必要在预测曲线上确定截断点L，预测曲线在该点截断，如图5所示。

图5 M⁃θ曲线的截断

Fig.5 Truncation of M⁃θ curve

从节点承载能力的角度出发，截断点的纵坐标M_lim不应大于按第3节计算方法得到的节点抗弯极限承载力M_max；现有研究认为，θ达到0.05 rad可视为节点拥有足够的转动能力^[25,37]，所以从变形的角度出发，截断点的横坐标θ_lim不应超过0.05 rad。

因此，本文提出预测曲线的截断方法如下：① 按第3节计算方法得到节点的抗弯极限承载力M_max；② 将θ=0.05 rad代入式(1)，求得节点弯矩M_0.05，若M_0.05≤M_max，则θ_lim=0.05 rad，M_lim=M_0.05；若M_0.05>M_max，则M_lim=M_max，将M_max代入式(1)，即可求得θ_lim；③ 由①和②确定了截断点L（θ_lim，M_lim），预测曲线在该点截断。

6 模型准确性验证

为验证所提出的k_ini、M_max、k_p和c计算方法的准确性，以及采用式(1)预测高强钢端板连接节点在常温、火灾下、火灾后M⁃θ曲线的合理性，与3组常温下^[26,27]、 3组火灾下^[30,31]及 4组火灾后^[26]的Q690和Q960高强钢端板连接节点足尺试验结果进行了对比。

6.1　试验的k_ini、M_max和k_p

采用提出的初始刚度k_ini和抗弯极限承载力M_max计算方法计算文献[26,27,28,29]各试验试件，得到计算值k_ini,cal和M_max,cal，如表1所示。表1同时列出节点初始刚度k_ini和抗弯极限承载力M_max的试验值k_ini,ep和M_max,ep，以及采用EN 1993—1‒8^[3]计算得到的节点初始刚度k_ini,EN。

表1 文献[26,27,28,29]各试验试件的k_ini和M_max试验和计算结果

Tab.1 Test and calculation results of kini and Mmax of various test specimens in Refs. [26-29]

试验状态	编号	Mmax/(kN ∙ m)			kini/(N ∙ mm-1)
		Mmax,ep	Mmax,cal	Mmax,ep/Mmax,cal	kini,ep	kini,EN	kini,cal	kini,ep/kini,EN	kini,ep/kini,cal
常温下	Q690A1	406.32	346.79	1.17	20 439	34 271	27 608	0.60	0.74
	Q690A2	394.98	346.79	1.14	22 142	34 271	27 608	0.65	0.80
	Q960A1	332.70	329.22	1.01	15 201	24 526	19 480	0.62	0.78
火灾下	Q690F1	207.39	161.63	1.28	12 121	17 733	13 463	0.68	0.90
	Q960F1	148.73	147.18	1.01	10 344	12 414	9 310	0.83	1.11
	Q960F2	152.44	147.18	1.04	9 375	12 414	9 310	0.76	1.01
火灾后	Q690P1	385.10	310.96	1.24	17 484	34 253	27 056	0.51	0.65
	Q690P2	387.56	310.96	1.25	18 632	34 253	27 056	0.54	0.69
	Q960P1	321.79	307.35	1.05	12 695	24 135	18 820	0.53	0.67
	Q960P2	311.46	307.35	1.01	13 209	24 135	18 820	0.55	0.70

从表1可看出，Q960高强钢端板连接节点的M_max,cal和M_max,ep较接近，相对误差在5%以内；Q690高强钢端板连接节点的M_max,ep大约是M_max,cal的1.2倍，这与ξ的计算结果有关，即与钢材本构模型的选取有关；除了Q960高强钢在火灾下的情况，其余节点的k_ini,cal高于试验值k_ini,ep，但都比k_ini,EN更接近试验值k_ini,ep。因此，提出的节点初始刚度k_ini和抗弯极限承载力M_max的计算方法能够较为准确地预测高强钢端板连接节点在常温、火灾下及火灾后的力学行为。同时由第4节可知，各试验试件的k_p,cal=0.02k_ini,cal，c_cal=0。

6.2　M⁃θ曲线的对比

将第6.1节得到的作者前期足尺试验研究中^[26,27,28,29]各试验试件的初始刚度计算值k_ini,cal、抗弯极限承载力计算值M_max,cal、屈服后刚度计算值k_p,cal和形状系数计算值c_cal代入式(1)，得到M与θ的表达式。然后，采用第5节的曲线截断方法确定截断点L，并在该点截断曲线，从而绘出各试验试件M⁃θ曲线的预测结果，其与试验结果的对比如图6所示。

图6 M⁃θ曲线预测模型与试验结果的对比

Fig.6 Comparisons of M⁃θ curves between prediction models and test results

由图6可知，高强钢端板连接节点的数学模型和试验结果较吻合；火灾后Q690高强钢端板连接节点的M_max,cal略低于M_max,ep，即数学模型预测的M⁃θ曲线峰值低于试验结果，因此模型对火灾后Q690的预测结果偏于安全，这可能与Q690钢材材性的离散性和Q690钢材火灾后材性试验（本文模型的材性输入数据）的准确性有关。综上所述，本文提出的节点屈服后刚度k_p和形状系数c的计算方法以及采用式(1)预测高强钢端板连接节点在常温、火灾下及火灾后的M⁃θ曲线是合理的。

7 结论

本文提出了高强钢端板连接节点在常温、火灾下及火灾后的初始刚度k_ini、屈服后刚度k_p和抗弯极限承载力M_max的计算方法，将k_ini、k_p、M_max和c代入Yee等提出的四参数指数模型，得到节点M⁃θ曲线的预测结果，将初始刚度k_ini和抗弯极限承载力M_max的计算结果以及节点M⁃θ曲线的预测结果，与作者前期的Q690与Q960节点足尺试验结果进行对比，研究结果表明：

（1） M_max的计算结果与试验值较吻合，其中Q960的相对误差在5%以内；k_ini的计算结果比EN 1993—1‒8的结果更接近试验值。所提出的节点初始刚度k_ini和抗弯极限承载力M_max的计算方法能够较为准确地预测Q690与Q960高强钢端板连接节点在常温、火灾下及火灾后的力学行为。

（2） 对于Q960高强钢端板连接节点，预测曲线与试验曲线较吻合；对于Q690高强钢端板连接节点，预测曲线的峰值略低于试验曲线的峰值，该预测结果偏于安全。由此可见，本文提出的屈服后刚度k_p的计算方法以及采用Yee等提出的四参数指数模型预测高强钢端板连接节点在常温、火灾下及火灾后的M⁃θ曲线是合理的。

本文研究成果可为Q690与Q960高强钢端板连接节点在常温下、火灾下及火灾后力学性能的预测与评定提供参考，为Q690与Q960高强钢梁柱节点的抗火设计、优化设计及火灾后评估，乃至高强钢在土木工程领域的合理应用提供理论依据。