一、寫在文前

隨著我國經(jīng)濟(jì)快速發(fā)展，全國各地的建筑迅速的拔地而起，魔都的上海環(huán)球金融中心是其著名標(biāo)志性建筑，以樓高492米直穿云霄。面對復(fù)雜的建筑的設(shè)計(jì)，它是由無數(shù)的梁、柱、板、墻最基本的單元組成，其中最基本的結(jié)構(gòu)是一榀門式框架，在我們對結(jié)構(gòu)仿真中，軟件其實(shí)只是工具，真正對我們設(shè)計(jì)負(fù)責(zé)的還是對力學(xué)、結(jié)構(gòu)變形機(jī)理的認(rèn)識，這是我們最應(yīng)該掌握的，這樣才能明白軟件是怎么算的，對我們的仿真結(jié)果才可以加以正確的判斷。

接 下來將以最簡單的一榀門式框架為例子，用不同軟件說明對門式框架模擬，并輔以手算過程。

二、 案例詳述—靜力計(jì)算

這是一個(gè)一榀門式框架，如圖1，左右兩側(cè)為柱，截面為A-A，寬高為均為5 feet (英尺)，柱高36 feet；上方為梁截面為B-B，寬高分別為5 feet，8 feet，梁長42 feet；梁受均布荷載，荷載總重為4000 kip (千磅力)；OpenSees軟件的建模示意圖，如圖2，右圖為桿件局部坐標(biāo)系示意圖；

圖1 一榀門式框架

圖2 OpenSees軟件的建模示意圖

圖3給出了一榀門式框架部分的命令，具體代碼在附錄中；圖4給出了全局坐標(biāo)系下，單元起始節(jié)點(diǎn)X向、Y向、MZ向的力和彎矩；圖5給出了全局坐標(biāo)系下，單元起始節(jié)點(diǎn)X向、Y向、MZ向的彎矩和轉(zhuǎn)角，具體見附錄；圖5給出了OpenSees算出來的彎矩圖結(jié)果和SAP2000給的彎矩圖，從彎矩圖上看，OpenSees和SAP2000計(jì)算結(jié)果一致。

圖3 部分OpenSees命令

在Sap2000中進(jìn)行建模分析，模型及荷載如圖所示，與OpenSees一樣，僅開啟3自由度(x/z/Roty)進(jìn)行分析計(jì)算。由于在OpenSees中，本例子采用的是歐拉伯努利梁設(shè)置，因此對Sap2000中的剪切截面屬性進(jìn)行調(diào)整，如下圖所示：

圖4 Sap2000模型

對模型進(jìn)行求得分析得到如下圖所示：

圖5 OpenSees中門式框架彎矩圖與SAP2000門式框架計(jì)算結(jié)果

下面我們采用位移法手算一下框架的彎矩，由于手算，我們假設(shè)框架軸向不可被壓縮，所以該框架是一個(gè)三自由度框架，該框架又是對稱結(jié)構(gòu)，可以簡化圖右的單自由度（B點(diǎn)轉(zhuǎn)角）結(jié)構(gòu)，i1 是桿件的彎曲線剛度。

查詢位移法相關(guān)結(jié)構(gòu)的彎矩：

假定無桿件軸向壓縮，因此只有B點(diǎn)轉(zhuǎn)角 ，其中：

框架各桿的截面面積、慣性矩、桿件彎曲線剛度計(jì)算如下，所用基本單位均為inch (英寸)、kip (英鎊)：

框架手算彎矩圖結(jié)果如下圖所示：

各節(jié)點(diǎn)彎矩、B點(diǎn)轉(zhuǎn)角對比表如下：

從對比結(jié)果來看，從圖中對比結(jié)果來看，手算和OpenSees計(jì)算有差別，但是很小，我們沒有考慮軸向壓縮，而OpenSees考慮了軸向壓縮，也可以從側(cè)面說明，在可以忽略軸向變形的影響的結(jié)構(gòu)中，手算的彎矩精度還是很高的。

三、案例詳述—?jiǎng)恿τ?jì)算

首先對模型進(jìn)行模態(tài)求解分析

• 通過OpenSees中讀取求解，可得到結(jié)構(gòu)周期T1=0.587s；

• 通過Sap2000進(jìn)行求解可得，結(jié)構(gòu)周期T1=0.58806s；

• 通過Matlab進(jìn)行求解可得，結(jié)構(gòu)周期T1=0.587s。

（可以看下： 推開土木工程振型求解值蘭索斯法大門 ）

動(dòng)力分析中，我們需要讀入地震波數(shù)據(jù)進(jìn)行時(shí)程分析，地震波（PEER格式）讀入OpenSees命令如下截圖，可在最后附件中下載；

下面是動(dòng)力分析命令相關(guān)記錄，并帶一部分解釋：

①將PEER地震波的時(shí)間間隔dt與加速度acc提取出來，

