1-鋼筋砼柱抗震性能模擬例子

1、柱截面彎矩-曲率關(guān)系的分析與梁相似（柱存在軸力）壓彎構(gòu)件壓彎構(gòu)件荷載荷載- -撓度撓度單調(diào)曲線(xiàn)的非線(xiàn)性全過(guò)程分析單調(diào)曲線(xiàn)的非線(xiàn)性全過(guò)程分析P效應(yīng)主要區(qū)別：pN12nPNmLMPLN二次矩分層法分層法程序結(jié)束計(jì)算開(kāi)始調(diào)用M計(jì)算程序并儲(chǔ)存M關(guān)系假設(shè)初始撓曲線(xiàn)1由M、1、N求水平力P由P、N、1計(jì)算各截面Mj 由各截面Mj根據(jù)M關(guān)系查各截面曲率j由j積分求2處理塑性鉸，調(diào)整各截面曲率j進(jìn)行下降段處理，調(diào)整各截面曲率j將2賦值給1記錄當(dāng)前的P和輸入基本信息=+1與2是否接近yuP是否降到極限承載力的80YesYesYesYesNoNoNoNo壓彎構(gòu)件壓彎構(gòu)件荷載荷載- -撓撓度度單調(diào)曲線(xiàn)非線(xiàn)性單調(diào)曲線(xiàn)

2、非線(xiàn)性全過(guò)程分析框圖全過(guò)程分析框圖下降段的處理塑性鉸的處理卸載加載P反向加載反向卸載壓彎構(gòu)件壓彎構(gòu)件荷載荷載- -撓度撓度滯回曲線(xiàn)模擬的目的滯回曲線(xiàn)模擬的目的試驗(yàn)滯回曲線(xiàn)試驗(yàn)滯回曲線(xiàn)滯回規(guī)則滯回規(guī)則滯回規(guī)則滯回規(guī)則+ +骨架曲線(xiàn)骨架曲線(xiàn)= =恢復(fù)力滯回模型恢復(fù)力滯回模型KKKd=K/=/yKdKdMy-yMyyMClough恢復(fù)力模型恢復(fù)力模型Clough模型是模型是雙線(xiàn)性模型雙線(xiàn)性模型，在這一模型中，構(gòu)件剛度不受混凝土裂縫對(duì)其帶來(lái)的影，在這一模型中，構(gòu)件剛度不受混凝土裂縫對(duì)其帶來(lái)的影響，當(dāng)所加荷載超過(guò)構(gòu)件屈服強(qiáng)度后，在卸載過(guò)程中，發(fā)生剛度降低的現(xiàn)象，這一響，當(dāng)所加荷載超過(guò)構(gòu)件屈服強(qiáng)度后，在卸

3、載過(guò)程中，發(fā)生剛度降低的現(xiàn)象，這一模型構(gòu)成簡(jiǎn)單，便于理解，在初期鋼筋混凝土結(jié)構(gòu)彈塑性分析中，經(jīng)常被使用。模型構(gòu)成簡(jiǎn)單，便于理解，在初期鋼筋混凝土結(jié)構(gòu)彈塑性分析中，經(jīng)常被使用。Ky=My/yKyKyMy-u-yMyyuMMuMu武藤武藤恢復(fù)力模型恢復(fù)力模型武藤模型是武藤模型是三線(xiàn)性模型三線(xiàn)性模型，以此來(lái)考慮構(gòu)件發(fā)生混凝土開(kāi)裂、鋼筋屈服等現(xiàn)象，以此來(lái)考慮構(gòu)件發(fā)生混凝土開(kāi)裂、鋼筋屈服等現(xiàn)象，模型體現(xiàn)的變化規(guī)律比較簡(jiǎn)單。此模型特點(diǎn)是卸載剛度恒定，卸載剛度不隨塑模型體現(xiàn)的變化規(guī)律比較簡(jiǎn)單。此模型特點(diǎn)是卸載剛度恒定，卸載剛度不隨塑性變形增大而發(fā)生變化，所以不能準(zhǔn)確地模擬較大程度的剛度降低現(xiàn)象。性變形增大而

