利用中心差分法计算习题.docVIP

4.2 利用中心差分法计算习题4-1所给系统的动力响应，注意考察时间积分步长对结果的影响。模型如下图所示： （一）运动方程为： （二）中心差分法（central difference）简介: 把位移函数按泰勒级数展开，得前差分公式： （1） 同样可得后差分公式： （2） 把以上两式相减和相加后分别可得： （3） （4） 在以上两式中忽略高阶小量，可得到用、、时刻的位移、、近似表示的t 时刻的速度和加速度： （5） （6） 为了求解时刻的位移，将t时刻的速度和加速度近似式（5）、（6）代入时刻t的运动方程 （7） 中得： （8） 其中： （9） （10） 如已经求得、时刻的位移和，则从式（8）式中可解的时刻的位移。由式（8）可见，只给定初值和无法得到，还必须确定，即该方法存在如何起步的问题。在后差分（2）式中取得到： 由上可构造中心差分法的计算算法，不再赘述，可参见文献[1]。下面再简要说一下中心差分法的稳定性： 对于系统中的步长必须满足 称为临界时间步长 对于无阻尼系统： 对于有阻尼系统： （为系统的固有最小的周期） 由于中心查分法的时间步长必须小于临界步长，因此它是条件稳定的。在用有限元法进行动力分析时，系统自由度数很高，很小，因此必须选的很小。另外在中心差分法中要求质量矩阵中的全部对角元素，如果有对角元素等于0或者接近于0，也将等于0 或者接近于0，中心差分法将不能使用。因此在用有限元进行动力分析时，如使用中心差分法积分，则不宜是个别的单元的尺寸比其它单元的尺寸小很多，否则会减少临界时间步长，因而大幅度增加计算量。 （三）系统的动响应 源程序见附件（4.2eig.m）（为便于进行比较各物理参数设置同4.1题） 首先给出运用振型叠加法算出的动响应时程图 图1.1 2．该系统无阻尼 Ⅰ.控制时间步长，得到的动响应图形如下图所示： 图2.1 图2.2 Ⅱ. 控制时间步长，得到的动响应图形如下图所示 图2.3 图2.4 图2.5 （四）结果分析 1.通过比较图1.1 和图2.1 图2.2 可知控制时间步长时，得到的动响应图形是一致的。 2.通过比较图2.1 图2.1 和图2.3、图2.4、图2.5可知 控制时间步长，得到的动响应图形不具有稳定性：由图2.3，可知此时时间步长只是大于系统的临界步长，对应的质量块3的响应曲线已具有有所不同（与图2.1比），不稳定；由图2.4知，此时只有质量块1的响应曲线大致与图2.1中的一致，但其余两质量块的响应曲线已明显发散；由图2.5知， 此时各质量块的响应曲线均发散。 （其中，为系统的相应的振动周期） 参考文献; [1] 张雄. 《计算动力学》 清华大学出版社：P147-151