风时程生成程序技术说明

更新时间：2023-12-01 22:44:01 阅读量：教育文库文档下载

说明：文章内容仅供预览，部分内容可能不全。下载后的文档，内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的，是否完整无缺。

风的生成原因推荐度：
相关推荐

编写：陈学伟 dinochen1983, WindHist说明书

1 程序原理........................................................................................... 3

1.1 风荷载动力分析方法简介........................................................................ 3 1.2 风速时程模拟的AR法 ............................................................................ 4

1.2.1 AR模型 ............................................................................................................ 4 1.2.2 AR模型模拟风速时程的基本过程 ................................................................ 5

1.3 风时程生成程序实现................................................................................ 7 1.4 风时程生成程序特点................................................................................ 9 1.5 风时程生成程序局限性说明.................................................................. 10

2 参数说明.......................................................................................... 11

2.1 顺向脉动风速功率谱密度函数Sv(n) .................................................... 11 2.2 脉动风空间相干函数rij .......................................................................... 13 2.3 地面粗糙系数k（紊流度）................................................................... 14 2.4 平均风速v .............................................................................................. 14 2.5 风压力时程Fw(x,y,z,t) ......................................................................... 16 2.6 数值计算的参数...................................................................................... 17

3 操作说明......................................................................................... 18

3.1 制作空间点信息表格（*.csv） ............................................................. 18 3.2 导入表格及输入参数.............................................................................. 19 3.3 计算风时程.............................................................................................. 20 3.4 显示计算结果.......................................................................................... 20 3.5 输出时程结果及分析代码...................................................................... 21 3.6 接力SAP2000进行时程分析 ................................................................ 21 3.7 接力ETABS进行时程分析 ................................................................... 22 3.8 SAP2000与ETABS的分析代码例子 ................................................... 23

3.8.1 ETABS分析代码 ........................................................................................... 23

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 1 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

3.8.2 SAP02000分析代码： .................................................................................. 24

4 计算实例......................................................................................... 25

4.1 操作步骤.................................................................................................. 25 4.2 24层框架风振分析结果分析................................................................. 29

4.2.1 风速时程结果 ................................................................................................ 29 4.2.2 风振分析计算结果与按现行《荷载规范》得出的结果对比..................... 31 4.2.3 风振分析的顶点加速度计算与按《高钢规》手算结果对比..................... 32

5 关于风振时程分析的若干建议 .................................................... 34

5.1 分析参数设置.......................................................................................... 34 5.2 输出结果处理.......................................................................................... 34

6 参考文献......................................................................................... 36

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 2 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

程序原理

风荷载动力分析方法简介

风荷载是作用在结构上的重要动力荷载之一，尤其对于高层、高耸及大跨结构来说，设计中必须考虑风荷载的作用。计算高层、大跨、悬索桥以及塔架结构的动力风振相应的一个有效方法是Monte Carlo法。即根据某些既定的统计参数产生一系列的时程样本，再对每个样本函数进行线性或非线性的结构分析。通过对结构不同单元在样本函数下的时程响应的统计分析，计算整个结构是否安全。在结构特定点上的风速通常认为是稳定的随机过程，其特性可完全由谱密度函数确定。目前，通常对于风敏感复杂结构的风振响应分析和风振系数的求解主要有以下4种方法:（1）频域分析法、（2）时域分析法、（3）风洞试验法、（4）数值风洞技术。其中，本程序采用时域分析法对结构进行风荷载的动力分析。

时域分析法直接运用风洞试验的风压时程或计算机模拟的风压时程作用于屋面结构进行风振响应时程分析。然后通过动力计算得到结构的动力响应，统计结构动力响应从而算得结构的风振系数。采用时域分析法，可以考虑自然风的时间相关性和结构非线性影响，更精确地反映结构的耦合风振情况。

由于风荷载可分解为平均风和脉动风，如图1-1，对于脉动风部分，借助经验公式给出的各种功率谱函数，可以从频域和时域两方面来分析。

图1-1 脉动风速vf、瞬时风速v和平均风速的关系v

