2. 哈尔滨工业大学 土木工程学院,黑龙江 哈尔滨 150090;
3. 扬州大学 建筑科学与工程学院,江苏 扬州 225127
2. School of Civil Engineering, Harbin Institute of Technology, Harbin 150090, China;
3. College of Civil Science and Engineering, Yangzhou University, Yangzhou 225127, China
目前,在实验室进行的结构抗震试验方法主要有3种:拟静力试验、地震模拟振动台试验、拟动力子结构试验(又称为联机试验或混合试验).拟静力试验由于加载制度简单、稳定, 而应用最为广泛,但是其缺点在于加载制度主观确定,与实际地震作用历程无关,不能直接考虑地震动对结构的影响,较多地应用于构件试验,很少进行整体结构模型试验.振动台试验能最真实地模拟地震对结构的影响,但是由于振动台尺寸及承载能力的限制,即使目前世界上最大的振动台——日本的E-Defense(台面尺寸20 m×15 m),对于很多结构模型来说,也往往只能进行缩尺结构试验.由于小尺度结构模型的动力相似率很难满足要求,且受尺寸效应的影响,尤其是进入非线性后,振动台试验不能很好地反映结构真实的抗震性能及抗震能力.虽然,近年来,有人将子结构技术引入到振动台试验中,提出了子结构振动台试验,为振动台试验发展注入了新的思路,然而其依然突破不了振动台台面尺寸的局限[1-3].子结构试验方法源于日本学者Hakuno[4]在1969年提出的拟动力试验,其将整体结构模型中非线性比较复杂的关键部位划分成试验子结构进行实际物理试验,而将其他部分通过数值模拟实现,大大降低了试验规模,使进行足尺模型抗震性能试验成为可能.此外,由于振动台昂贵的建设费用,目前,国内拥有振动台试验能力的院校及科研机构远远少于拥有子结构试验的硬件条件(反力墙和电液伺服加载系统)的单位.因此,采用子结构试验方法普遍被看作是一种经济合理、应用前景广泛的试验方法.
进入21世纪之后,随着网络技术及有限元的飞速发展,子结构试验方法出现了向网络化、精细化发展的趋势,多个国家开始逐步构建网络化结构试验室,其中规模最大的是美国1999年启动的NEES (network for earthquake engineering simulation)计划[5],并开发了多种拟动力子结构试验平台(包括OpenFresco[6]、UI-Simcor[7]、TUT[8]、NetSLab[9]、HyTest等).为了推动我国子结构试验方法的发展,许国山等[10]在OpenFresco基础上,开发了基于OpenSEES-OpenFresco-dSPACE和OpenSEES-OpenFresco-LabView-dSPACE两套混合模拟系统.其中,OpenSEES实现计算子结构模拟,OpenFresco实现软件接口,LabView/dSPACE实现试验子结构控制加载.这些复杂的子结构试验系统极大地推动了子结构试验方法的发展,但它们对于非抗震试验技术领域人员来说,需要熟悉多个软件的工作性能,并掌握各个软件之间的接口,尤其是上述系统多是基于MTS FlexTest平台基础的研发(如OpenFresco、TUT、NetSLab、HyTest等),程序的开发和扩展还涉及到试验设备控制系统的复杂底层技术,对非专业人员来说开发难度较大,从而浪费大量精力,而不能将主要精力集中在本身课题的研究.此外,很多子结构试验系统均需要使用人员开发相应软件,比如目前使用最多的基于OpenFresco的试验系统,只能采用其开发的几个常规的试验单元,如果是复杂特殊试验模型形式,比如本文为了减少加载控制难度,选用的在结构反弯点处作为界面,就没有相应试验单元可用.
基于上述考虑,本文提出一种不涉及任何试验设备控制系统的外部命令控制方法,采用目前学者所熟悉的开发软件Matlab实现运动方程的求解、试验子结构控制加载以及数据采集等功能,调用有限元软件OpenSEES,建立了通用子结构试验系统.其中,Matlab部分可以作为黑匣子提供给系统应用人员,其只用根据本身所需利用OpenSEES建立适合自身结构模型的计算子结构数值模型即可.当然,如果使用人员熟悉其他有限元软件,如ABAQUS和ANSYS等,只需修改Matlab与其接口程序即可.
1 系统组成及原理本文建立的子结构试验系统通过Matlab和OpenSEES混合编程实现,主要包括Matlab计算加载控制采集模块、OpenSEES数值子结构模拟模块以及两者的接口,如图 1所示.
由图 1可知,本文提出的子结构试验系统,是独立于MTS电液伺服设备控制系统的,其可以直接在外接的计算机上编程实现,再通过数据采集卡与MTS电液伺服设备控制系统硬件的I/O(输入/输出)通道连接即可,从而避开了设备控制系统复杂的底层编程技术,减少了对控制系统本身的要求(通常设备控制系统的计算通道数目和函数非常有限,且价格昂贵),降低了试验难度和费用.图 1中试验子结构控制加载系统是以MTS电液伺服系统为例,对于任何其他控制系统,只要找到与硬件系统对应的I/O通道,本文提出的子结构试验系统即可应用,具有很好的通用性.
1.1 Matlab模块Matlab软件具有强大的矩阵运算功能,并且含有丰富的应用工具箱功能,便于编程.因此,选择Matlab软件实现运动方程计算、试验子结构加载控制以及数据采集功能.
1.1.1 计算加载控制模块计算加载控制模块是本子结构试验系统的核心部分,主要实现运动方程的求解和试验子结构系统的控制加载,其流程如图 2所示.求解运动方程如下:
$ \mathit{\boldsymbol{M}}{\mathit{\boldsymbol{a}}_{i + 1}} + \mathit{\boldsymbol{C}}{\mathit{\boldsymbol{v}}_{i + 1}} + {\mathit{\boldsymbol{R}}_{\rm{N}}}({\mathit{\boldsymbol{d}}_{i + 1}}) + {\mathit{\boldsymbol{R}}_{\rm{E}}}({\mathit{\boldsymbol{d}}_{i + 1}}) = - \mathit{\boldsymbol{MI}}{\mathit{\boldsymbol{\overline x}} _{{\rm{g}}, i + 1}} $ | (1) |
式中:d、v、a分别为整体结构的位移、速度和加速度n×1阶向量(n为整体结构的自由度);M为n×n阶质量矩阵; C为n×n阶Rayleigh阻尼矩阵; xg, i+1为第i+1步地震动记录; I为n×n阶单位向量; RN为计算子结构n×1阶恢复力向量; RE为试验子结构n×1阶恢复力向量,其通常采用电液伺服作动器内置的力传感器测得其反力后换算得到.
算法采用平均加速度法,其算法系数选取a=0、β=1/4、γ=1/2,则第i+1步的加速度和速度表达式如下:
$ {\mathit{\boldsymbol{a}}_{i + 1}} = {\rm{ }}\frac{4}{{\Delta {t^2}}}\left( { - {\mathit{\boldsymbol{d}}_i} - \Delta t{\mathit{\boldsymbol{v}}_i} - \frac{{\Delta {t^2}}}{4}{\mathit{\boldsymbol{a}}_i} + {\mathit{\boldsymbol{d}}_{i + 1}}} \right) $ | (2) |
$ {\mathit{\boldsymbol{v}}_{i + 1}} = - \frac{2}{{\Delta t}}{\mathit{\boldsymbol{d}}_i} - {\mathit{\boldsymbol{v}}_i} + \frac{2}{{\Delta t}}{\mathit{\boldsymbol{d}}_{i + 1}} $ | (3) |
式中:Δt是积分时间间隔.
将式(2)、(3) 代入式(1) 中整理得到
$ {\mathit{\boldsymbol{R}}_{\rm{N}}}({\mathit{\boldsymbol{d}}_{i + 1}}) + {\mathit{\boldsymbol{K}}_{{\rm{PD}}}}{\mathit{\boldsymbol{d}}_{i + 1}} + {\rm{ }}{\mathit{\boldsymbol{R}}_{\rm{E}}}({\mathit{\boldsymbol{d}}_{i + 1}}) = {\mathit{\boldsymbol{F}}_{{\rm{EQ}}, i + 1}} $ | (4) |
$ {\mathit{\boldsymbol{K}}_{{\rm{PD}}}} = \frac{{4\mathit{\boldsymbol{M}}}}{{\Delta {t^2}}} + \frac{{2\mathit{\boldsymbol{C}}}}{{\Delta t}} $ | (5) |
式中:KPD为拟动力刚度.
$ \begin{array}{l} {\mathit{\boldsymbol{F}}_{{\rm{EQ}}, i + 1}} = - \mathit{\boldsymbol{MI}}{\mathit{\boldsymbol{\overline x}} _{{\rm{g}}, i + 1}} + \mathit{\boldsymbol{M}}{\mathit{\boldsymbol{a}}_i} + \\ \;\;\;\;\;\;\;\;\;\;\;\;\left( {\frac{{4\mathit{\boldsymbol{M}}}}{{\Delta t}} + \mathit{\boldsymbol{C}}} \right){\mathit{\boldsymbol{v}}_i} + \left( {\frac{{4\mathit{\boldsymbol{M}}}}{{\Delta {t^2}}} + \frac{{2\mathit{\boldsymbol{C}}}}{{\Delta t}}} \right){\mathit{\boldsymbol{d}}_i} \end{array} $ | (6) |
将FEQ, i+1称之为第i+1步的等效力命令,由式(6) 可知,在第i+1步时其为已知值,而且只有在式(4) 等式左边(称之为等效力反馈值)等于该步等效力命令时,即可得第i+1步的位移反应,从而利用式(2) 和(3) 计算出第i+1步的速度、加速度反应.
为使控制等效力平稳追踪等效力命令,本文采用PID(proportional integral derivative)控制,控制量为试验子结构的位移,采用闭环控制系统,即反馈控制,使等效力反馈值(等式(4) 左边)平稳渐近地趋于等效力命令(等式(4) 右边),其已被作者广泛应用于多自由度子结构位移控制[11-12].
1.1.2 采集模块采集模块主要实现试验模型位移采集、向电液伺服加载系统发送位移以及采集其反力等功能.本试验系统选用研华PCI 1712数据采集卡实现,当然可以选用研华其他系列采集卡、NI(national instruments)系列采集卡等.Matlab软件为研华系列采集卡提供了模拟输入、模拟输出的驱动控制模块,同时在研华公司主页上提供实时更新的各种情况下I/O实现的.m文件供用户免费下载使用,本文就是利用在该网站下载的程序AsynchronousOneBufferedAI.m和StaticAO.m修改实现的.目前,子结构试验系统选用的PCI 1712,有1 M转换速度的12位A/D (analog-to-digital)转换器,卡上带有FIFO (first in、first out)缓冲器(可存储1 k A/D采样值和32 k D/A转换数据),提供16路单端或8路差分的模拟量输入(也可单端差分混合使用),2路为D/A模拟量输出通道,本文中选用的均为差分输入.具体管脚线如图 3所示,具体输入、输出管脚线接法详见官网上提供的用户手册.
此外,为了方便子结构试验过程中随时观测、应变采集等,在每步发送采集完成后,在程序中设置暂停命令,等待用户反馈后进行下步运算.
1.2 计算子结构模块计算子结构模块为子结构试验系统的子程序,可供Matlab主程序调用,其可根据用户整体模型及子结构模型选择情况,设置单个或多个计算子结构模块,从而实现分布式子结构试验.
OpenSEES是Open System for Earthquake Engineering Simulation的简写,是一款面向对象的开源性有限元软件,初期主要由加州大学伯克利分校研发.其最大的特点是源代码的开放性,可为全世界相关研究人员所用,同时任何研发人员编写的被公认编码随时可加入到软件模块中,具有很好的应用和发展前景,得到越来越多研究人员的青睐.因此,本文选用OpenSEES软件进行计算子结构模拟.当然,用户可以根据自身情况,选用其他有限元软件,如ABAQUS、ANSYS等,只需修改Matlab与其接口程序即可.
1.3 接口 1.3.1 Matlab与OpenSEES接口子结构试验系统采用了两种接口方式实现编程软件Matlab和有限元软件OpenSEES的数据传输功能:① 读写文本文件,当需要向对方软件传输数据时,需要向两个文本文件写入数据,一个记录需传输的数据,另一个为对方软件的开关命令,比如当文件为0时,对方软件暂停等待读入数据,当文件变为1时,对方软件才开始读入数据;同理,当该软件输出数据后,就读对方软件给予的开关命令文件,等待对方软件传输数据;② Socket通信变量传输,通过Socket命令建立软件间的TCP/IP (transmission control protocol/internet protocol)通讯地址,实现软件间的变量直接传输.其中,读写文本方法比较繁琐,且读写文本增加了系统的风险,但其最大的优势就是可以实现任意软件之间的传输,也即采用此法用户目前就可用任意有限元软件.而对于Socket通信变量传输,直接实现两个软件间相应变量传输,大大提高了效率,减少系统风险,但目前仅实现Matlab和OpenSEES之间的传输,其他有限元软件和Matlab的接口还需进一步探索.
1.3.2 Matlab和MTS电液伺服设备控制系统接口Matlab和MTS电液伺服设备控制系统之间接口通过硬件数据采集卡实现,即将数据采集卡插入到Matlab实现的计算机主机内,通过屏蔽线将其与MTS电液伺服设备控制系统对应的I/O通道连接起来.数据采集卡的连线和驱动参考1.1中采集模块介绍.
2 系统验证在前面建立的系统基础上,目前完成了3个大比例尺模型的子结构混合模拟试验.其中,Matlab模块、OpenSEES数值子结构模拟模块以及二者的接口在如图 4所示的普通计算机上完成.该计算机实现本文所提出子结构混合模拟系统的程序计算和显示功能,可以在试验时实时追踪计算程序中的所有数据,从而监控试验进程,如图 5所示的等效力命令和等效力反馈值的比较曲线.同时,该计算机通过和图 6所示的数据卡PCI 1712连接,实现程序计算计算机和MTS电液伺服设备控制计算机之间的数据传输.图 7为MTS电液伺服设备控制计算机.该计算机实现从数据卡PCI 1712接收到数据通过MTS控制柜发送到电液伺服加载器上,实现模型的加载,同时将电液伺服作动器反馈的信号再传递给PCI 1712.文中地震动参数输入均为南北方向El Centro地震动.
该模型总层数为6层,平面尺寸为4.0 m×6.0 m,底部一层为框支层,层高为3.6 m,柱截面为400 mm×550 mm,梁截面为300 mm×600 mm,上部5层为配筋砌块砌体剪力墙层,层高为3 m,墙厚为190 mm.试验子结构模型选取底部3层,采用本文提出的子结构试验系统完成了峰值加速度为0.35g、0.70g、1.10g、2.20g、2.60g、3.10g、3.50g系统子结构试验,研究结构模型的抗震性能和抗震能力[13-14].试验子结构模型如图 8所示,其是在哈尔滨工业大学结构与抗震试验大厅完成,水平加载装置采用在每层放置1个电液伺服作动器控制加载实现.
图 9为峰值加速度为2.20g时的等效力命令和等效力反馈的比较时程曲线.从图 9可以看出,等效力反馈值可以很好地追踪到等效力命令,说明本文提出的子结构系统的控制模块是可行的.各工况下的位移时程及滞回曲线如图 10和图 11所示,这里仅给出试验子结构1层结果.从图 9~11可以看出,采用本文提出的子结构试验系统可以很好地评价多自由度结构体系抗震能力.
该模型总层数为12层,各层层高均为3.0 m,地梁高度为0.6 m,平面尺寸4.2 m×3.0 m,连梁截面均为200 mm×400 mm.试验子结构模型选取底部3层,采用本文提出的子结构试验系统完成了峰值加速度为0.35g、0.70g、1.10g、2.20g的系统子结构试验,研究装配式结构模型在地震作用下的水平、竖向连接性能[15].试验子结构模型如图 12所示,其是在哈尔滨工业大学结构与抗震试验大厅完成.水平加载装置采用4个电液伺服作动器加载控制:在1层和2层中轴线位置分别作用1个作动器,在3层中轴线两侧对称位置分别作用2个作动器.图 13为峰值加速度为2.20g时的等效力命令和等效力反馈的比较时程曲线.从图 13同样可以看出,等效力反馈值可以很好地追踪到等效力命令,说明本文提出的子结构系统的控制模块可行.试验子结构1层的滞回曲线如图 14所示.从图中可以看出,结构模型已经进入了下降段,说明本文提出的子结构试验系统可以得到复杂多自由度体系模型在地震作用下的下降段试验数据,具有很好的适用性.试验子结构1层在不同峰值加速度下的位移时程曲线见图 15.
该模型总层数为6层单榀单跨钢框架:首层层高为4.2 m,其他层高均为4.0 m,跨度为8.0 m,梁、柱均为H型钢,柱为HW 600 mm×600 mm×20 mm×30 mm,梁为HW 500 mm×300 mm×14 mm×22 mm.试验子结构选取为底层,上部5层为计算子结构,采用本文提出的子结构试验系统完成了峰值加速度为1.40g、3.10g、6.20g、8.00g、12.00g、16.00g的系统子结构试验.为实现子结构界面的准确模型,实际试验子结构选择在2层弹性反弯点处,即1层半模型,其中上部的半层作为加载层.试验模型采用1:2的缩尺模型,试验模型图如图 16所示.水平加载装置采用2个水平作动器加载控制,竖向力采用液压千斤顶控制.由于需要考虑倾覆力矩的影响,需要实时改变竖向力的大小,因此整个模型加载采用力位移混合控制的方法实现(水平采用位移控制, 竖向采用力控制).图 17为峰值加速度为16.00g时的等效力命令和等效力反馈的比较时程曲线.从图 17同样可以看出,等效力反馈值可以很好地追踪到等效力命令,说明本文提出的子结构系统的控制模块可行.图 18给出了试验子结构1层位移时程,从图中可以看出,在16.00g时经历了反向峰值后,进入正向后已经达不到12.00g时的正向峰值,说明结构已经进入了下降段,因此在正向峰值经历后回到零点即停止试验.由此可以说明,本文提出的试验系统可以得到力位移混合加载情况下的下降段数据,具有很好的鲁棒性.图 19给出了不同峰值加速度下的滞回曲线.随后,在完成子结构混合模拟试验后进行了单个作动器作用的低周往复试验,发现已经达不到16.00g时的峰值力,也说明了试验模型已经进入了下降段.
(1) 本文提出了基于Matlab和OpenSEES的通用子结构试验系统,主要包括Matlab计算加载控制采集模块、OpenSEES数值子结构模拟模块以及二者的接口.
(2) 本文提出的子结构试验系统,是独立于MTS电液伺服设备控制系统的,直接在外接的计算机上编程实现,避开了设备控制系统复杂的底层编程技术,减少了对控制系统本身的要求,降低了试验难度和费用,具有很好的通用性和推广度.
(3) 通过多个大比例尺结构试验模型的系统子结构试验,验证了本文提出的子结构试验系统的可行性.
[1] |
王向英, 田石柱. 子结构地震模拟振动台混合试验原理与实现[J]. 地震工程与工程振动, 2009, 29(4): 46 WANG Xiangying, TIAN Shizhu. Study on structural hybrid testing method with shaking table[J]. Journal of Earthquake Engineering and Engineering Vibration, 2009, 29(4): 46 |
[2] |
李振宝, 李晓亮, 唐贞云, 等. 基于振动台的动力子结构试验界面反力获取方法[J]. 地震工程与工程振动, 2011, 31(3): 65 LI Zhenbao, LI Xiaoliang, TANG Zhenyun, et al. Research on the methods of evaluating reaction force for dynamic sub-structure experiments using shaking table[J]. Journal of Earthquake Engineering and Engineering Vibration, 2011, 31(3): 65 |
[3] |
汪强, 王进廷, 金峰, 等. 基于振动台的RTDHT试验及其在双层钢架结构上的应用[J]. 工程力学, 2010, 27(10): 57 WANG Qiang, WANG Jinting, JIN Feng, et al. Real-time dynamic hybrid testing based on shaking-table and its application to a two storey rigid frame[J]. Engineering Mechanics, 2010, 27(10): 57 |
[4] |
HAKUNO M, SHIDOWARA M, HARA T. Dynamic destructive test of a cantilever beam, controlled by an analog-computer[J]. Transactions of the Japan Society of Civil Engineering, 1969, 171(12): 1 |
[5] |
BUCKLE I, REITHERMAN R, GEORGE E, et al. The consortium for the network for earthquake engineering simulation[C]//13th World Conference on Earthquake Engineering. Vancouver:[s.n.], 2004: 4016.
|
[6] |
TAKAHASHI Y, FENVES G L. Software framework for distributed experimental-computational simulation of structural systems[J]. Earthquake Engineering and Structural Dynamics, 2006, 35: 267 DOI:10.1002/(ISSN)1096-9845 |
[7] |
KWON O, ELNASHAI A S, SPENCER B F, et al. UI-SimCor : a global platform for hybrid distributed simulation[C]//Ninth Canadian Conference on Earthquake Engineering. Ottawa:[s.n.], 2007: 26-29.
|
[8] |
邱法维, 潘鹏, 钱稼茹, 等. 建筑结构拟动力实验软件的开发与应用[J]. 建筑结构学报, 2000, 21(5): 22 QIU Fawei, PAN Peng, QIAN Jiaru, et al. Development and application of pseudodynamic testing software[J]. Journal of Building Structures, 2000, 21(5): 22 |
[9] |
肖岩, 胡庆, 郭玉荣, 等. 结构拟动力远程协同试验网络平台的开发研究[J]. 建筑结构学报, 2005, 26(3): 122 XIAO Yan, HU Qing, GUO Yurong, et al. A network platform for remote pseudo-dynamic testing[J]. Journal of Building Structures, 2005, 26(3): 122 |
[10] |
许国山, 郝伟, 陈永盛, 等. 基于OpenFresco-LabVIEW-dSpace的混合试验系统研究[J]. 工程力学, 2013, 30(3): 417 XU Guoshan, HAO Wei, CHEN Yongsheng, et al. Experimental validation on hybrid testing system based on OpenFresco-LabVIEW-dSpace[J]. Engineering Mechanics, 2013, 30(3): 417 |
[11] |
陈再现, 吴斌, 王焕定, 等. 多自由度等效力控制拟动力子结构试验研究[J]. 工程力学, 2011, 28(8): 138 CHEN Zaixian, WU Bin, WANG Huanding, et al. Equivalent force control substructure pseudo-dynamic experimental study on a mdof model[J]. Engineering Mechanics, 2011, 28(8): 138 |
[12] |
CHEN Zaixian, XU Guoshan, WU Bin, et al. Equivalent force control method for substructure pseudo-dynamic test of a full-scale masonry structure[J]. Earthquake Engineering and Structural Dynamics, 2014, 43(7): 969 DOI:10.1002/eqe.v43.7 |
[13] |
王凤来, 陈再现, 王焕定, 等. 底部框支配筋砌块短肢砌体剪力墙结构抗震性能试验研究[J]. 土木工程学报, 2009, 42(11): 71 WANG Fenglai, CHEN Zaixian, WANG Huanding, et al. Experimental study on the seismic behavior of a frame-supported reinforced masonry short-limb shera wall structure[J]. China Civil Engineering Journal, 2009, 42(11): 71 DOI:10.3321/j.issn:1000-131X.2009.11.010 |
[14] |
陈再现, 王凤来, 杨同盖, 等. 底框砌体剪力墙震损房屋抗震加固拟动力子结构试验[J]. 土木建筑与环境工程, 2011, 33(6): 1 CHEN Zaixian, WANG Fenglai, YANG Tonggai, et al. Experimental analysis of seismic performance of strengthening frame supported reinforced concrete masonry short-leg shear wall structure damaged by earthquake using the substrcture[J]. Journal of Civil, Architectural & Environmental Engineering, 2011, 33(6): 1 |
[15] |
陈再现, 姜洪斌, 张家齐, 等. 预制钢筋混凝土剪力墙结构拟动力子结构试验研究[J]. 建筑结构学报, 2011, 32(6): 41 CHEN Zaixian, JIANG Hongbin, ZHANG Jiaqi, et al. Pseudo-dynamic substructure test on precast reinforced concrete shear wall structure[J]. Journal of Building Strutures, 2011, 32(6): 41 |