python时域转频域_理解波动（三）：从频域到时域

在理解波动(二)：以有限计算无限——人工边界条件中，我们说到了如何把“无限远”处理为一种特殊的边界条件。通过这样的处理，一个无阻尼的无界结构等效为一个在边界上具有强阻尼的有限结构。这样，原先由近场流向远场的能量现在流向了等价的边界阻尼器中——已经似乎都很完美。尤其是对于标准波动方程，其ABC既可以用于频域又可以用于时域。但对于弯曲波，还有一个小问题待讨论，那就它的ABC的等价阻尼系数与频率有关，不能处理为阻尼力。边界处的力与位移实际上是“分数阶导数”的关系。这导致它只能直接用于计算频域响应(或只含一个频率成分的激励力)，那么，如何计算无限梁对任意激励的时域响应呢？

无界梁的频域响应计算

欧拉伯努力梁的人工边界可以表示为两个传递函数，分别是：

位移-内力传递函数：

转角-弯矩传递函数：

对于一个无界梁，其近场用有限元建模，如果左边界节点编号是1，右边界是

号，那么可以将其动力学方程写为：

其中，

力向量也按照对应关系排序。

这样，边界传函矩阵写为：

在每个频率下生成边界传递函数矩阵&#