工字形截面轻钢结构构件弯扭耦合自振分析.docx

毕业论文

工字形截面轻钢结构构件弯扭耦合自振分析

专 业 工程力学
学 生
指导教师

河北工程大学土木工程学院
2013年6月10日

摘要
本文在轻钢结构满足弹性力学及薄壁杆件基本假定及工字形截面为开口薄壁结构前提下，运用所学过的材料力学、弹性力学和薄壁杆件结构力学相关的知识，采用薄壁杆件模型，通过线性插值分析方法，建立工字形截面轻钢结构构件弯扭耦合自振分析的哈密顿体系，并且采用拉格朗日能量法写出拉格朗日方程和哈密顿正则方程，最后运用精细积分法求解问题。
用MATLAB语言编制工字形截面轻钢结构构件弯扭耦合自振分析问题的计算程序，计算具体实例并对计算结果进行分析；求出问题的高精度数值解并验证本方法的可行性。运用此方法得出的结果精确度高、计算简便、计算操作易与程序化，为工程实例计算提供方便，为结构方案的确定和初步设计的完成提供了较为准确的依据。对于实际应用也有一定帮助。
关键词：工字形截面；轻钢结构；弯扭耦合自振分析；哈密顿对偶体系；线性插值

Abstract
We assume that the light steel structure is satisfied with the assumptions of the Elastic Mechanics and thin-walled structure and that the -beam is an opening thin-walled structure. Apply the knowledge which is related to the Material Mechanics， Elastic Mechanics and Thin-walled structure Mechanics. We use the model of thin-walled structure and the method of linear interpolation to build the Hamiltonian system of the -section and light steel structure member under the bending-torsional coupling vibration situation， and we can get the Lagrange equation and Hamilton’s canonical equation by applying Lagrange energy method， we solve the problem using the precise integration method of two end boundary at last.
We need to compile the procedures of the -section and light steel structure member under the bending-torsional coupling vibration situation by the MATLAB. Calculate the typical examples and analyze the computed result to obtain the highly precise numerical solution and verified this method’s feasibility. The results which obtained by this method have many advantages such as high accuracy， easy calculation， calculating operation easy to programmed and so on. It provides a new method for engineering case. And provide an accurate basis for the determination of the structural concept and the accomplishment of the initial design. It also benefits to some practical situation.
Key words: -section， light steel structure， bending-torsional coupling vibration analysis， Hamilton dual systems， linear interpolation

目 录
摘要
Abstract
1绪论 1
1.1 引言 1
1.2 轻钢结构的研究现状 1
1.3 弯扭耦合振动 2
1.4 本论文的研究内容和重点解决的问题 3
2 轻钢结构构件弯扭耦合自振的相关理论 4
2.1数值分析理论 4
2.2 哈密顿力学体系 5
2.3 哈密顿体系的求解方法 6
3 工字形截面轻钢结构构件弯扭耦合动力特性分析 8
3．1 基本假定和计算模型 8
3．2 工字形截面轻钢结构构件弯扭耦合动力特性分析的哈密顿对偶体系 9
3．3 工字形轻钢结构构件的截面特性 15
3．4 本章小结 17
4 算例计算及结果分析 18
4．1 算例验证 19
算例一 19
算例二 23
4．2 本章小结 27
5 结论与展望 28
5．1 结论 28
5．2 展望 28
鸣谢 29
参考文献 30
附录 34

