1 包含的內(nèi)容

（1）說明文本

（2）有限元模型及建模命令流

（3）模態(tài)分析全過程命令流

（4）EL Centro地震波詳細數(shù)據(jù)

（5）動力時程分析全過程命令流

（6）節(jié)點響應(yīng)后處理命令流

（7）完整算例文件

（8）《ANSYS結(jié)構(gòu)動力分析與應(yīng)用》

2 研究背景

在突如其來的地震面前，建筑結(jié)構(gòu)的每一次晃動，都是對工程師設(shè)計理念與分析方法的終極拷問。結(jié)構(gòu)是否具備足夠的延性？振動能否有效耗散？我們該如何預(yù)判這些動態(tài)響應(yīng)，做出科學(xué)決策？在現(xiàn)代結(jié)構(gòu)抗震設(shè)計中，有限元分析已成為工程師手中的核心工具。其中，ANSYS憑借其強大的建模能力與數(shù)值分析引擎，成為進行地震響應(yīng)模擬與結(jié)構(gòu)動力評估的主流平臺之一。然而，從構(gòu)建模型到輸入地震波、從模態(tài)分析到時程響應(yīng)，整個流程對初學(xué)者而言既嚴(yán)謹又復(fù)雜，亟需系統(tǒng)的操作指南。

作為一名科研博主，我希望通過這份教程，為你梳理出一條抗震建模之路。你將學(xué)到：如何搭建高層建筑的簡化有限元模型；如何進行模態(tài)分析與阻尼建模；如何輸入真實地震波并施加慣性力；如何提取關(guān)鍵節(jié)點的時程響應(yīng)數(shù)據(jù)；以及，如何一步步將“地震”變?yōu)椤皵?shù)據(jù)”，讓結(jié)構(gòu)的抗震能力變得可視、可量化、可優(yōu)化。無論你是結(jié)構(gòu)工程新手，還是希望將抗震仿真引入科研項目的研究者，這份教程都將成為你邁向工程抗震仿真實踐的重要起點。

3 研究的依據(jù)

[1] 王新敏. ANSYS結(jié)構(gòu)動力分析與應(yīng)用[M]. 人民交通出版社, 2014.

4 算例有限元模型

本模型采用ANSYS命令流構(gòu)建了一個典型的20層鋼筋混凝土高層框架結(jié)構(gòu)，旨在分析其在重力與地震荷載作用下的力學(xué)響應(yīng)。結(jié)構(gòu)主要特征如下：

（1）結(jié)構(gòu)形式：三維矩形平面框架，由梁柱構(gòu)件組成，不含剪力墻和樓板，以簡化分析。

（2）建模方法：使用ANSYS中的BEAM188單元模擬梁柱，具備考慮剪切變形與彎曲的能力，適合模擬細長框架構(gòu)件。

（3）空間布局：每層設(shè)4個節(jié)點（對應(yīng)柱腳和梁端），形成一個6m×6m的標(biāo)準(zhǔn)柱網(wǎng)，層高為3m，總高度為60m。

（4）材料參數(shù)：混凝土等效材料，彈性模量為30GPa，泊松比為0.2，密度為2500kg/m3。

（5）截面參數(shù)：梁柱截面均為0.4m×0.4m的矩形。

（6）邊界條件：底層所有節(jié)點完全固定，模擬剛性基礎(chǔ)。

（7）加載方式：重力荷載：通過加速度場模擬；地震荷載：在結(jié)構(gòu)上施加慣性力。

圖 1 鋼筋混凝土高層框架結(jié)構(gòu)有限元模型

5 模態(tài)分析

本分析采用ANSYS的命令流方式對結(jié)構(gòu)進行模態(tài)分析，以獲取其前10階固有頻率和振型。分析過程包括以下幾個步驟：

（1）設(shè)置分析類型：將分析類型指定為模態(tài)分析，以便求解結(jié)構(gòu)的固有頻率和振型。

（2）求解參數(shù)設(shè)定：

采用Lanczos方法，這是一種高效的特征值求解方法，適用于大型稀疏系統(tǒng)；指定提取前10階模態(tài)；設(shè)置頻率范圍，以排除非工程相關(guān)的高階模態(tài)；選用一致質(zhì)量矩陣，以確保模態(tài)結(jié)果的精度；采用PCG迭代求解器提高計算效率。

