建筑结构表面风压非高斯特性分析方法研究

发布时间：2020-11-18 07:17

　　 建筑结构抗风设计需要充分考虑风荷载的非高斯特征,不容忽视。本文首先回顾了非高斯风荷载研究领域的几个热点问题,综述了研究现状,指出了缺陷和不足。在此基础上,针对非高斯风压概率分布拟合、非高斯风压极值计算、非高斯风压过程模拟以及非高斯风压插值预测四个问题,展开深入研究。全文主要内容及结论如下:(1)现有基于穿越率理论的非高斯风压极值计算方法适用范围有限,且对于短尾侧极值的计算精度明显低于长尾侧。鉴于此,本文开发了“分离描述”(separate description,简称SD)算法以解决这些问题。SD算法采用Johnson变换作为潜在高斯过程与目标非高斯过程之间的转换工具,其适用范围覆盖了整个Pearson平面,适用于任意偏度、峰度组合的非高斯分布。由于短尾侧数据在极大似然函数的计算中所占权重更高,因而提出对母体风压分布拟合两次,以极大似然估计的拟合结果计算短尾侧极值,以矩估计拟合结果计算长尾侧极值,以针对性地改善短尾侧极值计算精度。经对比验证,SD算法的整体计算误差在4%以内,其精度高于传统算法,优势在短尾侧以及处理非高斯性较强的软过程时尤为明显。(2)现有基于极值理论的非高斯风压极值算法需要进行长时间风洞试验,其资源耗费较大。为解决这一问题,本文提出一种试验数据与数值模拟相结合(hybrid measurement and simulation-based,简称HMSB)的算法,通过数值模拟的手段,以前四阶统计量和功率谱密度(power spectral density,简称PSD)为口标,仿真生成大量风压时程,以获取足够多的极值样本,拟合极值分布。考虑到统计量对于极值计算精度影响更大,在保证统计量模拟精度的基础上,为了最大程度地提高模拟效率,提出了一种简化模拟方法,可以实现统计量的精确模拟和PSD的近似模拟。经验证,HMSB法计算精度高于传统方法,整体计算误差在4%左右。其计算效率极高,能在短时间内生成大量模拟时程,从而有效地节省了试验资源。(3)平稳非高斯过程模拟包含频率和概率两方面目标,传统方法大多基于先干涉频率再干涉概率的模拟思路。本文采用新的模拟思路,即先干涉概率再干涉频率,提出了一种新的模拟方法。在此过程中,推导了线性滤波系统输入过程与输出过程低阶统计量之间的转换关系,解决了统计量扭曲问题;探讨了新思路与传统思路所各自面临的不相容问题的成因,证明了两种不相容区间并不完全重叠,新方法可以解决某些传统方法的不相容问题;通过多个算例证明了新方法的有效性和准确性。此外,所提方法还具有模拟高阶相关非高斯过程的潜力。(4)现有非高斯非平稳过程模拟方法较少,且均需要迭代计算,模拟效率较低。本文基于线性滤波技术,提出了新的模拟方法,使用时变自回归(time-varying auto-regressive,简称TVAR)模型将潜在非高斯非平稳白噪声过滤成目标非高斯非平稳过程。新方法无需迭代,计算简便。在其开发过程中,推导了TVAR模型输入过程与输出过程时变低阶统计量之间的变换关系,从而可在模拟前根据目标时变低阶统计量计算出潜在非高斯非平稳白噪声的时变低阶统计量;将传统Johnson变换升级为时变版本,用于生成非平稳白噪声输入;提出了一种简便方法,用于确定TVAR模型阶数;通过算例证明了新方法的可行性和准确性。(5)将高斯过程回归(Gaussian process regression,简称GPR)技术引入到风压统计量的插值估计问题,相比于现有方法中使用的人工神经网络技术,GPR的优势在于超参数自适应选取、输出具有明确的概率意义等方面。对于低阶统计量估计问题,本文所提方法相对于传统方法精度更高,能够实现对未布测点位置和未测量风向角的风压信息的准确补充。对于风压时程估计问题,本文充分利用GPR善于处理小样本问题的特点,提出了时变GPR估计方法,其精度相对传统方法更高,且能够处理风场随时间变化的非平稳情形。尝试估计了高阶统计量和累积密度函数,其中三阶统计量和累积密度函数估计精度较高,可以满足使用要求。而四阶统计量由于其随机性更强、波动幅度更高,精度尚无法满足要求,需进一步提升算法。
【学位单位】：大连理工大学
【学位级别】：博士
【学位年份】：2019
【中图分类】：TU312.1
【部分图文】：

