[發明專利]一種相鄰樁基水平動力相互作用分析方法及系統

有效

• 申請號: 202011193171.2
• 申請日: 2020-10-30
• 公開（公告）號: CN112227434B
• 公開（公告）日: 2022-04-08
• 發明（設計）人: 崔春義;辛宇;梁志孟;王本龍;孟坤;劉海龍;裴華富
• 申請（專利權）人: 大連海事大學
• 主分類號: G06F30/20
• 分類號: G06F30/20;E02D33/00
• 代理公司: 大連至誠專利代理事務所(特殊普通合伙) 21242
• 代理人: 涂文詩;鄧珂

1.一種相鄰樁基水平動力相互作用分析方法，其特征在于，包括如下步驟：

S1、給定相鄰樁基水平動力相互作用簡化力學模型對應的假定條件，所述假定條件至少包括：樁身被簡化為圓形等截面、均質Euler梁，設定樁周土體為均質土層，并將其簡化為Pasternak地基模型，設定所述模型各部分均滿足小變形條件，樁土界面為完全接觸且無相對滑動；同時設定樁頂處僅發生水平位移，樁底處為固端約束；

S2、創建考慮軸、橫向力作用下，基于Euler梁和Pasternak地基模型的主動樁樁身單元的動力平衡方程；

S3、創建主動樁樁身水平振動位移和轉角關系方程并獲得對應的位移通解；

S4、計算主動樁的樁身水平振動位移、樁身彎矩以及剪力；

S5、計算由主動樁I引起的被動樁II的動態位移并計算出鄰樁相互作用因子；其中，所述S2中的所述模型的創建過程包括：

建立主動樁的單樁水平振動模型樁身單元的動力平衡方程，所述動力平衡方程如下式：

式(1)中：為主動樁I樁身質點的水平位移；E^p、I^p、m^p分別為樁體彈性模量、截面慣性矩和單位長度質量，N₀為作用在樁頂的軸向力；為樁周地基土剪切剛度；B₀為樁的計算寬度，B₀＝0.9(1.5d+0.5)；

其中，和則按如下公式確定：

式中：為樁周土的剪切波速；E^s、ρ^s、β^s和υ^s分別為樁周土的彈性模量、密度、阻尼系數及泊松比；a₀＝ωd/V^s為無量綱頻率；t^s為地基土的剪切層厚度，t^s＝11d，d為樁徑；其中，所述S3中的主動樁樁身水平振動位移和轉角關系方程表示為：

u^p(z，t)＝U^p(z)e^iωt (5)

式中，U^p(z)為樁身質點水平位移幅值；

所述方程的位移通解的獲取過程包括：

將式(5)分別代入式(1)，得到如下方程：

式中，W^p＝E^pI^p，

基于式(6)對應的4個特征根為則其方程位移通解為：

式(7)中，系數A₁₁、B₁₁、C₁₁、D₁₁的取值將由樁頂和樁底邊界條件確定；其中，在所述S4中，基于Euler梁理論，主動樁的樁身轉角、彎矩、剪力與樁身水平位移相互關系為：

其中，式(9)、(10)和(11)中各個系數對應的表達式：

系數A₁₂、A₁₃、A₁₄、B₁₂、B₁₃、B₁₄、C₁₂、C₁₃、C₁₄、D₁₂、D₁₃、D₁₄由樁頂和樁底邊界條件確定；

最后，基于上述主動樁的樁身轉角、彎矩、剪力與樁身水平位移相互關系表達式計算主動樁的樁身水平振動位移、樁身彎矩以及剪力；其中，在所述S4中，系數A₁₁、B₁₁、C₁₁、D₁₁的求解過程包括如下步驟：

基于樁頂和樁底邊界條件，具體邊界條件如下式：

令f₀(z)＝e^λzcosχz，f₁(z)＝e^λzsinχz，f₂(z)＝e^-λzcosχz，f₃(z)＝e^-λzsinχz，f₄＝λ³-3λχ²，f₅＝3λ²χ-χ³，并將上述各個表達式代入式(13)化簡得下式：

式中：[X₁]＝[A₁₁ B₁₁ C₁₁ D₁₁]^T

對式(14)化簡、整理得A₁₁、B₁₁、C₁₁、D₁₁的解析表達式如下：

