版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、.OpenSEES解題一般規(guī)律、技巧總結(jié)單位OpenSEES中是可以用公制單位(N,m )的(而并不是像某些文章中說的 OpenSees默 認(rèn)為英制單位”)。實際上我認(rèn)為OpenSEES中并沒有什么默認(rèn)單位,只要編程者自己保持單位一致就行;這點類似于SAP2000的風(fēng)格。建模順序做事要講究順序,OpenSEES建模亦如是:必須先定義材料才能離散截面(因為離散截面 時要對所劃分的截面指定材料屬性)。與之類似的,必須先定義(離散)截面,才能定義非線性梁柱單元(因為定義非線性梁柱單元時要指定單元截面)。關(guān)于BandSPD求解方式官網(wǎng)關(guān)于BandSPD方程形式的評價:This is a good ch
2、oice for most small size models.并且后面緊跟了一句:The equations have to be numbered so the widely used RCM (Reverse Cuthill-McKee) numberer is used. 可見numberer 類型不是隨便選,而是要根據(jù)方程類型來決定的!.(不過直到作業(yè)做完,我對 numberer, system, test, algorithm, analysis(還包括geomTransf, constraints )等求解控制命令還是一知半解!我覺得要想弄明白這些命令得先回頭好好翻翻有限元和數(shù)值
3、分析的書了 ! 自二)OpenSEES中默認(rèn)的計算精度比較高!“0.1000000000000001W0.1 ” :(自行總結(jié),未找到官方說明)這是一個真實的故事:我曾在程序中自以為是的將一連串相鄰均只有0.1左右的數(shù)的差強(qiáng)行賦值為0.1 ,而沒有采用循環(huán)命令將兩數(shù)作差并將結(jié)果賦給新變量一一其中即有這樣的強(qiáng)行截斷!我以為小數(shù)點后都n位了,即使我?guī)е詈笠矔幌到y(tǒng)截斷,還不如我直接預(yù)處理來得清爽!沒想到這樣做直接導(dǎo)致計算不收斂!真是失之毫厘謬以千里!可見在 OpenSEES中默認(rèn)的計算精度比較高!后來我還在老師給的一份范例程序(Silvia Mazzoni & Frank McKenna, 2
4、006)中發(fā)現(xiàn)了set Ubig 1.e10;# a really large number這么一段:set Usmall expr 1/$Ubig; # a really small number.可見系統(tǒng)并未認(rèn)為 Usmall=0 !再一次印證了這一點!劃分纖維截面時角點坐標(biāo)輸入的門道劃分纖維截面時角點坐標(biāo)輸入非常有講究!為了說的直白,我把要點放到下面這張圖中了:Avi(yi, zi)2(也,Z2)Opens EES中矩形橫截面的坐標(biāo)軸以及角點坐蝌述規(guī)則例:用鋼將上圖所示截面分為9 ( =3*3 )塊,則命令為:section Fiber $TAGsec ; #TAGse( 為纖維截面代號
5、patch quad STAGsteel 3 3 yl zl y2 z2 y3 z3 y4 z4; #TAG$tcel為鋼材料代號 這里弘z坐標(biāo)的頃序萬萬不能亂寫!否則運行結(jié)具會出造!.數(shù)據(jù)文件處理OpenSEES運行中是可以生成并讀寫txt文檔的!注意我說是“讀寫”哦!(生成 txt文檔的好處是方便運行完后雙擊生成的數(shù)據(jù)文件讀取數(shù)據(jù),你懂的。)Tcl編程語法(1)命令流中不能出現(xiàn)中文標(biāo)點(這一點和C語言編程類似)!(否則運行時 DOS窗會停住,給出警告,表明不識別命令流中的中文標(biāo)點。)if-else 語句中if和后面緊跟的大括號之間、else和前后大括號之間都要空一格。如:if $a0 se
6、t b 3 else set b -3; # 如果a大于0,則令b等于3,否則等于-3。類似的,相鄰的兩個大括號(一個反大括號和一個正大括號)之間也必須有一個空格。Tcl語言對命令名、變量名區(qū)分大小寫。一行一般只寫一條語句;若想寫多條,則各語句間應(yīng)用分號隔開一一當(dāng)一行只有一條語句時,句末分號可有可無。.同時還有一種特殊情況,就是當(dāng)在一條命令后(同一行中)加注釋時,該命令末尾必須有分號以告知編譯器該命令結(jié)束,否則編譯器會認(rèn)為該注釋也是前面的命令的一部分,導(dǎo)致 編譯出錯。(5)引用變量時,要在變量名前加上 $ (美元符號)!這個步驟非?,嵥?,不如 C語言編 程簡潔。大家就忍著點吧!運行方法目前我知
7、道至少有兩種:方法一:直接運行 OpenSEES,在OpenSees 提示符后輸入source *.tcl (*.tcl ”是提前編寫好的命令流),然后回車。優(yōu)點:個人認(rèn)為沒有二;缺點:命令流編輯時易犯格式錯誤,每次運行都需運行OpenSEES ,再在那個黑框里敲命令流,各種不方便?。ㄆ鋵嵾@個方法只是說說而已,實際我從來沒用過。)有的童鞋用Ultra Editor 之類的通用文本編輯軟件寫程序再導(dǎo)入OpenSEES運行,我沒試過,估計應(yīng)該沒有下面說的第二種方法好。方法二:借助第三方專用編譯環(huán)境。我目前一直用的是 Tcl Editor。它的優(yōu)點有:可以用不同顏色區(qū)分不同功能語句;還有 “查找”、
8、“加注釋”、“取消注釋”等基本常用功能;更好的是菜單欄有個按鈕直接與.OpenSEES關(guān)聯(lián),點擊就可以調(diào)用OpenSEES求解,比較方便(當(dāng)然,比起 VisualStidio之類的還是差遠(yuǎn)了 !可惜誰叫OpenSEES是這么小眾呢?)。但這個軟件有個非常大的缺點一一編程者無法獲知當(dāng)前所編輯文檔的路徑!如果你需要在編程時參考其他文件夾下同樣名字的文件(這樣的情況在我這次做作業(yè)時經(jīng)常發(fā)生,因為我編輯的文檔和模板文檔文件名相同),把兩個文件都用Tcl Editor打開后,你稍一不留神,就會忘了你當(dāng)前編輯的文檔到底是哪里的文件!那時可真是麻煩!所以我總結(jié),這就要求編程者:1、每次在Tcl Editor
9、里打開文件時,不要一看文件名對就急忙打開一一還要看看這個文件是不是在正確的文件夾下面;2、編程時最好一氣呵成;長時間休息時最好把Tcl Editor關(guān)掉。說句題外話,我認(rèn)為一個好的專業(yè)軟件應(yīng)該做到讓用戶大部分時間只需要考慮專業(yè)相關(guān)的東西,而不必操心其他。所以我想,如果以后自己需要經(jīng)常用OpenSEES的話,看能不能用其他的編譯器,不用這個Tcl Editor 了。計算不收斂,怎么辦?可以考慮如下幾點:1、材料本構(gòu)設(shè)定是否正確?材料本構(gòu)參數(shù)是否合理?而且有時steel01比steel02、concrete01 比concrete02 好收斂,如果可以的話不妨改改材料模型。.2-是否極限位移給的過
10、大,柱子已經(jīng)破壞?一一把極限位移改小一點試試。(這是針對我這次作業(yè)而言)3-是否收斂容差太苛刻?一一把容差改大一點試試。(盡管這也許并不是真正解決問題的辦法?。┱{(diào)試程序的技巧:控制變量法控制變量法大家應(yīng)該很熟悉了。調(diào)試程序中我的經(jīng)驗是:一次改動的參數(shù)不要太多,改動的是哪些變量自己要記得。最好一次只改一個變量。然后根據(jù)運行結(jié)果隨所作改動變化的規(guī)律,及時將變量修改到合適的值。這樣做看起來慢,其實我覺得是步步為營,效率比較高。(這些其實應(yīng)該是編程的通用技巧,經(jīng)常編程的朋友們應(yīng)該都有體會。)建議在程序中多用公式一個比較好的編程習(xí)慣是,程序里能輸公式的地方就輸公式,讓用戶只需要給定幾個基本參數(shù)。不要自己
11、事先把中間量在草稿紙上算出來然后輸?shù)匠绦蚶镆灰贿@樣不僅程序通用性不高,而且計算精度也沒有電腦算的高(我前面已經(jīng)說過, OpenSEES中默認(rèn)計算精度是非常高的?。?,真是“吃力不討好”!.關(guān)于wipe命令后面的分號上文說過,如果一條命令后面(同一行中)沒有其他命令或注釋,那么該命令末尾既可以帶分號,也可不帶。但是我發(fā)現(xiàn)對于 wipe命令則不然:因為在 Tcl Editor中可以發(fā)現(xiàn),如果wipe末尾帶了分號(該行再無其他字符),wipe這個單詞是黑色的;但若去掉該分號,wipe就變成了綠色一一從顏色變化上彳#測,莫非加了分號導(dǎo)致wipe命令不被識別?再考慮到下文將提到的“ OpenSEES似乎存
12、在計算不穩(wěn)定現(xiàn)象”與 wipe命令間千絲萬縷 的曖昧關(guān)系,對于 wipe這個“黑匣子”我還是保守處理一一去掉末尾的分號吧!事實證 明去掉分號后貌似有幾次曾經(jīng)不收斂的計算神奇般的收斂了!疑問:OpenSEES計算結(jié)果似乎不太穩(wěn)定?具體表現(xiàn)就是,你現(xiàn)在運行某個命令流算題,計算收斂,得到解了;然后你根本就不改程序,甚至連Tcl Editor都關(guān)了,更甚至連電腦都關(guān)了,等會再重新運行這個命令流,有可能不收斂!我遇到過很多次這個現(xiàn)象,還有同學(xué)出現(xiàn)剛開始算不通過,后來啥也沒改,重新運行一一竟然順利通過了!*口難道是內(nèi)存調(diào)用錯誤?可程序開頭不是由 wipe命令嗎?或許這個wipe根本就不像官網(wǎng)上介紹的那樣每
13、次運行都能徹底destory內(nèi)存中所有之前建的模型、對象?(當(dāng)然也不排除我們在兩次結(jié)果不同的計算中間無意改動了程序某個部分而自己又忘了一一畢竟調(diào)試程序很復(fù)雜,控制變量法調(diào)n個參數(shù),特別像我們這種初學(xué)者,一調(diào)就幾個小時,最后是頭暈眼花,腰酸背疼所以偶爾忘記自己對程序細(xì)微的改r t 口_八,動也是可能的).關(guān)于element recorder里軸力和剪力的正方向里的解釋是:element recorder里記錄的軸力和剪力的正方向是怎樣的?User ManualThese forces correspond to the global coordinate axes orientation.我的理
14、解圖解如下圖所示(圖中剪力和軸力都為正)其他兩個英制單位換算:.1、kipone thousand pounds force ,千磅力,約相當(dāng)于 4,445.205226 NkN。2、重力加速度 g=9.8 m/sA2 =385.8 inch/sA2。M.45在學(xué)習(xí)OpenSEES中常見的英文縮寫、專業(yè)英語術(shù)語:OpenSEES : Open System for Earthquake Engineering SimulationNSF : National Science FoundationPEER: Pacific Earthquake Engineering Research Cent
15、er(為什么不縮寫成PEERC呢?NEES: Network for Earthquake Engineering SimulationPBEE: Performance-Based Earthquake Engineeringportal frame :門式剛架,龍門架( planar portal frame:平面門式剛架)ndm : number of dimensions per nodendf : number of degrees of freedom per nodetranslation : uniform motion of a body in a straight line
16、剛體位移,線位移normal : 【數(shù)】 法線 rotation about the planes normal繞平面法線的轉(zhuǎn)動prompt :提示。(個人理解:指DOS窗中一閃一閃的光標(biāo),學(xué)名“命令提示符”。).a-priori :先驗的。(be generated a-priori ,個人理解:即“事先編輯好的“)geometry :幾何尺寸element :單元component :(地面運動的)分量drift (angle):位移角uniaxial : of or relating to only one axis , 單軸的time series : a set of data co
17、llected sequentially usually at fixed intervals of time時間序列argument : one of the independent variables upon whose value that of a functiondepends ,自變量,參數(shù)nodal load :節(jié)點荷載(這個詞組讀起來有點繞口A_A )flag :標(biāo)志變量(學(xué)過 C語言編程的同學(xué)應(yīng)該知道?。ヾiscretization :離散化offset :偏移(量)iteration :迭代(法)SOE : system of linear equations ,線性方程
18、系統(tǒng):線性搜索路線牛頓算法(這個翻譯可能不準(zhǔn)確)Newton with Line Search Algorithm.BandSPD : Banded Symmetric Positive Definitemap :映射關(guān)于OpenSEES的學(xué)習(xí)資料:來自官方:1、OpenSees User (里面有詳細(xì)的幫助文檔)2、Example Mannul (官網(wǎng)提供的幫助手冊)3、OpenSees官網(wǎng)提供的Tcl編程教程4、網(wǎng)頁版用戶手冊(09年版)5、OpenSeesNavigator(OpenSeesNavigator is a matlab interface for OpenSees. It
19、allows users to quickly create models, perform analysis, and look at the results. It runs on windows machines.)6、Basic Examples Manual(The objective of this primer is to provide new users of OpenSees (Open System for Earthquake Engineering Simulation) familiar structural engineering examples as a co
20、nvenient method for learning how to use the software.原創(chuàng)教程及其他:1、車t教程OpenSEES超簡單易懂的入門第一課2、OpenSEES最好的介紹(這個講了一些數(shù)值模擬的本質(zhì),比較理論化。)3、結(jié)構(gòu)藝術(shù)家劉金成的博客玩轉(zhuǎn) opensees e.g.1玩轉(zhuǎn) opensees- e.g.2 暨 sap2000 vs opensees (IV)4、(百度文庫)L形鋼管混凝土柱抗震性能非線性有限元分析5、(豆?。㎡penSees 使用與教程(一個豆單,里面有 5篇文章。其中:OpenSees 這篇是節(jié)選自一篇重慶大學(xué)碩士論文中關(guān)于OpenSees
21、的部分,個人感覺寫得很好。)6、(中華鋼結(jié)構(gòu)論壇)用過 OpenSees的前輩請進(jìn)來幫幫我?。槺阏f一下,這個論壇話題討論質(zhì)量很高,經(jīng)常有高手出沒! 、*)7、廈門大學(xué)建筑與土木工程學(xué)院古泉副教授創(chuàng)辦的OpenSees學(xué)習(xí)網(wǎng)站(以前印象中廈大只有文科,瀏覽了這個網(wǎng)站后我再不敢這么想了)8、推薦一個 QQ群:61667191(Opensees與地震工程2),感興趣的話可以加入;里面都是OpenSEES高手與愛好者,討論的基本都是學(xué)術(shù)方面的事情。.9、dinoChen 的博客(這是陳學(xué)偉博士的個人網(wǎng)站,他還有個QQ空間,兩個網(wǎng)站上都有非常豐富的優(yōu)質(zhì)原創(chuàng)資源!陳博士實乃青年才??!佩服!上面提到的 Q
22、Q群也是他創(chuàng)建的?。?0、TCL/TK編程實踐第四版 icelamp的專欄一一 CSDN博客11、Tcl簡介一一維基百科最后,謹(jǐn)摘錄兩段高手的OpenSEES的學(xué)習(xí)心得,與君共勉。感謝原作者!ocean2000:“我的當(dāng)初畢業(yè)論文也是用os做試驗仿真的,臺灣有一篇砌體結(jié)構(gòu)實驗用os來仿真的,可以google 之,os有一些可以模擬的 2D單元。這個程序只要過一遍mannual ,使用一點不難,而且資料也很多了,它的論壇不錯,問題的討論很積極。要加新單元和新材料也不難,接口都給大家提供好了。其實這個軟件的一大優(yōu)點是TCL/TK很容易上手,對于計算過程可以加入自己的判斷和控制條件,所以可以做出很好
23、的滯回曲線,還可以動態(tài)顯示曲線發(fā)展過程。非常難得的是在工作中居然看到單位的軟件有與opensees的數(shù)據(jù)轉(zhuǎn)化接口?!眃inochen1983:.“學(xué)習(xí)OPENSEES要求一定的有限元知識及非線性理論,最好會編程,因為建模需要用編程 的思想去簡化重復(fù)輸入,本人覺得OPENSEES值得大家女?好學(xué)習(xí),我學(xué)了半年,收獲很多。”(兩段話均來自:用過OpenSees的前輩請進(jìn)來幫幫我!一一中華鋼結(jié)構(gòu)論壇)(2012年4月19日11:26:49 更新)原文中所謂的“關(guān)于這次作業(yè)具體的心得”的博文鏈接已更新;本打算待該博文發(fā)表后及時回來更新的,結(jié)果忘了。抱歉!來源:顏冬啟全文鏈接10【轉(zhuǎn)載】OpenSEES
24、自學(xué)筆記(一) “博主按”:本文是我第一次用OpenSEES做仿真分析作業(yè)(基于 OpenSEES的方鋼管混凝土柱抗震性能分析)過程中點滴記錄的自學(xué)筆記,發(fā)表出來既是和各位(尤其是Open. HYPERLINK http:/163.fm/9qMUBJZ http:/163.fm/9qMUBJZ.全文鏈接10【轉(zhuǎn)載】OpenSEES自學(xué)筆記(一)“博主按”:本文是我第一次用 OpenSEES做仿真分析作業(yè)(基于 OpenSEES的方鋼管混凝土柱抗震性能分析)過程中點滴記錄的自學(xué)筆記,發(fā)表出來既是和各位(尤其是OpenSEES初學(xué)者)交流,同時也算作個自我小結(jié)以備日后查閱。盡管我力求完美,但這些習(xí)
25、得中仍然極有可能存在錯誤!請注意甄別!同時也衷心希望各位高手不吝賜教!另外,由于時間倉促,本人又是初學(xué)OpenSEES ,所以文章內(nèi)容上比較零散,見諒!.初識 OpenSEES我是在鋼筋混凝土結(jié)構(gòu)非線性分析這門課上第一次聽說這個軟件的。老師說(均為個人理解,可能不是老師原話)這個軟件能夠用纖維單元做有限元分析,在模擬大型結(jié)構(gòu)上 比ANSYS、SAP等利用實體單元的有限元程序有優(yōu)勢;經(jīng)常用于抗震分析科研中;不是一個“設(shè)計型”軟件(如 SAP、PKPM、橋博等);還要求我們用它做兩個大作業(yè)。在 Silvia Mazzoni, Frank McKenna, Michael H. Scott, Gre
26、gory L. Fenves等人編寫的OpenSEES 的 Users Manual (v2.0)開篇,是這樣回答 What is OpenSEES? 這個問題的:An object-oriented software framework for simulation applications in earthquake engineering using finite element methods. OpenSees is not a code.A communication mechanism within PEER for exchanging and building upon re
27、search accomplishments.As open-source software, it has the potential for a community code for earthquake engineering.好吧,既然是專業(yè)軟件,那咱就在接下來的使用中逐漸熟悉吧!軟件下載與安裝OpenSEES和Tcl的下載頁面鏈接在 OpenSEES官網(wǎng)首頁左側(cè)的欄目里,點擊aDownload ”即可進(jìn)入下載頁面(下載之前需要注冊(新用戶)或填寫電郵(已注冊用 戶)。我的電腦用的是 64位的 Win7 (2012年5月23日11:44:48更新:Win7 旗艦版 SP1) , Act
28、iveTcl.295402-win32-ix86-threaded”安裝成功。安裝路徑按照官網(wǎng)說明手動作了修改(即將安裝路徑由默認(rèn)的C:Tcl改為C:Program FilesTcl )。注意:必須以管理員身份運行ActiveTcl安裝程序(在安裝包上右擊,選擇“以管理員身份運行”),否則會安裝失?。penSEES的運行界面如下圖所示。.小技巧:可能由于 Tcl Editor (下文將介紹)是綠色軟件,系統(tǒng)默認(rèn)并沒有把 tcl文件和Tcl Editor關(guān)聯(lián)起來,而且一般也沒有建立tcl文件的其它關(guān)聯(lián)方式,所以如果直接雙擊tcl文件不僅無法打開它,還會彈出錯誤提示對話框(提示無法識別該tcl文
29、件頭幾行)!如果你想實現(xiàn)雙擊tcl文件調(diào)用Tcl Editor進(jìn)行編輯的話,可以自行更改文件關(guān)聯(lián)。具體操作就是在任意一個tcl文件上右擊,選擇打開方式,然后通過“瀏覽”找到下),并勾上始終用該程TclEditor.exe 這個可執(zhí)行文件(一般應(yīng)該是在 序打開此類文件選項,確認(rèn)。.以下大致總結(jié)下自行摸索的OpenSEES 一般編程規(guī)律、技巧。關(guān)于這次作業(yè)具體的心得在這篇博文里:OpenSEES自學(xué)筆記(二)。(兩篇寫一起既亂且長。)OpenSEES解題一般規(guī)律、技巧總結(jié)單位OpenSEES中是可以用公制單位(N,m )的(而并不是像某些文章中說的 OpenSees默 認(rèn)為英制單位”)。實際上我認(rèn)
30、為 OpenSEES中并沒有什么默認(rèn)單位,只要編程者自己保 持單位一致就行;這點類似于SAP2000的風(fēng)格。建模順序做事要講究順序,OpenSEES建模亦如是:必須先定義材料才能離散截面(因為離散截面 時要對所劃分的截面指定材料屬性)。與之類似的,必須先定義(離散)截面,才能定義非線性梁柱單元(因為定義非線性梁柱 單元時要指定單元截面)。關(guān)于BandSPD求解方式官網(wǎng)關(guān)于BandSPD方程形式的評價:This is a good choice for most small size models.并且后面緊跟了一句:The equations have to be numbered so th
31、e widely used RCM (Reverse Cuthill-McKee) numberer is used.可見numberer 類型不是隨便選,而是要根據(jù)方程類型來決定的!(不過直到作業(yè)做完,我對 numberer, system, test, algorithm, analysis(還包括geomTransf, constraints )等求解控制命令還是一知半解!我覺得要想弄明白這些命令得先回頭好好翻翻有限元和數(shù)值分析的書了 !起”)OpenSEES中默認(rèn)的計算精度比較高!“0.1000000000000001W0.1 ” :(自行總結(jié),未找到官方說明)這是一個真實的故事:我曾
32、在程序中自以為是的將一連串相鄰均只有0.1左右的數(shù)的差強(qiáng)行賦值為0.1 ,而沒有采用循環(huán)命令將兩數(shù)作差并將結(jié)果賦給新變量一一其中即有這樣的強(qiáng)行截斷!我以為小數(shù)點后都n位了,即使我?guī)е詈笠矔幌到y(tǒng)截斷,還不如我直接預(yù)處理來得清爽!沒想到這樣做直接導(dǎo)致計算不收斂!真是失之毫厘謬以千里!可見在 OpenSEES中默認(rèn)的計算精度比較高!后來我還在老師給的一份范例程序(Silvia Mazzoni & Frank McKenna, 2006)中發(fā)現(xiàn)了這么一段:.set Ubig 1.e10;# a really large numberset Usmall expr 1/$Ubig; # a rea
33、lly small number可見系統(tǒng)并未認(rèn)為 Usmall=0 !再一次印證了這一點!劃分纖維截面時角點坐標(biāo)輸入的門道劃分纖維截面時角點坐標(biāo)輸入非常有講究!為了說的直白,我把要點放到下面這張圖中 了:.Ar4(y Z4)3 仇 Z3)Zl)2(y2fz2)Opens EES中矩形橫截面的坐標(biāo)軸以及角點坐標(biāo)描述規(guī)則例:用鋼將上圖所示截面分為% =3*3 )塊,則命令為:section Fiber STAGsec ;#TAGse0 set b 3 else set b -3; # 如果a大于0,則令b等于3,否則等于-3。類似的,相鄰的兩個大括號(一個反大括號和一個正大括號)之間也必須有一個空
34、格。Tcl語言對命令名、變量名區(qū)分大小寫。一行一般只寫一條語句;若想寫多條,則各語句間應(yīng)用分號隔開一一當(dāng)一行只有一條語句時,句末分號可有可無。同時還有一種特殊情況,就是當(dāng)在一條命令后(同一行中)加注釋時,該命令末尾必須有分號以告知編譯器該命令結(jié)束,否則編譯器會認(rèn)為該注釋也是前面的命令的一部分,導(dǎo)致 編譯出錯。.(5)引用變量時,要在變量名前加上$ (美元符號)!這個步驟非常瑣碎,不如 C語言編程簡潔。大家就忍著點吧!運行方法目前我知道至少有兩種:方法一:直接運行 OpenSEES,在OpenSees 提示符后輸入source *.tcl (*.tcl ”是提前編寫好的命令流),然后回車。優(yōu)點:
35、個人認(rèn)為沒有三;缺點:命令流編輯時易犯格式錯誤,每次運行都需運行OpenSEES ,再在那個黑框里敲命令流,各種不方便?。ㄆ鋵嵾@個方法只是說說而已,實際我從來沒用過。)有的童鞋用Ultra Editor 之類的通用文本編輯軟件寫程序再導(dǎo)入OpenSEES運行,我沒試過,估計應(yīng)該沒有下面說的第二種方法好。方法二:借助第三方專用編譯環(huán)境。我目前一直用的是 Tcl Editor。它的優(yōu)點有:可以用不同顏色區(qū)分不同功能語句;還有“查找”、“加注釋”、“取消注釋”等基本常用功能;更好的是菜單欄有個按鈕直接與.OpenSEES關(guān)聯(lián),點擊就可以調(diào)用OpenSEES求解,比較方便(當(dāng)然,比起 VisualSt
36、idio之類的還是差遠(yuǎn)了 !可惜誰叫OpenSEES是這么小眾呢?)。但這個軟件有個非常大的缺點一一編程者無法獲知當(dāng)前所編輯文檔的路徑!如果你需要在編程時參考其他文件夾下同樣名字的文件(這樣的情況在我這次做作業(yè)時經(jīng)常發(fā)生,因為我編輯的文檔和模板文檔文件名相同),把兩個文件都用Tcl Editor打開后,你稍一不留神,就會忘了你當(dāng)前編輯的文檔到底是哪里的文件!那時可真是麻煩!所以我總結(jié),這就要求編程者:1、每次在Tcl Editor里打開文件時,不要一看文件名對就急忙打開一一還要看看這個文件是不是在正確的文件夾下面;2、編程時最好一氣呵成;長時間休息時最好把Tcl Editor關(guān)掉。說句題外話,
37、我認(rèn)為一個好的專業(yè)軟件應(yīng)該做到讓用戶大部分時間只需要考慮專業(yè)相關(guān)的東西,而不必操心其他。所以我想,如果以后自己需要經(jīng)常用OpenSEES的話,看能不能用其他的編譯器,不用這個Tcl Editor 了。計算不收斂,怎么辦?可以考慮如下幾點:1、材料本構(gòu)設(shè)定是否正確?材料本構(gòu)參數(shù)是否合理?而且有時steel01比steel02、concrete01 比concrete02 好收斂,如果可以的話不妨改改材料模型。.2-是否極限位移給的過大,柱子已經(jīng)破壞?一一把極限位移改小一點試試。(這是針對我這次作業(yè)而言)3-是否收斂容差太苛刻?一一把容差改大一點試試。(盡管這也許并不是真正解決問題的辦法?。┱{(diào)試程
38、序的技巧:控制變量法控制變量法大家應(yīng)該很熟悉了。調(diào)試程序中我的經(jīng)驗是:一次改動的參數(shù)不要太多,改動的是哪些變量自己要記得。最好一次只改一個變量。然后根據(jù)運行結(jié)果隨所作改動變化的規(guī)律,及時將變量修改到合適的值。這樣做看起來慢,其實我覺得是步步為營,效率比較高。(這些其實應(yīng)該是編程的通用技巧,經(jīng)常編程的朋友們應(yīng)該都有體會。)建議在程序中多用公式一個比較好的編程習(xí)慣是,程序里能輸公式的地方就輸公式,讓用戶只需要給定幾個基本參數(shù)。不要自己事先把中間量在草稿紙上算出來然后輸?shù)匠绦蚶镆灰贿@樣不僅程序通用性不高,而且計算精度也沒有電腦算的高(我前面已經(jīng)說過, OpenSEES中默認(rèn)計算精度是非常高的?。?,真
39、是“吃力不討好”!.關(guān)于wipe命令后面的分號上文說過,如果一條命令后面(同一行中)沒有其他命令或注釋,那么該命令末尾既可以帶分號,也可不帶。但是我發(fā)現(xiàn)對于 wipe命令則不然:因為在 Tcl Editor中可以發(fā)現(xiàn),如果wipe末尾帶了分號(該行再無其他字符),wipe這個單詞是黑色的;但若去掉該分號,wipe就變成了綠色一一從顏色變化上彳#測,莫非加了分號導(dǎo)致wipe命令不被識別?再考慮到下文將提到的“ OpenSEES似乎存在計算不穩(wěn)定現(xiàn)象”與 wipe命令間千絲萬縷 的曖昧關(guān)系,對于 wipe這個“黑匣子”我還是保守處理一一去掉末尾的分號吧!事實證 明去掉分號后貌似有幾次曾經(jīng)不收斂的計
40、算神奇般的收斂了!疑問:OpenSEES計算結(jié)果似乎不太穩(wěn)定?具體表現(xiàn)就是,你現(xiàn)在運行某個命令流算題,計算收斂,得到解了;然后你根本就不改程序,甚至連Tcl Editor都關(guān)了,更甚至連電腦都關(guān)了,等會再重新運行這個命令流,有可能不收斂!我遇到過很多次這個現(xiàn)象,還有同學(xué)出現(xiàn)剛開始算不通過,后來啥也沒改,重新運行一一竟然順利通過了!*口難道是內(nèi)存調(diào)用錯誤?可程序開頭不是由 wipe命令嗎?或許這個wipe根本就不像官網(wǎng)上介紹的那樣每次運行都能徹底destory內(nèi)存中所有之前建的模型、對象?(當(dāng)然也不排除我們在兩次結(jié)果不同的計算中間無意改動了程序某個部分而自己又忘了一一畢竟調(diào)試程序很復(fù)雜,控制變量
41、法調(diào)n個參數(shù),特別像我們這種初學(xué)者,一調(diào)就幾個小時,最后是頭暈眼花,腰酸背疼所以偶爾忘記自己對程序細(xì)微的改r t 口_八,動也是可能的).關(guān)于element recorder里軸力和剪力的正方向里的解釋是:element recorder里記錄的軸力和剪力的正方向是怎樣的?User ManualThese forces correspond to the global coordinate axes orientation.我的理解圖解如下圖所示(圖中剪力和軸力都為正)其他兩個英制單位換算:.1、kipone thousand pounds force ,千磅力,約相當(dāng)于 4,445.2052
42、26 NkN。2、重力加速度 g=9.8 m/sA2 =385.8 inch/sA2。M.45在學(xué)習(xí)OpenSEES中常見的英文縮寫、專業(yè)英語術(shù)語:OpenSEES : Open System for Earthquake Engineering SimulationNSF : National Science FoundationPEER: Pacific Earthquake Engineering Research Center(為什么不縮寫成PEERC呢?NEES: Network for Earthquake Engineering SimulationPBEE: Performan
43、ce-Based Earthquake Engineeringportal frame :門式剛架,龍門架( planar portal frame:平面門式剛架)ndm : number of dimensions per nodendf : number of degrees of freedom per nodetranslation : uniform motion of a body in a straight line剛體位移,線位移normal : 【數(shù)】 法線 rotation about the planes normal繞平面法線的轉(zhuǎn)動prompt :提示。(個人理解:指D
44、OS窗中一閃一閃的光標(biāo),學(xué)名“命令提示符”。).a-priori :先驗的。(be generated a-priori ,個人理解:即“事先編輯好的“)geometry :幾何尺寸element :單元component :(地面運動的)分量drift (angle):位移角uniaxial : of or relating to only one axis , 單軸的time series : a set of data collected sequentially usually at fixed intervals of time時間序列argument : one of the in
45、dependent variables upon whose value that of a functiondepends ,自變量,參數(shù)nodal load :節(jié)點荷載(這個詞組讀起來有點繞口A_A )flag :標(biāo)志變量(學(xué)過 C語言編程的同學(xué)應(yīng)該知道!)discretization :離散化offset :偏移(量)iteration :迭代(法)SOE : system of linear equations ,線性方程系統(tǒng):線性搜索路線牛頓算法(這個翻譯可能不準(zhǔn)確)Newton with Line Search Algorithm.BandSPD : Banded Symmetric Positive Definitemap :映射關(guān)于OpenSEES的學(xué)習(xí)資料:來自官方:1、OpenSees User (里面有詳細(xì)的幫助文檔)2、Example Mannul (官網(wǎng)提供的幫助手冊)3、OpenSees官網(wǎng)提供的Tcl編程教程4、網(wǎng)頁版用戶手冊(09年版)5、OpenSeesNavigator(OpenSeesNavigator is a
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 二零二五年版鋁合金散熱器研發(fā)與生產(chǎn)合同樣本2篇
- 引進(jìn)創(chuàng)新創(chuàng)業(yè)團(tuán)隊2025年度協(xié)議3篇
- 2025年度個人商鋪租賃合同裝修工程監(jiān)督合同4篇
- 二零二五年度地質(zhì)災(zāi)害風(fēng)險評估與防治合同3篇
- 2025年度個人閑置土地承包權(quán)置換合同范本4篇
- 2025年度大型活動策劃與執(zhí)行服務(wù)合同模板4篇
- 二零二五年度建筑廢棄物處理項目擔(dān)保合同范本3篇
- 二婚家庭專用2024離婚協(xié)議范本樣本版B版
- 2025年度趕工環(huán)保設(shè)備趕工銷售合作協(xié)議書
- 二零二五年度土地租賃權(quán)租賃與轉(zhuǎn)租合同
- 高考語文復(fù)習(xí)【知識精研】《千里江山圖》高考真題說題課件
- 河北省承德市2023-2024學(xué)年高一上學(xué)期期末物理試卷(含答案)
- 高中物理斜面模型大全(80個)
- 012主要研究者(PI)職責(zé)藥物臨床試驗機(jī)構(gòu)GCP SOP
- 農(nóng)耕研學(xué)活動方案種小麥
- 2024年佛山市勞動合同條例
- 污水管網(wǎng)規(guī)劃建設(shè)方案
- 城鎮(zhèn)智慧排水系統(tǒng)技術(shù)標(biāo)準(zhǔn)
- 采購管理制度及流程采購管理制度及流程
- 五年級美術(shù)下冊第9課《寫意蔬果》-優(yōu)秀課件4人教版
- 節(jié)能降耗課件
評論
0/150
提交評論