工形截面轻钢结构构件弯扭耦合自振分析
学生: 指导教师：
河北工程大学土木工程学院工程力学专业
1绪论
1.1 引言
轻型钢结构是指这样一种结构:围护结构自重轻，承重结构截面小，标准化、自动化、机械化快速制作安装，采用新结构钢材、新结构体系。它分为一般轻型和超轻型钢结构。一般轻型钢结构主要采用薄钢板焊接截面或冷弯薄壁型钢构件，典型的结构体系为门式刚架，也可采用轧制型钢板截面。超轻型轻钢结构主要采用压型钢板，冷弯薄壁构件和圆钢为承重构件，典型的结构体系是摺皱拱桥屋面。轻钢结构主要体系有:焊接(轧制)门式刚架结构体系、冷弯薄壁型钢结构体系、薄壁摺皱拱桥屋面体系、多层框架结构体系、空间和张拉结构体系。
工字形截面是一种特殊的截面形式，其结构形式上基本都是开口形式，一般情况下，这种形式的截面构件都是作为梁的受力构件。而轻钢结构有两种理解范畴：（1）现行的《钢结构设计规范》中第十一章圆钢、小角钢的轻型钢结构，是指用圆钢和小于和的角钢制作的轻型钢结构，主要在钢材缺乏年代时用于不宜用钢筋混凝土结构制造的小型结构，现已基本上不大采用；（2）另一种是《门式刚架轻型房屋钢结构技术规程》所规定的具有轻型屋盖和轻型外墙的单层实腹门式刚架结构，这里的轻型主要是指围护时用的轻质材料。
钢结构是以钢材制作为主的结构，是主要的建筑结构之一。钢材的特点是强度高、自重轻、整体刚性好、变形能力强、抗震性能优越；材料均匀性和各向同性好，属理想弹性体，符合一般工程力学的基本假定。其主要用于重型车间的承重结构、受动力荷载作用的厂房结构、高层和超高层建筑等。
1.2 轻钢结构的研究现状
现代轻钢结构是指以冷弯薄壁型钢、各种H 型钢、薄钢管、焊接薄钢板变截面梁、柱等构成的结构，它除了具有传统轻钢的优点外，还更加突出了“轻”的特征，而且其应用范围也不仅仅局限于小型厂房、仓储或一些临时性的工地住所，甚至也应用于单层或多层民用住宅，办公楼等各个领域，加之其无粉尘的干施工工艺及可回收并重复利用的特点，显示出非常广阔的应用前景。
冷弯薄壁型钢住宅结构体系是一种以冷弯薄壁型钢构件和轻型板材共同作为承重和维护结构的新型轻钢龙骨结构住宅，是轻钢结构住宅的一种主要形式。
与西方发达国家相比，我国的冷成型钢结构技术的发展起步较晚，上世纪 80年代初期，国内以西安冶金建筑科技大学(现西安建筑科技大学)的何保康教授、哈尔滨建筑工程学院(现哈尔冰工业大学)的张耀春教授为代表的一批钢结构领域的学者赴美国康乃尔大学师从乔治.温特教授，这是我国冷成型钢结构领域研究的开端。
自此北美先进的冷弯型钢结构技术被引入我国。结合我国钢结构技术的实际，1987 年我国发布了第一部冷成型钢结构设计标准《冷弯薄壁型钢结构技术规范》(GBJ18-87)，1998 年开始全面修订，并于 2002 年 9 月 27 日发布修订的《冷弯薄壁型钢结构技术规程》GB50018—2002，在新规范中还增加了单层房屋设计中考虑蒙皮作用的设计原则。这是我国冷成型钢结构技术领域最权威的标准，但它的适用范围受到的限制，主要作为门式钢架工业厂房的设计参考，还没有涉及钢板厚度 2.0mm 以下的范围，也就是说在(0.4mm～2.0mm)厚度范围内的冷成型钢结构我国还没有国家和行业标准可遵循，主要还靠从国外引入。
但是一些企业如北新从自身经营的考虑编写了自己的企业标准《薄板钢骨建筑体系技术规程》；上海绿筑公司编制了企业标准《低层冷弯薄壁型钢结构施工质量验收规程》(2005 年)；在行业标准方面，2004 年建设部开始组织编写《低层轻型钢结构装配式住宅技术要求》(属产品类标准)并于 2005 年发布；中国工程建设标准化协会(CECS)组织编制的冷成型钢结构房屋体系的设计标准《低层冷弯薄壁型钢建筑技术规程》已于 2011 年 1月得到建设部批准成为行业标准，编号为 JGJ227-2011，自 2011 年 12 月 1 日起实施。最近，据悉《冷弯薄壁型钢结构技术规程》GB50018—2002 会得到补充，将扩大构件壁厚的适用范围，这将是一个重大突破。对于 2.0mm 以下薄钢板冷成型结构设计到目前为止我国还没有可直接依据的国家标准和行业标准。
近几年，国内的有些科研单位也已经有了一定的进展，随着他们与国外的科学技术和学术交流的加强，不断的把国外在冷弯轻钢结构体系方面的最新研究成果及应用实践经验带回国内，其中以西安建筑科技大学、哈尔冰工业大学及长安大学的相关教授，以及中冶院等等工程技术方面的人员为主，他们参与了国内有关冷弯轻钢结构体系的相关研究，并且取得了显著的成绩，为我国冷弯轻钢结构体系的发展及应用推广起到了推动作用。
1.3 弯扭耦合振动
实际工程结构中的许多梁由于其横截面形心和剪切中心不重合，导致其振动是弯扭耦合振动。梁作为最基本的结构元件广泛应用于结构工程的各个领域，确保它们设计的可靠性和安全性是非常重要的。经典的梁的振动理论假定弯曲振动和扭转振动是互不耦合的，这种假定仅适用于横截面具有两个对称轴的梁。工程实际中存在着大量薄壁梁(如飞机机翼、螺旋桨叶片等)仅有一个对称轴，其形心和剪切中心的不重合将导致弯曲振动和扭转振动的耦合。因此，弯扭耦合梁的精确的振动分析将是非常有必要的。
迄今为止已有不少学者研究过薄壁梁的弯扭耦合振动：
Friberg基于Vlasov梁理论(考虑了转动惯量和翘曲刚度的影响，但忽略了剪切变形的影响)研究了薄壁梁的弯扭耦合振动。
Bishop等研究了薄壁梁的耦合振动，基于考虑了翘曲刚度影响的Bernoulli-Euler梁理论，数值结果表明翘曲刚度可能对固有频率产生显著的影响。
Banerjee等建立了考虑翘曲刚度的薄壁梁的动态刚度矩阵，但他们忽略了剪切变形和转动惯量的影响。
Bercin等通过考虑剪切变形和转动惯量的影响推广了Bishop等的工作，他们表明剪切变形和转动惯量可以显著减小开口截面薄壁梁的固有频率。
Ambrosini等针对变截面薄壁Timoshenko梁进行了动力分析。
Yaman和Ozdemir针对单跨和多跨周期薄壁Bernoulli-Euler梁的耦合振动，建立了精确的解析模型。
Arpaci等基于Vlasov梁理论，解析的研究了简支薄壁梁的耦合振动。
Prokic基于薄壁Timoshenko梁理论，解析的研究了简支薄壁梁的耦合振动。
Chen和Hsiao采用常规的有限元法研究了边界条件引起的薄壁Bernoulli-Euler梁的耦合振动。
李俊研究了受轴向载荷的 Bernoulli-Euler 梁和 Timoshenko 梁弯扭耦合动力响应，分别得到了轴向受载的 Bernoulli-Euler 梁和 Timoshenko 梁固有频率和模态形状，并计算了结构的弯扭耦合动力响应；
郁殿龙在前人的基础上，对周期结构梁的弯扭耦合振动特性进行了研究，计算出了周期结构梁的弯扭耦合振动能带结构。
1.4 本论文的研究内容和重点解决的问题
1.4.1 本论文研究内容
（1）提出工字形截面轻钢结构构件弯扭耦合自振分析计算模型；
（2）通过拉格朗日能量法和哈密顿原理导出弯扭耦合自振的哈密顿对偶方程；
（3）考虑振动作用，对工字形截面轻钢结构进行自振分析；
（4）运用MATLAB软件编程，对具体实例进行计算并分析计算结果。
1.4.2 本论文重点解决的问题
（1）如何建立工字形截面轻钢结构构件弯扭耦合自振分析的计算模型；
（2）如何将基于哈密顿体系及精细积分法，应用于工字形截面轻钢结构；
（3）如何用MATLAB程序来编写例题和计算。

