1.一种考虑风压关联性和风速方向性金属屋面风灾估计方法，其特征在于，考虑风荷载相关性和风速方向性，包括以下步骤：步骤1：采用POD插值得到屋面的外部风压设定测点处风压通过对已有测点处风压进行POD插值而得出，通过POD插值对低矮房屋金属板屋面进行风压重建；

为一组零均值的N元脉动风压系数向量，其中N表示风洞试验中测点的数量，利用POD找到一组最优正交基Θ＝[Θ1,Θ2,...,ΘN]，则CP(t)展开为其中，ai(t)是CP(t)在基向量Θi上的投影，i＝1,2,…,N；

基向量组Θ通过下列特征值方程确定

RpAΘ＝ΘΛ    (2)

其中，Rp是CP(t)的协方差矩阵；A＝diag(A1,A2,…,AN)，Ai是第i个测点处风压贡献面积；对角矩阵Λ＝diag(λ1,λ2,…,λN)为RpA的特征值矩阵，λ1、λ2、…、λN分别为RpA的N个特征值；当测点均匀分布且各测点风压贡献面积相等时，公式(2)写为RpΘ＝ΘΛ#    (3)

其中， 式中A表示测点处风压贡献面积；当测点分布不均匀时，将公式(2)乘以A1/2，则公式等效为其中， Θ*＝A1/2Θ，此时公式中的转换矩阵 是一个实对称矩阵，正交基向量组通过求逆得到

Θ＝A-1/2Θ*    (5)

通过对解得的基向量空间插值，得到设定测点处新向量分量，乘以投影a(t)得到这些设定测点处的脉动风压，进而得到屋面的外部风压；

步骤2：通过模拟得到内压

对一个单一开口，其控制方程为：

其中ρ是空气密度；le是“气塞”的有效长度；V是内部体积；γ是空气的比热容；a是开口面积；P0是大气压力； 是参考高度处的平均风速；κ是流量系数，弥补其他的能量损失；Cpe和Cpi分别表示外压和内压系数；

步骤3：螺钉上内力及其极值内力的确定通过影响系数来确定螺钉内力，获得螺钉内力的公式为：X(t)＝∫∫q(x,y,t)Ic(x,y)dxdy    (7)其中，Ic(x,y)是在(x,y)处螺钉内力的影响系数；q(x,y,t)表示相应的总压力，即外压和内压的总和；用 表示T1分钟过程内的极值内力，设它的累积分布函数服从Gumbel分布，即其中，w表示概率分布函数中的自变量；模态 是最可能的取值， 是散度；T分钟内极值W的累积分布函数为计算得到T1分钟内的累积分布函数为：其中，新的模态和散度为：

它们的乘积 是一个无量纲参数；

步骤4：求解螺钉内力间的相关性

采用Nataf转换来关联非高斯变量和高斯变量，用来求解极值内力的相关性，设W＝[W1,W2,…,Wn]T表示螺钉上极值内力分量的一个随机变量向量组，它的累积分布函数是n是螺钉的数量，此向量通过T

转换为相应的标准高斯随机向量Z＝[Z1,Z2,…,Zn] ，其中Φ(·)表示标准高斯累积分布函数；通过得到相应的Nataf转换；其中，fW(w)是极值内力的联合概率密度函数， 是相关矩阵RZ的n元标准高斯联合概率密度函数， 的表示Zj的高斯概率密度函数，j＝1,2,…,n；

是Wj和Wk间的相关系数，表示为

其中，μj和σj分别是Wj的均值和标准差，μk和σk分别是Wk的均值和标准差；

表示Zj和Zk的联合高斯概率密度函数；Zj和Zk之间的相关系数 是RZ中一个元素；若Wj和Wk均服从Gumbel分布，公式(14)近似表示为步骤5：金属屋面板失效概率和屋面损失率当使用高强度G550型号钢时，强度R的均值为其中，钢材的屈服应力为fy；螺钉帽的直径为dh；杨氏模量为E；金属板波峰高度为hc；波峰顶间距为hp；波谷宽度为Wt；金属板波峰宽度为Wc；金属板厚度为t；檩条之间的距离为L；

变异系数是0.12；

螺钉周围的金属板强度通过以下公式得到其中，fR(r)是螺钉周围金属板强度的概率密度函数；

金属板的失效概率由金属板上的极值内力和螺钉处的承载力决定，失效概率表示为其中， 是金属板强度的联合高斯概率密度函数；

设整个屋面金属板的数量是NC，引入损失率来描述屋面板的损坏程度，定义为失效金属板的百分比，即D＝MC/NC    (19)

其中，MC为失效金属板的数量；

运用蒙特卡洛模拟来估计金属板的失效概率和整个屋面的损失率；假设蒙特卡洛模拟重复了nt轮，在第m轮模拟中，首先通过Nataf转换来模拟所有相关螺钉的极值内力，相关高斯向量Z根据下述公式进行模拟：RZ＝LLT；Z＝L-1U    (20)其中，U是独立标准高斯向量，通过对RZ进行Cholesky分解得到下三角矩阵L；

对矩阵RZ进行重写：

RZ＝VTΩV    (21)

其中，V是特征向量矩阵，Ω是对角特征值矩阵；

用fl,m表示第l块金属板在第m次模拟是否破损坏，fl,m＝0或1分别表示未损坏或损坏，则第l块金属板的失效概率为：其中， 是第l块金属板在模拟中失效次数，第m次模拟的损失率为：其中， 是第m次模拟失效的金属板块数；由于随机变量D为近似高斯分布，损失率的均值和标准差为：步骤6：考虑风速方向性的金属屋面板的失效概率第一种情况：不考虑屋面损失率的变异性表示各向年最大风速联合分布，其中nd表示方向分段的总数，由多元极值理论得到联合累积分布函数，使用高斯Copula模型表示为其中Φ-1是标准高斯分布的累积分布函数的反函数； 是nd维零均值和协方差矩阵为Σ的正态分布，其中Σii＝1，是i方向上年最大风速的累积分布函数；

考虑各向风速屋面不超过损失等级d的概率按下式计算当各向风速相互独立时，则公式(26)简化为重现期为Y年的屋面损失率低于等级dY通过下式确定Y＝1/[1-P(D≤dY)]    (28)第二种情况：考虑屋面损失率的变异性

各向风速内屋面损失率D未超过损失等级d的概率表示为其中， 是在风速vi，i＝1,2,…,nd的条件下损失率的联合累积分布函数，屋面破坏在各方向认为是彼此独立的，则公式(29)写为其中， 是风速vi在第i个方向损失率的累积分布函数， 是的概率密度形式；若各向极值风速是相互独立的，则近似为其中， 是第i个方向上的极值风速的概率密度函数。