第五章 结构的强迫振动响应分析.docVIP

第五章 结构的强迫振动响应分析 §5.1 概述 如果结构已经用有限元方法进行了离散化，当一个结构系统受到外激励作用时，其响应就是一个多自由度系统的强迫振动问题的解。求解多自由度系统强迫振动响应的方法之一就是直接积分法。考虑到实际结构的高维数（自由度数很大）而给求解带来的困难，往往在实际求解中采用模态叠加法。直接积分法和模态叠加法这两种方法都可以得到具有相当精度的振动响应解，并且各有其特点。 §5.2 求解强迫振动响应的直接积分法 对动力学基本方程 （5－1） 进行直接积分，其含义是指在对方程进行积分之前，不对其进行任何形式的变换，在积分中，实际上是按时间步长逐步积分的。这样做的实质是基于如下考虑： 只在相隔的一些离散时间区间上、而不是在整个时间区间上的任一个时刻上满足方程，即平衡是在求解区间上的一些离散时刻上获得的。 假定位移、速度、加速度在每一个时间区间内按一定规律变化，也正是采用不同的变化形式，决定了各种直接积分解的精度、稳定性和求解速度。 首先，设表示初始时刻（）的位移、速度和加速度为已知向量，要求出从到的解，则把时间段均分为个间隔，所用的积分是在上求方程的近似解。即要在的解已知的情况下，求解时刻的解。 【中心差分法】 若基本方程式的平衡关系作为一个常系数微分方程组，则可以用任一种差分格式通过位移来表示速度和加速度。通常采用中心差分格式，这是一个行之有效的求解微分方程的格式。 （5－2） 假定 及前一时刻的位移已经求得，则将代入方程（5－1）得到： （5－3） 由此式求出 上述格式是一个显式格式。 具体计算时，还有一个步进递推格式启动的问题，即已知时，的求解问题。由的差分表达式，可求出： （5－4） 中心差分法的具体步骤为： 用有限元素法形成结构的刚度矩阵、质量矩阵、阻尼矩阵 计算初始值 选择步长，并计算积分常数 （5－5） 计算 （5－6） 形成 （5－7） 分解 （5－8） 对每一步长，进行如下计算： 求时刻的有效载荷 （5－9） 求解在时刻的位移 （5－10） 如果需要，计算t时刻的加速度和速度： （5－11） 中心差分格式使用中，一个重要的问题是步长必须小于临界步长，以保证步进递推的数值稳定性，这里，为系统的阶数，为系统最小自然周期，即： （5－12） 因此，中心差分格式是条件稳定的。 中心差分法作为显式算法的优点是，当质量阵为对角阵，阻尼阵也可以对角化时，可以避免矩阵求逆运算，而矩阵的分解运算非常简单，特别在进行非线性系统的响应分析时，由于需要在每个时间增量步修改刚度矩阵，采用中心差分法可以避免每一增量步对刚度矩阵的分解。 【Houbolt方法】 Houbolt方法也是一种差分方法，它是基于拉格朗日插值公式的步进方法，该方法利用向后差分，由位移导出速度和加速度的多步隐式公式。Houbolt方法得到的计算结果比较光滑。其差分的格式为： （5－13） 考虑时刻的平衡方程 （5－14） 从而有： （5－15） 显然，要求解必须知道 在使用Houbolt方法时，不是用此格式求初始两个时间步上的位移响应，而是用其它方法如中心差分法，步长取的几分之一来求得。 Houbolt方法是一个隐式差分格式，其步长可以取得比中心差分法大一些，不受的限制，但为了保证计算精度，步长也不能取得太大。 Houbolt方法的求解步骤总结如下： 1．用有限元素法形成结构的刚度矩阵、质量矩阵、阻尼矩阵 2．计算初始值 3．取步长，求 （5－16） 求（一般用一个特殊的启动格式如，中心差分法） 5．求解 （5－17） 6．矩阵分解 （5－18） 对每一步长，求 （1） （5－19） （2） （5－20） （3） 若需要，求出 （5－21） （5－22） 【Wilson－法】 Wilson－法假定加速度从时刻到时刻为线性变化，所以，可以认为它是线性加速度法的推广。在时，它是无条件稳定的。通常取。 具体方法为： 令为时间增量，其中，则在时刻到时刻的区间，有： （5－23） 积分上式：