地震波反應譜的案例
基于matlab編寫單自由度的地震波的反應譜曲線，及規范反應譜曲線 ￥9.9
程序的操作過程：
首先是地震波的反應譜
運行MATLAB程序之后會出現這個
輸入地震波后，如下輸入阻尼比，質量，加速度幅值后
自動生成圖形第二是規范反應譜的生成
如圖，運行程序，輸入設計的烈度，場地類型，設計地震分組，阻尼比
生成如圖，還可以在文件中提取橫坐標和縱坐標，方便在peer中選取地震波兩個程序具體如下
如何利用反應譜生成地震波
然后點擊SIMQKE;
生成如下的界面：
在上圖中的“SIMQKE”參數欄中輸入參數，各參數說明如下：
TS：最小的周期，導入的反應譜文件中周期比 TS 小的點都將被忽略
TL：最大的周期，導入的反應譜文件中周期比 TL 大的點都將被忽略，這里的值一定不要大于反應譜的周期(這里TL取值一定小于等于3，否則軟件會閃退)
TRISE：加速度平穩段的開始時間
TLVL：加速度平穩段的持續時間，最小值為 10s
DUR：輸出的時程總的持續時間，如 DUR=20，則輸出的時程為 20s.
NCYCLE：迭代的次數+1，如果 NCYCLE=1 則不進行迭代，使得在計算過
程中得到的反應譜曲線更加平滑
AGMAX：峰值加速度，由程序自動計算(這里的單位是g，跟我們輸入的是一樣的)
NPA：要生成的時程波條數，這里選擇 3 條
IIX：任意奇數，缺省值為 1235
AMOR：反應譜的阻尼比，這里選擇 0.05
===============================================================
然后點擊Simulated Earthquake,生成地震波如下圖所示：
查看相應的數據：
導出的數據如下所示：
這期的利用反應譜生成地震波的教程就結束了，這個軟件網上可以直接下載，也可以關注公眾號“生活中的力學仿真”，后臺回復“地震波生成軟件”獲取。
展開 大跨度鋼結構地震過程仿真分析
大跨度鋼結構地震過程仿真分析
近期，人們被日本九州地震和南美厄瓜多爾地震給震怕了，有人調侃“地球進入振動模式”，實際上從長期來看，世界范圍內的地震出現頻率和級別并無異常。在我們的科技水平還無法準確預測地震的時候，防震救災工作就顯得很重要，例如建筑設備的抗震能力能否抵抗住某次大級別地震，關系到人們生命安全問題。
以某個大跨度鋼結構為例，介紹一下CAE方法在地震分析中的應用。
地震分析方法一般分為底部剪力法、反應譜法和時程分析法。CAE領域常用反應譜法或時程分析法，時程分析法采用的載荷是某次具體的地震波，反應譜法是多個地震波的統計結果，由時域轉換成頻域而來，從計算的經濟性和通用性考慮，最常采用的是反應譜法。
地震波反應譜
學過力學的朋友都知道，地震波分為橫波和縱波，由于縱波的傳輸速度比橫波快，從振源首先到達地表，因此震中地區的人們首先感到的是上下震動，然后是橫向振動，豎向地震的加速度值比橫向要小，國標規定取橫向數值的65%，因此對地面設備造成最大破壞的是橫波的剪切作用。CAE分析過程也分為豎向地震和橫向地震。
在abaqus中采用反應譜法進行地震分析，需要先提取結構的頻率，然后再進行反應譜分析，同時輸入定義好的反應譜和阻尼，反應譜可以是加速度譜、位移譜、速度譜等內容：
提交計算后就可以得到結構在地震作用下的響應：
今年7月28日是家鄉唐山抗震40周年紀念日，在技術落后的年代，自然災害帶來的后果是慘痛的。
文章轉載自微信公眾號：SmartFEA
旨在分享知識，若侵即刪
展開 abaqus 自重+地震，其中地震用反應譜分析，如何做？
abaqus 自重+地震，其中地震用反應譜分析，如何做？
BAQUS反應譜法計算地震反應的簡單實例＋時程分析
補充時程分析cae操作
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析3.rar
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析1.rar
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析2.rar
ABAQUS反應譜法計算地震反應的簡單實例＋時程分析
琢磨了下ABAQUS如何進行地震反應譜計算。使用《有限元法及其應用》中的ansys算例的問題進行說明，供大家參考。
補充時程分析cae操作
ABAQUS地震時程分析小算例1-4.rar
demo-spc.rar
ABAQUS地震反應譜分析1-4.rar
框架結構 剪力墻結構 ABAQUS水平地震力計。算 ,push over ，動力時程計算 ，反應譜
框架結構 剪力墻結構 ABAQUS水平地震力計。算 ,push over ，動力時程計算 ，反應譜(二維建模，三維建模)
四層RC框架開洞，八度多遇地震下底部剪力頂層位移?模態分析，push over，時程分析,反應譜計算 ￥100
采用ABAQUS軟件計算RC框架水平地震。工程概況：四層RC框架，梁截面400*200，配筋4根HRB400直接20mm；柱截面600*600，配筋8根HRB400直接25mm；柱間距4200mm，樓板大開洞。計算八度多遇地震下底部剪力及頂層位移?(方法：push over，時程分析，反應譜計算作對比)
提共模型(DAT,CAE,INP,ODB文件)，計算小插件，計算結果，以及計算截圖!
PEER下載地震波方法介紹
下拉頁面，顯示搜索到的100條地震記錄，這些地震波的反應譜與剛才上傳的規范反應譜都存在不同程度的吻合，通過自行選擇，查找與規范反應譜吻合較好的地震波名稱。下拉網頁即可看到搜索的地震記錄。
經過篩選后的地震波(含有多條)，對于自己想要下載的地震波，不能直接在地震名稱前勾選下載，若這樣操作，則會將篩選的所有地震波的數據都下載下來，從下載的文件中也不易識別具體哪個是自己想要的數據；便捷的方法是明確自己做需要的地震波的名稱和測站，直接將這些信息輸入到上面的詳細篩選選項中，這樣重新篩選時，只會得到這一條地震波信息，勾選這一條地震波并下載它。
勾選地震波，點擊下載，下載的時候有兩種選擇，選擇后者，將地震波的時程下載下來。然后依次單個的下載地震波。
當然了，地震波下載后，還需要進行數據的轉換處理。本推文就不講解了，在下個推文中我再介紹相應的方法。
歡迎關注公眾號“土木愛研小站”
您的每一個贊和關注都是我前進的動力！！！
展開 【JY】結構動力學初步-單質點結構的瞬態動力學分析
數值分析算例
【例題】為了對比各軟件與理論計算結果的誤差和展示建模過程，以一個簡單的單質點體系作為算例，地震波選取三條人工波，單質點體系的自振周期為2.7s，阻尼比為0.05，提取質點的加速度及位移時程進行對比。
木質框架模型雙向地震仿真分析
4.3地震作用計算
采用隱式算法計算底部加速度時程作用下的結構響應，其基本思路如下：有限元分析分為兩個步驟，第一步采用Static General 步驟施加豎向重力荷載，模型底面采用固定約束；第二步采用Dynamic Implicit 步驟進行地震時程分析，釋放水平兩個方向的約束，并施加雙向地震波加速度時程。為驗證摩擦阻尼器的消能減震性能，進行了未設置阻尼器與設置阻尼器的兩個模型計算結果對比。兩個模型均作用峰值為250gal的雙向地震波時程曲線。圖8(a)為雙向250gal地震波，圖8(b)為6組地震波加速度反應譜。從圖9~10能夠看出，阻尼器的布置能明顯減小結構頂層位移和底部剪力。
五、阻尼力參數分析
為研究阻尼器屈服力對結構抗震性能的影響，分別計算了峰值加速度為250gal和500gal的地震荷載作用下，阻尼器屈服力Fd分別為24N（方案A）、30N（方案B）、36N（方案C）時結構各層加速度響應、位移響應、各層水平轉角響應。
從圖可以看出，在24N~36N范圍內，相同地震波作用下阻尼器屈服力的改變對各層峰值加速度影響不明顯。但是在大震下，結構的峰值加速度有所增加，并且WE方向的各層峰值加速度包絡圖的形狀發生了明顯的改變。
通過以上對框架模型的數值計算分析，可得到以下結論：
1）模型的數值計算及加載測試結果顯示，摩擦阻尼器能有效提高此結構模型的耗能能力，并且斜向棉繩拉索的設置，對結構的層間位移起到了良好的控制作用。
2）有限元計算結構各層的響應結果顯示，在一定范圍內隨著阻尼器屈服力增加，結構各層的位移響應降低。摩擦阻尼器屈服里的改變對結構加速度響應沒有明顯影響。
3）由于阻尼器的非對稱布置，數值計算顯示，阻尼器屈服力的增加，會造成樓層最大轉角增大。
展開 
【JY】動力學利器 —— JYdyn函數包分享與體驗
回調函數可以調用已經編寫好的m文件，這里將讀取地震波的主要代碼放在ReadData.m 文件中，同時將反應譜分析的函數放在SpectrumCal.m中，編寫好的回調函數如下圖所示。
一句命令做完了反應譜分析！
（嫌麻煩也可以直接用這個快速做反應譜出結果哦~）
第五步，點擊運行，可以得到如下界面，選擇一條地震波，運行便可得到地震波的反應譜，這樣一個簡單的繪制反應譜的APP便制作完成了。
其他更多更有趣的例子歡迎體驗和嘗試！