2 轻钢结构构件弯扭耦合自振的相关理论
本文中将轻钢结构假定为薄壁结构，所以对轻钢结构的分析可以由薄壁杆件的一些理论方法来代替和分析。
2.1数值分析理论
对于这一类问题，多年来人们寻找并发现了另一种求解途径和方法：数值解法；目前分析薄壁杆件常用的数值方法为：有限条法，有限单元法等，本文采用的是样条插值法。
2.1.1有限条法
有限条法是从有限元法发展而来的一种半解析方法，在土木工程领域得到广泛应用，其基本思路是：令求解域的一个方向为连续体，而将其沿其它方向离散为有限条元，然后选取条元的位移函数，利用最小势能原理导出有限条法的线性方程组进而得到结构的内力和位移响应。理论上只要选择合理的有限条数目和傅立叶级数函数中所保留的项数，就能达到足够的精确度，有限条法求解薄壁箱梁的一般步骤为：
（1）根据求解问题的边界条件，假定合理的有限条位移函数；
（2）按假设的位移函数导出单个有限条元在受力时的总势能，并将各条元能量的代数和叠加得到整个板的能量；
（3）由最小势能原理导出有限条的刚度矩阵和力矩阵，进而得到线性方程组；
（4）解线性方程组得到位移函数的闭合解，再由虎克定律得到薄壁箱梁翼板的应力。
有限条法能够分析薄壁箱梁所有的变形特征(弯曲、扭转、翘曲、畸变、横向弯曲等)，而且与有限元法相比，有限条法应用方便，精度高，节省大量计算时间和费用，其局限性在于目前对于解决变高度、变截面梁和不规则结构还有一定的困难。
2.1.2有限元法
近几十年来，有限单元法成为土木工程领域工程分析中最通用、最有效的方法，有时对某些结构是唯一适用的方法，用它可以分析形状复杂的各种工程结构，可以在计算中模拟各种复杂的材料本构关系、荷载和边界条件，通过完全交互式的前后处理技术实现结果的形象化显示，用有限单元法分析薄壁结构，通常采用变异的直接刚度法，此时的基本未知量为广义位移，这种方法本质上和其他力学体系的分析无多大区别，分析步骤一般为将结构纵向分为若干块，各块再分成很多单元，然后建立各单元的刚度矩阵，并变换到总体坐标系中，接着按照直接刚度法原理，将各刚度矩阵汇编为结构的刚度矩阵，建立荷载矩阵，应用边界条件，就能求解出节点的位移，最后将总体坐标系中的节点位移变换到单元坐标系中。
有限元法的计算精度依赖于单元网格的划分及单元模式的选择。目前，大型通用有限元程序(如ANSYS、ABAQUS等)非常流行，为有限单元法的应用提供了有利的条件。
2.1.3 插值法
插值法又称“内插法”，利用函数在某区间中若干点的函数值，作出适当的特定函数，在这些特定点上取已知值，在区间的其他点上用这特定函数的值作为函数的的近似值，这种方法称为插值法，插值法是函数逼近的一种重要方法。本文使用的就是这类方法。
2.2 哈密顿力学体系
力学求解的三大体系有牛顿力学、拉格朗日力学和哈密顿力学。
牛顿力学主要由三大定律组成，研究的对象是自由质点体系，可以解决基本力学问题。拉格朗日力学研究几何约束质点体系的问题，得出几何约束质点系统的拉格朗日第一类方程
(2-1)
上式是最早建立的比较完备的几何约束质点体系的数学模型。
在引入广义坐标的描述方法后，拉格朗日为理想的几何约束质点系统建立了拉格朗日第二类方程：
（2-2）
式中，——和所表示的系统动能，；
——对应于广义坐标的广义力；
——系统的自由度，；
——系统质点的个数；
——约束方程的个数。
对于有势系统，系统的势记作，当不显含时间时，特别称作系统的势能，此时的系统叫作保守系统。对于有势系统的拉格朗日函数为
(2-3)
包涵了有势系统的所有动力学行为，拉格朗日第二类方程可表示为
(2-4)
哈密顿力学在拉格朗日力学体系的基础上进行数学变换。
将拉格朗日函数积分
(2-5)
该式称为哈密顿作用量。有势系统的哈密顿变分原理可以描述为有势系统的哈密顿作用量在系统的正路上取驻值。即
(2-6)
系统的状态也可以用其他参数确定。通过勒让德变换引入变量
(2-7)
是广义动量。称为哈密顿变量。由上式可反解出
(2-8)
因此拉格朗日变量可由哈密顿变量表示。
拉格朗日函数对变量的勒让德变换为函数
(2-9)
其中需要用表示出来。函数称为哈密顿函数，也称作系统的广义能量。哈密顿建立了著名的有势系统的新数学模型哈密顿正则方程
； (2-10)
在弹性力学中，拉格朗日函数是弹性体的总势能密度，它的积分在动力学中称为哈密顿作用量，其物理意义是总势能。由最小总势能原理就可以导出相应的拉格朗日方程
(2-11)
由勒让德变换引入的对偶变量，在这里是广义力。得出的哈密顿函数叫做混合能密度。
2.3 哈密顿体系的求解方法
对于弹性体的线弹性系统，拉格朗日函数是和的函数，可写成矩阵形式
(2-12)
式中，广义位移，广义位移的导数，系统外荷载。系统的总势能为：
(2-13)
由最小势能原理，，即可以导出拉格朗日方程：
　　 　　 (2-14)
将以往的时间坐标换成空间坐标采用矩阵向量列式可写成：
(2-15)
式中。利用勒让德变换，引入广义位移的对偶变量广义力。
和的表达式分别为：
， (2-16)
导入哈密顿函数，消去，得到矩阵形式的表达式为：
(2-17)
式中，。
则哈密顿对偶方程成为：
(2-18)
(2-19)
将对偶向量合在一起组成为状态向量
(2-20)
(2-21)
其中，
，。 (2-22)
式(2-21)就是哈密顿对偶方程系统矩阵称为哈密顿矩阵，是一个辛矩阵。
对于初值问题和两端边值问题，可以结合边界条件，运用精细积分法求出。

