CN112854513B - 一种粘滞阻尼器力学性能系数及滑移识别方法 - Google Patents

一种粘滞阻尼器力学性能系数及滑移识别方法 Download PDF

Info

Publication number
CN112854513B
CN112854513B CN202110043789.9A CN202110043789A CN112854513B CN 112854513 B CN112854513 B CN 112854513B CN 202110043789 A CN202110043789 A CN 202110043789A CN 112854513 B CN112854513 B CN 112854513B
Authority
CN
China
Prior art keywords
viscous damper
hysteresis loop
slip
rigidity
damping
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
Application number
CN202110043789.9A
Other languages
English (en)
Other versions
CN112854513A (zh
Inventor
李照众
靳慧
陶雷
龚亮
冯岩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Jiangsu Engineering Co Ltd of China Railway 24th Bureau Group Co Ltd
Original Assignee
Southeast University
Jiangsu Engineering Co Ltd of China Railway 24th Bureau Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University, Jiangsu Engineering Co Ltd of China Railway 24th Bureau Group Co Ltd filed Critical Southeast University
Priority to CN202110043789.9A priority Critical patent/CN112854513B/zh
Publication of CN112854513A publication Critical patent/CN112854513A/zh
Application granted granted Critical
Publication of CN112854513B publication Critical patent/CN112854513B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E04—BUILDING
    • E04B—GENERAL BUILDING CONSTRUCTIONS; WALLS, e.g. PARTITIONS; ROOFS; FLOORS; CEILINGS; INSULATION OR OTHER PROTECTION OF BUILDINGS
    • E04B1/00—Constructions in general; Structures which are not restricted either to walls, e.g. partitions, or floors or ceilings or roofs
    • E04B1/62—Insulation or other protection; Elements or use of specified material therefor
    • E04B1/92—Protection against other undesired influences or dangers
    • E04B1/98—Protection against other undesired influences or dangers against vibrations or shocks; against mechanical destruction, e.g. by air-raids
    • E—FIXED CONSTRUCTIONS
    • E04—BUILDING
    • E04H—BUILDINGS OR LIKE STRUCTURES FOR PARTICULAR PURPOSES; SWIMMING OR SPLASH BATHS OR POOLS; MASTS; FENCING; TENTS OR CANOPIES, IN GENERAL
    • E04H9/00—Buildings, groups of buildings or shelters adapted to withstand or provide protection against abnormal external influences, e.g. war-like action, earthquake or extreme climate
    • E04H9/02—Buildings, groups of buildings or shelters adapted to withstand or provide protection against abnormal external influences, e.g. war-like action, earthquake or extreme climate withstanding earthquake or sinking of ground

Abstract

本发明提供一种粘滞阻尼器力学性能系数及滑移识别方法，包括以下步骤：首先将粘滞阻尼器试验滞回曲线离散为单个滞回环并连接首尾使其闭环得到单个闭合滞回环；由单个闭合滞回环卸载段得到卸载刚度/加载刚度的粗略计算结果，并确定卸载刚度/加载刚度的取值范围；再计算精确阻尼系数和阻尼指数取值范围；再确定粘滞阻尼器的取值范围；输入基于粘滞阻尼器耗能能力的目标函数，并进行基于NM单纯形法的参数识别；得到粘滞阻尼器具有初始滑移的恢复力模型；本发明考虑了由于试验机造成的误差同时为粘滞阻尼器力学性能参数及滑移识别提供了一种更加合理可靠方法，使其能够应用于粘滞阻尼器的参数分析和非线性分析。

Description

一种粘滞阻尼器力学性能系数及滑移识别方法
技术领域
本发明涉及土木工程领域，具体是一种粘滞阻尼器力学性能系数及滑移识别方法。
背景技术
粘滞流体阻尼器是一种减振装置，一般由缸体、活塞、阻尼孔(或间隙、或两者兼有)、粘滞流体阻尼材料和导杆等部分组成，活塞在缸筒内作往复运动.活塞上有适量小孔成为阻尼孔(和活塞与缸筒间配合间隙)，缸筒内装满粘滞流体阻尼材料。当活塞与缸筒之间发生相对运动时，由于活塞前后的压力差使流体阻尼材料从阻尼孔(或间隙)中通过，从而产生阻尼力，达到耗能的目的。流体阻尼器对结构进行振动控制的机理是将结构振动的部分能量通过阻尼器中粘滞流体阻尼材料的粘滞耗能耗散掉，达到减小结构振动(地震或风振)反应的目的。
在JG/T 209-2012《建筑消能阻尼器》中，粘滞阻尼器的阻尼系数和阻尼参数通常通过拟合的方法进行识别和分析，但拟合中同时存在两个未知数造成了该两个参数的拟合可调整范围较大、一旦出现较大误差无法解释。同时由于试验机的误差可能存在一定的滑移。以上原因均造成得到的参数无法精确反应粘滞阻尼器的恢复力性能。
发明内容
本发明的目的是为解决上述技术问题，提供一种粘滞阻尼器力学性能系数及滑移识别方法。
为达到上述目的，本发明的一种粘滞阻尼器力学性能系数及滑移识别方法，包括以下步骤：
S1：将粘滞阻尼器试验滞回曲线离散为单个滞回环并连接首尾使其闭环得到单个闭合滞回环；
S2：由步骤S1中得到的单个闭合滞回环卸载段得到卸载刚度/加载刚度的粗略计算结果，然后根据粗略计算结果确定卸载刚度/加载刚度的取值范围；
S3：计算出粘滞阻尼器最大阻尼与速度相关的参数：阻尼系数和阻尼指数，再计算精确阻尼系数和阻尼指数取值范围；
S4：在步骤S1中得到的单个闭合滞回环水平段开始和结束段分别取一个点，得到两者差值，根据两者的差值确定粘滞阻尼器的取值范围；
S5：输入基于粘滞阻尼器耗能能力的目标函数，并结合S1-S4得到的参数取值范围或初始值进行基于NM单纯形法或Genetic Algorithm的参数识别；
S6：得到精确的粘滞阻尼器的阻尼系数、阻尼指数和滑移的识别结果，得到粘滞阻尼器具有初始滑移的恢复力模型；
S7：通过S6中得到粘滞阻尼器恢复力模型进行阻尼器的参数分析和非线性分析。
进一步地，步骤S2中，根据S1中得到的单个闭合滞回环卸载段斜率得到卸载刚度/加载刚度的粗略计算结果Kr，确定粘滞阻尼器卸载刚度K的取值范围为[0.9Kr,1.1Kr]。
进一步地，步骤S3的具体方法是，根据S1中得到的单个闭合滞回环曲线，根据JG/T209-2012《建筑消能阻尼器》计算得到粘滞阻尼器的粗略的阻尼系数Cr和阻尼指数αr，由此确定精确的阻尼系数C0和阻尼指数α的取值范围分别为[0.95Cr，1.05Cr]和[0.95αr,1.05αr]。
进一步地，步骤S4的具体方法是，在S1得到的单个闭合滞回环曲线滑移段取其粗略的首尾位移数值，求得两个点的插值d；由此确定粘滞阻尼器滞回曲线真实的滑移值D的取值范围为[0.95d,1.05d]。
进一步地，步骤S5中，NM单纯形法/Genetic Algorithm的目标函数为：
Jobjective＝|∫C0|[Asin(2πft)]’|αdx-Sslip-S|/S
其中S为试验滞回环围成的面积，f为粘滞阻尼器试验过程中加载正弦波的频率，Sslip为由于滑移减小的面积，由滑移段积分得到，t为正弦波的加载试件，x＝Asin(2πft)；
根据S1-S4中得到的相关参数取值范围或初始值，调用NM单纯形法或GeneticAlgorithm进行参数识别，并输入S5中NM单纯形法/Genetic Algorithm的目标函数，优化得到卸载刚度K、最大位移A、阻尼系数C0、阻尼指数α和滑移值D。
本发明具有以下有益效果：
本发明考虑了由于试验机造成的误差同时为粘滞阻尼器力学性能参数及滑移识别提供了一种更加合理可靠方法，使其能够应用于粘滞阻尼器的参数分析和非线性分析。
附图说明
图1为本发明中带滑移的粘滞阻尼器滞回曲线；
图2为本发明中NM参数识别或Genetic Algorithm方法的流程图；
图3为JG/T 209-2012《建筑消能阻尼器》确定粘滞阻尼器力学参数的方法；
图4为本发明中参数识别得到的滞回曲线与原始滞回曲线比较图；
具体实施方式
下面结合具体实施方式和附图对本发明的技术方案作进一步的介绍。
一种粘滞阻尼器力学性能系数及滑移识别方法，包括以下步骤：
S1：将粘滞阻尼器试验滞回曲线离散为单个滞回环并连接首尾使其闭环得到单个闭合滞回环；
S2：由步骤S1中得到的单个闭合滞回环卸载段得到卸载刚度/加载刚度的粗略计算结果，然后根据粗略计算结果确定卸载刚度/加载刚度的取值范围；
S3：计算出粘滞阻尼器最大阻尼与速度相关的参数：阻尼系数和阻尼指数，再计算精确阻尼系数和阻尼指数取值范围；
S4：在步骤S1中得到的单个闭合滞回环水平段开始和结束段分别取一个点，得到两者差值，根据两者的差值确定粘滞阻尼器的取值范围；
S5：输入基于粘滞阻尼器耗能能力的目标函数，并结合S1-S4得到的参数取值范围或初始值进行基于NM单纯形法或Genetic Algorithm的参数识别；
S6：得到精确的粘滞阻尼器的阻尼系数、阻尼指数和滑移的识别结果，得到粘滞阻尼器具有初始滑移的恢复力模型；
S7：通过S6中得到粘滞阻尼器恢复力模型进行阻尼器的参数分析和非线性分析。
步骤S2中，根据S1中得到的单个闭合滞回环卸载段斜率得到卸载刚度/加载刚度的粗略计算结果Kr，确定粘滞阻尼器卸载刚度K的取值范围为[0.9Kr,1.1Kr]。
步骤S3的具体方法是，根据S1中得到的单个闭合滞回环曲线，根据JG/T 209-2012《建筑消能阻尼器》计算得到粘滞阻尼器的粗略的阻尼系数Cr和阻尼指数αr，由此确定精确的阻尼系数C0和阻尼指数α的取值范围分别为[0.95Cr，1.05Cr]和[0.95αr,1.05αr]。
步骤S4的具体方法是，在S1得到的单个闭合滞回环曲线滑移段取其粗略的首尾位移数值，求得两个点的插值d；由此确定粘滞阻尼器滞回曲线真实的滑移值D的取值范围为[0.95d,1.05d]。
步骤S5中，NM单纯形法/Genetic Algorithm的目标函数为：
Jobjective＝|∫C0|[Asin(2πft)]’|αdx-Sslip-S|/S
其中S为试验滞回环围成的面积，f为粘滞阻尼器试验过程中加载正弦波的频率，Sslip为由于滑移减小的面积，由滑移段积分得到，t为正弦波的加载试件，x＝Asin(2πft)；
根据S1-S4中得到的相关参数取值范围或初始值，调用NM单纯形法或GeneticAlgorithm进行参数识别，并输入S5中NM单纯形法/Genetic Algorithm的的目标函数，优化得到卸载刚度K、最大位移A、阻尼系数C0、阻尼指数α和滑移值D。
NM单纯形法就是使用单纯形的思想，对于N维空间间可以构成N+1顶点的多面体。亦即对于拥有N个变量的函数,可以构成N+1个点。在计算开始时,先计算出每个点的函数值,查找出最大值的点(G)和最小值的点(S),再计算除G点外的其他点的中点(M,一般取平均),在GM连线方向上,计算G的反射点、扩展点、压缩点、若没有找到这些点,则在G至S的线段上计算一个收缩点，若在这些点中,存在一个最佳点,则使用这个点作为一个新的有效点,替代G点,重现构成一个新的N+1多面体。如此进行,直到找到函数的最小值点。
NM单纯形法实现简单，与其他确定性局部优化技术相比，具有优化速度快，计算量小，抗干扰能力强等优点，搜索过程只计算目标函数值，无需任何倒数信息。
Genetic Algorithm(GA)是一种模拟自然选择和遗传机制的不依靠梯度信息的智能全局优化算法，其操作简单，优化不受限制性条件的约束，具有较强的鲁棒性，对初值不敏感，对于一些大型、复杂非线性系统的全局优化问题表现出良好的性能。GeneticAlgorithm通过模拟自然界生物适者生存、优胜劣汰的遗传法则，利用选择、交叉、变异及倒位等一系列运算将适应度搞的优良个体遗传给下一代，产生新的个体。基本步骤如下：
(1)确定编码的方案；(2)解码；(3)选择适应度函数(4)确定选择策略(5)选择算子(6)确定算法的终止准则。
在MATLAB中，可直接在工作区中输入gatool调用GA工具箱、fminsearch工具箱或直接通过编程的方法，来实现以上的操作。
如图4所示，通过S6中得到粘滞阻尼器恢复力模型进行阻尼器的参数分析和非线性分析。
以上所述的具体实施操作方法，对本发明的技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述为本发明具体实施方式，凡在本发明的精神和原则之内，所做的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种粘滞阻尼器力学性能系数及滑移识别方法，其特征在于，包括以下步骤：
S1：将粘滞阻尼器试验滞回曲线离散为单个滞回环并连接首尾使其闭环得到单个闭合滞回环；
S2：由步骤S1中得到的单个闭合滞回环卸载段得到卸载刚度/加载刚度的粗略计算结果，然后根据粗略计算结果确定卸载刚度/加载刚度的取值范围；
S3：计算出粘滞阻尼器最大阻尼与速度相关的参数：阻尼系数和阻尼指数，再计算精确阻尼系数和阻尼指数取值范围；
S4：在步骤S1中得到的单个闭合滞回环水平段开始和结束段分别取一个点，得到两者差值，根据两者的差值确定粘滞阻尼器的取值范围；
S5：输入基于粘滞阻尼器耗能能力的目标函数，并结合S1-S4得到的参数取值范围或初始值进行基于NM单纯形法或Genetic Algorithm的参数识别；
S6：得到精确的粘滞阻尼器的阻尼系数、阻尼指数和滑移的识别结果，得到粘滞阻尼器具有初始滑移的恢复力模型；
S7：通过S6中得到粘滞阻尼器恢复力模型进行阻尼器的参数分析和非线性分析；
步骤S4的具体方法是，在S1得到的单个闭合滞回环曲线滑移段取其粗略的首尾位移数值，求得两个点的差值d；由此确定粘滞阻尼器滞回曲线真实的滑移值D的取值范围为[0.95d,1.05d]；
步骤S5中，NM单纯形法/Genetic Algorithm的目标函数为：
Jobjective＝|∫C0|[Asin(2πft)]’|αdx-Sslip-S|/S
其中S为试验滞回环围成的面积，f为粘滞阻尼器试验过程中加载正弦波的频率，Sslip为由于滑移减小的面积，由滑移段积分得到，t为正弦波的加载试件，x＝Asin(2πft)；
根据S1-S4中得到的相关参数取值范围或初始值，调用NM单纯形法或GeneticAlgorithm进行参数识别，并输入S5中NM单纯形法/Genetic Algorithm的目标函数，优化得到卸载刚度K、最大位移A、阻尼系数C0、阻尼指数α和滑移值D。
2.根据权利要求1所述的一种粘滞阻尼器力学性能系数及滑移识别方法，其特征在于：步骤S2中，根据S1中得到的单个闭合滞回环卸载段斜率得到卸载刚度/加载刚度的粗略计算结果Kr，确定粘滞阻尼器卸载刚度K的取值范围为[0.9Kr,1.1Kr]。
3.根据权利要求1所述的一种粘滞阻尼器力学性能系数及滑移识别方法，其特征在于：步骤S3的具体方法是，根据S1中得到的单个闭合滞回环曲线，根据JG/T 209-2012《建筑消能阻尼器》计算得到粘滞阻尼器的粗略的阻尼系数Cr和阻尼指数αr，由此确定精确的阻尼系数C0和阻尼指数α的取值范围分别为[0.95Cr，1.05Cr]和[0.95αr,1.05αr]。
CN202110043789.9A 2021-01-13 2021-01-13 一种粘滞阻尼器力学性能系数及滑移识别方法 Active CN112854513B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110043789.9A CN112854513B (zh) 2021-01-13 2021-01-13 一种粘滞阻尼器力学性能系数及滑移识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110043789.9A CN112854513B (zh) 2021-01-13 2021-01-13 一种粘滞阻尼器力学性能系数及滑移识别方法

Publications (2)

Publication Number Publication Date
CN112854513A CN112854513A (zh) 2021-05-28
CN112854513B true CN112854513B (zh) 2021-11-30

Family

ID=76003509

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110043789.9A Active CN112854513B (zh) 2021-01-13 2021-01-13 一种粘滞阻尼器力学性能系数及滑移识别方法

Country Status (1)

Country Link
CN (1) CN112854513B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117251952A (zh) * 2023-09-08 2023-12-19 海南大学 基于多水准分级屈服阻尼器的减震结构的优化设计方法
CN117216846B (zh) * 2023-09-12 2024-04-19 华南理工大学 钢筋混凝土构件滞回曲线预测方法、系统、设备及介质

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002070943A (ja) * 2000-08-30 2002-03-08 Toyo Tire & Rubber Co Ltd 免震滑り支承装置
CN105865956B (zh) * 2016-04-01 2019-08-13 中国电力科学研究院有限公司 管母支撑伸缩线夹滞回特性模型预测装置及方法
CN110273555B (zh) * 2019-04-12 2021-06-15 同济大学建筑设计研究院(集团)有限公司 组合耗能减震结构设计方法及组合耗能减震高层建筑
CN111351622B (zh) * 2020-03-11 2020-11-27 东南大学 一种基于遗传算法的铅芯橡胶支座恢复力模型参数识别方法
CN111980195A (zh) * 2020-09-18 2020-11-24 苏州海德新材料科技股份有限公司 一种弹性滑板隔震支座

Also Published As

Publication number Publication date
CN112854513A (zh) 2021-05-28

Similar Documents

Publication Publication Date Title
CN112854513B (zh) 一种粘滞阻尼器力学性能系数及滑移识别方法
US7480750B2 (en) Optimization of buffer pool sizes for data storage
CN103425743A (zh) 基于贝叶斯神经网络算法的蒸汽管网预测系统
CN113654756B (zh) 一种海上浮式风机主动式实时混合模型试验方法
CN105354800A (zh) 基于图像结构的粒子群优化非凸压缩感知图像重构方法
CN110020721A (zh) 一种基于参数压缩的目标检测深度学习网络优化方法
CN105046203A (zh) 基于夹角dtw距离的卫星遥测数据自适应层次聚类方法
CN103207804B (zh) 基于集群作业日志的MapReduce负载模拟方法
Wang et al. An Improved Grey Prediction Model Based on Matrix Representations of the Optimized Initial Value.
CN114116778A (zh) 一种数据库查询优化方法
CN113821863B (zh) 一种桩基竖向极限承载力预测方法
JP2005050063A (ja) 流体の数値計算方法及び装置及び制御プログラム及び記憶媒体
CN103235879A (zh) 基于小生境主从式并行遗传算法的双相介质参数反演方法
Jones et al. A logarithmic-time solution to the point location problem for closed-form linear MPC
Pries A short guide to p-torsion of abelian varieties in characteristic p
CN106570220B (zh) 基于遗传算法的液压悬置参数辨识方法
CN114861984A (zh) 用于预测高含co2凝析气藏油环体积的方法及处理器
CN114799561A (zh) 基于案例推理和模型推理的激光切割工艺参数优化方法
CN113158297A (zh) 一种水工弧形钢闸门参数荷载识别方法
CN111159935A (zh) 基于lhs的bp神经网络参数标定方法
CN115146389B (zh) 一种永磁磁浮列车动力学特征建模方法
CN113673167B (zh) 基于高斯过程回归的风电机组输出功率的建模方法和系统
CN106383779A (zh) 一种数据库适配压力的测试方法及装置
CN115660488A (zh) 基于混沌粒子群遗传算法的区域水资源承载力评价方法
CN115563438A (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