隔震曲线连续梁桥粘滞阻尼器参数优化分析

第41卷第2期

哈　尔　滨　工　程　大　学　学　报JournalofHarbinEngineeringUniversity

Vol.41№.2

Feb.2020

隔震曲线连续梁桥粘滞阻尼器参数优化分析

(东南大学混凝土及预应力混凝土结构教育部重点实验室,江苏南京210096)

王浩,沈惠军,张寒,郑文智,沙奔,李爱群

摘　要:为了提升粘滞阻尼器在隔震曲线梁桥中的减隔震效果,本文对其力学性能参数进行优化分析。基于AN-SYS建立了某三跨隔震曲线连续梁桥的有限元模型,采用非线性动力时程法,分析了粘滞阻尼器速度指数和阻尼系数等参数对地震作用下支座位移和墩底剪力的影响,并确定了粘滞阻尼器参数的合理取值范围。在此基础上,基于零阶优化算法,以墩底剪力绝对值之和为目标函数,对粘滞阻尼器进行了参数优化,并对比分析了优化前后的减震效果。结果表明:阻尼系数对支座位移和墩底剪力的影响十分显著;速度指数对支座位移和墩底剪力的影响与阻尼系数有关;粘滞阻尼器参数优化后,支座切向和径向位移减震率分别在47%和22.2%以上,中墩支座的切向位移减震率高于边墩支座,而边墩支座的径向位移减震率高于中墩支座;边、中墩切向剪力差值由956kN减小为220kN,径向剪力差值由441kN减小为130kN,其受力相对更为均衡。

关键词:曲线连续梁桥;减隔震;粘滞阻尼器;敏感性分析;参数优化;零阶优化算法;隔震支座;地震响应DOI:10.11990/jheu.201906080

网络出版地址:http://www.cnki.net/kcms/detail/23.1390.u.20191021.1412.002.html中图分类号:U441　文献标志码:A　文章编号:1006-7043(2020)02-282-07

Parameteroptimizationanalysisofviscousdampersforisolated

continuouscurvedgirderbridges

(KeyLaboratoryofConcreteandPrestressedConcreteofMinistryofEducation,SoutheastUniversity,Nanjing210096,China)

WANGHao,SHENHuijun,ZHANGHan,ZHENGWenzhi,SHABen,LIAiqun

Abstract:Toimprovetheseismicperformenceofviscousdampersinisolatedcontinuouscurvedgirderbridges,themechanicalparametersofviscousdampersareoptimizedandanalyzed.Thefiniteelementmodelofathree-spani-solatedcontinuouscurvedgirderbridgeisperfermedthebasisonANSYS.Theinfluenceofthevelocityexponentanddampingcoefficientofviscousdampersareanalyzedusingthenonlineardynamictimehistorymethod.Therea-sonablerangeoftheparametersofviscousdampersaredeterminedusingthezero-orderoptimizationalgorithm,theparametersofviscousdampersarefurtheroptimized,withthesumoftheabsolutevaluesoftheshearforceatthebottomofpiersastheobjectivefunction.Resultsshowthatthedampingcoefficienthasasignificanteffectonthedisplacementofbearingsandthebaseshearforceofpiers.Theeffectofthevelocityexponentonthedisplacementofbearingsandthebaseshearforceofpiersisrelatedtothedampingcoefficient.Aftertheparametersofviscousdampersareoptimized,thereductionratiofortangentialandradialdisplacementofthebearingsreach47%and22.2%,respectively.Thereductionratiofortangentialdisplacementofthebearingsatthemiddlepiersarehigherthanthoseatthesidepiers.Meanwhile,thedampingratiosforradialdisplacementofthebearingsatthesideandmiddlepiersarehigherthanthoseatthemiddlepiers.Thedifferenceofthetangentialshearbetweenthesideandmiddlepiersisreducedfrom956kNto220kNandthedifferenceoftheradialshearbetweenthesideandmiddlepiersisreducedfrom441kNto130kN,whichensureabalanceddistributionofforceoneachpier.