基于线性迭加的频域分析方法概念清晰、简便，因此在工程中应用广泛。但不能给出反应的相关函数、瞬态反应，不能进行非线性结构分析，且对于大跨结构，由于其频率密集，因此在计算中应取多少阶模态，怎样补偿模态截断的能量损失，仍是一个问题。而时域分析法与之相比较，其优点表现为:时域法能进行较精确的非线性分析；响应量值如位移、力或加速度的最大值可以直接求出；在缺乏实测或试验资料的情况下，各种简化计算方法可以与

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 3 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

精确的时域分析方法进行比较验证；通过时域的信息可以获得幅值域及频域的信息。因此，研究时域模型一直是随机过程模拟的重要内容，时域模拟就是要通过已知的频域信息重现时程样本。

主要的时域模拟方法有：（1）线性滤波法模拟、（2）谐波叠加法模拟、（3）逆Fourier

变换法模拟、（4）小波分析。

风速时程模拟的AR法

近年来，线性滤波法中的自回归（Auto-Regressive，简称AR）模型因其计算量小、速度快，广泛用于随机振动和时间系列分析中。该模型将均值为零的白噪声随机系列通过线性滤波器，使其输出为具有指定谱特征的平稳随机过程。本程序采用线性滤波法的自回归模型计算。

AR模型

任何平稳随机信号x(n)都可以看一成由白噪声w(n)激励一个因果稳定的可逆系统H(z)产生的输出。AR模型将均值为零的白噪声随机系列通过线性滤波器，使其输出为具有指定谱特征的平稳随机过程。

自回归(AR)模型的系统函数H(z)为:

H(z)?G1??aiz?ii?1p （1-1）

式中，G为系统增益，可取为1； a1,a2,?,ap为待定参数。

AR模型的系统函数只有极点没有零点，故AR模型又称为全极点模型。模型的阶由分母多项式的阶p决定。如果在白噪声w(n)激励下的输出为x(n)，则模型输入、输出关系的时域表达式为:

x(n)??aix(n?i)?Gw(n) （1-2）

i?1p上式为AR模型的差分方程。

AR模型的一个重要特性是输出的自相关函数具有递推特性，这个递推关系便可表

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 4 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

示为:

?p???aiRx(m?i),m?1,2,?,p?i?1 （1-3） Rx(m)??p??aR(i)?G2,m?0?ix?i?1?上式就是AR模型的正则方程，也叫Yule-Walker方程。

AR模型模拟风速时程的基本过程

在满足工程计算精度要求的前提下，可对风速时程作以下假定:（1）任意一点处平均风速不随时间改变；（2）脉动风速时程是零均值平稳随机过程；（3）风速时程间具有空间相关性，AR法模拟风速时程的基本过程。

采用AR法推广到模拟多维风速时程的技术，M个点空间相关脉动风速时程V(x，y，z，t)列向量的AR模型可表示为：

V(x,y,z,t)???kV(x,y,z,t?k?t)?N(t) （1-4）

k?1p式中：x?[x1,x2,?,xM]T，y?[y1,y2,?,yM]T，z?[z1,z2,?,zM]T，(xi,yi,zi)为空间第i点坐标，i=1，2，…M；p为AR模型阶数；?t是模拟风速时程的时间步长；?k为AR模型自回归系数矩阵，为M?M阶方阵，k=1，…，p；N(t)为独立随机过程向量；

N(t)?L?n(t) （1-5）

式中，n(t)??n1(t),?nM(t)?,ni(t)为均值为0、方差为1的正态分布随机过程，i=1，…，M。

随机风过程的协方差R与回归系数?之间的关系可写成矩阵形式：

T?R?R????N? （1-6）

?OP????I,?1,??P? （1-7）

RN?R0???kR?k?t? （1-8）

k?1pT风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 5 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

参数说明

顺向脉动风速功率谱密度函数Sv(n)

国内外一些学者采用不同的方法，对强风记录资料进行了研究，得出了各种风速谱的公式。第一类风速谱是对强风观测记录进行相关的分析，获得相关曲线和相关函数，再通过超低频滤波器，直接测出风速的功率谱曲线，拟合出风速谱的数学表达式。该方法避免了第一类方法在理论计算过程中的误差，所求的风速谱公式里意义清晰，明确。风压谱是由风速谱与空气导纳函数之积换算出来的，它们都是风速脉动频率n的函数，脉动风与相关函数，风速谱、风压谱、风压谱及其概率分布统计量之间的关系，如图2-1所示。

