2014年11月29日,珠海生态树项目完成。上一篇文章未附带该程序,现附带如下:

%% 程序信息
% 程序内容:人行道风环境处理程序
% 参考资料:《建筑结构荷载规范》
% 编译单位:湖南大学建筑安全与节能教育部重点实验室
% 创建时间:2014年11月29日 
%% 初始化
clear,clc,close all
%% 读取风压数据
for m=8:-1:1
    for n=3:-1:1
        %读取风压数据
        press=load([pwd,'\风压数据\',num2str(n),'-1-5\',num2str(45*(m-1)),'.txt']);
        pra(:,m,5*n-4:5*n)=press(:,2:2:10)-press(:,3:2:11);
%         pra(:,m,5*n-4:5*n)=abs(press(:,2:2:10)-press(:,3:2:11));
    end
end
%% 风压转换风速
%已知拟合系数
p=load('p_code.txt');
speed_o=zeros(size(pra));
for i=5:-1:1
speed_o(:,:,i:5:end)=polyval(p(i,:),pra(:,:,i:5:end));
end
%% 风速风向插值
speed_c=(speed_o+[speed_o(:,2:end,:),speed_o(:,1,:)])/2;
speed_a=repmat(zeros(size(speed_o)),1,2);
speed_a(:,1:2:end,:)=speed_o;
speed_a(:,2:2:end,:)=speed_c;
%% 求解均值脉动
spe_mea=mean(speed_a);      
spe_mea=reshape(spe_mea(:,:),16,15)'; %风速均值
spe_std=std(speed_a);
spe_std=reshape(spe_std(:,:),16,15)'; %风速脉动
spe_max=max(speed_a);
spe_max=reshape(spe_max(:,:),16,15)'; %风速最值</pre>

标签: 风环境

评论已关闭