珠海项目的风压处理程式,特点在于对测压数据整合、预处理、插值点的添加和处理方式。

%% 程序信息
% 程序内容:风压系数的计算
% 参考资料:《建筑结构荷载规范》
% 编译单位:湖南大学建筑安全与节能教育部重点实验室
% 创建时间:2014年11月25日
%% 初始化
clear,clc,close all
%% 计算平均风压系数以及脉动风压系数
for i=35:-1:0;
    tic,disp(['当前计算风向角为',num2str(10*i),'°']);
    PP=load ([num2str(10*i),'.txt']);   %加载测点风压时程
    % 测压坏点修正
    % 测点数据整合
    pr=0.5*1.25*8.9583^2;               %求得皮托管(参考高度)的参考风压
    Pmru=PP(:,[258:278,322:342,2:22,66:86,130:150,194:214]);
    Pmrd=PP(:,[279:299,343:363,23:43,87:107,151:171,215:235]);               
    Pmro=Pmru-Pmrd;                     %屋面测点位置的净风压
    % 插值点计算
    Pmro1=Pmro(:,1:21);
    Pmro2=Pmro(:,22:42);
    Pmro3=Pmro(:,43:63);
    Pmro4=Pmro(:,64:84);
    Pmro5=Pmro(:,85:105);
    Pmro6=Pmro(:,106:126);
    Pmr1=[Pmro1, ...
        (Pmro1(:,2)+Pmro1(:,6)+Pmro1(:,7)+Pmro1(:,16))/4, ...
        (Pmro1(:,2)+Pmro1(:,3)+Pmro1(:,8)+Pmro1(:,9))/4, ...
        (Pmro1(:,3)+Pmro1(:,4)+Pmro1(:,10)+Pmro1(:,11))/4, ...
        (Pmro1(:,4)+Pmro1(:,5)+Pmro1(:,12)+Pmro1(:,13))/4, ...
        (Pmro1(:,5)+Pmro1(:,6)+Pmro1(:,14)+Pmro1(:,15))/4, ...
        (Pmro1(:,7)+Pmro1(:,16))/2*1.1, ...
        (Pmro1(:,8)+Pmro1(:,9))/2*1.1, ...
        (Pmro1(:,10)+Pmro1(:,11))/2*1.1, ...
        (Pmro1(:,12)+Pmro1(:,13))/2*1.1, ...
        (Pmro1(:,14)+Pmro1(:,15))/2*1.1];
    Pmr2=[Pmro2, ...
        (Pmro2(:,2)+Pmro2(:,6)+Pmro2(:,7)+Pmro2(:,16))/4, ...
        (Pmro2(:,2)+Pmro2(:,3)+Pmro2(:,8)+Pmro2(:,9))/4, ...
        (Pmro2(:,3)+Pmro2(:,4)+Pmro2(:,10)+Pmro2(:,11))/4, ...
        (Pmro2(:,4)+Pmro2(:,5)+Pmro2(:,12)+Pmro2(:,13))/4, ...
        (Pmro2(:,5)+Pmro2(:,6)+Pmro2(:,14)+Pmro2(:,15))/4, ...
        (Pmro2(:,7)+Pmro2(:,16))/2*1.1, ...
        (Pmro2(:,8)+Pmro2(:,9))/2*1.1, ...
        (Pmro2(:,10)+Pmro2(:,11))/2*1.1, ...
        (Pmro2(:,12)+Pmro2(:,13))/2*1.1, ...
        (Pmro2(:,14)+Pmro2(:,15))/2*1.1];
    Pmr3=[Pmro3, ...
        (Pmro3(:,2)+Pmro3(:,6)+Pmro3(:,7)+Pmro3(:,16))/4, ...
        (Pmro3(:,2)+Pmro3(:,3)+Pmro3(:,8)+Pmro3(:,9))/4, ...
        (Pmro3(:,3)+Pmro3(:,4)+Pmro3(:,10)+Pmro3(:,11))/4, ...
        (Pmro3(:,4)+Pmro3(:,5)+Pmro3(:,12)+Pmro3(:,13))/4, ...
        (Pmro3(:,5)+Pmro3(:,6)+Pmro3(:,14)+Pmro3(:,15))/4, ...
        (Pmro3(:,7)+Pmro3(:,16))/2*1.1, ...
        (Pmro3(:,8)+Pmro3(:,9))/2*1.1, ...
        (Pmro3(:,10)+Pmro3(:,11))/2*1.1, ...
        (Pmro3(:,12)+Pmro3(:,13))/2*1.1, ...
        (Pmro3(:,14)+Pmro3(:,15))/2*1.1];    
    Pmr4=[Pmro4, ...
        (Pmro4(:,2)+Pmro4(:,6)+Pmro4(:,7)+Pmro4(:,16))/4, ...
        (Pmro4(:,2)+Pmro4(:,3)+Pmro4(:,8)+Pmro4(:,9))/4, ...
        (Pmro4(:,3)+Pmro4(:,4)+Pmro4(:,10)+Pmro4(:,11))/4, ...
        (Pmro4(:,4)+Pmro4(:,5)+Pmro4(:,12)+Pmro4(:,13))/4, ...
        (Pmro4(:,5)+Pmro4(:,6)+Pmro4(:,14)+Pmro4(:,15))/4, ...
        (Pmro4(:,7)+Pmro4(:,16))/2*1.1, ...
        (Pmro4(:,8)+Pmro4(:,9))/2*1.1, ...
        (Pmro4(:,10)+Pmro4(:,11))/2*1.1, ...
        (Pmro4(:,12)+Pmro4(:,13))/2*1.1, ...
        (Pmro4(:,14)+Pmro4(:,15))/2*1.1];
    Pmr5=[Pmro5, ...
        (Pmro5(:,2)+Pmro5(:,6)+Pmro5(:,7)+Pmro5(:,16))/4, ...
        (Pmro5(:,2)+Pmro5(:,3)+Pmro5(:,8)+Pmro5(:,9))/4, ...
        (Pmro5(:,3)+Pmro5(:,4)+Pmro5(:,10)+Pmro5(:,11))/4, ...
        (Pmro5(:,4)+Pmro5(:,5)+Pmro5(:,12)+Pmro5(:,13))/4, ...
        (Pmro5(:,5)+Pmro5(:,6)+Pmro5(:,14)+Pmro5(:,15))/4, ...
        (Pmro5(:,7)+Pmro5(:,16))/2*1.1, ...
        (Pmro5(:,8)+Pmro5(:,9))/2*1.1, ...
        (Pmro5(:,10)+Pmro5(:,11))/2*1.1, ...
        (Pmro5(:,12)+Pmro5(:,13))/2*1.1, ...
        (Pmro5(:,14)+Pmro5(:,15))/2*1.1];
    Pmr6=[Pmro6, ...
        (Pmro6(:,2)+Pmro6(:,6)+Pmro6(:,7)+Pmro6(:,16))/4, ...
        (Pmro6(:,2)+Pmro6(:,3)+Pmro6(:,8)+Pmro6(:,9))/4, ...
        (Pmro6(:,3)+Pmro6(:,4)+Pmro6(:,10)+Pmro6(:,11))/4, ...
        (Pmro6(:,4)+Pmro6(:,5)+Pmro6(:,12)+Pmro6(:,13))/4, ...
        (Pmro6(:,5)+Pmro6(:,6)+Pmro6(:,14)+Pmro6(:,15))/4, ...
        (Pmro6(:,7)+Pmro6(:,16))/2*1.1, ...
        (Pmro6(:,8)+Pmro6(:,9))/2*1.1, ...
        (Pmro6(:,10)+Pmro6(:,11))/2*1.1, ...
        (Pmro6(:,12)+Pmro6(:,13))/2*1.1, ...
        (Pmro6(:,14)+Pmro6(:,15))/2*1.1];
    Pmr=[Pmr1,Pmr2,Pmr3,Pmr4,Pmr5,Pmr6];
    % 风压系数计算
    Cpall=Pmr./pr;                      %测点的风压系数时程
    Cpmean=mean(Cpall);                 %测点的平均风压系数
    Cpstd=std(Cpall);                   %测点的脉动风压系数
    % 存储各风向角下风压系数
    CpmeanALL(:,i+1)=Cpmean';
    CpstdALL(:,i+1)=Cpstd';
    disp(['风向角为',num2str(10*i),'°用时:']),toc
end
%% 存储变量到文件
save CpmeanALL.txt CpmeanALL -ascii
save CpstdALL.txt CpstdALL -ascii

标签: none

评论已关闭