系列文章的第1篇[1]已经对高强度Q460钢梁高温下的抗火性能进行了理论分析,采用有限差分法给出了高强度Q460钢梁的温度分布函数表达式,并结合Q460钢材的高温力学性能参数得到了温度均匀分布和不均匀分布下的Q460钢梁极限承载力、临界温度和稳定系数。本文将对理论结果进行有限元验证,并将普通Q235钢梁和高强度Q460钢梁的抗火性能进行对比,得到他们的区别。最后提出了温度不均匀分布的Q460钢梁简化设计方法。
系列文章的第1篇[1]给出了钢梁各个组件温度的计算方法和临界弯矩的计算理论,为了验证其正确性,采用有限元软件ANSYS对高强度Q460钢梁进行了有限元分析,包括热分析和结构分析两部分。在传热分析中,采用有限元分析的结果验证了有限差分法得到的钢梁各组件的温度分布,并对分析结果产生偏差的原因进行了分析;在结构分析中,采用有限元方法计算了高强度Q460钢梁的临界弯矩,并将有限元分析的结果与等效刚度法计算的结果进行了对比,据此判断该理论的可靠性。
钢梁的截面尺寸为HN400×200×8×13,温度计算考虑三面受火的情况,求解该钢梁的温度场分布。结构分析中,将热分析得到的结果作为体荷载施加到结构上。热分析单元类型选取SOLID70,该单元可用于三维稳态或瞬态热分析,每个单元有8个节点,每个节点仅有温度自由度。在对钢梁进行网格划分时,为了适应几何形状比较复杂的实体模型,该单元可以自动蜕化为四面体与棱柱体。
进行有限元热分析时,钢梁周围的空气温度按照ISO 834标准升温曲线确定,初始温度设定为20 ℃,钢梁的传热边界条件根据实际情况确定,即上翼缘上边界为绝热状态,其他受火的边界设置对流和辐射参数。分析中假定结构产生的变形对受热边界条件没有影响。有限元热分析中用到的钢梁的热工性能参数,例如:比热、导热系数、热膨胀系数和密度等,按照EC3[2]中的规定取值。
有限元分析的结果和有限差分法计算结果的对比见图 1。从图 1中可以看出,二者吻合较好,验证了采用有限差分法计算Q460钢梁各个部件温度的正确性。
为了进一步对热分析模型进行验证,将文献[3]抗火试验中的温度测量结果和该模型分析的结果进行了对比。试验中在3个不同的截面(1/4、1/2和3/4高度)测试了约束Q460钢柱的温度分布,约束钢柱截面尺寸为H200×195×8×8,空气升温为试验中实测的温度曲线,分析结果和试验结果的对比见图 2。从图 2中可以看出,二者吻合较好,最大误差控制在50 ℃以内。
钢梁承受横向荷载作用时,当荷载不大时,一般情况下可保持在平衡状态。随着荷载的增大并达到某一临界值时,如果此时在钢梁侧向施加一个微小的作用力,钢梁的平衡位置被打破,产生较大弯曲变形并伴随着扭转,这种现象称为梁的弯扭屈曲,该弯矩临界值定义为临界弯矩。采用有限元软件ANSYS对高强度Q460钢梁的整体稳定承载力进行计算。对构件进行分析时,采用的单元类型为SHELL181单元,该单元具有4个节点,每个节点有4个自由度,可以分析材料的塑性、大变形和大应变。在结构分析中,需要设置材料参数,主要有钢材高温屈服强度、高温弹性模量、泊松比、应力-应变关系等。其中高温下的强度和弹性模量取值按照系列文章第1篇中的式(1)和式(2)计算得到。应力应变关系采用双线性弹塑性模型。钢梁的有限元模型和网格划分见图 3(a)。
需要指出的是,当利用壳单元对H型钢梁进行有限元分析时,翼缘或腹板板件在钢梁发生整体屈曲之前可能会出现局部屈曲和畸变屈曲,这样就不能得到钢梁的弯扭屈曲临界弯矩。为了避免局部屈曲和畸变屈曲的出现,在有限元模型中,可以通过设置“特殊加劲肋”的方法来实现。一般来讲,沿梁的跨度方向设置5~7道特殊加劲肋即可避免出现局部屈曲和畸变屈曲。特殊加劲肋的设置方法是让加劲肋与钢梁的腹板相连,而与翼缘的连接仅耦合加劲肋平面内的位移,通过这种方法可以保证加劲肋支撑翼缘和腹板不发生屈曲,同时也不增加钢梁的抗扭刚度和抗弯刚度。钢梁端部的约束条件为夹支铰接,即约束端部截面内节点在平面内的位移和一个节点的纵向位移,加劲肋设置和约束条件见图 3(b)。
采用有限元分析钢梁的临界弯矩时,分两个步骤。第1步进行钢梁的静态分析,分析时需要设置激活预应力效应和设定两个荷载步,第1个荷载步中设置初始温度和参考温度,设定为20 ℃,并将前面分析得到的钢梁温度分布数据作为体荷载施加到钢梁上。第2个荷载步施加作用在钢梁上的横向荷载,考虑均布荷载和集中荷载两种情况,分析求解直至结束。第2步进行结构的特征值屈曲分析,首先设置屈曲模型和特征值提取方式,然后求解。求解完成后,可以在后处理中画出钢梁发生弯扭屈曲的屈曲模态。临界弯矩的数值即为施加的荷载和屈曲荷载系数的乘积。图 3(c)为钢梁的弯扭屈曲模态。
根据可以查到的文献,尚未发现有Q460钢梁的抗火性能试验。为了验证有限元模型的正确性,采用普通钢梁的抗火性能试验[4]对有限元模型进行验证。试验中,钢梁截面尺寸为H200×150×6×9,钢材为Q235B,测量的屈服强度为330 MPa,极限强度415 MPa。荷载为5分点加载,荷载大小为17.5 kN。构件按标准升温曲线ISO834进行升温,材料高温下的屈服强度和应力应变关系按照EC3[2]中的值给出,有限元分析得到的跨中挠度和试验结果的对比见图 4,从图 4可以看出,二者总体吻合较好,得到的临界温度相差很小,但升温过程梁的变形和达到临界温度后试件的变形有些差异,这主要是因为试验中随着挠度的增大,荷载很难保持恒定,通过手动调节油泵控制千斤顶荷载,其次是试验测量的试件温度分布和有限元分析得到的温度分布有差异。为了验证系列文章的第1篇中温度修正公式与等效刚度法用于分析温度不均匀分布的钢梁受力性能的可靠性,对截面为HN400×200×8×13,跨度分别为4、6、8 m的高强度度Q460钢梁进行了分析,该梁无侧向支撑,三面受火,考虑了温度的均匀分布和不均匀分布。采用有限元分析得到的温度分布作为温度荷载值施加到结构模型上,将ANSYS得到的结果与理论计算结果进行对比,如图 5所示,图 5中横坐标温度值,对于温度均匀分布的钢梁为截面平均温度,对于温度不均匀分布的钢梁,该数值是腹板的温度值。从图中可见,整体吻合较好,通过等效刚度法计算得到的高强度Q460钢梁的高温下承载力,其数值略大于ANSYS分析得到的结果。原因可能是由于温度修正公式具有一定误差,以及采用有限元模型与理想铰接模型具有一定差异等。总体来看,理论分析结果与ANSYS分析结果吻合较好。
为了验证系列文章的第1篇提出的等效刚度法计算温度不均匀分布的钢梁的临界弯矩的可靠性,考虑钢梁截面温度的不均匀分布,并选取腹板温度为300和600 ℃两个温度值以及4、6、8 m 3个跨度,分别采用等效刚度法和有限元法对临界弯矩进行计算,计算结果对比见表 1。从表 1可以看出,ANSYS分析计算结果与等效刚度法计算结果相差较小,不超出2.5%。
此外,还分别采用等效刚度法和有限元分析法对钢梁的临界温度结果进行了计算和对比,计算时考虑了钢梁截面温度的不均匀分布,选择两个荷载比R=0.7和0.9以及4、6、8 m 3个跨度,两种方法计算结果的对比见表 2。从表中可以看出,等效刚度法和ANSYS有限元分析结果吻合较好,二者差异均在6%以内。
从有限元分析过程可以看出,对考虑不均匀温度分布的高强度Q460钢梁进行承载力分析时,相比理论计算方法,有限元分析过程较复杂,理论方法更便于进行数据计算和分析,从而得到简化的设计方法。采用上文的理论分析方法对Q460钢梁是否考虑温度的不均匀分布的承载力和临界温度以及和普通Q235钢梁的稳定承载力和临界温度进行对比,从而得到高强度Q460钢梁和普通钢梁抗火性能的差异。
以H型高强度Q460钢梁和普通Q235钢梁为例,截面尺寸均为HN400×200×8×13,无侧向支撑;截面面积为85.76 cm2,绕强轴的惯性矩I=24 300 cm4,截面模量W=1 215cm3,考虑钢梁的三面受火,即上翼缘上表面绝热,假设常温下钢梁的整体稳定系数为φ′b,钢梁上的荷载比为R,考虑两种荷载类型:梁上翼缘作用有沿弱轴方向的均布荷载Q和集中荷载P。
系列文章第1篇中式(10)的温度修正根据《建筑钢结构防火技术规范》(CECS 200:2006)[5]的式(4.10)得到,为了便于对比,在对比温度均匀分布的钢梁和温度不均匀分布的钢梁时,钢梁的腹板温度保持相同,当考虑上翼缘和下翼缘的温度与腹板温度不同时,才能得到钢梁温度的不均匀分布对其稳定承载力的影响。所以算例中涉及的温度数值,对温度均匀分布的钢梁指的是截面温度,也是平均温度,对温度不均匀分布的钢梁指的是钢梁的腹板温度。
钢梁跨度选择了4、6、8 m共3种,荷载类型为均布荷载和集中荷载,在是否考虑温度不均匀分布的前提下,计算每种跨度每个作用荷载类型钢梁的临界弯矩随腹板温度的变化趋势。根据系列文章第1篇的计算方法,即采用等效刚度法得到不同温度不同跨度的两种荷载类型的临界弯矩,如图 6所示。由图 6可以看出,对于温度均匀分布和不均匀分布的高强度Q460钢梁,临界弯矩的大小都随温度的升高而逐渐降低。这是由于钢梁的弯扭屈曲稳定承载力主要由钢材的强度和刚度决定,而在高温下,二者皆受温度的影响,随着温度的升高而降低。从图 6还可以看出,当腹板温度相同时,对于温度均匀分布的高强Q460钢梁和温度不均匀分布的钢梁,临界弯矩随温度升高的变化具有差异,表现在温度不均匀分布时,腹板温度相同的条件下,临界弯矩高于温度均匀分布的钢梁。并且,腹板温度越高时,这种差异越明显。当腹板温度达到650 ℃左右时,差值达到最大值,随着腹板温度的继续增加,二者的差异开始降低。发生这种现象的原因是高温下截面的温度差不同,在温度较低时温差小,从而钢梁的临界弯矩差异较小。但随着温度的升高,对于温度不均匀分布的钢梁,上下翼缘的温度差增大,临界弯矩提高效应开始变得明显,造成与温度均匀分布的钢梁相比临界弯矩差值又逐渐增大。而当温度较高时,虽然温差仍然存在,但钢梁的材料力学性能下降严重,造成临界弯矩的差异又逐渐降低。
为了对比高强度Q460钢梁与普通Q235钢梁抗火性能的差异,采用系列文章的第1篇中提出的等效刚度法,分别对高强Q460钢梁和普通Q235钢梁在高温下的临界弯矩进行了计算,结果对比见图 7。从图 7可以看出,对于温度不均匀分布的钢梁,Q460钢梁比Q235钢梁在温度分布相同的条件下,前者的临界弯矩大于后者,在400 ℃之前,差异不是很明显,随着温度的升高,两者的临界弯矩都急剧降低,Q460钢梁下降的幅度明显低于Q235钢梁,温度越高二者的差异越明显。原因是高强度Q460钢材的高温力学性能优于普通钢,尤其是弹性模量数值,在同样的温度下,比普通钢高很多[6]。对不同跨度的钢梁,考虑温度的不均匀分布时,随着跨度的增大,腹板温度相同的高强度Q460钢梁或普通Q235钢梁,临界弯矩均降低,这是因为临界弯矩的大小和自由支撑的长度有关,该长度越大,钢梁越易于发生侧向弯扭屈曲。从图 7中还可以看出,承受集中荷载作用钢梁的临界弯矩要大于承受均布荷载作用的钢梁,这是因为钢梁的临界弯矩除了和最大弯矩有关,还与弯矩的分布有关,均布荷载下比集中荷载下分布更加饱满。
依据系列文章的第1篇中临界温度的计算方法得到不同稳定系数下高强Q460钢梁的临界温度随荷载比的变化情况,如图 8(a)所示。从该图中可以看出,对于相同的稳定系数,随温度的升高,Q460钢梁的临界温度逐渐降低。荷载比越大,临界温度降低得越多。在相同的荷载比下,随着稳定系数的增加,临界温度逐渐升高。稳定系数对临界温度的影响仅当荷载比较大时(R>0.9),才比较明显,当荷载比达到0.95时,稳定系数为0.5和1.0的钢梁相比,临界温度的差值已达120 ℃。
为了得到温度不均匀分布对高强度Q460钢梁临界温度的影响,考虑两个稳定系数:0.6和0.9,采用系列文章第1篇中给出的临界温度计算方法,计算高强度Q460钢梁的临界温度随荷载比变化的规律,如图 8(b)所示。从该图中可以看出,对于温度不均匀分布和温度均匀分布的Q460钢梁,随着钢梁上荷载比的增大,其临界温度均不同程度下降。对温度均匀分布的Q460钢梁,临界温度的下降较为平缓;但是对于温度不均匀分布的Q460钢梁,临界温度下降则相对较快。对于相同的稳定系数,温度不均匀分布的Q460钢梁临界温度要高于温度均匀分布的临界温度,并且临界温度的差值随荷载比的升高逐渐降低。主要是因为考虑钢梁的温度不均匀分布之后,当钢梁的腹板温度保持相同的条件下,钢梁的上下翼缘温温度较低,尤其是对于上翼缘,温度更低,从而造成钢梁上、下翼缘的弹性模量稍大,因而临界弯矩较高,相应的临界温度就稍高。在相同的荷载比下,温度均匀分布的Q460钢梁的临界温度随稳定系数的增加而降低;而温度不均匀分布的钢Q460梁的临界温度随稳定系数的提高而略微增加。
钢材种类可能会对钢梁的抗火性能产生影响,为了对比高强度Q460钢梁和普通Q235钢梁抗火性能的差异,根据提出的等效刚度法对两种钢梁的抗火性能进行了对比。对比结果见图 8(c)。从图中可以看出,对两种钢梁,稳定系数对临界温度的影响不大,普通Q235钢梁的临界温度随荷载比的提高而降低,荷载比小于0.8时近似线性相关。而对于Q460钢梁,临界温度随着荷载比的变化表现出折线性变化,而且其临界温度数值明显高于同等条件下的普通Q235钢梁。因而,高强度Q460钢梁的抗火性能与普通钢梁相比具有较大的差异。
为了简化计算过程并方便设计人员使用,根据前面推导的原理和方法,采用等效刚度法对多个不同截面的钢梁,考虑温度的不均匀分布后,计算了其稳定承载力和临界温度,计算结果以表格的形式给出,从而可以直接查找使用。
由系列文章的第1篇的理论推导易知,考虑温度不均匀分布的高强度Q460钢梁的整体稳定系数φbTeq′受钢梁上作用的荷载类型、钢梁的横截面尺寸等参数的影响。表 3给出了按Q460高强钢梁的自由长度、荷载情况和截面类型计算得到的φbTeq′,以直接查用,由于篇幅有限,仅给出了部分数据,对于其他荷载情况和截面类型以及其他温度下的整体稳定系数,可见参考文献[7]。
对温度不均匀分布的高强度Q460钢梁,采用数值计算方法,并考虑了荷载种类(均布荷载或集中荷载)与荷载作用点位置(上翼缘或下翼缘),计算得到多个截面尺寸的H型钢在多个荷载比及自由长度下的临界温度Tdeq,见表 4。对于表中未给出的其他截面和荷载比下的临界温度,可见参考文献[7]。
为了清楚的说明该理论的使用过程,设计了一个算例,对其进行抗火性能验算。该算例中,钢梁为Q460焊接H形简支梁,截面尺寸为H350×175×7×11,跨度为5 m,没有侧向支撑,按照三面受火确定钢梁各个板件的温度,并考虑温度不均匀分布的影响,已知常温下该梁的整体稳定系数为φ′b=0.401,梁上作用有均布荷载q=25 kN/m,作用位置在上翼缘;耐火极限为20 min,验算该梁是否满足抗火要求。
1) 计算钢梁的临界温度Tdeq
钢梁的荷载比:$R=\frac{M}{{{\varphi _{\rm{b}}}'Wf}}=\frac{1}{8} \cdot \frac{{q{l^2}{\gamma _{\rm{R}}}}}{{{\varphi _{\rm{b}}}'W{f_{\rm{y}}}}}=\frac{1}{8} \cdot \frac{{25 \times {5^2} \times {{10}^6} \times 1.1}}{{0.401 \times 782 \times {{10}^3} \times 503}}=0.545$
则根据表 4可得构件的临界温度:Tdeq=787℃
2) 计算钢梁受火20 min时刻的腹板温度Twb
三面受火下截面形状系数:F/V=190.23 m-1
由系列文章第1篇中式(3),得到钢梁受火20 min时刻的温度:Ts=721.9 ℃
再由系列文章第1篇中式(10)得此时刻腹板温度:Twb=742.7 ℃
故有Twb < Tdeq
因此,钢梁满足抗火要求。
1) 计算钢梁受火20 min时刻的腹板温度Twb
由4.1节中的计算结果可知Twb=742.7℃
2) 计算钢梁在腹板温度Twb时的稳定承载力McrTeq
根据表 3线性插值可得φbTeq′=0.838;
由系列文章第1篇中式(1)可计算腹板温度为Twb时Q460钢材的强度折减系数为${\eta _T} = \frac{{{f_{\rm{y}}}_T}}{{{f_{\rm{y}}}}} = 0.308$;
则钢梁腹板温度为Twb时的稳定承载力为:$M_{{\rm{cr}}T}^{{\rm{eq}}} = \varphi _{{\rm{b}}T}^{{\rm{eq}}'}{W^{{\rm{eq}}}}{f_{{\rm{y}}T}} \approx \varphi _{{\rm{b}}T}^{{\rm{eq}}'}W{\eta _T}{f_{\rm{y}}} = 0.838 \times 782 \times {10^3} \times 0.308 \times 503 \times {10^{ - 6}} = 101.52\left( {{\rm{kN}} \cdot {\rm{m}}} \right)$;
钢梁所承受的弯矩M为M=ql2/8=25×52/8=78.13(kN·m);
故有McrTeq>M,因此,钢梁满足抗火要求。
采用有限元分析对高强度Q460钢梁的温度分布和极限承载力进行了计算,将计算结果和系列文章第1篇的理论分析的结果进行对比,验证了理论分析的正确性。采用理论计算方法对高强度Q460钢梁和普通Q235钢梁的抗火性能进行了对比,得到以下结论:
1) 同一荷载下,温度均匀分布的Q460钢梁的临界温度随稳定系数的增加而降低;而温度不均匀分布的Q460钢梁的临界温度随稳定系数的增加而略微上升。
2) 在同样的稳定系数和荷载比下,高强度Q460钢梁的临界温度结果与普通Q235结构钢梁相差较大。