4、發(fā)生變化，所以不能準(zhǔn)確地模擬較大程度的剛度降低現(xiàn)象。Ky=My/yKdKdMy-u-yMyyuMMuMu武田武田恢復(fù)力模型恢復(fù)力模型武田模型是武田模型是三線(xiàn)性模型三線(xiàn)性模型，能夠較為精確地模擬鋼筋混凝土構(gòu)件在反復(fù)荷載作用下的，能夠較為精確地模擬鋼筋混凝土構(gòu)件在反復(fù)荷載作用下的彈塑性反應(yīng)，現(xiàn)階段在鋼筋混凝土結(jié)構(gòu)彈塑性反應(yīng)分析中，彈塑性反應(yīng)，現(xiàn)階段在鋼筋混凝土結(jié)構(gòu)彈塑性反應(yīng)分析中，這一模型得到了非常廣這一模型得到了非常廣泛的應(yīng)用。泛的應(yīng)用。在該模型中屈服前的卸載方向不指向原點(diǎn)，而是指向相反方向的開(kāi)裂點(diǎn)。在該模型中屈服前的卸載方向不指向原點(diǎn)，而是指向相反方向的開(kāi)裂點(diǎn)。屈服后卸載，剛度發(fā)生降低。屈服后

5、卸載，剛度發(fā)生降低。uyuP8 .0upypp1K2K3K帶軟化段的恢復(fù)力骨架曲線(xiàn)模型的提出根據(jù)影響荷載撓度關(guān)系的影響參數(shù)，如：混凝土強(qiáng)度、縱筋配筋率、體積配箍率、軸壓比、剪跨比等參數(shù)的大量計(jì)算結(jié)果多元回歸分析得到。121122123124YAAXBBXCCXDDX特征參數(shù)回歸系數(shù)1.9790.8621.0681.9862.096-0.0630.1892.275-0.4852.092-0.2042.0912.070-0.2150.1192.2220.0221.298-0.0971.2801.454-0.008-0.1121.1603.919-0.7563.125-1.0463.3540.161

6、0.3053.8484.657-1.9003.559-3.7754.6210.3330.0560.7651A2A1B2B1C2C1D2D1/KK2/KK3/KKy/uyuyuP8.0upypp1K2K3K三折線(xiàn)恢復(fù)力骨架曲線(xiàn)模型的驗(yàn)證40801201602000100200300400 全過(guò)程數(shù)值計(jì)算結(jié)果 回歸分析結(jié)果水平力P(kN)水平位移(mm)fcu=30MPa n=0.25 c=2 s=1.42%40801201602000100200300400 全過(guò)程分析結(jié)果 回歸分析結(jié)果水平力P(kN)水平位移(mm)fcu=30MPa; n=0.30 c=3 s=1.42%-60-40-200

7、204060-80-60-40-20020406080 計(jì)算曲線(xiàn)試驗(yàn)曲線(xiàn) P (kN)位移(mm)Z3-80-60-40-20020406080-100-80-60-40-20020406080100 試驗(yàn)曲線(xiàn)計(jì)算曲線(xiàn) P(kN)位移(mm)Z4 三折線(xiàn)恢復(fù)力滯回模型的驗(yàn)證第一次程序大作業(yè)1.自己查找柱構(gòu)件偽靜力試驗(yàn)的文獻(xiàn)(方柱、圓柱、組合柱)2.針對(duì)試驗(yàn)構(gòu)件建立分析模型3.在OpenSees中進(jìn)行滯回曲線(xiàn)的模擬4.完成研究報(bào)告（包括試驗(yàn)介紹，分析步驟，分析結(jié)果的討論Word文件）5.準(zhǔn)備PPT6.最終上交文檔內(nèi)容（PPT，Word文件及紙質(zhì)研究報(bào)告一份）7.時(shí)間：下周四（討論）鋼筋混凝土圓柱

