一种分区计算空间相关性的脉动风速时程输入方法与流程

本发明涉及非稳态数值风洞模拟技术领域，尤其是一种分区计算空间相关性的脉动风速时程输入方法。

背景技术：

建筑结构所受到的自然风主要包含长周期和短周期两种成分。在进行建筑结构风振响应分析计算时，需要考虑短周期成分对结构的脉动作用。短周期的脉动风一般被看作是平稳随机过程，当采用非稳态数值风洞模拟技术分析结构所受的带脉动成分的风荷载时程作用时，需要在数值风洞的入口边界输入带脉动成分的风速时程，该风速时程通常可采用人工模拟方法计算得到。为了确保输入的风速时程能满足自然风的随机性和统计特性，需要使模拟的风与自然风基本特性尽可能接近。其中，空间相关性是脉动风的基本特性之一，它是指当结构上一点风压达到最大值时，在一定范围内离该点越远处的风荷载同时达到最大值的可能性越小。我国现行《建筑结构荷载规范》采用了shiotani提出的相干函数来考虑空间相关性。考虑到数值风洞模拟中，建筑物最大迎风面积比流域横截面积一般不大于3％，因而对整个入口边界的风速输入考虑脉动分量会导致一定的电算资源浪费；同时，对每两个位置的风速时程考虑空间相关性也会大大增加风速时程模拟的工作量和求解时间，降低分析效率。

技术实现要素：

本发明的目的在于提供一种分区计算空间相关性的脉动风速时程输入方法，该方法不仅能使输入的风速时程能够满足规范建议的风速时程的相关性，而且能减少风速时程数据模拟的工作量和求解时间，从而提升整个结构表面风荷载模拟及风振响应分析的工作效率。

实现本发明目的的具体技术方案是：

一种分区计算空间相关性的脉动风速时程输入方法，该方法包括以下具体步骤：

步骤1：建立数值风洞计算域，并将数值风洞入口边界划分为平均风速区a和带脉动风速区d两部分，并根据带脉动风速区d的大小进行划分，形成d1,d2,……,dn共n(n>1)个带脉动风速区分区；

步骤2：对整个数值风洞计算域进行网格划分，使入口边界内形成多个网格，然后输出平均风速区a和带脉动风速分区di(i＝1,2,……,n)的网格中心点坐标数组pa[x,z]和pdi[x,z]；

步骤3：根据每个带脉动风速区分区di的网格中心点坐标数组pdi[x,z]，运用线性滤波法的ar模型，采用davenport顺风向脉动风速功率谱和shiotani提出的相干函数，仅对带脉动风速区分区di内的所有的网格中心点两两考虑空间相关性，生成与各个网格中心点对应的脉动风速时程vdi；

步骤3中，假设带脉动风速区分区di内的网格数为m，则所述ar模型生成的m个空间相关点脉动风速的时间序列函数定义为：

式中，p为ar模型阶数；△t为时间步长；pdi[x,z]为带脉动风速区分区di的网格中心点坐标数组，pdi[xjj,zjj]为带脉动风速区分区di第jj个网格中心点坐标，jj＝1,2,……,m；ψk为ar模型的m×m阶自回归系数矩阵，可由下式得到：

式中，r(j△t)是由时滞j△t的m×m的协方差矩阵，其中的元素为

spq(f)在p＝q时为脉动风速自谱密度函数，在p≠q时为脉动风速互谱密度函数，可由脉动风速自谱函数sn(f)与相干函数γpq(f)确定：

所述脉动风速功率谱sn(f)用函数定义为：

式中，f为脉动风频率；为标准高度为10m处的平均风速；k为与地面粗糙度有关的常数，计算公式为k＝0.0021522×35^36(α-0.16)，其中α为地面粗糙度系数。

所述shiotani提出的相干函数定义为：

式中，lx＝50m；lz＝60m；p＝1,2,……,m；q＝1,2,……,m。

n(t)为均值为0方差为1的均匀分布随机数，且与矩阵rn相关，rn可由下式求得：

求解出rn后可由下式解出n(t)：

n(t)＝ln(t)

式中，n(t)＝[n¹(t),…,n^m(t)]^t为m个独立的满足零均值、单位方差的随机变数。l满足rn＝ll^t，即l为rn经过乔累斯基分解后的下三角矩阵；

步骤4：根据建筑所处地貌确定平均风速剖面，结合网格中心点的z向坐标，计算得到平均风速区a和各个带脉动风速区分区di的平均风速数组和

步骤4中，所述平均风速剖面中某网格中心点的速度函数定义为：

式中，为标准高度为10m处的平均风速；α为地面粗糙度系数；

步骤5：对每个分区内所有网格逐个按时间步输入风速时程，平均风速区a输入各带脉动风速区分区di输入然后进行数值风洞模拟；

步骤6：完成求解后，查看平均风速区a控制计算域内的风速、风向，若方向与风速输入方向一致且大小不变，则表明平均风速区a的大小合适，不影响建筑在垂直于风速输入方向的风荷载的发展；若有变化，则需增大带脉动风速区d的大小，然后重复步骤2～6。

