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

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

姓名：胡斌

专业：结构工程

日期：2011-12-12

目录

题目要求：简述分布参数体系方程的建立、圆频率和振型的计算、振型正交性的理解、振型叠加法的建立及应用、找一算例编程并画出

2~3个图。

第一节

前言

第二节

考虑轴向力影响的梁的弯曲振动方程

第三节

梁的自振频率和振型

第四节振型的正交性

4.1分布质量的正交性

4.2分布刚度的正交性

第五节振型叠加法

5.1无阻尼体系振型叠加

5.2有阻力体系振型叠加

第六节算例

第七节matlab编程

7.1matlab程序

7.2图形

7.3分析

第八节总结

第一节前言

在前面的学习过程中，我们学习了多自由度离散体系的任意结构动力反应分析，但是对实际结构来说，本质上都是具有分布质量的弹性体，即分布参数体系。要描述这些弹性体系任意瞬时的空间位置，严格上说需要无限多个广义坐标，这样的体系称为无限自由度体系。要严格描述无限自由度体系的振动，需要建立位移关于空间位置坐标和时间两个独立变量的连续函数，因此，描述无限自由度体系的运动方程为偏微分方程。

连续结构体系可按描绘它们动力行为分布所需的独立变量数来分类。但本文讨论的梁结构或轴向变形的杆，属于一维结构，它们的物理性质和动力反应可用单独一个坐标，于是这种体系的偏微分方程只包含两个独立变量，即时间和沿轴的距离。

第二节考虑轴向力影响的梁的弯曲振动方程

已知：一非均匀简支梁，沿梁长度x方向变化的抗弯刚度EI(x)，单位长度的质量m（x），作用在梁上的横向荷载p（x，t），作用在梁端不随时间变化的轴向力N以及梁的横向位移u（x，t）。

求：弯曲梁的弯曲振动方程。

步骤：

1、假定梁的运动为平面弯曲，并假定变形前梁的横截面在变形后仍保持为平面，且垂直于变形后的梁轴线，即符合弯曲的平截面假定。

2、取梁上任一截面x处的微段dx为隔离体。作用在其两截面上有弯矩M，剪力Q,分布外荷载P（x，t）和假定的惯性力

dfI(x)=fI(x)dx=m(x)2,t)。

3、由竖向力平衡条件，得第一个平衡方程：

Q[P(x,t)m(x)δ2,t)]dx(Q+dx)=0 (1)

整理得：

=P(x,t)+m(x)δ2,t) (2)

4、由力矩平衡条件，对微段右截面和x轴的交点取矩，得到第二个平衡方程

M+Qdx+Nδx,t)[P(x,t)m(x)δ2,t)](dx)2(M+dx)=0.................(3)

整理得：=Q+N..................................(4)

将（2）代入（4）得

=P(x,t)+m(x)δ2,t)+Nδ2,t)...........(5)

根据梁的初等变形理论，梁的弯矩和曲率的关系式为：

M=EI(x)..........................(6)

将（6）代入（5）得：

m(x)δ2,t)+Nδ2,t)+2[EI(x)δ2,t)]=P(x,t) (7)

上式即为考虑轴向力影响的梁的弯曲振动方程。

..对于等截面梁，式（7）可以简化为：

..

mδ2,t)+Nδ2,t)+EIδ4,t)=P(x,t) (8)

第三节梁的自振频率和振型

分析轴向力影响欧拉梁的固有振动特性，这里仅讨论等截面直梁的情况。这

时，梁的自由振动运动方程为式（8）的齐次方程，即：

mδ2,t)+Nδ2,t)+EIδ4,t)=0 (9)

用“’”表示对位置x的导数，用“.”表示对时间t的导数，则（9）式可写成

..

mu+Nu+EIu=0 (10)

用分离变量法求解，假定解得形式为

u(x,t)=φ(x)q(t) (11)

式中，φ(x)表示振动的形状，它不随时间而变化，q(t)表示随时间变化的振

幅。

将式（11）代入（10），得到：

..

mq(t)φ(x)+Nφ(x)q(t)+EIφ(x)q(t)=0.......................(12)