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

结构系统动力分析可以采用总体结构有限元法，但该方法对于复杂大型结构进行分析存在计算规模大，计算时间长，所用的磁盘空间、计算机系统太庞大，如飞机、车辆、船舶、高层建筑等整体结构。特别是用有限元法进行较高频率振动分析时，要求结构被划分成非常多的单元数以便获得详细的位移和应力特性。这时结构模型的节点自由度可能达到几十万甚至上百万，直接求解如此庞大的模型是很困难。即使能够分析，也要耗费大量机时，效率极低。模态综合法（Component Mode Synthesis）就是在这样的背景下发展起来的一种缩减自由度方法。它可以将大模型化小，先进行各个子结构的模态分析，然后进行模态综合。由于仅采用了各个子结构的低阶模态，因而使所建立整体结构动力模型的自由度数大大降低，而且可以在不同的机器上对各子结构进行模态分析提高计算速度。 模态综合法的基本思想是根据复杂结构的特点将整体结构划分成若干子结构，对各个子结构分别进行模态分析，得到其动力特性。再利用子结构间力平衡条件及位移协调条件将各子结构部分低阶模态特性综合，由此得到整体结构的动力特性。 一.模态综合法基本原理 先将整个结构划分为s个子结构，每个子结构应该是容易分析的，子结构之间的连接尽可能要“弱”些，使子结构之间耦合较小，每个子结构在力学上有较大的独立性。 1.建立各个子结构的模态矩阵，其中；是交界面完全固定时所算出的一部分（参加综合的）固有模态，即主模态矩阵；是依次释放每一边界自由度，使其得到单位位移而产生的静位移分布所组成的约束模态矩阵。它可以由对子结构直接求解一系列静平衡方程（边界自由度发生单位强迫位移）而得到： 下标与分别代表该子结构的内部及界面自由度。结构作自由振动时，内部自由度上作用力为零，而界面自由度依次产生单位位移，则相当于是一个单位矩阵。于是由上式展开第一行可以得到： 从而得到约束模态矩阵，显然的列数等于界面自由度数。 2.用模态矩阵对子结构的刚度阵和质量阵做第一次坐标变换，得到各个子结构的缩减刚度矩阵和质量矩阵，再形式地拼装成不耦合的形式： ；； 其中， 一般来说，选定的前模态组成的模态矩阵，做坐标变化后得到的刚度矩阵，质量矩阵为阶矩阵，远比变换前的刚度矩阵，质量矩阵的阶数小得多。 3.第二次坐标变换，实现子结构的连接。 假如第号子结构与第号子结构有公共界面，则他们的位移向量分别写成： ， 下标B表示公共界面上的约束模态广义位移。由于在界面处接口点的位移必须相等，这就导致与该位移相对应大的约束模态广义位移亦必须相等，也就有： 同样的，对于整个结构可以形式地写出， 式中，是相应于所有子结构主模态的广义位移；是所有独立的约束模态广义位移的集合；是可以用表示的非独立约束模态广义位移，表示其间转换关系。 利用可以进一步将未耦合的总刚度阵及总质量阵缩减（耦合）为： ； 最后，求解下列广义特征值问题： 它的各低阶特征值与原结构低阶特征值的误差，跟子结构剖分是否得当，截取模态是否合理有很大关系。如果想要得到物理坐标下的结构模态向量（振型），需要逆向做两次坐标变换。 下面以悬臂板模态分析为例，分别用商用软件ANSYS和基于有限元算法的MATLAB程序，演示固定界面模态综合法的应用。 二.ANSYS模态综合法求解 ANSYS模态综合法是子结构在动力分析中的应用，其基本过程包括三方面：生成过程、使用过程、扩展过程。ANSYS提供了友好的模态综合法（Component Mode Synthesis）CMS向导（PreprocessorModelingCMS），可以方便的定义超单元和交界面，而且可以对模态综合法分析生成的文件进行管理和组织。 1. 创建超单元 2. 选择CMS求解方法(CMSOPT)：固定界面法或自由界面法。同时确定提取模态数、频率范围等。对于自由界面法还要确定刚体模态。 3. 命名超单元矩阵文件(SEOPT)。 4. 对于考虑阻尼时，输入集中质量矩阵公式。 5. 定义主自由度：在超单元的交界面定义主自由度。 6. 保持数据库：这是必须做的，因为在扩展模态时需要有相同的数据库文件。 7. 求解生成超单元矩阵文件。 2. CMS的使用和扩展过程 CMS的使用和扩展部分与子结构基本相同，但是CMS只支持模态分析。在子结构分析中，需要对整体结构中的每一个子结构进行生成和扩展，然后将各个子结构集合成整个模型。在自由界面模态综合法中，使用MODOPT扩展模态数，应小于每个超单元求解的模态数。若需要扩展更多的模态，需要CMS的超单元重新求解更多的阶数。 以长和宽2的正方形悬臂板为例： 材料参数：弹性模量E=2.1e11 泊松比?=0.3?密度??7.3e3 厚度t=0.05 网格划分：在长度和宽度方向上均划分20个单元（为了与后面程序一致） /PREP7 BLC