3 工字形截面轻钢结构构件弯扭耦合动力特性分析
结构的自振周期、振型和阻尼等动力特性只是和结构自身的质量和刚度有关，是结构本身所固有的属性，是结构振动的内因，也是结构动力学中的重要特性。阻尼的大小由试验测定，自振频率及振型可通过计算来确定。结构的动力特性分析是计算结构动力响应的前提条件，为进行复杂高层建筑结构抗震设计提供必要的理论依据，在实际的工程应用上有非常重要的意义。
目前常用的求解多质点弹性体系基本频率的方法有能量法、折算质量法；最常用的求解基本周期的近似方法是顶点位移法；求解多质点弹性体系多个频率和振型的近似计算方法有矩阵迭代法，Rayleigh-Ritz 法；求解多自由度弹性体系自振周期和振型的计算机方法有雅可比法。
3．1 基本假定和计算模型
3.1.1 基本假定
（1）将工字形轻钢结构构件从底到顶划分为若干个有限宽度的条形单元，沿条形单元的竖向是真实的连续函数(广义位移)，而沿横向用插值函数来模拟翘曲位移函数:
（3-1）
式中：：截面所分条元交界线处的纵向翘曲位移；
：关于各交界线处纵向翘曲位移的插值函数，我们选用较为简单的分段线性插值函数来描述闭口薄壁杆件的纵向翘曲位移。
（2）研究薄壁杆件绕z 轴扭转时，符合弗拉索夫的刚周边假定，假定截面的切向位移为：
（3-2）
式中：为截面上任一点的切线到扭心的距离；
为截面z 的扭转角。
（3）沿曲线坐标s 方向的环向应力和法向应力远比横截面的轴向应力小，可忽略不计。
3.1.2 计算模型
在薄壁杆件的弯扭作用问题中，为了更真实的表示薄壁杆件结构的实际变形，我们采取插值函数的思想对杆件的翘曲位移进行假定，以此达到逼近杆件实际变形的效果；将等效后的工字形薄壁结构依次划分为若干个有限宽度得到条形单元，沿条形单元的竖向真实的的连续函数（广义位移），而沿横向用插值函数来模拟翘曲位移函数：

