強化塑膠因其品質穩定及強度高，且有低壓力，低成本，速度快的優點，所以廣泛的被航太與汽車工業界所使用；而轉注成型便是應用於製造強化塑膠的一種方法。在樹脂轉注成型充填過程中，氣泡的生成對於機械性質的品質影響甚巨；當樹脂流經纖維束時，巨觀流（纖維束外區域）和微觀流（纖維束內區域）之流動消長關係，是因纖維蓆滲透率的差異，而造成流動面產生不一樣的流動速度與壓力梯度，進而影響氣泡的生成。
本文試以蒙地卡羅法模擬及實驗觀測轉注成型製程中充填階段毛細位移對流動面與初始氣泡消長之影響。分析在定流率下樹脂流經纖維束時，用達西定律來描述巨觀流及微觀流之流動情形，而以簡單正交法產生格點，用蒙地卡羅法來預測流動面前進的位置，而做出數值模擬分析及探討是否合乎物理性與可行性。其模擬結果顯示滲透度的比值、毛細係數與充填條件與流動機率的取決皆有密切的關連，並利用連續過度鬆弛法來增快程式執行的收斂速度；在實驗方面則發現不同黏度的流體在相同的滲透率之下，其形成的初始氣泡並不會相同，且崩裂時間也會有所差異；滲透率在一定的比值下時，對黏度較高的流體其初始氣泡變大；若壓力梯度大，初始氣泡溢出的微小氣泡則溢出速度快，崩裂時間變快。

摘要 ……………………………………………………………………………..Ⅰ
目錄 …………………………………………………………………………….Ⅱ
圖目錄 ………………………………………………………………………….Ⅴ
表目錄 ………………………………………………………………………….Ⅶ
符號說明 ……………………………………………………………………….Ⅷ
第一章 緒論………………………………………………………………………1
1-1 前言……………………………………………………………………1
1-2 文獻回顧………………………………………………………………3
1-3 研究方向………………………………………………………………8
第二章 理論模式……………………………………………………………….10
2-1 巨、微觀流動之物理模式………………………………………….10
2-2 毛細力對多孔性介質的影響……………………………………….13
2-3 流動機率的計算原則……………………………………………….15
第三章 數值方法與實驗原理………………………………………………….18
3-1 數值方法…………………………………………………………….18
3-1-1 蒙地卡羅法…………………………………………………..18
3-1-2 數值之定義方式………………………………………………18
3-1-3 數值之計算步驟……………………………………………..19
3-1-4 連續過度鬆弛法………………………………………………21
3-2 實驗原理…………………………………………………………….22
3-2-1 達西定律………………………………………………………22
3-2-2 空孔度…………………………………………………………23
3-2-3 滲透度…………………………………………………………23
3-2-4 相似理論分析…………………………………………………24
第四章 實驗裝備與步驟……………………………………………………….27
4-1 實驗材料…………………………………………………………….27
4-1-1 纖維蓆…………………………………………………………27
4-1-2 潤滑油…………………………………………………………28
4-2 實驗設備…………………………………………………………….29
4-2-1 精密天平………………………………………………………29
4-2-2 材料拉伸試驗機………………………………………………29
4-2-3 馬達、幫浦……………………………………………………30
4-2-4 高壓油管、節流閥……………………………………………31
4-2-5 油壓缸…………………………………………………………31
4-2-6 模具……………………………………………………………32
4-3 影像拍攝與處理…………………………………………………….33
4-4 實驗步驟…………………………………………………………….34
4-5 注意事項…………………………………………………………….35
第五章 結果與討論…………………………………………………………….37
5-1 數值方面…………………………………………………………….37
5-1-1 格點獨立………………………………………………………37
5-1-2 CPU Time……………………………………………………..38
5-1-3 程式之發展模式與其合理性…………………………………39
5-2 實驗方面…………………………………………………………….42
5-2-1 初始氣泡剛形成階段………………………………………..43
5-2-2 初始氣泡崩裂階段……………………………………………44
第六章 結論…………………………………………………………………….46
6-1 結論………………………………………………………………….46
6-2 未來發展方向……………………………………………………….47
參考文獻………………………………………………………………………..49

