張俊峰，　高健利

(鄭州大學 土木工程學院　河南 鄭州 450001)

?

基于OpenSees的空間桁架后屈曲材料本構(gòu)模型

張俊峰，高健利

(鄭州大學 土木工程學院河南 鄭州 450001)

摘要：空間桁架結(jié)構(gòu)在軸向壓力作用下容易發(fā)生失穩(wěn)，采用合理的鋼材本構(gòu)關(guān)系進行分析是準確追蹤桁架屈曲后行為的有效方法.借助OpenSees軟件平臺，基于桁架非彈性后屈曲單軸材料本構(gòu)模型，根據(jù)滯回準則，提出了該本構(gòu)模型新的實現(xiàn)方式；將模型分解為受拉屈服前卸載、受拉屈服后卸載、受壓屈服前卸載和受壓屈服后卸載4種加卸載路徑，基于材料基類編制了空間桁架彈性后屈曲、非彈性后屈曲本構(gòu)模型，通過算例驗證了這些本構(gòu)模型的正確性，并為如何采用這些模型進行空間桁架結(jié)構(gòu)分析提出了建議.

關(guān)鍵詞：空間桁架； 彈性后屈曲； 非彈性后屈曲； 滯回準則

0引言

空間桁架結(jié)構(gòu)廣泛應用于海洋平臺、輸電塔、大跨度穹頂?shù)冉Y(jié)構(gòu)中，這類結(jié)構(gòu)對非線性比較敏感，在軸向壓力作用下易發(fā)生失穩(wěn)，因此要想準確獲得空間桁架結(jié)構(gòu)的極限承載力，跟蹤關(guān)鍵節(jié)點的荷載-位移曲線，需要有效考慮桿件的后屈曲性能.在結(jié)構(gòu)非線性后屈曲分析中，有的學者[1]基于單元初彎曲進行研究，但大多數(shù)基于等直桿，如文獻[2]采用局部坐標系下的彈性與幾何剛度矩陣進行疊加來獲得等直桿單元切線剛度矩陣以考慮后屈曲性能，這類研究基于單元剛度矩陣來對桿件的后屈曲性能進行分析，單元剛度矩陣比較復雜，通用性相對較差；另一類研究方法基于材料本構(gòu)對桿件后屈曲性能進行分析，這類模型適用于各種能考慮幾何非線性的桁架單元，具有較好的通用性.目前關(guān)于桁架單元后屈曲的材料模型有兩類：彈性后屈曲材料本構(gòu)模型和非彈性后屈曲材料本構(gòu)模型.特別是非彈性后屈曲材料本構(gòu)模型能考慮構(gòu)件受壓屈服之后的剛度折減以及塑性變形，能夠很好地反映材料的非彈性行為，但這類模型如何在通用有限元上實現(xiàn)，以及模型間的對比研究還存在不足.

本文基于文獻[3]提出的本構(gòu)模型，在一些學者[4—10]對桁架單元材料本構(gòu)考慮后屈曲性能研究的基礎上，對桁架結(jié)構(gòu)的彈性后屈曲和非彈性后屈曲進行了對比分析，提出了非彈性后屈曲一種新的程序?qū)崿F(xiàn)方式，將本構(gòu)模型分解成4種加卸載路徑，并通過開關(guān)變量將它們銜接起來，程序更為方便簡潔，同時借助OpenSees軟件平臺，基于材料基類在OpenSees的材料庫中開發(fā)了空間桁架彈性后屈曲、非彈性后屈曲新的材料本構(gòu)模型類.

1考慮后屈曲的材料本構(gòu)模型

1.1彈性后屈曲本構(gòu)模型

彈性后屈曲本構(gòu)模型的加載和卸載規(guī)律也如圖1所示.

1.2非彈性后屈曲本構(gòu)模型

本文采用的非彈性后屈曲(IEPB)本構(gòu)模型[3]如圖2所示，該模型的骨架曲線方程的受壓段見式(1).

(1)

其中:σl為鋼材受壓下降段的臨界應力，X1、X2是基于構(gòu)件長細比的常數(shù)，ε′是距離非彈性開始點εib的長度.

非彈性后屈曲本構(gòu)模型的卸載規(guī)律遵循以下原則：

1) 受拉和受壓屈服前，均沿加載路徑卸載.

