结构动力方程的精细积分-FFT+方法.pdfVIP

2007 年第九届全国振动理论及应用学术会议论文集 杭州,2007.10.17-19 结构动力方程的精细积分-FFT 方法 张文志，富明慧，刘祚秋 (中山大学应用力学与工程系,广东 广州 510275） 摘 要：将结构动力方程离散，并利用假设始、末时刻位移为已知的方法，将初值问题形式上转化为边 值问题，然后利用快速傅立叶变换（FFT）求解此差分方程，从而得到非齐次结构动力方程的一个数值特 解。将上述非齐次方程的数值特解与精细积分法求得的齐次方程通解相结合，建立了求解结构动力方程的 一种新方法。该方法具有较高的精度和计算效率，算例的数值结果证明了该方法的有效性。 关键词：结构动力方程；精细积分法；差分方程；快速傅立叶变换 On precise integration-FFT method for structural dynamics ZHANG Wen-zhi, Fu Ming-hui, Liu Zuo-qiu (Department of Applied Mechanics and Engineering, Sun Yat-sen University, Guangzhou 510275, China) Abstract : By dispersing the structural dynamic equation and assuming the displacement of the initial and the final time is given, the initial issue form into a boundary problem and then using the fast Fourier transform (FFT) to solve the difference equation, hence a numerical special solution of non-homogeneous structural dynamic equation is obtained. Combining the special solution of the non-homogeneous equation with the homogeneous equation general solution which obtained by the precise integration method (PIM), a new method to solve the structural dynamic equation is established. The method has higher accuracy and computational efficiency. Numerical examples show the validity of the present method. Key words : structural dynamic equation ；precise integration method (PIM) ；difference equation ；fast Fourier transform (FFT) 引言 钟万勰提出的精细积分法[1]用于求解线性定常结构动力方程，可得到逼近精确解的数 值结果。但对于非齐次结构动力方程，会遇到矩阵求逆的问题，由于逆阵有可能不存在或 不稳定，这在一定程度上限制了精细积分法的应用范围。林家浩等[2]将非齐次项化为 Fourier 级数，虽然避免了对方程的系数矩阵求逆，但需要对另外的矩阵求逆。顾元宪等 [3]提出了对结构方程增维的方法，这种方法无需矩阵求逆，但其适用范围有限。汪梦甫[4] 采用将 Pade逼近与指数矩阵运算结合起来，用高斯积分计算积分项的方法，建立了无条件 稳定的精细积分法，但这种方法在精度和计算效率上的矛盾十分突出。 本文利用快速傅立叶变换法求出非齐次结构动力方程的数值特解，并将其与用精细积 分法求得的齐次方程通解结合起来，建立了一种精度好、效率高且普遍适用的结构动力分 析方法。 4 －173