《圆管层流》PPT课件.ppt

《《圆管层流》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《圆管层流》PPT课件.ppt（51页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、主讲教师机械工程学院 王开松第五章第五章 圆管层流圆管层流安徽理工大学机械学院安徽理工大学机械学院l流体流动的四种形式：流体流动的四种形式：l流体在固体中间的管中流动或缝隙中流动（机械常见）流体在固体中间的管中流动或缝隙中流动（机械常见）l流体在固体外部的绕流流体在固体外部的绕流l流体在固体一侧的明渠流动流体在固体一侧的明渠流动l流体与固体不接触的孔口出流或射流流体与固体不接触的孔口出流或射流l 主要讨论流动中的能量损失；流动的形式是层流主要讨论流动中的能量损失；流动的形式是层流还是湍流，及其形成的原因；沿程阻力和局部阻力。还是湍流，及其形成的原因；沿程阻力和局部阻力。安徽理工大学机械学院安徽

2、理工大学机械学院l5.1 5.1 雷诺实验和流态判据雷诺实验和流态判据l5.2 5.2 圆管的层流运动圆管的层流运动l5.3 5.3 椭圆管层流椭圆管层流l5.4 5.4 圆管中流体的湍流运动圆管中流体的湍流运动l5.5 5.5 流体运动的两种流动阻力流体运动的两种流动阻力l5.6 5.6 圆管湍流运动的沿程损失圆管湍流运动的沿程损失l5.7 5.7 管流的局部损失管流的局部损失l5.8 5.8 复杂管路计算复杂管路计算安徽理工大学机械学院安徽理工大学机械学院5.1 雷诺实验和流态判据雷诺实验和流态判据1、雷诺实验、雷诺实验 流体的阻力特性直接影响到流体流动时的能量损失，为探索流体摩擦阻力的规

3、律，人们进行了长期研究。1883年，雷诺(Osborne Reynolds)通过大量实验，终于发现了液体在管道中流动时有着两种不同的流动状态，阻力特性也不相同。这种现象可用图5-1所示的雷诺实验装置观测出来。图 5-1 雷诺实验装置1 水龙头；2容器；3水管；4容器；5控制阀安徽理工大学机械学院安徽理工大学机械学院 打开水龙头1，使容器2保持溢流状态；再打开管3的控制阀5使水处于连续滴出状态。为观察流动状态，在容器4中加入染色（如红色）液体。逐步打开阀门5，使管3中流体流速变大，可以观察到：(1)当水平管3中流体流速较小时，染色流体呈一条鲜明的细流（线），非常平稳，染色线与水平管轴线平行或重合

4、（图51(a)）；(2)当管中流速增大到某定值时，染色线开始弯曲颤动，这表明管内流体不再保持安定，不仅有横向脉动速度，而且纵向速度脉动（图5-1(b)）；(3)继续增大流速，染色液体不再保持完整形状而是破裂成杂乱无章、瞬息变化的状态。当使管内流速下降到一定程度时又重复前述状态。这就是著名的雷诺实验。1 水龙头；2容器；3水管；4容器；5控制阀安徽理工大学机械学院安徽理工大学机械学院 vv层流：流体质点无横向脉动，质点互不混杂，层次分明，层流：流体质点无横向脉动，质点互不混杂，层次分明，层流：流体质点无横向脉动，质点互不混杂，层次分明，层流：流体质点无横向脉动，质点互不混杂，层次分明，稳定安详的

5、流动状态。稳定安详的流动状态。稳定安详的流动状态。稳定安详的流动状态。vv紊流紊流紊流紊流 （湍流）：流体质点不仅在轴（横）向而且在纵向（湍流）：流体质点不仅在轴（横）向而且在纵向（湍流）：流体质点不仅在轴（横）向而且在纵向（湍流）：流体质点不仅在轴（横）向而且在纵向(径径径径向向向向)均有不规则脉动速度，流体质点杂乱交错的混沌流动现象。均有不规则脉动速度，流体质点杂乱交错的混沌流动现象。均有不规则脉动速度，流体质点杂乱交错的混沌流动现象。均有不规则脉动速度，流体质点杂乱交错的混沌流动现象。vv实验表明，实验表明，实验表明，实验表明，流体流动具有两种形态，并且可以相互转变流体流动具有两种形态，

