受轴向力影响欧拉梁动力反应分析.doc

考虑轴向力影响欧拉梁
动力反应分析

姓名: 胡斌
专业:结构工程
日期:2011-12-12
目录
题目要求:简述分布参数体系方程的建立、圆频率和振型的计算、振型正交性的理解、振型叠加法的建立及应用、找一算例编程并画出2~3个图。
第一节前言
第二节考虑轴向力影响的梁的弯曲振动方程
第三节梁的自振频率和振型
第四节振型的正交性
分布质量的正交性
分布刚度的正交性
第五节振型叠加法
无阻尼体系振型叠加
有阻力体系振型叠加
第六节算例
第七节 matlab编程
matlab程序
图形
分析
第八节总结
第一节前言
在前面的学****过程中,我们学****了多自由度离散体系的任意结构动力反应分析,但是对实际结构来说,本质上都是具有分布质量的弹性体,即分布参数体系。要描述这些弹性体系任意瞬时的空间位置,严格上说需要无限多个广义坐标,这样的体系称为无限自由度体系。要严格描述无限自由度体系的振动,需要建立位移关于空间位置坐标和时间两个独立变量的连续函数,因此,描述无限自由度体系的运动方程为偏微分方程。
连续结构体系可按描绘它们动力行为分布所需的独立变量数来分类。但本文讨论的梁结构或轴向变形的杆,属于一维结构,它们的物理性质和动力反应可用单独一个坐标,于是这种体系的偏微分方程只包含两个独立变量,即时间和沿轴的距离。
第二节考虑轴向力影响的梁的弯曲振动方程

已知:一非均匀简支梁,沿梁长度x方向变化的抗弯刚度EI(x),单位长度的质量m(x),作用在梁上的横向荷载p(x,t),作用在梁端不随时间变化的轴向力N以及梁的横向位移u(x,t)。
求:弯曲梁的弯曲振动方程。
步骤:
1、假定梁的运动为平面弯曲,并假定变形前梁的横截面在变形后仍保持为平面,且垂直于变形后的梁轴线,即符合弯曲的平截面假定。
2、取梁上任一截面x处的微段dx为隔离体。作用在其两截面上有弯矩M,剪力Q,分布外荷载P(x,t)和假定的惯性力。
3、由竖向力平衡条件,得第一个平衡方程:
4、由力矩平衡条件,对微段右截面和x轴的交点取矩,得到第二个平衡方程

整理得:
将(2)代入(4)得
根据梁的初等变形理论,梁的弯矩和曲率的关系式为:

将(6)代入(5)得:
上式即为考虑轴向力影响的梁的弯曲振动方程。
对于等截面梁,式(7)可以简化为:
第三节梁的自振频率和振型
分析轴向力影响欧拉梁的固有振动特性,这里仅讨论等截面直梁的情况。这时,梁的自由振动运动方程为式(8)的齐次方程,即:
用“’”表示对位置x的导数,用“”表示对时间t的导数,则(9)式可写成

用分离变量法求解,假定解得形式为
式中,表示振动的形状,它不随时间而变化,表示随时间变化的振幅。
将式(11)代入(10),得到:
公式两端同时除以得:
式中C=m
由此得到两个独立的常微分方程:
方程(15)是单自由度体系无阻尼自由振动方程,其解为
式中的系数可以根据初始位移q(0)和初始速度确定,即
方程(14)是四阶微分方程,设其解得形式为
将式(18)代入式(14),得
解得
其中
将(20)代入(18)式,并用三角函数和双曲函数等式代替指数函数得
式中四个常数A~D决定梁振动的形状和振幅,它们可以利用梁端的边界条件确定。对于每一个可能指定的轴向力值,方程(22)给出了相应梁的振动形式。
第四节振型的正交性
分布质量的正交性
根据功的互等定理,第n阶振型的惯性力在第m阶振型位移上所做的功等于第m阶振型的惯性力在第n阶振型位移上所作的功。用数学表达式可表示为:
当梁以某种振型振动时,其各点的位移可表示为:
由振型引起的相应的分布惯性力为
将以上四式的幅值代入式(23)得
即
对于一般工程结构,,则有
这就是分布参数简支梁关于分布质量的正交条件。
分布刚度的正交性
如果用分布刚度作为加权函数,可以得到分布参数体系关于振型的第二个正交条件
对于有轴向力的变截面梁,自由运动方程为
左边第一项
左边第二项
左边第三项
则式(31)可写成
将式(35)代入(30)得:
式(36)即是受轴向力影响的简支梁关于分布刚度的正交条件。
对式(36)进行两次分布积分,得:
对于简支梁,边界处的位移弯矩都等于零,则式(37)中得第一项和第二项都等于零,仅剩下第三项,即:
第五节振型叠加法
无阻尼体系振型叠加
振型叠加法的基本运算就是把几何位移坐标变换为用振型幅值表示的广义坐标。对于梁这样的一维连续体,这个变换的表达式为
………………………………………….(39)
式中,u(x, t)是体系的几何位移坐标;qn(t)是