（3）模態(tài)歸一化與輸出控制：

設(shè)置MXPAND命令，展開模態(tài)振型并進行質(zhì)量歸一化，便于后續(xù)分析；不進行模態(tài)后應(yīng)力計算，以簡化分析輸出；不保存模態(tài)結(jié)果文件，用于快速試算階段。

（4）求解與完成：通過SOLVE命令啟動計算，最終輸出結(jié)構(gòu)的固有頻率和對應(yīng)振型信息，用于后續(xù)動力響應(yīng)分析。

圖 2 模態(tài)分析結(jié)果

6 動力時程分析

6.1 地震波選取

在PEER強震數(shù)據(jù)庫中選取典型的EL Centro地震波作為動力輸入，加速度時程如下圖。

圖 3 EL Centro地震波

6.2 在ANSYS中施加地震慣性力

本分析采用ANSYS平臺進行結(jié)構(gòu)的地震動響應(yīng)時程分析，模擬結(jié)構(gòu)在地震波作用下的動態(tài)響應(yīng)特性。分析流程如下：

（1）地震波數(shù)據(jù)讀取

從外部文件ElCentro.txt中讀取實測或合成的地震加速度時程數(shù)據(jù)，數(shù)據(jù)格式為四列（時間，X、Y、Z方向加速度），共2688個時間步。通過*VREAD語句將其導(dǎo)入數(shù)組，并自動計算時間間隔dt。

（2）瑞利阻尼系數(shù)設(shè)定

采用經(jīng)典的雙頻點瑞利阻尼方法，根據(jù)結(jié)構(gòu)第一階（0.708577Hz）與第二階（7.63773Hz）固有頻率計算出阻尼系數(shù)α與β，對質(zhì)量項與剛度項進行阻尼控制，阻尼比設(shè)定為5%。

（3）分析類型與控制參數(shù)設(shè)置

分析類型為瞬態(tài)動力分析，使用直接積分法進行時程積分。啟用集中質(zhì)量矩陣以提高慣性力計算效率。設(shè)定自動時間步長、強制階躍荷載輸入，并采用PCG迭代求解器以提升求解速度。

（4）慣性力施加與求解循環(huán)

使用ACEL命令在每個時間步中施加地震加速度（X/Y/Z方向），通過循環(huán)控制結(jié)構(gòu)響應(yīng)的積分計算，并以等效慣性力的形式參與系統(tǒng)平衡方程的求解，模擬結(jié)構(gòu)在整個地震作用過程中的動力響應(yīng)。

圖 4 時程分析計算完成

6.3 時程分析結(jié)果后處理

為提取結(jié)構(gòu)在地震作用下的動力響應(yīng)特征，本命令流使用ANSYS的/POST26時程后處理模塊，對結(jié)構(gòu)關(guān)鍵節(jié)點（節(jié)點編號201）在地震時程分析過程中的位移、速度與加速度響應(yīng)進行了提取與計算。

（1）模塊切換與變量預(yù)設(shè)

進入時程分析專用的后處理模塊/POST26，并預(yù)設(shè)了最多20個變量存儲空間。

（2）提取節(jié)點位移數(shù)據(jù)

使用NSOL命令分別提取節(jié)點201在X、Y、Z 三個方向的位移時程（UX、UY、UZ），作為基礎(chǔ)響應(yīng)量。

（3）計算速度響應(yīng)

利用DERIV命令對位移曲線進行一階時間導(dǎo)數(shù)計算，獲得各方向上的速度響應(yīng)（VX、VY、VZ）。

（4）計算加速度響應(yīng)

對速度時程繼續(xù)求導(dǎo)，獲得加速度時程（AX、AY、AZ），用于進一步評估結(jié)構(gòu)的地震響應(yīng)。

7 計算結(jié)果分析

圖5展示了結(jié)構(gòu)頂部在地震作用下三個方向的加速度響應(yīng)時程曲線。

圖 5 結(jié)構(gòu)頂部加速度響應(yīng)

自此，基于ANSYS的工程結(jié)構(gòu)抗震分析全過程結(jié)束，感興趣的小伙伴可以私信聯(lián)系。