论文部分内容阅读
本文主要研究由理想可压的流体区域和线弹性的结构区域所构成的流体与结构相耦合系统,对其系统中的若干基本问题进行了分类处理,在此基础上通过粘性项来考虑流体阻尼,而结构阻尼则考虑为Rayleigh阻尼。在经过对用于耦合系统分析的三种基本有限元格式进行讨论和比较,指出了它们各自的优缺点之后,选择位移——位移格式有限元法作了较深入的分析研究。 在本文中,给出了程序直接实现法、罚函数单元法和混合单元法等三种界面对接实现方法,它们能够有效地实现耦合系统中的交界面对接;给出了预条件处理的Lanczos算法,明显地提高了在耦合系统分析过程中可能出现的病态矩阵迭代求解的收敛速度;由罚函数原理引入无旋约束,提出在系统模态中识别流体漩涡模态的差分方法,并用外推法求出系统的真实频率值,还讨论了无旋罚数值在系统分析中的具体选择及其对数值解的影响。 通过对流体位移单元特性的分析,指出一部分常用单元存在仿形能力差以致在求解耦合系统时流体域中出现变形“伪能”,提出了改善单元的插值函数和引进单元降阶/选择技巧等克服这种缺陷的措施,并给出了一种经过改进的具有高仿形能力的流体四节点非协调单元,用位移型有限元法高精度地求得了在流体域中结构振动频率和模态与在平面和回转壳柔性腔内流体晃动频率和模态。 对于流体有限扰动时的非线性耦合系统的有限元分析作了一些探讨,给出了非线性系统的基本方程、变分原理及有限元方程,并具体地导出了三维几何非线性流体单元的切线刚度公式,采用增量形式的Newmark预报——校正逐步积分法计算系统的动力响应,由Newton-Raphson迭代法及其修正方法求解系统的非线性方程组,还列出了不平衡力余量检查和位移迭代增量检查等迭代收敛判断准则。 本文最后列举了由作者编制的FSAP所完成的一部分计算实例,包括静力耦合问题分析、耦合系统的固有摸态分析、线性与非线性的有阻尼耦合系统的动力响应计算等,并同有关的解析解及其它方法解作了详细的对比,结果表明本文中给出的数值计算方法和所编程序都是可靠的和行之有效的。