8、抗震性能試驗(yàn)和有限元模擬鋼筋混凝土圓柱抗震性能試驗(yàn)和有限元模擬 為研究鋼筋混凝土圓柱的抗震性能，對(duì)軸壓比為為研究鋼筋混凝土圓柱的抗震性能，對(duì)軸壓比為0.45的鋼筋混凝土圓的鋼筋混凝土圓柱進(jìn)行了柱進(jìn)行了偽靜力實(shí)驗(yàn)偽靜力實(shí)驗(yàn)。基于?；贠penSEES的纖維模型，對(duì)柱的纖維模型，對(duì)柱水平力水平力-位移曲位移曲線(xiàn)線(xiàn)進(jìn)行了有限元模擬。進(jìn)行了有限元模擬。1.1 試驗(yàn)工況及材料的力學(xué)性能試驗(yàn)工況及材料的力學(xué)性能1. 試驗(yàn)研究試驗(yàn)研究試驗(yàn)工況試驗(yàn)工況砼強(qiáng)度砼強(qiáng)度f(wàn)c（Mpa)軸壓比軸壓比材料性能材料性能箍筋箍筋縱筋縱筋H-156.20.45屈服屈服/抗拉強(qiáng)度抗拉強(qiáng)度（MPa)402353.4H-272.90

9、.45彈性模量（彈性模量（MPa)2.011052.01 1051.2 試件尺寸及配筋試件尺寸及配筋 圓柱高為圓柱高為1260mm，直徑為，直徑為180mm。具體尺寸和配筋參數(shù)如圖所示。具體尺寸和配筋參數(shù)如圖所示。1.3 加載方案加載方案建研式加載裝置示意建研式加載裝置示意先將豎向軸力加至預(yù)定軸壓比先將豎向軸力加至預(yù)定軸壓比, 然后開(kāi)始水平加載試驗(yàn)然后開(kāi)始水平加載試驗(yàn), 試驗(yàn)過(guò)程中隨時(shí)調(diào)整油壓以保持軸力的試驗(yàn)過(guò)程中隨時(shí)調(diào)整油壓以保持軸力的穩(wěn)定。水平加載制度穩(wěn)定。水平加載制度: 模型柱受拉縱筋屈服前模型柱受拉縱筋屈服前,采用采用力控制加載力控制加載, 每級(jí)荷載循環(huán)一次每級(jí)荷載循環(huán)一次; 縱筋受拉

10、屈縱筋受拉屈服后改為服后改為位移控制加載位移控制加載, 每級(jí)加載位移為屈服位移的每級(jí)加載位移為屈服位移的0.5倍倍, 各級(jí)位移下循環(huán)兩次以考察模型柱各級(jí)位移下循環(huán)兩次以考察模型柱承載力及剛度的退化情況。當(dāng)試件承載力下降到極限承載力的承載力及剛度的退化情況。當(dāng)試件承載力下降到極限承載力的80%以下時(shí)以下時(shí), 認(rèn)為試件已喪失承認(rèn)為試件已喪失承載能力而結(jié)束試驗(yàn)。載能力而結(jié)束試驗(yàn)。1.4 試驗(yàn)結(jié)果試驗(yàn)結(jié)果H-1的水平力和位移滯回曲線(xiàn)的水平力和位移滯回曲線(xiàn)H-2的水平力和位移滯回曲線(xiàn)的水平力和位移滯回曲線(xiàn)2. OpenSEES有限元分析有限元分析2.1 OpenSEES分析介紹分析介紹 建模要求先將所分

11、析的對(duì)象劃分結(jié)點(diǎn)和單元，再定義荷載，施加約建模要求先將所分析的對(duì)象劃分結(jié)點(diǎn)和單元，再定義荷載，施加約束，建模命令包括：束，建模命令包括：結(jié)點(diǎn)結(jié)點(diǎn)(node)、質(zhì)量質(zhì)量(mass)、材料材料(material)、截面截面(section)、約束約束(constraints)、單元單元(element)、荷載荷載(load pattern)等。通等。通過(guò)這些命令建立所分析對(duì)象的結(jié)點(diǎn)位置、約束形式、材料特性、截面恢過(guò)這些命令建立所分析對(duì)象的結(jié)點(diǎn)位置、約束形式、材料特性、截面恢復(fù)力模型、外加荷載等。復(fù)力模型、外加荷載等。2.1.1 OpenSEES建立模型建立模型2.1.2 非線(xiàn)性分析非線(xiàn)性分析 建模

