黏彈性邊界插件的案例
適用于ABAQUS的黏彈性邊界(粘彈性邊界)及等效地震荷載施加插件程序 v3.2.1 ￥9999
本程序已停售，由于之前有人購買所以無法刪除帖子，價格設置為防拍價，請勿購買，謝謝
粘彈性邊界等效節點力公式的推導(黏彈性邊界)
等效節點力的計算在粘彈性邊界的地震動輸入中至關重要，公式的最終表達式很多論文中都有，但是對于初學者來說，直接使用可能會有些吃力。筆者在前不久發表的論文中對其進行了細致的推導，現在正式版（印刷版）已經刊出，正式版參考文獻鏈接如下，直接點擊文章標題即可：
黏彈性人工邊界在ABAQUS中的實現及地震動輸入方法的比較研究
DOI: 10.13722/j.cnki.jrme.2019.1068
這里將正式版文獻中，正確完整的粘彈性邊界等效節點力公式推導放在下面以供大家參考（公式5-24），希望能及時地給大家帶來一些幫助，相信大家能成功實現粘彈性邊界的地震動輸入。
展開 abaqus黏彈性人工邊界的實現
網上的文章都挺多的，肝了一個月終于肝出matlab程序了，還是挺有成就感的，以后抗震分析可以一鍵生成了哈哈
新論文：黏彈性邊界中靜動力邊界轉化方法（地震靜動力耦合分析）
進行地震靜動力耦合計算時，如果采用黏彈性邊界作為動力邊界條件，則會面臨靜動力邊界轉化的問題。而靜動力邊界條件與地應力平衡橫容易混淆，地應力平衡應該包含在該過程中，許多文獻描述很模糊。新論文：《靜動力邊界轉換及其合理性驗證方法的研究》給出了在進行地震靜動力耦合計算時，靜力邊界條件（固定邊界）向動力邊界條件（黏彈性邊界）轉換的詳細步驟及檢驗方法。論文鏈接：https://doi.org/10.11939/jass.20220136
標簽：粘彈性邊界 黏彈性邊界 等效節點力 靜動力耦合模擬 靜動力邊界條件轉換 黏彈性邊界 疊加原理 地震反應 ABAQUS
展開 
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之后自由地表停止振動（圖中藍色部分為數值震蕩），說明入射波在底部黏彈性邊界處被吸收，沒有二次反射。
展開 