ansys 地震計算的案例
地震波的90%的python腳本持時計算 ￥2
話不多說，先展示計算后的結(jié)果示意圖，包括下面給的代碼，運行計算后都會給出相應(yīng)的持時圖！
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
# 設(shè)置中文字體和負號顯示
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
# 檢查數(shù)據(jù)采樣是否均勻
dt = np.diff(time)
if not np.allclose(dt, dt[0], atol=1e-5):
raise ValueError("時間列采樣不均勻，需預(yù)處理為均勻采樣數(shù)據(jù)")
dt = dt[0] # 時間步長
展開 【CAE案例】支持地震風(fēng)險評估的土壤-結(jié)構(gòu)相互作用的概率計算
1979年，第一次完成了地震的風(fēng)險評估，在1990年法國也完成了第一次EPS。
EPS方法需要的輸入為：
1、初始事件（電路故障、地震等）；
2、系統(tǒng)的組成部分；
3、他們的不確定性。
輸出為：
1、最終的風(fēng)險；
2、各組成部分對最終風(fēng)險的影響。
圖一 EPS流程示意圖
02 研究方法
易損性曲線不僅隨著平板的加速度的變化而變化，也和土壤的加速度有關(guān)，如下圖二所示。為了求出不確定性曲線，研究人員使用結(jié)構(gòu)仿真進行計算。有兩種計算方法，第一，通過ISS（Interaction Soil-Structure）和各組成部分進行易損性的完整計算；第二，分開計算，首先通過ISS計算得到平板的頻譜，其次進行各組分的易損性計算。在這兩種情況下都通過結(jié)構(gòu)仿真和Miss3D進行ISS的概率計算。
圖二 物體不確定性組成
拉丁超立方抽樣作為一種不確定性傳播的概率方法有以下幾個優(yōu)勢：首先對參數(shù)空間有良好的代表性，因為其以等概率間隔繪制并根據(jù)超立方體隨機組合；其次它對線性情況的收斂性較好，只需要大約三十次循環(huán)計算就足夠。圖三為拉丁超立方抽樣的示意圖。
圖三 拉丁超立方抽樣示意圖
研究人員得到了以下兩個模型，鏈狀模型和板狀模型，示意圖如下圖四。
圖四 鏈狀模型（左）和板狀模型（中、右）
計算前準備五個隨機變量，類型有三種：1、土壤的加速度；2、結(jié)構(gòu)的楊氏模量和阻尼；3、土壤的楊氏模量和阻尼，其中土壤的楊氏模量和阻尼是相關(guān)變量，他們的關(guān)系如下圖五所示。
展開 ABAQUS UEL 二次開發(fā)（Koyna混凝土壩地震響應(yīng)計算）
概述：以Koyna混凝土壩為對象進行地震響應(yīng)計算。將自編的八節(jié)點UEL和二十節(jié)點UEL應(yīng)用到計算中。分別進行了混凝土壩模態(tài)計算和地震時程計算。
其中，在模態(tài)計算中共設(shè)置四種計算工況，分別為：ABAQUS-C3D8、UEL-C3D8、ABAQUS-C3D20、UEL-C3D20。
在地震時程計算中設(shè)置兩種計算工況，分別為：ABAQUS-C3D8、UEL-C3D8。
計算結(jié)果表明，自編UEL與ABAQUS自帶單元結(jié)果一致。
（）模型信息
Koyna混凝土重力壩位于印度孟買東南200 km處，1967年12月11日，Koyna混凝土重力壩遭受里氏6.5級的地震（Koyna地震），該地震給大壩和水電站等水工建筑帶來了巨大的損壞，給下游數(shù)十萬居民的人身生命財產(chǎn)安全造成了巨大損失，該大壩地震案例事后成為諸多學(xué)者進行地震作用下壩基動力相互作用、混凝土材料的動態(tài)力學(xué)性能等領(lǐng)域研究的對象。
Koyna混凝土重力壩的橫斷面尺寸如下圖：
該壩壩高113 m，壩頂寬度14.8 m，壩底寬度70 m，正常蓄水位91.75 m，壩基尺寸分別取向上游、下游延伸兩倍壩高，向地基方向同樣延伸兩倍壩高，向橫河向延伸20 m，最終的三維幾何模型示意圖如下圖：
采用六面體單元離散，有限元計算模型如下圖：
該模型的壩體和壩基共計用17950個六面體單元，其中壩體5000個，壩基12950個，壩基的網(wǎng)格采用疏密漸進過渡的方式避免計算結(jié)果在網(wǎng)格突變處不連續(xù)。
展開 ANSYS知識普及系列16——在ANSYS里施加地震慣性力的方法
本人準備出一個ANSYS知識普及系列，將有用的網(wǎng)上資料歸攏，由于知識水平有限，不對之處請諒解。也歡迎各位網(wǎng)友提供好的資料分享，讓我們共同完成這個ANSYS知識普及系列。
編輯人：技術(shù)鄰ANSYS專家
業(yè)務(wù)咨詢網(wǎng)址：http://www.yqgqt.org.cn/content/other/402981
（打個小廣告）
聲 明:1、ANSYS知識普及系列中所有資料均來自網(wǎng)上；
2、如侵犯知識產(chǎn)權(quán)，請聯(lián)系ANSYS專家本人或者技術(shù)鄰，我將第一時間刪除。
小技巧：加本人關(guān)注，可以及時觀看本人發(fā)布的技術(shù)貼
在ANSYS里做地震分析時，需要對結(jié)構(gòu)施加地震慣性荷載，地震慣性力是通過加速度的方式輸入進結(jié)構(gòu)的，然后與結(jié)構(gòu)的質(zhì)量一起形成動力計算時的慣性荷載，下面說一下在ANSYS里施加地震慣性力的方法。
展開 
BAQUS反應(yīng)譜法計算地震反應(yīng)的簡單實例＋時程分析
使用《有限元法及其應(yīng)用》中的ansys算例的問題進行說明，供大家參考。
補充時程分析cae操作
ABAQUS反應(yīng)譜法計算地震反應(yīng)的簡單實例＋時程分析3.rar
ABAQUS反應(yīng)譜法計算地震反應(yīng)的簡單實例＋時程分析1.rar
ABAQUS反應(yīng)譜法計算地震反應(yīng)的簡單實例＋時程分析2.rar
ABAQUS反應(yīng)譜法計算地震反應(yīng)的簡單實例＋時程分析
琢磨了下ABAQUS如何進行地震反應(yīng)譜計算。使用《有限元法及其應(yīng)用》中的ansys算例的問題進行說明，供大家參考。
補充時程分析cae操作
ABAQUS地震時程分析小算例1-4.rar
demo-spc.rar
ABAQUS地震反應(yīng)譜分析1-4.rar
ansys地震分析圖片
ansys地震分析圖片
ANSYS地震分析算例
一個四層彈簧－質(zhì)點模型的地震ANSYS分析
! Example of seismic analysis in ANSYS
! 包括頻率分析，譜分析和時程分析
! With model & frequency analysis, spectrum analysis and transient analysis
! 作者：陸新征，清華大學(xué)土木工程系
! Author: Lu Xinzheng Dept. Civil Engrg. of Tsinghua University
[Money=20]
FINI
/CLEAR
/UNITS,SI
/PREP7
DAMPRATIO=0.02 !
展開 ANSYS workbench石油井架地震響應(yīng)分析 ￥10
本案例適合哪些人學(xué)習(xí)：
1、學(xué)習(xí)型仿真工程師
2、理工科院校學(xué)生
你會得到什么：
1、學(xué)習(xí)石油井架模型的三維模型處理
2、學(xué)習(xí)地震響應(yīng)分析相關(guān)的分析步的建立
3、學(xué)習(xí)地震響應(yīng)分析相關(guān)的約束條件的建立
4、學(xué)習(xí)地震響應(yīng)分析的載荷施加
案例介紹：
所使用軟件為ANSYS workbench2020R2.
案例介紹了ANSYS workbench 石油井架地震響應(yīng)分析。
本案例完整得提供了分析相關(guān)所有分析文件。
?
ANSYS強度折減法邊坡穩(wěn)定性分析及地震荷載分析 ￥30
采用ANSYS有限元強度折減方法對滑坡穩(wěn)定系數(shù)進行求解，通過有限元強度折減方法對不同工況下滑坡穩(wěn)定系數(shù)進行計算，并將模擬計算值與極限平衡方法進行對比，驗證了強度折減方法的有效性。
有限元強度折減法是20世紀70年代末由英國科學(xué)家Zienkiewicz提出的，是通過不斷提高強度折減系數(shù)來降低坡體巖土抗剪強度參數(shù)，并反復(fù)試算，直到達到極限破壞狀態(tài)，程序自動根據(jù)彈塑性有限元計算結(jié)果得到滑動破壞面，同時得到滑坡的強度儲備安全系數(shù)。該方法在理論體系上比極限平衡法更嚴格，它全面滿足了靜力許可、應(yīng)變相容以及土體的非線性應(yīng)力-應(yīng)變關(guān)系。
地震荷載加載前需要對模型進行模態(tài)分析求解，來獲得固有頻率及瑞麗阻尼系數(shù)，然后再對模型進行動態(tài)加載。
第一步：模型建立、施加邊界條件、自重工況下強度折減
第二步：模態(tài)分析求解
第三步：求解瑞麗阻尼系數(shù)、地震波加載
展開 新論文：新型地震和連續(xù)倒塌綜合防御韌性PC框架承載力計算方法
對于土木工程中最常見的混凝土框架結(jié)構(gòu)而言，諸多研究表明由地震引起的倒塌和結(jié)構(gòu)局部破壞導(dǎo)致的連續(xù)倒塌是混凝土框架結(jié)構(gòu)的兩大主要破壞形式。因此我們提出了一種新型的地震和連續(xù)倒塌綜合防御韌性PC框架（multi-hazard resilient prefabricated concrete frame ，MHRPC框架），詳情請參考：
新論文：這個混凝土框架能抗震，能防連續(xù)倒塌，還功能可恢復(fù)，您不進來看看么？
地震與連續(xù)倒塌綜合防御韌性PC框架結(jié)構(gòu)，如圖1所示：
(1) 框架梁和框架柱通過剪力傳遞板傳遞剪力，通過可更換耗能裝置和預(yù)應(yīng)力筋傳遞彎矩。
(2) 預(yù)應(yīng)力筋可以同時作為抗震的自復(fù)位鋼筋和抗連續(xù)倒塌的拉結(jié)配筋；
(3) 剪力傳遞板保障大變形下剪力的可靠傳遞；
(4) 可更換耗能裝置可以消耗地震和連續(xù)倒塌作用下的動能。
圖1 MHRPC框架結(jié)構(gòu)
為了驗證該體系的效果，我們同時開展了抗震和防連續(xù)倒塌子結(jié)構(gòu)試驗，簡化版試驗結(jié)果（與常規(guī)RC框架試件對比）
圖2. 抗震性能試驗
（MHRPC框架：承載力高，二階剛度穩(wěn)定，殘余變形小）
圖3.
展開 
ANSYS橋梁地震時程分析
在對斜拉橋施加300gal的ELCentro水平方向的地震波后，提取了跨中位置點的水平和垂直位移，可是水平方向的位移幾乎為0，最大位移只有4.14775e-5m，而垂直方向上也產(chǎn)生了位移并且位移最大值達到2m，并且最終跨中點停留在了1m處，沒有恢復(fù)到原來的起始位置（0m），提取了橋墩和橋塔的點也發(fā)現(xiàn)了同樣的問題，提取點水平方向位移很小，垂直方向位移很大，并且最終停留的位置也各不相同，水平方向的地震波按理不會對豎直方向產(chǎn)生這么大的影響。顯然是計算出現(xiàn)了問題，可是現(xiàn)在找不到具體是哪里的問題，軟件沒有報錯，約束應(yīng)該也沒有少，請各位大佬幫忙看看是哪里出現(xiàn)了問題，謝謝了！
跨中位置水平方向位移
跨中位置垂直方向位移
梁柱交接點水平方向位移
梁柱交接點垂直方向位移
斜拉橋模型
展開 ANSYS框架結(jié)構(gòu)地震時程分析
主要闡述了地震波選波的一些關(guān)鍵點，如何根據(jù)設(shè)計反應(yīng)譜人工生成地震波，ANSYS讀入地震波的方法以及計算結(jié)果的輸出方法以及其他的一些相關(guān)技巧。
ansys之——地震波的輸入和求解
對于地震波的輸入，可以把荷載記錄做成文件，利用apdl的讀取功能讀入倒數(shù)據(jù)庫中。下面的例子是自己編的一個小文件。修改一下可以更簡潔。有用到的朋友自己作一下把。
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數(shù)據(jù)庫中作為數(shù)組。
也可以把加速度記錄做成一個文件，這樣程序就簡單多了。大家可以試看看修改一下。
下面是計算部分語句：
/SOLU
ANTYPE,trans
！
展開 框架結(jié)構(gòu) 剪力墻結(jié)構(gòu) ABAQUS水平地震力計。算 ,push over ，動力時程計算 ，反應(yīng)譜
框架結(jié)構(gòu) 剪力墻結(jié)構(gòu) ABAQUS水平地震力計。算 ,push over ，動力時程計算 ，反應(yīng)譜(二維建模，三維建模)