Keywords:continuouscurvedgirderbridge;seismicisolation;viscousdamper;sensitivityanalysis;parameterop-timization;zero-orderoptimizationalgorithm;isolationbearings;seismicresponse

　　曲线梁桥因其外形美观、适应性强的特点,而被

收稿日期:2019-06-24.网络出版日期:2019-10-21.基金项目:国家自然科学基金项目(51578151,57978155);国家“万

人计划”青年科技人才项目(W03070080);江苏省重点研发计划(BE2018120).

作者简介:王浩,男,教授,博士生导师,“长江学者奖励计划”青年

学者.

通信作者:王浩,E-mail:wanghao1980@seu.edu.cn.

广泛应用于各类山区公路和城市道路的桥梁工程。由于其平面不规则,地震作用下的受力更为复杂[1]。当采用隔震体系时,为了满足其位移需求,边墩处通常采用四氟滑板橡胶支座,但已有研究表明,地震作用下,由于边墩四氟滑板橡胶支座可发生相对滑移,上部结构的地震作用将主要由中墩来承

第2期王浩,等:隔震曲线连续梁桥粘滞阻尼器参数优化分析

·283·

担,从而会导致边墩、中墩受力不均衡。针对此问题,目前通常采用隔震支座加阻尼器的组合方法对结构的地震响应进行控制。粘滞阻尼器的力-位移滞回曲线非常饱满,耗能限位能力强,能在不增加地震惯性力的情况下,有效地控制桥梁结构的位移响应[2-3]。此外,粘滞阻尼器还具有构造简单、便于施工、受激励频率和温度影响小的特点,因此很受工程技术人员的青睐,成为减小桥梁结构关键部位地震响应的主要手段之一[4-7]。

为了使粘滞阻尼器达到预期的效果,有必要对其参数进行优化分析。对此,国内外诸多学者开展了相关研究。文献[8]通过建立与粘滞阻尼器参数相关的拉格朗日方程,利用顺序搜索法求解阻尼器参数的最优解。文献[9]为优化粘滞阻尼器对双塔斜拉桥的减震效果,应用最小二乘回归分析法建立关键截面参数与阻尼参数之间的数学模型,以控制截面内力和变形最小为原则,通过求解拟合方程的极值得到最优设计参数。文献[10]针对现有阻尼器参数优化方法的缺陷,基于斜拉桥的随机地震响应,通过最小化粘滞阻尼器的阻尼力得到阻尼器参数的最优解。以上关于粘滞阻尼器参数优化的研究

对象均是斜拉和悬索等大跨度柔性桥梁,针对连续梁桥的相关研究还不够系统,特别是对曲线连续梁桥则更为少见[11-12]。显然,对隔震曲线连续梁桥采用粘滞阻尼器进行减震控制优化研究有待进一步加强。

为此,本文以某三跨隔震曲线连续梁桥为研究对象,采用非线性动力时程法,分析了粘滞阻尼器阻尼系数和速度指数对其地震响应的影响。在此基础上,基于零阶优化算法,以各墩墩底剪力之和为目标函数,进行了粘滞阻尼器参数优化分析,旨在为隔震曲线连续梁桥的减震设计与研究提供参考。

1　隔震曲线连续梁桥分析模型

该隔震曲线连续梁桥跨度为3×35m,位于曲率半径为200m的平曲面上。主梁采用现浇预应力混凝土箱梁,其截面为单箱双室,桥墩采用圆形双柱墩,其具体尺寸分别见图1(a)、(b)。中墩、边墩处支座类型分别为铅芯橡胶支座、四氟滑板橡胶支座。该隔震曲线连续梁桥构造图如图1所示,其中,桥墩编号为P1~P4,支座编号为B1~B4。由设计资料可知,该桥址区的抗震设防烈度为9度,场地类别为二类。

　　基于ANSYS建立该桥的空间有限元模型,如图2所示。主梁、横梁和桥墩采用弹性梁单元模拟,中墩铅芯橡胶支座、边墩四氟滑板橡胶支座水平切向和径向均采用COMBIN40单元模拟,竖向采用COMBIN14单