本发明的有益效果是：可在实现建筑表面脉动风荷载模拟的基础上，减少考虑脉动成分的入口大小，从而减少分析的复杂性，提升数值风洞计算效率；对带脉动风速区进行分区并对各分区进行风速时程模拟，可大量减少计算空间相关性的矩阵大小，节省风速时程模拟时间；对每个分区中的网格点进行风速时程模拟和输入，脉动风的分布可更接近实际情况。

附图说明

图1是发明流程图；

图2是本发明实施例的轴测图；

图3是本发明实施例的左视图；

图4是图2的i节点详图；

图5是本发明实施例的平均风速区控制计算域示意图。

图中：1.数值风洞计算域，2.入口边界，3.平均风速区a，4.带脉动风速区d，5.建筑，6.带脉动风速区分区，7.网格，8.网格中心点，9.平均风速区控制计算域

具体实施方式

下面结合附图对本发明作进一步的说明。本发明的一种分区计算空间相关性的脉动风速时程输入方法，包括：

1)如图2、3所示，首先建立数值风洞计算域1，并将数值风洞入口边界2划分为平均风速区a3和带脉动风速区d4两部分，并根据带脉动风速区d4的大小进行划分，形成d1,d2,……,dn共n(n>1)个带脉动风速区分区6。

2)如图4所示，对整个数值风洞计算域1进行网格划分，使入口边界2内形成多个网格7，然后输出平均风速区a3和带脉动风速分区di(i＝1,2,……,n)6的网格中心点8坐标数组pa[x,z]和pdi[x,z]。

3)根据每个带脉动风速区分区di6的网格中心点8坐标数组pdi[x,z]，运用线性滤波法的ar模型，采用davenport顺风向脉动风速功率谱和shiotani提出的相干函数，仅对带脉动风速区分区di6内的所有的网格中心点8两两考虑空间相关性，生成与各个网格中心点8对应的脉动风速时程vdi。

步骤3)中，假设带脉动风速区分区di6内的网格7数为m，则所述ar模型生成的m个空间相关点脉动风速的时间序列函数定义为：

式中，p为ar模型阶数；△t为时间步长；pdi[x,z]为带脉动风速区分区di6的网格中心点8坐标数组，pdi[xjj,zjj]为带脉动风速区分区di6第jj个网格中心点8坐标，jj＝1,2,……,m；ψk为ar模型的m×m阶自回归系数矩阵，可由下式得到：

式中，r(j△t)是由时滞j△t的m×m的协方差矩阵，其中的元素为

spq(f)在p＝q时为脉动风速自谱密度函数，在p≠q时为脉动风速互谱密度函数，可由脉动风速自谱函数sn(f)与相干函数γpq(f)确定：

所述脉动风速功率谱sn(f)用函数定义为：

式中，f为脉动风频率；为标准高度为10m处的平均风速；k为与地面粗糙度有关的常数，计算公式为k＝0.0021522×35^36(α-0.16)，其中α为地面粗糙度系数。

所述shiotani提出的相干函数定义为：

式中，lx＝50m；lz＝60m；p＝1,2,……,m；q＝1,2,……,m。

n(t)为均值为0方差为1的均匀分布随机数，且与矩阵rn相关，rn可由下式求得：

求解出rn后可由下式解出n(t)：n(t)＝ln(t)

式中，n(t)＝[n¹(t),…,n^m(t)]^t为m个独立的满足零均值、单位方差的随机变数。l满足rn＝ll^t，即l为rn经过乔累斯基分解后的下三角矩阵。

4)根据建筑所处地貌确定平均风速剖面，结合网格中心点8的z向坐标，计算得到平均风速区a3和各个带脉动风速区分区di6平均风速数组和

步骤4)中，所述平均风速剖面中某网格中心点8的速度函数定义为：

式中，为标准高度为10m处的平均风速；α为地面粗糙度系数。

5)对各个区域内所有网格7逐个按时间步输入风速时程，平均风速区a3输入各带脉动风速区分区di6输入然后进行数值风洞模拟。

6)完成求解后，如图5所示，查看平均风速区控制计算域9内的风速、风向，若方向与风速输入方向一致且大小不变，则表明平均风速区a3的大小合适，不影响建筑5在垂直于风速输入方向的风荷载的发展；若有变化，则需增大带脉动风速区d4的大小，然后重复步骤2～6。

本发明在实现建筑表面脉动风荷载模拟的基础上，减少考虑脉动成分的入口大小，从而减少分析的复杂性，提升数值风洞计算效率；对带脉动风速区进行分区并对各分区进行风速时程模拟，可大量减少计算空间相关性的矩阵大小，节省风速时程模拟时间。例如，不进行分区时，对整个区域的n个网格需要求解的矩阵维数为n，若分成3个网格数量一致的分区，则矩阵维数为n/3，对3个n/3维的矩阵进行求解所消耗的时间将远少于对1个n维矩阵求解的时间，即风速时程的模拟时间将大大减少；对每个分区中的网格进行风速时程模拟和输入，由于生成的网格尺寸一般相对较小，输入的脉动风的分布可更接近实际情况。