12、完成后，程序進(jìn)入非線(xiàn)性工作階段。建模完成后，程序進(jìn)入非線(xiàn)性工作階段。OpenSEES有豐富的有豐富的模塊設(shè)置并完成非線(xiàn)性分析過(guò)程。這些模塊包括：非線(xiàn)性方程組的模塊設(shè)置并完成非線(xiàn)性分析過(guò)程。這些模塊包括：非線(xiàn)性方程組的約束處理方式約束處理方式(constraints)、結(jié)點(diǎn)自由度編號(hào)優(yōu)化、結(jié)點(diǎn)自由度編號(hào)優(yōu)化(numberer)、非線(xiàn)、非線(xiàn)性方程存取計(jì)算方法性方程存取計(jì)算方法(system)、積分法則、積分法則(integrator)、迭代準(zhǔn)則、迭代準(zhǔn)則(algorithm)、容差判斂精度、容差判斂精度(test)等。每個(gè)模塊又包含了多種選項(xiàng)供等。每個(gè)模塊又包含了多種選項(xiàng)供用戶(hù)靈活調(diào)用。用戶(hù)靈活調(diào)

13、用。2.1.3 結(jié)果輸出結(jié)果輸出 在結(jié)果記錄和輸出階段，在結(jié)果記錄和輸出階段，OpenSEES程序主要提供的是程序主要提供的是Recorder模塊。模塊?？蛇x的輸出、記錄選項(xiàng)包括：非線(xiàn)性時(shí)程分析過(guò)程中各時(shí)刻結(jié)點(diǎn)的可選的輸出、記錄選項(xiàng)包括：非線(xiàn)性時(shí)程分析過(guò)程中各時(shí)刻結(jié)點(diǎn)的位移、位移、速度、加速度、位移增量速度、加速度、位移增量；時(shí)程分析過(guò)程中各時(shí)刻單元的；時(shí)程分析過(guò)程中各時(shí)刻單元的桿端力、桿端桿端力、桿端變形，截面抗力、變形和剛度的變化情況。變形，截面抗力、變形和剛度的變化情況。2.2 OpenSEES TCL語(yǔ)言編程語(yǔ)言編程 基于基于OpenSEES中的纖維模型，對(duì)中的纖維模型，對(duì)1.1-1.

14、4所介紹的試驗(yàn)進(jìn)行了所介紹的試驗(yàn)進(jìn)行了OpenSEES的有限的有限元模擬，具體的編程介紹如下。元模擬，具體的編程介紹如下。2.2.1 定義量綱定義量綱(units)OpenSEES中的量綱定義分為基本量綱和其它量綱。中的量綱定義分為基本量綱和其它量綱。例：例：set NT 1.0; # 基本量綱力基本量綱力NT set kN expr 1000*$NT ; # 其它量綱其它量綱kN=1000NT定義的量綱或參數(shù)若為表達(dá)式，則用定義的量綱或參數(shù)若為表達(dá)式，則用expr $x*$y表示。表示。expr=expression注意注意: 編程中編程中#后語(yǔ)句起解釋說(shuō)明作用，為不執(zhí)行語(yǔ)句，相當(dāng)于后語(yǔ)句起

15、解釋說(shuō)明作用，為不執(zhí)行語(yǔ)句，相當(dāng)于Matlab中中 的的%；所有定義的參數(shù)如在編程中引用，前面須加所有定義的參數(shù)如在編程中引用，前面須加$號(hào)號(hào)；除特殊說(shuō)；除特殊說(shuō) 明，例子均為有限元分析中的語(yǔ)句明，例子均為有限元分析中的語(yǔ)句 。 2.2.2 建立模型建立模型wipe; #清除之前定義的所有模型清除之前定義的所有模型1. 定義模型的維度和自由度定義模型的維度和自由度model BasicBuild ndm $ndm 例：例：model BasicBuild ndm 2 ndf 3; #模型為模型為2維維3個(gè)自由度個(gè)自由度2. 定義結(jié)點(diǎn)定義結(jié)點(diǎn)node $nodeTag (ndm $coords)