风速的频率分布风压的概率分布风速谱Sv(n) 空气导纳风压谱Sw(n) 低频滤波相关函数R(τ) 风速风速图2-1 脉动风与相关函数，风速谱、风压谱、风压谱及其概率分布统计量之间的关系

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 11 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

表2-1、两类风速功率谱密度函数表达式

类别发表者 A.G.Davenport Harris A类西安热工所风速谱的数学表达式 x21200n Sv(n)?4kv,x?n(1?x2)4/3v102102Sv(n)?4kv10x1800n ,x?n(2?x2)5/6v10?(lgx?0.61)2?1.611kv1021200n Sv(n)?exp??,x??n0.5408?v10?v3x21800nSv(n)??,x?24/33n(1?x)v x1400n或 Sv(n)?0.238?v2,x?n(1?x2)5/6vxny,x?,x?0.2 n(1?50x)5/3v1nySv(n)?0.26kv10?2/3,x?,x?0.2nxvSv(n)?200kv10?Shiotani Emil Simiu B类 Hino yk3/2()1?4??Wx210 Sv(?)?18.0kv10?,x?850.57825/63?(1?x)?v10Sv(n),Sv(?)———风速谱

n,?———分别为频率、圆频率

k———表征地面粗糙度的系数

v10———为离地面10m高度处的平均风速值

y、x———高度值、无量纲频率

表2-1中A类风速谱的特点为：谱形式与高度无关，紊流尺度沿高度不变；B类风速谱则相反，其特点是：不同高度处风速谱形式不同，随着高度的增加，谱的峰值减小，峰值频率下降。图2-2所示为不同脉动风速谱的比较。

Davenport（译名：达文波特）根据世界上不同地点、不同高度测得到90多次的强风记录，并假定水平阵风谱中的湍流积分尺度L沿高度不变，取常数值1200m，并取脉动风速谱为不同离地高度实测值的平均值，建立了经验数学表达式如表2-1所示。

我国规范及在风工程应用中一般采用Davenport脉动风速谱。由图2-2看出，Davenport谱比其它谱偏大，而谱值偏大的范围正好是风频率与结构物自振频率接近的地方，影响较大，

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 12 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

故Davenport风速谱可能会高估结构的动力响应，其结果可能会偏于保守，但是在结构抗风的设计角度而言，却提高了结构的安全度。

图2-2不同风速功率谱曲线对比图

脉动风空间相干函数rij

当结构上一点i的脉动动风压达到最大时，与i点距离为r的j点的脉动风压一般不会同时达到最大值，在一定的范围内，离开i点越远，脉动风压同时达到最大值的可能性越小，这种性质称为脉动风的空间相关性。

程序采用以下相关系数公式：

在顺风向，对于像高层建筑那样的高度和宽度方向的尺度需同时考虑的建筑物，一般考

虑水平(x，y方向)和竖向（z方向)的相关，对此，Davenport提出了指数形式的经验公式：

rij?Coh(r,n)?Rxyz(xi,xj,yi,yj,zi,zj,n)?e?c （2-1）

式（2-1）中

c?2?2nCx2(xi?xj)2?Cy(yi?yj)2?Cz2(zi?zj)2v(z)?v(z?) （2-2）

其中Cx，Cy，Cz分别是x，y，z三个方向的空间衰减系数，Davenport建议取值为

Cx?Cy?16,Cz?10

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 13 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

其它形式的相关系数表达式，这里不一一列举。

地面粗糙系数k（紊流度）

根据我国可靠指标的规定的数值，我国规范保证系数（峰因子）μ的取值在2.2（保证率

在98.61%）左右，并常以带有保证系数（峰因子）μ和脉动风压实测数据的脉动系数来表达。

如果采用Davenport谱，由于它是由实测得来的，因而根据Davenport谱可以求出脉动

系数、紊流度和风速根方差之间的关系式（详细推导请参考《结构风工程理论.规范.实践》）：

