第43卷第3期 西 安 建 筑 科技 大 学 学 报(自然科学版) J.Xi an Univ.of Arch.&Tech.(Natural Science Edition) Vo1.43 NO.3 2011年6月 Jun.2011 门式钢框架稳定分析的精细传递矩阵法 孙建鹏 ,李青宁 (1.西安建筑科技大学土木工程学院,陕西西安710055; 2.西部建筑科技国家重点实验室(筹),陕西西安710055) 摘 要:将精细传递矩阵法应用于门式钢框架的稳定分析,研究门式钢框架在考虑构件轴向变形及不考虑轴 向变形时的屈曲荷载.根据精细传递矩阵法原理和压杆微弯平衡状态下的平衡微分方程,建立门式钢框架稳 定分析的精细传递矩阵式.并运用精细传递矩阵法对铰支约束下的门式钢框架进行了屈曲荷载的计算.计算 结果与有限元软件ANSYS的计算结果基本~致,从而验证了本文方法是正确、有效的.运用以上理论对于不 同边界条件的门式钢框架进行了稳定分析.分析表明:门式钢框架的屈曲荷载在考虑轴向变形时与不考虑轴 向变形时是一致的.对门式钢钢框架进行稳定分析时,可不考虑构件轴向变形的影响. 关键词:精细传递矩阵法;门式钢框架;屈曲荷载;平衡微分方程 中图分类号:Tu213.2;TU451.4 文献标志码:A 文章编号:1006—7930(2011)03—0374 05 门式钢框架是工业厂房中一种常见的结构形式.由于其承受竖向荷载的框架柱均为细长杆件,当其 承受的压力达到一定值时,易产生失稳现象 引.尽管不少学者对门式钢框架的稳定问题进行了研 究[4-63,但均未分析构件轴线变形对门式钢框架稳定的影响. 本文基于精细传递矩阵原理和压杆微弯平衡状态下的微分平衡方程,建立了考虑构件轴向变形及 不考虑轴向变形时门式钢架稳定分析的精细传递矩阵式 ;并根据以上理论编制了相应的MATLAB 计算程序,对不同边界条件下的门式钢框架进行了屈曲荷载计算.计算结果表明:门式钢框架的屈曲荷 载在考虑轴向变形时与不考虑轴向变形时是一致的,从而验证了轴向变形对门式钢框架的稳定不产生 影响.表明在对门式钢框架进行稳定分析时,可以不考虑该类构件轴向变形的影响. 酣 y1 基本原理 1.1构件的变形几何关系 轴心受压杆的计算坐标系及受力图如图1所示.假定压杆的轴线是直的,材料质地均匀,且横截面沿杆长方向保持不变. 压杆在平面内仅包括6个状态矢量: S一{N,Q,M, ,口, } l/2 埘 埘 (1) 1/2 分析压杆几何物理关系. 以指向压杆的变形为正,则其轴向应变为: £一一du 一一dx (2) (z J 图1 中心压杆 当不考虑轴向变形式压应变为零,即e一0. Fig.1 The central compression bar 收稿日期:2010—07—29 修改稿日期:201卜O4—12 基金项目:国家自然科学基金资助项目(51078306);西安建筑科技大学人才科技基金项目(RC1027);高等学校博士学科点专项科研 基金资助项目(20106120110004);西安建筑科技大学重大科技项目创新基金资助项目(ZXO901);陕西省重点学科建设专 项资金资助项目(E01004) 作者简介:孙建鹏(1981一),男,山东乐陵人,博士,讲师,主要从事桥梁结构分析及抗震. 第3期 孙建鹏等:门式钢框架稳定分析的精细传递矩阵法 375 布日对位移引起的转角为: dv 一。dx (3) 挠曲率K为: K一 d2v (4) 1.2 内力和变形的物理关系 轴力与轴向应变间的关系为: N==——EAe (5) 弯矩与曲率间的关系为: Nv=一E K一一M (6) 式中:M是弯矩;E是弹性模量;A是构件的截面面积; 是截面抗弯惯性矩. o o o 一o o 1.3构件的平衡方程 0 O —O O O dN一一0 (7) d.Z O O O O 0 dQ _0 l—d(8) 一一Q+NO (9) dZ 令s==={N,Q,M,U, ,0} ,对式(2)~式(9)进行运算整理并写成矩阵的形式为: dS—As (10) O O 0 0 0 0 0 O N 其中, A== 0 0 O O 0 1 O O O 式(10)的一般解为: S( )一P S (11) 将杆件总长离散成步长为Ax的距离间隔 m],则任一位移 一kAx(k一1,2,3…),而X + 一 + Ax,根据(11)式可得S(x…)和S(x )之间的转换关系: S( +1)一Tf(Ax)S(x ) (12) 式中T (△ )一eA . 根据指数矩阵的精细算法即可以求得指数矩阵T (Ax)一 简记为T ,为门式钢框架稳定分析的 传递矩阵. 1.4转折处的状态传递矩阵 设转折处的转折角为 ,转折处两侧状态向量的传递关系为: S 一TS r (13) 一 E 一丁i 376 西安建筑科技大学学报(自然科学版) 第43卷 cos ̄‘——sine 0 0 0 0 0 sine cos ̄O 0 0 1 0 0 0 0 0 0 0 0 其中, T:: 0 0 0 0 0 0 0 0 0 cos ̄’——sine 0 0 0 sine cos ̄0 O 0 0 1 0 0 】 0 O 0 0 T是转折处的状态传递矩阵. 1.5总传递矩阵 门式钢框架上各段的场矩阵确定以后,其上第i截面的状态向量s 与初始端状态向量s o之间的传 递关系为: S 一T T …T1S。 (14) 若整个钢框架划分成n段,则有 S 一T T 一l…TlS。一TS。 (15) 式中T是整个钢框架的总体传递矩阵. 1.6边界条件 门式钢框架任一边界处截面有6个状态向量,边界处必有3个状态向量为零. 一般铰支: 』 【M一0 悬臂端: N 0 , Q:==0 l M一0 (17) 固定端: u=O 叫一0 (18) I 一。 将压杆两端的边界条件引入式(14),就可以得到压杆关于轴向荷载N的 方程厂(N)一0,运用Newton—Raphson法 n 即可求出临界荷载N… 图2 门式钢框架的 2 算例分析 算例1某两端铰支的门式钢框架在轴向荷载N作用下的计算简图如图2 所示.杆件的基本参数见表1. 计算简图 Fig.2 The plane frame of stee[portal frame z为构件的长度.考虑构件的轴向变形的影响.根据本文的 精细传递矩阵理论,编制了相应的Matlab计算程序,计算的屈曲 荷载值列于表2.作为比较,表2还列出了有限元软件ANSYS的 计算结果. 表1基本参数 Tab.1 The parameters 由表2可以看出,门式钢框架在轴向荷载N作用下,精细传递矩阵法和有限元软件ANSYS的计算 第3期 孙建鹏等:门式钢框架稳定分析的精细传递矩阵法 377 结果完全一致.由此可见,本文的精细传递矩阵法是正确、有效的. 表2结果对比 表3结果对比 Tab.2 The comparison of results Tab.3 The comparison of results /kN Considering axial deformation 0.084 0.596 0.781 Unconside ring axial deformation 0.084 0.596 0.781 为了比较轴向变形对门式钢框架屈曲荷载的影响,运用本文方法对其在 考虑及不考虑轴向变形时的屈曲荷载进行了计算.计算结果列于表3. 由表3可以看出,门式钢框架在轴向荷载N作用下的屈曲荷载,考虑轴 向变形与不考虑轴向变形是一致的. 算例2某两端固定的门式钢框架在轴向荷载N作用下的计算简图如图3 所示.杆件的基本参数与算例1相同. 同样运用本文的方法对其在考虑及不考虑轴向变形时屈曲荷载进行了计 算,结果列于表4. 由表4可以看出,在轴向荷载N作用下两端固定门式钢框架的屈曲荷 图3 门式钢框架的 载,考虑轴向变形与不考虑轴向变形是一致的. 计算简图 算例3某一端固定一端铰接的门式钢框架在轴向荷载N作用下的计算 Fig.3 The plane frame of 简图如图4所示.杆件的基本参数与算例1相同. steei portal frame 对其在考虑及不考虑轴向变形时屈曲荷载进行了计算,结果列于表5. 由表5可以看出,在轴向荷载N作用下一端固定一端铰接门式钢框架的 屈曲荷载,在考虑轴向变形与不考虑轴向变形时也是一致的. 表4结果对比 Tab.4 The comparison of results /kN 图4 门式钢框架的 计算简图 Fig.4 The plane frame of steel portal frame 3 结 语 精细传递矩阵法兼有传递矩阵法和指数矩阵精细算法的优点.精细传递矩阵法是一种正确、有效的 计算分析方法,可以用于门式钢框架的稳定性分析.算例表明:构件的轴向变形对结构的屈曲荷载没有 影响,在计算细长结构构件的屈曲荷载时不需考虑. 参考文献References I-1-1 LOANNIS G.The effect of semi—rigid joints and an elastic bracing system on the buckling load of simple rectangular steel framesEJ-].Journal of Constructional Steel Research,2005,61:1205—1225. 378 西安建筑科技大学学报(自然科学版) 第43卷 [23 ERMOPOULOS J.Buckling length of framed compression members with semi—rigid connections ,lJ].Journal of Con-structionaI Steel Research,1991,18:139—154. [3] WANG C M.Exact Solutions for Buckling of Structural Members[J].CRC Press,Boca Raton,2005. [4] 韦成龙,李斌,刘小燕.槽型宽翼受压构件剪力滞分析的传递矩阵法[J].力学与实践,2008,30(6):39—42. WEI Cheng—long,L1 Bin,Llu xiaoyan.Transfer matrix method for shear lag analysis of a beam-and-slab structure under axial loadsl,J].Mechanics in Engineering,2008,30(6):39—42. ,ls-1 MORTEZA A,T0RKAMANI M.Second—Order Elastic Plane-Frame Analysis Using Finite-Element Method[J]. Journal of Structural Engineering,1997,123(9):1225—1235. [6] ZHOU Z H。CHAN S L.Second—Order Analysis of Slender Steel Frames under Distributed Axial and Member Loads [J].Journal of Structural Engineering,1997,123(9):1187—1193. [7] 孙建鹏,李青宁.多点地震输入下结构地震反应的频域精细传递矩阵法[J].建筑结构学报,2010,31(2):48—54. SUN Jian-peng,LI Qing—ning.Precise frequency domain transfer matrix method for seismic response analysis of structures under multi—support excitationsl,J].Journal of Building Structures,2010,31(2):48—54. [8] 孙建鹏,李青宁.求解结构自振频率的精细传递矩阵法EJ].世界地震工程,2009,25(2):140—145. SUN Jian-peng,LI Qing-ning.Precise transfer matrix method for resolving natural frequencies of structuresl,J,1。 World Earthquake Engineering,2009,25(2):140—145. [9] 段忠东,沈洪宇.非粘滞阻尼系统时程响应分析的精细积分方法[J].计算力学学报,2009,46(5):638—644. DUAN Zhong—dong,SHEN Hong~yu.Time-history analysis of a non-viscous damped system using precise integra— tionEJ].Chinese Journal of Computational Mechanics,2009,46(5):638—644. [1O] 孙建鹏,李青宁.结构动力方程的离散精细积分格式I-J].西安建筑科技大学学报:自然科学版,2010,42(1):42—46. SUN Jian-peng,LI Qing—ning.Discrere precise time-integration method for structural dynamic analysisl,J].Journal of Xi an University of Architecture&Technology(Natural Science Edition),2010,42(1):42—46. [11] HE Ji—Huan.A modified Newton-Raphson method[J].Communications in Numerical Methods in Engineering 2004,2O(10):801-805. ・ Precise transfer matrix method for stability analysis of steeel portal frame SUN Jian—peng ~。LI Qing—ning (1.School of Civil Engineering,Xi an University of Architecture and Technology,Xi an 710055 China: 2.State Key Laboratory of Architecture Science and Technology in West China(XAUAT),Xi an 710055,Chin) Abstract:Considering or neglecting the axial deformation of components.the buckling loads of steel portal frame are stud— led by the precise transfer matrix rnethod(PTMM).According to the principle and slightly bending equilibrium differential equation of compression bar,the PTMM for stability analysis of steel portal frame is built under axial force.and the buck— ling loads of steel portal frame under hinge support are solved.The results were basically identical with the result by finite element method.This shows that the method is right and effedtive.Based on the algorithm mentioned above,the stability analysis of steel portal frame with dfifferent boundary conditionis has been done.Some examples indicate that the buckling loads of steel portal frame considered are identical with that which is neglected.The suggestion is giving no consideration to the axial deformation of components to analysis the stability of steel portal frame. Key words:precise transfer matrix method;steel portal frame;buckling load;equilibrium diJ‘ferential equation Biography;SUN Jian-peng,Candidate for Ph.n,Xi an 710055,P.R.China,Tel:0086—13484666479,E-mail:suNianpeng2001@163.com