[发明专利]一种多维多变量非高斯空间随机场模拟方法

说明书

技术领域

本发明涉及土木工程可靠度分析方法，尤其涉及一种多维多变量空间随机场模拟方法。

背景技术

材料参数的空间随机场模拟是土木工程可靠度研究与设计领域的重要工作之一，由于土木工程材料在制造过程中制造工艺的缺陷常导致材料的物理力学性质在材料空间内存在一定的空间变异性，而非一个常数；同样，天然土体的工程性质由于受沉积、后沉积、化学风化和搬运等作用以及不同荷载历史的影响而呈现一定的空间变异性。因此，在对土木工程结构物进行安全评估时应当考虑材料、岩土体物理力学参数的不确定性，应当采用基于可靠度理论的土木工程设计理念，而提供材料参数的空间随机场是可靠度设计与研究的基础工作。

目前多维多变量非高斯空间随机场模拟方法方面的研究，主要是先通过给定的表征参数空间变化特性的相关距离，使用已有相干函数模型，建立参数的多维功率谱模型，再运用谱表示法生成多维多变量高斯空间随机场，然后根据材料参数的非高斯分布特性，运用非线性转换模型，得到符合给定分布特性的材料参数的多维多变量非高斯空间场。其核心步骤为：通过实验数据，求得参数的相关距离；建立各参数之间的目标多维非高斯功率谱矩阵；使用某个多维功率谱矩阵，采用谱表示法生成潜在的多维多变量高斯空间随机场；通过非线性转换，将潜在的多维多变量非高斯空间随机场转换为多维多变量非高斯空间随机场。

其中，最为困难的是“某个多维功率谱矩阵”的确定。传统的方法主要是直接使用目标多维非高斯功率谱矩阵，由于从高斯空间场转换到非高斯空间场采用了非线性转换，这样导致最终得到的非高斯空间随机的功率谱矩阵与目标功率谱矩阵不一，使得模拟结果存在误差。

发明内容

发明目的：本发明的目的是提供一种采用迭代程序求解潜在多维高斯功率谱矩阵，从而使得所模拟的多维多变量非高斯空间随机场的功率谱与目标功率谱一致的多维多变量空间随机场模拟方法。

技术方案：一种多维多变量非高斯空间随机场模拟方法，包括如下步骤：

(A)通过参数的相关距离，建立参数之间的目标多维非高斯功率谱矩阵，给定一个目标多维非高斯功率谱矩阵与参数的分布特性；

(B)假定一个初始的多维高斯功率谱矩阵为目标多维非高斯功率谱矩阵；

(C)由上述多维高斯功率谱矩阵求解多维非高斯功率谱矩阵；

(D)比较多维非高斯功率谱矩阵与目标多维非高斯功率谱矩阵之间的误差；

(E)对目标多维非高斯功率谱矩阵、多维高斯功率谱矩阵以及多维非高斯功率谱矩阵进行分解；

(F)对分解后矩阵的元素采用指数型迭代公式进行迭代；

(G)通过对迭代后得到的元素求解新的多维高斯功率谱矩阵，并对矩阵元素进行标准化；

(H)对各个变量进行随机排序；

重复步骤(C)至步骤(H)，直到误差矩阵的平均误差大于前一次迭代的平均误差，得到潜在多维高斯功率谱矩阵；

(I)使用迭代得到的潜在多维高斯功率谱矩阵模拟多维多变量高斯空间随机场，再通过非线性变换，得到所要模拟的多维多变量非高斯空间随机场。

步骤(C)中，采用多维傅里叶变化公式，由所述多维高斯功率谱矩阵求解得到高斯互相关函数再根据下述公式求解得到非高斯互相关函数：

其中，φ[g_j1,g_k2；ρ_jk(ξ)]是第j个元素与第k个元素的在相关系数ρ_jk(ξ)下的联合正态分布函数，表示第j个函数概率分布函数的逆函数；

由所述非高斯互相关函数求解得到多维非高斯功率谱矩阵。

步骤(I)具体为：首先使用谱表示法将迭代得到的潜在多维高斯功率谱矩阵模拟多维多变量高斯空间随机场，其次通过基于高斯分布与目标分布之间的非线性转换关系，得到所要模拟的多维多变量非高斯空间随机场。

步骤(F)中，采用了以下指数型迭代公式：

其中，和为目标多维非高斯功率谱矩阵多维高斯功率谱矩阵以及多维非高斯功率谱矩阵分解后矩阵的元素，参数β用来控制收敛速率；参数β选取为0.5-2，优选为0.8-1.2。

工作原理：本发明提出了一种确定潜在高斯空间场模拟所用的多维高斯功率谱矩阵的迭代求解方法，得到能够使得最终的多维非高斯功率谱矩阵与目标多维非高斯功率谱矩阵一致的多维高斯功率谱矩阵；通过迭代求解得到的潜在多维高斯功率谱矩阵后，使用谱表示法模拟多维多变量高斯空间随机场，然后再通过目标概率分布函数与高斯分布之间的非线性转换，将其转化为多维变量非高斯空间随机场。