②設(shè)置accelSeries的內(nèi)容是“Series -dt $dt -filePath BM68elc.acc -factor 1”，時(shí)間間隔$dt，加速度是BM68elc.acc，調(diào)整系數(shù)是1（根據(jù)自己需要的加速度值進(jìn)行調(diào)整），這個(gè)內(nèi)容格式來源于OpenSees讀入時(shí)程數(shù)據(jù)的要求，

③標(biāo)識是2，“1”表示X方向，“$accelSeries”表示②中設(shè)置的加速度值，

④設(shè)置瑞利阻尼，rayleigh阻尼需要兩個(gè)參數(shù)，a0和a1是由下圖計(jì)算出，a0=0時(shí)，a1=(2*ζ)/(sqrt(ω))，a1在OpenSees表達(dá)為“[expr 2*0.02/pow([eigen 1],0.5)]”，下圖中寫出OpenSees中rayleigh阻尼的相關(guān)要求。

a 0 和 a 1 可以分別根據(jù)給定的第i階振型阻尼比 和第j階振型 確定 ：

⑤設(shè)置“總分析步數(shù)”和“單步時(shí)間間隔”，要注意，官網(wǎng)上“analyze [expr 10/0.02] 0.02”命令會(huì)報(bào)錯(cuò)，“總步數(shù)”數(shù)據(jù)格式錯(cuò)誤，會(huì)報(bào)錯(cuò)，改為下圖中即會(huì)正常運(yùn)行。

Sap2000中，對模型定義了線性直接分析方法，也采用Newmark-β法，以及瑞利阻尼（采用剛度阻尼同OpenSees）。

下面給出SAP2000、OpenSees、Matlab對應(yīng)的分析方向的位移時(shí)程對比，吻合良好！

在本分析中，麻雀雖小，五臟俱全，在線彈性體靜力分析與動(dòng)力分析中，很好的表達(dá)的平面計(jì)算的機(jī)理，靜力分析下，采用位移法手算與矩陣位移法電算對比；動(dòng)力分析中在考慮慣性力下，考慮的平衡是形式上的平衡是瞬時(shí)的平衡，很好的闡述了概念為先，機(jī)理為本的本質(zhì)。

往期精彩

#性能分析

【JY】基于性能的抗震設(shè)計(jì)淺析（一）

【JY】基于性能的抗震設(shè)計(jì)淺析（二）

【JY】淺析消能附加阻尼比

【JY】近斷層結(jié)構(gòu)設(shè)計(jì)策略分析與討論

【JY】淺析各動(dòng)力求解算法及其算法數(shù)值阻尼(人工阻尼)

理念

【JY|體系】結(jié)構(gòu)概念設(shè)計(jì)之(結(jié)構(gòu)體系概念)

【JY|理念】結(jié)構(gòu)概念設(shè)計(jì)之(設(shè)計(jì)理念進(jìn)展)

【JY】有限單元分析的常見問題及單元選擇

【JY】結(jié)構(gòu)動(dòng)力學(xué)之顯隱式

【JY】淺談結(jié)構(gòu)設(shè)計(jì)

【JY】淺談混凝土損傷模型及Abaqus中CDP的應(yīng)用

#概念機(jī)理

【JY】推開土木工程振型求解之蘭索斯法(Lanczos法)的大門

【JY】基于Ramberg-Osgood本構(gòu)模型的雙線性計(jì)算分析

【JY】結(jié)構(gòu)動(dòng)力學(xué)初步-單質(zhì)點(diǎn)結(jié)構(gòu)的瞬態(tài)動(dòng)力學(xué)分析

【JY】從一根懸臂梁說起

【JY】反應(yīng)譜的詳解與介紹

【JY】結(jié)構(gòu)瑞利阻尼與經(jīng)濟(jì)訂貨模型

【JY】主成分分析與振型分解

【JY】淺談結(jié)構(gòu)多點(diǎn)激勵(lì)之概念機(jī)理（上）

【JY】淺談結(jié)構(gòu)多點(diǎn)激勵(lì)之分析方法（下）

【JY】板殼單元的分析詳解

【JY】橡膠支座的簡述和其力學(xué)性能計(jì)算

【JY】振型求解之子空間迭代

【JY】橡膠支座精細(xì)化模擬與有限元分析注意要點(diǎn)

#軟件討論

【JY】復(fù)合材料分析利器—內(nèi)聚力單元

【JY】SDOF計(jì)算教學(xué)軟件開發(fā)應(yīng)用分享

【JY】Abaqus案例—天然橡膠隔震支座豎(軸)向力學(xué)性能

【JY】Abaqus6.14-4如何關(guān)聯(lián)fortran？

【JY】如何利用python來編寫GUI？

【JY】如何解決MATLAB GUI編程軟件移植運(yùn)行問題？

【JY】淺談結(jié)構(gòu)分析與設(shè)計(jì)軟件

【JY|STR】求解器之三維結(jié)構(gòu)振型分析

【JY】SignalData軟件開發(fā)應(yīng)用分享

【JY】基于Matlab的雙線性滯回代碼編寫教程

#其他

【JY】位移角還是有害位移角？

【JY】如何利用python來編寫GUI？

【JY】今日科普之BIM

~關(guān)注未來更精彩~