一、前言

    粘彈性動力邊界是工程仿真中比較常用，效果也不錯的局部時域人工動力邊界條件，目前已經在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之后自由地表停止振動（圖中藍色部分為數值震蕩），說明入射波在底部黏彈性邊界處被吸收，沒有二次反射。與文獻上的結果也是一樣的。

圖4 頂部位移檢測點時程曲線

圖5 底部監測點時程曲線

文獻結果

   從上述結果可看出，本貼方法對比ABAQUS模擬結果效果相當。當然該方法在數值結果中精確度也有待提高，歡迎給位聯系我互相探討！想要學習本算例的歡迎各位前來咨詢。