CN110287637B - 一种弹塑性屈曲承载力计算方法 - Google Patents
一种弹塑性屈曲承载力计算方法 Download PDFInfo
- Publication number
- CN110287637B CN110287637B CN201910596944.2A CN201910596944A CN110287637B CN 110287637 B CN110287637 B CN 110287637B CN 201910596944 A CN201910596944 A CN 201910596944A CN 110287637 B CN110287637 B CN 110287637B
- Authority
- CN
- China
- Prior art keywords
- truss structure
- elastic
- bearing capacity
- buckling bearing
- unfavorable
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Architecture (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
本发明公开了一种弹塑性屈曲承载力计算方法,包括:S1、对在受力荷载状况下的桁架结构进行线弹性分析,确定桁架结构在该荷载状况下的最不利杆件;S2、计算并得到桁架结构最不利杆件线弹性屈曲承载力和最不利杆件弹塑性屈曲承载力;S3、根据桁架结构最不利杆件弹塑性屈曲承载力计算得到桁架结构的弹塑性屈曲承载力;S4、采用有限元软件Abaqus中的特征值屈曲分析方法和riks分析方法分别计算桁架结构的弹塑性屈曲承载力,并将得到的桁架结构的弹塑性屈曲承载力结果与计算所得的桁架结构的弹塑性屈曲承载力结果进行对比分析,以验证计算所得桁架结构的弹塑性屈曲承载力的正确性和可行性。
Description
技术领域
本发明属于土木工程的技术领域,具体涉及一种弹塑性屈曲承载力计算方法。
背景技术
桁架结构因其特殊的结构形式,其杆件主要受轴向拉力或压力,能充分发挥结构材料的作用,减轻结构的自重,故而在工程实际中得到了广泛的应用。而一般钢构件的长细比较大,故桁架结构稳定问题是桁架结构工程应用中面临的一个主要问题,同时大量学者对桁架结构的稳定性问题作了大量的研究分析[1-8]。孙焕纯等(2005)[1]通过对桁架结构屈曲稳定性理论进行讨论,提出了桁架结构线弹性屈曲承载力的计算方法。柯秋鸿等(2006)[2]对无支撑和有支撑的两种立体桁架结构体系进行了线弹性屈曲承载力分析。郭彦霖等(2010)[3]针对钢管桁架拱平面内的弹性屈曲进行了分析。乐慈等(2015)[4]结合天津梅江会展中心工程,对张玄桁架结构的线性屈曲承载力进行了分析。杜兆宇(2017)[5]应用ABAQUS有限元软件对某悬挑双向折形斜交平面桁架结构进行了线弹性屈曲承载力分析。侯亚委等(2018)[6]对某电厂大跨度干煤棚钢结构进行了线性特征值屈曲分析。Dou等(2013)[7]对圆钢管桁架拱结构采用静态平衡法,研究分析了的圆钢管桁架拱结构弹性面外屈曲荷载。Madah等(2017)[8]基于梯度的移动渐近线方法,在考虑几何缺陷的情况下利用共旋转梁方程探讨了结构的局部屈曲和桁架的整体屈曲问题。
然而,目前大量研究只研究了桁架结构的线弹性屈曲承载力或者弹性屈曲承载力,少有研究桁架结构的弹塑性屈曲承载力。
参考文献:
[1]孙焕纯,王跃方.对桁架结构稳定分析经典理论的讨论[J].计算力学学报,2005,22(3):316-319.
Sun Huanchun,Wang Yuefang.Discussion on the Classical Theory of TrussStructural Stability Analysis[J].Chinese Journal of Computational Mechanics,2005,22(3):316-319.
[2]柯秋鸿,刘锋,李丽娟,等.立体桁架拱结构的稳定性分析[J].空间结构,2006,12(2):44-48.
Ke Qiuhong,Liu Feng,Li Lijuan,et al.Stability Analysis of a Three-dimensional Truss Arch Structure[J].Spatial Structure,2006,12(2):44-48.
[3]郭彦林,郭宇飞,窦超.纯压圆弧形钢管桁架拱平面内稳定性能及设计方法[J].建筑结构学报,2010,31(8):45-53.
GuoYanlin,GuoYufei,DouChao.In-plane buckling and design of two-hingedsteel tube circular truss-arches under pure compression[J].Journal ofBuilding Structures,2010,31(8):45-53.
[4]乐慈,袁海峰,杨洁,等.天津梅江会展中心张弦桁架稳定性分析[J].建筑结构,2015,45(14):72-76.
Le Ci,Yuan Haifeng,Yang Jie,et al.Stability Analysis on Truss StringStructure of Tianjin Meijiang Convention and Exhibition Center[J].BuildingStructure,2015,45(14):72-76.
[5]杜兆宇.瑞金体育中心体育场钢屋盖整体稳定性能分析[J].建筑结构,2017,47(S1):707-711.
Du Zhaoyu.Stability Analysis of Ruijin Sports Center Stadium Steelroof[J].Building Structure,2017,47(s1):707-711.
[6]侯亚委,赵迎九,李庆建,等.某电厂大跨度干煤棚钢结构整体稳定分析[J].建筑结构,2018,48(S1):432-434.
Hou Yawei,Zhao Yingjiu,Li Qingjian,et al.Global Stability Analysis onthe Large-span Dry-coal-shed of a Power Station[J].Building Structure,48(S1):432-434.
[7]Dou C,Guo Y L,Zhao S Y et al.Elastic out-of-plane buckling load ofcircular steel tubular truss arches incorporating shearing effects[J].Engineering Structures,2013,52(9):697-706.
[8]Madah H,Amir O.Truss optimization with buckling considerationsusing geometrically nonlinear beam modeling[J].Computers&Structures,2017,192,233-247.
[9]龙驭球,包世华,袁驷,等.结构力学Ⅱ[M].第四版.北京:高等教育出版社.2001,169-197.
Long Yuqiu,Bao Shihua,Yuan Si et al.Structural mechanicsⅡ[M].thefourth edition.Beijing,China:Higher Education Press,2001,169-197.
[10]Kato S.Guide to Buckling Load Evaluation of Metal ReticulatedRoof Structures[M].International Association for Shell and SpatialStructures,2014:23-44.
[11]T.Ogawa,T.Kumagai,S.Kuruma,K.Minowa.Buckling Load of Elliptic andHyperbolic Paraboloidal Steel Single-Layer Reticulated Shells of RectangularPlan[J].IASS Journal,2008,49(1):31-36.
[12]E.Dulacska,L.Kollar.Buckling Analysis of Reticulated Shells[J].International Journal of Space Structures,2000,15(3&4):195-203.
[13]ABAQUS 6.13.(2013).Theory Reference,ABAQUS Inc..
发明内容
本发明的目的在于针对现有技术中的上述不足,提供一种弹塑性屈曲承载力计算方法,以解决或改善上述的问题。
为达到上述目的,本发明采取的技术方案是:
一种弹塑性屈曲承载力计算方法,其包括:
S1、对在受力荷载状况下的桁架结构进行线弹性分析,确定桁架结构在该荷载状况下的最不利杆件;
S2、计算并得到桁架结构最不利杆件线弹性屈曲承载力和最不利杆件弹塑性屈曲承载力;
S3、根据所述简单桁架结构最不利杆件弹塑性屈曲承载力计算得到桁架结构的弹塑性屈曲承载力;
S4、采用有限元软件Abaqus中的特征值屈曲分析方法和riks分析方法分别计算桁架结构的弹塑性屈曲承载力,并将得到的桁架结构的弹塑性屈曲承载力结果与计算所得的桁架结构的弹塑性屈曲承载力结果进行对比分析,以验证计算所得桁架结构的弹塑性屈曲承载力的正确性和可行性。
优选地,步骤S1中确定桁架结构在该荷载状况下的最不利杆件的方法为:
计算桁架结构在该荷载状况下的荷载向量{P}:
[K]{D}={P}
其中,[K]为桁架结构线弹性刚度矩阵,{D}为桁架结构在荷载向量{P}作用下桁架结构的位移向量;
对在受力荷载状况下的桁架结构进行线弹性分析,得到桁架结构各杆件在该荷载状况下的应力状态:
其中,Ni为桁架结构受压区第i根杆件的轴向力,Ai为桁架结构受压区第i根杆件的横截面面积,计算得到桁架结构受压区中杆件的应力绝对值σb最大的杆件即为该桁架结构在该荷载状况下的最不利杆件;同时,定义最不利杆件在该荷载状况下的轴向力为N0,荷载向量{P}为参考荷载P。
优选地,步骤S2中计算得到最不利杆件线弹性屈曲承载力的方法为:
([Kr]-[S])·{Δ}=0
[S]=Ple·[s]
|[Kr]-Ple·[s]|=0
其中,[Kr]为桁架结构最不利杆件的线弹性刚度矩阵,[S]为桁架结构最不利杆件在线弹性屈曲承载力作用下的几何刚度矩阵,Ple为桁架结构最不利杆件的线弹性屈曲承载力,[s]为杆件单元的几何刚度矩阵,{Δ}为桁架结构最不利杆件在其线弹性屈曲承载力作用下的位移向量。
优选地,步骤S2中计算得到简单桁架结构最不利杆件弹塑性屈曲承载力Pcr_R的方法为:
Ny=fy·A
其中,Ple为桁架结构最不利杆件线弹性屈曲承载力,Ny为桁架结构最不利杆件的屈服承载力,Λ为最不利杆件的广义长细比,fy为最不利杆件材料的屈服应力,A为最不利杆件的横截面面积。
优选地,步骤S3中计算得到桁架结构的弹塑性屈曲承载力的方法为:
根据简单桁架结构最不利杆件的屈曲承载力Pcr_R、桁架结构弹塑性屈曲承载力的比例因子和桁架结构的参考荷载,计算得到桁架结构的弹塑性屈曲承载力Pcr:
Pcr=λ·P
λ=Pcr_R/N0
其中,λ为桁架结构弹塑性屈曲承载力的比例因子,P为桁架结构的参考荷载。
本发明提供的弹塑性屈曲承载力计算方法,具有以下有益的效果:
本发明确定桁架结构最不利杆件,并计算简单桁架结构最不利杆件弹塑性屈曲承载力,最终基于简单桁架结构最不利杆件的弹塑性屈曲承载力获得桁架结构的弹塑性屈曲承载力;同时通过与有限元软件Abaqus中的特征值屈曲分析方法和riks分析方法计算得到桁架结构的弹塑性屈曲承载力进行对比验证本发明计算得到的方法的合理性和正确性。
本发明的桁架结构最不利杆件的确定充分考虑了结构布置形式、荷载分布、杆件与节点之间的连接以及材料非线性和几何非线性的影响。因此,基于简单桁架结构最不利杆件弹塑性屈曲承载力计算桁架结构弹塑性屈曲承载力的计算方法是合理的、可行的。
根据本发明所提方法计算的桁架结构弹塑性屈曲承载力与根据有限元模特征值屈曲分析方法和riks分析方法模拟的屈曲承载力结果较为接近,验证了本发明方法计算桁架结构弹塑性屈曲承载力的正确性。
桁架结构的弹塑性屈曲承载力与结构的边界条件、荷载分布以及结构形式有关,进而对于复杂桁架结构的弹塑性屈曲承载力分析还有待更进一步的分析研究。
附图说明
图1为B31单元几何图。
图2为模型-a、模型-b、模型-c和模型-d简单桁架结构边界、布置、及荷载分布。
图3为钢材应力-应变曲线。
图4为桁架结构特征值分析位移云图。
图5为桁架结构riks分析位移云图。
图6为桁架结构荷载加载曲线。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
根据本申请的一个实施例,参考图1,本方案的弹塑性屈曲承载力计算方法,包括:
S1、对在受力荷载状况下的桁架结构进行线弹性分析,确定桁架结构在该荷载状况下的最不利杆件;
S2、计算并得到桁架结构最不利杆件线弹性屈曲承载力和最不利杆件弹塑性屈曲承载力;
S3、根据所述桁架结构最不利杆件弹塑性屈曲承载力计算得到桁架结构的弹塑性屈曲承载力;
S4、采用有限元软件Abaqus中的特征值屈曲分析方法和riks分析方法分别计算桁架结构的弹塑性屈曲承载力,并将得到的桁架结构的弹塑性屈曲承载力结果与计算所得的桁架结构的弹塑性屈曲承载力结果进行对比分析,以验证计算所得桁架结构的弹塑性屈曲承载力的正确性和可行性。
以下对上述各个步骤进行详细描述
S1、确定桁架结构在该荷载状况下的最不利杆件;
对桁架结构进行弹塑性屈曲承载力分析,其首要任务是找到桁架结构在其已知荷载状况下的最不利受压杆件。因为,最不利杆件是桁架结构在该荷载状况下最早发生屈曲的杆件,进而可能导致整个桁架结构的屈曲破坏:
[K]{D}={P} (1)
其中,[K]为桁架结构线弹性刚度矩阵,{P}为桁架结构在该荷载状况下的荷载向量,{D}为桁架结构在荷载向量{P}作用下桁架结构的位移向量。
通过对桁架结构在其受力荷载状况下进行线弹性分析,可获得该桁架结构各杆件在该荷载状况下的应力状态,在其受压杆件中,根据式(2)可定义桁架结构在该荷载状况下的最不利杆件:
其中,Ni为桁架结构受压区第i根杆件的轴向力,Ai为桁架结构受压区第i根杆件的横截面面积。
计算得到桁架结构受压区中杆件的应力绝对值σb最大的杆件即为该桁架结构在该荷载状况下的最不利杆件。同时,定义最不利杆件在该荷载状况下的轴向力为N0,荷载向量{P}为参考荷载P。
分析最不利杆件与桁架结构的关系
本发明桁架结构最不利杆件的确定是综合考虑了结构的布置形式、荷载分布、杆件与节点之间的连接以及材料非线性,这些因素体现在最不利杆件的确定及最不利杆件线弹性屈曲承载力计算部分。
因此,桁架结构最不利杆件是桁架结构中的关键构件,因为最不利杆件的屈曲可能会导致桁架结构的局部屈曲、进而可能引起结构连续倒塌,甚至引起整个结构的倒塌。同时,结构的破坏是由于结构杆件破坏逐渐积累导致。综上所述,可以将最不利杆件的屈曲承载力作为分析桁架结构屈曲承载力的重要参考指标,即可通过参考杆件的屈曲承载力估算桁架结构屈曲承载力的下限值。
S2、计算并得到桁架结构最不利杆件线弹性屈曲承载力和最不利杆件弹塑性屈曲承载力;
S2.1、计算最不利杆件线弹性屈曲承载力:
根据结构力学中稳定计算理论,计算桁架结构最不利杆件的线弹性屈曲承载力:
([Kr]-[S])·{Δ}=0 (3)
[S]=Ple·[s] (4)
其中,[Kr]为桁架结构最不利杆件的线弹性刚度矩阵,[S]为桁架结构最不利杆件在线弹性屈曲承载力作用下的几何刚度矩阵,Ple为桁架结构最不利杆件的线弹性屈曲承载力,[s]为杆件单元的几何刚度矩阵,{Δ}为桁架结构最不利杆件在其线弹性屈曲承载力作用下的位移向量。因此,通过求解等式(5)即可求得桁架结构最不利杆件的线弹性屈曲承载力Ple:
|[Kr]-Ple·[s]|=0 (5)
S2.2、计算简单桁架结构最不利杆件弹塑性屈曲承载力:
计算得到桁架结构最不利杆件线弹性屈曲承载力Ple后,简单桁架结构最不利杆件的弹塑性屈曲承载力Pcr_R即可根据Dunkerley提出的压杆强度曲线计算得到:
其中:Ny为桁架结构最不利杆件的屈服承载力,可根据式(8)计算;Λ为最不利杆件的广义长细比,其中fy为最不利杆件材料的屈服应力,A为最不利杆件的横截面面积:
Ny=fy·A (8)
S3、计算得到桁架结构的弹塑性屈曲承载力;
计算得到简单桁架结构最不利杆件的屈曲承载力Pcr_R后,根据式(9)和式(10)可得到桁架结构的弹塑性屈曲承载力[10]。
Pcr=λ·P (9)
λ=Pcr_R/N0 (10)
式中:Pcr为桁架结构的弹塑性屈曲承载力,λ为桁架结构弹塑性屈曲承载力的比例因子,P为桁架结构的参考荷载。
S4、采用有限元软件Abaqus中的特征值屈曲分析方法和riks分析方法分别计算桁架结构的弹塑性屈曲承载力,并将得到的桁架结构的弹塑性屈曲承载力结果与计算所得的桁架结构的弹塑性屈曲承载力结果进行对比分析,以验证计算所得桁架结构的弹塑性屈曲承载力的正确性和可行性。
其中,桁架结构有限元模拟分析的屈曲承载力是通过采用有限元软件Abaqus进行模拟分析。桁架结构模型中各杆件均采用Abaqus中B31单元进行模拟,模型节点均采用固结;采用两种方法进行对比分析:特征值屈曲分析方法和riks分析方法,模型网格划分方式采用自由网格划分。
根据本申请的一个实施例,参考图1,梁单元B31的几何图,假定:构件的变形可以由沿长度变化的位移函数确定,并且适合于以剪切变形为主的构件建模,如深梁和细长梁进行建模。梁单元B31在每个节点处具有六个自由度:沿x向,y向以及z向的平移,以及关于x轴,y轴以及z轴的旋转。
本发明所选的桁架结构模型为单杆结构模型和单层单跨结构模型,其结构布置及荷载分布如图2所示。桁架结构采用材料为低碳钢Q345,其弹性模量E为206GPa,屈服强度σy为354N/mm2,钢材泊松比ν为0.3,各杆件杆长l为2m,杆件截面尺寸为φ100×5(mm×mm)。钢材Q345材料性质假定为理想弹塑性,其应力-应变曲线如图3所示。
桁架结构弹塑性屈曲承载力分析
表1最不利杆件轴向力及参考荷载
Table 1 The axial force of the most unfavorable member and itsreference value
如表1所示,为了便于分析,对所有分析的桁架结构模型的荷载P的大小均假设为单位荷载1N。为确定桁架结构模型的最不利杆件,可利用有限元软件Abaqus对桁架结构模型-a、模型-b、模型-c和模型-d分别进行线弹性分析。桁架结构模型-a、模型-b、模型-c和模型-d的最不利杆件在荷载{P}作用下的轴向力N0及参考荷载P见表1。
最不利杆件线弹性屈曲承载力分析
根据结构力学中结构稳定计算理论,简单桁架结构模型-a、模型-b、模型-c和模型-d的最不利杆件的弹性刚度矩阵[Kr]和单元几何刚度矩阵[s]见表2。因此,根据式(5)可分别计算得到桁架结构模型-a、模型-b、模型-c和模型-d的最不利杆件的线弹性屈曲承载力,计算结果见表3。
表2桁架结构模型最不利杆件[Kr]和[s]
Table 2 The[Kr]and[s]of the most unfavorable member of the trussstructure models
表3 桁架结构最不利杆件线弹性屈曲承载力
Table 3 The linear elastic buckling strength of the most unfavorablemember of the truss structures
最不利杆件弹塑性屈曲承载力分析
在根据式(5)计算得到桁架结构最不利杆件的线弹性屈曲承载力Ple后,根据式(6)即可计算得到桁架结构最不利杆件的弹塑性屈曲承载力Pcr_R,具体相关参数的计算及结果见表4。
表4 桁架结构最不利杆件弹塑性屈曲承载力
Table 4 The elastic-plastic buckling strength of the most unfavorablemember of the truss structures
桁架结构弹塑性屈曲承载力分析
根据式(6)计算得到桁架结构最不利杆件的弹塑性屈曲承载力Pcr_R后,桁架结构的弹塑性屈曲承载力Pcr即可根据式(9)和式(10)计算得到,其具体相关参数的计算及结果见表5。
表5 桁架结构弹塑性屈曲承载力
Table 5 The elastic-plastic buckling strength of the truss structures
桁架屈曲承载力有限元模拟分析
桁架结构模型-a、模型-b、模型-c以及模型-d应用有限元软件Abaqus特征值屈曲分析方法分析得到的桁架结构屈曲承载力Pcr(B)和应用riks分析方法得到的桁架结构屈曲承载力Pcr(R)见表6。
为便于对比分析,表中还给出了基于最不利杆件桁架结构的弹塑性屈曲承载力估算值Pcr。另外,桁架结构模型-a、模型-b、模型-c以及模型-d应用特征值屈曲分析得到的位移云图和应用riks分析得到的位移云图分别如图4和图5所示。桁架结构模型-a、模型-b、模型-c以及模型-d应用riks分析方法的荷载加载曲线如图6所示。
表6 桁架结构屈曲承载力
Table 6 The buckling strength of truss structure
从表6中可见,桁架结构模型-a、模型-c和模型-d根据本阀门给出的方法计算的弹塑性屈曲承载力与应用有限元软件Abaqus特征值屈曲分析方法和riks分析方法得到的屈曲承载力较为接近。由于特征值屈曲分析是假定结构在整个分析过程中均为弹性(未考虑材料的塑性),桁架结构模型-b在荷载达到在弹性屈曲荷载前,结构已经发生强度破坏,故基于本文方法估算的桁架模型-b的弹塑性屈曲承载力与有限元软件Abaqus特征值屈曲分析方法得到的屈曲承载力相差较大。而riks分析的结果更接近压溃荷载,其与基于最不利杆件方法估算的桁架结构屈曲承载力较为接近。
研究发现,当构件长细比较大时,结构失效以屈曲破坏为主,此时基于最不利杆件方法估算的桁架屈曲承载力与利用特征值屈曲分析所获得的屈曲承载力较为接近;当构件的长细比较小时,构件以强度破坏为主,此时基于最不利杆件方法估算的桁架结构屈曲承载力与利用riks分析结果较为接近。
对比桁架结构模型-a、模型-b、模型-c以及模型-d应用特征值屈曲分析得到的位移云图和应用riks分析得到的位移云图,由于两种分析方法的原理不同,故在模型屈曲破坏时的破坏形态并不相同,但分析得到的模型屈曲承载力较为接近,也与本发明所提方法计算桁架结构弹塑性屈曲承载力较为接近。结合图5中桁架结构模型-b应用riks分析得到的位移云图和图6中桁架结构荷载加载曲线,可发现桁架结构模型-b的破坏形态更接近杆件的受压强度破坏。同时,从图6中,可发现桁架结构模型-a、模型-c和模型-d的破坏形态以屈曲破坏为主。
本发明的桁架结构最不利杆件的确定充分考虑了结构布置形式、荷载分布、杆件与节点之间的连接以及材料非线性和几何非线性的影响。因此,基于简单桁架结构最不利杆件弹塑性屈曲承载力计算桁架结构弹塑性屈曲承载力的计算方法是合理的、可行的。
根据本发明所提方法计算的桁架结构弹塑性屈曲承载力与根据有限元模特征值屈曲分析方法和riks分析方法模拟的屈曲承载力结果较为接近,验证了本发明方法计算桁架结构弹塑性屈曲承载力的正确性。
桁架结构的弹塑性屈曲承载力与结构的边界条件、荷载分布以及结构形式有关,进而对于复杂桁架结构的弹塑性屈曲承载力分析还有待更进一步的分析研究。
虽然结合附图对发明的具体实施方式进行了详细地描述,但不应理解为对本专利的保护范围的限定。在权利要求书所描述的范围内,本领域技术人员不经创造性劳动即可做出的各种修改和变形仍属本专利的保护范围。
Claims (5)
1.一种弹塑性屈曲承载力计算方法,其特征在于,包括:
S1、对在受力荷载状况下的桁架结构进行线弹性分析,根据线弹性分析结果将受压区中杆件的应力绝对值σb最大的杆件确定为桁架结构在该荷载状况下的最不利杆件;同时,定义该荷载状况下荷载向量{P}为参考荷载P,最不利杆件在参考荷载作用下的轴向力为N0;
S2、根据最不利杆件的线弹性刚度矩阵[Kr]、几何刚度矩阵[S]以及方程|([Kr]-[S])·{Δ}=0、[S]=Ple·[s]和|[Kr]-Ple·[s]|=0,其中,{Δ}为桁架结构最不利杆件在其线弹性屈曲承载力作用下的位移向量,计算并得到桁架结构最不利杆件线弹性屈曲承载力Ple;并根据桁架结构最不利杆件线弹性屈曲承载力Ple,最不利杆件屈服承载力Ny,确定最不利杆件弹塑性屈曲承载力Pcr_R;
S3、根据所述桁架结构最不利杆件弹塑性屈曲承载力Pcr_R以及最不利杆件在参考荷载P作用下的轴向力为N0,计算得到桁架结构的弹塑性屈曲承载力Pcr,Pcr=P·Pcr_R/N0;
S4、采用有限元软件Abaqus中的特征值屈曲分析方法计算桁架结构的弹塑性屈曲承载力或采用有限元软件Abaqus中的Riks分析方法计算桁架结构的弹塑性屈曲承载力,将得到的桁架结构的弹塑性屈曲承载力结果与计算所得的桁架结构的弹塑性屈曲承载力结果进行对比分析,以验证计算所得桁架结构的弹塑性屈曲承载力的正确性和可行性。
2.根据权利要求1所述的弹塑性屈曲承载力计算方法,其特征在于,所述S1中确定桁架结构在该荷载状况下的最不利杆件的方法为:
计算桁架结构在该荷载状况下的荷载向量{P}:
[K]{D}={P}
其中,[K]为桁架结构线弹性刚度矩阵,{D}为桁架结构在荷载向量{P}作用下桁架结构的位移向量;
对在受力荷载状况下的桁架结构进行线弹性分析,得到桁架结构各杆件在该荷载状况下的应力状态:
其中,Ni为桁架结构受压区第i根杆件的轴向力,Ai为桁架结构受压区第i根杆件的横截面面积,计算得到桁架结构受压区中杆件的应力绝对值σb最大的杆件即为该桁架结构在该荷载状况下的最不利杆件;同时,定义最不利杆件在该荷载状况下的轴向力为N0,荷载向量{P}为参考荷载P。
3.根据权利要求1所述的弹塑性屈曲承载力计算方法,其特征在于,所述S2中计算得到最不利杆件线弹性屈曲承载力的方法为:
([Kr]-[S])·{Δ}=0
[S]=Ple·[s]
|[Kr]-Ple·[s]|=0
其中,[Kr]为桁架结构最不利杆件的线弹性刚度矩阵,[S]为桁架结构最不利杆件在线弹性屈曲承载力作用下的几何刚度矩阵,Ple为桁架结构最不利杆件的线弹性屈曲承载力,[s]为杆件单元的几何刚度矩阵,{Δ}为桁架结构最不利杆件在其线弹性屈曲承载力作用下的位移向量。
5.根据权利要求1所述的弹塑性屈曲承载力计算方法,其特征在于,所述S3中计算得到桁架结构的弹塑性屈曲承载力的方法为:
根据桁架结构最不利杆件的屈曲承载力Pcr_R、桁架结构弹塑性屈曲承载力的比例因子和桁架结构的参考荷载,计算得到桁架结构的弹塑性屈曲承载力Pcr:
Pcr=λ·P
λ=Pcr_R/N0
其中,λ为桁架结构弹塑性屈曲承载力的比例因子,P为桁架结构的参考荷载。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910596944.2A CN110287637B (zh) | 2019-07-03 | 2019-07-03 | 一种弹塑性屈曲承载力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910596944.2A CN110287637B (zh) | 2019-07-03 | 2019-07-03 | 一种弹塑性屈曲承载力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110287637A CN110287637A (zh) | 2019-09-27 |
CN110287637B true CN110287637B (zh) | 2020-03-10 |
Family
ID=68021880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910596944.2A Expired - Fee Related CN110287637B (zh) | 2019-07-03 | 2019-07-03 | 一种弹塑性屈曲承载力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110287637B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112347686B (zh) * | 2020-11-30 | 2024-03-26 | 中国运载火箭技术研究院 | 一种仪器安装结构的螺栓校核方法 |
CN112668083A (zh) * | 2020-12-30 | 2021-04-16 | 中国电力工程顾问集团东北电力设计院有限公司 | 基于等稳定原则的板式连接支撑计算长度控制方法 |
CN112784460B (zh) * | 2021-01-28 | 2023-05-16 | 北京理工大学重庆创新中心 | 一种力学超材料压杆稳定性分析方法 |
CN113536646B (zh) * | 2021-08-12 | 2022-03-08 | 西南交通大学 | 一种单层球壳地震失效荷载计算方法 |
CN114036677B (zh) * | 2021-11-15 | 2023-04-28 | 中国空气动力研究与发展中心超高速空气动力研究所 | 对板钢结构承载力进行分析的方法 |
CN114048607B (zh) * | 2021-11-15 | 2023-04-25 | 中国空气动力研究与发展中心超高速空气动力研究所 | 对弹性边界板极限承载力进行分析的方法 |
CN117057023B (zh) * | 2023-10-10 | 2024-01-19 | 中国农业大学 | 空间钢结构构件统一计算长度确定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101566013A (zh) * | 2009-06-03 | 2009-10-28 | 西安建筑科技大学 | 大型火力发电厂主厂房结构体系 |
CN105544725A (zh) * | 2015-11-22 | 2016-05-04 | 东北农业大学 | 一种新型拉杆式单层柱面温室网壳结构体系及其应用 |
CN106599489A (zh) * | 2016-12-16 | 2017-04-26 | 广西大学 | 空间圆管结构极限承载力分析的一次线弹性估算方法 |
CN106777694A (zh) * | 2016-12-16 | 2017-05-31 | 广西大学 | 平面圆管结构极限承载力分析的一次线弹性估算方法 |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4777040B2 (ja) * | 2005-10-24 | 2011-09-21 | 住友ゴム工業株式会社 | タイヤの走行シミュレーション方法 |
US20090254286A1 (en) * | 2005-12-02 | 2009-10-08 | Pratt & Whitney | Systems and Methods for Modeling Surface Properties of a Mechanical Component |
CN100489256C (zh) * | 2008-02-13 | 2009-05-20 | 福建省九龙建设集团有限公司 | 宽截面梁模板支撑系统及施工方法 |
CN102277899B (zh) * | 2011-04-29 | 2014-01-29 | 深圳市建筑设计研究总院有限公司 | 一种大跨空间结构复杂边界条件杆件计算长度的确定方法 |
US8781795B2 (en) * | 2011-12-12 | 2014-07-15 | King Fahd University Of Petroleum And Minerals | Strut and tie method for waffle slabs |
CN102720295B (zh) * | 2012-04-04 | 2013-07-31 | 中国航空规划建设发展有限公司 | 一种基于索穹顶张拉和承载全过程分析的预应力确定方法 |
US9639638B2 (en) * | 2014-07-10 | 2017-05-02 | Livermore Software Technology Corp. | Methods and systems for numerically simulating physical behaviors of a string drawn out of a yarn feeder |
CN104881576B (zh) * | 2015-05-15 | 2017-12-26 | 中国电力科学研究院 | 轴心受压构件的弹塑性屈曲荷载的计算方法 |
CN105045977A (zh) * | 2015-07-01 | 2015-11-11 | 许昌学院 | 一种研究抗滑桩位的三维边坡模型建立方法 |
CN105117536B (zh) * | 2015-08-12 | 2018-01-19 | 苏州热工研究院有限公司 | 一种rpv含裂纹类缺陷的简化弹塑性断裂力学分析方法 |
CN106248502A (zh) * | 2016-07-12 | 2016-12-21 | 西南交通大学 | 悬臂梁弯曲获取材料弹塑性力学性能的方法 |
CN106503400B (zh) * | 2016-11-23 | 2019-05-10 | 福州大学 | 考虑内力分项贡献度的框架结构构件重要性系数计算方法 |
CN106978873A (zh) * | 2017-05-11 | 2017-07-25 | 四川省建筑科学研究院 | 一种可拆底模钢筋桁架楼承板结构及其施工方法 |
CN109214041B (zh) * | 2018-07-19 | 2019-05-31 | 东南大学 | 一种考虑力载荷的板结构屈曲温度分析方法 |
CN109583116B (zh) * | 2018-12-10 | 2020-12-08 | 西南交通大学 | 基于多源实测信息的地下工程初始应力场动态反演方法 |
-
2019
- 2019-07-03 CN CN201910596944.2A patent/CN110287637B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101566013A (zh) * | 2009-06-03 | 2009-10-28 | 西安建筑科技大学 | 大型火力发电厂主厂房结构体系 |
CN105544725A (zh) * | 2015-11-22 | 2016-05-04 | 东北农业大学 | 一种新型拉杆式单层柱面温室网壳结构体系及其应用 |
CN106599489A (zh) * | 2016-12-16 | 2017-04-26 | 广西大学 | 空间圆管结构极限承载力分析的一次线弹性估算方法 |
CN106777694A (zh) * | 2016-12-16 | 2017-05-31 | 广西大学 | 平面圆管结构极限承载力分析的一次线弹性估算方法 |
Non-Patent Citations (4)
Title |
---|
"Topological sensitivity to diagonal member flips of two-layered statically determinate trusses under worst loading";Z. Rychter .etal;《ELSEVIER International Journal of Solids and Structures》;20071231;第4942–4957页 * |
"基于应变能密度的网壳结构抗震性能参数分析";张明 等;《土木工程学报》;20161130;第49卷(第11期);第11-18页 * |
"基于能量的网壳结构协调抗震性能分析";张明 等;《西南交通大学学报》;20180430;第53卷(第2期);第234-243页 * |
"大跨异形钢连廊连体结构振动台试验研究";潘毅 等;《土木工程学报》;20190228;第52卷(第2期);第66-77页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110287637A (zh) | 2019-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110287637B (zh) | 一种弹塑性屈曲承载力计算方法 | |
Xiong et al. | Experimental and numerical studies on lateral torsional buckling of welded Q460GJ structural steel beams | |
Huang et al. | Experimental investigation of the flexural behavior of CFST trusses with interfacial imperfection | |
Shi et al. | Experimental investigation on bending and shear performance of two-way aluminum alloy gusset joints | |
Totaro et al. | Local buckling modelling of anisogrid lattice structures with hexagonal cells: An experimental verification | |
US10824778B2 (en) | Computer-implemented method for simulating behavior of a screw connection between metal components of engineering structures | |
Basaglia et al. | Buckling analysis of thin-walled steel structural systems using generalized beam theory (GBT) | |
CN112528415B (zh) | 一种复合材料轴结构宏-细观失效模式分析方法 | |
Ghannam | Axial load capacity of cold-formed steel built-up stub columns | |
Ma et al. | Experimental and theoretical analysis on static behavior of bolt-column joint under in-plane direction bending in single-layer reticulate shells | |
Collins et al. | Dynamic relaxation modelling of braced bending active gridshells with rectangular sections | |
Guner et al. | Pushover analysis of shear-critical frames: formulation | |
Rezaiee-Pajand et al. | Computing the structural buckling limit load by using dynamic relaxation method | |
Chen et al. | Static and dynamic analysis of Isogrid Stiffened Composite Plates (ISCP) using equivalent model based on variational asymptotic method | |
Tounsi et al. | Dynamic stiffness formulation for circular rings | |
Zhong et al. | 6A13-T6 high-strength aluminium alloy columns under axial compression: Experiments, finite element analysis and design recommendations | |
Alanjari et al. | Nonlinear joint flexibility element for the modeling of jacket-type offshore platforms | |
Wang et al. | Buckling behavior of a circular steel tube with a bolt–ball joint under installation eccentricity | |
Garcia et al. | Design charts for the local buckling analysis of integrally web-stiffened panels with filleted junctions subjected to uniaxial compressive loads | |
Al-Mosawi et al. | Optimum shape design of cold-formed thin-walled steel sections | |
Wang et al. | Static behavior of steel tubular structures considering local joint flexibility | |
Yang et al. | Analysis of the Mechanical Behavior of Bolted Beam‐Column Connections with Different Structural Forms | |
Wang et al. | Deformation performance analysis of thin plates based on a deformation decomposition method | |
Lu et al. | Mixed-dimensional coupling method for box section member based on the optimal stress distribution pattern | |
Guo et al. | Bearing capacity of bolted ball-cylinder joint under uniaxial tensile force |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200310 Termination date: 20200703 |
|
CF01 | Termination of patent right due to non-payment of annual fee |