的所有分布类型。第三阶和第四阶统计量，即偏度（ｓｋｅｗｎｅｓｓ，简称从）和峰度（ｋｕｒｔｏｓｉｓ，??简称ｈ?）是用来区分高斯分布与非高斯分布的重要参数，也是用来描述非高斯分布概率??性质的重要指标。图１．１展示了偏度、峰度与随机变量概率密度函数之间的联系。由图??１．１（ａ）可见，偏度描述分布的偏斜程度，正偏度向左，负偏度向右，绝对值越大表示越??向一侧倾斜，偏度为０则为对称分布。峰度描述分布的尖削程度，峰度越大表示分布中??心越高耸，尾部越长。由于高斯分布的偏度为０，峰度为３，因此对于任意偏度不等于０??或峰度不等于３的分布都可归类为非高斯分布，而其偏度与峰度的大小可以用来衡量其??偏离高斯分布的程度。??绝大多数非高斯风压分布是非对称的，相较于高斯分布其一侧尾部较长，一侧尾部??较短。体现在风压时程上则具有明显的非对称特征，一侧伴有大幅度的风压脉冲。这些??大幅脉冲更容易使结构发生破坏，是导致局部结构构件失效的主要原因。许多知名建筑??曾在风荷载作用下发生过不同程度的毁坏

，。分布也有类似的规律，在此不一一列举，可查阅文献Ｉ４７］。这些传，当待分析数据与某一分布适用范围更接近时，该分布的拟合研究者们在寻找最适合描述非高斯风压的分布时，由于他们的率特性，因而得出了不同的结论。??分析非高斯风压数据时传统指数类分布已较少使用，然而该类于目前常用的Ｈｅｒｍｉｔｅ转换分布和非参数分布，传统指数类分便。只是应当摒弃使用一类分布描述所有数据的习惯，应熟悉据数据类型选择合适的分布。??ｉｔｅ转换分布??理论人为任意非高斯过程都可以由一个与之对应的标准高斯调非线性转换函数变换得到，即??ｙ?＝?ｇ｛ｕ＼?ｕ?＝?ｇ￣＇（ｙ）高斯过程为均值等于０、方差等于１的标准高斯过程；

建筑结构表面风压非高斯特性分析方法研究?－３和％?．Ｃ７２?．Ｃ７３＝４￣Ｋｓ／：２的关系，进而可以推出三个根必是一正１＞０，所以Ｃ７的值应该与正根的值相等。将极带入式（１．１２）可Ａｔ／，如果ｈ?＜知＇则选择＆系统，即式（〗．９ｂ）；如果ｔｏ?＞?Ａ＾１．９ｃ）。??２５?■?Ｉ?？?？?？?？?Ｉ?＿?Ｉ?，?Ｉ?，?Ｉ?■?Ｉ?■?Ｉ?１?／?Ｉ?■．??、＼?、？?／?／??
【参考文献】

相关期刊论文 前10条

1 李寿科;李寿英;陈政清;孙洪鑫;;基于非高斯仿真的风压系数极值计算方法[J];振动与冲击;2014年24期

2 李锦华;吴春鹏;陈水生;;矩形结构非高斯风荷载特性研究[J];振动.测试与诊断;2014年05期

3 罗俊杰;宋雄彬;李璟;;考虑随机风压场非高斯特性的大跨度膜屋盖风振响应[J];土木建筑与环境工程;2013年06期

4 林巍;黄铭枫;楼文娟;;大跨屋盖脉动风压的非高斯峰值因子计算方法[J];建筑结构;2013年15期

5 王飞;全涌;顾明;;基于广义极值理论的非高斯风压极值计算方法[J];工程力学;2013年02期

6 陈伏彬;李秋胜;卢春玲;黄生洪;傅学怡;郭明;;复杂大跨结构屋盖风荷载特性的试验与计算研究[J];空气动力学学报;2012年05期

7 韩宁;顾明;;方形高层建筑风压脉动非高斯特性分析[J];同济大学学报(自然科学版);2012年07期

8 李秋胜;胡尚瑜;李正农;;低矮房屋风荷载实测研究(Ⅱ)——双坡屋面风压特征分析[J];土木工程学报;2012年04期

9 胡尚瑜;李秋胜;;低矮房屋风荷载实测研究(Ⅰ)——登陆台风近地边界层风特性[J];土木工程学报;2012年02期

10 李寿科;李寿英;陈政清;黄韬;;大跨开合式屋盖峰值风压的试验研究[J];振动与冲击;2010年11期

相关博士学位论文 前1条

1 孙瑛;大跨屋盖结构风荷载特性研究[D];哈尔滨工业大学;2007年

本文编号：2888465