abaqus地震計算的案例
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個，壩基的網格采用疏密漸進過渡的方式避免計算結果在網格突變處不連續。
展開 ABAQUS反應譜法計算地震反應的簡單實例＋時程分析
琢磨了下ABAQUS如何進行地震反應譜計算。使用《有限元法及其應用》中的ansys算例的問題進行說明，供大家參考。
補充時程分析cae操作
ABAQUS地震時程分析小算例1-4.rar
demo-spc.rar
ABAQUS地震反應譜分析1-4.rar
框架結構 剪力墻結構 ABAQUS水平地震力計。算 ,push over ，動力時程計算 ，反應譜
框架結構 剪力墻結構 ABAQUS水平地震力計。算 ,push over ，動力時程計算 ，反應譜(二維建模，三維建模)
abaqus 自重+地震，其中地震用反應譜分析，如何做？
abaqus 自重+地震，其中地震用反應譜分析，如何做？
地震波的90%的python腳本持時計算 ￥2
話不多說，先展示計算后的結果示意圖，包括下面給的代碼，運行計算后都會給出相應的持時圖！
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
# 設置中文字體和負號顯示
plt.rcParams['font.sans-serif'] = ['SimHei'] # 顯示中文
plt.rcParams['axes.unicode_minus'] = False # 正常顯示負號
# 讀取 Excel 文件
file_path = 'x.xlsx' # 替換為你的文件路徑
df = pd.read_excel(file_path, engine='openpyxl')
# 提取時間和加速度列
time = df['Time'].values
acc = df['Acceleration'].values
# 檢查數據采樣是否均勻
dt = np.diff(time)
if not np.allclose(dt, dt[0], atol=1e-5):
raise ValueError("時間列采樣不均勻，需預處理為均勻采樣數據")
dt = dt[0] # 時間步長
展開 【CAE案例】支持地震風險評估的土壤-結構相互作用的概率計算
1979年，第一次完成了地震的風險評估，在1990年法國也完成了第一次EPS。
EPS方法需要的輸入為：
1、初始事件（電路故障、地震等）；
2、系統的組成部分；
3、他們的不確定性。
輸出為：
1、最終的風險；
2、各組成部分對最終風險的影響。
圖一 EPS流程示意圖
02 研究方法
易損性曲線不僅隨著平板的加速度的變化而變化，也和土壤的加速度有關，如下圖二所示。為了求出不確定性曲線，研究人員使用結構仿真進行計算。有兩種計算方法，第一，通過ISS（Interaction Soil-Structure）和各組成部分進行易損性的完整計算；第二，分開計算，首先通過ISS計算得到平板的頻譜，其次進行各組分的易損性計算。在這兩種情況下都通過結構仿真和Miss3D進行ISS的概率計算。
圖二 物體不確定性組成
拉丁超立方抽樣作為一種不確定性傳播的概率方法有以下幾個優勢：首先對參數空間有良好的代表性，因為其以等概率間隔繪制并根據超立方體隨機組合；其次它對線性情況的收斂性較好，只需要大約三十次循環計算就足夠。圖三為拉丁超立方抽樣的示意圖。
圖三 拉丁超立方抽樣示意圖
研究人員得到了以下兩個模型，鏈狀模型和板狀模型，示意圖如下圖四。
圖四 鏈狀模型（左）和板狀模型（中、右）
計算前準備五個隨機變量，類型有三種：1、土壤的加速度；2、結構的楊氏模量和阻尼；3、土壤的楊氏模量和阻尼，其中土壤的楊氏模量和阻尼是相關變量，他們的關系如下圖五所示。
展開 BAQUS反應譜法計算地震反應的簡單實例＋時程分析
補充時程分析cae操作
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析3.rar
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析1.rar
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析2.rar
新論文：新型地震和連續倒塌綜合防御韌性PC框架承載力計算方法
對于土木工程中最常見的混凝土框架結構而言，諸多研究表明由地震引起的倒塌和結構局部破壞導致的連續倒塌是混凝土框架結構的兩大主要破壞形式。因此我們提出了一種新型的地震和連續倒塌綜合防御韌性PC框架（multi-hazard resilient prefabricated concrete frame ，MHRPC框架），詳情請參考：
新論文：這個混凝土框架能抗震，能防連續倒塌，還功能可恢復，您不進來看看么？
地震與連續倒塌綜合防御韌性PC框架結構，如圖1所示：
(1) 框架梁和框架柱通過剪力傳遞板傳遞剪力，通過可更換耗能裝置和預應力筋傳遞彎矩。
(2) 預應力筋可以同時作為抗震的自復位鋼筋和抗連續倒塌的拉結配筋；
(3) 剪力傳遞板保障大變形下剪力的可靠傳遞；
(4) 可更換耗能裝置可以消耗地震和連續倒塌作用下的動能。
圖1 MHRPC框架結構
為了驗證該體系的效果，我們同時開展了抗震和防連續倒塌子結構試驗，簡化版試驗結果（與常規RC框架試件對比）
圖2. 抗震性能試驗
（MHRPC框架：承載力高，二階剛度穩定，殘余變形小）
圖3.
展開 ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬 ￥70
ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬
一、建模技術
地震工況下邊坡可能失穩進而出現滑坡現象，為避免模擬滑坡時網格產生的畸變問題，采用耦合歐拉拉格朗日法(CEL)進行滑坡的大變形模擬；土體本構采用摩爾庫倫模型；采用模型底部小范圍內的周期性荷載模擬地震荷載。
二、模型及部分結果展示
圖1：藍色為邊坡；紅色為空氣層
圖2：網格的劃分
圖3：賦予模型初始應力
圖4：土體達到地應力平衡時的應力分布
圖5：土體底部的地震荷載施加區域
圖6：所施加的周期性荷載（地震荷載）
圖7：邊坡因地震荷載產生的位移
圖8：地震波產生的區域
展開 abaqus地震方面的資料以及總結 ￥2
abaqus地震方面的資料以及總結
abaqus地震分析關于模型漂移問題
關于解決地震分析模型漂移的方法，某前輩在博客已經做了比較完美的解答，即基線修正法，http://blog.sina.com.cn/s/blog_417a7d700100z77g.html但這個方法只適用于隱式分析。如果模型是顯式分析該怎么解決呢？在顯式分析中如果添加了基線修正，就會提示keyword在顯式分析中不可獲得。自己的模型屬于車橋耦合振動，比較復雜，隱式計算太耗費時間，所以采用了顯式計算。有沒有在顯式分析的前提下解決模型漂移的辦法呢？求助.........
也歡迎大家在此帖下討論
四層RC框架開洞，八度多遇地震下底部剪力頂層位移?模態分析，push over，時程分析,反應譜計算 ￥100
采用ABAQUS軟件計算RC框架水平地震。工程概況：四層RC框架，梁截面400*200，配筋4根HRB400直接20mm；柱截面600*600，配筋8根HRB400直接25mm；柱間距4200mm，樓板大開洞。計算八度多遇地震下底部剪力及頂層位移?(方法：push over，時程分析，反應譜計算作對比)
提共模型(DAT,CAE,INP,ODB文件)，計算小插件，計算結果，以及計算截圖!
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混凝土壩地震動力響應分析）
摘要：以平面線性四節點單元為例，分別采用有限元法和比例邊界有限元法（SBFEM）在ABAQUS提供的UEL子程序接口進行二次開發，編寫的UEL均包含動力計算部分，即采用HHT隱式時程積分法求解動力方程。將ABAQUS自帶的CPS4單元、自編四節點等參單元和自編SBFEM的UEL三者進行對比。將以上三種單元應用到Koyna混凝土壩地震動力響應分析中，對比壩體關鍵點數據，驗證三種單元的計算結果吻合良好。
（一）模型基本信息
（1）材料信息
壩體彈性模量：E = 31027 MPa，泊松比：0.15，密度：2643 kg/m3。
壩體尺寸
（2）網格信息
模型網格
采用四節點單元離散壩體，共計1891個節點，1800個單元。
（二）Koyna混凝土壩模態分析
模態分析時將壩體底部設置固定邊界，約束雙向位移。
“ABAQUS”代表軟件自帶的四節點單元計算結果，“ABAQUS-CPS4”指的是用UEL實現的四節點單元計算結果，“UEL-SBFEM”指的是用UEL實現的SBFE單元計算結果，與“Chopra and Chakrabarti (1973)”的結果進行對比，可以看到三者計算精度基本保持一致。
展開 ABAQUS中粘彈性人工邊界及地震力的施加
Add選項定義阻尼Damping，力F為水平方向F1，阻尼系數C11由之前 計算確定。
③賦予屬性：菜單欄-Connector-assignment-create，選取①中定義的set和②中定義的屬性，進行賦值。賦值后模型底部如：
7、創建彈簧
彈簧的創建根據不同的土體性質而不同，假設全是融土，則分地基左右側彈簧和填土右側彈簧兩種。
菜單欄-special-springs/dashpots-create選擇接地彈簧
點Mesh選擇兩側節點，并定義彈簧剛度，計算由 確定。自由度選擇水平方向1方向（彈簧水平方向震動）
定義之后：
同樣方式定義填土右側彈簧剛度。
8、邊界條件與荷載
關掉填土和地基左右側水平方向位移約束（因為三個面已定義彈簧，另一面與墻背摩擦），但底部豎向位移邊界條件不可關閉，否則模型會在地震力作用下飛走……
定義地震加速度積分得到的速度，并創建Amp-02g-v（最大0.2g加速度為例）。定義地震力為集中力，施加到底部節點上，數值大小由 定義。【地震過程中節點力是不同時刻速度v的 倍】
9、創建job并提交
10、后處理提取墻頂加速度，位移，墻背土壓力（創建path），總土壓力（）
展開 