连续体振动.doc

Vibration of Countinuous System 1. Euler—Bernoulli Beam Euler—Bernoulli运动过程中没有考虑剪切效应的影响。 Euler—Bernoull Beam：变形前垂直于梁中心线的截面在变形后仍保持垂直于梁的中心线。 Timoshenko Beam：Euler—Bernoull 梁中并没有考虑梁的剪切变形，在实际工程中，会存在梁的剪切变形，变形后截面与中心线存在一个夹角，截面的转角变为 中性面定义：+++ 变形的几何关系： 假设距离中性层距离为的层为。根据平面假设，单元体变形后层面为 其中，为变形的角度 应变的表达式为 弯矩的表达式为 其中，为截面的惯性矩。 转角的表达式，A点的转角为，B点的转角为 对于B点，假定转角对位置坐标线性变化，有 因此，弯曲的角度表示为 由于梁弯曲变形为小变形，有如下 得到 得到弯矩的表达式 1.1 Newton 2th law Equation. 取长度为的梁中的微元体研究，单元体的长度为。假定受到与位置坐标相关的载荷的作用，考虑到变截面梁，假定截面面积为。梁的密度表示为位置坐标的函数。单元体受力情况如图++所示。其中为弯矩(bending)，为剪力(shear force) 这里认为弯矩与剪力在位置坐标系下随线性变化，因此单元体右端的弯矩和剪力分别表示为如下形式 (1) (2) 假设梁在方向的位移(挠度)为时间与位置坐标的函数，写成的表达形式。在方向合力为零，即满足，有 (3) 式中，为作用在单元体上的载荷；等号右端表示惯性力项。 公式(3)略去，化简为 (4) 考虑到Euler—Bernoulli Beam假设，忽略剪切力的影响，弯矩平衡。假设单元体左端固定，有弯矩平衡关系 (5) 整理公式(5)，有 (6) 略去项有 (7) 将式(7)带入式(4)有 (8) 弯矩与挠度的关系为 (9) 式中，为杨氏模量；为梁截面对轴的惯性矩，对于矩形截面梁，表达式为。 将弯矩表达式(9)带入式(8)，有 (10) 公式(10)即为梁强迫振动的微分方程形式，若梁受到的外力同时随时间与位置坐标的影响，方程右侧的力可以写成的形式。 方程(10)总阶次为6，其中包含对空间坐标的4阶导数和对时间的2阶导数，求解需要6个条件。在动力学领域，采用边界条件提供求解条件。提供四个边界条件和两个初始条件，可以确定方程(10)的解。 对于求解梁的固有频率问题和自由振动问题，方程(10)右端的力。对于方程(10)的求解采用一些列正交基函数序列进行求解，即所谓的Garlerkin法。求解微分方程采用变量分离法，基于这种方法，将式（10）的解写成如下形式 (11) 其中，表示与位置坐标有关的函数；表示与时间相关的函数。 将式(11)带入式(10)并令，有 (12) 式(12)可以写成如下形式 (13) 公式（13）中，方程左侧与时间有关，方程右侧与位置坐标有关。若要求方程成立，方程两端等于同一常数，令常数为。则方程(13)有 (14) 进一步写成如下形式 (15) 观察方程（15-a），发现方程的解应为三角函数组合形式 (16) 这里进行分析时，假定梁是各项同性材料，密度为常数，梁截面积不发生变化，即，方程（15-b）可以写成如下形式 (17) 其中， 方程(17)的解可以表达为 (18) 带入方程（17）得到特征方程为 (19) 求解方程(19)得到4个特征根为 (20) 因此，方程(17)的解的表达式为 (21) 根据Euler公式如下关系 (22) 将欧拉公式（22）带入方程（21）得到梁模型的振型函数为 (23) 式中，表示常数，可以通过几何边界条件确定。 因此，梁模型的振动方程的解的形式为 可以通过方程（16）和方程（23）确定梁的固有频率和振型。 1.2 梁模型的几何边界条件 Boundary condition 常见的梁模型的支撑方式有：两端固支、两端自由、两端简支、固支自由、固支简支等，这些支撑方式可以简化为几何边界条件，梁模型常见的几何边界条件有：固定边界条件、自由边界条件、铰支边界条件等。 (1) 固定端： 位移和转角为零 (2) 自由端： 弯矩和剪力等于零 (3) 简支端： 位移和弯矩为零 1.3 悬臂梁cantilever beam 悬臂梁是指一端固定，一端自由支撑的梁模型，如图所示 图 悬臂梁示意图 悬臂梁的边界条件为： 固定端位移和转角为零，自由端弯矩和剪力为零。假设惯性矩，有边界条件表达式