第个薄壁杆件经过沿平动以及绕扭心的三个过程以后，其上任意一点在整体坐标系下的切向位移可表示为:
（3-3）
式中：；；
；；
―分别为坐标原点沿轴、轴方向的位移；
―分别为第个截面的扭心在整体坐标系下沿轴方向的位移；
―分别为第个截面的扭心在局部坐标系下沿轴方向的位移；
―分别为第个截面的任意一点在整体坐标系下沿轴方向的位移；
―分别为第个截面任意一点在局部坐标系下沿轴方向的位移；
―为第个断面上任一点的切线到其扭心的距离；
―为第个断面绕轴的扭转角；
―分别为第个断面扭心的坐标；
3．2 工字形截面轻钢结构构件弯扭耦合动力特性分析的哈密顿对偶体系
3.1.1 工字形轻钢结构构件的线性插值函数

…
1
2
…

S
图 3-1 截面分段插值图

按上图依次写出：
，
（3-4）

式中： 为截面上第点处的坐标，设工字形截面上左上部第一个点为起点，按汉字中工的写法依次标号从1~n，上翼缘交点处的点设为，最右边的点设为；下翼缘最左边第一个点设为，交点设为，最右边点设为，如上图所示。
d：表示薄壁结构划分区段的长度，其值为：。
对上式求导可得：
，
（3-5）

