华科流固耦联悬臂输液管大作业微分求积法详细推导与源代码探析.docxVIP

华中科技大学流固耦联悬臂输液管动态特性微分求积法附源代码研究一端固支一端悬臂输液管道的动态特性。管道的流固质量比为，无量纲流速为，最大值为14，悬臂输液管道的无量纲控制方程为为管道的横向位移，为无量纲的轴向位置，范围为0~1，为时间项图1 一端固支一端悬臂输液管模型示意图1、问题分析与方法选择对于一边固支，一边悬臂形式的输液管道，边界条件为当时挠度（横向位移对轴向位移零阶导数）与转角（横向位移对轴向位移一阶导数）均为0，即当时，曲率（横向位移对轴向位移二阶导数）为0，而且曲率的空间变化率（横向位移对轴向位移三阶导数）为0，即最近 30 年，微分求积法（DQM,Differential Quadrature Method）逐渐成为一种求解边值和初值问题的有效方法，逐渐被人们认可，广泛地应用于工程结构的振动分析之中。该方法的本质是把函数在给定网点处的导数值近似地用域上全部网点处函数值的加权和表示对于输液管道模型,传统的作法是采用Galekrin法将控制方程进行离散化,仅提取有限的低阶振型函数作近似处理；DQM则克服了Galekrin法的这种局限性,本质上考虑了所有振型的综合贡献。而且DQM的实施过程避免了Galekrin的大量积分运算,其计算精度和效率较好。所以本文最终选取DQM方法。2、结构离散、计算权系数矩阵微分求积法的核心是将给定网格节点的导数值由由网格域上所有节点的函数值的加权