1. Dave, R. "A Unified Approach to Modeling Resin Flow During Composite Processing," Journal of Composite Materials., 24, pp. 22-41 （1990）
2. F. Trochu, and R. Gauvin, "Limitations of a Boundary-Fitted Finite Difference Method for the Simulation of the Resin Transfer Molding Process," Journal of Reinforced Plastics and Composites, 11, pp.772-786 （1992）
3. L. Fong, J. Xu and L. J. Lee, "Performing Analysis of Thermoformable Glass Fiber Mats-Deformation Modes and Reinforcement Characterization," Polymer Composites, 15, pp. 134-146 （1994）
4. T. J. Wang, C. H. Wu, and L. J. Lee, "In-Plane Permeability Measurement and Analysis in Liquid Composite Molding," Polymer Composites, 15, pp. 278-288 （1991）
5. K. Han, L. Trevino, W. B. Young, L. J. Lee and M. J. Liou, "Fiber Mat Deformation During Mold Filling in Structural RIM," The 46th Annual Conference of The Composites Institute of the The Society of The Plastics Industry Inc., 4, pp. 18-21（1991）
6. C. J. Wu, L. W. Hourng and J. C. Liao, " Numerical and Experimental Study on the Edge Effect of Resin Transfer Molding," J. of Resinforced Plastics and Composites, 14, pp. 692-722 （1995）
7. C. J. Wu, and L. W. Hourng, "Permeable Boundary Condition For Numerical Simulation in Resin Transfer Molding," Polymer Engineering and Science, 35, pp. 1272-1281 （1995）
8. J. H. Eckier and D. A. Rust, "Proceedings of The ASM Internal/EDS," Advanced Composites Conference., 109, pp. 1147-1160 （1987）
9. W. B. Young, K. Rupel, K. Han, L. James Lee, * and Ming J. Liou*, "Analysis of Resin Molding in With Preplaced Fiber Mats. II：Numerical Simulation and Experiment of Mold Filling," Polymer Composites, 12, pp. 30-38 （1991）
10. M. V. Bruschke and S. G. Advani, "A Finite Element/Control Volume Approach to Mold Filling in Anisotropic Porous Media," Polymer Composites, 11, pp. 398-405 （1990）
11. A. W. Chan, and S. T. Hwang, "Molding Nonisothermal Impregnation of Fibrous Media With Reactive Polymer Resin," Polymer Composites, 32, pp. 310-317 （1992）
12. A. W. Chan, and S. T. Hwang, "Modeling Resin Transfer Molding of Polyimide（PMR-15）/Fiber Composites," Polymer Composites, 14,pp. 524-528 （1993）
13. W. B. Young, "Tree-Dimensional Nonisothermal Mold Filling Simulation in Resin Transfer Molding," Polymer Composites, 15, pp. 118-127 （1994）
14. S. A. Wymer, and R. S. Engel, "A Numerical Study of Nonisothermal Resin Flow in RTM with Heated Uniaxial Fibers," J. Composite Materials, 28, pp. 53-65 （1994）
15. D. A. Weitz, J. P. Stokes, R. C. Ball, and A. P. Kushnick, "Dynamic Capillary Pressure in Porous Media：Orgin of the Viscous-Fingering Length Scale," Physical Review Letters, 59, pp. 2967-2970 （1992）
16. G. L. Batch, Y. T. Chen, and C. W. Macosko, "Capillary Impregnation of Aligned Fibrous Beds：Experiments and Mold," J. of Reinforced Plastics and Composites, 15, pp. 1027-1051（1996）
17. N. Patel, V. Rohatgi and L. J. Lee, "Micro Scale Flow Behavior and Void Formation Mechanism During Impregnation Through a Undirection Stitched Fiberglass Mat," Polymer Engineering and Science, 35, pp. 837-851 （1995）
18. D. G. Kiriakidis, G. H. Neale, and E. Mitsoulis, "The Effects of Capillary Forces on the Island Size Distribution in Two-Fluid Immiscible Displacement Flow in Porous Media," J. Phys. A：Math. Gen., 24, pp. 3797-3805 （1991）
19. R. Chandler, J. Koplik, K. Lerman, and J. F. Willemsen, "Capillary Displacement and Percolation in Porous Media," J. Fluid Mech., 119, pp. 249-267 （1982）
20. V. Yamakov* and A. Milchev, "Polymer Chain in a Flow Through a Porous Medium：A Monte Carlo Simulation," Physical Review E , 56, pp. 7043-7052 （1997）
21. Binder, K. "Application of Monte Carlo Methods to Statistical Physics" Reports on Progress in Physics, 60, pp. 487-559 （1997）
22. Binder, K. and Paul, W. "Monte Carlo Simulation of Polymer Dynamics-Recent Advances," Journal of Polymer Science Part B-Polymer Physics, 35, pp. 1-31 （1997）
23. A. D. Mahale, R. K. Prud'homme and L. Rebenfeld, "Quantitative Measurement of Voids Formed During Liquid Impregnation of Nonwoven Multifilament Glass Networks Using an Optical Visualization Technique," Polymer Engineering and Science, 32, 319-326 （1992）
24. P. G. De Gennes, "Wetting：Static and Dynamics," Reviews of Modern Physics, 57, Part 1, pp. 827-863 （1985）
25. P. G. De Gennes, "Dynamics Capillary Pressure in Porous Media," Europhysics Letters, 5, pp. 689-691 （1988）
26. R. Lenormand, "Invasion Percolation in an Etched Network：Measurement of a Fractal Dimension," Physical Review Letters, 54, pp. 2226-2229 （1985）
27. D. G. Kiriakidis, and G. H. Neale, "The Effects of Medium Thickness on the Island Size Distribution in Immiscible Displacement Flow in Porous Media," Chem. Eng. Comm., 123, pp. 127-134 （1993）
28. K. Han, L. J. Lee* and S. Nakamura, "Dry Spot Formation and Changes in Liquid Composite Molding：II-Modeling and Simulation," Journal of Composite Materials, 30, pp.1475-1492 （1996）
29. E. W. Washburn, "The Dynamics of Capillary Flow," The Physical Review 2nd series, 12, pp. 273-283 （1991）
30. Shillcock, J. C. and Seifert, U. "Monte Carlo Simulation of Multiple Pores in a Model Fluid Membrane," Biophysical Journal, 72, pp. 319-329 （1997）
31. Riechelmann, D. and Nanbu, K. ?-Dimensional Simulation of Wavy Taylor Vortex Flow by Direct Simulation Monte-Carlo Method," Physical Fluids, 9, pp. 811-813 （1997）
32. 黃文淵, "蒙地卡羅法在熱傳導逆運算問題的應用" 國立交通大學機械研究所碩士論文（1988）
33. 洪明宏, "樹脂轉注成形製程中氣泡消長機制之分析" 國立中央大學機械研究所碩士論文（1998）
34. 張修銘, "轉注成形製成中氣泡形成機制之可視化探討" 國立中央大學機械研究所碩士論文（1997）
35. C. J. Wu and L. W. Hourng, "Permeability Boundary Condition For Numerical Simulation in Resin Transfer Molding," Polymer Engineering and Science, 35, pp. 1272-1281 （1995）
36. Albert W. Chan and Sun-Tak Hwang, "Anisotropic In-Plane Permeability of Fabric Media," Polymer Engineering and Science, 31, pp.1233-1238 （1991）
37. J. Bear, Dynamics of Fluids in Porous Media, Dover Publication, New York （1972）
38. F. A. L. Dullien, Porous Media Fluid Transport and Pore Structure, Academic Press, New York （1979）
39. G. P. Androutsopoulos and R. Mann, "Evaluation of Mercury Porosimeter Experiments Using a Network Pore Structure Model," Chemical Engineering Science, 34, pp. 1203-1212 （1979）
40. D. A. Anderson, J. C. Tannehill, R. H. Pletcher, "Computational Fluid Mechanics and Heat Transfer," Hemisphere Publishing Corporation, New York, （1987）
41. J. C. Tannehill, D. A. Anderson, R. H. Pletcher, "Computational Fluid Mechanics and Heat Transfer," Taylor and Francis, Washington, （1997）
42. R. Gauvin, M. Chibani and P. Lanfntaine, "The Molding Pressure Distribution in Resin Transfer Molding," Journal of Reinforced Plastics and Composites, 6, pp. 367-377 （1987）
43. 廖致欽, "轉注成形之滲透率的量測與流場之觀察" 國立中央大學機械工程研究所碩士論文 （1994）
44. B. G. Gebart, P. Gudmundson and C. Y. Lundemo, "An Evaluation of Alternative Injection Strategies in RTM," The 47th Annual Conference of Composites Institute of the Plastics Industry Inc., Session 16D/1-Session 16D/8 （1992）
45. C. Y. Chang and L. W. Hourng, "Numerical Simulation for the Transverse Impregnation in Resin Transfer Molding," J. of Reinforced Plastics and Composites, Accepted （1997）
46. David Bensimon, * Leo P. Kadanoff, S. Liang, B. I. Shraiman, and Chao Tang, "Viscous Flows in Dimensions," Reviews of Modem Physics, 58, pp. 977-999. （1986）