基于圖像處理的細集料休止角研究

2019-10-29 08:48:48單佳欣馮勇劉婷婷姜旭東李雪婷陳嘉樂陳維翰

中國建材科技 2019年4期

單佳欣馮勇劉婷婷姜旭東李雪婷陳嘉樂陳維翰

（長安大學，陜西西安 710021）

0 前言

在實際工程應用中，集料堆積場設計、集料運輸轉場等均涉及顆粒堆積問題。

國內學者對休止角的實驗測試主要是利用裝置法。通常是在測試過程中將待測顆粒經下料裝置以一定的高度自由落在方格紙上堆積成錐，多次測量其底面直徑，求得其半徑的平均值，再測得其錐頂高度，計算休止角，此方法操作麻煩，底面半徑和圓錐高度存在較大誤差，所得測量結果精度低，參考價值不大。

國內學者對細集料的研究主要集中在顆粒級配[1]、顆粒孔隙特征分析[2]、路用性能[3]和內摩擦角[4]等方面，而較少對細集料的休止角進行研究。細集料是工程實際中一種重要的材料，在道路橋梁、房屋建筑、水利水電等工程中都有重要的應用，在集料堆積場建設、集料運輸轉場等工程實踐中，都需要科學的設計，提高場地利用率，預防安全事故。另外，采砂場砂石采集、臨時堆放等，目前都沒有系統的理論作為支撐，多靠經驗設計，存在較大風險。因此，細集料休止角的研究對集料堆積場設計、集料運輸轉場、采砂場等工程實踐有重要的指導意義。

國內外對休止角的研究主要集中于天然砂。劉建寶等[5]對沙丘沙的粒徑、磨圓度、表面粗糙度、分選性、顆粒之間的分子力、礦物成分以及含水率等影響因素對休止角的影響；張紅武等[6]通過試驗測得砂與小礫石的水下休止角隨著粒徑增大由31°增大至40°，水下摩擦系數與粒徑的平方根近似成正比，而且當泥沙相對密度減少時，水下休止角將增大；孟震等[7]進行了無黏性泥沙休止角與表層沙摩擦角試驗；王成華等[8]通過對溜沙坡天然休止角的測定，分析了影響休止角的因素；金臘華等[9]認為顆粒級配、粒徑、粒面磨蝕程度是影響散粒體模型沙水下休止角的主要因素。國外早在20世紀40年代，就已經展開了對各種不同材料的休止角廣泛研究[10，11]。一批早期國外學者總結了包括粒徑、磨圓度、表面粗糙度等影響因子在內的休止角推測方程或經驗公式用于估算休止角的大小；1953年，Lane[12]通過試驗發現，泥沙水下休止角與砂粒粒徑成正比，并隨顆粒形狀而變化。

綜上所述，國內學者對細集料的研究主要集中在常規性能的研究，本研究則對細集料的休止角進行研究。國內外對休止角的研究主要集中于天然砂，本研究則選擇從細集料出發進行休止角的探究。在探究影響因素上，前人對天然砂的天然休止角的影響因素進行了大量研究，主要側重于粒徑、重度、顆粒級配、表面粗糙度等對天然休止角的影響，本文則著重探究粒徑和傾倒條件對休止角的影響。

1 休止角測試方法分析

1.1 傳統測試方法理論分析[13]

休止角的傳統測試方法具有沖擊法的缺點，結果都存在偶然誤差之外的系統誤差。

在對成堆現象的分析中，沙粒收到自身重力而產生的沖力作用會沖擊已形成的堆頂，頂部會出現削峰現象；同時，由于慣性作用顆粒下滑到接近底盤的位置時會繼續向外部蔓延一部分，導致錐形堆并不是理想的圓錐形狀，堆積模型如圖1所示。

圖1 堆積模型

其中：tanθ為真實休止角的正切值；tanθ′為實測休止角的正切值；S為錐形堆的表觀高度；R為錐形堆的表觀半徑；μ為高度偏差；ρ半徑偏差。

由于存在高度偏差μ和半徑偏差ρ，因此會導致tanθ≥tanθ′，即休止角的實測值小于真實值，這就是傳統方法最主要的誤差來源。

在堆積過程中，一般μ和ρ不隨S的變化而變化，當S→∞時，tanθ=tanθ′。

1.2 基于圖像處理的休止角測量方法

本測量方法由CCD圖像采集模塊、信息提取模塊、信息處理模塊三大模塊集合而成（具體流程示意圖如圖2）。

圖2 基于圖像處理的細集料休止角測量方法流程示意圖

具體步驟如下：

1）通過CCD圖像采集模塊采集視頻資料（如圖3示），調整拍攝高度和角度，使得鏡頭中心與自制休止角測試儀底座上表面在同一水平面上；

圖3 數據采集模塊示意圖

2）通過信息提取模塊抽取符合休止角定義的關鍵幀；

3）將（2）中所得幀導入到信息處理模塊中，根據圖片擬合錐形堆的母線和底，使用計算機測角工具，完成休止角測量；

圖4 信息處理模塊示意圖

4）同一幀重復步驟3）三次，取其平均值。

1.3 新方法的優點

1）本方法消除傳統方法高度偏差μ和半徑偏差ρ產生的系統誤差；

2）本方法減少了傳統方法由于人工操作產生的偶然誤差，運用圖像處理技術的休止角處理方法對數據進行詳細的分析，極大提高了測量精度；

3）本方法對于特殊情況，比如有粘性的粉體或者吸濕的粉體，還可能形成第二休止角（如圖5所示），即崩塌的正在流動的粒子斜面與底面的夾角。對此，傳統方法實測休止角θ′既不是第一休止角也不是第二休止角，而新的測量方法只需要在信息處理模塊再擬合一條直線，用測角工具測得即可，突破了傳統方法的禁區。

圖5 特殊情況下第二休止角示意圖

2 試驗

2.1 試驗材料

本次試驗所用細集料為河沙（已除雜、烘干），經標準方孔砂石篩（如圖6所示）長時間逐級篩分后，得到①-⑤五組實驗樣品（如圖7所示），其粒徑見表1。

圖6 標準方孔砂石篩

圖7 五組試驗樣品

表1 各組別樣品粒徑區間

2.2 試驗方案設計

本文通過自制的砂土休止角測定儀，研究不同粒徑、不同傾倒條件下，細集料休止角的變化規律；具體試驗方案如表2所示。為了使試驗數據更加具有準確性，每個試驗砂樣進行了3組試驗，當3組試樣結果偏差在規定值內時，取3次測量結果的平均值。

表2 試驗方案

2.3 試驗方法

顆粒材料休止角常用的試驗方法有：注入法、排出法和旋轉容器法[14]。為了模擬工程實際情況，本試驗選取注入法，注入法是通過漏斗裝置將所選的沙土顆粒材料傾倒到一個光滑的板面上形成圓錐體，圓錐體的母線與水平面的夾角即為所測角（如圖8所示）。

其中：d為漏斗內徑；H為漏斗下緣距離底座上表面的豎直距離。

具體休止角測量方法為本研究提出的基于圖像處理的休止角測量方法。

3 分析與討論

3.1 休止角與顆粒粒徑的關系

表3為不同粒徑區間細集料在不同內徑漏斗下休止角的試驗結果。由試驗結果可知，當漏斗內經分別為11.60mm、12.40mm和17.90mm時，粒徑區間為0.16-0.315mm的休止角均最小，分別為32.5°、32.8°和30.5°，粒徑區間為1.25-2.5mm的休止角均最高，達38.5°、41.5°和37.3°。

圖8 試驗裝置示意圖

表3 不同粒徑區間細集料在不同內徑漏斗下的休止角

圖9為細集料的休止角θ與粒徑D的關系，從中可以看出，細集料休止角與粒徑的關系曲線分為兩段：當粒徑D＜0.315mm時，細集料休止角與粒徑成反比關系；而當粒徑D＞0.315mm時，細集料休止角與粒徑成正比關系。

原因分析：當粒徑D＜0.2mm時，范德華力與靜電力發揮作用，粒子容易發生聚集，內聚力超過細集料重力，妨礙了細集料的重力行為，所以此區間的細集料主要為黏附性粉狀粒子，隨其粒徑的增大，粒子間的接觸點數減少，粒子間的范德華力、靜電力降低，休止角減小。而當粒徑D＞0.2mm時，主要為散粒體，休止角的大小主要受顆粒的球度和顆粒表面的磨蝕程度的影響，隨著細集料粒徑的增大，球形度減小，粒面粗糙度減小，但休止角卻隨著粒徑的增大而變大，可見對于不同粒徑散粒體細集料的休止角大小主要取決于顆粒的球形度，粒徑主要是通過影響細集料的形體規則度（球形度）來間接影響細集料休止角的大小[15-17]。

圖9 細集料的休止角θ與粒徑D的關系

3.2 休止角與傾倒條件的關系

休止角本應為顆粒物質的固有屬性，但經過對多次堆積試驗觀察分析，隨著傾倒條件的改變，實測休止角也會隨之改變，因此，為了探究傾倒條件是如何影響休止角的，首先對顆粒下落成堆進行建模分析。

基本假定：

1）顆粒流是均勻連續的；

2）顆粒流之間互不接觸。

設顆粒流下落到錐形堆頂部時對頂部的沖量為dI；顆粒流末段的動量為dp；

其中，ρ為顆粒的表觀密度；A漏斗內壁橫截面積；g為重力加速度；dx為漏斗口處距離微元；其他符號同前。

綜上所述，傾倒條件對實測休止角的影響主要通過傾倒流量和傾倒流速兩個影響因素來實現的。

3.2.1 休止角與傾倒流量的關系

傾倒口徑是指本休止角測試方法中漏斗的口徑d。為探究傾倒口徑對細集料休止角的影響，設置在相同的粒徑（①-⑤五組）和傾倒高度（H=10.2cm）下，改變漏斗口徑，分別為11.60mm、12.40mm、17.90mm。得到如圖10所示結果，可以得到，隨著漏斗口徑的增大，細集料的休止角總體上呈現下降趨勢。

原因分析：由建模分析得到，隨著漏斗口徑d的增大，顆粒流對錐形堆頂部的沖量會越大，導致削峰現象更加嚴重，錐角處的慣性作用也更大，綜合使得錐形堆坍落更加嚴重，從而導致實測休止角呈現整體下落趨勢。

3.2.2 休止角與傾倒流速的關系

傾倒高度是指本休止角測試方法中漏斗下緣距離底座的豎直距離H。為探究傾倒高度對細集料休止角的影響，設置在相同的粒徑（0.315-0.63mm）和傾倒口徑（d=12.40mm）下，改變漏斗放置高度，分別為10.2cm、15.7cm、22.7cm。實驗數據如表4得到如圖11所示結果，可以得到，細集料的休止角與傾倒高度成正比。