1?v10?f(z)??24k??24k?z2(z)?2?I(z) （2-3）

v1?v10I(z)?6k?6k?z2(z) （2-4）

v参考我国《荷载规范》附录G的风荷载说明（第166页），脉动系数?f是根据国内实测数据，并参考国外规范资料取：

?f?0.5?351.8(??0.16)1?z??1.8(??0.16)()?0.5?35?z2(z) （2-5） 10因此，按我国规范公式推导出地面粗糙系数k为：

k?0.002152?353.6(??0.16) （2-6）

式中，?为地面粗糙度系数，按《荷载规范》由A、B、C、D类地面可得?为0.12、0.16、0.22、0.30。按式可得各种地面的地面粗糙系数k：

表2-2 不同类地面的地面粗糙系数k取值

地面粗糙度类别地面粗糙系数k

A 0.00129 B 0.00215 C 0.00464 D 0.01291 平均风速v

根据《荷载规范》所述，在大气边界内，风速随离地面高度而增大。当气压场随高度不变是，风速随高度增大的规律，主要取决于地面粗糙度和温度垂直梯度。通常认为在离地面高度为300~500m时，风速不再受地面粗糙度的影响，也即达到所谓“梯度风速”，该高度称

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 14 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

之梯度风高度。地面粗糙度等级低的地区，其梯度高度比等级高的地区低。下表是不同地面情况的梯度风高度：

表2-3 不同类地面的梯度风高度HG（m）取值

地面粗糙度类别 A B C D 梯度风高度HG（m） 300 350 400 450

根据地面粗糙度指数及梯度风高度，即可得到出风压高度变化系数如下：