2) 受拉屈服后卸載按照式σ=E·(ε-εP)進行，其中,εP=εD-εy是在D點的塑性應變.

3) 受壓非彈性后屈曲段(ε<εib)按照直線BC段卸載，其中B點為卸載點，C點坐標為(0.5εy，0.5σy)[10]，如圖2所示.

圖1　彈性后屈曲本構(gòu)模型Fig.1　Elastic post-buckling constitutive model

圖2　非彈性后屈曲本構(gòu)模型Fig.2　Inelastic post-buckling constitutive model

2基于OpenSees的本構(gòu)模型實現(xiàn)

2.1OpenSees程序開發(fā)

OpenSees是在美國自然科學基金資助下，并由太平洋地震工程研究中心組織，十余所美國著名高校在加州大學伯克利分校主持下開發(fā)的地震工程數(shù)值模擬開源軟件.OpenSees中種類豐富的單元庫和材料類型、程序的開放性、計算方法的高效性和程序的持續(xù)集成等諸多優(yōu)點，吸引了國內(nèi)外眾多學者的注意，被科研機構(gòu)、眾多高校等學習和開發(fā)研究[11].本文基于單軸材料類編寫了IEPB材料類和EPB材料類.

2.2程序?qū)崿F(xiàn)方法

EPB材料類的實現(xiàn)比較簡單，本文重點對IEPB材料本構(gòu)的實現(xiàn)進行了分析，將其卸載規(guī)律分為了初次加載和卸載及再加載兩部分.

初次加載按照骨架曲線進行加載，初次加載以后均按照卸載處理，分為以下4種情況：受拉屈服前卸載、受拉屈服后卸載、受壓屈服前卸載和受壓屈服后卸載(圖3～6).

σr和εr分別代表開始卸載點的應力和應變.在編制程序時該點需要及時更新并儲存，在編制程序時設置了幾個開關(guān)變量，用來處理卸載過程中重新加載的情況.

圖3　受拉屈服前卸載Fig.3　Unloadingpath before tensile yield

圖4　受拉屈服后卸載Fig.4　Unloading path after tensile yield

圖5　受壓屈服前卸載Fig.5　Unloading path before compressive yield

圖6　非彈性后屈曲段卸載Fig.6　Unloading path after nonlinear compressive yield

3算例驗證

為了驗證提出的后屈曲材料本構(gòu)實現(xiàn)方法和基于OpenSees開發(fā)的新材料類的正確性，對兩個標準算例進行了彈性、彈性后屈曲、非彈性后屈曲3種情況的對比分析，其中進行非彈性后屈曲分析時標準算例[10]假定構(gòu)件的長細比足夠小，沒有受壓屈服平臺.

3.1六角星桁架

六角星桁架共24根桿件，所有桿件截面面積A和彈性模量E都相同，A=0.1 cm2，E=2.034×107N/cm2，σy=4×104N/cm2，幾何及邊界條件如圖7所示.分別進行彈性、彈性后屈曲和非彈性后屈曲分析，得到節(jié)點1和2的荷載位移曲線與文獻[10]的對比結(jié)果如圖8～10所示.

圖7　六角星桁架Fig.7　Star dome truss

圖8　節(jié)點1的荷載位移曲線Fig.8　Load-displacement curve of node 1

圖9　節(jié)點2的水平荷載位移曲線Fig.9　Horizontal load-displacement curve of node 2

圖10　節(jié)點2的豎向荷載位移曲線Fig.10　Vertical load-displacement curve of node 2

3.2測地圓頂桁架

測地圓頂桁架的幾何構(gòu)型[10]如圖11所示.該桁架包含156根桿件，所有桿件的截面相同且面積A=650 mm2，材料彈性模量E=6.895×104N/mm2.在節(jié)點1處承受集中荷載P，測地圓頂桁架的邊緣假定為鉸接，桁架的高度由方程x2+y2+(z+7.2)2=60.84確定.分別進行彈性、彈性后屈曲和非彈性后屈曲分析，得到節(jié)點1和2的荷載-位移曲線與文獻[10]的對比結(jié)果如圖12～14所示.

圖11　測地圓頂桁架Fig.11　Geodesic dome truss

圖12　節(jié)點1的荷載位移曲線Fig.12　Load-displacement curve of node 1

