• ISSN 1673-5722
  • CN 11-5429/P

基于小比例缩尺模型结构试验的小型钢架反力墙优化设计与力学性能分析

孙柏涛 姜琪 闫培雷

孙柏涛, 姜琪, 闫培雷. 基于小比例缩尺模型结构试验的小型钢架反力墙优化设计与力学性能分析[J]. 震灾防御技术, 2018, 13(4): 869-877. doi: 10.11899/zzfy20180414
引用本文: 孙柏涛, 姜琪, 闫培雷. 基于小比例缩尺模型结构试验的小型钢架反力墙优化设计与力学性能分析[J]. 震灾防御技术, 2018, 13(4): 869-877. doi: 10.11899/zzfy20180414
Sun Baitao, Jiang Qi, Yan Peilei. Optimization Design and Mechanical Property Analysis of Small Steel Frame Reaction-wall Based on Large Scale Model Structure Test[J]. Technology for Earthquake Disaster Prevention, 2018, 13(4): 869-877. doi: 10.11899/zzfy20180414
Citation: Sun Baitao, Jiang Qi, Yan Peilei. Optimization Design and Mechanical Property Analysis of Small Steel Frame Reaction-wall Based on Large Scale Model Structure Test[J]. Technology for Earthquake Disaster Prevention, 2018, 13(4): 869-877. doi: 10.11899/zzfy20180414

基于小比例缩尺模型结构试验的小型钢架反力墙优化设计与力学性能分析

doi: 10.11899/zzfy20180414
基金项目: 

黑龙江省应用技术研究与开发计划 国家项目省级资助GX16C009

国家自然科学青年基金 51608495

中国地震局创新团队发展计划 中国大陆地区地震灾害模拟与评估

详细信息
    作者简介:

    孙柏涛, 男, 生于1961年。研究员。主要从事地震灾害评估等方面的研究。E-mail:sunbt@iem.cn

    通讯作者:

    姜琪, 男, 生于1993年。硕士研究生。主要从事结构试验方面的研究。E-mail:jiangqi_hit@163.com