Fig.1　Thestructureoftheisolatedcontinuouscurvedgirderbridge(unit:m)

图1　隔震曲线连续梁桥构造(单位:m)

元模拟。铅芯橡胶支座和四氟滑板橡胶支座的力学模型分别如图2(a)和(b)所示,其优化前的力学模型参数如表1所示。主梁和支座之间采用刚性单元连接,桥墩底部固接,未考虑桩土相互作用。

k0/(kN·mm-1)k1/(kN·mm-1)

18.3

2.8

铅芯橡胶支座(Y4Q920∗315G1.2)(P2和P3墩)初始刚度屈服后刚度屈服力屈服位移

Fy/kN323

y/mm22.7

表1　隔震支座简化力学模型参数

Table1　Simplifiedmechanicalparametersofisolationbearings

竖向刚度/1422

四氟滑板橡胶支座(LNR(H)-d670∗188)(P1和P4墩)

滑移前刚度/(kN·mm-1)

2.84

滑移系数0.03

屈服位移y/mm14.5

kN

竖向刚度/

(kN·mm-1)

931

·284·

哈　尔　滨　工　程　大　学　学　报第41卷

图2　隔震曲线连续梁桥有限元模型

Fig.2　Finiteelementmodeloftheisolatedcontinuouscurvedgirderbridge

2　粘滞阻尼器参数敏感性分析

2.1　地震动输入

根据场地剪切波速,从PEER地震动数据库中选取了3条地震波,其信息如表2所示。根据《公路桥梁抗震设计细则》[13],将地震动主、次分量的加速度峰值分别调整为0.4、0.34g。本文仅研究了水平双向地震动输入的情况,未考虑竖向地震动的影响。

Table2　Accelerationrecordsoftheearthquake

表2　地震加速度记录

记录台站

地震事件

时间震级

NorthernCalif-03FerndaleCityHall19546.50SuperstitionHills-02E1CentroImp.Co.Cent19876.54

LomaPrieta

HollisterCityHall

19896.93

2.2　参数敏感性分析

地震作用下,当该桥边墩四氟滑板橡胶支座发生相对滑移时,该支座提供的恢复力不再增加,上部结构的地震作用将主要通过铅芯橡胶支座传递给中墩来承担,从而导致边、中墩受力不均衡。因此,在两边墩处沿径向和切向各设置一个粘滞阻尼器。以期增大结构阻尼,减小支座位移,同时有效调节结构内力在各墩之间的分配。下面以内侧四氟滑板橡胶支座(B1)和内侧桥墩(P1、P2)为研究对象,分析粘滞阻尼器参数对其地震响应的影响。其中,阻尼系数c取值范围为0~3MN·(m/s)-α,速度指数α分别取0.3、0.5、0.7、1.0。NorthernCalif-03地震动作用下的计算结果如图3、图4所示。

图3　c和α对B1支座位移的影响

Fig.3　TheeffectofcandαonthedisplacementoftheB1bearing

　　由图3可知,当α一定时,支座切向和径向位移随着c的增大而减小,且逐渐趋于平缓,c取值在0~1MN·(m/s)-α内,对支座切向和径向位移的影响

较为显著;当c一定时,支座切向和径向位移随着α

的减小而减小。当c和α分别取3MN·(m/s)-α和0.3时,支座切向和径向位移减幅均高达88.89%。

第2期王浩,等:隔震曲线连续梁桥粘滞阻尼器参数优化分析

·285·

由此说明,当c和α取值合理时,粘滞阻尼器可以对隔震曲线连续梁桥的位移响应进行较好的控制。由图4可知,当α一定时,边墩墩底切向和径向剪力随着c的增大而增大,而中墩的变化规律相反;当c一定时,α越小,对中、边墩墩底切向和径向剪力的影响越大。

图4　c和α对桥墩墩底剪力的影响

Fig.4　Theeffectofcandαontheshearforceatthebottomofpiers