3.2.2 工字形轻钢结构构件的总势能
引用向量，则式（3-1）可写成向量形式：
（3-7）
将式（3-6）代入应变表达式（3-3），即可得轴向应变和剪应变，分别为：
（3-7）
（3-8）
式中：；。
则体系总势能为：
（3-9）
3.2.3 工字形轻钢结构构件的拉格朗日方程和哈密顿对偶方程
（1）由体系总势能公式（3-9）可求得势能密度，即拉格朗日函数：
（3-10）
由于；
（3-11）
式中：；此为等效正交异形板所产生的轴向应应变能。
；；；。
，绕扭心的极惯性矩；
，圣维南扭转常数；
，沿x方向的剪切面积；
，沿y方向的剪切面积；
，混合剪切面积；
，y方向的剪切静矩；
，x方向的剪切静矩。
将上式写成拉格朗日矩阵形式：
（3-12）
式中：，，
.
，，，

将上式（3-12）代入拉格朗日静力平衡方程：， （3-13）
因为 ，所以
，，则：
（3-14）
式中；
结构进行自由振动时，在静力平衡方程式（3-13）的基础上，加上结构自由振动时的惯性力，并将相关方程改为运动方程便可得到结构的运动方程：
（3-15）
其中：
，，。
，。
为结构单元的密度，其他几何尺寸如上。
在分析动力问题时，对时间t 常采用化为频域的方法，此时采用：
（3-16）
将频域化的的式（3-16）代入式（3-15）可得结构自振问题的动力方程。
（3-17）
对比系数可得：
，，。
式中： （3-18）
（2）引入变量的对偶变量
（3-19）
由上式可求得 为：
（3-20）
进而可引入哈密顿函数：
（3-21）
将式（3-20）代入式（3-21），消去，整理后可得哈密顿函数新表达式：
（3-22）
其中：
，，，，。
由上式可导出结构的哈密顿对偶方程即哈密顿正则方程：
（3-23）
把哈密顿正则方程写成矩阵形式：
（3-24）
其中：
，， （3-25）
上式中系统矩阵称为哈密顿矩阵，它是一个辛矩阵。哈密顿的物理意义是动能和势能的和。运用MATLAB编制程序，并应用两端边值问题精细积分法进行求解计算，从而得到频率。
3．3 工字形轻钢结构构件的截面特性
对于工字形截面，截面上任意点到截面扭心的距离为：
（3-26）
S
2b
2h
2b

图3-2 工型截面图
1

（3-27）
（3-28）
（3-29）
（3-30）
（3-31）
（3-32）
（3-33）
（3-34）
（3-35）
（3-36）
（3-37）
（3-38）
（3-39）
（3-40）
3．4 本章小结
本章通过对工字形截面轻钢结构构件的分析，建立了开口薄壁杆件的受力模型。采用线性插值的方法无限逼近真实结果，把等效后的薄壁结构划分为若干个有限宽度的条形单元，并以最小总势能原理为基础，通过引入原变量的对偶变量，运用哈密顿体系导出哈密顿对偶方程，从而推导出组合薄壁杆件的弯扭耦合分析的对偶方程。本论文的重点的工作之一是对截面特性进行计算，对此本章给出了详细的介绍。

4 算例计算及结果分析
（4-1）
将式（3-4）代入式（4-1）中，可求得：
（4-2）
，
（4-2a）
（4-2b）
（4-2c）
（4-3）
4．1 算例验证
算例一
如图4-1a所示，一悬臂工字形断面薄壁杆件断面尺寸如图4-1a所示（单位m），杆长，材料弹性模量，泊松比，，求杆件约束扭转振动的自振频率。
1.6
1.6
0.18
0.14
0.14
0.94
图4-1a 截面尺寸图
x
y

1
2
4
5
9
10
21
13
14
17
3
6
7
8
11
12
15
16
18
19
20
X
Y
图 4-1b 分段插值
O