式中：T＝χ(e^2λl-e^-2λl)+(2sin2χl)

同時，為便于后續分析，引入如下無量綱參數如下：

式中，u_max(z)、m_max(z)、q_max(z)分別為樁基水平振動位移、彎矩和剪力最大值；其中，在所述S5中，計算由主動樁I引起的被動樁II的動態位移的過程包括：

S51、假定主動樁I與被動樁II中各樁的幾何尺寸和材料性質均相同；

S52、基于土體水平位移的衰減函數，獲取由主動樁I引起的場地振動位移，其中，土體水平位移的衰減函數為：

式(17)中：θ為兩樁連線與振動方向x的夾角，S為兩樁的間距；

對應的由主動樁I引起的場地振動位移為：

S53、考慮樁與土體之間的動力相互作用，被動樁II的動力平衡方程為：

由于被動樁II樁身內各單元水平位移和轉角表示為

對式(19)進行化簡進一步得到：

式中：式(20)方程的解由通解和特解兩部分組成，其相應齊次方程的通解為：

式中：λ、χ的表達式A₂₁、B₂₁、C₂₁、D₂₁為待定系數；

式(20)的特解為：

式中：γ₁＝λ+χi；γ₂＝λ-χi

將式(22)代入式(20)分別求得：

則式(20)方程的解為：

S54、確定系數A₂₁、B₂₁、C₂₁、D₂₁，并計算被被動樁II的位移，具體包括：

基于Euler梁理論，被動樁II樁身轉角、彎矩、剪力與樁身水平位移相互關系為：

式中：

式(22)～式(24)中系數由邊界條件確定，對應的表達式A₂₂、B₂₂、C₂₂、D₂₂、A₂₃、B₂₃、C₂₃、D₂₃、A₂₄、B₂₄、C₂₄、D₂₄為：

考慮邊界條件樁頂約束轉角，樁底固定端的情況，則所述邊界條件為：

令F₅＝F₁+F₃；F₆＝F₂+F₄將被動樁II樁身轉角、彎矩、剪力與樁身水平位移相互關系表達式分別代入式(26)得：

式中：[X₂]＝[A₂₁ B₂₁ C₂₁ D₂₁]^T；

求解式(27)得未知變系數表達式A₂₁、B₂₁、C₂₁、D₂₁，并解得被動樁的位移；鄰樁相互作用因子的計算公式為：

2.一種相鄰樁基水平動力相互作用分析系統，其特征在于，包括：

第一數據獲取單元，其用于給定相鄰樁基水平動力相互作用簡化力學模型對應的假定條件，所述假定條件至少包括：樁身被簡化為圓形等截面、均質Euler梁，設定樁周土體為均質土層，并將其簡化為Pasternak地基模型，設定所述模型各部分均滿足小變形條件，樁土界面為完全接觸且無相對滑動；同時設定樁頂處僅發生水平位移，樁底處為固端約束；

第一模型創建單元，其用于創建考慮軸、橫向力作用下，基于Euler梁和Pasternak地基模型的主動樁樁身單元的動力平衡方程；

第二模型創建單元，其用于創建主動樁樁身水平振動位移和轉角關系方程并獲得對應的位移通解；

第一數據獲取單元，其用于計算主動樁的樁身水平振動位移、樁身彎矩以及剪力；

第二數據獲取單元，其用于計算由主動樁I引起的被動樁II的動態位移并計算出鄰樁相互作用因子；其中，所述第一模型創建單元中的所述模型的創建過程包括：

建立主動樁的單樁水平振動模型樁身單元的動力平衡方程，所述動力平衡方程如下式：

式(1)中：為主動樁I樁身質點的水平位移；E^p、I^p、m^p分別為樁體彈性模量、截面慣性矩和單位長度質量，N₀為作用在樁頂的軸向力；為樁周地基土剪切剛度；B₀為樁的計算寬度，B₀＝0.9(1.5d+0.5)；

其中，和則按如下公式確定：

式中：為樁周土的剪切波速；E^s、ρ^s、β^s和υ^s分別為樁周土的彈性模量、密度、阻尼系數及泊松比；a₀＝ωd/V^s為無量綱頻率；t^s為地基土的剪切層厚度，t^s＝11d，d為樁徑；

