CN1851436A - 大跨度屋盖和超高层建筑结构风振响应的检测计算方法 - Google Patents
大跨度屋盖和超高层建筑结构风振响应的检测计算方法 Download PDFInfo
- Publication number
- CN1851436A CN1851436A CN 200610035732 CN200610035732A CN1851436A CN 1851436 A CN1851436 A CN 1851436A CN 200610035732 CN200610035732 CN 200610035732 CN 200610035732 A CN200610035732 A CN 200610035732A CN 1851436 A CN1851436 A CN 1851436A
- Authority
- CN
- China
- Prior art keywords
- response
- displacement
- truncated
- modal
- function
- 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.)
- Pending
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种用于大跨度屋盖和超高层建筑结构风振响应的检测计算方法,本发明利用风洞多点同步测压技术测得大跨度屋盖和超高层建筑结构的时变气动荷载分布再采用随机振动理论计算得到结构的风振响应以及等效静力风荷载的技术领域;其特点在于它直接采用模态分析方法在频域计算结构在随机气动荷载结尾函数作用下的模态坐标和物理坐标的响应,其计算结果具有和传统精确的完全二次型组合(CQC)方法同样的精度,但计算速度要大大地快于CQC方法,且由于不需要直接计算激振力的功率谱密度矩阵而降低了对内存的需求,特别适合在中档PC机上实现,为大跨度屋盖和超高层结构的风振响应和等效静力风荷载的检测计算解决了两个关键技术问题。
Description
技术领域:
本发明涉及大跨度屋盖和超高层建筑结构的结构抗风分析领域,具体涉及建筑模型的风洞同步风压测量并计算此类结构的风致振动响应和结构的等效静力风荷载。
背景技术:
大跨度屋盖和超高层建筑结构受到平稳的随机风荷载的作用,在理论上是可以计算其任意位置上的响应的,在本发明之前,所采用的精确方法便是完全二次型组合(CQC)方法,结构风振响应分析的精确计算是进一步计算可用于结构设计的等效静力风荷载的关键。然而,由于大型的工程结构均有非常大的自由度,直接采用CQC方法存在计算量大和占用内存空间巨大的问题而影响了它的工程适应性,这些问题在计算大跨度屋盖结构在多点随机激励的风振响应中尤为明显,它最终导致的结果是计算效率低下,数据分析计算时间长,且一旦在结构参数调整后再计算又要花费同样的计算时间,往往难以满足结构设计的需要。因此,在保证计算精度不变的条件下,寻找一种既节省内存又高效快速的方法对于此类结构的工程抗风的快速和优化设计无疑具有非常重大的意义。
发明内容:
本发明所要解决的技术问题是大跨度屋盖和超高层建筑风振响应和等效静力风荷载的检测计算中的两个关键问题:CQC方法计算量太大和占有存储空间太大的问题。
本发明是通过下面的技术方案解决上述技术问题,其特征在于它包括以下步骤:
1)利用高速同步测压设备在风洞中对建筑模型表面的风压进行高速数据采样,并进行模/数转换(A/D转换),得到Np×1维的数字化离散风压时程数据样本矢量{PT(t)},其中t=iΔt,Δt为采样间隔,i=1,2,...N,j=1,2,...Np;N为样本点数,T=NΔt为采样时间,Np为同步测压测点数,采样的结果被存储在计算机里;本质上{PT(t)}相当于实际风压时程{p(t)}的截尾样本序列,即:
2)计算该截尾函数向量的傅立叶变换:
3)用模态分析法可得结构在以上截尾函数向量作用下的模态坐标和结点位移响应的傅立叶变换:
{Q(f,T)}=[hp(f)]*[Φ]T[R]{P(f,T)}
{Y(f,T)}=[Φ]]{Q(f,T)}
其中:[hp(ω)]为模态坐标下的频率响应函数矩阵,[Φ]为振型矩阵,[R]为荷载指示矩阵或扩展荷载分配矩阵,Q、Y分别表示模态坐标q(t)和结点位移y(t)的截尾函数的傅立叶变换;由以上结果可以分别计算模态响应和结点位移响应的功率谱密度函数的估计:
4)采用交叠取样方式可以得到不同时间段的截尾样本函数,重复过程2)3)对所计算出的[SqTqT(f)]和[SyTyT(f)]进行算术平均以提高谱估计的精度,最终得到模态响应和结点位移响应的功率谱密度函数:
相应模态响应的协方差矩阵为:
而结点位移响应的协方差矩阵为:
以上两个矩阵是进一步计算等效静力风荷载的关键。
本发明具有以下技术效果:
1、避免了直接计算随机荷载激励的功率谱密度矩阵Spp(f)可以大大地节省计算机内存的占用。
在采用以风洞试验的多点同步测压数据计算结构的风振响应中,激励的功率谱密度矩阵通常转换为测点风压的功率谱密度函数矩阵。由于需要计算不同频率值上的功率谱密度计算,[Spp(f)]实际上是个三维的复数矩阵,其长度取决于试验的测点数以及数据采用的样本长度,一般大跨度屋盖结构的同步测压通常需要400个测点以上,假定最终做功率谱密度计算的谱线数为1024(取决于数据采样的样本长度和谱分析中的分段数),这样以上测点风压所需要的功率谱密度矩阵所需要的内存量就为:16×400×400×1024/10243=2.44GB字节,这在一般的高档微机或工作站上也是难以实现的,而存放同样长度的风压时程函数傅立叶变换所需要的内存仅为6.25MB字节。
大量内存的占用使得即使是高性能的工作站都无法采用CQC方法计算,而本发明方法则使得复杂结构的风振分析问题可以在中档微机上实行。
2、事实上,正是由于避免直接对三维功率谱密度矩阵Spp(f)的运算以及此过程中蕴涵的运算简化使得后续模态响应和结点位移响应的功率谱密度计算所需的计算工作量也远比传统的CQC方法少。
以上解决了大型复杂结构尤其是大跨度屋盖风振响应的检测计算和等效静力风荷载计算中的两个关键问题,在一般的中档微机上即可实现此类结构的快速计算,可以满足结构设计分析的需要。
附图说明:
图1本发明的算法流程图;
图2本算法和CQC算法计算出某超高层建筑顶部动态位移随不同风向变化的计算结果比较;
图3某大跨度屋盖风振响应计算所需时间的对比;
图4某大跨度屋盖风振响应计算中数组所需要内存的对比。
具体实施方式:
应用多点同步测压设备测量结构表面的风压信号得到作用于建筑物的激励时程,直接由激励截尾函数的傅立叶变换计算模态坐标以及结点位移响应的傅立叶变换,然后计算模态响应和主要结点位移响应的功率谱密度函数,并由此得到用于等效静力风荷载计算模态坐标响应和结点位移响应的协方差。本发明算法流程见图1。
实际应用结果:
图2为采用本方法和CQC计算不同风向下某超高层建筑的最大动态位移的比较,可见本方法和CQC方法的计算结果是一致的。
图3和图4分别是计算中CPU时间和数组内存占用情况的比较,其计算对象是一个包括屋盖在内的大型体育中心结构,从图3、图4可看出,本发明占用CPU计算时间不到CQC方法的十分之一,内存占用只有CQC方法的五分之一左右,从中可以看到本发明的优越之处,应当指出的是,对于该算例,由于内存问题普通中档微机已不能满足其计算要求,该算例的计算是在DELL 670工作站上完成的,时间对比只具有相对意义;同时为了对比两种算法,对解题规模做了适当的限制,如果不做限制,有可能在很好的工作站也难以采用CQC方法进行计算,从而无法进行对比。
对于较为复杂的工程结构,节省内存空间的意义在于它解决了即使是在很好的工作站上CQC方法都不能计算的问题,本发明使复杂的大跨度屋盖结构的风致响应和等效静力风荷载的计算的实用化和工程应用成为可能。在结构设计的深化过程中,结构的参数会一直处在调整之中,一旦结构参数做了调整,风振分析和等效荷载计算就必须重新进行,本发明可以对结构修正作出最快速度的响应。
Claims (1)
1、一种用于大跨度屋盖和超高层建筑结构风振响应的检测计算方法,是通过上述建筑模型的风洞试验测试结构表面的风压分布,并据此进行计算实现的,其特征在于:
1)风洞中利用高速同步测压设备对建筑模型表面的风压进行高速的模/数转换,即A/D转换,得到Np×1维的数字化离散风压时程数据样本矢量{pT(t)},其中t=iΔt,Δt为采样间隔,i=1,2,...N,j=1,2,...Np;N为样本点数,T=NΔt为采样时间,Np为同步测压测点数,采样的结果被存储在计算机里;本质上{pT(t)}相当于实际风压时程{p(t)}的截尾样本序列,即:
2)计算该截尾序列向量的傅立叶变换:
3)用模态分析法得到结构在以上截尾函数向量作用下的模态坐标和结点位移响应的傅立叶变换:
{Q(f,T)}=[hp(f)]*[Φ]T[R]{P(f,T)}
{Y(f,T)}=[Φ]]{Q(f,T)}
其中:[hp(ω)]为模态坐标下的频率响应函数矩阵,[Ф)为振型矩阵,[R]为荷载分配矩阵,Q、Y分别表示模态坐标q(t)和结点位移y(t)的截尾函数的傅立叶变换;由以上结果分别计算模态响应和结点位移响应的功率谱密度函数的估计:
相应模态响应的协方差矩阵为:
而结点位移响应的协方差矩阵为:
[Cyy]=[Ф][Cqq][Ф]T
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200610035732 CN1851436A (zh) | 2006-05-31 | 2006-05-31 | 大跨度屋盖和超高层建筑结构风振响应的检测计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200610035732 CN1851436A (zh) | 2006-05-31 | 2006-05-31 | 大跨度屋盖和超高层建筑结构风振响应的检测计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1851436A true CN1851436A (zh) | 2006-10-25 |
Family
ID=37132923
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200610035732 Pending CN1851436A (zh) | 2006-05-31 | 2006-05-31 | 大跨度屋盖和超高层建筑结构风振响应的检测计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1851436A (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102426628A (zh) * | 2011-10-21 | 2012-04-25 | 上海大学 | 非同轴刚度偏心结构平扭耦合风振响应的分析方法 |
CN102589838A (zh) * | 2012-02-15 | 2012-07-18 | 浙江大学 | 单层建筑风洞试验模型内部气承刚度模拟装置 |
CN102663246A (zh) * | 2012-04-06 | 2012-09-12 | 上海大学 | 考虑ssi一般平面不规则建筑地震反应的分析方法 |
CN104535254A (zh) * | 2014-12-23 | 2015-04-22 | 太原科技大学 | 一种建筑物外表面风压测量方法 |
CN104866688A (zh) * | 2015-06-09 | 2015-08-26 | 同济大学建筑设计研究院(集团)有限公司 | 一种基于加速度敏感性分析的结构构件尺寸修正方法 |
CN104964806A (zh) * | 2015-06-25 | 2015-10-07 | 哈尔滨工业大学 | 一种风洞试验中高层建筑模型的吸吹气控制装置 |
CN105466661A (zh) * | 2016-01-08 | 2016-04-06 | 武汉理工大学 | 基于改进卡尔曼滤波的超高层建筑风荷载反分析方法 |
CN105479491A (zh) * | 2016-01-20 | 2016-04-13 | 蔡权 | 一种精确智能的机械手 |
CN106096102A (zh) * | 2016-06-02 | 2016-11-09 | 东南大学 | 基于fem‑kl的非平稳随机动态响应分析方法 |
CN106679798A (zh) * | 2016-12-15 | 2017-05-17 | 中国机械工业集团有限公司 | 一种高耸型建筑结构异层振动多点同步测试方法 |
CN107194036A (zh) * | 2017-04-24 | 2017-09-22 | 浙江大学 | 一种大跨屋盖结构非高斯风荷载数值模拟方法 |
CN109060292A (zh) * | 2018-06-01 | 2018-12-21 | 华南理工大学 | 考虑高频底座力天平试验的双耦合***的风振实现方法 |
CN115455685A (zh) * | 2022-09-05 | 2022-12-09 | 大连理工大学 | 一种二维紊流风激励下长索结构抖振响应的计算方法 |
CN116186859A (zh) * | 2023-03-01 | 2023-05-30 | 广州广检建设工程检测中心有限公司 | 一种大跨索屋盖面组分区的风振动力计算方法及其*** |
CN116467878A (zh) * | 2023-04-21 | 2023-07-21 | 广州非攻智造建筑科技有限公司 | 一种装配式建筑安全健康检测方法、***、设备及介质 |
-
2006
- 2006-05-31 CN CN 200610035732 patent/CN1851436A/zh active Pending
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102426628A (zh) * | 2011-10-21 | 2012-04-25 | 上海大学 | 非同轴刚度偏心结构平扭耦合风振响应的分析方法 |
CN102589838A (zh) * | 2012-02-15 | 2012-07-18 | 浙江大学 | 单层建筑风洞试验模型内部气承刚度模拟装置 |
CN102589838B (zh) * | 2012-02-15 | 2014-04-16 | 浙江大学 | 单层建筑风洞试验模型内部气承刚度模拟装置 |
CN102663246A (zh) * | 2012-04-06 | 2012-09-12 | 上海大学 | 考虑ssi一般平面不规则建筑地震反应的分析方法 |
CN102663246B (zh) * | 2012-04-06 | 2016-03-30 | 上海大学 | 考虑ssi一般平面不规则建筑地震反应的分析方法 |
CN104535254A (zh) * | 2014-12-23 | 2015-04-22 | 太原科技大学 | 一种建筑物外表面风压测量方法 |
CN104535254B (zh) * | 2014-12-23 | 2017-03-15 | 太原科技大学 | 一种建筑物外表面风压测量方法 |
CN104866688B (zh) * | 2015-06-09 | 2017-11-24 | 同济大学建筑设计研究院(集团)有限公司 | 一种基于加速度敏感性分析的结构构件尺寸修正方法 |
CN104866688A (zh) * | 2015-06-09 | 2015-08-26 | 同济大学建筑设计研究院(集团)有限公司 | 一种基于加速度敏感性分析的结构构件尺寸修正方法 |
CN104964806A (zh) * | 2015-06-25 | 2015-10-07 | 哈尔滨工业大学 | 一种风洞试验中高层建筑模型的吸吹气控制装置 |
CN105466661A (zh) * | 2016-01-08 | 2016-04-06 | 武汉理工大学 | 基于改进卡尔曼滤波的超高层建筑风荷载反分析方法 |
CN105479491A (zh) * | 2016-01-20 | 2016-04-13 | 蔡权 | 一种精确智能的机械手 |
CN106096102A (zh) * | 2016-06-02 | 2016-11-09 | 东南大学 | 基于fem‑kl的非平稳随机动态响应分析方法 |
CN106679798A (zh) * | 2016-12-15 | 2017-05-17 | 中国机械工业集团有限公司 | 一种高耸型建筑结构异层振动多点同步测试方法 |
CN107194036A (zh) * | 2017-04-24 | 2017-09-22 | 浙江大学 | 一种大跨屋盖结构非高斯风荷载数值模拟方法 |
CN109060292A (zh) * | 2018-06-01 | 2018-12-21 | 华南理工大学 | 考虑高频底座力天平试验的双耦合***的风振实现方法 |
CN109060292B (zh) * | 2018-06-01 | 2019-10-18 | 华南理工大学 | 考虑高频底座力天平试验的双耦合***的风振实现方法 |
CN115455685A (zh) * | 2022-09-05 | 2022-12-09 | 大连理工大学 | 一种二维紊流风激励下长索结构抖振响应的计算方法 |
CN115455685B (zh) * | 2022-09-05 | 2023-12-19 | 大连理工大学 | 一种二维紊流风激励下长索结构抖振响应的计算方法 |
CN116186859A (zh) * | 2023-03-01 | 2023-05-30 | 广州广检建设工程检测中心有限公司 | 一种大跨索屋盖面组分区的风振动力计算方法及其*** |
CN116186859B (zh) * | 2023-03-01 | 2023-09-08 | 广州广检建设工程检测中心有限公司 | 一种大跨索屋盖面组分区的风振动力计算方法及其*** |
CN116467878A (zh) * | 2023-04-21 | 2023-07-21 | 广州非攻智造建筑科技有限公司 | 一种装配式建筑安全健康检测方法、***、设备及介质 |
CN116467878B (zh) * | 2023-04-21 | 2024-01-26 | 广州非攻智造建筑科技有限公司 | 一种装配式建筑安全健康检测方法、***、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1851436A (zh) | 大跨度屋盖和超高层建筑结构风振响应的检测计算方法 | |
CN108345741B (zh) | 基于无网格rkpm各向异性材料二维热变形和热应力分析方法 | |
Bediz | A spectral-Tchebychev solution technique for determining vibrational behavior of thick plates having arbitrary geometry | |
Noor et al. | Recent advances in reduction methods for instability analysis of structures | |
CN100468422C (zh) | 利用交叉模型交叉模态的结构模型修正方法 | |
EP2866019A1 (en) | Air-water-soil skeleton coupled calculation device, coupled calculation method, and coupled calculation program | |
CN106096105A (zh) | 输电线路风偏瞬态响应计算方法 | |
CN102254068A (zh) | 大跨度桥梁抖振响应的多尺度分析方法 | |
CN113607373B (zh) | 基于离散多点测量数据的大跨度屋面风压分布快速重构算法 | |
CN105954802A (zh) | 一种岩性数据体的转换方法及装置 | |
Song et al. | Investigation on wind-induced aero-elastic effects of tall buildings by wind tunnel test using a bi-axial forced vibration device | |
CN104281730A (zh) | 一种大转动变形的板壳结构动响应的有限元分析方法 | |
Jian et al. | Buffeting response of cable-stayed bridge during construction under skew winds and pylon interference | |
CN103018729A (zh) | 金属圆柱定标体雷达散射截面的计算方法 | |
CN115358026A (zh) | 一种基于多元线性回归与曲面拟合的五孔探针数据处理方法 | |
CN115993649A (zh) | 基于等效方位杨氏模量的裂缝参数预测方法及*** | |
CN102567627A (zh) | 基于卫星重力梯度观测数据的圆环面调和分析方法 | |
CN111400898A (zh) | 基于主模态法和应变的阵列天线振动变形预测方法及设备 | |
Kardani et al. | On the application of the maximum entropy meshfree method for elastoplastic geotechnical analysis | |
CN111651920B (zh) | 塔线体系下大跨越输电塔风致机械能的传递分析方法 | |
Huang et al. | Evolutionary spectra estimation of a nonstationary process based on a single sample | |
Choi et al. | Efficient calculation of statistical moments for structural health monitoring | |
da Costa et al. | Bridge buffeting by skew winds: A quasi-steady case study | |
CN114818550B (zh) | 一种飞机振动试验中时变气动载荷地面等效模拟方法 | |
CN114492012A (zh) | 一种考虑带宽影响的多轴随机振动寿命预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |