1、第六章第六章 流動阻力和水頭損失流動阻力和水頭損失 張炯張炯2022-4-91主要內容61 61 流動阻力和水頭損失的分類流動阻力和水頭損失的分類62 62 黏性流體的兩種流態黏性流體的兩種流態63 63 沿程水頭損失與剪應力的關系沿程水頭損失與剪應力的關系64 64 圓管中的層流運動圓管中的層流運動65 65 紊流運動紊流運動66 66 紊流的沿程水頭損失紊流的沿程水頭損失67 67 局部水頭損失局部水頭損失68 68 邊界層概念與繞流阻力邊界層概念與繞流阻力1、研究內容、研究內容:恒定不可壓縮流體中的機械能損失。2、流動阻力及其分類：、流動阻力及其分類： 由于流體存在粘性（內 因）及由固體

2、邊壁發生變化（外因）所產生的阻礙流體運動的力。 61 流動阻力和水頭損失的分類一、水頭損失的分類一、水頭損失的分類按固體邊壁情況的不同，分為：按固體邊壁情況的不同，分為：（1）沿程阻力沿程阻力（2）局部阻力局部阻力由流體粘性所產生的阻礙流體運動的力。由流體粘性所產生的阻礙流體運動的力。 在邊壁沿程無變化（邊壁形狀、尺寸、過流方向均在邊壁沿程無變化（邊壁形狀、尺寸、過流方向均無變化）的均勻流流段上，產生的流動阻力。無變化）的均勻流流段上，產生的流動阻力。由固體邊壁發生改變所產生的阻礙流體運動的力。由固體邊壁發生改變所產生的阻礙流體運動的力。 在邊壁沿程急劇變化，流速分布發生變化的局部在邊壁沿程急

3、劇變化，流速分布發生變化的局部區段上（如管道入口、異徑管、彎管、三通、閥區段上（如管道入口、異徑管、彎管、三通、閥門等），集中產生的流動阻力。門等），集中產生的流動阻力。3、水頭損失的分類：、水頭損失的分類：（2）局部損失）局部損失 hj （1）沿程損失）沿程損失 hf總水頭總水頭損失損失hw流體克服流體克服沿程沿程阻力所損失掉的能量。阻力所損失掉的能量。流體克服流體克服局部局部阻力所損失掉的能量。阻力所損失掉的能量。兩者不相互干擾時兩者不相互干擾時hw = hf + hj注：沿程水頭損失均勻分布在整個流段上，與流段的長度成正比。常見的發生局部水頭損失區域常見的發生局部水頭損失區域 jcjbj

4、afcdfbcfabjfwhhhhhhhhh總水頭損失：總水頭損失：二、水頭損失的計算公式二、水頭損失的計算公式1.1.沿程阻力沿程阻力沿程損失沿程損失( (長度損失、摩擦損失長度損失、摩擦損失) )gvdlhf22達西達西- -魏斯巴赫公式魏斯巴赫公式 沿程阻力系數沿程阻力系數 d d 管徑管徑 v v 斷面平均流速斷面平均流速 g g 重力加速度重力加速度2.2.局部阻力局部阻力局部損失局部損失gvhj22 局部阻力系數局部阻力系數 v v 對應的斷面平均速度對應的斷面平均速度（2）紊流）紊流流體質點在流動過程中發生相互混摻，流體質點在流動過程中發生相互混摻， 流體質點的軌跡與其流向不平行

5、。流體質點的軌跡與其流向不平行。6-2 粘性流體的兩種流態粘性流體的兩種流態流體質點作規則運動，相互不干擾，流體流體質點作規則運動，相互不干擾，流體 質點的運動軌跡與流向平行。質點的運動軌跡與流向平行。（1）層流）層流一、雷諾實驗（一、雷諾實驗（1880 1883年）年）1、實驗裝置：Qhf實驗證明：沿程損失實驗證明：沿程損失h hf f與與V V的關系與流態有關。的關系與流態有關。1122顏色水K2K1LhfBDCAEVK1VK2LghfLgV2、實驗方法、實驗方法：使水流的速度分別使水流的速度分別由小到大由小到大由大到小由大到小改變改變。觀測現象，并測出相應的數值（觀測現象，并測出相應的數

6、值（ v、 hf ）。）。3、實驗結果與分析：、實驗結果與分析：（1）實驗現象：）實驗現象：1 流速流速 v 由小由小大：大：當當 v vcr 時，時， 玻璃管中的紅線消失；玻璃管中的紅線消失；2 流速流速 v 由大由大?。盒。寒敭?v vcr 時，時， 玻璃管中的紅線又重新出現。玻璃管中的紅線又重新出現。vcr上臨界流速； vcr下臨界流速。（2）流態的劃分：）流態的劃分：vcr v vcrv vcr層流；紊流；可為層流也可為湍流，保持原有流態。（3）流速）流速 v 與沿程損失與沿程損失 hf 的關系的關系 ： lg vcr lg vlg hflg vcr ABDCE在雷諾實驗中，在雷諾實驗

7、中，測得多組測得多組 hf 與與 v 的值，得到的值，得到 v hf 的對應關的對應關系，在對數紙上系，在對數紙上點繪出點繪出 v hf關系曲線關系曲線. 如圖所示。如圖所示。k2=1.752.0k1=1.01 當流速由小到大時當流速由小到大時曲線沿曲線沿AEBCD 移動；移動；2 當流速由大到小時當流速由大到小時曲線沿曲線沿DCEA 移動移動 lg vcr lg vlg hflg vcr ABDCE分析：分析：1 AE 段： lg vcr lg vlg hflg vcr ABDCE層流v CD段：紊流v vcr，為直線段，hf = kv 1.752.0 m2=1.752.0， 直線的斜率 3

8、 EC 段：vcr v vcr ，為折線段。屬過渡區 ，狀態取決于原流動狀態。 lg vcr lg vlg hflg vcr ABDCE由于沿程損失與流態有關，故計算hf 時，應先判斷流體的流動型態。二、流態的判別標準臨界雷諾數 Rec1、圓管：實驗發現： 上臨界流速vcr不穩定，受起始擾動的影響很大；下臨界流速vcr穩定，不受起始擾動的影響。一般取Rec=2300vdvdRe雷諾數雷諾數23002000Redvcc用臨界雷諾數作為流態判別標準，只需計算用臨界雷諾數作為流態判別標準，只需計算出流管的雷諾數出流管的雷諾數將Re值與Rec=2300比較，便可判別流態： ReRec，則vRec，則v

9、vc，流動是紊流； Re=Rec，則v=vc，流動是臨界流。vdReAR hbR水力半徑 ; 濕周,為過流斷面與固體邊壁相接觸的周界。x = b +2hR=Ax2、非圓管：、非圓管：3.用量綱分析說明雷諾數的物理意義 Lv dnduAam粘性力慣性力 LvLLvL223 Re慣性力與粘性力作用之比判斷流態圓滿管流（如圖右）圓滿管流（如圖右）4412dddR以水力半徑以水力半徑R R為特征長度，相應的為特征長度，相應的臨界雷諾數臨界雷諾數575vRRRec 例：某段自來水管，d=100mm，v=1.0m/s。水溫10，（1）試判斷管中水流流態？（2）若要保持層流，最大流速是多少？解:（1）水溫為

10、10時，水的運動粘度，由下式計算得：則： 即：圓管中水流處在紊流狀態。 （2）要保持層流，最大流速是0.03m/s。63 沿程水頭損失與剪應力的關系沿程水頭損失與剪應力的關系本節只對簡單均勻流作分析，找出本節只對簡單均勻流作分析，找出 hf 與與 的關系。的關系。以圓管為例以圓管為例一、均勻流基本方程一、均勻流基本方程1、沿程損失（、沿程損失（伯努利方程）伯努利方程）gvgv22222211因為流體的流動是恒定、均勻流，因為流體的流動是恒定、均勻流，所以有：所以有：故有故有)()(2211pzpzhf2、均勻流基本方程（力的平衡）：、均勻流基本方程（力的平衡）：如果流體的流動為均勻流，則流體的

11、受力應平衡。如果流體的流動為均勻流，則流體的受力應平衡。（1）分析受力，如圖：）分析受力，如圖：重力： ()慣性力： 01 質量力：質量力：lgAG2 表面力：側面所受切力：兩斷面所受壓力：p2 A2 (-) p1 A1 (+)soFl 濕周0壁面剪應力0fhgRgRJl120cos0p Ap AgAll 21zz 01212flppzzhgggA （2）基本方程：）基本方程：Rlhf0J R0適用于層流與紊流，只要是均勻流即可。gJRv0*v*動力速度、阻力速度、剪切速度。3、圓管過流斷面上剪應力分布規律、圓管過流斷面上剪應力分布規律表明有壓圓管均勻流過流斷面上切應力呈直線分布。表明有壓圓管

12、均勻流過流斷面上切應力呈直線分布。如圖右所示一水平恒定圓管均勻流，R=r0/2，則由上式可得 同理可得： 所以圓管均勻流切應力分布為 或 0 xy二、沿程損失的普遍表達式二、沿程損失的普遍表達式達西公式達西公式 gvRhf2412gvdlhf22適用于圓形管路適用于非圓形管路 適用于層流與 紊 流。64 圓管中的層流運動圓管中的層流運動一、流動特征一、流動特征由于層流各流層質點互不摻混，對于圓管來說，各層質點沿平行管軸線方向運動。與管壁接觸的一層速度為零，管軸線上速度最大，整個管流如同無數薄壁圓筒一個套著一個滑動。二、流速分布、流量、平均流速二、流速分布、流量、平均流速 1、流速分布： 積分

13、當r=r0時， u=0由drduJrg2)(2204rrgJu204rgJumax（2）管軸中心處流速最大，為：注：（1）圓管層流過水斷面上流速分布呈旋轉拋物面分布；（3）管壁處流速最小，為： umin= 0 umax0ArrgJrdrrrgJudAQ4002208240)(2、流量：3、斷面平均流速：208rgJAQvmaxu21 即圓管層流的平均流速是最大流速的一半。 gvdlgvdlRdvlhef226432222三、沿程水頭損失的計算三、沿程水頭損失的計算 208rJv且：且：1、在雷諾實驗中，已知如果流體的流態為、在雷諾實驗中，已知如果流體的流態為層流層流， 則有：則有：hf =k

14、v 。 而由以上理論也證明而由以上理論也證明 hf 與與 v 的一次方成正比。的一次方成正比。2、在圓管、在圓管層流層流中，中，只與只與Re有關。有關。即：即：=f ( Re )64Re例例1 1 =0.85g/cm=0.85g/cm3 3的油在管徑的油在管徑100mm,100mm,v v=0.18cm=0.18cm2 2/s/s的管中以的管中以v v=6.35cm/s=6.35cm/s的速度作層流運動，求的速度作層流運動，求 （1 1）管中心處的最大流速；（）管中心處的最大流速；（2 2）在離管中心）在離管中心r=20mmr=20mm處的流速；處的流速；（3 3）沿程阻力系數）沿程阻力系數

15、；（；（4 4）管壁切應力）管壁切應力0 0及每及每kmkm管長的水頭損失。管長的水頭損失。 解 :（1）求管中心最大流速 （2）離管中心r=20mm處的流速2220r44gJurrgJumax)(寫成 當r=50mm時，管軸處u=0，則有0=12.7-K52，得K=0.51，則r=20mm在處的流速 （3 3）沿程阻力系數）沿程阻力系數 先求出Re （層流） 則 （4）切應力及每千米管長的水頭損失 紊流運動較為復雜，到目前尚處于半經驗紊流運動較為復雜，到目前尚處于半經驗階段，此處只介紹與流動阻力損失有關的階段，此處只介紹與流動阻力損失有關的理論。理論。65 紊流運動紊流運動一、紊流的特點一、

16、紊流的特點紊流的流動參數在時間、空間上的分布具有隨機性，同時又具有統計上的確定性。脈動脈動流體的運動要素在某段時間內以一定值為中流體的運動要素在某段時間內以一定值為中 心，隨時間不斷改變的現象，稱脈動。心，隨時間不斷改變的現象，稱脈動。二、紊流脈動與時均化二、紊流脈動與時均化1、脈動現象、脈動現象由于在紊流運動中，各流體質點間相互混摻，由于在紊流運動中，各流體質點間相互混摻， 使流體各運動要素發生脈動現象。使流體各運動要素發生脈動現象。2、時均化概念、時均化概念脈動流速時均值：010dtuTuTxx運動要素在一定時段內時大時小，但總運動要素在一定時段內時大時小，但總圍繞一定值（平均值）上下波動

17、。圍繞一定值（平均值）上下波動。（1）設某點的瞬時流速為 ux（以 x 方向為例）：ppp瞬時壓強：瞬時流速：xxxuuuTxxdtuTu01時均流速：常用紊流度常用紊流度N來表示紊動的程度來表示紊動的程度xzyxuuuuN)(22231想一想：紊流的瞬時流速、時均流速、脈動流速、斷面平均流速有何聯系和區別？ 3、脈動現象及時均化的意義：、脈動現象及時均化的意義：（1）由于紊流中各點的運動要素隨時間改變，所以應屬）由于紊流中各點的運動要素隨時間改變，所以應屬 于非恒定流動，但若其時均值不隨時間改變，即可于非恒定流動，但若其時均值不隨時間改變，即可 將其視為恒定流（時均恒定流），故所有關于恒定將

18、其視為恒定流（時均恒定流），故所有關于恒定 流的公式均可用于此。即采用了時均化概念，紊流流的公式均可用于此。即采用了時均化概念，紊流 脈動有可能按恒定流來處理。脈動有可能按恒定流來處理。（2）由于存在脈動現象，故紊流與層流相比，其速度分）由于存在脈動現象，故紊流與層流相比，其速度分 布、溫度分布、懸浮物分布都更趨平均化。表現在布、溫度分布、懸浮物分布都更趨平均化。表現在 動能修正系數和動量修正系數上，紊流時近似為動能修正系數和動量修正系數上，紊流時近似為 1， 而層流時差別較大。而層流時差別較大。三、紊流的切應力三、紊流的切應力1、紊流切應力表達式：、紊流切應力表達式：在紊流中的水流阻力除了粘

19、性阻力外，液體質點混摻和動量交換還將產生附加的切應力，簡稱紊流附加應力。因此，紊流的水流阻力可表達為:yuxdd 121時均流速產生的粘滯切應力時均流速產生的粘滯切應力脈動流速產生的附加切應力脈動流速產生的附加切應力（1）粘滯切應力可由牛頓內摩擦定律解決；）粘滯切應力可由牛頓內摩擦定律解決；（2）附加切應力主要依靠紊流半經驗理論解決。）附加切應力主要依靠紊流半經驗理論解決。普朗特混合長度理論； 卡門理論； 泰勒理論。ReRe數較小時，占主導地位數較小時，占主導地位ReRe數很大時，數很大時，121紊流附加切應力的半經驗理論：紊流附加切應力的半經驗理論：22_dyudl vu四、粘性底層四、粘性

20、底層粘性底層：圓管作紊流運動時，靠近管壁處存在著一薄層，粘性底層：圓管作紊流運動時，靠近管壁處存在著一薄層，該層內流速梯度較大，粘性影響不可忽略，紊流該層內流速梯度較大，粘性影響不可忽略，紊流附加切應力可以忽略，速度近似呈線性分布，附加切應力可以忽略，速度近似呈線性分布， 這這一薄層就稱為粘性底層。一薄層就稱為粘性底層。 特點：1 流動近似層流； 2 時均流速為線性分布； 3 紊流附加切應力可忽略。 紊流核心：粘性底層之外的液流統稱為紊流核心。紊流核心：粘性底層之外的液流統稱為紊流核心。 （1）粘性底層的厚度：）粘性底層的厚度：dRe8 .320粘性底層的厚度為：66 紊流的沿程水頭損失紊流的

21、沿程水頭損失任務： 確定紊流流動中確定紊流流動中的值。的值。確定方法： 以紊流的半經驗理論為基礎，結合實驗結果，以紊流的半經驗理論為基礎，結合實驗結果，整理成整理成的半經驗公式；的半經驗公式； 直接根據實驗結果，綜合成直接根據實驗結果，綜合成的經驗公式。的經驗公式。gvdlhf22gvRlhf242一、尼古拉茲實驗一、尼古拉茲實驗 尼古拉茲進行了人工粗糙管的阻力實驗。尼古拉茲進行了人工粗糙管的阻力實驗。他用經過仔細篩分后，粒徑為他用經過仔細篩分后，粒徑為的沙子粘的沙子粘貼在管壁上形成人工粗糙管。貼在管壁上形成人工粗糙管。 并對并對d/=30d/=30、6161、120120、252252、50

22、4504和和10141014六種管道進行實驗，得到六種管道進行實驗，得到-Re-Re對數曲線，對數曲線，稱為尼古拉茲實驗曲線。稱為尼古拉茲實驗曲線。一、尼古拉茲實驗一、尼古拉茲實驗（1）絕對粗糙度）絕對粗糙度粗糙凸出固體壁面的平均高度。粗糙凸出固體壁面的平均高度。（2）相對粗糙度）相對粗糙度絕對粗糙度與過流斷面上某一特絕對粗糙度與過流斷面上某一特性性 幾何尺寸的比值。幾何尺寸的比值。層流沿程阻力系數層流沿程阻力系數只是只是Re的函數，紊流中沿程阻力系數的函數，紊流中沿程阻力系數除和流動狀況除和流動狀況(Re)有關外，由于壁面粗糙是對流動的一種有關外，由于壁面粗糙是對流動的一種擾動，因此壁面粗糙

23、是影響沿程阻力系數的另一個重要因擾動，因此壁面粗糙是影響沿程阻力系數的另一個重要因素。素。1、實驗方法：、實驗方法：（1）選擇一組不同相對粗糙度的人工粗糙管。）選擇一組不同相對粗糙度的人工粗糙管。 壁面粗糙一般包括粗糙突起的高度、形狀以及疏密和壁面粗糙一般包括粗糙突起的高度、形狀以及疏密和排列等許多因素。尼古拉茲將經過篩選的均勻砂粒緊密地排列等許多因素。尼古拉茲將經過篩選的均勻砂粒緊密地粘在管壁表面，做成人工粗糙。粘在管壁表面，做成人工粗糙。 用砂粒的突出高度（砂粒直徑）表示壁面的絕對粗用砂粒的突出高度（砂粒直徑）表示壁面的絕對粗糙。糙。/d表示相對粗糙。表示相對粗糙。分析得出：)/,(dRf

24、e（2）實驗裝置：）實驗裝置：gvdlhf22)100lg(lgeR（3）具體實驗內容：）具體實驗內容： 對不同相對糙度的管路，分別測得一系列 Q， hf ，t， l ，d 。1 計算：2 做做 Re /d曲線（如圖）曲線（如圖）尼古拉茲圖可分為五個區域：尼古拉茲圖可分為五個區域：I.層流區層流區 II. 過渡區過渡區 III.紊流光滑區紊流光滑區 IV. 紊流過渡區紊流過渡區 V. 紊流粗糙區紊流粗糙區2、實驗分析、實驗分析（據的變化特征分為五個區）：（1）層流區（）層流區（區區）：）：ab線，lgRe 3.36, Re 2300，=f(Re )（2）層流）層流紊流過渡區（紊流過渡區（區區）

25、bc線，3.36 lgRe 3.6 2300 Re 3.6，Re 4000，=f(Re )（4）紊流過渡區（）紊流過渡區（區）：區）：=f(Re ，/d ) cd，ef之間的曲線族，不同的相對粗糙管的實驗點分布落在不同的曲線上。表明即與Re有關，又與/d有關。范圍越大。越小，則其水力光滑區，dfRe3、實驗意義：（1）較具體的揭示了影響）較具體的揭示了影響的各因素之間的具體關系。的各因素之間的具體關系。（2）補充了普朗特半經驗理論，為推導紊流半經）補充了普朗特半經驗理論，為推導紊流半經 驗公式提供了依據。驗公式提供了依據。（5）紊流粗糙區、阻力平方區（區）：=f(/d )ef右側水平的直線簇，

26、不同的相對粗糙管的實驗點分別落在不同的的水平直線上。表明只與/d有關，與Re無關。1 (Re , /d )， 2 各流區流速分布資料。（2）普朗特動量傳遞理論 （紊流速度對數表達式）； （1）尼古拉茲實驗資料：二、人工粗糙管路二、人工粗糙管路和和v計算計算公式推導依據：公式推導依據：1、紊流光滑區、紊流光滑區2、紊流粗糙區、紊流粗糙區d7321.lg51221.lgeR3、紊流光滑區向紊流過渡區：、紊流光滑區向紊流過渡區： 因人工管路與工業管路在此區存在較大因人工管路與工業管路在此區存在較大差異，研究此區無實際意義，故在此不做介紹。差異，研究此區無實際意義，故在此不做介紹。另：由于層流紊流過渡

27、區 （區）極窄，無實際 意義，故不作研究。4、光滑區、過渡區、粗糙區判別標準、光滑區、過渡區、粗糙區判別標準 （2）過渡區：）過渡區： 0.40 60 5 Re* 70 （1）光滑區：）光滑區： 0. 40 Re* 60 Re* 70*vRe 粗糙雷諾數粗糙雷諾數當量粗糙度當量粗糙度 指和工業管路粗糙區指和工業管路粗糙區相等的、相等的、 同粒徑的、人工粗糙的砂粒高度。同粒徑的、人工粗糙的砂粒高度。人工管路為等粒徑，即人工管路為等粒徑，即相等。相等。工業管路粒徑不等，故常用工業管路粒徑不等，故常用當量粗糙度當量粗糙度表示。表示。三、實用（工業）管路三、實用（工業）管路的確定的確定實際工程中常用的

28、管道，稱為工業管道。實際工程中常用的管道，稱為工業管道。64工業管道流區的劃分及工業管道流區的劃分及 的的計算公式：計算公式： 1. 1.紊流流區的劃分紊流流區的劃分： . .過渡粗糙區過渡粗糙區： . .水力光滑區水力光滑區： . .紊流粗糙區紊流粗糙區：。，則水力光滑區不存在。若范圍：40008080Re4000dd85. 024160Red范圍：85. 024160Re80dd范圍：652. 2. 的的計算公式：計算公式：層流：層流：Re64水力光滑區：水力光滑區：，布拉修斯經驗公式。25. 0Re3164. 0紊流粗糙區紊流粗糙區(阻力平方區阻力平方區)：，普朗特公式。8 . 0Rel

29、og21274. 12log21d66過渡粗糙區：過渡粗糙區：柯列勃洛克公式。,Re51. 27 . 3log21d 柯列勃洛克公式適用于紊流的三個流區，故它又柯列勃洛克公式適用于紊流的三個流區，故它又被稱為紊流沿程損失系數被稱為紊流沿程損失系數 的綜合的綜合計算公式計算公式。 謝才公式謝才公式 通常用于均勻流。將達西魏斯巴赫公式變形 lhdgvf22lhJf以d=4R， ，代入上式，整理得 RJCRJgv8式中： v 斷面平均流速 R 水力半徑 J 水力坡度 C 謝才系數（1）曼寧公式：611Rnc （2）巴甫洛夫公式：適用于紊流粗糙區 謝才系數謝才系數式中： n 綜合反映壁面對水流阻滯作用

30、的系數，稱為粗糙系數。見表6-3。 R 水力半徑，單位為m。穆迪圖(Moody)gvRhf2412R水力半徑x濕周A過流斷面面積R =Ax四、非圓管沿程損失計算四、非圓管沿程損失計算 例 有一新的給水管道，管徑d=400mm，管長 l =100m，糙率n =0.011，沿程水頭損失hf=0.4m，水流屬于紊流粗糙區，問通過的流量為多少？ 解 管道過水斷面面積 水力半徑 利用曼寧公式計算C值，則 所以流量 例. 給水管為長40m，直徑d75mm的新鑄鐵管，流量Q7.25L/s, 水溫10C。試求該管段的沿程水頭損失。6-7 局部水頭損失局部水頭損失 局部水頭損失與沿程水頭損失一樣，不同的局部水頭

31、損失與沿程水頭損失一樣，不同的流態所遵循的規律也不同。流態所遵循的規律也不同。 目前能用理論公式推導出的只有突然擴大的局部損失計算，目前能用理論公式推導出的只有突然擴大的局部損失計算，其它計算均由實驗而得。因大多數的流動為其它計算均由實驗而得。因大多數的流動為紊流紊流，故在此，故在此只研究紊流時的局部損失計算。只研究紊流時的局部損失計算。層流區： hj v紊流區： hj v 2gvhj22 一、局部水頭損失的一般分析一、局部水頭損失的一般分析 局部水頭損失產生的原因：邊壁的變化導致漩渦區的產生，引起能量損失。突然擴大突然縮小閘閥三通匯流管道彎頭管道進口分離區分離區分離區分離區分離區分離區分離區分離區分離區分離區分離區分離區分離區分離區gvhj22 二、局部損失計算公式 =f (局部阻礙的形狀 )，