abaqus 地震輸入的案例
ansys之——地震波的輸入和求解
對于地震波的輸入，可以把荷載記錄做成文件，利用apdl的讀取功能讀入倒數據庫中。下面的例子是自己編的一個小文件。修改一下可以更簡潔。有用到的朋友自己作一下把。
fini
/config,nres,1000
*dim,aceX,TABLE,3000,1
*dim,aceY,TABLE,3000,1
*dim,aceZ,TABLE,3000,1
*creat,ff
*vread,aceX(1,1),acex,txt,,1
(e16.6)
*vread,aceX(1,0),ACETT,,,1
(e17.6)
ACEX(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceY(1,1),acey,txt,,1
(e16.6)
*vread,aceY(1,0),ACETT,,,1
(e17.6)
ACEY(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceZ(1,1),acez,txt,,1
(e16.6)
*vread,aceZ(1,0),ACETT,,,1
(e17.6)
ACEZ(0,1)=1
*end
/input,ff
！地震波時程記錄分成了3個文件，每個文件是一列。分別記錄x,y,z方向的加速度。acett是時間記錄。
這樣就可以把加速度記錄讀取倒ansys數據庫中作為數組。
也可以把加速度記錄做成一個文件，這樣程序就簡單多了。大家可以試看看修改一下。
下面是計算部分語句：
/SOLU
ANTYPE,trans
！
展開 多層復雜土體地形上部既有斜拉橋地震波輸入模型
圖片內容為多層復雜土體地形上部既有斜拉橋地震波輸入模型以及附有對應腳本程序和使用教程。 提供了一套強大而高效的MATLAB自動化腳本，它能一鍵生成精確的粘彈性人工邊界和復雜的等效節點荷載，徹底解決因邊界處理不當導致的波形反射與結果失真問題。通過我這套教程，能讓學習者學會如何將其應
abaqus 自重+地震，其中地震用反應譜分析，如何做？
abaqus 自重+地震，其中地震用反應譜分析，如何做？
ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬 ￥70
ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬
一、建模技術
地震工況下邊坡可能失穩進而出現滑坡現象，為避免模擬滑坡時網格產生的畸變問題，采用耦合歐拉拉格朗日法(CEL)進行滑坡的大變形模擬；土體本構采用摩爾庫倫模型；采用模型底部小范圍內的周期性荷載模擬地震荷載。
二、模型及部分結果展示
圖1：藍色為邊坡；紅色為空氣層
圖2：網格的劃分
圖3：賦予模型初始應力
圖4：土體達到地應力平衡時的應力分布
圖5：土體底部的地震荷載施加區域
圖6：所施加的周期性荷載（地震荷載）
圖7：邊坡因地震荷載產生的位移
圖8：地震波產生的區域
展開 
abaqus地震方面的資料以及總結 ￥2
abaqus地震方面的資料以及總結
ABAQUS UEL 二次開發（Koyna混凝土壩地震動力響應分析）
第1階模態云圖（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
第2階模態云圖（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
第3階模態云圖（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
第4階模態云圖（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
（三）Koyna混凝土壩地震動力響應分析
在壩體底部邊界輸入加速度時程（地表水平向、豎直向加速度地震動記錄見附件）。
（1）時程數據對比
壩頂-壩踵水平向相對位移時程
壩頂-壩踵豎直向相對位移時程
壩頂-壩踵水平向相對速度時程
壩頂-壩踵豎直向相對速度時程
（2）云圖對比
第4s水平向位移云圖對比（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
第4s水平向位移云圖對比（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
（3）峰值數據對比
統計三種情況的壩頂-壩踵相對數據峰值，比較表格如下。
(四)附件
附件包括包含兩個文件夾，分別為abaqus自帶單元計算文件和自編uel計算文件（for文件加密后的obj文件）。不包含sbfem的相關計算文件。
Koyna dam dynamic analysis.zip
?注：帖子不定時更新，也可能永遠不更新，慎重參考，如給您帶來誤導，深感抱歉。
展開 ABAQUS中粘彈性人工邊界及地震力的施加
實現方式：模型底部邊界節點施加地震力 （集中力）
1、線性攝動分析提取模型固有頻率f
設置材料參數和邊界條件，不需輸入荷載，分析步時間任意。
編輯step，設置提取最大階數位3：
計算結束提取模型一階固有頻率（Mode-1，Freq=2.3946）
局部放大：
2、創建Mesh-Part和Instance
在Part模塊建立part之后，不對其賦予屬性，直接到Mesh模塊進行網格劃分，之后對劃分好網格的part創建mesh-part【菜單欄-Mesh-Create Mesh Part】，之后再回到Property模塊對mesh-part賦予屬性，再到Assembly模塊用mesh-part創建Instance。
3、地應力平衡
運用生死單元將墻體殺死，對填土和地基進行平衡。約束所有土體左右側水平位移（尤其是填土左側，第二部地震時需關掉此邊界條件）。平衡后的模型位移可達e-17次方數量級。
4、創建地震分析部并編輯
分析步類型為動態隱式（Dynamic，Implicit）【動態顯示Dynamic/Explicit不支持彈簧剛度的定義】，開啟大變形。
分析步時間為地震響應時間，最大增量調到很大以保障足夠多的increment，同時初始盡可能地小以保證迭代收斂，最小也更小；最大increment不能選為默認，最大設為0.1（最好是地震時間間隔0.02s，但計算代價太大）。
5、激活上一步殺死的墻體，并設置墻體和填土與地基間的接觸屬性，其中墻-地基不可分離（Rough），墻-填土可分離（Penalty）。在第二步創建接觸屬性：
6、創建底部邊界阻尼（Interaction模塊）
共分三個步驟：創建幾何、創建屬性、賦予屬性。
展開 abaqus用戶材料參數輸入
abaqus安裝目錄下有一個文件夾，其中有一行代碼是控制用戶材料模塊中參數輸入的問題。
有誰知道這行代碼在哪里嗎？
應用ABAQUS進行復雜建筑結構的彈塑性地震反應分析
Abaqus上海土木研討會上的演講ppt.
Abaqus.rar
abaqus地震分析關于模型漂移問題
關于解決地震分析模型漂移的方法，某前輩在博客已經做了比較完美的解答，即基線修正法，http://blog.sina.com.cn/s/blog_417a7d700100z77g.html但這個方法只適用于隱式分析。如果模型是顯式分析該怎么解決呢？在顯式分析中如果添加了基線修正，就會提示keyword在顯式分析中不可獲得。自己的模型屬于車橋耦合振動，比較復雜，隱式計算太耗費時間，所以采用了顯式計算。有沒有在顯式分析的前提下解決模型漂移的辦法呢？求助.........
也歡迎大家在此帖下討論
基于python編程操作ABAQUS輸入文件生成PD3D單元顆粒
在空間中生成剛性顆粒（注意是剛性顆粒）有下列幾種方法：
1.修改關鍵字，構建粒子生成器模型生成隨機分布剛性顆粒
2.使用python語言直接在ABAQUS中生成顆粒，并進行剛體綁定，使其成為剛性顆粒，或者直接生成解析剛體或離散剛體。
方法1生成顆粒的隨機性較好，操作簡單。方法2直接在ABAQUS界面生成顆粒，當所需顆粒數量以萬為計量單位時，在前處理界面時就會卡死，對顯卡要求極高。因此，在僅考慮到這些弊端情況下，就已經使研究人員頭皮發麻，無從下手。
在一些特定應用場合下，比如所需顆粒數量數以萬計，我們只能采用方法1生成顆粒，但我們不僅僅是需要顆粒，還需將這些顆粒與其它模型進行耦合求解計算，這個時候粒子生成器就會有局限性。此外，考慮到顆粒在空間中排布的多樣性，比如最典型的高斯分布，那么粒子生成器很難做到一步到位生成所需分布特征的顆粒。
我們今天介紹的通過python編程操作ABAQUS輸入文件生成PD3D單元顆粒，其可操作性更強，我們可以不采用粒子生成器內部定義的隨機算法生成顆粒，用戶可以根據需求自定義顆粒分布算法，以契合實際工況。此外，可省去粒子生成顆粒的分析步，直接進行工況建模求解計算。
本貼只是個人興趣，只提供思路，不提供源碼，用戶需了解ABAQUS的inp文件的書寫規則、python操作文件語法和生成顆粒的底層邏輯(分布模型)。感興趣的可以私信，提供編寫思路。
下面我們采用這一方法生成直徑2mm、3mm、4mm和5mm的混合顆粒，數量為1000。具體生成結果如下圖所示。
展開 
Abaqus模擬 | 響應譜法求解結構地震響應
反應譜按7度多遇地震，取地震影響系數為0.08，第一組，III類場地卓越周期Tg=0.45s，求解結構的最大位移和內力。
步驟1：根據條件計算對應的設計反應譜
在這里，我們采用網上現有的小程序進行生成，創建QM-KZZY文件，生成反應譜的相關數據，得到了周期—加速度反應譜關系中的600個數據點。
我們選取重要的14個數據，滿足曲線的完整性，同時對前兩列數據進行處理、單位換算得到Abaqus中需要輸入的頻率和加速度，其中頻率為周期的倒數，加速度為Sa(g)*9.86。
步驟2：創建模型，材料屬性和截面及集中質量
本次統一量綱位m（米）級別，我們采用線單元建立懸壁柱；建立鋼材，輸入密度7800，楊氏模量2.1e11，泊松比0.3；創建工字鋼梁截面，賦予給懸壁柱，注意要指派截面方向。
關于集中質量的指定，我們在特殊設置—慣性—管理器中創建類型為點質量的慣性，即3m處為160kg，其他均為120kg。
步驟3：設置響應譜分析步
反應譜分析要創建兩個分析步，第一個分析步是模態分析，是響應譜分析的基礎，選擇線性攝動—頻率，設置特征值個數為10，即只關心前10階模態。
展開 ABAQUS UEL 二次開發（Koyna混凝土壩地震響應計算）
概述：以Koyna混凝土壩為對象進行地震響應計算。將自編的八節點UEL和二十節點UEL應用到計算中。分別進行了混凝土壩模態計算和地震時程計算。
其中，在模態計算中共設置四種計算工況，分別為：ABAQUS-C3D8、UEL-C3D8、ABAQUS-C3D20、UEL-C3D20。
在地震時程計算中設置兩種計算工況，分別為：ABAQUS-C3D8、UEL-C3D8。
計算結果表明，自編UEL與ABAQUS自帶單元結果一致。
（）模型信息
Koyna混凝土重力壩位于印度孟買東南200 km處，1967年12月11日，Koyna混凝土重力壩遭受里氏6.5級的地震（Koyna地震），該地震給大壩和水電站等水工建筑帶來了巨大的損壞，給下游數十萬居民的人身生命財產安全造成了巨大損失，該大壩地震案例事后成為諸多學者進行地震作用下壩基動力相互作用、混凝土材料的動態力學性能等領域研究的對象。
Koyna混凝土重力壩的橫斷面尺寸如下圖：
該壩壩高113 m，壩頂寬度14.8 m，壩底寬度70 m，正常蓄水位91.75 m，壩基尺寸分別取向上游、下游延伸兩倍壩高，向地基方向同樣延伸兩倍壩高，向橫河向延伸20 m，最終的三維幾何模型示意圖如下圖：
采用六面體單元離散，有限元計算模型如下圖：
該模型的壩體和壩基共計用17950個六面體單元，其中壩體5000個，壩基12950個，壩基的網格采用疏密漸進過渡的方式避免計算結果在網格突變處不連續。
展開 c30~c45塑性損傷本構ABAQUS輸入 ￥1
快速查詢混凝土彈模和拉壓強度標準值.xlsx
以下內容也可在我發布的文檔中免費下載
abaqus鋼管混凝土的本構關系輸入
各位大佬，請問abaqus里面關于鋼管混凝土的核心混凝土的本構關系怎么輸入？是按照一般的分成彈性和塑性分別輸入嗎，但是彈性模量不知道怎么取。本構關系是根據韓林海的鋼管混凝土結構計算的。