6、并且可以相互转变流体流动具有两种形态，并且可以相互转变流体流动具有两种形态，并且可以相互转变。安徽理工大学机械学院安徽理工大学机械学院 2 2、雷诺数、雷诺数流态判别准则流态判别准则 雷诺经过大量实验发现，流体的流动状态与流速、管径、动力粘度和雷诺经过大量实验发现，流体的流动状态与流速、管径、动力粘度和流体密度有关系，可归结为一个无因数流体密度有关系，可归结为一个无因数雷诺数（就是雷诺数（就是惯性力与粘性力惯性力与粘性力的比值的比值）。）。式中：式中：流体密度，流体密度，kg/m3；管内平均流速，管内平均流速，m/s；动力粘度，动力粘度，Pa.s；运动粘度，运动粘度，m2/s；圆管直径，对异形

7、管，则为水力直径，圆管直径，对异形管，则为水力直径，m水力直径可表示为水力直径可表示为 湿周长度（与液面接触的壁面长度）。湿周长度（与液面接触的壁面长度）。式中式中:过流断面面积。过流断面面积。安徽理工大学机械学院安徽理工大学机械学院雷诺数的物理性质及其成因：水流的流态由惯性力和粘性力所起的作用大小所决定。紊流：惯性力起主导作用，质点受约束降低，无规则的脉动增加，产生可见尺度的漩涡，即紊流。雷诺数大；层流：粘性力起主导作用，质点受其约束只能沿运动方向运动，不会偏向，属于层流状态。雷诺数小。安徽理工大学机械学院安徽理工大学机械学院 vv Re 2320Re 2320Re 2320时为紊流；时为紊

8、流；时为紊流；时为紊流；vv 雷诺数表示作用于流体微团的雷诺数表示作用于流体微团的雷诺数表示作用于流体微团的雷诺数表示作用于流体微团的惯性力惯性力惯性力惯性力与粘性力之比。两个几何与粘性力之比。两个几何与粘性力之比。两个几何与粘性力之比。两个几何相似流场的雷诺数相等，则对应微团的惯性力与粘性力之比相等。相似流场的雷诺数相等，则对应微团的惯性力与粘性力之比相等。相似流场的雷诺数相等，则对应微团的惯性力与粘性力之比相等。相似流场的雷诺数相等，则对应微团的惯性力与粘性力之比相等。vv 雷诺数越小，意味着粘性力影响越显著，越大则惯性力影响越雷诺数越小，意味着粘性力影响越显著，越大则惯性力影响越雷诺数越

9、小，意味着粘性力影响越显著，越大则惯性力影响越雷诺数越小，意味着粘性力影响越显著，越大则惯性力影响越显著。雷诺数很小的流动（如润滑膜内的流动），其粘性影响遍及显著。雷诺数很小的流动（如润滑膜内的流动），其粘性影响遍及显著。雷诺数很小的流动（如润滑膜内的流动），其粘性影响遍及显著。雷诺数很小的流动（如润滑膜内的流动），其粘性影响遍及全流场。雷诺数很大的流动（如一般飞行器绕流），其粘性影响仅全流场。雷诺数很大的流动（如一般飞行器绕流），其粘性影响仅全流场。雷诺数很大的流动（如一般飞行器绕流），其粘性影响仅全流场。雷诺数很大的流动（如一般飞行器绕流），其粘性影响仅在物面附近的边界层或尾迹中才是重要的

10、。在涉及粘性影响的流体在物面附近的边界层或尾迹中才是重要的。在涉及粘性影响的流体在物面附近的边界层或尾迹中才是重要的。在涉及粘性影响的流体在物面附近的边界层或尾迹中才是重要的。在涉及粘性影响的流体力学实验中，雷诺数是主要的相似准数力学实验中，雷诺数是主要的相似准数力学实验中，雷诺数是主要的相似准数力学实验中，雷诺数是主要的相似准数圆管中流体的流态判别：圆管中流体的流态判别：安徽理工大学机械学院安徽理工大学机械学院1、圆管层流时的运动微分方程（牛顿力学分析法）、圆管层流时的运动微分方程（牛顿力学分析法）5.2 5.2 圆管中的层流运动圆管中的层流运动 取长为dx 半径为r 的圆柱体，不计质量力和

