065基于NASTRAN的贮箱结构与液体耦合.docVIP

基于NASTRAN的贮箱结构与液体耦合 模态分析方法 杨琼梁 唐国安? （复旦大学力学与工程科学系，上海，200433） 摘要：在理想流体及小幅度晃动的假设条件下，表征液体晃动的压力函数的控制方程可以线性化，经有限元方法离散后，便能导出结构振动与液体晃动耦合问题的二阶线性微分方程组，其中系数矩阵是非对称的。采用静态缩聚技术，方程中的液体非自由面上的节点压力可被凝聚掉。进一步用模态坐标表示液体自由面的压力向量，得到了具有对称系数矩阵的耦合方程组。最终，液体是以超单元形式出现在方程中，从而可以用NASTRAN程序求解结构与液体的耦合动力学问题。 关键词：流固耦合；模态分析；有限元 Modal Analysis to the Tank Structure Coupled with liquid Based on NASTRAN YANG Qiong-liang, TANG Guo-an (Department of Mechanics and Engineering Science, Fudan University, Shanghai 200433) Abstract: Based on hypothesis of ideal fluid and small shaking amplitude, the equations of the pressure functions can be linearized. The second order differential equations of the structures coupled with fluid can be obtained through the FEM methods, in which the coefficient matrix is asymmetric. The node pressure on the non- free surface can be reduced using the static condensation technique. The pressure vector on the free surface is expressed in the modal coordinate which makes the coupled equations symmetrical. The fluid is the super-element in equations which can be solved by the NASTRAN. 1 引言 内的压力满足Laplace方程 当满足小幅晃动条件时，自由面上液体压力满足 在液体与结构接触面上液体压力满足 在边界条件和中，表示液体表面的法向单位矢量(指向液体外部)，为结构沿方向的加速度，为液体密度，为重力加速度。 用有限元方法对液体压力进行插值，将压力函数表示成 其中是液体有限元模型的节点压力向量，为插值型函数。根据插值公式式可以推导出液体压力梯度的表达式 于是，方程和,边界条件可以写成离散形式 其中 结构有限元模型的控制方程为 式中为结构的节点位移向量，和分别是结构的质量和刚度矩阵，为结构除液体压力之外所受的外力向量。 将和写成联立形式，就得到结构振动与液体晃动的耦合方程 3 液体节点压力自由度缩聚 将液体的节点压力向量划分成自由面压力和非自由面压力两部分，即，液体质量矩阵和刚度矩阵也相应地写成分块形式。从公式可知，计算矩阵的面积分仅限于自由面，因此矩阵只有与自由面相关的元素非零，因此 将分块形式的和代入方程，即有 其中。从方程的第二式可以确定与约束关系 将上式代回到从方程的第一式，得到仅由自由面节点压力表示的液体控制方程 其中 再将代入方程，则可得到结构的控制方程 其中 将和写成联立形式后得到 4 液体压力主坐标表示的耦合方程 求解广义特征值问题 可以得到个特征值和对应的线性无关的特征向量。若液体自由面上有个节点，则线性空间中向量一定可以有这一组特征向量线性表示： 式中，。 实际上广义特征值问题就是当结构为刚性时液体自由晃动的特征方程，特征值就是刚壁容器内液体的自由晃动频率，且。主座标表示液体自由面的均匀压力，即自由面平直升降。特征向量满足正交性条件 其中 将主坐标表示的液体自由面压力向量代入到方程中，利用正交性条件可以得到 式中。 在零初始条件下，从方程的第一式可解出 这个式子表示了液体自由面平直升降完全取决于结构的变形，是液体不可压缩假设下的必然结果。方程第二式可得 再将和代入到方程中，又可得到 其中 用同乘以方程第二式两端，并于方程联立后最