工程振动基础 教学课件 ppt 作者 任兴民 秦卫阳第6章.ppt

模态迭加法是建立在展开定理基础上的一种求解动力响应的近似方法。对于一个n维自由度系统，可以通过模态坐标的解耦，确定模态坐标响应，然后通过线性变换得到物理坐标响应。 模态迭加法的计算步骤可概括如下： ①求解系统特征值问题，形成模态变换矩阵； ②通过模态坐标变换，解耦系统方程； ③求解解耦后单自由度系统的模态坐标响应解； ④进行模态坐标反变换，将模态坐标下的解变换到物理坐标。这一过程也是将各模态坐标下的响应解进行迭加，得到结构系统总的响应解。 6.3.2 中心差分法 将加速度 和速度 用中心差分法代替 将式（6-87）和式（6-88）分别代人式（6-89）中， 可得 由式（6-90）可知，只要利用t时刻的状态变量计算出 ，即可直接求解 。该积分过程称为显式 积分过程，有时也称为显式直接积分法。 ? 中心差分积分法在计算机上实现的时间积分格式 ? B. 对每一时间步长： 1. 计算在时刻t的有效荷载： 2. 求解在时刻 的位移: 3. 如果需要，计算在时刻t的加速度和速度: 6.3.3 Houbolt法 Houbolt方法与中心差分方法的区别是差分形式不同，它采用的是向后差分格式 ? Houbolt积分方法在计算机上实现的时间积分格式 ? B. 对于每一个时间步长： 1. 计算在时刻 的有效荷载: 2. 求解在时刻 的位移: 3. 若需要，计算在时刻 的加速度和速度: 6.3.4 Wilson-θ法 ? Wilson-θ积分方法在计算机上的实现 ? B. 对于每一个时间步长： 1. 计算在时刻 的有效荷载： 2. 求解在时刻 的位移： 3. 若需要，计算在时刻 的加速度和速度： 6.3.5 Newmark法 Newmark法也是线性加速度法的推广，采用如下形式的近似表达式 其中α和δ是参数，根据积分精度和稳定性的要求来确定这两个参数。当δ=1/2和α=1/6时，上面两式相应线性加速度法。 一般取δ=1/2，α=1/4，在这种情况下，它是一种定常平均加速度法（如图6-2），是无条件稳定的。 ? Newmark积分法在计算机上的实现 ? B. 对于每一个时间步长： 1. 计算在时刻 的有效荷载： 2. 求解在时刻 的位移： 3. 若需要，计算在时刻 的加速度和速度： 6.3.6 方法的选择 响应求解方法的选择取决的因素有：载荷、结构、精度要求、非线性影响程度、方法的稳定性等。 有限元法是一种将连续系统离散化的方法。这种方法先把复杂结构分割成若干个彼此之间只在结点处相互连接的单元，每个单元都是一个弹性体，单元内位移用节点位移插值函数来表示。对每个单元，由位移插值函数和动力学基本原理确定刚度矩阵、质量矩阵及其它特征矩阵。 6.4.1 结构的离散化 先把梁分成S个单元，再分别对单元和结点进行编号（见图6-3）。对于平面梁，结点的位移为结点所在截面的挠度和转角，并被取作结点的广义坐标。每个结点有两个自由度 6.4.2 单元特性分析 在梁上任取一单元，如图6-4所示。图中坐标原点取在单元左端，x轴沿单元轴线。这种与单元相联系的坐标系也称为局部坐标系。设单元长度为l，抗弯刚度为El(x)，单位长度的质量为ρ，由于单元长度一般较小，所以 El(x) 、ρ均当作常量。 梁的挠曲线微分方程为 将边界条件代入式（6-120）有 即梁的位移 是4个结点位移的线性组合，其系数由下式确定 单元内位移的表达式为 建立了单元的位移模式后，其动能、势能均可用结 点的位移表示。单元的动能为 6.4.3 结构的综合 例6-1 如图6-6所示，试用有限元分析长为2l的均匀悬臂梁的固有频率。 解：为简单起见，把梁分成两个长度均为l的单元。全梁共有3个结点，取总结点位移列阵q为 6.4.4 建立分析模型考虑的主要因素 总之，一个好的分析模型的确定，依赖于对研究对象的综合考虑分析、采用的求解方法以及精度要求。综合平衡有效性与经济性，经过实验验证逐步修正，才能完善、合理地解决。对于一个结构物，在不同的阶段，任务要求的不同，