Optimization Design and Mechanical Property Analysis of Small Steel Frame Reaction-wall Based on Large Scale Model Structure Test

  • 摘要: 为了降低模型制作加工成本和节约试验空间,开展结构缩尺比例小于1:10的小模型结构抗震试验。目前,国内相关的试验设备较少,本文将参考大型反力墙的结构形式,研究小型钢架反力墙,并对其进行优化设计和力学性能分析,为后续的制作加工提供理论依据。
  • 结构抗震试验是一门试验科学和技术,是研究和发展土木工程结构新材料、新结构、新工艺以及检验结构计算分析和设计理论的重要手段(马永欣等,2001)。在结构抗震试验中,为了使模型的几何形状、边界条件、材料与原型结构相似且具有等效荷载,缩尺比例一般不小于1:10,这也导致试验需要足够大的加载设备、试验空间以及较高的试验成本。

    如果仅为验证和研究结构体系、构件的抗震性能、破坏机理和破坏特征,如结构体系的薄弱环节等,对缩尺比例小于1:10的微模型结构进行抗震试验也可以达到目的。而结构抗震试验主要包括拟静力试验、拟动力试验和振动台试验,针对微模型结构,只有振动台试验有相应的试验设备——小型振动台,而拟静力和拟动力试验则缺少相应的设备。通过研究针对大型结构拟静力和拟动力试验的试验设备——大型反力墙,其结构形式可为小型反力墙的设计提供启示。

    本文对小型反力墙进行逐步优化,设计了3种方案,分别对其进行力学性能分析,从结构变形、功能性、实际操作性、经济性、美观性等角度选出最优设计方案,并对最优方案进行模态分析和疲劳验算。

    参考大型反力墙L型主、副反力墙+反力台的构造类型,本文设计的小型反力墙为三面相互垂直的形式。除了满足微模型结构试验精度外,小型反力墙及台座还应满足:①具有足够的抗压、抗拉、抗弯强度;②具有足够的刚度,满足“大荷载,小变形”;③具备良好的整体性,不因局部破坏而丧失整体使用功能;④方便反力架等其他配套设备的安装和拆卸,易于操作。本文针对上述要求设计了3种方案。

    设计方案1的具体尺寸为:长2640mm、宽1560mm、高2440mm,底部台座离地660mm;加载孔(板面上连接固定其它加载设备及模型的螺栓孔)为通透孔,直径18mm,孔中心间距的模数为180mm。

    反力墙和反力台的板面采用30mm厚的单层钢板,背面采用T形连接、点焊(极大避免通长焊缝带来的残余应力变形问题)的方式焊接20mm厚的钢板做加劲肋,加劲肋间距340mm,以增强方案1的整体刚度。反力墙长、短轴和反力台两两之间依靠直径为20mm的高强度螺栓连接。方案1的三维模型及侧视图见图 12

    图 1  方案1三维模型
    Figure 1.  A three-dimensional model of design#1
    图 2  方案1侧视图
    Figure 2.  The side-map of design#1

    设计方案2的具体尺寸为:长2400mm、宽1900mm、高1900mm,底部台座离地300mm;加载孔为通透孔,直径18mm,孔中心距模数为100mm。外侧支撑腿长100mm,内侧支撑腿长1000mm。

    反力墙的长、短轴和反力台均由截面为200mm×200mm的H型钢与水平向T形加劲肋焊接而成,H型钢间距300mm,加劲肋间距400mm,腹板与翼缘厚度均为20mm(加劲肋根据整体性要求、刚度要求以及美观程度自行设计)。短、长肢墙H型钢柱高1800mm,底部台座H型钢柱长2200mm。加载孔均匀对称地分布在一侧翼缘。外侧支撑腿置于地面之上,内侧支撑腿深埋入地下,以增强整体稳定性,减少不均匀沉降。方案2的三维模型及俯视图见图 34

    图 3  方案2三维模型
    Figure 3.  A three-dimensional model of design# 2
    图 4  方案2俯视图
    Figure 4.  The overlook-map of design#2

    设计方案3的具体尺寸为:长2380mm、宽1880mm、高2100mm,底部台座离地300mm;加载孔为通透孔,直径18mm,孔中心距模数为100mm,支撑腿长100mm。

    方案3对方案2钢架的H型钢柱进行了削减。短、长肢墙钢柱前翼缘尺寸由200mm减小到160mm,后翼缘尺寸由200mm减小到80mm,底部台座钢柱前后翼缘尺寸均由200mm减小到160mm,腹板与翼缘厚度依然为20mm。对比方案2,方案3的长、宽外廓尺寸基本不变,主、副反力墙的高度增加200mm;为了方便实际安装,支撑腿也全部修改为地上支撑,同时在钢架底部制作钢筋混凝土基础,后期安装时将支撑腿与基础中的钢预埋件焊接在一起。方案3的三维模型及正视图见图 56

    图 5  方案3三维模型
    Figure 5.  A three-dimensional model of design#3
    图 6  方案3正视图
    Figure 6.  The positive-map of design#3

    利用有限元软件ABAQUS对3种方案进行时程分析,求得各模型的变形数据,从应力、变形及经济的角度分析对比,选出最优设计方案。

    ABAQUS中只有相对单位制,首先需要统一尺寸单位,而模型体量较小,所以选取mm级的单位制,如表 1所示。本文研究的小型反力墙通体材料一致,适宜选取实体单元,网格划分类型为C3D8R(8节点的六面体),相应的3D实体模型与真实结构比较接近。

    表 1  ABAQUS统一单位
    Table 1.  Unified unit of ABAQUS
    长度 时间 荷载 质量 密度 应力
    m s N kg kg/m3 Pa
    mm s N 103kg 103kg/mm3 MPa
    下载: 导出CSV 
    | 显示表格

    大型反力墙为预应力钢筋混凝土结构,而本文研究对象体量小、受荷大、要求变形小,故选取金属材料进行制作。通过对中国金属市场的调研以及强度、刚度、经济等方面的考虑,选取市场常见的Q235B级钢。

    小型反力墙的设计前提是“大荷载,小变形”,而Q235B钢材在小应变情况下表现出近似弹性的性质,所以默认小型反力墙的使用一直处于线弹性范围。按照表 2输入参数。

    表 2  ABAQUS输入的钢材参数
    Table 2.  Input steel parameters in ABAQUS
    材料 密度/kN·m-3 杨氏模量/N·mm-2 泊松比
    Q235B钢 78.5 2.06×105 0.3
    下载: 导出CSV 
    | 显示表格

    为了精确模拟小型反力墙模型节点之间的传递运动关系,需严格定义各加劲肋与H型钢之间的约束,有限元软件中的约束类型包括绑定约束、刚体约束、耦合约束等,而模型中零部件的连接基本为焊接,符合绑定约束的主从面无滑移、分离情况,因此采用绑定约束。

    针对不同试验目的,得出的试验加载制度方案也不尽相同(王墩等,2014)。中国《建筑抗震试验方法规程》中指出,加载制度包括变位移控制加载、变荷载控制加载以及变荷载-变位移混合控制加载等方式,而本文模型所研究的内容是模拟实际试验时作动器的反力对小型反力墙的变形影响,因此适宜采用变荷载控制的加载制度。由于微模型试验的模型质量较小(小于500kg),经折算取荷载峰值为5kN,加载制度如图 7所示。为了横向对比3种方案,在其相同位置(图 810中的红色位置)施加同样的荷载,并对主、副反力墙角点提取变形数据进行分析。

    图 7  加载制度
    Figure 7.  Rules for loading
    图 8  方案1施加荷载位置
    Figure 8.  Load location of design#1
    图 9  方案2施加荷载位置
    Figure 9.  Load location of design#2
    图 10  方案3施加荷载位置
    Figure 10.  Load location of design#3

    有限元分析主要研究设计方案的变形情况,应力分析则是为了探究各设计方案的薄弱部位,进而对其进行加强,因此为了便于计算,将3种方案的短肢墙和长肢墙三维模型简化为杆系模型,将均布荷载简化为集中荷载,计算误差忽略不计。以方案3为例,简化杆系受力状态如图 11所示。

    图 11  简化杆系受力状态
    Figure 11.  The stress state of simplified rod system

    从变形挠度限值来看,小型反力墙的设计没有相应的参考文献以及规范,但其与大型反力墙的使用功能一致,因此两者对挠度限值的控制要求是相同的。参照新疆大学大型结构试验室大型反力墙的设计任务书中规定的最大挠度不超过L/1800,L为反力墙的有效高度(万宇,2008),方案1、2、3的有效高度分别为1800mm、1500mm、1700mm。

    经过对3种设计模型进行简化计算,得到了实际荷载值、挠度限值及有限元分析的最大变形值,见表 3

    表 3  实际荷载值、最大挠度限值及实际最大变形值
    Table 3.  Actual load value, maximum deflection limit and actual maximum deformation value
    方案1 方案2 方案3
    根部最大剪力/kN 5 5 5
    根部最大弯矩/kN·m 7.25 6.25 7.25
    最大挠度限值/mm 1 0.83 0.94
    短肢墙最大变形/mm 0.525 0.499 0.187
    长肢墙最大变形/mm 0.734 0.540 0.441
    下载: 导出CSV 
    | 显示表格

    分析表 3和3种方案的应力云图及变形云图(图 1213)发现,3种设计方案的最大变形都符合挠度限值的要求,但方案3的变形最小。除方案1反力墙的根部和底部支撑腿有明显的应力集中现象外,方案2、3均无明显的应力薄弱部位,但是方案2实际埋置、安装复杂。

    图 12  3种设计方案的应力云图
    Figure 12.  Stress cloud diagram of three design schemes
    图 13  3种设计方案的变形云图
    Figure 13.  The deformation cloud diagram of three design schemes

    经计算,在满足变形要求的前提下,3种方案反力墙根部的最大允许剪力和最大允许弯矩见表 4

    表 4  根部最大允许剪力和最大允许弯矩
    Table 4.  Maximum allowable shear and maximum allowable bending moment of root
    方案1 方案2 方案3
    根部允许最大剪力/kN 6.8 7.7 10.7
    根部允许最大弯矩/kN·m 12.24 11.55 18.19
    下载: 导出CSV 
    | 显示表格

    表 5给出了各方案的经济性指标,其中市场估价是指试件的材料费用,不包括试件加工费用和人工费用。从中可以看出方案3的重量相对较轻,且成本最低。因此,在满足微模型结构试验精度的前提下,从节约成本考虑,选取方案3进行加工制作。

    表 5  3种方案的自重及市场估价
    Table 5.  The weights and market valuation of the three designs
    设计方案 类型 材料 重量/kg 市场估价/万元
    方案1 板材满布形式 Q235B钢 5470 20
    方案3 型钢钢架 Q235B钢 4125 14
    方案3 型钢钢架 Q235B钢 3348 10
    下载: 导出CSV 
    | 显示表格

    结构的模态指其振动特性,包括频率、振型和阻尼比等,求解模态是为了获取结构的振动参数(张敏政,2015)。为了避免试验时作动器的输出频率与小型反力墙的固有频率重合而造成的共振现象,对方案3进行模态分析,求其固有频率。前3阶振型图及各参数见图 14表 6

    图 14  方案3的前3阶振型图
    Figure 14.  The first three-order vibration mode of design#3
    表 6  方案3的前3阶振型参数
    Table 6.  The first three-order vibration mode parameters of design#3
    振型 方向 特征值 固有频率/Hz 自振周期/s
    1阶振型 X轴正向平动+扭转 28072 26.666 0.0375
    2阶振型 X轴负向平动+扭转 46798 34.430 0.0290
    3阶振型 扭转 68393 41.622 0.0240
    下载: 导出CSV 
    | 显示表格

    考虑到小型反力墙将应用于动力试验,并结合表 6的数据,提出2种方法可以有效避免共振问题:①避开其固有频率;②输入荷载方向与其主要平动方向X轴垂直(即沿Y轴输入)。

    钢结构在连续交变荷载作用下会逐渐累计损伤、产生裂缝及裂缝逐渐扩展,直至最后破坏,因此要对其进行疲劳验算(张耀春,2011)。本文中的疲劳验算包含2种意义:①验算多循环荷载下结构不会破坏;②验算多循环荷载下结构的变形量不会增大。

    《钢结构设计规范GB 50017—2003》(中华人民共和国建设部,2003)中6.1.1节指出,直接承受动力荷载重复作用的钢结构构件及其连接,当应力变化的循环次数大于5×104时,应进行疲劳计算。假设方案3的循环次数能够达到规范中的限值要求,针对这一假设进行疲劳验算。

    根据《钢结构设计规范》中的6.2.1节,对常幅(所有应力循环内的应力幅值保持常量)疲劳按式(1)进行计算:

    $$ \Delta \sigma \le {\rm{[}}\Delta \sigma {\rm{ ]}} $$ (1)

    其中,Δσ为应力幅;[Δσ]为常幅疲劳的容许应力幅(N/mm2),按式(2)计算。对焊接部位,应力幅Δσ=σmax-σmin;对非焊接部位,折算应力幅Δσ=σmax-0.7σminσmax为计算部位每次应力循环中的最大拉应力(取正值);σmin为计算部位每次应力循环中的最小拉应力或压应力(拉应力取正值,压应力取负值)。

    $$[\Delta \sigma ] = {\left({\frac{C}{n}} \right)^{1/\beta }}$$ (2)

    其中,n为应力循环次数;Cβ为参数,根据《钢结构设计规范》附录E中的构件和连接类别取值(表 7)。

    表 7  参数值
    Table 7.  Related parameters
    参数 构件和连接类别
    1 2 3 4 5 6 7 8
    C 1940×1012 861×1012 3.26×1012 2.18×1012 1.47×1012 0.96×1012 0.65×1012 0.41×1012
    β 4 4 3 3 3 3 3 3
    下载: 导出CSV 
    | 显示表格

    根据制作类型,查询《钢结构设计规范》附录E中的构件和连接类别,方案3属于类别4,故根据表 6中4类的具体数值和公式(2),计算其常幅疲劳的容许应力幅。经计算,常幅疲劳的容许应力幅为351.96N/mm2,根据规范中对焊接部位应力幅公式的规定,计算出焊接面积,求得方案3中设计荷载5kN时的应力幅为0.14205N/mm2。从结果可以看出,即使假设方案3能达到疲劳循环次数,计算应力幅也远小于容许应力幅,因此方案3在多循环荷载下不至于发生破坏。

    为了探究方案3在不同循环次数下的最大变形情况有无增大,分别设置时长1s、10s和100s的3种工况,同时为了放大变形数值、便于宏观判断,将输入荷载值扩大为20kN(若变形值太小,即使最大变形值有变化也不易进行宏观判断,增大输入荷载仅仅是为了观察变形情况)。反力墙长轴变形量最大,故以其变形数据为例,验证常幅多循环荷载对变形情况的影响,3种分析部工况下的最大变形值见表 8

    表 8  3种分析部的最大变形值
    Table 8.  The maximum deformation value of the three types of analysis
    分析部时长 1s 10s 100s
    长肢墙最大变形值/mm 1.286 1.286 1.286
    下载: 导出CSV 
    | 显示表格

    表 8图 15可以看出,施加荷载的循环次数对最大变形值没有影响。综合2个角度的疲劳验算结果来看,应力循环次数对小型反力墙的影响基本可以忽略不计。

    图 15  3种分析部下长肢墙的变形曲线
    Figure 15.  The deformation curve of the long limb wall under the three types of analysis

    本文参照大型反力墙的结构形式、布置样式设计了3种小型反力墙的方案,并对它们进行了优化和ABAQUS有限元软件力学性能分析,从经济、应力和变形的角度比对选取出最优设计方案,并进行制作加工。再对最优设计方案进行模态分析和疲劳验算,针对实际情况可能会出现的共振现象,给出了具体的解决办法。优化后的小型反力墙可用于小模型结构抗震试验,为研究结构体系的抗震性能、破坏机理和破坏特征等提供试验支撑。

  • 图  1  方案1三维模型

    Figure  1.  A three-dimensional model of design#1

    图  2  方案1侧视图

    Figure  2.  The side-map of design#1

    图  3  方案2三维模型

    Figure  3.  A three-dimensional model of design# 2

    图  4  方案2俯视图

    Figure  4.  The overlook-map of design#2

    图  5  方案3三维模型

    Figure  5.  A three-dimensional model of design#3

    图  6  方案3正视图

    Figure  6.  The positive-map of design#3

    图  7  加载制度

    Figure  7.  Rules for loading

    图  8  方案1施加荷载位置

    Figure  8.  Load location of design#1

    图  9  方案2施加荷载位置

    Figure  9.  Load location of design#2

    图  10  方案3施加荷载位置

    Figure  10.  Load location of design#3

    图  11  简化杆系受力状态

    Figure  11.  The stress state of simplified rod system

    图  12  3种设计方案的应力云图

    Figure  12.  Stress cloud diagram of three design schemes

    图  13  3种设计方案的变形云图

    Figure  13.  The deformation cloud diagram of three design schemes

    图  14  方案3的前3阶振型图

    Figure  14.  The first three-order vibration mode of design#3

    图  15  3种分析部下长肢墙的变形曲线

    Figure  15.  The deformation curve of the long limb wall under the three types of analysis

    表  1  ABAQUS统一单位

    Table  1.   Unified unit of ABAQUS

    长度 时间 荷载 质量 密度 应力
    m s N kg kg/m3 Pa
    mm s N 103kg 103kg/mm3 MPa
    下载: 导出CSV

    表  2  ABAQUS输入的钢材参数

    Table  2.   Input steel parameters in ABAQUS

    材料 密度/kN·m-3 杨氏模量/N·mm-2 泊松比
    Q235B钢 78.5 2.06×105 0.3
    下载: 导出CSV

    表  3  实际荷载值、最大挠度限值及实际最大变形值

    Table  3.   Actual load value, maximum deflection limit and actual maximum deformation value

    方案1 方案2 方案3
    根部最大剪力/kN 5 5 5
    根部最大弯矩/kN·m 7.25 6.25 7.25
    最大挠度限值/mm 1 0.83 0.94
    短肢墙最大变形/mm 0.525 0.499 0.187
    长肢墙最大变形/mm 0.734 0.540 0.441
    下载: 导出CSV

    表  4  根部最大允许剪力和最大允许弯矩

    Table  4.   Maximum allowable shear and maximum allowable bending moment of root

    方案1 方案2 方案3
    根部允许最大剪力/kN 6.8 7.7 10.7
    根部允许最大弯矩/kN·m 12.24 11.55 18.19
    下载: 导出CSV

    表  5  3种方案的自重及市场估价

    Table  5.   The weights and market valuation of the three designs

    设计方案 类型 材料 重量/kg 市场估价/万元
    方案1 板材满布形式 Q235B钢 5470 20
    方案3 型钢钢架 Q235B钢 4125 14
    方案3 型钢钢架 Q235B钢 3348 10
    下载: 导出CSV

    表  6  方案3的前3阶振型参数

    Table  6.   The first three-order vibration mode parameters of design#3

    振型 方向 特征值 固有频率/Hz 自振周期/s
    1阶振型 X轴正向平动+扭转 28072 26.666 0.0375
    2阶振型 X轴负向平动+扭转 46798 34.430 0.0290
    3阶振型 扭转 68393 41.622 0.0240
    下载: 导出CSV

    表  7  参数值

    Table  7.   Related parameters

    参数 构件和连接类别
    1 2 3 4 5 6 7 8
    C 1940×1012 861×1012 3.26×1012 2.18×1012 1.47×1012 0.96×1012 0.65×1012 0.41×1012
    β 4 4 3 3 3 3 3 3
    下载: 导出CSV

    表  8  3种分析部的最大变形值

    Table  8.   The maximum deformation value of the three types of analysis

    分析部时长 1s 10s 100s
    长肢墙最大变形值/mm 1.286 1.286 1.286
    下载: 导出CSV
  • 马永欣, 郑山锁, 2001.结构试验.北京:科学出版社.
    万宇, 2008.新疆大学结构实验室X方向反力墙有限元分析.乌鲁木齐: 新疆大学.
    王墩, 赵海琼, 吕西林, 2014.建筑结构拟静力试验方法的加载制度.四川建筑科学研究, 40(3):54-59. doi: 10.3969/j.issn.1008-1933.2014.03.014
    张敏政, 2015.地震工程的概念和应用.北京:地震出版社.
    张耀春, 2011.钢结构设计原理.北京:高等教育出版社.
    中华人民共和国建设部, 2003.钢结构设计规范GB 50017-2003.北京: 中国计划出版社.
  • 期刊类型引用(3)

    1. 余佳科,景立平,王展,陆新宇,齐文浩. 桩-土-钢、砼结构动力相互作用试验对比研究. 震灾防御技术. 2024(01): 151-159 . 本站查看
    2. 张尹,杨嘉胤. 反力墙结构受力性能有限元分析. 低温建筑技术. 2022(08): 67-70 . 百度学术
    3. 谭晓晶,文龙,许国山,吴斌,张九峰,宋军. 6层装配式钢筋混凝土盒子结构的子结构拟动力试验与数值模拟分析. 建筑结构. 2020(S1): 539-544 . 百度学术

    其他类型引用(3)

  • 加载中
图(15) / 表(8)
计量
  • 文章访问数:  95
  • HTML全文浏览量:  39
  • PDF下载量:  12
  • 被引次数: 6
出版历程
  • 收稿日期:  2018-03-15
  • 刊出日期:  2018-12-01

目录

/

返回文章
返回