Matlab谐波叠加法风速时程生成代码

【资源介绍】

很多同学在做结构风致响应时，需要根据结构项目情况生成脉动风风速时程曲线，然后在转为风压时程曲线，最后转为节点力时程曲线，导入ANSYS或者ABAUQUS等商用有限元软件进行计算。

本代码基于Matlab,采用谐波叠加法法生成脉动风时程曲线，功率谱Davenport谱，采用的理论在任何关于风振的文献中均可查阅，这里便不在给出具体文献。

该代码需要用户确认的参数有如下：

1、风速生成点的离地高度，单位m

2、采样的点数，通常可取500~1000

3、时间间隔，通常可取0.1~0.2s

4、粗糙度类别，根据《荷载规范》取A、B、C、D

5、基本风压值，根据《荷载规范》取值，其中如果要生成竖向风速，一般取竖向基本风压为水平基本风压的0.18倍。

这里以某项目某个节点的结果进行演示。

项目情况：

1、生成点位离地标高26.5m；2、采样点数为500；3、时间间隔为0.1s；4、地面粗糙度为C类；5、基本风压为 0.6 KN/m^2

采用Matlab生成的风速曲线以及对应的功率谱如下所示：

从图中可见，本代码生成的风速 功率谱和目标功率谱 契合度较高~

【资源预览】

为方便各位同学方便调用，本资源含有三个代码，如下所示：

其中第一个代码为复核单个节点风速与功率谱用；

第二个代码则直接是一个matlab函数，直接使用函数即可批量生成风速时程，函数参数如下所示。

第三个则为一个通用的批量生成节点风速时程的matlab代码，采用的前处理软件为ANSYS，其他软件皆可通用。

【资源下载】