　　由此可知,存在一个最优的c和α,可以使得边、中墩受力更为均衡,同时有效减小支座位移。以上粘滞阻尼器的参数敏感性分析提供了一个参数取值范围,为了进一步得到最优参数的具体值,以下进行优化分析。

利用罚函数将上述有约束问题的目标函数转化为无约束问题的目标函数[15],其拟合公式为:minF(X,pk)=^f+f0pk[∑X(xi)+

m2

i=1

m3

n

3　粘滞阻尼器参数优化

3.1　零阶优化算法

零阶优化算法是在一定次数的抽样基础上,拟合设计变量、状态变量和目标函数的响应函数,从而寻求最优解,故又可称其为子问题方法。函数曲线的形式可采用线性拟合、平方拟合或平方交叉项拟合。对于含设计变量和状态变量的约束问题,可采用罚函数将其转化为无约束的最小值问题[14]。

该方法是建立在目标函数及状态变量近似的基础上,若采用平方交叉项拟合,则目标函数的拟合公式为:

^f=a0+

aixi+∑∑bijxixj∑iij

n

n

n

式中:xi、xj为设计变量;a0、ai、bij为拟合系数,拟合

系数随迭代过程而变,由加权最小二乘法确定。

(1)

式中:xi为设计变量;gi、hi、wi为状态变量;X、G、H、W为对应的罚函数;f0为目标函数参考值;pk为响应面参数;无约束目标F(X,pk)为响应面函数,随设计变量X及响应面参数pk而变化。3.2　优化模型的建立

基于零阶优化算法,对粘滞阻尼器参数进行优化分析。优化过程中,阻尼系数c和阻尼指数α为设计变量,其取值范围为:

lu-+ρ1cρ≤c≤cì1111ï

ïcl-ρ≤c≤cu+ρï21221

(3)ílu

-+ïα1ρ1≤α1≤α1ρ1

ïluïα2

î-ρ1≤α2≤α2+ρ1

ll

=c2=式中:上标u、l代表上下限,其中,c1

ul

=1.8MN·(m/s)-α、α1=0.6MN·(m/s)、c=c2

H(hi)+∑W(wi)]∑=i1i=1

G(gi)∑i=1

m1

+(2)

-α

u

1

·286·

哈　尔　滨　工　程　大　学　学　报第41卷

luu=0.3、α1=α2=1.0,ρ1=2×10-2。设计变量初α2

0.6MN·(m/s)-α、α2=0.5。这些值的选用参考了以上参数敏感性分析的结果。

将支座位移作为含有约束条件的状态变量,以各墩底切向剪力绝对值之和为目标函数fob:

4Ui-T≤0.2+ρ2

=fobmin(∑Fi),(4)

Ui-R≤0.15+ρ2i=1

始值设定为c1=0.8MN·(m/s)-α、α1=0.5和c2=

3.3　优化结果分析

响应最大,基于上述优化模型,开展粘滞阻尼器参数优化分析。优化后的切向粘滞阻尼器参数c和α分别为0.775MN·(m/s)-α、0.411;径向粘滞阻尼器参

由于NorthernCalif-03地震作用下的结构地震

