基于比例邊界有限元法的結(jié)構(gòu)地基動力相互作用時(shí)域算法的改進(jìn)

1、 基金項(xiàng)目 作者簡介 河南新密人 博士生 主要從事結(jié)構(gòu) 基于比例邊界有限元法的結(jié)構(gòu)土木水利學(xué)院 遼寧 大連摘要 通過選擇適當(dāng)?shù)幕瘮?shù)數(shù)目 縮減結(jié)構(gòu) 以減小卷積積 分所帶來的計(jì)算工作量結(jié)果表明 基函數(shù)的縮減可使計(jì)算效率明顯 提高 但精度損失不大 當(dāng)采用 計(jì)算效率提高因此 該算法比較適合于大型結(jié)構(gòu)關(guān)鍵詞 結(jié)構(gòu) 比例邊界有限元法 縮減基函數(shù)法 時(shí)域算法 結(jié)構(gòu) 尤其對于大型結(jié) 但是有限元 條件 其次 它只需離散部分邊界此外 近幾年這種方法已經(jīng)應(yīng)用到了結(jié)構(gòu) 不僅能夠用于大型工程問題的數(shù)值計(jì)算時(shí)仍然存在比較大的困難特別是對于大型結(jié)構(gòu) 域算法做了改進(jìn) 式中 表示維數(shù) 二維和三維時(shí)分別等于 和 其中系數(shù)矩陣表

2、示為 為邊界 上的 對于結(jié)構(gòu) 可推導(dǎo)出無限地基位移動力剛度的 則含有 個(gè)未知數(shù)的一階常微分方 程 其中系數(shù)矩陣為 和 則由于 并且存在關(guān)系 與位移動力剛度之間關(guān)系為 時(shí)域中加速度 得到 結(jié)構(gòu) 時(shí)域中結(jié)構(gòu) 其中子塊的上標(biāo) 表示結(jié) 和 分別表示結(jié)構(gòu)的自由度和結(jié)構(gòu)與地基接觸面上的自由度 式中 所以 時(shí)刻 和 對于 時(shí)刻 算例驗(yàn)證 二維重力壩 以 高 的 重力壩為例 斷面如圖 壩體 材料彈性模量為 從接觸面上均勻 輸入水平向和豎向的單位加速度脈沖 變換后得到壩頂加 本文分別選取 圖 當(dāng)自由度縮減至總自由度的約 水平向和豎向所產(chǎn)生的最大誤差分別 為 圖 重力壩斷面 水平向和豎向壩頂加速度頻響曲線三維拱壩

3、 影響 可虛擬一個(gè)結(jié)構(gòu) 以 高的大崗山雙曲 地基材料彈性模量為 壩體及近場地基離散如圖圖 壩體及近場地基的離散觸面上有 個(gè)節(jié)點(diǎn) 各方向上從接觸面上均勻輸入單位加速度脈沖 計(jì)算了拱冠梁頂部的加速度頻 個(gè)基函數(shù) 與不進(jìn)行自由度縮減計(jì)算的結(jié)果進(jìn)行對比 如圖 所 響曲線 本文分別選取 示 當(dāng)自由度縮減至總自由度的約 ? 時(shí) 個(gè)基函數(shù) 橫河向 順河向和豎向所產(chǎn)生的最大誤差分 別為 和 圖 各個(gè)方向的拱冠梁處加速度頻響曲線 分析與結(jié)論 在原來的時(shí)域算法中 卷積積分的運(yùn)算在整個(gè)分析中所占的計(jì)算量是比較大的 并且隨著所分析時(shí) 程總時(shí)間的增長而增大 在改進(jìn)的時(shí)域算法中 卷積積分并沒有避免 分別出現(xiàn)在由式 計(jì)算加

4、速度 單位脈沖響應(yīng)矩陣時(shí)和式 計(jì)算相互作用力時(shí) 珟 珟 珟 個(gè)自由 假設(shè)所分析時(shí)程 被離散為 個(gè)時(shí)間步 接觸面上有 度 則卷積積分中矩陣為 向量為 如果記一次乘法和加 法為一個(gè)計(jì)算量單位的話 則算法中僅卷積積分所引起的計(jì)算量為 ? ? 當(dāng)自由度縮減為 ? 時(shí) 通過計(jì)算可知計(jì)算量為原來計(jì)算量的 倍 從算例的結(jié)果可以 看出 此時(shí)計(jì)算精度已經(jīng)達(dá)到很高 圖 給出了與未縮減的算法相 圖 自由度縮減倍數(shù)與卷積積分 計(jì)算效率的關(guān)系 比的自由度縮減倍數(shù)與卷積積分計(jì)算量的關(guān)系 由于在整個(gè)分析過程中 卷積積分運(yùn)算所占整個(gè)計(jì)算量的比重 很大 所以以上只是討論了卷積積分部分所節(jié)約的計(jì)算量 同樣 對 于此算法中其它計(jì)算

5、 如特征值求解 求解方程組等 效率同樣能夠得到不同程度的提高 例如對于本文 三維拱壩算例 改進(jìn)后的算法總耗時(shí)僅約為原算法總耗時(shí)的 從算例的結(jié)果中看出 為了保證計(jì)算結(jié)果的可靠性 采用的基函數(shù)數(shù)目不能太少 一般要達(dá)到接觸 面總自由度的 左右 同時(shí) 本文提出的改進(jìn)算法與原有算法相比 計(jì)算效率得到很大程度的提高 博士提供相關(guān)的資料 并提出寶貴意見 并且保證了計(jì)算精度 比較適合于大型結(jié)構(gòu) 地基動力相互作用的時(shí)域分析 致謝 感謝澳大利亞南新威爾士大學(xué) 參 考 文 獻(xiàn) 劉云賀 張伯艷 陳厚群 拱壩地震輸入模型中黏彈性邊界與黏性邊界的比較 水利學(xué)報(bào) 吳健 金峰 張楚漢 等 無限地基輻射阻尼對溪洛渡拱壩地震響應(yīng)的影響 巖土工程學(xué)報(bào) 杜建國 林皋 地基剛度和不均勻性對混凝土大壩地震響應(yīng)的影響 巖