CN113065208A - 调谐液体阻尼器设计方法、装置、介质和设备 - Google Patents
调谐液体阻尼器设计方法、装置、介质和设备 Download PDFInfo
- Publication number
- CN113065208A CN113065208A CN202110334937.2A CN202110334937A CN113065208A CN 113065208 A CN113065208 A CN 113065208A CN 202110334937 A CN202110334937 A CN 202110334937A CN 113065208 A CN113065208 A CN 113065208A
- Authority
- CN
- China
- Prior art keywords
- tld
- damping ratio
- response
- design
- fitting
- 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.)
- Granted
Links
- 238000013461 design Methods 0.000 title claims abstract description 138
- 239000007788 liquid Substances 0.000 title claims abstract description 120
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000013016 damping Methods 0.000 claims abstract description 161
- 230000004044 response Effects 0.000 claims abstract description 108
- 230000008878 coupling Effects 0.000 claims abstract description 31
- 238000010168 coupling process Methods 0.000 claims abstract description 31
- 238000005859 coupling reaction Methods 0.000 claims abstract description 31
- 230000009471 action Effects 0.000 claims abstract description 26
- 238000012360 testing method Methods 0.000 claims abstract description 24
- 238000004458 analytical method Methods 0.000 claims abstract description 21
- 238000004088 simulation Methods 0.000 claims abstract description 17
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 38
- 230000005284 excitation Effects 0.000 claims description 22
- 238000010206 sensitivity analysis Methods 0.000 claims description 21
- 238000005265 energy consumption Methods 0.000 claims description 16
- 230000001133 acceleration Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 229920002430 Fibre-reinforced plastic Polymers 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 239000011151 fibre-reinforced plastic Substances 0.000 claims description 8
- 238000012986 modification Methods 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 6
- 230000004048 modification Effects 0.000 claims description 5
- 238000000926 separation method Methods 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 abstract description 7
- 230000035945 sensitivity Effects 0.000 abstract 1
- 230000003595 spectral effect Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 239000000463 material Substances 0.000 description 7
- 238000004519 manufacturing process Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000007797 corrosion Effects 0.000 description 3
- 238000005260 corrosion Methods 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 238000013398 bayesian method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000007769 metal material Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 description 1
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- 229910052796 boron Inorganic materials 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000003365 glass fiber Substances 0.000 description 1
- JEIPFZHSYJVQDO-UHFFFAOYSA-N iron(III) oxide Inorganic materials O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000012779 reinforcing material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000005654 stationary process Effects 0.000 description 1
- 238000011410 subtraction method Methods 0.000 description 1
- 229920003002 synthetic resin Polymers 0.000 description 1
- 239000000057 synthetic resin Substances 0.000 description 1
- 230000009466 transformation 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/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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
- Vibration Prevention Devices (AREA)
Abstract
本发明公开了调谐液体阻尼器设计方法、装置、介质和设备，首先TLD初步设计，得到TLD初步设计参数；风荷作用下结构‑TLD耦合系统的TLD阻尼比敏感性分析，并拟合TLD响应与TLD阻尼比关系曲线；基于TLD当前设计参数，进行TLD减振性能振动台试验或数值模拟分析，得到液体晃荡频率ft′，并拟合TLD阻尼比与TLD响应关系曲线；得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
判断TLD当前的ft′和
能否使受控结构响应满足控制目标；若否，则调整TLD设计参数，并且通过上述重计算液体晃荡频率ft′和阻尼比值
直到获取到的ft′和
使受控结构响应满足控制目标，得到TLD的最终设计结果。本发明使得TLD的设计更加细致和精确。
Description
技术领域
本发明属于超高层建筑风致振动控制领域，特别涉及一种调谐液体阻尼器设计方法、装置、介质和设备。
背景技术
调谐液体阻尼器(TLD)是一种常用在超高层建筑风振控制中的被动控制装置，因构造相对简单、造价低，且可兼做消防水箱，而受到越来越多的关注。对于结构-TLD耦合系统而言，质量比、液体晃荡频率、TLD阻尼比是与控制效果息息相关的三个重要设计参数；其中，质量比为TLD中液体有效质量与结构受控模态质量的比值，取值范围通常为1％～5％；液体晃荡频率与水箱长度和液体深度有关，该频率与结构受控模态频率相近时控制效果最佳；而TLD阻尼主要来自于自由液面破碎、液体粘度以及液体与水箱壁面的摩擦，但仅依靠上述几种耗能方式，TLD所能够提供的阻尼比远小于最优值，难以达到有效控制结构风致振动的效果，因此常采取在TLD内设置耗能装置如阻尼网、格栅等措施增大液体晃荡耗能。
结构-TLD耦合系统中TLD的设计通常以受控结构响应为控制目标，先估计结构所需有效阻尼，有效阻尼是质量比的函数，进而确定满足条件的质量比；然后由确定的质量比计算得到最优频率比(液体晃荡频率与结构受控模态频率的比值)与最优TLD阻尼比；最后根据最优频率比计算满足TLD所在楼面空间大小的水箱长度和液深，根据最优TLD阻尼比计算TLD内置耗能装置设计参数。
该设计方法通过等效线性化的方式将TLD中液体晃荡问题化简，从而能够快速获得结构-TLD耦合系统的重要设计参数进而确定TLD外形尺寸、液深与内置耗能装置设计参数。研究表明在风荷载作用下TLD响应呈明显非线性，且带有内置耗能装置的TLD阻尼比会随着TLD响应的变化而改变，因此采用线性化的设计方法不能准确评估TLD控制效果。
发明内容
本发明的第一目的在于克服现有技术的缺点与不足，提供一种调谐液体阻尼器设计方法，该方法能够考虑TLD响应非线性特性，使得TLD的设计更加细致和精确，基于此方法设计的TLD能够准确的评估出其控制效果。
本发明的第二目的在于提供一种调谐液体阻尼器设计装置。
本发明的第三目的在于提供一种存储介质。
本发明的第四目的在于提供一种计算设备。
本发明的目的通过下述技术方案实现：一种调谐液体阻尼器设计方法，包括步骤：
步骤1、获取设计要求的结构-TLD耦合系统的质量比μ、液体晃荡频率ft和TLD阻尼比
通过等效线性化方法对TLD进行初步设计，得到初步设计的TLD的初步设计参数；
步骤2、构建结构-TLD耦合系统运动方程，对结构-TLD耦合系统进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，并拟合TLD响应与TLD阻尼比关系曲线；
步骤3、基于TLD当前设计参数，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析，对测得的TLD响应进行解耦分离，再进行模态参数识别，得到液体晃荡频率ft′，并拟合TLD阻尼比与TLD响应关系曲线；
步骤4、基于步骤2和步骤3拟合TLD响应与TLD阻尼比关系曲线，计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
步骤5、判断步骤3中得到液体晃荡频率ft′和步骤4中计算得到的阻尼比值
能否使受控结构响应满足控制目标；
若否，则针对TLD设计参数进行修改，然后返回步骤3；
若是，则将最终获取到的TLD设计参数作为TLD的最终设计结果。
优选的，步骤2中构建的结构-TLD耦合系统运动方程为：
式中：
Ms、Cs和Ks分别为结构质量、阻尼和刚度矩阵，由设计信息可得；
mt、kt和ct分别为TLD质量、刚度和阻尼，β为质量参与系数，由TLD初步设计结果得到；
假设TLD与结构第n个自由度耦合，则X、xn和ut分别为结构位移向量、结构第n个自由度位移、TLD位移；
H为TLD位置向量；F为风荷载向量；
优选的，步骤3中，基于TLD当前设计参数，包括外形尺寸、液深以及内置耗能装置设计参数，构建TLD缩尺模型，针对上述构建的TLD缩尺模型，进行有色噪声激励下的TLD减振性能振动台试验，在得到试验结果后换算到TLD缩尺模型对应的原型，从而得到液体晃荡频率ft′；
或者步骤3中，基于TLD当前设计参数，包括外形尺寸、液深以及内置耗能装置设计参数，采用有限元软件建立TLD数值模型，根据实际试验条件在软件中设置相应初始条件、边界条件与外部激励，再通过求解流体运动方程获得流体晃动响应，从而实现数值模拟分析。
优选的，步骤2中，根据结构-TLD耦合系统运动方程进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，得到若干TLD响应对应的TLD阻尼比数据点，对这些数据点进行拟合得到TLD响应与TLD阻尼比关系曲线：
式中，
为敏感性分析得到的TLD阻尼比拟合结果，σr1为相应的TLD响应，α1、β1为拟合常数；
步骤3中，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析后，对测得的TLD响应进行解耦分离再进行模态参数识别，得到若干TLD响应对应的TLD阻尼比数据点，对这些数据点进行拟合可得TLD响应与TLD阻尼比关系曲线：
式中，
为减振性能分析得到的TLD阻尼比拟合结果，σr2为相应的TLD响应，α2、β2为拟合常数；
步骤4中，基于步骤2中拟合得到的TLD响应与TLD阻尼比关系曲线
以及步骤3中拟合得到的TLD响应与TLD阻尼比关系曲线
计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
为：
优选的，在步骤5中：将步骤5中得到液体晃荡频率ft′与液体晃荡频率ft进行比较，若两者之间的差值超过阈值F1，则判定TLD当前设计使得使受控结构响应不满足控制目标，对TLD的设计参数，包括TLD的外形尺寸长度或液深进行修改；
步骤5中，将步骤4中得到阻尼比值
与TLD阻尼比
进行比较，若两者之间的差值超过阈值F1，则判定TLD当前设计使受控结构响应不满足控制目标，对TLD的设计参数，包括TLD内置耗能装置设计参数进行修改。
更进一步的，在步骤5中：
当ft′与ft的差值超过阈值F1，且ft′小于ft，则减小TLD的长度或增大TLD的液深；
当ft′与ft的差值超过阈值F1，且ft′大于ft，则增大TLD的长度或减小TLD的液深；
当
与
的差值超过阈值F2，且
小于
则增大TLD内置耗能装置的尺寸、增加TLD内置耗能装置的数量或将TLD内置装置摆放的位置更加靠近TLD的中间位置处；
当
与
的差值超过阈值F2，且
大于
则减小TLD内置耗能装置的尺寸、减少TLD内置耗能装置的数量或将TLD内置装置摆放的位置更加远离TLD的中间位置处。
优选的，TLD内置耗能装置采用纤维增强塑料制作。
本发明的第二目的通过以下技术方案实现：一种调谐液体阻尼器设计装置，包括：
TLD初步设计模块，用于获取设计要求的结构-TLD耦合系统的质量比μ、液体晃荡频率ft和TLD阻尼比
通过等效线性化方法对TLD进行初步设计，得到初步设计的TLD的初步设计参数；
构建模块，用于结构-TLD耦合系统运动方程；
敏感分析和拟合模块，用于对结构-TLD耦合系统进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，并拟合TLD响应与TLD阻尼比关系曲线；
模态参数识别和拟合模块，在基于TLD当前设计，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析的情况下，用于对测得的TLD响应进行解耦分离再进行模态参数识别，得到液体晃荡频率ft′，并拟合TLD阻尼比与TLD响应关系曲线；
阻尼比计算模块，用于基于敏感分析和拟合模块与模态参数识别和拟合模块拟合得到TLD响应与TLD阻尼比关系曲线，计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
判定模块，用于判断模态参数识别和拟合模块得到液体晃荡频率ft′与阻尼比计算模块计算得到的阻尼比值
能否使受控结构响应满足控制目标；
修改模块，用于在判定模块判定出液体晃荡频率ft′和阻尼比值
不能使得受控结构响应满足控制目标时，则针对设计参数进行修改；修改后由模态参数识别和拟合模块以及阻尼比计算模块重新计算液体晃荡频率ft′和阻尼比值
并通过判定模块再进行判定，直到通过判定模块判定出得到的液体晃荡频率ft′和阻尼比值
使得受控结构响应满足控制目标，将最终获取到的TLD设计参数作为TLD的最终设计结果。
本发明的第三目的通过以下技术方案实现：一种存储介质，存储有程序，所述程序被处理器执行时，实现实施例所述的调谐液体阻尼器设计方法。
本发明的第四目的通过以下技术方案实现：一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，其特征在于，所述处理器执行存储器存储的程序时，实现实施例1所述的调谐液体阻尼器设计方法。
本发明相对于现有技术具有如下的优点及效果：
(1)本发明调谐液体阻尼器设计方法，在TLD初步设计基础上对结构-TLD耦合系统设计参数进行校验，为TLD设计方案的调整提供数据支撑，使得TLD的设计更加细致、精确。
(2)本发明调谐液体阻尼器设计方法，直接通过TLD减振性能振动台试验或数值模拟分析的方式获取液体晃荡频率与TLD阻尼比随TLD响应变化曲线，能够反映TLD真实晃动特性，考虑了TLD响应的非线性特征，基于此方法设计的TLD能够准确的评估出其控制效果。
(3)本发明调谐液体阻尼器设计方法，能够很好地体现了风荷载、结构响应、TLD响应以及TLD阻尼比间的相互联系，考虑了TLD阻尼比随TLD响应变化的特性，获得更准确的TLD阻尼比。
(4)本发明调谐液体阻尼器设计方法，TLD内置耗能装置采用纤维增强塑料制作，克服现有技术中TLD内置耗能装置采用金属材料(如钢材)制作，由于金属材料易被腐蚀且装置连接处存在动荷载作用下的疲劳问题，纤维增强塑料作为TLD内置耗能装置的制作材料，其具有轻质高强、防锈、耐腐蚀、工艺性优良，且总材料成本与采用钢材制作成本相差不多的优点，能够有效延长所设计TLD的使用寿命。
附图说明
图1是本发明调谐液体阻尼器设计方法的流程图。
图2是本发明中分别采用敏感性分析与减振性能分析得到的TLD响应对TLD阻尼比的影响示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
实施例1
本实施例公开了一种调谐液体阻尼器设计方法，如图1中所示，具体步骤包括：
步骤1、获取设计要求的结构-TLD耦合系统的质量比μ、液体晃荡频率ft和TLD阻尼比
通过等效线性化方法对TLD进行初步设计，得到初步设计的TLD外形尺寸、液深以及内置耗能装置设计参数。
结构-TLD耦合系统的质量比μ、液体晃荡频率ft和TLD阻尼比
是根据项目设计要求信息得到的，基于上述参数，通过等效线性化方法可以对TLD进行初步设计。
在本实施例中，TLD内置耗能装置采用纤维增强塑料制作，纤维增强塑料(FRP)是以玻璃纤维、碳纤维或硼纤维等作为增强材料，以合成树脂作基体材料的一种复合材料。FRP轻质高强、防锈、耐腐蚀、工艺性优良，且总材料成本与采用钢材制作成本相差不多，是一种适合应用于TLD内置耗能装置制作的材料。
步骤2、构建结构-TLD耦合系统运动方程，对结构-TLD耦合系统进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，并拟合TLD响应与TLD阻尼比关系曲线。
本实施例中，构建的结构-TLD耦合系统运动方程为：
式中：
Ms、Cs和Ks分别为结构质量、阻尼和刚度矩阵，由设计信息可得；
mt、kt和ct为TLD质量、刚度和阻尼，β为质量参与系数，由TLD初步设计结果得到；
假设TLD与结构第n个自由度耦合，则X、xn和ut分别为结构位移向量、结构第n个自由度位移、TLD位移；
H为TLD位置向量；F为风荷载向量；
在本步骤中，根据结构-TLD耦合系统运动方程进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，得到若干TLD响应对应的TLD阻尼比数据点，对这些数据点进行拟合得到TLD响应与TLD阻尼比关系曲线，如图2中曲线1所示：
式中，
为敏感性分析得到的TLD阻尼比拟合结果，σr1为相应的TLD响应，α1、β1为拟合常数；
步骤3、基于TLD当前设计，包括TLD当前外形尺寸、液深以及内置耗能装置设计参数，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析，对测得的TLD响应进行解耦分离再进行模态参数识别，得到液体晃荡频率ft′，并拟合TLD阻尼比与TLD响应关系曲线。
在本实施例中，基于TLD当前外形尺寸、液深以及内置耗能装置设计参数，构建TLD缩尺模型，针对上述构建的TLD缩尺模型，进行有色噪声激励下的TLD减振性能振动台试验，在得到试验结果后换算到TLD缩尺模型对应的原型，从而得到液体晃荡频率ft′。或者，在本实施例中，基于TLD当前外形尺寸、液深以及内置耗能装置设计参数，采用有限元软件建立TLD数值模型，根据实际试验条件在软件中设置相应初始条件、边界条件与外部激励，再通过求解流体运动方程获得流体晃动响应，从而实现数值模拟分析，得到液体晃荡频率ft′。
上述TLD缩尺模型为依据TLD当前设计结构为原型结构，按照一定相似准则进行设计制作；具体的TLD缩尺模型和TLD原型结构相似准则为：
长度缩尺比为：
质量缩尺比为：
时间缩尺比为：
频率缩尺比为：
加速度缩尺比为：
其中，n为常数，代表长度缩尺比值，Lm为TLD缩尺模型长度，Lp为TLD原型结构长度，Mm为TLD缩尺模型质量，Mp为TLD原型结构质量，ρ为液体密度，Vm为TLD缩尺模型体积，Vp为TLD原型结构体积，Tm为TLD缩尺模型周期，Tp为TLD原型结构周期，fm为TLD缩尺模型频率，fp为TLD原型结构频率，am为TLD缩尺模型加速度，ap为TLD原型结构加速度。
本实施例中，本步骤进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析后，对测得的TLD响应进行解耦分离再进行模态参数识别，得到若干TLD响应对应的TLD阻尼比数据点，对这些数据点进行拟合可得TLD响应与TLD阻尼比关系曲线，如图2中所示曲线2所示：
式中，
为减振性能分析得到的TLD阻尼比拟合结果，σr2为相应的TLD响应，α2、β2为拟合常数。
本实施例中，本步骤中所提到的模态参数识别方法可以是曲线拟合法、随机减量法、改进的贝叶斯方法等参数识别方法。其中，通过改进的贝叶斯方法对解耦分离得到的模特响应信号进行模态参数识别的过程如下：
S1、构造模态坐标下的激励功率谱密度函数
其中：
式中，S0为输入激励在固有频率处的功率谱密度值，β为荷载幂指数，fj为TLD的第j阶模态频率，f为频率变量；
S2、计算TLD模态响应信号y(t)功率谱密度Sy,N(f')的期望；其中，TLD模态响应y(t)为：
其中
为解耦后的第j阶模态响应，η(t)为测量噪声和模型误差等造成的实际响应和测量系统所获得响应间的预测误差；
计算TLD模态响应y(t)功率谱密度Sy,N(f')的期望：
式中，Sη为η(t)的功率谱密度分布；f'＝k△f，k＝0,1,…,int(N/2)，int表示取实数的整数部分，
T为采样时间；θ为需要识别的TLD模态参数，包括TLD模态频率和阻尼比参数值；
为TLD模态响应的功率谱密度；N为采样总数；
其中，
根据随机振动理论，得TLD模态响应的功率谱密度为：
式中，
为激励功率谱密度，
为TLD第j阶阻尼比；
S3、计算Sy,N(f')的概率密度函数；具体为：对于平稳过程，在特定的频率区间段，Sy，N(f')的概率密度函数p(Sy,N(f')|θ)近似为Chi-square分布：
当k1≠k2时，Sy,N(k1△f)与Sy,N(k2△f)的联合概率密度分布函数可表示为：
式中，
为[k1△f,k2△f]频率区间段内TLD模态响应的功率谱密度，其中k1，k2＝0,1,…,int(N/2))；
S4、计算TLD模态参数最优解，得到TLD的多阶频率f1、f2、f3...fN'，阻尼比ζ1、ζ2、ζ3...ζN'，N′为TLD所求频率和阻尼的总阶数；本步骤，通过求解函数
的最小值来确定TLD模态参数最优解，从而得到TLD的多阶频率f1、f2、f3...fN'，阻尼比ζ1、ζ2、ζ3...ζN'；其中，p(θ)为先验概率密度函数；上述计算得到的f1即为液体晃荡频率ft′。
另外，曲线拟合法与上述改进的贝叶斯方法采用相同的模态激励功率谱密度函数假设，根据模态响应信号的功率谱密度可以采用最小二乘法局部拟合出TLD频率和阻尼比。随机减量法将模态响应信号视为由两部分线性叠加而成，即初始条件决定的自由振动响应以及外荷载激励产生的随机强迫振动部分，采用相同的阈值按照相等长度截取信号，并对所得到的一系列等长信号进行平均，得到自由振动衰减曲线，对该曲线按照有阻尼系统的自由振动方程进行拟合，可得TLD频率、阻尼比。
步骤4、基于步骤2和步骤3拟合TLD响应与TLD阻尼比关系曲线，计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
在本实施例中，联立步骤2和步骤3拟合得到的TLD响应与TLD阻尼比关系曲线，具体如下：
两公式联立求交点，此时
σr1＝σr2；则
得到
变换公式进一步得到：
将该公式代入
则计算得到：
步骤5、判断步骤3中得到液体晃荡频率ft′和步骤4中计算得到的阻尼比值
能否使受控结构响应满足控制目标，若否，进入步骤6，若是，进入步骤7；在本实施例中，将步骤3中得到液体晃荡频率ft′与液体晃荡频率ft进行比较，若两者之间的差值超过阈值F1，则判定TLD当前设计使受控结构响应不满足控制目标；
步骤6、针对TLD外形尺寸、液深、内置耗能装置设计参数中的其中一个或多个进行修改，然后返回步骤3；其中TLD内置耗能装置包括装置的外形尺寸大小、装置数量以及装置摆放位置，TLD内置耗能装置包括TLD内部的格栅、可旋转隔板、桨柱等。
具体的：
当ft′与ft的差值超过阈值F1，且ft′小于ft，则减小TLD的长度或增大TLD的液深；
当ft′与ft的差值超过阈值F1，且ft′大于ft，则增大TLD的长度或减小TLD的液深；
当
与
的差值超过阈值F2，且
小于
则增大TLD内置耗能装置的尺寸、增加TLD内置耗能装置的数量或将TLD内置装置摆放的位置更加靠近TLD的中间位置处；
当
与
的差值超过阈值F2，且
大于
则减小TLD内置耗能装置的尺寸、减少TLD内置耗能装置的数量或将TLD内置装置摆放的位置更加远离TLD的中间位置处。
步骤7、将最终获取到的TLD外形尺寸、液深和内置耗能装置设计参数作为TLD的最终设计结果。
本实施例方法在TLD初步设计基础上，计算结构-TLD耦合系统在不同重现期风荷载作用下的响应，并进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析，得到更为准确的结构-TLD耦合系统设计参数，以此为依据进行TLD设计参数的调整，直到获取到的TLD设计参数ft′、
能否使受控结构响应满足控制目标。
本领域技术人员可以理解，实现本实施例方法中的全部或部分步骤可以通过程序来指令相关的硬件来完成，相应的程序可以存储于计算机可读存储介质中。应当注意，尽管在附图中以特定顺序描述了本实施例1的方法操作，但是这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。相反，描绘的步骤可以改变执行顺序，有些步骤也可以同时执行。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。
实施例2
本实施例公开了一种调谐液体阻尼器设计装置，包括如下功能模块：
TLD初步设计模块，用于获取设计要求的结构-TLD耦合系统的质量比μ、液体晃荡频率ft和TLD阻尼比
通过等效线性化方法对TLD进行初步设计，得到初步设计的TLD外形尺寸、液深以及内置耗能装置设计参数；
构建模块，用于结构-TLD耦合系统运动方程；
敏感分析和拟合模块，用于对结构-TLD耦合系统进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，并拟合TLD响应与TLD阻尼比关系曲线；
模态参数识别和拟合模块，在基于TLD当前设计，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析的情况下，用于对测得的TLD响应进行解耦分离再进行模态参数识别，得到液体晃荡频率ft′，并拟合TLD阻尼比与TLD响应关系曲线；
阻尼比计算模块，用于基于敏感分析和拟合模块与模态参数识别和拟合模块拟合得到TLD响应与TLD阻尼比关系曲线，计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
判定模块，用于判断模态参数识别和拟合模块得到液体晃荡频率ft′与阻尼比计算模块计算得到的阻尼比值
能否使受控结构响应满足控制目标；
修改模块，用于在判定模块判定出液体晃荡频率ft′和阻尼比值
不能使得受控结构响应满足控制目标时，则针对TLD外形尺寸、液深、内置耗能装置设计参数中的其中一个或多个进行修改；修改后由模态参数识别和拟合模块以及阻尼比计算模块重新计算液体晃荡频率ft′和阻尼比值
并通过判定模块再进行判定，直到通过判定模块判定出得到的液体晃荡频率ft′和阻尼比值
使得受控结构响应满足控制目标，将最终获取到的TLD外形尺寸、液深和内置耗能装置设计参数作为TLD的最终设计结果。
本实施例上述各个模块的具体实现可以参见上述实施例1，在此不再一一赘述。需要说明的是，本实施例提供的装置仅以上述各功能模块的划分进行举例说明，在实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。
实施例3
本实施例公开了一种存储介质，存储有程序，所述程序被处理器执行时，实现实施例1所述的调谐液体阻尼器设计方法，如下：
步骤1、获取设计要求的结构-TLD耦合系统的质量比μ、液体晃荡频率ft和TLD阻尼比
通过等效线性化方法对TLD进行初步设计，得到初步设计的TLD外形尺寸、液深以及内置耗能装置设计参数。
步骤2、构建结构-TLD耦合系统运动方程，对结构-TLD耦合系统进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，并拟合TLD响应与TLD阻尼比关系曲线。
步骤3、基于TLD当前设计，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析，对测得的TLD响应进行解耦分离再进行模态参数识别，得到液体晃荡频率ft′，并拟合TLD阻尼比与TLD响应关系曲线。
步骤4、基于步骤2和步骤3拟合TLD响应与TLD阻尼比关系曲线，计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
步骤5、判断步骤3中得到液体晃荡频率ft′和步骤4中计算得到的阻尼比值
能否使受控结构响应满足控制目标；
若否，则针对TLD外形尺寸、液深、内置耗能装置设计参数中的其中一个或多个进行修改，然后返回步骤3；
若是，则将最终获取到的TLD外形尺寸、液深和内置耗能装置设计参数作为TLD的最终设计结果。
在本实施例中，存储介质可以是磁盘、光盘、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、U盘、移动硬盘等介质。
实施例4
本实施例公开了一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，其特征在于，所述处理器执行存储器存储的程序时，实现实施例1所述的调谐液体阻尼器设计方法，如下：
步骤1、获取设计要求的结构-TLD耦合系统的质量比μ、液体晃荡频率ft和TLD阻尼比
通过等效线性化方法对TLD进行初步设计，得到初步设计的TLD外形尺寸、液深以及内置耗能装置设计参数。
步骤2、构建结构-TLD耦合系统运动方程，对结构-TLD耦合系统进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，并拟合TLD响应与TLD阻尼比关系曲线。
步骤3、基于TLD当前设计，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析，对测得的TLD响应进行解耦分离再进行模态参数识别，得到液体晃荡频率ft′，并拟合TLD阻尼比与TLD响应关系曲线。
步骤4、基于步骤2和步骤3拟合TLD响应与TLD阻尼比关系曲线，计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
步骤5、判断步骤3中得到液体晃荡频率ft′和步骤4中计算得到的阻尼比值
能否使受控结构响应满足控制目标；
若否，则针对TLD外形尺寸、液深、内置耗能装置设计参数中的其中一个或多个进行修改，然后返回步骤3；
若是，则将最终获取到的TLD外形尺寸、液深和内置耗能装置设计参数作为TLD的最终设计结果。
本实施例中，计算设备可以是台式电脑、笔记本电脑、PDA手持终端、平板电脑等终端设备。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。
Claims (10)
1.一种调谐液体阻尼器设计方法，其特征在于，包括步骤：
步骤1、获取设计要求的结构-TLD耦合系统的质量比μ、液体晃荡频率ft和TLD阻尼比
通过等效线性化方法对TLD进行初步设计，得到初步设计的TLD的初步设计参数；
步骤2、构建结构-TLD耦合系统运动方程，对结构-TLD耦合系统进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，并拟合TLD响应与TLD阻尼比关系曲线；
步骤3、基于TLD当前设计参数，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析，对测得的TLD响应进行解耦分离，再进行模态参数识别，得到液体晃荡频率ft′，并拟合TLD阻尼比与TLD响应关系曲线；
步骤4、基于步骤2和步骤3拟合TLD响应与TLD阻尼比关系曲线，计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
步骤5、判断步骤3中得到液体晃荡频率ft′和步骤4中计算得到的阻尼比值
能否使受控结构响应满足控制目标；
若否，则针对TLD设计参数进行修改，然后返回步骤3；
若是，则将最终获取到的TLD设计参数作为TLD的最终设计结果。
2.根据权利要求1所述的调谐液体阻尼器设计方法，其特征在于，步骤2中构建的结构-TLD耦合系统运动方程为：
式中：
Ms、Cs和Ks分别为结构质量、阻尼和刚度矩阵，由设计信息可得；
mt、kt和ct分别为TLD质量、刚度和阻尼，β为质量参与系数，由TLD初步设计结果得到；
假设TLD与结构第n个自由度耦合，则X、xn和ut分别为结构位移向量、结构第n个自由度位移、TLD位移；
H为TLD位置向量；F为风荷载向量；
3.根据权利要求1所述的调谐液体阻尼器设计方法，其特征在于，步骤3中，基于TLD当前设计参数，包括外形尺寸、液深以及内置耗能装置设计参数，构建TLD缩尺模型，针对上述构建的TLD缩尺模型，进行有色噪声激励下的TLD减振性能振动台试验，在得到试验结果后换算到TLD缩尺模型对应的原型，从而得到液体晃荡频率ft′；
或者步骤3中，基于TLD当前设计参数，包括外形尺寸、液深以及内置耗能装置设计参数，采用有限元软件建立TLD数值模型，根据实际试验条件在软件中设置相应初始条件、边界条件与外部激励，再通过求解流体运动方程获得流体晃动响应，从而实现数值模拟分析。
4.根据权利要求1所述的调谐液体阻尼器设计方法，其特征在于，步骤2中，根据结构-TLD耦合系统运动方程进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，得到若干TLD响应对应的TLD阻尼比数据点，对这些数据点进行拟合得到TLD响应与TLD阻尼比关系曲线：
式中，
为敏感性分析得到的TLD阻尼比拟合结果，σr1为相应的TLD响应，α1、β1为拟合常数；
步骤3中，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析后，对测得的TLD响应进行解耦分离再进行模态参数识别，得到若干TLD响应对应的TLD阻尼比数据点，对这些数据点进行拟合可得TLD响应与TLD阻尼比关系曲线：
式中，
为减振性能分析得到的TLD阻尼比拟合结果，σr2为相应的TLD响应，α2、β2为拟合常数；
步骤4中，基于步骤2中拟合得到的TLD响应与TLD阻尼比关系曲线
以及步骤3中拟合得到的TLD响应与TLD阻尼比关系曲线
计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
为：
5.根据权利要求1所述的调谐液体阻尼器设计方法，其特征在于，在步骤5中：将步骤5中得到液体晃荡频率ft′与液体晃荡频率ft进行比较，若两者之间的差值超过阈值F1，则判定TLD当前设计使受控结构响应不满足控制目标，对TLD的设计参数，包括TLD的外形尺寸长度或液深进行修改；
步骤5中，将步骤4中得到阻尼比值
与TLD阻尼比
进行比较，若两者之间的差值超过阈值F1，则判定TLD当前设计使受控结构响应不满足控制目标，对TLD的设计参数，包括TLD内置耗能装置设计参数进行修改。
6.根据权利要求5所述的调谐液体阻尼器设计方法，其特征在于，在步骤5中：
当ft′与ft的差值超过阈值F1，且ft′小于ft，则减小TLD的长度或增大TLD的液深；
当ft′与ft的差值超过阈值F1，且ft′大于ft，则增大TLD的长度或减小TLD的液深；
当
与
的差值超过阈值F2，且
小于
则增大TLD内置耗能装置的尺寸、增加TLD内置耗能装置的数量或将TLD内置装置摆放的位置更加靠近TLD的中间位置处；
当
与
的差值超过阈值F2，且
大于
则减小TLD内置耗能装置的尺寸、减少TLD内置耗能装置的数量或将TLD内置装置摆放的位置更加远离TLD的中间位置处。
7.根据权利要求1至6中任一项所述的调谐液体阻尼器设计方法，其特征在于，TLD内置耗能装置采用纤维增强塑料制作。
8.一种调谐液体阻尼器设计装置，其特征在于，包括：
TLD初步设计模块，用于获取设计要求的结构-TLD耦合系统的质量比μ、液体晃荡频率ft和TLD阻尼比
通过等效线性化方法对TLD进行初步设计，得到初步设计的TLD的初步设计参数；
构建模块，用于结构-TLD耦合系统运动方程；
敏感分析和拟合模块，用于对结构-TLD耦合系统进行不同重现期风荷载作用下的TLD阻尼比敏感性分析，并拟合TLD响应与TLD阻尼比关系曲线；
模态参数识别和拟合模块，在基于TLD当前设计，进行有色噪声激励下的TLD减振性能振动台试验或数值模拟分析的情况下，用于对测得的TLD响应进行解耦分离再进行模态参数识别，得到液体晃荡频率ft′，并拟合TLD阻尼比与TLD响应关系曲线；
阻尼比计算模块，用于基于敏感分析和拟合模块与模态参数识别和拟合模块拟合得到TLD响应与TLD阻尼比关系曲线，计算得到不同重现期风荷载作用下的TLD所能提供的阻尼比值
判定模块，用于判断模态参数识别和拟合模块得到液体晃荡频率ft′与阻尼比计算模块计算得到的阻尼比值
能否使受控结构响应满足控制目标；
修改模块，用于在判定模块判定出液体晃荡频率ft′和阻尼比值
不能使得受控结构响应满足控制目标时，则针对设计参数进行修改；修改后由模态参数识别和拟合模块以及阻尼比计算模块重新计算液体晃荡频率ft′和阻尼比值
并通过判定模块再进行判定，直到通过判定模块判定出得到的液体晃荡频率ft′和阻尼比值
使得受控结构响应满足控制目标，将最终获取到的TLD设计参数作为TLD的最终设计结果。
9.一种存储介质，存储有程序，其特征在于，所述程序被处理器执行时，实现权利要求1～6中任一项所述的调谐液体阻尼器设计方法。
10.一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，其特征在于，所述处理器执行存储器存储的程序时，实现权利要求1～6中任一项所述的调谐液体阻尼器设计方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110334937.2A CN113065208B (zh) | 2021-03-29 | 2021-03-29 | 调谐液体阻尼器设计方法、装置、介质和设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110334937.2A CN113065208B (zh) | 2021-03-29 | 2021-03-29 | 调谐液体阻尼器设计方法、装置、介质和设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113065208A true CN113065208A (zh) | 2021-07-02 |
| CN113065208B CN113065208B (zh) | 2021-11-23 |
Family
ID=76564282
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110334937.2A Active CN113065208B (zh) | 2021-03-29 | 2021-03-29 | 调谐液体阻尼器设计方法、装置、介质和设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113065208B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113673012A (zh) * | 2021-07-29 | 2021-11-19 | 华南理工大学 | 基于风洞试验的结构-tld系统的耦合振动实现方法 |
| CN113686528A (zh) * | 2021-07-28 | 2021-11-23 | 华南理工大学 | 一种结构-tld系统的子系统动力特性检测方法 |
| CN114528618A (zh) * | 2022-01-24 | 2022-05-24 | 华南理工大学 | 一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法 |
| CN114528652A (zh) * | 2022-01-13 | 2022-05-24 | 华南理工大学 | 一种内置桨柱调谐液体阻尼器的设计方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109139769A (zh) * | 2018-07-23 | 2019-01-04 | 大连理工大学 | 基于特高压复合套管的多环形调谐液体阻尼器及其方法 |
| CN109299571A (zh) * | 2018-10-25 | 2019-02-01 | 内蒙古科技大学 | 一种基于建筑消防水箱的调谐液体阻尼器的设计方法 |
| CN110823537A (zh) * | 2019-10-16 | 2020-02-21 | 华南理工大学 | 一种超高层建筑风振控制阻尼器性能的检测评价方法 |
| CN111751070A (zh) * | 2020-07-08 | 2020-10-09 | 华南理工大学 | 调谐液体阻尼器性能参数的检测方法、系统、设备和介质 |
| CN112000144A (zh) * | 2020-08-14 | 2020-11-27 | 武汉理工大学 | 基于涡流调谐质量阻尼器的应力带桥车致振动控制方法 |
-
2021
- 2021-03-29 CN CN202110334937.2A patent/CN113065208B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109139769A (zh) * | 2018-07-23 | 2019-01-04 | 大连理工大学 | 基于特高压复合套管的多环形调谐液体阻尼器及其方法 |
| CN109299571A (zh) * | 2018-10-25 | 2019-02-01 | 内蒙古科技大学 | 一种基于建筑消防水箱的调谐液体阻尼器的设计方法 |
| CN110823537A (zh) * | 2019-10-16 | 2020-02-21 | 华南理工大学 | 一种超高层建筑风振控制阻尼器性能的检测评价方法 |
| CN111751070A (zh) * | 2020-07-08 | 2020-10-09 | 华南理工大学 | 调谐液体阻尼器性能参数的检测方法、系统、设备和介质 |
| CN112000144A (zh) * | 2020-08-14 | 2020-11-27 | 武汉理工大学 | 基于涡流调谐质量阻尼器的应力带桥车致振动控制方法 |
Non-Patent Citations (4)
| Title |
|---|
| AL JAMAL等: "Numerical Investigation of Tuned Liquid Damper Performance Attached to a Single Degree of Freedom Structure", 《MECHANICAL ENGINEERING》 * |
| 宗刚等: "减震调谐液体阻尼器参数优化设计", 《同济大学学报(自然科学版)》 * |
| 嵇学培等: "广州西塔风致响应和气动阻尼特性的试验研究", 《土木工程学报》 * |
| 顾家扬等: "建筑结构阻尼比参数的优化识别方法", 《同济大学学报(自然科学版)》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113686528A (zh) * | 2021-07-28 | 2021-11-23 | 华南理工大学 | 一种结构-tld系统的子系统动力特性检测方法 |
| CN113673012A (zh) * | 2021-07-29 | 2021-11-19 | 华南理工大学 | 基于风洞试验的结构-tld系统的耦合振动实现方法 |
| CN113673012B (zh) * | 2021-07-29 | 2023-08-22 | 华南理工大学 | 基于风洞试验的结构-tld系统的耦合振动实现方法 |
| CN114528652A (zh) * | 2022-01-13 | 2022-05-24 | 华南理工大学 | 一种内置桨柱调谐液体阻尼器的设计方法 |
| CN114528618A (zh) * | 2022-01-24 | 2022-05-24 | 华南理工大学 | 一种连体塔楼和多个调谐液体阻尼器的风致振动实现方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113065208B (zh) | 2021-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113065208B (zh) | 调谐液体阻尼器设计方法、装置、介质和设备 | |
| Oberkampf et al. | Measures of agreement between computation and experiment: validation metrics | |
| CN111751070B (zh) | 调谐液体阻尼器性能参数的检测方法、系统、设备和介质 | |
| Karsh et al. | A hybrid stochastic sensitivity analysis for low-frequency vibration and low-velocity impact of functionally graded plates | |
| CN110110437B (zh) | 一种基于相关区间不确定性理论的汽车高频噪声预测方法 | |
| Clark Jr et al. | Engineering design exploration using locally optimized covariance kriging | |
| CN107676273B (zh) | 一种低振高效离心泵叶轮数值优化方法 | |
| CN108256693B (zh) | 一种光伏发电功率预测方法、装置及系统 | |
| CN111783351A (zh) | 针对结构系统不确定性参数的非概率可信集合定量化方法 | |
| Huang et al. | A Kriging-based non-probability interval optimization of loading path in T-shape tube hydroforming | |
| Khodaparast et al. | Fuzzy model updating and its application to the DLR AIRMOD test structure | |
| CN107704672B (zh) | 一种主动减振的船用离心泵浮筏系统多工况数值优化方法 | |
| Ling et al. | An efficient method for estimating time-dependent global reliability sensitivity | |
| Bedonian et al. | An Adaptive Sampling Method for Gradient-Enhanced Surrogate Modeling | |
| CN114528652A (zh) | 一种内置桨柱调谐液体阻尼器的设计方法 | |
| Droste et al. | Parametric model order reduction for vibroacoustic metamaterials based on modal superposition | |
| Chen et al. | Research on prediction and control method of car exterior noise | |
| Dong et al. | A Novel Adaptive Bayesian Model Averaging-Based Multiple Kriging Method for Structural Reliability Analysis | |
| Papadimitriou et al. | Information-Driven Modeling of Structures Using a Bayesian Framework | |
| Papadioti | Management of uncertainties in structural response and reliability simulations using measured data | |
| Koo et al. | Optimal dynamic vibration absorber design for minimizing the band-averaged input power using the residue theorem | |
| Marulo et al. | Probabilistic analysis and experimental results of sound transmission loss of composite panels | |
| Setúbal et al. | Force Identification from Vibration Data by Response Surface and Random Forest Regression Algorithms. Energies 2022, 15, 3786 | |
| CN118070431B (zh) | 一种飞行器轻质结构安全评估方法 | |
| Kim et al. | Active learning-based structural design optimization under constraints on first-passage probability |
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 |