ansys計算地震的案例
地震波的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、土壤的楊氏模量和阻尼，其中土壤的楊氏模量和阻尼是相關變量，他們的關系如下圖五所示。
展開 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個，壩基的網格采用疏密漸進過渡的方式避免計算結果在網格突變處不連續。
展開 ANSYS知識普及系列16——在ANSYS里施加地震慣性力的方法
本人準備出一個ANSYS知識普及系列，將有用的網上資料歸攏，由于知識水平有限，不對之處請諒解。也歡迎各位網友提供好的資料分享，讓我們共同完成這個ANSYS知識普及系列。
編輯人：技術鄰ANSYS專家
業務咨詢網址：http://www.yqgqt.org.cn/content/other/402981
（打個小廣告）
聲 明:1、ANSYS知識普及系列中所有資料均來自網上；
2、如侵犯知識產權，請聯系ANSYS專家本人或者技術鄰，我將第一時間刪除。
小技巧：加本人關注，可以及時觀看本人發布的技術貼
在ANSYS里做地震分析時，需要對結構施加地震慣性荷載，地震慣性力是通過加速度的方式輸入進結構的，然后與結構的質量一起形成動力計算時的慣性荷載，下面說一下在ANSYS里施加地震慣性力的方法。
展開 
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析
琢磨了下ABAQUS如何進行地震反應譜計算。使用《有限元法及其應用》中的ansys算例的問題進行說明，供大家參考。
補充時程分析cae操作
ABAQUS地震時程分析小算例1-4.rar
demo-spc.rar
ABAQUS地震反應譜分析1-4.rar
BAQUS反應譜法計算地震反應的簡單實例＋時程分析
使用《有限元法及其應用》中的ansys算例的問題進行說明，供大家參考。
補充時程分析cae操作
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析3.rar
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析1.rar
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析2.rar
ANSYS地震分析算例
一個四層彈簧－質點模型的地震ANSYS分析
! Example of seismic analysis in ANSYS
! 包括頻率分析，譜分析和時程分析
! With model & frequency analysis, spectrum analysis and transient analysis
! 作者：陸新征，清華大學土木工程系
! Author: Lu Xinzheng Dept. Civil Engrg. of Tsinghua University
[Money=20]
FINI
/CLEAR
/UNITS,SI
/PREP7
DAMPRATIO=0.02 !
展開 ANSYS workbench石油井架地震響應分析 ￥10
本案例適合哪些人學習：
1、學習型仿真工程師
2、理工科院校學生
你會得到什么：
1、學習石油井架模型的三維模型處理
2、學習地震響應分析相關的分析步的建立
3、學習地震響應分析相關的約束條件的建立
4、學習地震響應分析的載荷施加
案例介紹：
所使用軟件為ANSYS workbench2020R2.
案例介紹了ANSYS workbench 石油井架地震響應分析。
本案例完整得提供了分析相關所有分析文件。
?
ansys地震分析圖片
ansys地震分析圖片
ANSYS強度折減法邊坡穩定性分析及地震荷載分析 ￥30
采用ANSYS有限元強度折減方法對滑坡穩定系數進行求解，通過有限元強度折減方法對不同工況下滑坡穩定系數進行計算，并將模擬計算值與極限平衡方法進行對比，驗證了強度折減方法的有效性。
有限元強度折減法是20世紀70年代末由英國科學家Zienkiewicz提出的，是通過不斷提高強度折減系數來降低坡體巖土抗剪強度參數，并反復試算，直到達到極限破壞狀態，程序自動根據彈塑性有限元計算結果得到滑動破壞面，同時得到滑坡的強度儲備安全系數。該方法在理論體系上比極限平衡法更嚴格，它全面滿足了靜力許可、應變相容以及土體的非線性應力-應變關系。
地震荷載加載前需要對模型進行模態分析求解，來獲得固有頻率及瑞麗阻尼系數，然后再對模型進行動態加載。
第一步：模型建立、施加邊界條件、自重工況下強度折減
第二步：模態分析求解
第三步：求解瑞麗阻尼系數、地震波加載
展開 ANSYS框架結構地震時程分析
主要闡述了地震波選波的一些關鍵點，如何根據設計反應譜人工生成地震波，ANSYS讀入地震波的方法以及計算結果的輸出方法以及其他的一些相關技巧。
新論文：新型地震和連續倒塌綜合防御韌性PC框架承載力計算方法
對于土木工程中最常見的混凝土框架結構而言，諸多研究表明由地震引起的倒塌和結構局部破壞導致的連續倒塌是混凝土框架結構的兩大主要破壞形式。因此我們提出了一種新型的地震和連續倒塌綜合防御韌性PC框架（multi-hazard resilient prefabricated concrete frame ，MHRPC框架），詳情請參考：
新論文：這個混凝土框架能抗震，能防連續倒塌，還功能可恢復，您不進來看看么？
地震與連續倒塌綜合防御韌性PC框架結構，如圖1所示：
(1) 框架梁和框架柱通過剪力傳遞板傳遞剪力，通過可更換耗能裝置和預應力筋傳遞彎矩。
(2) 預應力筋可以同時作為抗震的自復位鋼筋和抗連續倒塌的拉結配筋；
(3) 剪力傳遞板保障大變形下剪力的可靠傳遞；
(4) 可更換耗能裝置可以消耗地震和連續倒塌作用下的動能。
圖1 MHRPC框架結構
為了驗證該體系的效果，我們同時開展了抗震和防連續倒塌子結構試驗，簡化版試驗結果（與常規RC框架試件對比）
圖2. 抗震性能試驗
（MHRPC框架：承載力高，二階剛度穩定，殘余變形小）
圖3.
展開 ANSYS橋梁地震時程分析
在對斜拉橋施加300gal的ELCentro水平方向的地震波后，提取了跨中位置點的水平和垂直位移，可是水平方向的位移幾乎為0，最大位移只有4.14775e-5m，而垂直方向上也產生了位移并且位移最大值達到2m，并且最終跨中點停留在了1m處，沒有恢復到原來的起始位置（0m），提取了橋墩和橋塔的點也發現了同樣的問題，提取點水平方向位移很小，垂直方向位移很大，并且最終停留的位置也各不相同，水平方向的地震波按理不會對豎直方向產生這么大的影響。顯然是計算出現了問題，可是現在找不到具體是哪里的問題，軟件沒有報錯，約束應該也沒有少，請各位大佬幫忙看看是哪里出現了問題，謝謝了！
跨中位置水平方向位移
跨中位置垂直方向位移
梁柱交接點水平方向位移
梁柱交接點垂直方向位移
斜拉橋模型
展開 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
！
展開 框架結構 剪力墻結構 ABAQUS水平地震力計。算 ,push over ，動力時程計算 ，反應譜
框架結構 剪力墻結構 ABAQUS水平地震力計。算 ,push over ，動力時程計算 ，反應譜(二維建模，三維建模)