11、惯性力，仅考虑压力和剪应力，则有根据牛顿粘性定律再考虑到压力的线性分布，有：则有安徽理工大学机械学院安徽理工大学机械学院2、速度分布规律与流量、速度分布规律与流量对上式作不定积分，边界条件：r=R，u=0；则可得定积分常数则图 5-4 圆管层流的速度和剪应力分布圆圆管管层层流的速度分布是以流的速度分布是以圆圆管的管的轴线为轴线为中心中心线线的二次抛物面。的二次抛物面。安徽理工大学机械学院安徽理工大学机械学院2、速度分布规律与流量、速度分布规律与流量图 5-4 圆管层流的速度和剪应力分布 在半径在半径r处处取壁厚取壁厚为为dr的微的微圆环圆环，在，在dr上的速度上的速度为为u，微，微圆环圆环界面

12、上的界面上的微流量微流量dq为为：安徽理工大学机械学院安徽理工大学机械学院3 最大流速与平均流速最大流速与平均流速 由知，r=0时有最大流速 u max，且平均流速4 剪应力分布规律剪应力分布规律 根据牛顿内摩擦定律可求剪应力，它服从线性分布。在r=R也即在管壁上的剪应力最大：安徽理工大学机械学院安徽理工大学机械学院5.4 5.4 圆管中流体的湍流运动圆管中流体的湍流运动1.紊流（湍流）概念及研究方法紊流（湍流）概念及研究方法 紊流特征：紊流是指流体从一种稳定状态向另一种紊流特征：紊流是指流体从一种稳定状态向另一种稳定状态变化过程中的一种无序状态。具体是指稳定状态变化过程中的一种无序状态。具体

13、是指流体流体流流动时，各质点作无规则的运动，各质点间的动时，各质点作无规则的运动，各质点间的惯性力惯性力占主占主要地位，流体作不规则地流动。要地位，流体作不规则地流动。特征有：各层的质点相互掺混；运动要素主特征有：各层的质点相互掺混；运动要素主要是脉动。要是脉动。安徽理工大学机械学院安徽理工大学机械学院5.4 5.4 圆管中流体的湍流运动圆管中流体的湍流运动1.紊流（湍流）概念及研究方法紊流（湍流）概念及研究方法 质点运动轨迹曲折无序，研究紊流的运动规律只能采用时均质点运动轨迹曲折无序，研究紊流的运动规律只能采用时均法，因为时均法的平均值有一定的规律可循。将紊流中各物理量法，因为时均法的平均值

14、有一定的规律可循。将紊流中各物理量的瞬时值看做为：时均值和脉动值之和。的瞬时值看做为：时均值和脉动值之和。将运动要素瞬时值看成将运动要素瞬时值看成时间平均值时间平均值与与脉动值叠加脉动值叠加的方法的方法叫运动要素叫运动要素时均化处理时均化处理。以后就以。以后就以时均值时均值替代瞬时值研究。替代瞬时值研究。紊流的瞬时流速紊流的瞬时流速紊流的瞬时流速紊流的瞬时流速=时均流速时均流速时均流速时均流速+脉动流速脉动流速脉动流速脉动流速安徽理工大学机械学院安徽理工大学机械学院分析时均流速与脉动速度分析时均流速与脉动速度 在时间间隔在时间间隔T内，速度内，速度u随时间进行变化，但时均流速随时间进行变化，但

15、时均流速不变，它只是空间点的函数不变，它只是空间点的函数 瞬时流速瞬时流速u与时均流速的差与时均流速的差值，叫做脉动流速。即值，叫做脉动流速。即安徽理工大学机械学院安徽理工大学机械学院2.紊流层次结构和光滑管概念紊流层次结构和光滑管概念2.1 紊流结构紊流结构 (1)紊流核心区紊流核心区 (2)层流边层层流边层 (3)过渡层过渡层 湍流越剧烈，层流边层越薄。湍流越剧烈，层流边层越薄。层流边层的厚度：层流边层的厚度：阻力系数，与管径、流速和：阻力系数，与管径、流速和管壁的光滑程度有关。下面讨管壁的光滑程度有关。下面讨论光滑管和粗糙管。论光滑管和粗糙管。安徽理工大学机械学院安徽理工大学机械学院层流