圖13　節(jié)點2的水平荷載位移曲線Fig.13　Horizontal load-displacement curve of node 2

圖14　節(jié)點2的豎向荷載位移曲線Fig.14　Vertical load-displacement curve of node 2

4結(jié)論

本文針對非彈性后屈曲本構(gòu)模型提出了一種新的程序?qū)崿F(xiàn)方式.基于OpenSees軟件，開發(fā)了彈性后屈曲EPB和非彈性后屈曲IEPB新的材料本構(gòu)類，通過對兩個經(jīng)典算例的對比分析，可以看出本文提出的程序?qū)崿F(xiàn)方式簡便有效，程序的收斂性好；編制的材料類能夠準確地模擬空間桁架結(jié)構(gòu)的彈性、非彈性后屈曲行為.

參考文獻：

[1]周臻, 吳京, 孟少平. 基于初彎曲單元的某弦支穹頂非線性穩(wěn)定承載力分析[J]. 計算力學學報, 2010, 27(4):771—726.

[2]劉樹堂. 等直桿單元切線剛度矩陣的精確分析方法[J]. 計算力學學報, 2014, 31(1):48—53.

[3]HILL C,BLANDFORDGG E,WANG S. Post-buckling analysis of steel space trusses[J].Journal of structural engineering,1989,115(4):900—919.

[4]PAPADRAKAKIS M. Inelastic post-buckling analysis of trusses[J]. Journal of structural engineering,ASCE, 1983, 109(9):2129—2147.

[5]BLANDFORD G E. Progressive failure analysis of inelastic space trussstructures[J].Computers & structures, 1996, 58(5):981—990.

[6]RAMESH G,KRISHNAMOORTHY C.Inelastic post-buckling analysis of truss structures by dynamic relaxation method[J]. International journal for numerical methods in engineering,1994,37(21):3633—3657.

[7]BLANDFORD G E. Large deformation analysis of inelastic space truss structures[J].Journal of structural engineering, ASCE, 1996,122(4):407—415.

[8]THAI H T,KIM S E. Nonlinear inelastic time-history analysis of truss structures[J]. Constructional steel research,2011, 67(12):1966—1972.

[9]YU Y,ZHAO X H,LUOY Z.Multi-snap-through and dynamic fracture based on finite particle method[J].Journal of constructional steel research,2013, 82(1):142—152.

[10]THAI H T, KIM S E. Large deflection inelastic analysis of space trusses using generalized displacement control method[J]. Journal of constructional steel research,2009, 65(10/11):1987—1994.

[11]FRANK M, GREGORY L F. Introducing a newmaterial into OpenSees[R]. Berkeley: University of California, 2000.

(責任編輯：王浩毅)

Post-buckling Material Constitutive Model of Space Truss in OpenSees Platform

ZHANG Junfeng,GAO Jianli

(SchoolofCivilEngineering,ZhengzhouUniversity,Zhengzhou450001,China)

Abstract:It is possible to occur buckling for space truss structures. The constitutive model of structural steel is the basis for tracking the buckling behavior of truss members accurately. Therefore, a new method of inelastic post-buckling(IEPB) constitutive model was proposed in OpenSees platform according to the rule of hysteresis. The model was divided into four loading and unloading paths including unloading before yield in tension,unloading after yield in tension,unloading before yield in compression and unloading after yield in compression. A procedure of elastic post-buckling (EPB) and inelastic post-buckling constitutive model of space truss was complied based on material base class.Numerical examples were presented to check the accuracy of the new material constitutive model.Suggestions of using these models to analysis was proposed.

Key words:space truss; elastic post-buckling; inelastic post-buckling; hysteresis criteria

收稿日期:2015-07-12

基金項目：河南省科技攻關(guān)計劃項目(142102310345).

作者簡介:張俊峰(1981—)，男，河南洛陽人，講師，博士，主要從事空間鋼結(jié)構(gòu)非線性研究，E-mail：zhangjunfeng@zzu.edu.cn.

中圖分類號：O212.1

文獻標志碼：A

文章編號：1671-6841(2016)01-0116-05

DOI:10.3969/j.issn.1671-6841.201507019

引用本文:張俊峰，高健利．基于OpenSees的空間桁架后屈曲材料本構(gòu)模型[J]．鄭州大學學報(理學版)，2016，48(1)：116—120.