CN114528618B - 一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法 - Google Patents
一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法 Download PDFInfo
- Publication number
- CN114528618B CN114528618B CN202210077429.5A CN202210077429A CN114528618B CN 114528618 B CN114528618 B CN 114528618B CN 202210077429 A CN202210077429 A CN 202210077429A CN 114528618 B CN114528618 B CN 114528618B
- Authority
- CN
- China
- Prior art keywords
- tld
- tower
- conjoined
- wind
- kth
- 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
- 238000000034 method Methods 0.000 title claims abstract description 71
- 239000007788 liquid Substances 0.000 title claims abstract description 37
- 238000013016 damping Methods 0.000 claims abstract description 70
- 230000004044 response Effects 0.000 claims abstract description 41
- 238000012360 testing method Methods 0.000 claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 16
- 230000008878 coupling Effects 0.000 claims abstract description 13
- 238000010168 coupling process Methods 0.000 claims abstract description 13
- 238000005859 coupling reaction Methods 0.000 claims abstract description 13
- 230000001360 synchronised effect Effects 0.000 claims abstract description 13
- 230000000694 effects Effects 0.000 claims abstract description 12
- 238000003860 storage Methods 0.000 claims abstract description 10
- 230000009471 action Effects 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 24
- 230000001133 acceleration Effects 0.000 claims description 22
- 239000002245 particle Substances 0.000 claims description 18
- 238000013461 design Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 12
- 239000012530 fluid Substances 0.000 claims description 10
- 230000003595 spectral effect Effects 0.000 claims description 9
- 230000005284 excitation Effects 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000009434 installation Methods 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 238000005316 response function Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 5
- 230000009977 dual effect Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/30—Wind power
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Optimization (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Architecture (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
本发明公开了一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法,包括以下步骤:S1、进行大气边界层流场模拟,通过风洞试验获取连体塔楼的风荷载时程;S2、建立连体塔楼的运动方程;S3、根据连体塔楼的动力特性计算每个TLD达到最优控制效果时的频率和阻尼比,设计TLD的尺寸和储水高度;S4、通过对比TLD数值模型的阻尼比和最优阻尼比的差别,确定内部阻尼构件的方案;S5、建立多个TLD同步振动的数值模型;S6、建立连体塔楼在TLD控制下的耦合运动方程;S7、根据耦合运动方程计算不同重现期风荷载作用下,连体塔楼与多个TLD系统的风致振动响应。本发明不再局限于计算单体塔楼和单个TLD的振动,而是可以计算连体多塔结构和多个TLD的耦合风致振动。
Description
技术领域
本发明属于结构振动控制技术领域,具体涉及一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法。
背景技术
随着城市发展和技术进步,超高层建筑结构的形式日趋多样化,大型连体多塔结构层出不穷。由于连体塔楼的外形多变且结构间往往存在气动干扰现象,此类型超高建筑的风效应成为控制其安全性和舒适性的主要因素之一。
调谐液体阻尼器(TLD)是高层建筑风振控制中一种非常有效的被动动力减振器,可以显著降低结构顶部加速度,提高居住者的舒适性。TLD减振的基本原理是利用液体在晃荡过程中吸收和耗散振动能量,为结构提供附加阻尼。已有的风荷载作用下结构-TLD耦合振动的实现方法主要有:1)通过风洞试验获取作用在建筑上的风荷载时程,并将TLD简化为调谐质量阻尼器(TMD),然后进行结构风振计算便可以确定建筑受控后的风振响应。2)建立单栋高层建筑的理论模型和单个TLD数值模型,通过结构响应和TLD控制力的实时数据交互实现系统的耦合振动。
对于以上方法,方法1)将TLD等效为TMD来简化非线性的液体晃荡问题,但该做法对于形状不规则或者设置阻尼构件的复杂TLD可能会造成较大的误差。方法2)只能考虑基本的单体建筑结构和单个TLD的振动,然而连体塔楼是相互关联的,将塔楼隔开单独分析不能全面准确地反映其动力特性,从而造成计算结果失真。
发明内容
本发明的主要目的在于克服现有技术的缺点与不足,提出一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法。
为了达到上述目的,本发明采用以下技术方案:
一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法,包括以下步骤:
S1、进行大气边界层流场模拟,通过风洞试验获取连体塔楼的风荷载时程;
S2、由连体塔楼的集中质量模型组装质量矩阵和刚度矩阵,并以此建立连体塔楼的运动方程;
S3、根据连体塔楼的动力特性计算每个TLD达到最优控制效果时的频率和阻尼比,然后针对塔楼平面构造设计TLD的尺寸和储水高度;
S4、通过对比TLD数值模型的阻尼比和最优阻尼比的差别,确定内部阻尼构件的方案;
S5、在确定每个TLD的设计方案后,建立多个TLD同步振动的数值模型;
S6、根据多个TLD在连体塔楼的楼层位置,建立连体塔楼在TLD控制下的耦合运动方程;
S7、根据耦合运动方程计算不同重现期风荷载作用下,连体塔楼与多个TLD系统的风致振动响应。
进一步的,步骤S1具体为:
基于连体塔楼周围的地貌环境判断地面粗糙度类型,根据建筑风洞试验相关标准进行大气边界层流场模拟;
由塔楼高度确定合适的几何缩尺比λL并根据结构外形制作刚性试验缩尺模型;
假定将连体塔楼的各楼层简化为集中质量模型,则通过风洞试验获得各质点处的合力,然后根据相似理论转化为原型的风荷载,具体公式为:
其中,λF表示气动力缩尺比,λL表示几何缩尺比,λV风速比。
进一步的,步骤S2具体为:
假定连体部分与塔楼的连接形式满足强连接的要求,将连体所在楼层看作一个集中质量;各质点的质量和侧向刚度组装为总体质量矩阵M和刚度矩阵K,组装方法由每个质点的受力状态确定,然后阻尼矩阵C采用瑞利阻尼形式,具体公式为:
C=a0M+a1K (2)
其中,a0和a1两个常数通过给定结构任意两个模态的动力特性来确定,具体公式为:
其中,ωn和ωm为结构第n和m阶模态频率,ζn和ζm为结构第n和m阶模态阻尼比;
建立连体塔楼的运动方程,具体为:
将不同重现期的风荷载时程Fw代入公式(4),通过数值积分方法计算得到连体塔楼的风致振动响应。
进一步的,步骤S3具体为:
根据连体塔楼的基阶固有频率fs和第k个TLD的质量比μk,计算每个TLD达到最优控制效果时的晃动频率fk和阻尼比ζk,具体公式为:
TLD采用矩形水箱的形式,考察TLD安装楼层的平面构造,在空间允许情况下设计第k个TLD长度为Lk,然后由基阶频率理论公式计算得到TLD的储水高度hk,具体公式为:
其中,g为重力加速度;
储水高度hk需满足TLD所在楼层的层高限制;
通过连体塔楼的模态质量ms计算第k个TLD的宽度Wk,具体公式为:
其中,ρ表示水的密度,βk表示第k个TLD的质量参与系数。
进一步的,若TLD的设计尺寸不满足楼面的空间限制,则调整长度Lk重新进行设计。
进一步的,步骤S4具体为:
在TLD内部设置阻尼构件来增大阻尼比,提高TLD减振控制能力;
通过改变阻尼构件的数量、尺寸以及安装方位得到不同的阻尼比;
根据步骤S3得到TLD的尺寸和最优阻尼比情况,对每个TLD内部阻尼构件的配置进行设计,然后建立相应的数值模型,并计算其在有限带宽白噪声激励下的液面晃动响应;
通过引入流体体积分数α来表示模型单元内液体的成分,从而实时监测TLD液面的变化,具体公式为:
其中,u为流体速度矢量,当α=0.5时表示自由液面;
采用贝叶斯谱密度法或曲线拟合法由第k个TLD的液面响应获取其基阶阻尼比大小,具体公式为:
Sk(f)=|hk(f)|S0 (10)
其中,Sk(f)为第k个TLD的基阶响应功率谱密度,hk(f)为第k个TLD的频率响应函数,S0为白噪声激励的功率谱密度值。
进一步的,贝叶斯谱密度法基于Sk(f)在特定频率区段为Chi-square分布的假设,通过令模态参数的后验概率密度函数最大化得到第k个TLD的基阶阻尼比ζk *;
曲线拟合法通过对Sk(f)进行最小二乘法拟合得到第k个TLD的基阶阻尼比ζk *;
将第k个TLD数值模型的阻尼比ζk *与最优阻尼比ζk进行对比,如果ζk *与ζk的相差不超过5%,则表明第k个TLD的构件方案符合要求;
如果ζk *与ζk的相差超过5%,则改变构件方案并重新进行数值模拟,直至数值模型的ζk *与ζk相差不超过5%。
进一步的,步骤S5具体为:
基于步骤S3和步骤S4确定每个TLD的设计方案后,建立它们的同步振动数值模型;
同步振动模型将所有TLD建在同一个网格文件内,同时为了使得TLD晃动过程中互不影响,将TLD模型设置为互不接触或者接触面设置为壁面边界,而且每个TLD都设置为不同的流体域。
进一步的,步骤S6具体为:
假设TLD在连体塔楼的某一楼层,按照集中质量模型,相当于在该楼层质点施加控制力,因此得到连体塔楼在TLD控制下的运动方程,具体公式为:
其中,Fk为第k个TLD的控制力,Hk为反映TLD所在楼层的列向量;Fk由TLD所有壁面边界的压强积分得到,具体公式为:
其中,nk为第k个TLD的壁面边界数量,Ski、Pki和θki分别为第k个TLD中第i个边界、边界处压强以及边界与振动方向的夹角,ds代表边界微元。
进一步的,步骤S7具体为:
由于步骤S1得到的不同重现期下连体塔楼风荷载是各个质点的荷载时程数据,对于每一个时刻都将对应的各质点荷载和所有TLD控制力代入公式(11),并通过数值积分方法计算结构的响应,获取每个TLD所在楼层的结构加速度响应,具体公式为:
通过连体塔楼加速度和TLD控制力的实时数据交互实现了耦合振动的过程,从而得到不同重现期风荷载作用下整个系统的风致振动响应。
本发明与现有技术相比,具有如下优点和有益效果:
1、本发明通过建立多个不同尺寸和阻尼构件配置的TLD数值模型,能够充分反映TLD液体晃荡过程中的非线性特征,有效避免了TMD简化模型造成的计算误差。
2、本发明通过集中质量法建立具有强连接的连体塔楼理论模型,可以较为准确的反映其动力特性,有效解决现有方法只能考虑单体结构的问题。
3、本发明具有普适性,适用于形状不规则或者设置阻尼构件的复杂TLD,也适用于安装在任意塔楼和任意楼层的多个TLD。
附图说明
图1是本发明方法的流程图;
图2是本发明实施例中的连体塔楼风洞试验示意图;
图3是本发明实施例中的连体塔楼集中模型示意图;
图4是本发明实施例中的双TLD设计方案俯视图;
图5是本发明实施例中的双TLD数值模型示意图;
图6a是本发明实施例中塔楼1屋顶层的加速度响应示意图;
图6b是本发明实施例中塔楼2屋顶层的加速度响应示意图;
图7a是本发明实施例中塔楼1屋顶层的加速度响应功率谱示意图;
图7b是本发明实施例中塔楼2屋顶层的加速度响应功率谱示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例
如图1所示,本发明,一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法,包括以下步骤:
S1、根据相关标准进行大气边界层流场模拟,然后通过风洞试验获取连体塔楼的风荷载时程;具体为:
基于连体塔楼周围的地貌环境判断地面粗糙度类型,根据建筑风洞试验相关标准进行大气边界层流场模拟。由塔楼高度确定合适的几何缩尺比λL,并根据结构外形制作刚性试验缩尺模型。假定连体塔楼的各楼层可简化为集中质量模型,则通过风洞试验可以获得各质点处的合力,然后根据相似理论转化为原型的风荷载,具体公式为:
其中,λF表示气动力缩尺比,λL表示几何缩尺比,λV风速比。
在本实施例中,连体塔楼由塔楼1和塔楼2构成,如图2所示,对塔楼进行多点同步测压风洞试验得到风荷载时程Fw。
S2、由连体塔楼的集中质量模型组装质量矩阵和刚度矩阵,并以此建立塔楼的运动方程;具体为:
假定连体部分与塔楼的连接形式满足强连接的要求,也就是说连接部分将两侧塔楼连成一个整体来协调受力和变形,此时就可以将连体所在楼层看作一个集中质量。各质点的质量和侧向刚度组装为总体质量矩阵M和刚度矩阵K,组装方法由每个质点的受力状态确定。然后阻尼矩阵C采用瑞利阻尼形式,具体公式为:
C=a0M+a1K (2)
其中,a0和a1两个常数通过给定结构任意两个模态的动力特性来确定,具体公式为:
其中,ωn和ωm为结构第n和m阶模态频率,ζn和ζm为结构第n和m阶模态阻尼比。然后建立连体塔楼的运动方程,具体公式为:
在本实施例中,连体塔楼的集中质量模型如图3所示,1至6层为连体层,7-11层分离,12至14层为连体层,右边的15至55层为塔楼1,左边的15至65层为塔楼2。连体塔楼一共有111个质点,因此位移xs是111个元素组成的列向量,质量矩阵M和刚度矩阵K是111行×111列的矩阵。阻尼矩阵C通过设置第1和2阶的模态阻尼比为0.02计算得到。
S3、根据连体塔楼的动力特性计算每个TLD达到最优控制效果时的频率和阻尼比,然后针对塔楼平面构造设计TLD的尺寸和储水高度;具体为:
根据连体塔楼的基阶固有频率fs和第k个TLD的质量比μk,然后计算每个TLD达到最优控制效果时的晃动频率fk和阻尼比ζk,具体公式为:
TLD采用矩形水箱的形式,考察TLD安装楼层的平面构造,在空间允许情况下设计第k个TLD长度为Lk,然后由基阶频率理论公式计算得到TLD的储水高度hk,具体公式为:
其中,g为重力加速度。储水高度hk需要满足TLD所在楼层的层高限制。然后通过连体塔楼的模态质量ms计算第k个TLD的宽度Wk,具体公式为:
其中,ρ表示水的密度,βk表示第k个TLD的质量参与系数。若TLD的尺寸不满足楼面的空间限制,需要调整长度Lk重新进行设计。
在本实施例中,设计了两个不同尺寸的TLD1和TLD2来对连体塔楼进行风振控制,其中TLD1放置在塔楼1屋顶层,基本尺寸为:L1=15.3m、W1=9.85m和h1=3.08m。TLD2放置在塔楼2屋顶层,基本尺寸为:L2=18.2m、W2=3.6m和h2=4.6m。这两个TLD方案均满足所在楼面的空间限制。
S4、通过对比TLD数值模型的阻尼比和最优阻尼比的差别,确定内部阻尼构件的方案;具体为:
在TLD内部设置阻尼构件能有效增大阻尼比,提高TLD减振控制能力。常用的构件类型包括:阻尼网、格栅、挡板、桨柱等,通过改变构件的数量、尺寸和安装方位可以得到不同的阻尼比。根据步骤S3得到TLD的基本尺寸和最优阻尼比情况后,对每个TLD内部阻尼构件的配置进行设计,然后建立相应的数值模型,并计算其在有限带宽白噪声激励下的液面晃动响应。通过引入流体体积分数α来表示模型单元内液体的成分,从而实时监测TLD液面的变化,具体公式为:
其中,u为流体速度矢量,当α=0.5时表示自由液面。
然后采用贝叶斯谱密度法或曲线拟合法由第k个TLD的液面响应获取其基阶阻尼比大小,具体公式为:
Sk(f)=|hk(f)|S0 (10)
其中,Sk(f)为第k个TLD的基阶响应功率谱密度,hk(f)为第k个TLD的频率响应函数,S0为白噪声激励的功率谱密度值。贝叶斯谱密度法基于Sk(f)在特定频率区段近似为Chi-square分布的假设,通过令模态参数的后验概率密度函数最大化得到第k个TLD的基阶阻尼比ζk *。曲线拟合法通过对Sk(f)进行最小二乘法拟合得到第k个TLD的基阶阻尼比ζk *。
将第k个TLD数值模型的阻尼比ζk *与最优阻尼比ζk进行对比,如果ζk *与ζk的相差不超过5%,则表明第k个TLD的构件方案符合要求。如果ζk *与ζk的相差超过5%,则改变构件方案并重新进行数值模拟,直至数值模型的ζk *与ζk相差不超过5%。
在本实施例中,在TLD1的1/4L1、1/2L1、3/4L1位置分别安装一列等间距桨柱,每列3个,桨柱构件的宽度ap1为0.6m。在TLD2的1/4L2、1/2L2、3/4L2位置分别安装一列等间距桨柱,每列2个,桨柱构件的宽度ap2为0.43m。TLD1和TLD2的俯视图如图4所示。这两个TLD的内部阻尼构件方案满足与最优阻尼比相差不超过5%的要求。
S5、在确定每个TLD的设计方案后,建立多个TLD同步振动的数值模型;具体为:
基于步骤S3和步骤S4确定每个TLD的设计方案后,建立它们的同步振动数值模型。与各个TLD独立振动数值模型不同的是,同步振动模型将所有TLD建在同一个网格文件内,同时为了使得TLD晃动过程中互不影响,应将TLD模型设置为互不接触或者接触面设置为壁面边界,而且每个TLD都设置为不同的流体域。
在本实施例中,TLD1和TLD2的同步振动数值模型如图5所示,两者没有直接接触,属于两个不同的流体域。
S6、根据多个TLD在连体塔楼的楼层位置,建立连体塔楼在TLD控制下的耦合运动方程;具体为:
假设TLD在连体塔楼的某一楼层,按照集中质量模型,相当于在该楼层质点施加控制力,因此得到连体塔楼在TLD控制下的运动方程,具体公式为:
其中,Fk为第k个TLD的控制力,Hk为反映TLD所在楼层的列向量。Fk可由TLD所有壁面边界的压强积分得到,具体公式为:
其中,nk为第k个TLD的壁面边界数量,Ski、Pki和θki分别为第k个TLD中第i个边界、边界处压强以及边界与振动方向的夹角,ds代表边界微元。
S7、根据耦合运动方程计算不同重现期风荷载作用下,连体塔楼与多个TLD系统的风致振动响应;具体为:
由于步骤S1得到的不同重现期下连体塔楼风荷载是各个质点的荷载时程数据,所以对于每一个时刻都将对应的各质点荷载和所有TLD控制力代入公式(11),并通过数值积分方法计算结构的响应,然后获取每个TLD所在楼层的结构加速度响应,具体公式为:
其中,为第k个TLD所在楼层的结构加速度响应,Hk T表示Hk的转置。然后将/>作为外激励施加到第k个TLD上模拟液体晃荡过程,进而获取下一时刻计算所需的TLD控制力。通过连体塔楼加速度和TLD控制力的实时数据交互实现了耦合振动的过程,从而得到不同重现期风荷载作用下整个系统的风致振动响应。
在本实施例中,基于公式(4)和公式(11)分别计算10年重现期风荷载作用下,连体双塔在受到双TLD控制之前和之后的振动响应,其中塔楼1和塔楼2的屋顶层加速度响应如图6a和6b所示,对应的加速度响应功率谱如图7a和7b所示,充分说明这两个TLD对连体双塔结构的减振效果明显。
本发明首先对连体塔楼进行风洞试验获取风荷载时程数据,然后分别通过集中质量法和有限元法建立连体塔楼的理论模型和多个TLD的同步振动数值模型,最后计算得到不同重现期风荷载作用下连体结构和多个TLD的振动响应。本发明与现有技术相比,避免了采用TLD简化模型造成的误差,并且不再局限于计算单体塔楼和单个TLD的振动,而是可以计算连体多塔结构和多个TLD的耦合风致振动。本发明具有精度高、稳定性好、适用性广的特点,而且实现形式清晰简洁,适于工程应用。
还需要说明的是,在本说明书中,诸如术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其他实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (8)
1.一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法,其特征在于,包括以下步骤:
S1、进行大气边界层流场模拟,通过风洞试验获取连体塔楼的风荷载时程;
S2、由连体塔楼的集中质量模型组装质量矩阵和刚度矩阵,并以此建立连体塔楼的运动方程;
S3、根据连体塔楼的动力特性计算每个TLD达到最优控制效果时的频率和阻尼比,然后针对塔楼平面构造设计TLD的尺寸和储水高度;具体为:
根据连体塔楼的基阶固有频率fs和第k个TLD的质量比μk,计算每个TLD达到最优控制效果时的晃动频率fk和阻尼比ζk,具体公式为:
TLD采用矩形水箱的形式,考察TLD安装楼层的平面构造,在空间允许情况下设计第k个TLD长度为Lk,然后由基阶频率理论公式计算得到TLD的储水高度hk,具体公式为:
其中,g为重力加速度;
储水高度hk需满足TLD所在楼层的层高限制;
通过连体塔楼的模态质量ms计算第k个TLD的宽度Wk,具体公式为:
其中,ρ表示水的密度,βk表示第k个TLD的质量参与系数;
S4、通过对比TLD数值模型的阻尼比和最优阻尼比的差别,确定内部阻尼构件的方案;具体为:
在TLD内部设置阻尼构件来增大阻尼比,提高TLD减振控制能力;
通过改变阻尼构件的数量、尺寸以及安装方位得到不同的阻尼比;
根据步骤S3得到TLD的尺寸和最优阻尼比情况,对每个TLD内部阻尼构件的配置进行设计,然后建立相应的数值模型,并计算其在有限带宽白噪声激励下的液面晃动响应;
通过引入流体体积分数α来表示模型单元内液体的成分,从而实时监测TLD液面的变化,具体公式为:
其中,u为流体速度矢量,当α=0.5时表示自由液面;
采用贝叶斯谱密度法或曲线拟合法由第k个TLD的液面响应获取其基阶阻尼比大小,具体公式为:
Sk(f)=|hk(f)|S0 (10)
其中,Sk(f)为第k个TLD的基阶响应功率谱密度,hk(f)为第k个TLD的频率响应函数,S0为白噪声激励的功率谱密度值;
S5、在确定每个TLD的设计方案后,建立多个TLD同步振动的数值模型;
S6、根据多个TLD在连体塔楼的楼层位置,建立连体塔楼在TLD控制下的耦合运动方程;
S7、根据耦合运动方程计算不同重现期风荷载作用下,连体塔楼与多个TLD系统的风致振动响应。
3.根据权利要求1所述的一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法,其特征在于,步骤S2具体为:
假定连体部分与塔楼的连接形式满足强连接的要求,将连体所在楼层看作一个集中质量;各质点的质量和侧向刚度组装为总体质量矩阵M和刚度矩阵K,组装方法由每个质点的受力状态确定,然后阻尼矩阵C采用瑞利阻尼形式,具体公式为:
C=a0M+a1K(2)
其中,a0和a1两个常数通过给定结构任意两个模态的动力特性来确定,具体公式为:
其中,ωn和ωm为结构第n和m阶模态频率,ζn和ζm为结构第n和m阶模态阻尼比;
建立连体塔楼的运动方程,具体为:
将不同重现期的风荷载时程Fw代入公式(4),通过数值积分方法计算得到连体塔楼的风致振动响应。
4.根据权利要求1所述的一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法,其特征在于,若TLD的设计尺寸不满足楼面的空间限制,则调整长度Lk重新进行设计。
5.根据权利要求1所述的一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法,其特征在于,贝叶斯谱密度法基于Sk(f)在特定频率区段为Chi-square分布的假设,通过令模态参数的后验概率密度函数最大化得到第k个TLD的基阶阻尼比ζk *;
曲线拟合法通过对Sk(f)进行最小二乘法拟合得到第k个TLD的基阶阻尼比ζk *;
将第k个TLD数值模型的阻尼比ζk *与最优阻尼比ζk进行对比,如果ζk *与ζk的相差不超过5%,则表明第k个TLD的构件方案符合要求;
如果ζk *与ζk的相差超过5%,则改变构件方案并重新进行数值模拟,直至数值模型的ζk *与ζk相差不超过5%。
6.根据权利要求5所述的一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法,其特征在于,步骤S5具体为:
基于步骤S3和步骤S4确定每个TLD的设计方案后,建立它们的同步振动数值模型;
同步振动模型将所有TLD建在同一个网格文件内,同时为了使得TLD晃动过程中互不影响,将TLD模型设置为互不接触或者接触面设置为壁面边界,而且每个TLD都设置为不同的流体域。
8.根据权利要求7所述的一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法,其特征在于,步骤S7具体为:
由于步骤S1得到的不同重现期下连体塔楼风荷载是各个质点的荷载时程数据,对于每一个时刻都将对应的各质点荷载和所有TLD控制力代入公式(11),并通过数值积分方法计算结构的响应,获取每个TLD所在楼层的结构加速度响应,具体公式为:
通过连体塔楼加速度和TLD控制力的实时数据交互实现了耦合振动的过程,从而得到不同重现期风荷载作用下整个系统的风致振动响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210077429.5A CN114528618B (zh) | 2022-01-24 | 2022-01-24 | 一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210077429.5A CN114528618B (zh) | 2022-01-24 | 2022-01-24 | 一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114528618A CN114528618A (zh) | 2022-05-24 |
CN114528618B true CN114528618B (zh) | 2023-06-20 |
Family
ID=81620697
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210077429.5A Active CN114528618B (zh) | 2022-01-24 | 2022-01-24 | 一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114528618B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115130225B (zh) * | 2022-08-10 | 2023-02-17 | 华南理工大学 | 一种结构风振控制的mtld系统频率设计方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106948256B (zh) * | 2017-04-26 | 2019-06-21 | 中铁大桥科学研究院有限公司 | 一种超低频液体质量调谐阻尼器及设计方法 |
CN113065208B (zh) * | 2021-03-29 | 2021-11-23 | 华南理工大学 | 调谐液体阻尼器设计方法、装置、介质和设备 |
CN113673012B (zh) * | 2021-07-29 | 2023-08-22 | 华南理工大学 | 基于风洞试验的结构-tld系统的耦合振动实现方法 |
-
2022
- 2022-01-24 CN CN202210077429.5A patent/CN114528618B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114528618A (zh) | 2022-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chan | Optimal lateral stiffness design of tall buildings of mixed steel and concrete construction | |
Yalla et al. | Optimum absorber parameters for tuned liquid column dampers | |
CN114528618B (zh) | 一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法 | |
CN113673012B (zh) | 基于风洞试验的结构-tld系统的耦合振动实现方法 | |
CN109299571B (zh) | 一种基于建筑消防水箱的调谐液体阻尼器的设计方法 | |
CN103294857B (zh) | 以频响位移为优化目标的约束阻尼板的拓扑优化方法 | |
CN108416083B (zh) | 一种高耸电视塔结构二维动力模型分析方法及系统 | |
CN111859741A (zh) | 一种加筋板壳结构阻尼层拓扑与加筋肋布局协同优化方法 | |
CN110717216A (zh) | 不规则波下带柔性气囊直升机横摇响应预报方法 | |
CN110826132A (zh) | 一种结构分散振动控制系统设计方法 | |
CN114491748B (zh) | 一种基于oc-pso的超高层建筑抗风性能设计优化方法 | |
CN115062500A (zh) | 分布式随机激励下结构振动响应分析方法 | |
CN112883328A (zh) | 群桩基础水平动力响应时域非线性计算方法 | |
CN115130225B (zh) | 一种结构风振控制的mtld系统频率设计方法 | |
CN104820762A (zh) | 含有钢管混凝土柱的高层建筑框架结构优化设计方法 | |
Wang et al. | Relationships between vibrations of main device and mechanical network in a classical vibration control system | |
Qu et al. | A reduced‐order modeling technique for tall buildings with active tuned mass damper | |
CN110704922A (zh) | 一种基于ansys软件的环状水箱tld振动分析方法 | |
CN114528652B (zh) | 一种内置桨柱调谐液体阻尼器的设计方法 | |
CN108647383A (zh) | 一种结构增强调谐质量阻尼器优化设计方法 | |
Chen | A more precise computation of along wind dynamic response analysis for tall buildings built in urban areas | |
CN114676543B (zh) | 一种用于两相邻建筑结构的黏弹性阻尼器设计方法 | |
CN116150842B (zh) | 基于iwl法计算弯扭柱螺旋梁景观塔设计风荷载的方法 | |
CN117236218A (zh) | 一种超高层风作用下超高速电梯轿厢的模拟方法 | |
CN113282995A (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 |