CN107194036A - 一种大跨屋盖结构非高斯风荷载数值模拟方法 - Google Patents
一种大跨屋盖结构非高斯风荷载数值模拟方法 Download PDFInfo
- Publication number
- CN107194036A CN107194036A CN201710270627.2A CN201710270627A CN107194036A CN 107194036 A CN107194036 A CN 107194036A CN 201710270627 A CN201710270627 A CN 201710270627A CN 107194036 A CN107194036 A CN 107194036A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- simulation
- msub
- wind
- 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.)
- Withdrawn
Links
Classifications
-
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种大跨屋盖结构非高斯风荷载数值模拟方法。本发明所提供的方法改进了传统的单点非高斯风压模拟算法,将传统的单点非高斯风压模拟与本征正交分解(POD)技术、Hermite矩阵相结合,通过POD技术将风洞试验所得的风压场数据分解为主坐标和本征向量的组合,选取能量贡献率较大的前N阶POD本征向量对应的主坐标时程进行模拟,并与原本征向量组合形成模拟风压场。一方面,本发明通过POD技术不仅保持了风压场在时间与空间上的相关性,而且大大减小了模拟矩阵阶数,从而减小计算量,提高计算效率;另一方面,本发明采用Hermite矩阵将非高斯过程转化为高斯过程,并引入三个模拟误差评价参数修正模拟结果,更加符合实际大跨屋盖脉动风压场的特性。
Description
技术领域
本发明属于大跨屋盖结构抗风设计领域,涉及一种随机脉动风压场的模拟方法,特别涉及一种大跨屋盖结构非高斯风荷载数值模拟方法。
背景技术
风荷载是控制结构安全的重要荷载,特别是对于柔度较大的大跨度结构和高耸结构,有时甚至起决定性的作用。在考虑分离流作用的局部重要区域,风荷载常表现出强烈的脉动和非高斯特性。与高斯过程相比,在相同均值和方差的情况下,非高斯随机过程的峰值更大,会对结构产生更不利的效应,因而对非高斯风荷载的研究具有重要意义。对于体育场馆和干煤棚等大跨建筑,结构的几何非线性程度较高,相比于频域法,时域法可以考虑结构的非线性因素,能较精确地反映结构的风振耦合情况。然而,由于各种因素的影响,大跨屋盖结构风洞测压试验给出的风荷载时程数据都是有限的,难以满足研究非高斯风荷载和结构动力风效应的要求。传统的谐波叠加法虽然可以用于风压模拟,但随着测点数量的增加,计算量会大大增加,且无法模拟非高斯过程。因此,利用数值方法模拟获得满足特定条件的长时距非高斯风荷载是非常必要的。
发明内容
本发明的目的在于提供一种新的非高斯风压场模拟方法,尤其是提供一种大跨屋盖结构非高斯风荷载数值模拟方法,该方法采用本征正交分解(POD)技术将风洞试验所获取的目标零均值风压场分解为仅依赖时间的主坐标和仅随空间变化的本征向量的组合,根据各阶本征向量本征值大小,选取前N阶POD本征向量所对应的主坐标时程进行模拟,将前N阶POD主坐标时程用Hermite矩阵转化为高斯过程,用谐波叠加法(WAWS)模拟并修正,再与原本征向量组合形成模拟零均值风压场。计算得到的测点模拟脉动风压时程保持了风压场在时间与空间上的相关性,更加符合实际大跨屋盖脉动风压场的特性,可推广到其他屋盖结构抗风设计。
为此,本发明的上述目的通过以下技术方案来实现:
本发明提供一种大跨屋盖结构非高斯风荷载数值模拟方法,该方法包括以下步骤:
第一步:通过风洞测压试验,获取风压场:
P(x,y,t)={p1(t)、p2(t)...pn(t)},计算各测点风压均值:得到目标零均值风压场P*(x,y,t)={p1 *(t);p2 *(t)...pn *(t)};
第二步:进行本征正交分解(POD),目标零均值风压场P*(x,y,t)可以展开为:
式中:为第i阶POD主坐标时程,为第i阶本征向量,相应的本征值则记为
第三步:根据各阶本征向量的本征值大小,选取模拟的前N阶POD主坐标时程 及对应的本征向量
第四步:对选取后的前N阶POD主坐标时程进行模拟,模拟算法如下:
(1)计算第i阶POD主坐标时程的目标非高斯过程的功率谱、偏度和峰度,分别记为:
(2)根据Hermite矩变换关系,推导出第i阶POD主坐标时程自相关函数对应的高斯过程自相关函数
(3)对作快速傅里叶变换得到初始模拟高斯过程的功率谱采用谐波叠加法得到第j次模拟的高斯过程计算其功率谱,记为
(4)采用Hermite矩变换方法,得到第j次模拟非高斯过程计算其功率谱、偏度和峰度,分别记为:
(5)引入第j次模拟误差评价参数
为目标非高斯过程功率谱密度函数与第j次模拟得到的非高斯过程功率谱密度函数之间的相对误差,fn为离散的频率点,为目标高阶统计量与模拟得到的非高斯过程高阶统计量之间的误差;
(6)若或或其中,目标误差和则对(3)中模拟高斯过程功率谱进行修正,方法如下:
式中:β是修正常数,这里取值1.3,修正后返回至(2),进行下一步迭代计算;
若且且其中,目标误差和或者迭代步数大于50000步,停止迭代,得到的模拟非高斯过程即为第i阶POD主坐标时程模拟结果,记为并返回至(1),模拟第i+1阶POD主坐标时程;
(7)重复上述(1)至(6),直至完成前N阶POD主坐标时程模拟,模拟主坐标时程分别为:
第五步:将模拟得到的主坐标时程与原本征向量组合形成模拟零均值风压场:
本发明所提供的方法将谐波叠加法(WAWS)与本征正交分解(POD)技术、Hermite矩阵相结合,通过POD技术将风洞试验所得的风压场数据分解为主坐标和本征向量的组合,选取能量贡献率较大的前N阶POD本征向量对应的主坐标时程进行模拟,并与原本征向量组合形成模拟风压场。一方面,本发明通过POD技术不仅保持了风压场在时间与空间上的相关性,而且大大减小了模拟矩阵阶数,从而减小计算量,提高计算效率;另一方面,本发明采用Hermite矩阵将非高斯过程转化为高斯过程,并引入三个模拟误差评价参数修正模拟结果,更加符合实际情况。本发明所提供的方法可以推广到其他屋盖结构设计。
附图说明
图1为本发明的流程图。
图2为风洞试验模型测点布置图。
图3为前五阶本征向量等值线云图。
图4为前五阶模拟主坐标时程与目标时程功率谱密度对比图。
图5为代表性测点脉动风压模拟值与目标值的功率谱密度对比图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
本发明提供一种大跨屋盖结构非高斯风荷载数值模拟方法,其流程如图1所示,该方法包括以下步骤:
第一步:通过风洞测压试验,获取风压场:
P(x,y,t)={p1(t)、p2(t)...pn(t)},计算各测点风压均值:得到目标零均值风压场P*(x,y,t)={p1 *(t);p2 *(t)...pn *(t)};
第二步:进行本征正交分解,目标零均值风压场P*(x,y,t)可以展开为:
式中:为第i阶POD主坐标时程,为第i阶本征向量,相应的本征值则记为
第三步:根据各阶本征向量的本征值大小,选取模拟的前N阶POD主坐标时程 及对应的本征向量
第四步:对选取后的前N阶POD主坐标时程进行模拟,模拟算法如下:
(1)计算第i阶POD主坐标时程的目标非高斯过程的功率谱、偏度和峰度,分别记为:
(2)根据Hermite矩变换关系,推导出第i阶POD主坐标时程自相关函数对应的高斯过程自相关函数
(3)对作快速傅里叶变换得到初始模拟高斯过程的功率谱采用谐波叠加法得到第j次模拟的高斯过程计算其功率谱,记为
(4)采用Hermite矩变换方法,得到第j次模拟非高斯过程计算其功率谱、偏度和峰度,分别记为:
(5)引入第j次模拟误差评价参数
为目标非高斯过程功率谱密度函数与第j次模拟得到的非高斯过程功率谱密度函数之间的相对误差,fn为离散的频率点,为目标高阶统计量与模拟得到的非高斯过程高阶统计量之间的误差;
(6)若或或其中,目标误差和则对(3)中模拟高斯过程功率谱进行修正,方法如下:
式中:β是修正常数,这里取值1.3,修正后返回至(2),进行下一步迭代计算;
若且且其中,目标误差和或者迭代步数大于50000步,停止迭代,得到的模拟非高斯过程即为第i阶POD主坐标时程模拟结果,记为并返回至(1),模拟第i+1阶POD主坐标时程;
(7)重复上述(1)至(6),直至完成前N阶POD主坐标时程模拟,模拟主坐标时程分别为:
第五步:将模拟得到的主坐标时程与原本征向量组合形成模拟零均值风压场:
下面例举一个实施例:
第一步:图2为苍南发电厂干煤棚网架风洞试验模型测点布置图,根据该风洞试验数据,选取具有代表性的0°(无煤堆工况)荷载进行非高斯风压场模拟。截取中间30秒时程,减去各测点风压均值,得到各测点目标零均值风压场数据。
第二步:对所得目标零均值风压场进行本征正交分解,得到各阶主坐标时程和本征向量。
第三步:根据各阶本征向量对应的本征值大小,选取模拟的前100阶主坐标时程及对应的本征向量,前五阶本征向量等值线云图如图3所示。
第四步:对选取后的前100阶主坐标时程进行模拟:
(1)计算第i阶主坐标时程的目标非高斯过程功率谱、偏度和峰度,前5阶主坐标偏度和峰度如表1所示;
表1
(2)根据Hermite矩变换关系推导出第i阶主坐标时程自相关函数;
(3)对作快速傅里叶变换得到初始模拟高斯过程功率谱采用谐波叠加法得到一个模拟高斯过程计算其功率谱;
(4)采用Hermite矩变换方法,得到模拟非高斯过程计算其功率谱、偏度和峰度;
(5)引入第j次模拟误差评价参数
(6)若或或则对(3)中模拟高斯过程功率谱进行修正,返回至(2),进行下一步迭代计算;
若且且或者迭代步数大于50000步,停止迭代,得到的模拟非高斯过程即为第i阶主坐标时程模拟结果;返回至(1),模拟第i+1阶主坐标时程;
(7)重复上述(1)至(6),直至完成前100阶主坐标时程模拟,前五阶模拟主坐标时程与目标时程对比如图4所示。
第五步:将模拟得到的主坐标时程与原征向量组合形成模拟零均值风压场。脉风压系数均方根值较大的代表性测点G3(见图2中标注所示)模拟值与目标值的功率谱密度对比如图5所示。
如上所述,尽管参照特定的优选实施例已经表示和表述了本发明,但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下,可对其在形式上和细节上作出各种变化。
Claims (1)
1.一种大跨屋盖结构非高斯风荷载数值模拟方法,其特征在于,该方法包括以下步骤:
第一步:通过屋盖刚性模型风洞测压试验,获取风压场数据样本:
P(x,y,t)={p1(t)、p2(t)...pn(t)},计算各测点风压均值:得到目标零均值风压场P*(x,y,t)={p1 *(t);p2 *(t)...pn *(t)};
第二步:进行本征正交分解,目标零均值风压场P*(x,y,t)可以展开为:
式中:为第i阶POD主坐标时程,为第i阶本征向量,相应的本征值则记为
第三步:根据各阶本征向量的本征值大小,选取模拟的前N阶POD主坐标时程 及对应的本征向量
第四步:对选取后的前N阶POD主坐标时程进行模拟,模拟算法如下:
(1)计算第i阶POD主坐标时程的目标非高斯过程的功率谱、偏度和峰度,分别记为:
(2)根据Hermite矩变换关系,推导出第i阶POD主坐标时程自相关函数对应的高斯过程自相关函数
(3)对作快速傅里叶变换得到初始模拟高斯过程的功率谱采用谐波叠加法得到第j次模拟的高斯过程数据计算其功率谱,记为
(4)采用Hermite矩变换方法,第j次模拟得到非高斯过程数据计算其功率谱、偏度和峰度,分别记为:
(5)引入对第j次模拟结果的误差评价参数
<mrow>
<msubsup>
<mi>&epsiv;</mi>
<mn>1</mn>
<mi>j</mi>
</msubsup>
<mo>=</mo>
<mn>100</mn>
<msqrt>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
<mi>j</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>S</mi>
<mi>a</mi>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>S</mi>
<mi>a</mi>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</msqrt>
</mrow>
<mrow>
<msubsup>
<mi>&epsiv;</mi>
<mn>2</mn>
<mi>j</mi>
</msubsup>
<mo>=</mo>
<mn>100</mn>
<mo>|</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&gamma;</mi>
<mn>3</mn>
<mi>T</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&gamma;</mi>
<mn>3</mn>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msubsup>
<mi>&gamma;</mi>
<mn>3</mn>
<mi>T</mi>
</msubsup>
<mo>|</mo>
</mrow>
<mrow>
<msubsup>
<mi>&epsiv;</mi>
<mn>3</mn>
<mi>j</mi>
</msubsup>
<mo>=</mo>
<mn>100</mn>
<mo>|</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&gamma;</mi>
<mn>4</mn>
<mi>T</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&gamma;</mi>
<mn>4</mn>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msubsup>
<mi>&gamma;</mi>
<mn>4</mn>
<mi>T</mi>
</msubsup>
<mo>|</mo>
</mrow>
为目标非高斯过程功率谱密度函数与第j次模拟得到的非高斯过程功率谱密度函数之间的相对误差,fn为离散的频率点,为目标高阶统计量与模拟得到的非高斯过程高阶统计量之间的误差;
(6)若或或其中,目标误差和则对(3)中模拟高斯过程功率谱进行修正,方法如下:
<mrow>
<msubsup>
<mi>S</mi>
<mi>g</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>S</mi>
<mi>g</mi>
<mi>j</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<msubsup>
<mi>S</mi>
<mi>a</mi>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
<mi>j</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mi>&beta;</mi>
</msup>
</mrow>
式中:β是修正常数,这里取值1.3,修正后返回至(2),进行下一步迭代计算;
若且且其中,目标误差和或者迭代步数大于50000步,停止迭代,得到的非高斯过程数据即为第i阶POD主坐标时程模拟结果,记为并返回至(1),模拟第i+1阶POD主坐标时程;
(7)重复上述(1)至(6),直至完成前N阶POD主坐标时程模拟工作完成,模拟得到的主坐标时程数据分别为:
第五步:将模拟得到的主坐标时程与原风压场本征向量组合得到零均值风压场的模拟结果:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710270627.2A CN107194036A (zh) | 2017-04-24 | 2017-04-24 | 一种大跨屋盖结构非高斯风荷载数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710270627.2A CN107194036A (zh) | 2017-04-24 | 2017-04-24 | 一种大跨屋盖结构非高斯风荷载数值模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107194036A true CN107194036A (zh) | 2017-09-22 |
Family
ID=59872220
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710270627.2A Withdrawn CN107194036A (zh) | 2017-04-24 | 2017-04-24 | 一种大跨屋盖结构非高斯风荷载数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107194036A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108804838A (zh) * | 2018-06-15 | 2018-11-13 | 辽宁工程技术大学 | 一种复杂大跨度双曲屋盖的抗风设计方法 |
CN111027261A (zh) * | 2019-11-15 | 2020-04-17 | 四川大学 | 一种用于研究结构风激励响应的混合模拟试验方法 |
CN112131638A (zh) * | 2020-09-09 | 2020-12-25 | 石家庄铁道大学 | 大跨屋盖结构的风致动力特性类型判定方法及终端设备 |
CN112749476A (zh) * | 2020-11-26 | 2021-05-04 | 重庆交通大学 | 基于Piecewise-Johnson变换的非高斯风压模拟方法及***和存储介质 |
CN115034154A (zh) * | 2022-06-09 | 2022-09-09 | 北京航空航天大学 | 一种基于欠均匀离散测点数据的屋面风压重构算法 |
CN115935485A (zh) * | 2022-12-29 | 2023-04-07 | 北京建筑大学 | 矩形高层建筑非平稳横风向风荷载条件模拟方法 |
WO2023229115A1 (ko) * | 2022-05-27 | 2023-11-30 | 서울대학교 산학협력단 | 비대칭도를 이용한 시간이력풍하중 생성 방법 |
CN115034154B (zh) * | 2022-06-09 | 2024-07-09 | 北京航空航天大学 | 一种基于欠均匀离散测点数据的屋面风压重构方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1851436A (zh) * | 2006-05-31 | 2006-10-25 | 汕头大学 | 大跨度屋盖和超高层建筑结构风振响应的检测计算方法 |
CN103020471A (zh) * | 2012-12-27 | 2013-04-03 | 黑龙江大学 | 一种大跨度屋盖结构脉动风致响应计算的块里兹向量生成方法 |
US20150071978A1 (en) * | 2013-09-06 | 2015-03-12 | Alice Chang | Clothing and covering system with various functions |
-
2017
- 2017-04-24 CN CN201710270627.2A patent/CN107194036A/zh not_active Withdrawn
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1851436A (zh) * | 2006-05-31 | 2006-10-25 | 汕头大学 | 大跨度屋盖和超高层建筑结构风振响应的检测计算方法 |
CN103020471A (zh) * | 2012-12-27 | 2013-04-03 | 黑龙江大学 | 一种大跨度屋盖结构脉动风致响应计算的块里兹向量生成方法 |
US20150071978A1 (en) * | 2013-09-06 | 2015-03-12 | Alice Chang | Clothing and covering system with various functions |
Non-Patent Citations (1)
Title |
---|
冯鹤: "大跨开敞式干煤棚风荷载研究及非高斯风压场模拟", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108804838A (zh) * | 2018-06-15 | 2018-11-13 | 辽宁工程技术大学 | 一种复杂大跨度双曲屋盖的抗风设计方法 |
CN111027261A (zh) * | 2019-11-15 | 2020-04-17 | 四川大学 | 一种用于研究结构风激励响应的混合模拟试验方法 |
CN111027261B (zh) * | 2019-11-15 | 2023-09-26 | 四川大学 | 一种用于研究结构风激励响应的混合模拟试验方法 |
CN112131638A (zh) * | 2020-09-09 | 2020-12-25 | 石家庄铁道大学 | 大跨屋盖结构的风致动力特性类型判定方法及终端设备 |
CN112131638B (zh) * | 2020-09-09 | 2022-03-29 | 石家庄铁道大学 | 大跨屋盖结构的风致动力特性类型判定方法及终端设备 |
CN112749476A (zh) * | 2020-11-26 | 2021-05-04 | 重庆交通大学 | 基于Piecewise-Johnson变换的非高斯风压模拟方法及***和存储介质 |
CN112749476B (zh) * | 2020-11-26 | 2022-09-30 | 重庆交通大学 | 基于Piecewise-Johnson变换的非高斯风压模拟方法及***和存储介质 |
WO2023229115A1 (ko) * | 2022-05-27 | 2023-11-30 | 서울대학교 산학협력단 | 비대칭도를 이용한 시간이력풍하중 생성 방법 |
CN115034154A (zh) * | 2022-06-09 | 2022-09-09 | 北京航空航天大学 | 一种基于欠均匀离散测点数据的屋面风压重构算法 |
CN115034154B (zh) * | 2022-06-09 | 2024-07-09 | 北京航空航天大学 | 一种基于欠均匀离散测点数据的屋面风压重构方法 |
CN115935485A (zh) * | 2022-12-29 | 2023-04-07 | 北京建筑大学 | 矩形高层建筑非平稳横风向风荷载条件模拟方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107194036A (zh) | 一种大跨屋盖结构非高斯风荷载数值模拟方法 | |
CN107194625B (zh) | 基于神经网络的风电场弃风电量评估方法 | |
CN104102836B (zh) | 一种电力***快速抗差状态估计方法 | |
CN107422376B (zh) | 一种非平稳非高斯地震动时程模拟方法 | |
CN113189561B (zh) | 一种海杂波参数估计方法、***、设备及存储介质 | |
Abramov et al. | A new algorithm for low-frequency climate response | |
El Gharamti et al. | Comparing adaptive prior and posterior inflation for ensemble filters using an atmospheric general circulation model | |
CN109800449A (zh) | 一种基于神经网络的航空发动机压缩部件特性修正方法 | |
Do et al. | Demand estimation in water distribution systems: solving underdetermined problems using genetic algorithms | |
Li et al. | Simulation of stationary non-Gaussian stochastic vector processes using an eigenvalue-based iterative translation approximation method | |
CN114442153A (zh) | 一种近断层地震动拟合方法 | |
Tian et al. | Calibrating dynamic building energy models using regression model and Bayesian analysis in building retrofit projects | |
CN108169794A (zh) | 近断层人工地震波与规范反应谱匹配的调整方法 | |
CN104105049A (zh) | 一种减少传声器使用数量的房间冲激响应函数测量方法 | |
CN117933095A (zh) | 一种基于机器学习的地表发射率实时反演及同化方法 | |
CN103414671A (zh) | 基于单位基函数的2d-dpd迭代削减算法及应用 | |
CN105488351A (zh) | 一种生成移动心电信号的噪声模型的方法 | |
Tafazoli et al. | Performance improvement of GPS GDOP approximation using recurrent wavelet neural network | |
CN111277950B (zh) | 一种基于到达时差和到达频差的定位方法、装置及设备 | |
Lee et al. | Wind power scenario generation for stochastic wind power generation and transmission expansion planning | |
CN112749476B (zh) | 基于Piecewise-Johnson变换的非高斯风压模拟方法及***和存储介质 | |
Rios et al. | PEM-based Generation Adequacy Assessment of Power Systems with Non-Conventional Renewable Energy Sources | |
CN107888272A (zh) | 一种随机分布鲁棒波束成形器的构建及其计算方法 | |
CN104037756B (zh) | 一种含复杂电力设备模型的电力***稳定评估方法 | |
Ali et al. | Statistical downscaling of precipitation and temperature using gene expression programming |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20170922 |
|
WW01 | Invention patent application withdrawn after publication |