z0.24)10z?zB?1.000()0.3210 （2-7） z?zC?0.616()0.4410z?zD?0.318()0.6010?zA?1.379(由风压与风速的关系式如式（2-8）得平均风速高度变化系数为公式（2-9）：

0.0120182v2??v?v?kN/m2 （2-8）

2g2?9.81600?2z0.12?1.379()v10z?zBv?1.000()0.1610 （2-10） z?zCv?0.616()0.2210z?zDv?0.318()0.3010?zA算例2-1：结构位于地面粗糙度等级为C的地区，基本风压为0.5kN/m，位于10m和20m处的空间点的平均风速为多少？

解：v10?1600?10?1600?0.616?0.5?22.2m/s

2?z??20?v20?v10???22.2????10??10??0.22?25.86m/s。

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 15 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

风压力时程Fw(x,y,z,t)

通过AR法计算可得到空间点的模拟总风速时程V(x,y,z,t)?v(z)?v(x,y,z,t)，按《荷载规范》引入空间点的体型系数可得风压动力荷载计算公式：

????V2(x,y,z,t)FW(x,y,z,t)??(??sxAx???syAy???sxAz) （2-11）

1600式中，?sx，?sy，?sz分别为空间点yz受风面、xz受风压和xy受风面体形系数，Ax，Ay，

Az分别为空间点yz受风面、xz受风面和xy受风面的面积。如图2-3所示。

图2-3 空间点的体型系数、受风面积与风压力的关系

算例2-2：已知空间点i处x方向受风面积为Ax?36m2，y方向受风面积为Ay?27m，z方向没有受风面，x方向风荷载作用下，体型系数为?x风荷载作用下，体型系数为?x风速为V(t)，风压力为多少？

windy2windx?1.4,?ywindx?0.2，y方向

?0.1,?ywindy?1.3，现风荷载角度是与x方向成30°，

V2(t)V2(t)F(t)??(??sxAx)??(36?1.4?cos30??36?0.1?sin30?)解：Wx 16001600?0.0284?V2(t)V2(t)V2(t)FWy(t)??(??syAy)??(27?1.3?sin30??27?0.2?cos30?) 16001600?0.0139?V2(t)风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 16 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

???V2(t)程序将部分作为风压时程文本输出，将(??sxAx???syAy???sxAz)作为静

1600力荷载输入SAP2000或ETABS作用于结构，再将这个静力荷载通过时程放大，实现风压力时程的加载。

数值计算的参数

运算阶次p是AR模型的阶次，详细请查看1.2.2的说明。

计算时间步长?t是风速时程的时间系列的单位，为了使计算稳定建议取值为0.1~0.25s。计算步数tsteps，风时程总时间为tsteps??t，自回归模型中建议取值不大于3000步。缓冲步数是在输入时程文本时为了使风压力时程开始阶段不产生荷载的突然施加使结构的加速度产生突变，在施加风压力时，荷载从零增加到初始风压力时程荷载的步数，建议缓冲步数大于30步。

在AR模型计算时，产生R矩阵之前必须通过式（1-10）对自功率谱密度进行数值积分，程序积分方法采用可控的复化梯形积分公式。其中积分的积分域与积分段数可供输入可以控制精度与计算速度。

fa：频区始值（单位是Hz）； fb：频始终值（单位是Hz）； n：复化梯形积分的分段数。

建议取值为：fa=0.001Hz，fb=10Hz，n=500。积分参数示意图如下：

FF?Sij(f)cos?2?f???

n图2-4 积分参数示意图

fbf

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 17 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

操作说明

制作空间点信息表格（*.csv）

打开程序目录下的sample.xls表格文件，这是个空间点信息的标准输入格式如表3-1：

表3-1标准空间点信息输入格式表格

Num 54 55 56 57 58

表格表头的信息内容如下：

Num——SAP2000点号；ETABS层号； Pnt——SAP2000点号；ETABS平面点号； x，y，z——空间点坐标，单位是（mm）；

Ax，Ay，Az——空间点三个方向受风面积，单位是（m2）；

μxx——x方向风荷载作用下，x方向受风压的体型系数，（表示风压与风速方向一致）； μxy——y方向风荷载作用下，x方向受风压的体型系数，（表示风压与风速方向一致）； μyx——x方向风荷载作用下，y方向受风压的体型系数，（表示风压与风速方向一致）； μyy——y方向风荷载作用下，y方向受风压的体型系数，（表示风压与风速方向一致）； μz——z方向受风压的体型系数，表示风压与风速方向一致。体型系数定义请参考图3-1所示。

Csv表格可以通过Microsoft Office Excel来编辑制作后，另存为*.csv文件即可。

x 0 0 0 0 0 y 0 0 0 z Ax Ay μxx μyx μxy μyy Az 36 36 36 36 36 1.3 1.3 1.3 1.3 1.3 0 0 0 0 0 0 0 0 0 0 1.3 1.3 1.3 1.3 1.3 0 0 0 0 0 μz 0 0 0 0 0 Pnt 54 55 56 57 58 3000 36 6000 36 9000 36 0 12000 36 0 15000 36 风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 18 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

μsyx=-0.6μsyy=+0.8Windxμsyx=+0.8yxμsyx=-0.6图3-1 体型系数取值示意图

μsyy=+0.8μsyy=-0.5

导入表格及输入参数

图3-2程序主界面图

图3-2为程序主界面。分为菜单栏、工具栏、参数栏、进度条、表格和图形区。

制作好空间点的信息表格后，按工具栏中按钮，导入CSV表格。导入成功后，

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 19 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

表格区自动更新CSV表格。

然后在参数栏输入计算参数、空间相关性参数和风荷载信息，其中参数的意义请参加第

2部分的参数说明。其中风荷载水平风角度是与x方向的顺时针夹角，单位是°。

计算风时程

在检查输入信息无误后，按工具栏中按钮，程序自动批量n个风时程，在计算

过程中可以通过观察进度条来检查程序计算时度，如图3-3。

图3-3计算进度显示

显示计算结果

在计算完成后，选择参数栏底部的显示图像的选框如图3-4，可实时查看全部风时程的计算结果，包括脉动风速时程、总风速时程和总风压时程。

图3-4图像显示对话框

图3-5风时程显示图像

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 20 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

如图3-5，图形显示区显示内容是红色线代码时程值，蓝色线是平均值线，左边轴线上从上往下分别是最大正值，中值，最大负值或零。蓝线的右边数据是平均值，水平刻度是时间（单位是sec），纵向刻度是风速（单位是m/s）或风压（单位是kN/m2）。

点击

可将图形区显示的曲线保存为bmp文件。

输出时程结果及分析代码

在计算时程完成后，点击工具栏中的按钮，弹出保存文件对话框，确定文件保

存文件夹以后，填写文件名后，程序自动保存全部风速时程，风压时程，ETABS时程分析代码与SAP2000时程分析代码在同一文件夹里。文件命名为：

etabs.e2k——ETABS时程分析代码 sap2000.s2k——SAP2000时程分析代码

VV（a）_b.dat——第a批点号（层号）为b的风速时程 W（a）_b.dat——第a批点号（层号）为b的风压时程

注：在SAP2000默认时程分析的模态分析工况取“Modal”工况。

接力SAP2000进行时程分析

打开SAP2000软件。

点击菜单中的File/Imporr/sap2000 v8/v9/v10.s2k Text File。弹出对话框如图3-6所示。

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 21 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

图3-6 SAP2000导入对话框

选择Add to existing model这一项，选择3.5步导出的sap2000.s2k文件即可。导入成功后便可以在sap2000下进行结构批量风振时程分析。

接力ETABS进行时程分析

1．打开ETABS软件；

2．点击菜单中的File/Export/Save model as Etabs.e2k Text file，成功导出未加风时程的

ETABS的e2k文本文件；

3．打开输出的e2k文本文件，在文件中找到代码如下：

$ LOG STARTCOMMENTS 4．这代码前，空白处粘贴风时程输出的ETABS.e2k的全部代码进去，保存文件退出； 5．点击菜单中的File/Import/Etabs.e2k Text file，把修改后的e2k文件导入ETABS，检查模

型无误后可以进行结构的批量风振时程分析。

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 22 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

SAP2000与ETABS的分析代码例子 ETABS分析代码

$ FUNCTIONS FUNCTION \ FUNCTYPE \ FILE \风速时程研究_0314\\haha\\W_(1)54.dat\ DATATYPE \ DT 0.10 FUNCTION \ POINTSPERLINE 1 FORMAT \FUNCTION \ FUNCTYPE \ FILE \风速时程研究_0314\\haha\\W_(1)55.dat\ DATATYPE \ DT 0.10 FUNCTION \ POINTSPERLINE 1 FORMAT \FUNCTION \ FUNCTYPE \ FILE \风速时程研究_0314\\haha\\W_(1)56.dat\ DATATYPE \ DT 0.10 FUNCTION \ POINTSPERLINE 1 FORMAT \ FUNCTION \ FUNCTYPE \ FILE \风速时程研究_0314\\haha\\W_(1)57.dat\ DATATYPE \ DT 0.10 FUNCTION \ POINTSPERLINE 1 FORMAT \FUNCTION \ FUNCTYPE \ FILE \风速时程研究_0314\\haha\\W_(1)58.dat\ DATATYPE \ DT 0.10 FUNCTION \ POINTSPERLINE 1 FORMAT \FUNCTION \ FUNCTYPE \ FILE \风速时程研究_0314\\haha\\W_(1)59.dat\ DATATYPE \ DT 0.10 FUNCTION \ POINTSPERLINE 1 FORMAT \ $ STATIC LOADS LOADCASE \ TYPE \ SELFWEIGHT 0 LOADCASE \ TYPE \ SELFWEIGHT 0 LOADCASE \ TYPE \ SELFWEIGHT 0 LOADCASE \ TYPE \ SELFWEIGHT 0 LOADCASE \ TYPE \ SELFWEIGHT 0 LOADCASE \ TYPE \ SELFWEIGHT 0 $ POINT OBJECT LOADS POINTLOAD \ \ TYPE \ LC \ FX 46.8000 FY 0.0000 FZ 0.0000 POINTLOAD \ \ TYPE \ LC \ FX 46.8000 FY 0.0000 FZ 0.0000 POINTLOAD \ \ TYPE \ LC \ FX 46.8000 FY 0.0000 FZ 0.0000 POINTLOAD \ \ TYPE \ LC \ FX 46.8000 FY 0.0000 FZ 0.0000 POINTLOAD \ \ TYPE \ LC \ FX 46.8000 FY 0.0000 FZ 0.0000 POINTLOAD \ \ TYPE \ LC \ FX 46.8000 FY 0.0000 FZ 0.0000 $ TIME HISTORY CASES THCASE \ TYPE \ NSTEPS 1024 DTOUT 0.10 DAMP 0 THCASE \ RFTOL .00001 RETOL .00001 MAXITER 100 MINITER 2 CONVFACT 1 THCASE \ LOAD \ FUNC \ SF 1 THCASE \ LOAD \ FUNC \ SF 1 THCASE \ LOAD \ FUNC \ SF 1 THCASE \ LOAD \ FUNC \ SF 1 THCASE \ LOAD \ FUNC \ SF 1 THCASE \ LOAD \ FUNC \ SF 1

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 23 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

SAP02000分析代码

TABLE: \ ProgramName=SAP2000 Version=10.0.7 ProgLevel=\ LicenseOS=Yes LicenseSC=Yes LicenseBR=Yes LicenseHT=No CurrUnits=\ SteelCode=\ ConcCode=\ AlumCode=\ ColdCode=AISI-ASD96 StiffCase=None TABLE: \ Name=W(1)54 Time=0.00 Value=0.000000 HeaderLines=0 PrefixChars=0 PtsPerLine=1 DataType=\ FormatType=Free Interval=0.10 FileName=\风速时程研究_0314\\haha\\W(1)_54.dat\ Name=W(1)55 Time=0.00 Value=0.000000 HeaderLines=0 PrefixChars=0 PtsPerLine=1 DataType=\ FormatType=Free Interval=0.10 FileName=\风速时程研究_0314\\haha\\W(1)_55.dat\ Name=W(1)56 Time=0.00 Value=0.000000 HeaderLines=0 PrefixChars=0 PtsPerLine=1 DataType=\ FormatType=Free Interval=0.10 FileName=\风速时程研究_0314\\haha\\W(1)_56.dat\ Name=W(1)57 Time=0.00 Value=0.000000 HeaderLines=0 PrefixChars=0 PtsPerLine=1 DataType=\ FormatType=Free Interval=0.10 FileName=\风速时程研究_0314\\haha\\W(1)_57.dat\ Name=W(1)58 Time=0.00 Value=0.000000 HeaderLines=0 PrefixChars=0 PtsPerLine=1 DataType=\ FormatType=Free Interval=0.10 FileName=\风速时程研究_0314\\haha\\W(1)_58.dat\ Name=W(1)59 Time=0.00 Value=0.000000 HeaderLines=0 PrefixChars=0 PtsPerLine=1 DataType=\ FormatType=Free Interval=0.10 FileName=\风速时程研究_0314\\haha\\W(1)_59.dat\TABLE: \ LoadCase=F54 DesignType=Other SelfWtMult=0 LoadCase=F55 DesignType=Other SelfWtMult=0 LoadCase=F56 DesignType=Other SelfWtMult=0 LoadCase=F57 DesignType=Other SelfWtMult=0 LoadCase=F58 DesignType=Other SelfWtMult=0 LoadCase=F59 DesignType=Other SelfWtMult=0 TABLE: \ Joint=54 LoadCase=F54 CoordSys=GLOBAL F1=46.8000 F2=0.0000 F3=0.0000 M1=0 M2=0 M3=0 Joint=55 LoadCase=F55 CoordSys=GLOBAL F1=46.8000 F2=0.0000 F3=0.0000 M1=0 M2=0 M3=0 Joint=56 LoadCase=F56 CoordSys=GLOBAL F1=46.8000 F2=0.0000 F3=0.0000 M1=0 M2=0 M3=0 Joint=57 LoadCase=F57 CoordSys=GLOBAL F1=46.8000 F2=0.0000 F3=0.0000 M1=0 M2=0 M3=0 Joint=58 LoadCase=F58 CoordSys=GLOBAL F1=46.8000 F2=0.0000 F3=0.0000 M1=0 M2=0 M3=0 Joint=59 LoadCase=F59 CoordSys=GLOBAL F1=46.8000 F2=0.0000 F3=0.0000 M1=0 M2=0 M3=0 TABLE: \ Case=WIND_HIST(1) Type=LinModHist InitialCond=Zero ModalCase=MODAL RunCase=Yes TABLE: \ Case=WIND_HIST(1) HistoryType=Transient OutSteps=1024 StepSize=0.10 DampingType=Constant ConstDamp=.05 TABLE: \ Case=WIND_HIST(1) LoadType=LoadPattern LoadName=F54 Function=W(1)54 LoadSF=1 TimeFactor=1 ArrivalTime=0 Case=WIND_HIST(1) LoadType=LoadPattern LoadName=F59 Function=W(1)59 LoadSF=1 TimeFactor=1 ArrivalTime=0 Case=WIND_HIST(1) Type=LinModHist InitialCond=Zero ModalCase=MODAL RunCase=Yes END TABLE DATA 风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 24 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

计算实例

下面以如图4-1所示的24层框架为实例，介绍用WINDHIST PRODUCER与SAP2000、ETABS分析结构的风振反应。详细参数如表4-1所示。

图4-1 24层框架结构平面布置图

表4-1 24层框架结构各种参数层数层高跨长梁截面混凝土材料

24 3000 6000 B300X800 C30 楼板 SLAB100 柱截面 C800X800 梁截面 B300X800 恒载活载 3.0kN/m2 3.5kN/m2 操作步骤

1．打开SAP2000进行结构的建模及恒活载的施加，完成建模后的截图如图4-2所示。

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 25 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

图4-2 结构未加风时程的模型

2．点选每层加载点，模型取中心点为加载点，选取后按Ctrl+C，复制结点到新建的Excel

表格，会得到如下图所示的信息。

图4-3 粘贴到Excel表的sap2000空间点的信息

3．参考3.1所述的方法编辑如图4-3的Excel表格，得到空间点信息输入表格形式，在这

个过程中需要计算空间点的受风面，输入查阅规范或参考风洞试验得到的体型系数。通过Excel的计算与编辑，输出*.csv文件，如下表所示。

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 26 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

表4-2 空间点信息输入表格

Num 302 303 304 305 306 307 308 309 310 x 0 0 0 0 0 0 0 0 0 y 0 0 0 0 0 0 0 0 0 z 3000 6000 9000 12000 15000 18000 21000 24000 27000 Ax 72 72 72 72 72 72 72 72 72 Ay 72 72 72 72 72 72 72 72 72 μxx 1.3 1.3 1.3 1.3 1.3 1.3 1.3 1.3 1.3 μyx 0 0 0 0 0 0 0 0 0 μxy 0 0 0 0 0 0 0 0 0 μyy 1.3 1.3 1.3 1.3 1.3 1.3 1.3 1.3 1.3 Az 0 0 0 0 0 0 0 0 0 μz 0 0 0 0 0 0 0 0 0 Pnt 302 303 304 305 306 307 308 309 310 ? ? 324 325 ? ? 0 0 ? ? 0 0 ? ? 69000 72000 ? ? 72 36 ? ? 72 36 ? ? 1.3 1.3 ? ? 0 0 ? ? 0 0 ? ? 1.3 1.3 ? ? 0 0 ? ? 0 0 ? ? 324 325

4．打开WINDHIST PRODUCER V1.0程序，按第3部分的操作说明所示，计算风时程。计

算风时程的相关参数如下表所示，其它按默认。

表4-3 风时程生成计算的相关参数批数阶数步长 10 4 2基本风压 0.5kN/m 地面分类 C类 1024 100 0.1s 计算步数频区终值 5Hz 缓冲步数

5．点击工具栏的计算时程按钮，计算完成后，程序自动绘制风时程曲线如下图。

风时程生成程序WINDHIST PRODUCER V1.0技术说明书第 27 页共 36 页

编写：陈学伟 dinochen1983, WindHist说明书

图4-4计算完成时程序的截图

6．计算完成后，导出计算结果与SAP2000的分析代码，点击工具栏中的导出文本，程序自

动生成SAP2000和ETABS的分析代码。创建目录如名为“windhist”，将计算结果存放于文件夹里。

7．打开SAP2000的计算模型，按3.6所述方法导入SAP2000的时程分析代码，注意需选

择Add to existing model这一项。导入完成后，弹出如图4-5对话框表示导入成功。

图4-5 SAP2000导入报告对话框