截面分为20段，每段长度为：；扭心坐标为：
，

，，，
，，。
， ，，
，，

表4-1 不同长度悬臂梁自振频率对比
长度
（）
一阶频率（）
二阶频率（）
三阶频率（）
四阶频率（）
五阶频率（）
六阶频率（）
10
189.43
289.93
490.43
590.93
791.43
1191.9
15
79.543
180.04
304.82
405.32
505.82
707.32
20
55.4
155.9
228.4
266.65
367.15
467.65
30
34.019
112.33
116
152.12
217.75
318.25

算例二
工字形薄壁断面悬臂梁的自由端承受均布荷载q，截面尺寸如图4-3b所示，，，杆长，弹性模量，泊松比，荷载值，，求杆件约束扭转振动的自振频率。
q
图 4-3a 荷载分布图

t
x
b
b
d
t
h
图4-3b 截面尺寸图
y

图 4-3c 分段插值图
1
2
3
5
6
7
12
13
15
17
…
4
14
16

截面分为16段，每段长度为：，。
，

，，，
，，。
， ，，
，，

表4-2 不同长度悬臂梁自振频率对比
长度
（）
一阶频率（）
二阶频率（）
三阶频率（）
四阶频率（）
五阶频率（）
六阶频率（）
5
83.5
393
837.5
858
1005.5
1857
10
31.5
121
282.5
419
428.5
520
15
19.5
66
140.5
249
279.5
285
20
13.5
45
88.5
151
209.5
214

4．2 本章小结
本章通过对悬臂轻钢工型截面梁的算例分析，可以看出：二阶频率比一阶频率的数值大，三阶的比二阶大，并且随着阶数的增大，其数值逐渐增大；另外，对于不同杆长的梁，在杆长分段数、截面插值数、步长值相等的情况下，频率大小随杆长的增大而减小。MATLAB的编程计算可以很大的简化计算利于分析，不过上述数值的结果尽可供参考。
5 结论与展望
5．1 结论
随着建筑行业的发展，轻钢结构在其中发辉的作用越来越显著。钢结构有着强度大和重量轻等众多优点，本文中所讨论的工型轻钢结构构件在简化为薄壁结构模型的情况下，利用薄壁杆件的相关知识推到出了拉格朗日方程和哈密顿对偶方程，再通过精细积分法求出未知量。虽然在matlab程序运行中出现了一点问题，但认可得出结论如下：
（1）通常情况下，结构不仅受到弯矩或者某一个单独广义力的作用，而是受到复合力的综合作用，本文考虑弯扭耦合自振情况，这就更加符合真实效果。
（2）将工型轻钢结构简化为薄壁结构后，求出哈密顿对偶方程，利用精细积分法计算，（3）本文中考虑了结构的轴向翘曲位移，这更加符合真实情况，也进一步增加了本文所采取方法的准确度。
5．2 展望
本文通过对工字形轻钢结构构件的简化，将其看成一开口薄壁结构。然后采用线性插值的方法来逼近结构真实的效果，拉格朗日方程和哈密顿对偶体系及两端边值精细积分法的运用也有很高的精度，我国的轻钢结构的发展也有一些问题：
（1） 轻钢结构的技术标准和规范有待进一步形成和完善。我国的标准规范是针对几十年来大量使用的结构体系编制的，近几年，根据建设部技术标准编制工作计划，虽然有一些轻钢结构技术规范出版，但在全国还尚未形成强制性规范条文。在建造试验工程并取得经验的基础上，就应不失时机地编制相应的轻钢结构住宅技术标准和规范。
（2） 轻钢建筑设计理念与国外存在一定的差距。轻钢结构建筑由于其特殊的材料和先进的设计软件，得以使设计程序实现建筑结构一体化设计，即建筑设计与结构设计同时完成，并不像混凝土建筑的设计是按着先建筑设计后结构设计的理念进行设计，国外轻钢结构企业都是遵循这样一条理念，而国内的大多数企业仍执行混凝土建设设计的原则。另外，我国轻钢结构的主体结构设计、维护系统的设计分开来进行，而国外普遍的做法是统一进行设计，甚至有些企业仍然执行建筑模数制，这就保证了建筑整体安装的精确度。
（3） 缺少研究轻钢结构的专业工程技术人员。我国在较长的一段时间内，大多从事单一的混凝