{

数c和α分别为0.386MN·(m/s)-α、0.633。

为了分析隔震曲线连续梁桥粘滞阻尼器优化后

的减震效果,基于上述优化后的参数,采用非线性动力时程法,进行了有、无粘滞阻尼器时支座位移和墩底剪力的对比分析,其结果分别如表3、4所示,相应的支座位移和墩底剪力时程曲线如图5、6所示。

式中:Fi为单墩墩底切向剪力;Ui-T为支座切向位

移;Ui-R为支座径向位移;ρ2为容差,其值为1×10-3。

Table3　Thedisplacementofbearingsanditsdampingrate

表3　支座位移及其减震率

支座位移

无阻尼器/m有阻尼器/m减震率/%

0.360.1947

B1-T

0.370.1948.6

B2-T

0.370.1948.6

B3-T

0.380.2047.4

B4-T

0.190.1426.3

B1-R

0.180.1422.2

B2-R

0.170.1323.5

B3-R

0.170.1229.4

B4-R

注:T代表切向、R代表径向。

Table4　Theshearforceatthebottomofpiersanditsdampingrate

表4　墩底剪力及其减震率

无阻尼器/kN有阻尼器/kN减震率/%

1001-106.8

484

F1-T

1331

F2-T

39.4

806

1385

F3-T

38.3

854

1026-139.1

429

F4-T

墩底剪力F1-R610-82.1335

F2-R

6956812.0

F3-R

7767059.1

F4-R

740-77.5

417

目标函数3628

3687-1.6

注:T代表切向、R代表径向。

图5　支座位移时程

Fig.5　Thedisplacementtimehistoryofthebearing

　　由表3和图5可知,阻尼器优化后,支座切向和径向位移减震率分别在47%和22.2%以上,切向位移和径向位移分别在0.2m和0.15m以内。两中墩支座和两边墩支座的位移减震率分别保持一致,而中墩支座与边墩支座间的位移减震率相差较大。中墩支座的切向位移减震率高于边墩支座,而边墩支座的径向位移减震率高于中墩支座。

由表4和图6可知,阻尼器优化后,边墩墩底

切向剪力明显增大,最大增加了139.1%,最小也达到了106.8%,中墩墩底切向剪力减震率在38.3%以上。边墩墩底径向剪力增大的幅度相对较小,但也在77.5%以上,而中墩墩底径向剪力的减震率最大仅9.1%。值得注意的是,设置阻尼器后,虽然各墩墩底剪力之和仅变化了1.6%,但中墩与边墩的墩底剪力差值大幅度减小,其中,切向剪力差值由956kN减小为220kN,径向剪力差值

第2期

王浩,等:隔震曲线连续梁桥粘滞阻尼器参数优化分析

·287·

由441kN减小为130kN,边、中墩受力更为均衡。总体而言,在边墩处设置粘滞阻尼器,对边墩墩底

剪力的影响大于中墩,对墩底切向剪力的影响大于径向。

图6　墩底剪力时程

Fig.6　Thesheartimehistoryatthebottomofpiers

4　结论

bridge[J].JournalofSouthChinaUniversityofTechnology

1)合理的粘滞阻尼器参数可以使隔震曲线连续梁桥边、中墩在地震作用下的受力更为均衡,同时有效控制桥梁位移。

2)采用优化后的粘滞阻尼器参数,支座位移减幅明显,切向和径向位移的减幅分别在47%和22.2%以上。中墩支座的切向位移减幅高于边墩支座,而边墩支座的径向位移减幅高于中墩支座。

3)设置粘滞阻尼器前后的各墩墩底剪力之和基本无变化,但墩底剪力差值大幅度减小,边、中墩受力更为均衡。

[3]马长飞,胡可,汪正兴,等.芜湖长江公路二桥斜置粘

(naturalscienceedition),2019,47(4):90-98.

65-68.

滞阻尼器动力性能监测[J].世界桥梁,2019,47(5):MAChangfei,HUKe,WANGZhengxing,etal.Dynamic

behaviormonitoringforinclinedviscousdampersofSecond[4]SINGHMP,MORESCHILM.Optimalplacementof

2019,47(5):65-68.

WuhuChangjiangRiverhighwaybrige[J].Worldbridges,dampersforpassiveresponsecontrol[J].Earthquakeengi-[5]IMPOLLONIAN,PALMERIA.Seismicperformanceof

neering&structuraldynamics,2002,31(4):955-976.

参考文献:

[1]王丽,周锡元,闫维明.曲线梁桥地震响应的简化分析

方法[J].工程力学,2006,23(6):77-84.

buildingsretrofittedwithnonlinearviscousdampersandad-jacentreactiontowers[J].Earthquakeengineering&struc-WANGLi,ZHOUXiyuan,YANWeiming.Simplifieda-[6]BOUGTEBY,RAYT.Choicebetweenseriesandparallel

turaldynamics,2018,47(5):1329-1351.

nalysismethodforseismicresponseofcurvedbridges[J].

connectionsofhystereticsystemandviscousdamperfor

