高等学校教学用书 流体力学与流体机械 习题参考答案 主讲：陈庆光 中国矿业大学出版社 张景松编.流体力学与流体机械, 徐州：中国矿业大学出版社，2001.6（2005.1重印）
删掉的题目：1-14、2-6、2-9、2-11、2-17、3-10、3-19、4-5、4-13 《流体力学与流体机械之流体力学》 第一章 流体及其物理性质 1-8 1.53m 的容器中装满了油。已知油的重量为12591N 。求油的重度γ和密度ρ。 解：312591 856.5kg/m 9.8 1.5 m V ρ= ==?；38394N/m g γρ== 1-11 面积20.5m A =的平板水平放在厚度10mm h =的油膜上。用 4.8N F =的水平力拉它以0.8m/s U =速度移动（图1-6）。若油的密度3856kg/m ρ=。求油的动力粘度和运动粘度。 解：29.6N/m F A τ= =，U h τμ=， 所以，0.12Pa s h U τμ==g ，42/0.12/856 1.410m /s νμρ-===? 1-12 重量20N G =、面积20.12m A =的平板置于斜面上。其间充满粘度0.65Pa s μ=g 的油液（图1-7）。当油液厚度8mm h =时。问匀速下滑时平板的速度是多少。
解：sin 20 6.84F G N ==o ，57Pa s F A τ==g ， 因为U h τμ =，所以570.0080.7m/s 0.65h U τμ?=== 1-13 直径50mm d =的轴颈同心地在50.1mm D =的轴承中转动（图1-8）。间隙中润滑油的粘度0.45Pa s μ=g 。当转速950r/min n =时，求因油膜摩擦而附加的阻力矩M 。 解：将接触面沿圆柱展开，可得接触面的面积为： 20.050.10.016m A dL ππ==??= 接触面上的相对速度为：2 2.49m/s 2260d d n u πω=== 接触面间的距离为：0.05mm 2D d δ-== 接触面之间的作用力：358.44N du F A A dy u δ μμ=== 则油膜的附加阻力矩为：8.9N m 2 d M F ==g 1-14 直径为D 的圆盘水平地放在厚度为h 的油膜上。当驱动圆盘以转速n 旋转时，试证明油的动力粘度μ与驱动力矩M 的关系为： 24 960hM nD μπ= 证明：26030n n ππω= = ，30 nr v r πω== 2dA rdr π=，2215v nr dr dF dA h h μπμ== ，2315nr dr dM dFr h μπ== /2 2324 15960D nr dr nD h M h μπμπ= =? 所以：24 960hM nD μπ=
算法分析与设计题目：最大流算法 院系：软件工程 班级：软件11-2班 姓名：慕永利 学号：23 号
目录 1算法提出背景............................................................................................................................- 3 - 2 问题实例及解决.......................................................................................................................- 3 - 3算法论述....................................................................................................................................- 4 - 3.1、可行流..........................................................................................................................- 4 - 3.2 最大流..........................................................................................................................- 5 - 3.3最大流算法.....................................................................................................................- 6 - 3.3.1 增广路径.......................................................................................................- 6 - 3.3.2沿增广路径增广..................................................................................................- 7 - 3.3.3样例：..................................................................................................................- 8 - 3.3.4定理：............................................................................................................... - 13 - 3.3.5算法的实现：................................................................................................... - 13 - 3.3.6 优化.................................................................................................................. - 16 - 4算法应用................................................................................................................................. - 18 -
工程流体力学与流体机械(二) (总分：54.00，做题时间：90分钟) 一、{{B}}单项选择题{{/B}}(总题数：54，分数：54.00) 1.水力最优梯形断面的水力半径为水深的{{U}} {{/U}}。 ? A. 2倍 ? B. 1倍 ? C. 一半 ? D. 0.3倍 （分数：1.00） A. B. C. √ D. 解析：梯形断面的水力半径为[*]，水力最优条件为： [*]，代入得：[*]。2.运动黏度的量纲为{{U}} {{/U}}。 ? A. [L1T2] ? B. [L-1T-2] ? C. [L-2T-1] ? D. [L2T-1] （分数：1.00） A. B. C. D. √ 解析： 3.超音速一元等熵气流在渐扩管道中随截面积增大，温度T{{U}} {{/U}}。 ? A. 减小 ? B. 增大 ? C. 不变 ? D. 不确定 （分数：1.00） A. √ B. C. D. 解析： 4.在圆管流动中，层流时切应力在断面分布为{{U}} {{/U}}。
? A. 常数 ? B. 管轴处为零，与半径成正比 ? C. 管壁处为零，向管轴线形增大 ? D. 按抛物线分布 （分数：1.00） A. B. √ C. D. 解析：参考[例4]。 5.气蚀产生的主要原因是{{U}} {{/U}}。 ? A. 压力降低 ? B. 压力升高 ? C. 温度升高 ? D. 速度较大 （分数：1.00） A. √ B. C. D. 解析：气蚀产生的主要原因是在较低压力下，水在常温下汽化及原溶于水中气体大量析出，又经加压后气泡破裂。这一过程伴随局部高温，逸出氧气等对金属有腐蚀的过程。 6.可压缩气体一元恒定流绝热过程的能量方程是{{U}} {{/U}}。 （分数：1.00） A. B. C. D. √ 解析：可压缩气体一元恒定流的能量方程是：[*]，将[*]代入可得到结论(D)。 7.管嘴的流速系数φn为( )。 ? A. 0.62 ? B. 1.32 ? C. 0.75 ? D. 0.82 （分数：1.00） A. B. C. D. √ 解析：见[例9]。
1459:Power Network Time Limit: 2000MS Memory Limit: 32768K Total Submissions: 17697 Accepted: 9349 Description A power network consists of nodes (power stations, consumers and dispatchers) connected by power transport lines. A node u may be supplied with an amount s(u) >= 0 of power, may produce an amount 0 <= p(u) <= p max(u) of power, may consume an amount 0 <= c(u) <= min(s(u),c max(u)) of power, and may deliver an amount d(u)=s(u)+p(u)-c(u) of power. The following restrictions apply: c(u)=0 for any power station, p(u)=0 for any consumer, and p(u)=c(u)=0 for any dispatcher. There is at most one power transport line (u,v) from a node u to a node v in the net; it transports an amount 0 <= l(u,v) <= l max(u,v) of power delivered by u to v. Let Con=Σu c(u) be the power consumed in the net. The problem is to compute the maximum value of Con. An example is in figure 1. The label x/y of power station u shows that p(u)=x and p max(u)=y. The label x/y of consumer u shows that c(u)=x and c max(u)=y. The label x/y of power transport line (u,v) shows that l(u,v)=x and l max(u,v)=y. The power consumed is Con=6. Notice that there are other possible states of the network but the value of Con cannot exceed 6. Input There are several data sets in the input. Each data set encodes a power network. It starts with four integers: 0 <= n <= 100 (nodes), 0 <= np <= n (power stations), 0 <= nc <= n (consumers), and 0 <= m <= n^2 (power transport lines). Follow m data triplets (u,v)z, where u and v are node identifiers (starting from 0) and 0 <= z <= 1000 is the value of l max(u,v). Follow np doublets (u)z, where u is the identifier of a power station and 0 <= z <= 10000 is the value of p max(u). The data set ends with nc doublets (u)z, where u is the identifier of a consumer and 0 <= z <= 10000 is the value of
流体力学与流体机械 模拟题A 一 、 填空题 1、 往复式压缩机容积系数的含义是（ ），表达式（ ）； 2、大、中型压缩机级数的选择，一般以（ ）为原则； 小型压缩机级数选择多取决于每级允许的（ ）； 3、对于理想气体，多变压缩理论循环指示功为L=（ ）； 4、压缩机正常运转时，产生的作用力有四类，一类是质量力，另外三类 分别为（ ）；（ ）；（ ）。 5、离心泵扬程的概念是（ ）； 6、影响离心泵汽蚀的因素至少应包括以下内容（ ），（ ），（ ），（ ）。 7、离心泵相似包含（ ），（ ），（ ），（ ）； 8、离心泵不稳定工况产生的原因有（ ），（ ）。 9、对同一台压缩机，在等温、绝热和多变循环中耗功最小的为（ ）。 二、一个充满水的密闭容器以等角速度ω绕一水平轴旋转，同时考虑重力的影响，试 证明其等压面为圆柱面，且等压面的中心轴线比容器的转动轴线高2ωg 。 三、流场中任意点的速度为：2 2y x ky u x +- =，2 2y x kx u y += ，0=z u 。求流场流线方程和迹线方程。 四、沿喷嘴中心线速度为2 12??? ??+=L x t v ，求流体的位变和时变加速度。 五、判断题：
1、实际压缩循环指压缩机实际运转的气体循环，而理论循环则指在设计阶段理论上的论证循环。 2、实际循环功可用示功图表示，而理论循环则不能。 3、温度系数是压缩机冷却系统的换热系数。 4、在使用铅基润滑油润滑的压缩机中，排气温度不亦太高，否则会使润滑油粘度降低，造成积炭现象。 5、水析是指压缩机工艺管道中，混在原料气中的水分在油水分离器中分离出来。 6、活塞是利用阻塞和节流的作用以达到密封原理的。 7、离心泵串联是为了增大泵的流量，而并联则提高泵的扬程。 8、比转数相同，则离心泵一定相似。 9、离心泵叶轮由于切割量较小，则切割后的泵与原泵相似。 10、离心泵进口速度三角形决定了流体在进口处的能量。 六、简述题 1、压缩机复算性热力计算的步骤有哪些？ 2、简要说明离心泵H —Q 线的作图步骤 流体力学与流体机械A 卷参考答案 二、证： 对旋转圆筒建立oxyz 坐标系，取y 轴为旋转轴。 （1’） 单位质量的重力分量为0=x g ，0=y g ，g g z -=。 （3’） 单位质量流体所受的惯性力为x a x 2ω=-，0=y a ，z a z 2ω=-。 （4’）
最大流算法 clc,clear,M=1000; c(1,2)=3;c(1,4)=3; c(2,3)=1;c(2,4)=20; c(3,6)=3; c(4,5)=10; c(5,1)=4;c(5,3)=2;c(5,6)=13; n=length(u); list=[]; maxf=zeros(1:n);maxf(n)=1; while maxf(n)>0 maxf=zeros(1,n);pred=zeros(1,n); list=1;record=list;maxf(1)=M; while (~isempty(list))&(maxf(n)==0) flag=list(1);list(1)=[]; index1=(find(u(flag,:)~=0)); label1=index1(find(u(flag,index1)... -f(flag,index1)~=0)); label1=setdiff(label1,record); list=union(list,label1); pred(label1(find(pred(label1)==0)))=flag; maxf(label1)=min(maxf(flag),u(flag,label1)... -f(flag,label1)); record=union(record,label1); label2=find(f(:,flag)~=0); label2=label2'; label2=setdiff(label2,record); list=union(list,label2); pred(label2(find(pred(label2)==0)))=-flag; maxf(label2)=min(maxf(flag),f(label2,flag)); record=union(record,label2); end if maxf(n)>0 v2=n; v1=pred(v2); while v2~=1 if v1>0
一、判断题： 1、温度对液体和气体粘性的影响是截然不同的。当温度升高时，液体的粘性增加；而气体的粘性降低。 A、对 B、错 2、在平衡的流体中，通过任意一点的等压面，必与该点所受的质量力相垂直。 A、对 B、错 3、连续性方程的实质是质量守恒定律在流体力学中的表现形式。 A、对 B、错 4、动能修正系数表示截面上以平均流速表示的单位重量流体的动能与实际的平均单位重量流体的动能之比。 A、对 B、错 5、沿程阻力一般发生在沿流程边界形状（过流断面）变化不大的区域，一般在急变流区域。 A、对 B、错 6、在一层流运动的管流中，管流中心处的最大流速是其断面上的平均流速的二倍。 A、对 B、错 7、对于几条管段组成的串联管路，它的水头损失等于这些管路中各管路的水头损失之和。 A、对 B、错 8、按能量传递的方式划分，流体机械可划分为离心式、轴流式与混流式。 A、对 B、错 9、扬程是指单位体积流体通过水泵时，自水泵获得的能量。 A、对 B、错 10、通风机扩散器的作用是减速增压，使气流中一部分动压转变为静压，减少损失，以提高风机的静压与静压效率。 A、对 B、错 二、问答题： 1、流体静压强的两个重要特性是什么?
2、何为过流断面、湿周和水力半径？ 3、何为水力光滑管? 4、什么是物理方程量纲的和谐性原理？ 5、何为经济流速与经济管径？ 6、水泵(合理工况的)正常工作条件是什么？ 三、计算题： 1、用如图的双U 形管来测定比重比水小的液体 的比重。试用液柱高差来确定未知液体的相对密度（比 重）。（管中的水是在标准大气压下，C 04的纯水。忽 略U 形管中气柱的影响。） 题1图
删掉的题目：1-14、2-6、2-9、2-11、2-17、3-10、3-19、4-5、4-13 《流体力学与流体机械之流体力学》 第一章 流体及其物理性质 1-8 1.53m 的容器中装满了油。已知油的重量为12591N 。求油的重度γ和密度ρ。 解：312591 856.5kg/m 9.8 1.5 m V ρ= ==?；38394N/m g γρ== 1-11 面积20.5m A =的平板水平放在厚度10mm h =的油膜上。用 4.8N F =的水平力拉它以0.8m/s U =速度移动（图1-6）。若油的密度3856kg/m ρ=。求油的动力粘度和运动粘度。 解：29.6N/m F A τ= =，U h τμ=， 所以，0.12Pa s h U τμ==，42/0.12/856 1.410m /s νμρ-===? 1-12 重量20N G =、面积2 0.12m A =的平板置于斜面上。其间充满粘度 0.65Pa s μ=的油液（图1-7）。当油液厚度8mm h =时。问匀速下滑时平板的速度是多少。
解：sin 20 6.84F G N ==，57Pa s F A τ==， 因为U h τμ =，所以570.0080.7m/s 0.65 h U τμ?=== 1-13 直径50mm d =的轴颈同心地在50.1mm D =的轴承中转动（图1-8）。间隙中润滑油的粘度0.45Pa s μ=。当转速950r/min n =时，求因油膜摩擦而附加的阻力矩M 。 解：将接触面沿圆柱展开，可得接触面的面积为： 20.050.10.016m A dL ππ==??= 接触面上的相对速度为：2 2.49m/s 2260d d n u πω= == 接触面间的距离为：0.05mm 2 D d δ-== 接触面之间的作用力：358.44N du F A A dy u δ μμ=== 则油膜的附加阻力矩为：8.9N m 2 d M F == 1-14 直径为D 的圆盘水平地放在厚度为h 的油膜上。当驱动圆盘以转速n 旋转时，试证明油的动力粘度μ与驱动力矩M 的关系为： 24 960hM nD μπ= 证明：26030n n ππω= =，30 nr v r πω== 2dA rdr π=，2215v nr dr dF dA h h μπμ==，2315nr dr dM dFr h μπ==
~ 高等学校教学用书 》 流体力学与流体机械习题参考答案 : 主讲：张明辉 中国矿业大学出版社 &
张景松编.流体力学与流体机械, 徐州：中国矿业大学出版社，（重印） 删掉的题目：1-14、2-6、2-9、2-11、2-17、3-10、3-19、4-5、4-13 《流体力学与流体机械之流体力学》 第一章 流体及其物理性质 1-8 3m 的容器中装满了油。已知油的重量为12591N 。求油的重度γ和密度ρ。 解：312591 856.5kg/m 9.8 1.5 m V ρ= ==?；38394N/m g γρ== 1-11 面积20.5m A =的平板水平放在厚度10mm h =的油膜上。用 4.8N F =的水平力拉它以0.8m/s U =速度移动（图1-6）。若油的密度3856kg/m ρ=。求油的动力粘度和运动粘度。 . 解：29.6N/m F A τ= =，U h τμ=， 所以，0.12Pa s h U τμ==，42/0.12/856 1.410m /s νμρ-===? 1-12 重量20N G =、面积20.12m A =的平板置于斜面上。其间充满粘度0.65Pa s μ=的油液（图1-7）。当油液厚度8mm h =时。问匀速下滑时平板的速度是多少。
解：sin 20 6.84F G N ==，57Pa s F A τ==， 因为U h τμ =，所以570.0080.7m/s 0.65h U τμ?=== 1-13 直径50mm d =的轴颈同心地在50.1mm D =的轴承中转动（图1-8）。间隙中润滑油的粘度0.45Pa s μ=。当转速950r/min n =时，求因油膜摩擦而附加的阻力矩M 。 / 解：将接触面沿圆柱展开，可得接触面的面积为： 20.050.10.016m A dL ππ==??= 接触面上的相对速度为：2 2.49m/s 2260d d n u πω=== 接触面间的距离为：0.05mm 2D d δ-== 接触面之间的作用力：358.44N du F A A dy u δ μμ=== 则油膜的附加阻力矩为：8.9N m 2 d M F == 1-14 直径为D 的圆盘水平地放在厚度为h 的油膜上。当驱动圆盘以转速n 旋转时，试证明油的动力粘度μ与驱动力矩M 的关系为： — 24 960hM nD μπ= 证明：26030n n ππω==，30 nr v r πω== 2dA rdr π=，2215v nr dr dF dA h h μπμ== ，2315nr dr dM dFr h μπ== /2 2324 15960D nr dr nD h M h μπμπ= =?
《流体力学与流体机械》复习考试资料 仅供部学习交流使用安全131班编制 绪论： 1.流体力学是以研究流体（包括液体和气体）为研究对象，研究其平衡和运动基本规律的科学。主要研究流体在平衡和运动时的压力分布、速度分布、与固体之间的相互作用以及流动过程中的能量损失。 2.流体力学的主要研究方法：实验研究、理论分析、数值计算。第一章流体及其物理性质 1.流体：在任何微小剪切力下能产生连续变形的物质即为流体。 主要特征：流动性 2.连续介质假说：质点（而不是分子）是组成宏观流体的最小基元，质点与质点之间没有间隙其物理性质各向同性，且在空间和时间上具有连续性。 3.流体的粘性 （1）流体产生粘性的原因：流体的聚力；动量交换；流体分子和固体壁面之间的附着力。 （2）流层之间的摩擦力:带动力和阻力（一对大小相等、方向相反的作用力） （3）流体摩擦切应力:τ=μ·（du/dy) (N/m2) τ=F/A=μ·U/h (N/m2) （4）相对运动的结果使流体产生剪切变形。流体的粘性就是阻止发生剪切变形的一种特性，而摩擦力则是粘性的动力表现。
（5）粘性的度量：动力粘度μ＝τ/（du/dy) (pa·s) 运动粘度ν＝μ/ρ (m2/s) 温度升高时，流体的粘性降低，气体的粘性增加。 4.课后习题答案 第二章流体静力学 1.作用在流体上的力
（1）表面力：作用在被研究流体的表面上，其大小与被作用的面积成正比，如法向压力和切向摩阻力。（平衡流体不存在表面切向力，只有表面法向力） （2）质量力：作用在被研究流体的每个质点上，其大小与被研究流体的质量成正比，如重力和惯性力。质量力常用单位质量力表示，所谓单位质量力，是指作用在单位质量流体上的质量力。 2.流体静压力及其特性 流体处于平衡状态时，表面力只有压力，称其为静压力，单位面积上作用的静压力称为静压强。 静压力有两个重要特性： ①静压力垂直于作用面，并沿着作用面法线方向； ②平衡流体中任何一点的静压力大小与其作用面的方位无关，其值均相等。 3.流体平衡微分方程式（压力差公式） dp=ρ(Xdx+Ydy+Zdz) 4.等压面：平衡流体中压力相等的点所组成的平面或曲面称为等压面。等压面的两个性质：（1）平衡流体中，任一点的等压面恒与质量力正交；（2）当两种互不相混的液体处于平衡时，它们的分界面必为等压面。 5.重力作用下流体静压力的分布规律 （1）静压强分布规律
从一道题目的解法试谈网络流的构造与算法 福建师大附中江鹏 1. 引论 A. 对网络流算法的认识 网络流算法是一种高效实用的算法，相对于其它图论算法来说，模型更加复杂，编程复杂度也更高，但是它综合了图论中的其它一些算法（如最短路径），因而适用范围也更广，经常能够很好地解决一些搜索与动态规划无法解决的，看似NP的问题。 B. 具体问题的应用 网络流在具体问题中的应用，最具挑战性的部分是模型的构造。这没用现成的模式可以套用，需要对各种网络流的性质了如指掌（比如点有容量、容量有上下限、多重边等等），并且归纳总结一些经验，发挥我们的创造性。
2. 例题分析 【问题1】项目发展规划（Develop） Macrosoft?公司准备制定一份未来的发展规划。公司各部门提出的发展项目汇总成了一张规划表，该表包含了许多项目。对于每个项目，规划表中都给出了它所需的投资或预计的盈利。由于某些项目的实施必须依赖于其它项目的开发成果，所以如果要实施这个项目的话，它所依赖的项目也是必不可少的。现在请你担任Macrosoft?公司的总裁，从这些项目中挑选出一部分，使你的公司获得最大的净利润。 ●输入 输入文件包括项目的数量N，每个项目的预算Ci和它所依赖的项目集合Pi。格式如下：第1行是N； 接下来的第i行每行表示第i个项目的信息。每行的第一个数是Ci，正数表示盈利，负数表示投资。剩下的数是项目i所依赖的项目的编号。 每行相邻的两个数之间用一个或多个空格隔开。 ●输出 第1行是公司的最大净利润。接着是获得最大净利润的项目选择方案。若有多个方案，则输出挑选项目最少的一个方案。每行一个数，表示选择的项目的编号，所有项目按从小到大的顺序输出。 ●数据限制 0≤Ｎ≤1000 -1000000≤Ci≤1000000 ●输入输出范例
第一部分：流体力学 1.力按物理性质的不同分类：重力、摩擦力、惯性力、弹性力、表面力等。 2. 力按作用方式分：质量力和面积力。 3 质量力是指作用于隔离体每一流体质点上的力，它的大小与质量成正比 4 最常见的质量力有：重力、惯性力。 5 比较重力场（质量力只有重力）中，水和水银所受的单位质量力f水和f水银的大小？ A. f水
《流体力学与流体机械》复习考试资料 仅供内部学习交流使用安全131班编制 绪论： 1.流体力学是以研究流体（包括液体和气体）为研究对象，研究其平衡和运动基本规律的科学。主要研究流体在平衡和运动时的压力分布、速度分布、与固体之间的相互作用以及流动过程中的能量损失。 2.流体力学的主要研究方法：实验研究、理论分析、数值计算。第一章流体及其物理性质 1.流体：在任何微小剪切力下能产生连续变形的物质即为流体。 主要特征：流动性 2.连续介质假说：质点（而不是分子）是组成宏观流体的最小基元，质点与质点之间没有间隙其物理性质各向同性，且在空间和时间上具有连续性。 3.流体的粘性 （1）流体产生粘性的原因：流体的内聚力；动量交换；流体分子和固体壁面之间的附着力。 （2）流层之间的内摩擦力:带动力和阻力（一对大小相等、方向相反的作用力） （3）流体内摩擦切应力:τ=μ·（du/dy) (N/m2) τ=F/A=μ·U/h (N/m2) （4）相对运动的结果使流体产生剪切变形。流体的粘性就是阻止发生剪切变形的一种特性，而内摩擦力则是粘性的动力表现。
（5）粘性的度量：动力粘度μ＝τ/（du/dy) (pa·s) 运动粘度ν＝μ/ρ (m2/s) 温度升高时，流体的粘性降低，气体的粘性增加。 4.课后习题答案
第二章流体静力学 1.作用在流体上的力 （1）表面力：作用在被研究流体的表面上，其大小与被作用的面积成正比，如法向压力和切向摩阻力。（平衡流体不存在表面切向力，只有表面法向力） （2）质量力：作用在被研究流体的每个质点上，其大小与被研究流体的质量成正比，如重力和惯性力。质量力常用单位质量力表示，所谓单位质量力，是指作用在单位质量流体上的质量力。 2.流体静压力及其特性 流体处于平衡状态时，表面力只有压力，称其为静压力，单位面积上作用的静压力称为静压强。 静压力有两个重要特性： ①静压力垂直于作用面，并沿着作用面内法线方向； ②平衡流体中任何一点的静压力大小与其作用面的方位无关，其值均相等。 3.流体平衡微分方程式（压力差公式） dp=ρ(Xdx+Ydy+Zdz) 4.等压面：平衡流体中压力相等的点所组成的平面或曲面称为等压面。等压面的两个性质：（1）平衡流体中，任一点的等压面恒与质量力正交；（2）当两种互不相混的液体处于平衡时，它们的分界面必为等压面。
传感技术学报 CHINESE JOURNAL OF SENSORS AND ACTUATORS 第29卷第11期2016年11月Vol 29No.11 Nov.2016Repairing Barrier Gaps in WSN Using Set -Based Max -Flow Algorithm * DAI Guanglin ，FANG Kai ，FANG Fei ，DAI Guoyong ，XIA Ming ， HUAN Ruohong ，MAO Keji * （College of Computer Science and Technology ，Zhejiang University of Technology ，Hangzhou 310023，China ） Abstract ：In wireless sensor networks （WSNs ），barrier coverage is a typical coverage model and is of paramount im?portance for intrusion detection.A barrier divides the area of interest into two regions such that any intruder pene?trates from one region to another is guaranteed to be detected by one or more sensor nodes in the barrier.The emer?gence of a gap in the barrier ，which is caused by running out of energy or other reasons ，may leads to the penetration without being detected.Therefore ，the functions of WSNs will be seriously influenced.A mobile nodes based gap healing method is proposed in this paper.The set -based max -flow algorithm is introduced to efficiently carry out the number of existed gaps and then the mobile nodes are scheduled to the right position to heal the gap under the con?dition of minimizing the total moving distance.Some experiments are conducted and shows the effectiveness of the proposed method.Key words ：WSN ；barrier repairing ；set -based Max -flow algorithm ；efficiency EEACC ：6150P ；6210；7230doi ：10.3969/j.issn.1004-1699.2016.11.019 基于集合最大流算法的WSN 栅栏修复方法研究 *戴光麟，方 凯，方飞，戴国勇，夏明，宦若虹，毛科技*（浙江工业大学计算机科学与技术学院，杭州310023）摘要：无线传感器网络栅栏覆盖在入侵检测方面发挥着重要作用，如何修复栅栏间隙是该领域重点研究问题之一。栅栏将监测区域划分为二部分，任何入侵目标从一个区域穿越到另外一个区域都会被栅栏中至少一个传感器节点监测到。栅栏中的节点由于某些原因过早死亡导致栅栏出现间隙，监测目标可以通过间隙而不被栅栏监测到。提出一种利用移动节点修复栅栏间隙的方法，该方法采用基于集合的最大流算法计算出能修复间隙的数量并且具有较高的效率，然后利用移动节点修复栅栏，修复过程中，移动节点的总移动距离最短。最后仿真实验验证了该方法的有效性。 关键词：无线传感器网络；栅栏修复；集合最大流算法；效率 中图分类号：TP393文献标识码：A 文章编号：1004-1699（2016）11-1742-06 栅栏覆盖是无线传感器网络领域主要的覆盖模型之一，是覆盖控制研究的热点，主要考察监测目标穿越传感器网络时被检测的情况［1］。无线传感器网络栅栏覆盖有着广泛的用途，如在国防应用中，将栅栏部署在边境线可以探测非法越境者。在环保方面，将栅栏部署在污染源周围可检测污染物的扩散情况。在林业保护方面，将栅栏部署在森林火灾现场可检测火灾蔓延情况等［2-4］。 目前国内外对无线传感器网络栅栏覆盖的相关 研究已经取得了很大的成果，栅栏覆盖的概念最早 出现在机器人传感器中［5］。舒坚等人在栅栏覆盖中 合理的引入了移动模型，能够保证以较高的概率发 现入侵者以及提早发现入侵者［6］。Anwar Saipulla 等 人提出了line -based 部署方法，该方法形成栅栏的概 率比均匀部署等方式更高，然后提出2阶段算法修 补栅栏间隙，第一阶段，FIND -GAPS 算法发现栅栏间———————————— 项目来源：国家自然科学基金项目（61379023，61401397，61302129）；浙江省公益性技术应用研究计划项目（2015C31066）；浙 江省安全生产科技计划项目（2013A1001，2013A1002） 收稿日期：2016-03-26修改日期：2016-07-10 万方数据
%最大流算法源程序 clc,clear,M=1000; u(1,2)=1;u(1,3)=1;u(1,4)=2; u(2,3)=1;u(2,5)=2; u(3,5)=1; u(4,3)=3;u(4,5)=3; f(1,2)=1;f(1,3)=0;f(1,4)=1; f(2,3)=0;f(2,5)=1; f(3,5)=1; f(4,3)=1;f(4,5)=0; n=length(u); list=[]; maxf=zeros(1:n);maxf(n)=1; while maxf(n)>0 maxf=zeros(1,n);pred=zeros(1,n); list=1;record=list;maxf(1)=M; while (~isempty(list))&(maxf(n)==0) flag=list(1);list(1)=[]; index1=(find(u(flag,:)~=0)); label1=index1(find(u(flag,index1)... -f(flag,index1)~=0)); label1=setdiff(label1,record); list=union(list,label1); pred(label1(find(pred(label1)==0)))=flag; maxf(label1)=min(maxf(flag),u(flag,label1)... -f(flag,label1)); record=union(record,label1); label2=find(f(:,flag)~=0); label2=label2'; label2=setdiff(label2,record); list=union(list,label2); pred(label2(find(pred(label2)==0)))=-flag; maxf(label2)=min(maxf(flag),f(label2,flag)); record=union(record,label2); end if maxf(n)>0 v2=n; v1=pred(v2); while v2~=1 if v1>0 f(v1,v2)=f(v1,v2)+maxf(n); else v1=abs(v1); f(v2,v1)=f(v2,v1)-maxf(n);
实用 文档 流力与流机课后作业题！ 流体力学部分 第一章 1-8解：由于油膜厚度较小，假设板间流体中的速度分布是线性的 由牛顿内摩擦定律：u F A h μ= 得 4.80.01 0.120.50.8a F h P S A u μ= ==? 420.12 1.410/856 m s μυρ-===? 1-13解：轴承与轴颈间间隙：550.150 0.0551022 D d h mm m ---====? 轴颈相对轴承的线速度：2 3.140.05950/min v rn m π==?? 由于间隙较小，由牛顿内摩擦定律： 50.45 3.140.050.1 3.140.05950 351.751060 u F A N h μ-??????===?? 因油膜摩擦而附加的阻力矩M ：0.05 351.78.7922 d M F N m ==?=? 第二章 2-7解：123234()a a p h h h h h h p γγγγγ+---+-=水水水水 123432 h h h h h h γγ-+-?= -水（） 即123432 h h h h h h ρρ-+-=-水（） 2-10解：321a p h h h p γγγ+-+=汞汞气 555=1.01310 1.3332100.7=3.0777710a p p ??+????气（0.8+0.8）-9800 2-16解：（1）2()tan H h a l g θ-= =- 20.4 1.307/3 a g m s ?=-= （2） 1.29800 1.5 1.21058422c H P h A b H N γγ=?=??=???=