16、边层厚度层流边层厚度，与管径、流，与管径、流速和管壁光滑程度有关。速和管壁光滑程度有关。管壁绝对粗糙度管壁绝对粗糙度/d相对粗糙度相对粗糙度 水力水力光滑管光滑管（图（图a)水力水力粗糙管粗糙管(图图b）同一个管，到底是什么管，取同一个管，到底是什么管，取决于层流边层的厚度，或者决于层流边层的厚度，或者雷诺数。雷诺数。安徽理工大学机械学院安徽理工大学机械学院3.紊流的剪应力情况紊流的剪应力情况 紊流的剪应力由两个部分组成：一个是时均流层相对运紊流的剪应力由两个部分组成：一个是时均流层相对运动产生的黏性剪应力动产生的黏性剪应力 ，另外就是上下层质点之间互相掺，另外就是上下层质点之间互相掺混而引起

17、的附加剪应力混而引起的附加剪应力 。当内诺数较小时，紊流运动不很激烈，占主导作用；随当内诺数较小时，紊流运动不很激烈，占主导作用；随着雷诺数不断增大，紊流加剧，附加剪应力增大。着雷诺数不断增大，紊流加剧，附加剪应力增大。安徽理工大学机械学院安徽理工大学机械学院4.紊流基本理紊流基本理论论普朗特混合普朗特混合长长度理度理论论（略略讲讲）两条假设：两条假设：(1)类似于分子的平均自由行程，类似于分子的平均自由行程，紊流流体微团有一个紊流流体微团有一个“混合长度混合长度”。如图，对于某一给定的如图，对于某一给定的y点，点，和和 的流体微团各以时间间隔的流体微团各以时间间隔 dt到达到达y点，在此之前

18、，保持原来的时均速度点，在此之前，保持原来的时均速度 和和不变；一旦达到不变；一旦达到y点，就与该处原流体微团发生碰撞而产生动量点，就与该处原流体微团发生碰撞而产生动量交换。交换。(2)x和和y向的速度涨落（脉动）量向的速度涨落（脉动）量 和和 为同阶量。为同阶量。安徽理工大学机械学院安徽理工大学机械学院5.紊流的速度分布紊流的速度分布 在粘性底层，无流体质点混杂，附加切应力在粘性底层，无流体质点混杂，附加切应力2可略去；可略去；在层流条件下，速度梯度在层流条件下，速度梯度 为常数，则剪应力为常数，则剪应力1为常数，为常数，即即 根据边界条件；根据边界条件；y=0，u=0，可知层流区域的速度分

19、布规，可知层流区域的速度分布规律为律为：安徽理工大学机械学院安徽理工大学机械学院5.紊流的速度分布紊流的速度分布 当紊流充分发展时，流体质点充分混杂，附加切应力当紊流充分发展时，流体质点充分混杂，附加切应力2占主导地位；黏性剪应力占主导地位；黏性剪应力1可忽略不计可忽略不计 尼克拉德塞等人实验证明，对湍流的三个边界层，速度尼克拉德塞等人实验证明，对湍流的三个边界层，速度分布满足经验公式分布满足经验公式5.4-2628 安徽理工大学机械学院安徽理工大学机械学院5.55.5流体运动的两种流动阻力流体运动的两种流动阻力 流体存在粘性，在管路中流动就要受到阻力作用。根据成流体存在粘性，在管路中流动就要

20、受到阻力作用。根据成因不同，分为因不同，分为沿程阻力沿程阻力和和局部阻力局部阻力。沿程阻力沿程阻力hf，是流体在过流断面沿程不变的均匀流中受到的是流体在过流断面沿程不变的均匀流中受到的阻力，主要由流体与管壁面得摩擦引起的，可表示为阻力，主要由流体与管壁面得摩擦引起的，可表示为：安徽理工大学机械学院安徽理工大学机械学院5.55.5流体运动的两种流动阻力流体运动的两种流动阻力 局部阻力局部阻力hj，是流体流过局部装置，因为流体与这些装置内是流体流过局部装置，因为流体与这些装置内部的冲击以及流体质点流速大小和方向的急剧变化引起碰撞引起部的冲击以及流体质点流速大小和方向的急剧变化引起碰撞引起的阻力，可

21、表示为的阻力，可表示为：安徽理工大学机械学院安徽理工大学机械学院5.6 5.6 圆管紊流的沿程损失圆管紊流的沿程损失1.圆管紊流的沿程损失圆管紊流的沿程损失 由于紊流的复杂性，目前工程上确定沿程阻力系数由于紊流的复杂性，目前工程上确定沿程阻力系数的方的方法是：法是：以半经验半理论为基础，整理出来的经验公式；以半经验半理论为基础，整理出来的经验公式；根据实验结果，综合出来的经验公式根据实验结果，综合出来的经验公式对于圆管紊流，沿程损失对于圆管紊流，沿程损失hf 可表示为可表示为 安徽理工大学机械学院安徽理工大学机械学院5.6 圆管紊流的沿程损失绝对粗糙度绝对粗糙度/d相对粗糙度相对粗糙度1.1.

22、尼古拉兹实验尼古拉兹实验（1933-1934）（1 1）实验曲线）实验曲线人工的尼古拉兹粗糙管，表明人工的尼古拉兹粗糙管，表明与与/d。Re的关系的关系安徽理工大学机械学院安徽理工大学机械学院区区-层流区（层流区（ab线，线，lgRe3.3，Re3.6，Re4000）紊流光滑区的紊流光滑区的=f（Re）与相对粗糙度与相对粗糙度/d无关无关 /d小的管子在小的管子在Re较大时，会有稍微的偏离直线较大时，会有稍微的偏离直线 的现象的现象安徽理工大学机械学院安徽理工大学机械学院区区-紊流过渡区紊流过渡区（cd、ef之间的曲线族）之间的曲线族）=f（Re，/d），不同相对粗糙度管有不同的曲线），不同相

23、对粗糙度管有不同的曲线区区-紊流粗糙区（阻力平方区）紊流粗糙区（阻力平方区）（ef右侧水平的直线族）右侧水平的直线族）=f（/d），），不同相对粗糙度管有不同的曲线，但与不同相对粗糙度管有不同的曲线，但与Re无关无关安徽理工大学机械学院安徽理工大学机械学院2.卡门普朗特卡门普朗特(Karman-Prandtl)经验公式经验公式光滑管光滑管 粗糙管粗糙管 上两式有一定理论基础，又有实验资料确定系数，比较上两式有一定理论基础，又有实验资料确定系数，比较精确，缺点是计算不方便。精确，缺点是计算不方便。3.布拉休斯布拉休斯(Blasins)公式公式安徽理工大学机械学院安徽理工大学机械学院4.莫迪莫迪(

24、Moody)图图 实际工业管道难以用相对粗糙度来表征，因而尼古拉兹采用均匀沙粒管道做的实验与自然管道还是有区别的，实际管材内壁的凸凹度，与均匀沙粒的粗糙度有很大的区别，因而其实验结果不能直接应用。1940年美国的莫迪对工业用管进行了大量的实验，得出与Re、/d 的关系图，简便而准确。安徽理工大学机械学院安徽理工大学机械学院4.莫迪莫迪(Moody)图图 安徽理工大学机械学院安徽理工大学机械学院5.7 5.7 管流局部损失管流局部损失 工业管道中，由于设置进出口、弯头、三通、水表等，使得流体流经工业管道中，由于设置进出口、弯头、三通、水表等，使得流体流经此处是，速度、方向均有变化，产生流体质点之

25、间的撞击、分离和漩涡等此处是，速度、方向均有变化，产生流体质点之间的撞击、分离和漩涡等现象，带来附加阻力，造成能里的损失，即为局部损失。现象，带来附加阻力，造成能里的损失，即为局部损失。由于过流断面变化引起的局部损失由于过流断面变化引起的局部损失流动方向的变化引起的局部损失。流动方向的变化引起的局部损失。安徽理工大学机械学院安徽理工大学机械学院5.7 5.7 管流局部损失管流局部损失1管道进口处的损失管道进口处的损失 管道进口处的损失很复杂，经过实验总结得到下面的局管道进口处的损失很复杂，经过实验总结得到下面的局部损失系数部损失系数为：为：局部损失系数局部损失系数=0.5 安徽理工大学机械学院

26、安徽理工大学机械学院5.7 5.7 管流局部损失管流局部损失2.管径突然扩大的局部阻力系数管径突然扩大的局部阻力系数 当管道截面突然扩大时，如图所示，由于流线不能折转，管道截面由A1突然扩大到A2时，管中的流线是逐渐扩散的，因而在管壁的拐角处形成旋涡，由于旋祸要靠主流带动旋转，因此旋涡运动必然要消耗流体的能量。由于细管流速高，粗管流速低，因此从细管流出的流体微团必然要和粗管的流体微团发生碰撞，碰撞和旋祸均会引起流体的能量损失，变成热量耗散。下面用分析法来推导因管道截向突然扩大形成的局部损失。截面11处流体的压强为p1，流速为v1，截面积为A1；截面22处流体的压强为P2，流速为v2，截面面积为

27、A2。列11至22截面的伯努利方程:安徽理工大学机械学院安徽理工大学机械学院5.7 5.7 管流局部损失管流局部损失2.管径突然扩大的局部阻力系数管径突然扩大的局部阻力系数安徽理工大学机械学院安徽理工大学机械学院5.7 5.7 管流局部损失管流局部损失安徽理工大学机械学院安徽理工大学机械学院4出口处的损失出口处的损失安徽理工大学机械学院安徽理工大学机械学院3.线性渐扩和渐缩管局部损失线性渐扩和渐缩管局部损失图 3(a)线性渐扩管图 如图所示，线性扩散或收缩角为如图所示，线性扩散或收缩角为，这时局部损失比较，这时局部损失比较复杂，与复杂，与A1/A2的比值和的比值和 角相关。对于渐扩管，局部阻力

28、角相关。对于渐扩管，局部阻力系数系数 可表示为可表示为 式中 沿程阻力系数。k和扩张角有关的系数。图3(b)线性渐缩管 安徽理工大学机械学院安徽理工大学机械学院用速度用速度U2表示水头损失表示水头损失 :式中式中 管径突然扩大时的局部阻力系数，管径突然扩大时的局部阻力系数，,表 管径突扩时局部损失阻力系数A1/A210.90.80.70.60.50.40.30.20.100.01230.06250.1840.44412.255.441681安徽理工大学机械学院安徽理工大学机械学院局部局部损损失也可用速度失也可用速度 和和 表示为表示为 同样对于管径突然收缩情况（图同样对于管径突然收缩情况（图2

29、），局部损失），局部损失 可表示为可表示为 式中：u 突变面的下游流速。安徽理工大学机械学院安徽理工大学机械学院5收缩管道的局部损失收缩管道的局部损失收缩管道分突然缩小和逐渐缩小。收缩管道分突然缩小和逐渐缩小。突然缩小管道的局部损失发生在细管界面的突然缩小管道的局部损失发生在细管界面的C处附近的漩处附近的漩涡区域，其损失系数涡区域，其损失系数取决于收缩的面积比取决于收缩的面积比A2/A1表 管径突缩时局部损失阻力系数 A1/A20.010.10.20.30.40.50.60.70.80.910.500.470.450.380.340.300.250.200.150.090 逐渐缩小的管道不会出

30、现流线脱离壁面的现象，其局部水逐渐缩小的管道不会出现流线脱离壁面的现象，其局部水头损失头损失取决于收缩的面积比取决于收缩的面积比A2/A1和收缩叫和收缩叫a决定。决定。安徽理工大学机械学院安徽理工大学机械学院6弯管出的水头损失弯管出的水头损失 圆滑弯管和折角弯管的管径不变化，流速大小不变，方向圆滑弯管和折角弯管的管径不变化，流速大小不变，方向改变，会造成能量损失。改变，会造成能量损失。弯管的局部损失主要包括两部分；弯管的局部损失主要包括两部分；旋涡损失；旋涡损失；二次流二次流损失。损失。安徽理工大学机械学院安徽理工大学机械学院7附件处的流动损失附件处的流动损失三通处的损失三通处的损失闸板阀与截止阀处的损失闸板阀与截止阀处的损失液压附件处的损失：如表液压附件处的损失：如表5-4安徽理工大学机械学院安徽理工大学机械学院8安徽理工大学机械学院安徽理工大学机械学院8安徽理工大学机械学院安徽理工大学机械学院8安徽理工大学机械学院安徽理工大学机械学院8安徽理工大学机械学院安徽理工大学机械学院9管路计算管路计算安徽理工大学机械学院安徽理工大学机械学院9管路计算管路计算