所述第二模型創建單元中的主動樁樁身水平振動位移和轉角關系方程表示為：

u^p(z，t)＝U^p(z)e^iωt (5)

式中，U^p(z)為樁身質點水平位移幅值；

所述方程的通解的獲取過程包括：

將式(5)分別代入式(1)，得到如下方程：

式中，W^p＝E^pI^p，

基于式(6)對應的4個特征根為則其方程位移通解為：

式(7)中，系數A₁₁、B₁₁、C₁₁、D₁₁的取值將由邊界條件確定；其中，基于Euler梁理論，主動樁的樁身轉角、彎矩、剪力與樁身水平位移相互關系為：

其中，式(9)、(10)和(11)中各個系數對應的表達式：

系數A₁₂、A₁₃、A₁₄、B₁₂、B₁₃、B₁₄、C₁₂、C₁₃、C₁₄、D₁₂、D₁₃、D₁₄由式(13)邊界條件確定；

最后，基于上述主動樁的樁身轉角、彎矩、剪力與樁身水平位移相互關系表達式計算主動樁的樁身水平振動位移、樁身彎矩以及剪力其中，在所述第一數據獲取單元中，系數A₁₁、B₁₁、C₁₁、D₁₁的求解過程包括如下步驟：

基于樁頂和樁底邊界條件，具體邊界條件如下式：

令f₀(z)＝e^λzcosχz，f₁(z)＝e^λzsinχz，f₂(z)＝e^-λzcosχz，f₃(z)＝e^-λzsinχz，f₄＝λ³-3λχ²，f₅＝3λ²χ-χ³并將上述表達式代入式(13)化簡得：

式中：[X₁]＝[A₁₁ B₁₁ C₁₁ D₁₁]^T

對式(14)化簡、整理得A₁₁、B₁₁、C₁₁、D₁₁的解析表達式如下：

式中：T＝χ(e^2λl-e^-2λl)+λ(2sin2χl)

同時，為便于后續分析，引入如下無量綱參數如下：

式中，u_max(z)、m_max(z)、q_max(z)分別為樁基水平振動位移、彎矩和剪力最大值；在所述第一數據獲取單元中，計算由主動樁I引起的被動樁II的動態位移的過程包括：

S51、假定主動樁I與被動樁II中各樁的幾何尺寸和材料性質均相同；

S52、基于土體水平位移的衰減函數，獲取由主動樁I引起的場地振動位移，其中，土體水平位移的衰減函數為：

式(17)中：θ為兩樁連線與振動方向x的夾角，S為兩樁的間距；

對應的由主動樁I引起的場地振動位移為：

S53、考慮樁與土體之間的動力相互作用，被動樁II的動力平衡方程為：

由于被動樁II樁身內各單元水平位移和轉角表示為

對式(19)進行化簡進一步得到：

式中：式(20)方程的解由通解和特解兩部分組成，其相應齊次方程的通解為：

式中：λ、χ的表達式A₂₁、B₂₁、C₂₁、D₂₁為待定系數；

式(20)的特解為：

式中：γ₁＝λ+χi；γ₂＝λ-χi

將式(22)代入式(20)分別求得：

則式(20)方程的解為：

S54、確定系數A₂₁、B₂₁、C₂₁、D₂₁，并計算被被動樁II的位移，具體包括：

基于Euler梁理論，被動樁II樁身轉角、彎矩、剪力與樁身水平位移相互關系為：

式中：

式(22)～式(24)中系數由邊界條件確定，對應的表達式A₂₂、B₂₂、C₂₂、D₂₂、A₂₃、B₂₃、C₂₃、D₂₃、A₂₄、B₂₄、C₂₄、D₂₄為：

考慮邊界條件樁頂約束轉角，樁底固定端的情況，則所述邊界條件為：

令F₅＝F₁+F₃；F₆＝F₂+F₄將被動樁II樁身轉角、彎矩、剪力與樁身水平位移相互關系表達式分別代入式(26)得：

式中：[X₂]＝[A₂₁ B₂₁ C₂₁ D₂₁]^T；

求解式(27)得未知變系數表達式A₂₁、B₂₁、C₂₁、D₂₁，并解得被動樁的位移；

在所述S5中，鄰樁相互作用因子的計算公式為：