CN111351622B - 一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法 - Google Patents
一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法 Download PDFInfo
- Publication number
- CN111351622B CN111351622B CN202010165052.XA CN202010165052A CN111351622B CN 111351622 B CN111351622 B CN 111351622B CN 202010165052 A CN202010165052 A CN 202010165052A CN 111351622 B CN111351622 B CN 111351622B
- Authority
- CN
- China
- Prior art keywords
- rubber support
- max
- restoring force
- lead core
- core rubber
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0041—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/086—Learning methods using evolutionary algorithms, e.g. genetic algorithms or genetic programming
Abstract
本发明提供一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法,包括以下步骤:首先输入试验得到的铅芯橡胶支座第三圈滞回曲线,平均非对称的两峰值点使其关于原点中心对称;在第三圈滞回曲线加载段和卸载段各取两个点,得到支座恢复力模型参数的取值范围;输入遗传算法的终端性能指标函数,并结合得到的参数取值范围,通过遗传算法进行参数识别;得到再加载刚度、卸载刚度和屈服力的识别结果,构造出其恢复力模型;通过恢复力模型得到的滞回环与测得的铅芯橡胶支座第一圈滞回曲线对比,可以得到试验过程中的滑移。本发明有效简化了铅芯橡胶支座恢复力模型参数识别方法,且可以计算出铅芯橡胶支座在试验过程中的滑移。
Description
技术领域
本发明涉及土木工程领域,具体是一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法。
背景技术
隔震建筑是将地震作用产生的能量与隔震层以上的结构隔离,从而减弱地震作用对上部结构的破坏,提高建筑物的稳定性,已经成为建筑设计的新理念。其中铅芯橡胶隔振支座以其耐久性好、性能稳定、构造便捷等特别备受青睐。
在有限元模拟中,铅芯橡胶隔振支座的剪切性能常采用无下降段的双折线恢复力模型,其恢复力参数主要包括水平等效刚度Kreload、屈服后刚度Kdeload、屈服荷载(剪力)Vyield等。GB/T 20688.1-2007《橡胶支座第一部分:隔震橡胶支座试验方法》中对橡胶隔震支座(天然橡胶支座LNR和铅芯橡胶支座LRB和高阻尼橡胶支座HDR)竖向压缩刚度和剪切性能的测试方法、计算方法均做了详细的规定。通过试验得到的滞回曲线特征点值,便可以得到竖向刚度、水平等效刚度、屈服后刚度等恢复力特征参数。然而在实际数据处理中不可避免地存在采样频率、读数等系统误差和偶然误差,导致试验中得到的滞回曲线与规范中相差较大,计算得到的恢复力特征参数不能十分有效地反应隔震支座真实的滞回性能。
发明内容
本发明的目的是为解决上述技术问题,提供一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法。
为达到上述目的,本发明采用的方法是:一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法,包括以下步骤:
S1:输入试验得到的铅芯橡胶支座的100%剪切性能曲线;
S2:根据测得的铅芯橡胶支座第三圈滞回曲线,平均非对称的两峰值点使其关于原点中心对称;
S3:在第三圈滞回曲线加载段和卸载段各取两个点,得到支座恢复力模型参数(再加载刚度、卸载刚度和屈服力)的取值范围;
S4:输入遗传算法的终端性能指标函数,并结合S3得到的取值范围,通过遗传算法进行参数识别;
S5:得到再加载刚度、卸载刚度和屈服力的识别结果,得到其恢复力模型;
S6:通过恢复力模型得到的滞回环与测得的铅芯橡胶支座第三圈滞回曲线对比,可以得到试验过程中的滑移。
进一步地,步骤S2中,根据测得的铅芯橡胶支座第三圈滞回曲线,识别出滞回环的最大值(最大位移Xmax和Fmax)和最小值(Xmin和Fmin),根据最值结果,分别将非对称的滞回环的最大值和最小值的绝对值平均,得到对称的正负位移(Dmax和Dmin)和剪力(Vmax和Vmin),从而将非对称的滞回环转化为对称滞回环:
Vmax=-Vmin=(Fmax-Fmin)/2 式(1)
Dmax=-Dmin=(Xmax-Xmin)/2 式(2)
进一步地,步骤S3中,在再加载段和卸载段任取两个相差较远的点,求得其斜率为分别为K1和K2;由此确定其两折线恢复力模型参数中再加载刚度Kreload和卸载刚度Kdeload取值范围分别为[0.9K1,1.1K1]和[0.9K2,1.1K2],屈服荷载Fyield取值范围为[0,Vmax],最值点分别为(Dmax,Vmax)和剪力(Dmin,Vmin)。
进一步地,步骤S4中,遗传算法的终端性能指标函数为:
Fitnessfunction=ABS(S-△w)/△w 式(3)
其中S为简化滞回环围成的面积,△w为试验滞回环围成的面积。
S=∫V(x)dx 式(4)
△w=∫F(x)dx 式(5)
进一步地,步骤S5中,在MATLAB中调用遗传算法进行参数识别或优化,代码如下输入:
A=[1 0 0;0 1 0;0 0 1];
B=[1.1K1;1.1K2;Vmax];
Lowerbound=[0.9K1,0.9K2,0];
Uowerbound=[1.1K1,1.1K2,Vmax];
X=ga(@Fitnessfunction,3,A,B,[],[],Lowerbound,Uowerbound);
得到Kreload=X(1),Kdeload=X(2),Vyield=X(3)。
进一步地,步骤S6中,根据S5中得到再加载刚度Kreload、卸载刚度Kdeload和屈服力Vyield,得到其恢复力模型,另外两个特征点(D1,V1)和(D2,V2)计算如下:
D1=-D2=[Dmax(Kreload+Kdeload)-2Vmax]/(Kdeload-Kreload) 式(6)
V1=-V2=Vmax+Kdeload(D1-Dmax) 式(7)
进一步地,步骤S7中,根据测得的铅芯橡胶支座第一圈滞回曲线,得到再加载点位移值△可以得到试验过程中的滑移△slip:
△slip=△-Vyield/Kdeload 式(8)
本发明具有以下有益效果:
本发明有效简化了铅芯橡胶支座恢复力模型参数识别方法,且可以计算出铅芯橡胶支座在试验过程中的滑移,具有较高的精度。
附图说明
图1为本发明基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法的流程图;
图2为本发明铅芯橡胶支座第三圈滞回曲线示意图;
图3为本发明中恢复力模型示意图;
图4为本发明铅芯橡胶支座第一圈滞回曲线示意图。
具体实施方式
下面结合具体实施方式和附图对本发明的技术方案作进一步的介绍。
如图1所示,本发明的一种应用于新型耗能节点的简化三折线恢复力模型,包括以下步骤:
S1:输入试验得到的铅芯橡胶支座的100%剪切性能曲线;
S2:根据测得的铅芯橡胶支座第三圈滞回曲线,平均非对称的两峰值点使其关于原点中心对称;
S3:在第三圈滞回曲线加载段和卸载段各取两个点,得到支座恢复力模型参数(再加载刚度、卸载刚度和屈服力)的取值范围;
S4:输入遗传算法的终端性能指标函数,并结合S3得到的取值范围,通过遗传算法进行参数识别;
S5:得到再加载刚度、卸载刚度和屈服力的识别结果,得到其恢复力模型;
S6:通过恢复力模型得到的滞回环与测得的铅芯橡胶支座第三圈滞回曲线对比,可以得到试验过程中的滑移。
步骤S2中,根据测得的铅芯橡胶支座第三圈滞回曲线,如图2所示,根据测得的铅芯橡胶支座第三圈滞回曲线,识别出滞回环的最大值(最大位移Xmax和Fmax)和最小值(Xmin和Fmin),根据最值结果,分别将非对称的滞回环的最大值和最小值的绝对值平均,得到对称的正负位移(Dmax和Dmin)和剪力(Vmax和Vmin),从而将非对称的滞回环转化为对称滞回环::
Vmax=-Vmin=(Fmax-Fmin)/2 式(1)
Dmax=-Dmin=(Xmax-Xmin)/2 式(2)
步骤S3中,在再加载段和卸载段任取两个相差较远的点,求得其斜率为分别为K1和K2;由此确定其两折线恢复力模型参数中再加载刚度Kreload和卸载刚度Kdeload取值范围分别为[0.9K1,1.1K1]和[0.9K2,1.1K2],屈服荷载Fyield取值范围为[0,Vmax],最值点分别为(Dmax,Vmax)和剪力(Dmin,Vmin)。
步骤S4中,遗传算法的终端性能指标函数为:
Fitnessfunction=ABS(S-△w)/△w 式(3)
其中S为简化滞回环围成的面积,△w为试验滞回环围成的面积。
S=∫V(x)dx 式(4)
△w=∫F(x)dx 式(5)
步骤S5中,在MATLAB中调用遗传算法进行参数识别或优化,代码如下输入:
A=[1 0 0;0 1 0;0 0 1];
B=[1.1K1;1.1K2;Vmax];
Lowerbound=[0.9K1,0.9K2,0];
Uowerbound=[1.1K1,1.1K2,Vmax];
X=ga(@Fitnessfunction,3,A,B,[],[],Lowerbound,Uowerbound);
得到Kreload=X(1),Kdeload=X(2),Vyield=X(3)。
步骤S6中,根据S5中得到再加载刚度Kreload、卸载刚度Kdeload和屈服力Vyield,得到其恢复力模型,如图3所示,另外两个特征点(D1,V1)和(D2,V2)计算如下:
D1=-D2=[Dmax(Kreload+Kdeload)-2Vmax]/(Kdeload-Kreload) 式(6)
V1=-V2=Vmax+Kdeload(D1-Dmax) 式(7)
步骤S7中,根据测得的铅芯橡胶支座第一圈滞回曲线,如图4所示,得到试验屈服点位移值△可以得到试验过程中的滑移△slip:
△slip=△-Vyield/Kdeload 式(8)。
Claims (5)
1.一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法,其特征在于,包括以下步骤:
S1:输入试验得到的铅芯橡胶支座的100%剪切性能曲线;
S2:根据测得的铅芯橡胶支座第三圈滞回曲线,平均非对称的两峰值点使其关于原点中心对称;
S3:在第三圈滞回曲线加载段和卸载段各取两个点,得到支座恢复力模型参数的取值范围,所述的支座恢复力模型参数包括再加载刚度、卸载刚度和屈服力;
其中步骤S3中,在再加载段和卸载段任取两个相差较远的点,求得其斜率为分别为K 1 和K 2 ;由此确定其两折线恢复力模型参数中再加载刚度K reload 和卸载刚度K deload 取值范围分别为[0.9 K 1 ,1.1 K 1 ]和[0.9 K 2 ,1.1 K 2 ],屈服荷载V yield 取值范围为[0,V max ],最值点分别为(D max ,V max )和(D min ,V min );
S4:输入基于铅芯橡胶支座滞回环面积相对误差的终端性能指标函数,并结合S3得到的取值范围,通过遗传算法进行参数识别;其中遗传算法的终端性能指标函数为:
Fitnessfunction=ABS(S-△ w )/△ w 式(1)
其中S为简化滞回环围成的面积,△w为试验滞回环围成的面积;
S=∫V(x) dx 式(2)
△ w =∫F(x) dx 式(3)
其中V为简化滞回环的剪力,F为试验滞回环的剪力,x为剪切位移;
S5:得到再加载刚度、卸载刚度和屈服力的识别结果,得到其恢复力模型;
S6:通过恢复力模型得到的滞回环与测得的铅芯橡胶支座第三圈滞回曲线对比,可以得到试验过程中的滑移。
2.根据权利要求1所述的一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法,其特征在于:步骤S2中,根据测得的铅芯橡胶支座第三圈滞回曲线,识别出滞回环的位移、剪力的最大值和最小值,分别将滞回环的位移、剪力的最大值和最小值绝对值平均,将非对称的滞回环转化为对称滞回环:
V max =-V min =(F max -F min )/2 式(4)
D max =-D min =(X max -X min )/2 式(5)
其中Xmax、Xmin为初始滞回环的最大位移和最小位移; Fmax、Fmin为初始滞回环的最大剪力和最小剪力, Dmax、Dmin为处理后滞回环的最大位移和最小位移; Vmax和Vmin为处理后滞回环的最大剪力和最小剪力。
3.根据权利要求1所述的一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法,其特征在于:步骤S5中,是通过在MATLAB中调用遗传算法进行参数识别或优化实现的。
4.根据权利要求3所述的一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法,其特征在于:根据S5中得到再加载刚度K reload 、卸载刚度K deload 和屈服荷载V yield ,得到其恢复力模型,另外两个特征点(D 1 ,V 1)和(D 2 ,V 2)计算如下:
D 1 =- D 2 =[D max (K reload +K deload )-2V max ]/(K deload -K reload ) 式(6)
V 1 =-V 2 =V max +K deload (D 1 - D max ) 式(7) 。
5.根据权利要求4所述的一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法,其特征在于:根据测得的铅芯橡胶支座第一圈滞回曲线,得到再加载点位移值△可以得到试验过程中的滑移△slip:
△slip=△-V yield / K deload 式(8)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010165052.XA CN111351622B (zh) | 2020-03-11 | 2020-03-11 | 一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010165052.XA CN111351622B (zh) | 2020-03-11 | 2020-03-11 | 一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111351622A CN111351622A (zh) | 2020-06-30 |
CN111351622B true CN111351622B (zh) | 2020-11-27 |
Family
ID=71194450
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010165052.XA Active CN111351622B (zh) | 2020-03-11 | 2020-03-11 | 一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111351622B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112854513B (zh) * | 2021-01-13 | 2021-11-30 | 东南大学 | 一种粘滞阻尼器力学性能系数及滑移识别方法 |
CN114444191B (zh) * | 2022-02-11 | 2022-12-09 | 四川大学 | 一种重型储仓设备的铅芯橡胶隔震支座参数设计方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09217784A (ja) * | 1996-02-14 | 1997-08-19 | Nitta Ind Corp | 免震装置 |
JP2001074625A (ja) * | 1999-09-07 | 2001-03-23 | Showa Electric Wire & Cable Co Ltd | せん断変形量測定装置 |
WO2006088600A2 (en) * | 2005-01-20 | 2006-08-24 | Bae Systems Information And Electronic Systems Integration Inc. | Microradio design, manufacturing method and applications for the use of microradios |
KR20130043855A (ko) * | 2011-10-21 | 2013-05-02 | 김재욱 | 복원성을 조절할 수 있는 교량용 납 면진받침 |
CN103440383A (zh) * | 2013-09-02 | 2013-12-11 | 东南大学 | 一种自复位屈曲约束支撑的滞回性能分析方法 |
CN104677587A (zh) * | 2015-03-25 | 2015-06-03 | 重庆邮电大学 | 基于智能磁控橡胶支座的桥墩缓冲隔振多功能试验平台 |
CN104879437A (zh) * | 2015-03-27 | 2015-09-02 | 中国电子工程设计院 | 一种隔振支撑基座的滑移复合耗能隔振方法 |
CN106639024A (zh) * | 2016-10-20 | 2017-05-10 | 兰州理工大学 | 限位及自复位橡胶‑滑移隔震支座的施工方法 |
CN109933855A (zh) * | 2019-02-18 | 2019-06-25 | 北京交通大学 | 基于遗传算法的自复位支撑结构抗震性能优化方法与系统 |
CN110414161A (zh) * | 2019-08-01 | 2019-11-05 | 广州建筑股份有限公司 | 一种梭形滞回曲线的参数确定方法及系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102539236B (zh) * | 2011-12-13 | 2013-11-13 | 北京工业大学 | 一种检测橡胶隔震支座定剪切竖向拉伸刚度的方法 |
WO2014193328A1 (en) * | 2013-05-31 | 2014-12-04 | Toker Mehmet | Natural rubber or synthetic rubber elastomer-based earthquake isolator with rigid polyurethane core |
CN103512740B (zh) * | 2013-09-27 | 2016-04-13 | 中交公路长大桥建设国家工程研究中心有限公司 | 一种对支座进行静动力多功能测试的装置 |
CN204694465U (zh) * | 2015-06-23 | 2015-10-07 | 宁波大学 | 橡胶支座在复合受力状态下抗压和抗剪强度测试装置 |
CN105865766B (zh) * | 2016-04-01 | 2020-04-17 | 中国电力科学研究院有限公司 | 滞回特性模型预测装置及方法 |
-
2020
- 2020-03-11 CN CN202010165052.XA patent/CN111351622B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09217784A (ja) * | 1996-02-14 | 1997-08-19 | Nitta Ind Corp | 免震装置 |
JP2001074625A (ja) * | 1999-09-07 | 2001-03-23 | Showa Electric Wire & Cable Co Ltd | せん断変形量測定装置 |
WO2006088600A2 (en) * | 2005-01-20 | 2006-08-24 | Bae Systems Information And Electronic Systems Integration Inc. | Microradio design, manufacturing method and applications for the use of microradios |
KR20130043855A (ko) * | 2011-10-21 | 2013-05-02 | 김재욱 | 복원성을 조절할 수 있는 교량용 납 면진받침 |
CN103440383A (zh) * | 2013-09-02 | 2013-12-11 | 东南大学 | 一种自复位屈曲约束支撑的滞回性能分析方法 |
CN104677587A (zh) * | 2015-03-25 | 2015-06-03 | 重庆邮电大学 | 基于智能磁控橡胶支座的桥墩缓冲隔振多功能试验平台 |
CN104879437A (zh) * | 2015-03-27 | 2015-09-02 | 中国电子工程设计院 | 一种隔振支撑基座的滑移复合耗能隔振方法 |
CN106639024A (zh) * | 2016-10-20 | 2017-05-10 | 兰州理工大学 | 限位及自复位橡胶‑滑移隔震支座的施工方法 |
CN109933855A (zh) * | 2019-02-18 | 2019-06-25 | 北京交通大学 | 基于遗传算法的自复位支撑结构抗震性能优化方法与系统 |
CN110414161A (zh) * | 2019-08-01 | 2019-11-05 | 广州建筑股份有限公司 | 一种梭形滞回曲线的参数确定方法及系统 |
Non-Patent Citations (3)
Title |
---|
"A simplified simulation method of friction pendulum bearings";Wang Yanhua、Feng Yan、Wu Jing;《Journal of Southeast University ( English Edition)》;20181231;第34卷(第4期);第480-487页 * |
"铅芯叠层橡胶支座恢复力模型研究";王建强、王丽娟;《世界地震工程》;20051231;第21卷(第2期);第151-153页 * |
"铅芯橡胶隔震支座研究";赵毅 等;《科技与创新》;20181231(第1期);第33页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111351622A (zh) | 2020-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111351622B (zh) | 一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法 | |
CN103308381B (zh) | 一种疲劳裂纹扩展速率归一化预测方法 | |
CN104236889B (zh) | 橡胶密封圈的贮存寿命优化设计方法 | |
CN109212413B (zh) | 线性马达带宽测量方法 | |
CN112528411B (zh) | 一种基于模态减缩的几何非线性结构噪声振动响应计算方法 | |
CN110909413A (zh) | 一种混凝土面板堆石坝水平位移预测精确方法及系统 | |
CN111428361A (zh) | 一种适用于多种循环载荷的寿命预测方法 | |
CN111881564B (zh) | 一种关于机械结构变幅疲劳寿命预测方法 | |
CN110414161B (zh) | 一种梭形滞回曲线的参数确定方法及系统 | |
CN108110336B (zh) | 电池包均衡控制方法、装置和设备 | |
CN110687458A (zh) | 终端电池电量确定方法及装置 | |
CN110472315A (zh) | 基于ersa算法的结构参数识别方法 | |
CN110618029A (zh) | 极限载荷极短时间工况下滚珠丝杠副疲劳弹性寿命测试方法 | |
CN112100806A (zh) | 一种表征复合材料本构曲线的获取方法 | |
CN111487011B (zh) | 压力芯片的烘烤方法及设备、可读存储介质 | |
CN114182713B (zh) | 一种基于含水率的承载力预测方法 | |
CN110147643B (zh) | 车钩钩体剩余寿命确定方法和装置 | |
CN110580391B (zh) | 一种柔性结构的基频模态测量方法 | |
CN110706761A (zh) | 一种预估金属材料的疲劳强度退化的方法 | |
CN109856794B (zh) | 一种强过载环境下透镜安装结构最佳参数设计方法 | |
CN115628998B (zh) | 节点疲劳试验的可靠性确定方法、装置及存储介质 | |
CN108870677A (zh) | 控制盒的电流检测方法、装置及计算机可读存储介质 | |
CN109241614B (zh) | 一种纤维与土界面强度参数确定方法 | |
CN116577716B (zh) | 一种电流传感器振动特性测试方法、相关设备及相关系统 | |
CN111595542B (zh) | 一种航天器动力学试验数据分析系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |