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

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