粘彈性邊界及地震動輸入的案例
COMSOL黏彈性動力邊界及地震動輸入
一、前言
粘彈性動力邊界是工程仿真中比較常用，效果也不錯的局部時域人工動力邊界條件，目前已經在ANASYS、ABAQUS和Fssicas等通用有限元軟件中有了較為通用的使用方法，但是在COMSOL這款以多物理場和PDE建模為特色的通用軟件中卻比較少見。因此本帖展示的是本人在COMSOL有限元平臺實現的粘彈性邊界的施加以及地震動輸入的介紹。
本貼采用的驗證算例引用于文獻《黏彈性人工邊界在ABAQUS中的實現及地震動輸入方法的比較研究》-巖土力學與工程學報-馬笙杰等。
下面是建模介紹和模擬結果與文獻結果的對比驗證。
二、模型建立
通過場外垂直入射sv波算例來驗證黏彈性邊界設置和地震動輸入的準確性。在二維無限彈性空間中截取長50m，高50m的有限元區域作為計算區域，設置模型的頂部中點和底部中點作為監測點，如圖1所示，模型材料參數如下：密度為2000kg/m^3,彈性模量：2e8[Pa],泊松比0.25，剪切波速為200m/s，采用四邊形網格單元，網格尺寸為0.5m×0.5m，在模型底部垂直輸入sv波，波形和速度圖像如圖2、3所示。持續時間為0.2s，計算時長為1s，計算時間步為0.001s，瞬態隱式求解，時間進步方法為向后差分。
圖1 二維土體計算模型
圖2 入射波位移時程曲線圖
圖3 入射波速度時程曲線圖
計算結果如圖4、5所示，入射波在經過0.25s之后到達自由表面與反射波疊加，變成入射波位移的2倍，0.4s之后自由地表停止振動（圖中藍色部分為數值震蕩），說明入射波在底部黏彈性邊界處被吸收，沒有二次反射。
展開 粘彈性人工邊界及地震動輸入方法
最近看到大家有不少人想用abaqus做粘彈性人工邊界，一般情況下粘彈性人工邊界搭配等效節點力地震動輸入方法效果更好，所以在半無限域地基和地下結構的模擬中，經常采用粘彈性邊界和等效節點力輸入方法。
在這里我給大家推薦網上一位大神做的算例，非常不錯，應該可以解決大家的問題。因為大神非常慷慨地在網上公布了算例和代碼，我姑且認為大神是希望與大家分享，因而我本著讓大家早日解決問題的心態把他的算例網址公布在這里，如有侵權我會立即刪除。
網址：https://www.cnblogs.com/w-tao13614/p/10498697.html
下面是大神做出的模擬效果：
最后感謝大神的無私奉獻
展開 ABAQUS中粘彈性邊界的實現（二維+三維，均質+多層介質）
本教程深入講解了粘彈性邊界理論及地震動轉換為等效節點力的理論基礎，并通過實際編程演示，詳細展示了如何在ABAQUS軟件中實現粘彈性邊界和節點力地震動的輸入。
針對均質土體，教程介紹了使用MATLAB軟件計算彈簧、阻尼文件及等效節點力文件的全過程，并在ABAQUS中構建模型。通過添加關鍵字的方式將這些文件整合至模型中，再次導入ABAQUS進行地震響應計算。針對多層土體，教程基于波動理論和斯奈爾定律，推導出粘彈性人工邊界的分層土地震等效節點力計算公式，并編寫了相應的MATLAB程序。在完成彈簧、阻尼文件及等效節點力文件的計算后，這些文件將被導入ABAQUS進行驗證。此外，教程還詳細探討了地下工程中土體初始地應力的影響，即粘彈性邊界中靜—動力邊界的轉化問題。
本教程適用于均質介質和多層介質中的橫波與縱波計算，涵蓋二維和三維模型。如有需求或疑問，請聯系 QQ: 2636336968。
展開 錯過等一年！技術鄰雙十一精選課程案例6折起！更有VIP等你來拿！
www.yqgqt.org.cn/video/c15288
六折
ABAQUS建立CRTSⅠ型無砟軌道(梁-殼)模型【以漁計劃】
https://www.yqgqt.org.cn/video/c15291
六折
ABAQUS-Simpack車輛-柔性軌-浮置板耦合動力學20講【完結】
https://www.yqgqt.org.cn/video/c15344
六折
Abaqus-CEL方法在爆炸中的應用
https://www.yqgqt.org.cn/video/c15594
六折
UMESHMOTION子程序
https://www.yqgqt.org.cn/video/c15742
六折
CAE車輛-軌道耦合動力學仿真
https://www.yqgqt.org.cn/video/c15775
六折
abaqus切削仿真（切屑調試不求人）
https://www.yqgqt.org.cn/video/c15950
六折
ABAQUS Python二次開發工程實例教程
https://www.yqgqt.org.cn/video/c16214
六折
ABAQUS模擬混凝土水化熱溫度場、熱應力裂縫擴展（XFEM）
https://www.yqgqt.org.cn/video/c16449
六折
粘彈性邊界及節點力地震動輸入在
展開 
