數值模擬顆粒堆積體行為分析方法及程式效率改良之研究.pptVIP

南亞技術學院97年度教師專題研究計畫成果簡報 摘要 大自然界中有許多具有粒狀特性的材料，工業界也有許多具有粒狀特性的原料或成品，在處理這些材料時常發生無法精確掌握粒狀材料行為特性的困擾。在欠缺足夠參考資料情況下，一般常採用試誤法進行多次實驗來取得所需之基本資料，而以試誤法進行之實驗過程往往耗時費力，成果也常因設想不夠周全而有所遺漏，甚至造成後續實際工作運轉時的重大損失。隨著計算機功能逐漸強大，近年來有相當多研究人力投入數值模擬顆粒堆積體行為程式之開發，如DEM、DDA等，希望藉由計算機高速的計算能力來建立可靠且涵蓋面廣的資料，提供工程界及工業界作為後續相關設計之參考。但目前所開發之分析程式普遍碰到因顆粒數量大、顆粒間關係複雜及耗時過長三項基本問題而導致分析效率不佳的困擾，影響實用之價值。本文著眼於第三項基本問題－耗時過長進行分析效率改善之研究，針對現有分析之方法及計算過程進行流程解析，在釐清程式之分析模式及流程後，著手去除分析程式中空轉計算之部分，並探討改良分析效率或重新建立更高效率分析流程的可行性。研究成果顯示透過分析流程改善確實可以大幅縮短分析所需時間，增加數值分析顆粒堆積體行為方法的實用價值，將可提供作為後續相關研究及設計之參考。 關鍵詞：粒狀材料、顆粒堆積體、模擬堆積 * 前言 近幾年有關數值模擬顆粒堆積及相關後續之力學及堆積體行為分析工作停滯不前，其中最主要因素為現行模式在分析顆粒堆積體時，當遭遇顆粒外形較複雜（參數較多）、數量較多之情況時，往往由於計算所耗時間也隨之大幅增加，且增加之幅度通常為幾何增加模式，造成研究及分析工作上之困擾，因此大部分之文獻仍係分析較少數量並配合較少簡單外形顆粒之模式，無法再進一步進行與實際顆粒較為接近之模擬分析，更遑論進入三維顆粒之模擬，即使是商用軟體也面臨相同發展困境。本文著眼於堆積及動靜態行為分析過程之檢討及改善工作，先行蒐集相關文獻所用之堆積方法，將其中幾種較為典型之方法改寫成相同語法程式，藉以進行分析方式研究，探討改善執行效率之可能性，嘗試進行各種程式分析流程效能探討，再依效能探討成果進行程式改寫工作。探索在相同條件下，執行各種外形顆粒配合各種分析方法改寫前後堆積之耗時比較分析。本文所分析之內容包括二維圓形、橢圓形、類多邊形、三維圓球形等形狀顆粒，所探討之分析方法如表1所示。隨著計算機運算速度不斷增加及其儲存容量也不斷提昇，數值分析方法在未來相信會成為分析之主流，此一趨勢將對顆粒堆積體之堆積模擬及動靜態行為之分析工作有很大的影響，因此分析效率的提昇是現階段最重要的課題，成果將影響使用者對於相關分析法之使用意願，以及數值模擬顆粒堆積體相關行為持續發展的動能。 表1 本文研究改善之堆積方法名稱及採用之改善策略 霣降堆積代表 改良式隨機霣降堆積法 隨機散佈霣降堆積法 垂直向壓縮 隨機霣降堆積法 水平重疊分區制 隨機搜尋堆積法 重疊分區制 隨機散佈堆積法 動態外環區域 支援盤堆積法 備註 改善策略或方法 堆積方法 成果與討論 支援盤堆積法－由於此法係由中心處向外擴散堆積，後續堆積之顆 粒不大可能堆積至中心區域，因此在判斷可放置地點時，可以將內 側區域顆粒排除在外，避免無用之判斷程序，以縮短整體之堆積時 間。本研究採用動態外環區域方式進行堆積，排除與內側區域顆粒 之重疊判斷程序，堆積耗時之改善程度相當顯著，隨顆粒數之增 加，堆積效率亦較之未改善前佳，惟實際改善效率仍會受顆粒粒徑 組合及動態外環區域設定之影響。本文以顆粒數相同之五種粒徑 (D=1.2,1.4,1.6,1.8,2.0)模擬堆積，當達600顆時，改善前後之堆 積時間比約可降至15%左右，如圖11所示。 隨機散佈堆積法－由於此法之關鍵在於盲目隨機放置新顆粒，造成 失敗率隨總堆積密度增加而增加，若能增加考量已堆積處密度之機 制，應能改善堆積之效率。本研究採用重疊分區制，放置新顆粒時 優先以密度最低分區為目標區並依隨機散佈堆積法原則進行堆積， 放置新顆粒後立即推算因堆積而改變之分區密度，依此方式持續進 行堆積。本文以顆粒數相同之五種粒徑(D=1.2,1.4,1.6,1.8,2.0) 模擬堆積，堆積空間100×100，當達600顆時，改善前後之堆積時間 比約可降至33%左右，成果顯示確實可以改善堆積效率，改善情形 如圖12所示。採用此法仍須注意當密度到增加到相當程度後又會面 臨可堆積空間相當少之情況，仍然會有無法完成堆積之問題。 隨機搜尋堆積法－此法雖已經多次改善，但如堆積之空間水平方向 較大時仍會有耗時稍多之問題。本研究採用重疊分區制，由於計算 新顆粒堆放位置時水平方向可大幅縮短搜尋時間，整體堆積效率仍 會進一步提昇，惟實際改善效