结构的强迫振动响应分析
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第五章 结构的强迫振动响应分析
§5.1 概述
如果结构已经用有限元方法进行了离散化，当一个结构系统受到外激励作用时，其响应就是一个多自由度系统的强迫振动问题的解。求解多自由度系统强迫振动响应的方法之一就是直接积分法。考虑到实际结构的高维数（自由度数很大）而给求解带来的困难，往往在实际求解中采用模态叠加法。直接积分法和模态叠加法这两种方法都可以得到具有相当精度的振动响应解，并且各有其特点。
§5.2 求解强迫振动响应的直接积分法
对动力学基本方程
)}({}]{[}]{[}]{[t P U K U C U M =++ （5－1）
进行直接积分，其含义是指在对方程进行积分之前，不对其进行任何形式的变换，在积分中，实际上是按时间步长逐步积分的。这样做的实质是基于如下考虑：
（1） 只在相隔t ∆的一些离散时间区间上、而不是在整个时间区间上的任一个
时刻t 上满足方程，即平衡是在求解区间上的一些离散时刻上获得的。 （2） 假定位移、速度、加速度在每一个时间区间t ∆内按一定规律变化，也正
是采用不同的变化形式，决定了各种直接积分解的精度、稳定性和求解速度。
首先，设}{}{}{0
00U U U 表示初始时刻（0=t ）的位移、速度和加速度为已知向量，要求出从0=t 到T t =的解，则把时间段T 均分为n 个间隔n T t /=∆，所用的积分是在T t t ,2,∆∆上求方程的近似解。即要在t t t ,2,∆∆的解已知的情况下，求解t t ∆+时刻的解。
【中心差分法】
若基本方程式的平衡关系作为一个常系数微分方程组，则可以用任一种差分格式通过位移来表示速度和加速度。通常采用中心差分格式，这是一个行之有效的求解微分方程的格式。
}){}({21}{}){}{2}({1
}{2
t t t t t
t t t t t t
U U t U U U U t U ∆∆∆∆∆∆-++--=+-= （5－2）
假定}{t U 及前一时刻的位移}{t t U ∆-已经求得，则将}{t U }{t U 代入方程（5－1）得到：
}]){[21][1(}]){[2]([}{}]){[21][1(
22
2t t t t t t U C t
M t U M t K P U C t M t ∆∆∆∆∆∆∆-+----=+ （5－3）
由此式求出}{t t U ∆+
上述格式是一个显式格式。
具体计算时，还有一个步进递推格式启动的问题，即}{}{}{0
00U U U 已知时，}{t U ∆-的求解问题。由}{t U }{t
U 的差分表达式，可求出： }{2
}{}{}{0
2
00U t U t U U t ∆∆∆+-=- （5－4） 中心差分法的具体步骤为：
1． 用有限元素法形成结构的刚度矩阵][K 、质量矩阵][M 、阻尼矩阵][C
2． 计算初始值}{}{}{000U U U 3． 选择步长t ∆，并计算积分常数
2
3
021201
,2,21,1a a a a t a t a ====
∆∆ （5－5） 4． 计算
}{}{}{}{0
300U a U t U U t +-=-∆∆ （5－6） 5． 形成
][][]ˆ[1
0C a M a M += （5－7） 6． 分解
T L D L M
]][][[]ˆ[= （5－8） 对每一步长，进行如下计算： 1． 求t 时刻的有效载荷
}]){[][(}]){[]([}{}ˆ{102
t t t t t U C a M a U M a K P P ∆-----= （5－9）
2． 求解在时刻t t ∆+的位移
}ˆ{}{]][][[t
t t T P U L D L =+∆ （5－10） 3． 如果需要，计算t 时刻的加速度和速度：
})
{}{(}{}){}{2}({}{1
0t
t t
t t
t
t t t t t U U a U U U U a U ∆∆∆∆+-+-+-=+-= （5－11）
中心差分格式使用中，一个重要的问题是步长必须小于临界步长n
T t n
cr =
∆，以保证步进递推的数值稳定性，这里，n 为系统的阶数，n T 为系统最小自然周期，即：
n
T t t n
cr =
≤∆∆ （5－12） 因此，中心差分格式是条件稳定的。
中心差分法作为显式算法的优点是，当质量阵为对角阵，阻尼阵也可以对角化时，可以避免矩阵求逆运算，而矩阵的分解运算非常简单，特别在进行非线性系统的响应分析时，由于需要在每个时间增量步修改刚度矩阵，采用中心差分法可以避免每一增量步对刚度矩阵的分解。 【Houbolt 方法】
Houbolt 方法也是一种差分方法，它是基于拉格朗日插值公式的步进方法，该方法利用向后差分，由位移导出速度和加速度的多步隐式公式。Houbolt 方法得到的计算结果比较光滑。其差分的格式为：
}){2}{9}{18}{11(61}{}){}{4}{5}{2(1
}{222t t t t t t t t
t t t t t t t t
U U U U t U U U U U t
U ∆∆∆∆∆∆∆∆--+--+-+-=-+-= （5－13）
考虑t t ∆+时刻的平衡方程
}{}]{[}]{[}]{[t t t t t t t t P U K U C U M ∆∆∆∆++++=++ （5－14）
从而有：