16、例：例：node 1 0 0; #定義結(jié)點(diǎn)定義結(jié)點(diǎn)1的坐標(biāo)的坐標(biāo)x、y為為(0,0)xy123. 定義結(jié)點(diǎn)質(zhì)量定義結(jié)點(diǎn)質(zhì)量mass $nodeTag (ndf $massValus)例：例：mass 1 2 0; #定義結(jié)點(diǎn)定義結(jié)點(diǎn)1的的x方向質(zhì)量為方向質(zhì)量為2 (此命令在此程序中未應(yīng)用此命令在此程序中未應(yīng)用)當(dāng)進(jìn)行時(shí)程反應(yīng)分析時(shí)需要指定質(zhì)量。當(dāng)進(jìn)行時(shí)程反應(yīng)分析時(shí)需要指定質(zhì)量。4. 定義約束（邊界條件）定義約束（邊界條件）fix $nodeTag (ndf $ConstValues)例：例：fix 1 1 1 1; #定義結(jié)點(diǎn)定義結(jié)點(diǎn)1的的x、y、z三向約束，即為固定端。三向約束，即為固定端。

17、 0表示未約束，表示未約束，1表示約束表示約束。5. 定義材料的應(yīng)力定義材料的應(yīng)力-應(yīng)變關(guān)系應(yīng)變關(guān)系uniaxialMaterial Command鋼筋混凝土纖維模型示意圖鋼筋混凝土纖維模型示意圖纖維模型須分別定義截面各部分的本構(gòu)關(guān)系：纖維模型須分別定義截面各部分的本構(gòu)關(guān)系：(1)、箍筋約束砼的箍筋約束砼的本構(gòu)關(guān)系本構(gòu)關(guān)系；(2)、保護(hù)層素砼的本構(gòu)關(guān)系保護(hù)層素砼的本構(gòu)關(guān)系；(3)、鋼筋的本構(gòu)關(guān)系鋼筋的本構(gòu)關(guān)系。箍筋約束混凝土采用箍筋約束混凝土采用Mander模型模型 .colcolcoccffffff2947125422541ccssspyhldssdAff 12122rcccxrxrff 1

18、cccx secccErEE 式中各參數(shù)含義請(qǐng)查閱式中各參數(shù)含義請(qǐng)查閱 Mander 1998原始文獻(xiàn)。原始文獻(xiàn)。Concrete 01單調(diào)與滯回本構(gòu)關(guān)系單調(diào)與滯回本構(gòu)關(guān)系uniaxialMaterial Concrete01 $matTag $fpc $epsc0 $fpcu $epsu $lambda $ft $Ets單調(diào)應(yīng)力單調(diào)應(yīng)力-應(yīng)變曲線(xiàn)應(yīng)變曲線(xiàn)滯回應(yīng)力滯回應(yīng)力-應(yīng)變曲線(xiàn)應(yīng)變曲線(xiàn)(1) 定義箍筋約束混凝土的應(yīng)力定義箍筋約束混凝土的應(yīng)力-應(yīng)變關(guān)系應(yīng)變關(guān)系uniaxialMaterial Concrete01 $matTag $fpc $epsc0 $fpcu $epsu $lambda

19、$ft $Etsset IDreinf 2 ; #約束混凝土區(qū)編號(hào)約束混凝土區(qū)編號(hào)set fc1f expr ($fc0*(-1.245+2.245*sqrt(1+7.94*$y2/$fc0)-2*$y2/$fc0; #峰值應(yīng)力峰值應(yīng)力set eps1f expr 2.0*$fc1f/$Ec ; #峰值應(yīng)變峰值應(yīng)變set eps2f expr 3.0*$eps1f ; #極限應(yīng)變極限應(yīng)變set fc2f expr ($fc1f*($eps2f/$eps1f)*$R)/($R-1+pow($eps2f/$eps1f,$R) ;#極限應(yīng)力極限應(yīng)力set lambda 0.15; #極限應(yīng)變的卸載斜