[2]徐艳,童川,李建中.设置纵桥向粘滞阻尼器的斜拉桥

Engineeringmechanics,2006,23(6):77-84.

seismicprotectionofstructures[J].Earthquakeengineering

的简化动力模型[J].华南理工大学学报(自然科学版),2019,47(4):90-98.

[7]POLLININ,LAVANO,AMIRO.Minimum-costoptimiza-&structuraldynamics,2018,47(1):237-244.

tionofnonlinearfluidviscousdampersandtheirsupportingmembersforseismicretrofitting[J].Earthquakeengineer-XUYan,TONGChuan,LIJianzhong.Seismicmitigation

analysisofviscousdampersforcurvedcontinuousgirder

[8]AYDINE,BODUROGLUMH,GUNEYD.Optimal

ing&structuraldynamics,2017,46(12):1941-1961.

·288·

哈　尔　滨　工　程　大　学　学　报

第41卷

damperdistributionforseismicrehabilitationofplanar

buildingstructures[J].Engineeringstructures,2007,29[9]巫生平,张超,房贞政.斜拉桥粘滞阻尼器设计方案及

参数回归分析[J].桥梁建设,2014,44(5):21-26.

analysisofseismiccontrolschemesforcontinuouscurved

(2):176-185.

[13]中华人民共和国交通运输部.JTG/TB02-01-2008,公

35(10):157-161,173.

girderbridges[J].Journalofvibrationandshock,2016,

WUShengping,ZHANGChao,FANGZhenzheng.Designschemesandparameterregressionanalysisofviscousdamp-2008.

路桥梁抗震设计细则[S].北京:人民交通出版社,MinistryofTransportofthePeople′sRepublicofChina.

[10]王波,马长飞,刘鹏飞,等.基于随机地震响应的斜拉

44(5):21-26.

ersforcable-stayedbridge[J].Bridgeconstruction,2014,桥粘滞阻尼器参数优化[J].桥梁建设,2016,46(3):WANGBo,MAChangfei,LIUPengfei,etal.Parameter

JTG/TB02-01-2008,GuidelinesforseismicdesignofPress,2008.

highwaybridges[S].Beijing:ChinaCommunications[14]王浩,王付全,李爱群,等.大跨度缆索支撑桥梁分阶

17-22.

optimizationofviscousdamperforcable-stayedbridge

111-116.

段有限元模型修正[J].工程力学,2009,26(10):

[11]LIUYi,QIXingjun,WANGYijian,etal.Seismicmiti-basedonstochasticseismicresponse[J].Bridgeconstruc-tion,2016,46(3):17-22.

WANGHao,WANGFuquan,LIAiqun,etal.Multi-phaseFEmodelupdatingonlong-spancable-supported

gationanalysisofviscousdampersforcurvedcontinuous

[12]李正英,蒋林均,李正良.曲线连续梁桥不同减隔震方

2011,90-93:1230-1233.

girderbridge[J].Appliedmechanicsandmaterials,案对比分析[J].振动与冲击,2016,35(10):157-[15]郭彤,李爱群,费庆国,等.零阶与一阶优化算法在悬

bridges[J].Engineeringmechanics,2009,26(10):111-116.

2007,26(4):35-38.

索桥模型修正中的应用对比分析[J].振动与冲击,GUOTong,LIAiqun,FEIQingguo,etal.Application

LIZhengying,JIANGLinjun,LIZhengliang.Comparative

161,173.

comparisonbetweenzero-orderandfirst-orderoptimization

methodsinmodelupdatingofsuspensionbridges[J].

Journalofvibrationandshock,2007,26(4):35-38.

本文引用格式:

王浩,沈惠军,张寒,等.隔震曲线连续梁桥粘滞阻尼器参数优化分析[J].哈尔滨工程大学学报,2019,41(2):282-288.

WANGHao,SHENHuijun,ZHANGHan,etal.Parameteroptimizationanalysisofviscousdampersforisolatedcontinuouscurvedgirderbridges[J].

JournalofHarbinEngineeringUniversity,2019,41(2):282-288.