[發(fā)明專利]一種混凝土重力壩材料參數(shù)空間變異性的模擬方法

有效

• 申請?zhí)枺?/td>: 201810410957.1
• 申請日: 2018-05-02
• 公開（公告）號: CN108629111B
• 公開（公告）日: 2019-09-10
• 發(fā)明（設(shè)計）人: 吳震宇;盧祥;陳建康;李艷玲;裴亮;張瀚
• 申請（專利權(quán)）人: 四川大學(xué)
• 主分類號: G06F17/50
• 分類號: G06F17/50
• 代理公司: 北京眾合誠成知識產(chǎn)權(quán)代理有限公司 11246
• 代理人: 夏艷

1.一種混凝土重力壩材料參數(shù)空間變異性的模擬方法，其特征在于，包括以下步驟：

(1)建立重力壩模型，對所述模型進(jìn)行網(wǎng)格劃分，獲取每個網(wǎng)格單元的編號和形心坐標(biāo)，組成矩陣其中n為網(wǎng)格單元的個數(shù)，x_i,y_i分別為第i個網(wǎng)格單元形心的橫、縱坐標(biāo)，1≤i≤n；

(2)利用MATLAB軟件生成隨機矩陣A_m×n，其中m為每個網(wǎng)格單元參數(shù)的隨機抽樣次數(shù)，n為網(wǎng)格單元的個數(shù)；

(3)利用自相關(guān)函數(shù)計算所述矩陣E_n×3的相關(guān)系數(shù)矩陣其中元素ρ_i,j為第i和j個網(wǎng)格單元中參數(shù)的相關(guān)系數(shù)，1≤i≤n，1≤j≤n；

(4)對所述矩陣ρ_n×n進(jìn)行Cholesky分解，即ρ_n×n＝B_n×n^T·B_n×n，得到上三角矩陣B_n×n；

(5)利用所述矩陣B_n×n對所述矩陣A_m×n進(jìn)行線性變換，即D_m×n＝A_m×n·B_n×n，得到矩陣D_m×n；

(6)相關(guān)隨機場X_m×n的計算：通過公式lnX_m×n＝σ_lnx×D_m×n+μ_lnx計算得到lnX_m×n，相關(guān)隨機場其中μ_lnx和σ_lnx為滿足對數(shù)正態(tài)分布的參數(shù)。

2.根據(jù)權(quán)利要求1所述的一種混凝土重力壩材料參數(shù)空間變異性的模擬方法，其特征在于，步驟(1)中所述劃分網(wǎng)格的尺度為1～1.5m。

3.根據(jù)權(quán)利要求1所述的一種混凝土重力壩材料參數(shù)空間變異性的模擬方法，其特征在于，步驟(3)中所述自相關(guān)函數(shù)為負(fù)指數(shù)自相關(guān)函數(shù)，所述元素其中θ_x和θ_y分別為橫、縱坐標(biāo)方向上混凝土材料的相關(guān)距離，x_i,y_i分別為第i個網(wǎng)格單元形心的橫、縱坐標(biāo)，x_j,y_j分別為第j個網(wǎng)格單元形心的橫、縱坐標(biāo)，1≤i≤n，1≤j≤n。

4.根據(jù)權(quán)利要求1所述的一種混凝土重力壩材料參數(shù)空間變異性的模擬方法，其特征在于，步驟(6)中所述滿足對數(shù)正態(tài)分布的參數(shù)μ_lnx和σ_lnx通過下式計算，其中μ_x為混凝土材料參數(shù)的均值，α_x為混凝土材料參數(shù)的變異系數(shù)。