20、率和初始斜率比極限應(yīng)變的卸載斜率和初始斜率比set ftu expr -$fc1u/10.0; #抗拉強(qiáng)度抗拉強(qiáng)度set Ets expr -$Ec/2 ; #抗拉軟化剛度抗拉軟化剛度uniaxialMaterial Concrete02 $IDreinf 2 $fc1f $eps1f $fp2f $eps2f $lambda$ftu $Ets;注意：模型中的各參數(shù)設(shè)置需按照命令格式中的順序，而參數(shù)符號(hào)可自注意：模型中的各參數(shù)設(shè)置需按照命令格式中的順序，而參數(shù)符號(hào)可自行定義。實(shí)際上是主程序與子程序之間的行定義。實(shí)際上是主程序與子程序之間的實(shí)參虛參實(shí)參虛參調(diào)用調(diào)用(2) 定義未約束混凝土的應(yīng)力定

21、義未約束混凝土的應(yīng)力-應(yīng)變關(guān)系應(yīng)變關(guān)系uniaxialMaterial Concrete01 $matTag $fpc $epsc0 $fpcu $epsu $lambda$ft $Ets;set IDconcu 1; #未約束砼編號(hào)未約束砼編號(hào)set fc1u $fc; #峰值應(yīng)力峰值應(yīng)力set eps1u -0.002; #峰值應(yīng)變峰值應(yīng)變set fc2u expr 0.2*$fc ; #極限應(yīng)力極限應(yīng)力set eps2u -0.005; #極限應(yīng)變極限應(yīng)變uniaxialMaterial Concrete01 $IDconcu $fc1u $eps1u $fc2u $eps2u $lam

22、bda $ft $Ets;Steel 02單調(diào)與滯回本構(gòu)關(guān)系單調(diào)與滯回本構(gòu)關(guān)系uniaxialMaterial Steel02 $matTag $Fy $E $b $R0 $cR1 $cR2 $a $b $c $d;Steel 02單調(diào)應(yīng)力單調(diào)應(yīng)力-應(yīng)變曲線(xiàn)應(yīng)變曲線(xiàn)Steel 02滯回應(yīng)力滯回應(yīng)力-應(yīng)變曲線(xiàn)應(yīng)變曲線(xiàn)(3) 定義鋼筋的應(yīng)力定義鋼筋的應(yīng)力-應(yīng)變關(guān)系應(yīng)變關(guān)系uniaxialMaterial Steel02 $matTag $Fy $E $b $R0 $cR1 $cR2 $a $b $c $d;set IDsteel 3; #鋼筋區(qū)編號(hào)鋼筋區(qū)編號(hào)set Fy expr 353.4*$M

23、Pa; #鋼筋屈服應(yīng)力鋼筋屈服應(yīng)力set Es expr 2.0e5*$MPa; #初始彈性模量初始彈性模量set Bs 0.00085; #應(yīng)變強(qiáng)化比應(yīng)變強(qiáng)化比set R0 18.5; set cR1 0.925;set cR2 0.15; #R0、cR1、cR2是鋼筋由彈性階段向塑性階段的控制點(diǎn)是鋼筋由彈性階段向塑性階段的控制點(diǎn)6. 纖維截面單元的劃分纖維截面單元的劃分section Fiber $secTag ;patch circ $matTag $numSubdivCirc $numSubdivRad $yCenter $zCenter $intRad $extRad ;layer

24、circ $matTag $numBar $areaBar $yCenter $zCenter $radius ;砼截面單元?jiǎng)澐质疽鈭D砼截面單元?jiǎng)澐质疽鈭D鋼筋單元分布示意圖鋼筋單元分布示意圖纖維單元的劃分纖維單元的劃分section Fiber $ColSecTag ; #纖維截面纖維截面patch circ $matTag $numSubdivCirc $numSubdivRad $yCenter $zCenter $intRad $extRad ;set nfcoreC 25; #核心砼沿周長(zhǎng)方向的纖維數(shù)核心砼沿周長(zhǎng)方向的纖維數(shù)set nfcoreR 21; #核心砼沿半徑方向的纖維數(shù)核心

