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

考虑轴向力影响欧拉梁 动力反应分析 姓 名： 胡斌 专 业：结构工程 日 期：2011-12-12  目 录 题目要求：简述分布参数体系方程的建立、圆频率和振型的计算、振型正交性的理解、振型叠加法的建立及应用、找一算例编程并画出2~3个图。 第一节 前言 第二节 考虑轴向力影响的梁的弯曲振动方程 第三节 梁的自振频率和振型 第四节 振型的正交性 4.1 分布质量的正交性 4.2 分布刚度的正交性 第五节 振型叠加法 无阻尼体系振型叠加 5.2 有阻力体系振型叠加 第六节 算例 第七节 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）给出了相应梁的振动形式。 第四节 振型的正交性 4.1 分布质量的正交性 根据功的互等定理，第n阶振型的惯性力在第m阶振型位移上所做的功等于第m阶振型的惯性力在第n阶振型位移上所作的功。用数学表达式可表示为： 当梁以某种振型振动时，其各点的位移可表示为： 由振型引起的相应的分布惯性力为 将以上四式的幅值代入式（23）得 即 对于一般工程结构，，则有 这就是分布参数简支梁关于分布质量的正交条件。 4.2 分布刚度的正交性 如果用分布刚度作为加权函数，可以得到分布参数体系关于振型的第二个正交条件 对于有轴向力的变截面梁，自由运动方程为 左边第一项 左边第二项 左边第三项 则式（31）可写成 将式（35）代入（30）得： 式（36）即是受轴向力影响的简支梁关于分布刚度的正交条件。 对式（36）进行两次分布积分，得： 对于简支梁，边界处的位移弯矩都等于零，则式（37）中得第一项和第二项都等于零，仅剩下第三项，即： 第五节 振型叠