结构系统动力分析可以采用总体结构有限元法但该方法对于.docVIP

结构系统动力分析可以采用总体结构有限元法，但该方法对于复杂大型结构进行分析存 在计算规模大，计算时间长，所用的磁盘空间、计算机系统太庞大，如飞机、车辆、船舶、 高层建筑等整体结构。 特别是用有限元法进行较高频率振动分析时， 要求结构被划分成非常 多的单元数以便获得详细的位移和应力特性。 这时结构模型的节点自由度可能达到几十万甚 至上百万， 直接求解如此庞大的模型是很困难。即使能够分析，也要耗费大量机时，效率极 低。模态综合法（ Component Mode Synthesis ）就是在这样的背景下发展起来的一种缩减自 由度方法。 它可以将大模型化小，先进行各个子结构的模态分析， 然后进行模态综合。 由于 仅采用了各个子结构的低阶模态， 因而使所建立整体结构动力模型的自由度数大大降低， 而 且可以在不同的机器上对各子结构进行模态分析提高计算速度。 模态综合法的基本思想是根据复杂结构的特点将整体结构划分成若干子结构，对各个子 结构分别进行模态分析， 得到其动力特性。 再利用子结构间力平衡条件及位移协调条件将各 子结构部分低阶模态特性综合，由此得到整体结构的动力特性。 一.模态综合法基本原理 先将整个结构划分为 s 个子结构，每个子结构应该是容易分析的，子结构之间的连接尽 可能要“弱”些，使子结构之间耦合较小，每个子结构在力学上有较大的独立性。 F B 1.建立各个子结构的模态矩阵 [ ]r (r 1,2, s) ，其中 [ ] [ ] r ； r r F r 是交界面完 全固定时所算出的一部分（参加综合的）固有模态，即主模态矩阵； B r 是依次释放每一边 界自由度， 使其得到单位位移而产生的静位移分布所组成的约束模态矩阵。 它可以由对子结 构直接求解一系列静平衡方程（边界自由度发生单位强迫位移）而得到： k kii u 0 ij i k k u F ji jj j 下标 i 与 j 分别代表该子结构的内部及界面自由度。结构作自由振动时，内部自由度上 作用力为零，而界面自由度依次产生单位位移，则相当于 u 是一个单位矩阵。于是由上式 j 展开第一行可以得到： 1 kii ui kij I 0 ui kii kij I 从而得到约束模态矩阵，显然 B r 的列数等于界面自由度数。 B r ui I 2.用模态矩阵 [ i ]对子结构的刚度阵 [k ]和质 [mi ] 量阵做第一次坐标变换，得到各个子 i 结构的缩减刚度矩阵 [k ]和质量矩阵 [mi ] ，再形式地拼装成不耦合的形式： i [ ] 0 0 [ 1 k ] 0 0 1 [m ] 0 1 0 [ ] 0 [ ] 0 i ； [ k] 0 [k ] 0 ；[m] i 0 [m i ] 0 0 0 [ ] 0 0 s [k ] 0 s 0 [m s ] T T 其中 [ki ] [ ] [ki ][ i ] ， [mi ] [ ] [mi ][ i ] i i 一般来说，选定的前 N 模态组成的模态矩阵，做坐标变化后得到的刚度矩阵，质量矩阵 为 N N 阶矩阵，远比变换前的刚度矩阵，质量矩阵的阶数小得多。 3.第二次坐标变换，实现子结构的连接。 假如第 r 号子结构与第 q号子结构有公共界面，则他们的位移向量分别写成： p rF [ p ] ， r p rB [ p qF p ] q p qB 下标 B 表示公共界面上的约束模态广义位移。由于在界面处接口点的位移必须相等，这 就导致与该位移相对应大的约束模态广义位移亦必须相等，也就有： [ prB ] [ pqB ] [ p r p q ] p rF p rB p qF p qB I 0 0 0 0 I 0 I 0 0 I 0 p rF p rB p qF 同样的，对于整个结构可以形式地写出， p F p B p B [ I 0 0 ] [ [ 0 I ] ] [ p F p B [ p] ] [ ] [ p] 式中， [ pF ] 是相应于所有子结构主模态的广义位移； [ pB ] 是所有独立的约束模态广义 位移的集合； [ p ] 是可以用 [ pB ] 表示的非独立约束模态广义位移， [ ] 表示其间转换关系。 B 利用 [ ] 可以进一步将未耦合的总刚度阵 [k] 及总质量阵 [ m]缩减（耦合）为： T [K ] [ ] [ k][ ] T ；[ M ] [ ] [ m][ ] 最后，求解下列广义特征值问题： [ K ] [ x] [ ][ ]2 M x 2 M x 它的各低阶特征值与原结构低阶特征值的误差，跟子结构剖分是否得当，截取模态是否 合理有很大关系。如果想要得到物理坐标下的结构模态向量（振型） ，需要逆向做两次坐标 变换。 下面以悬臂板模态分析为例， 分别用商用软件 ANSYS 和基于有限元算法的