25、砼沿半徑方向的纖維數(shù)set nfcoverC 25; #保護(hù)層砼沿周長(zhǎng)方向的纖維數(shù)保護(hù)層砼沿周長(zhǎng)方向的纖維數(shù)set nfcoverR 2; #保護(hù)層砼沿半徑方向的纖維數(shù)保護(hù)層砼沿半徑方向的纖維數(shù)patch circ $IDconcu $nfcoverC $nfcoverR 0. 0. $Rcore $Rcol 0. 360.;#對(duì)未約束砼區(qū)進(jìn)行纖維元?jiǎng)澐謱?duì)未約束砼區(qū)進(jìn)行纖維元?jiǎng)澐謕atch circ $Idreinf $nfcoreC $nfcoreR 0. 0. 0. $Rcore 0. 360.;#對(duì)約束砼區(qū)進(jìn)行纖維元?jiǎng)澐謱?duì)約束砼區(qū)進(jìn)行纖維元?jiǎng)澐謑ayer circ $IDsteel $n

26、umBarsCol $Abar 0. 0. $Rcore 0. expr 360-360/$numBarsCol; #箍筋單元?jiǎng)澐止拷顔卧獎(jiǎng)澐?. 幾何坐標(biāo)轉(zhuǎn)換幾何坐標(biāo)轉(zhuǎn)換geomTransf PDlter $transfTag ;set ColTransfTag 1 ; #坐標(biāo)轉(zhuǎn)換編號(hào)坐標(biāo)轉(zhuǎn)換編號(hào)geomTransf PDlter $ColTransfTag ; #考慮考慮P-Dlter效應(yīng)的坐標(biāo)轉(zhuǎn)換效應(yīng)的坐標(biāo)轉(zhuǎn)換 8. 定義單元模型定義單元模型element nonlinearBeamColumn $eleTag $iNode $jNode $numIntgrPts $secTag $tr

27、ansfTag ; set numIntgrPts 5 ; #取取5個(gè)積分點(diǎn)個(gè)積分點(diǎn)element nonlinearBeamColumn 1 1 2 $numIntgrPts $ColSecTag $ColTransfTag ; #截面單元為非線(xiàn)性梁柱單元截面單元為非線(xiàn)性梁柱單元9. 結(jié)果輸出結(jié)果輸出recorder Node -dof ($dof1 $dof 2 ) $respType ;recorder Node -file $dataDir/Node2.dat -time -node 2 -dof 1 disp ;#記錄記錄2結(jié)點(diǎn)結(jié)點(diǎn)x方向的位移方向的位移10. 豎向軸力的加載豎向軸力

28、的加載pattern Plain $patternTag (TimeSeriesType arguments) load (load-command arguments)set ZYB $n ; #軸壓比為軸壓比為nset FN expr $ZYB*$fc*$PI*pow($RCol ,2)*$kN/1000 ; #軸壓比軸壓比n的豎向軸力的豎向軸力pattern Plain 1 Linear load $IDctrlNode 0 $FN 0 ; #在控制節(jié)點(diǎn)逐步增加豎向軸力至軸壓比在控制節(jié)點(diǎn)逐步增加豎向軸力至軸壓比n下的軸力下的軸力11. 軸向力的分析參量軸向力的分析參量力控制力控制set Tol 1.0e-8 ; #收斂精度收斂精度constraints Plain ; #非線(xiàn)性方程組的約束處理方式非線(xiàn)性方程組的約束處理方式numberer Plain ; #節(jié)點(diǎn)自由度編號(hào)優(yōu)化節(jié)點(diǎn)自由度編號(hào)優(yōu)化system BandGeneral ; #非線(xiàn)性存取計(jì)算方法非線(xiàn)性存取計(jì)算方法test NormDispIncr $Tol 6 0 ; #容差判斂精度容差判斂精度algorithm Newton