近年來因節能與環保的需求，輕量化已成為各國主要政策發展之一，尤其是利用纖維強化塑膠(Fiber Reinforce Thermoplastics, FRT)來取代部份金屬零組件，其具有高強度、重量輕之材料特性，使汽車重量能大大減輕。然而應用射出成型製備纖維強化塑膠時，由於新世代產品之幾何越變越複雜，加上材料本身特性與操作條件複合後，造就了熔膠在模穴中的流動行為變得難以掌握及預測，因而纖維如何透過其微結構之變化，進而達成巨觀之強化效果，至今仍未完全地明瞭。因此，本研究之目的為利用三種不同澆口型態的標準拉伸試片模型(定義為Model I, II, III)的射出成型系統，探索因為設計、材料差異及操作條件變化下，纖維微結構如何受到影響，並進而造成巨觀形狀之變異，其中Model I, II, III分別是側邊單點進澆、直接單點進澆、雙點進澆。具體執行方法是先採用Moldex3D進行之模擬分析，之後再進行實際之試模驗證與研究。CAE模擬分析結果顯示，由於系統幾何設計之效應，三個標準拉伸試片模型呈現非常不同之翹曲變形之行為，此時我們雖然預期纖維強化了塑料強度，但卻導致Model I 往下方翹曲並且往內彎；在此同時Model II 卻往上方翹曲，而Model III 變形相對不明顯。為探索此項差異之成因，我們深入觀察纖維排向分佈(Fiber Orientation Distribution, FOD)發現，由於熔膠進入不同產品幾何時之入口效應(Entrance Effect) 引發不對稱的流場，造成不對稱之纖維排向分佈所導致特殊之翹曲變形。另外針對纖維長度對產品翹曲品質之模擬分析，可以發現纖維長度越長，對翹曲變形的改善效果更佳的明顯。此等翹曲變形改善行為也利用實驗驗證，趨勢及量值皆吻合。再則，針對纖維長度對樣品拉伸強度的變化之研究，隨著添加的纖維長度增加，拉伸強度有顯著的提升，以Model I為例，從純PP拉伸強度約為20 N/mm2，至纖維長度為25 mm約為140 N/mm2，增強了7倍。但有經過幾何入口效應之影響，不同區域所呈現之拉伸強度有明顯之差異。

Due to its great potential and capability, the fiber reinforced thermoplastics (FRT) material and technology have been applied into industry recently. However, due to the microstructures of fiber inside plastic matrix are very complex, they are not easy to be visualized. The connection from microstructures to the final shrinkage/warpage is far from our understanding.   
In this study, we have performed a benchmark with three standard specimens based on ASTM D638 where those specimens have different gate designs. Due to the geometrical effect, the warpage behaviors are quite different for those three specimens. Although we expect long fiber reinforced to enhance strength, it causes one specimen warped downward and bended inward, another warped upward, and the other slightly upward at the same time. The difference might be due to the interaction of the entrance effect of molten plastic with fiber content to cause high asymmetrical fiber orientation distribution (FOD). We further studied that the asymmetrical FOD is introduced by the entrance effect from the sprue entrance of the part to the gate of each standard specimen. This entrance effect will cause the variation of FOD and further influence the fiber distribution and the bundle phenomena. Moreover, the experimental study is also performed to validate the simulation results. From short shot testing to the warpage and bending measurement for each individual model, overall, the tendency for both numerical simulation and experimental results is in a reasonable agreement. However, some deviation still existed which needs for further study. Furthermore, in the fiber length effect study experimentally, the longer fiber reinforced the less warpage. Also, the elongation and strength of injected parts can be enhanced significantly as fiber length increased. For example, from pure PP to FRP with 25 mm fiber inside, the strength can be increased from 20 N/mm2 to 140 N/mm2 .

目錄
致謝	I
中文摘要	II
英文摘要	III
目錄	V
圖目錄	IX
表目錄	XV
符號說明	XVI
第一章 緒論	1
1.1 前言	1
1.2 文獻回顧	3
1.2.1 纖維強化塑膠	3
1.2.2 CAE模擬分析與纖維強化塑膠之關係	6
1.2.3 產品翹曲品質與強度	7
1.2.4 文獻總結	8
1.3 研究動機與目的	10
1.4 論文架構	12
第二章 射出成型纖維複合材料之微結構特性及產品巨觀之翹曲變形	14
2.1 塑膠射出成型製程簡介	15
2.2 高分子材料介紹	18
2.2.1 聚丙烯(Polypropylene)	19
2.2.2 纖維強化塑膠	20
2.3 纖維強化結構之機理介紹	23
2.3.1 纖維排向	23
2.3.2 纖維長度	24
2.3.3 纖維濃度	26
第三章 研究方法與流程	27
3.1 研究方法與流程	27
3.2 數值模擬分析與系統資訊	29
3.2.1 基本理論	29
3.2.1.1 連續方程式(Continuity Equation)	29
3.2.1.2 動量方程式(Momentum Equation)	29
3.2.1.3 能量方程式(Energy Equation)	30
3.2.1.4 本質黏度方程式	31
3.2.1.5 纖維排向模型	32
3.2.1.6 纖維長度斷裂模型	35
3.2.1.7 纖維濃度模型	36
3.2.2 產品幾何模型與模具設計	37
3.2.3 網格模型建立	39
3.2.4 材料選擇	45
3.2.5 加工參數設定	49
3.2.6 量測位置之目的與選定	51
3.2.7 數值模擬分析所使用之硬體	54
3.2.8 數值模擬分析所使用之軟體	54
3.2.9 模擬分析專案建立	54
3.3 實驗研究與相關資訊	56
3.3.1 實際射出成型之流程	56
3.3.2 射出成型系統與相關設備	58
3.3.3 射出成型成品巨觀翹曲變形量測方法	63
3.3.4 射出成型成品巨觀強度量測方法	67
3.3.5 射出成型成品纖維長度量測方法	67
3.3.6 射出成型成品微觀品質量測方法	68
第四章 結果與討論	71
4.1 理論模擬分析	71
4.1.1 巨觀翹曲變形	71
4.1.2 微觀結構	79
4.1.3 巨觀翹曲與微觀結構之關係	85
4.1.4 不同操作參數對巨觀翹曲變形之影響	86
4.1.5 入口效應	90
4.2 實驗驗證與研究	93
4.2.1 充填短射結果之比對	93
4.2.2 不同試片收縮翹曲模擬與實驗比對	96
4.2.3 纖維長度	103
4.2.4 拉伸試驗	107
4.2.5 光學方式檢視纖維微結構	113
第五章 結論	120
第六章 未來研究方向之建議	122
第七章 參考文獻	123
第八章 附錄	128
8.1 微觀結構之纖維濃度	128
8.2 單一因子之纖維排向	132
作者簡歷	142
 
圖目錄
圖 1.2 1 噴泉效應示意圖[14]	5
圖 1.2 2 模擬平板之纖維排向	5
圖 1.2 3 纖維排向之五層結構	6
圖 2.1 1 射出成型製程循環圖[39]	17
圖 2.2 1 熱塑性高分子之分子鏈	19
圖 2.2 2 熱固性高分子之分子鏈	19
圖 2.2 3 聚丙烯結構示意圖	20
圖 2.2 4 含玻璃纖維之塑料流動方向與強度之關係[43]	22
圖 2.3 1 模穴內之厚度方向纖維排向示意圖[47]	24
圖 2.3 2 纖維長度對品質的影響[48]	25
圖 3.1 1 本研究之流程圖	28
圖 3.2 1 單根纖維在空間中的方向[51]	32
圖 3.2 2 幾何模型尺寸大小(單位：mm)	38
圖 3.2 3 拉伸試片尺寸大小(單位：mm)	38
圖 3.2 4 冷卻水管配製圖	39
圖 3.2 5 模座配製圖(單位：mm)	39
圖 3.2 6 網格類型種類[52]	40
圖 3.2 7 Tetra網格解剖[52]	40
圖 3.2 8 網格層數差異之溫度分佈結果[53]	41
圖 3.2 9 不同網格尺寸之進澆口壓力曲線圖	43
圖 3.2 10 整體模穴及流道實體網格圖	44
圖 3.2 11 模穴系統網格品質圖	45
圖 3.2 12 純PP ST868M 材料黏度對剪切率之關係圖	46
圖 3.2 13 純PP ST868M 材料比容對溫度在不同壓力之關係圖	47
圖 3.2 14 純PP ST868M 材料比熱對溫度之關係圖	47
圖 3.2 15 純PP ST868M 材料熱傳導係數對溫度之關係圖	48
圖 3.2 16 比對三種不同含纖PP材料黏度對剪切速率之關係圖	48
圖 3.2 17 巨觀翹曲變形之觀測位置	52
圖 3.2 18 匯出翹曲變形後(X1)之數值方法	53
圖 3.2 19 纖維微觀變化之觀測位置	53
圖 3.2 20 本研究進行模擬分析之流程	55
圖 3.3 1 射出實驗流程圖	57
圖 3.3 2 外掛式微射出系統主要結構(專利號碼：I 258419)	59
圖 3.3 3 三柱式長纖複合材料射出機示意圖[54]	59
圖 3.3 4 射出成型機 CLF-180TXL	60
圖 3.3 5 實驗用模具	62
圖 3.3 6 (a)公模仁 (b)母模仁	62
圖 3.3 7 模具設計圖 (a)母模面 (b)公模面	63
圖 3.3 8 整個射出產品內三種不同拉伸試片之收縮翹曲量測位置	64
圖 3.3 9 量測收縮示範圖	64
圖 3.3 10 將四種材料之產品於量測位置做記號	65
圖 3.3 11 量測整個產品翹曲之示意圖	66
圖 3.3 12 電子式游標尺	66
圖 3.3 13 萬能拉力機	67
圖 3.3 14 高溫燒結爐	68
圖 3.3 15 鍍金機	69
圖 3.3 16 掃描式電子顯微鏡(TESCAN Vega II)	70
圖 4.1 1 幾何系統之收縮行為預測	72
圖 4.1 2 從Side view (1) 觀察之Model I下翹	72
圖 4.1 3 從Side view (2) 觀察之Model II上翹	73
圖 4.1 4 Model I收縮與翹曲行為	73
圖 4.1 5 Model II收縮與翹曲行為	73
圖 4.1 6 Model III收縮與翹曲行為	74
圖 4.1 7 整體試片四邊之翹曲量值模擬結果	76
圖 4.1 8 Model II短邊之P16量測節點	77
圖 4.1 9 三種材料模擬分析之翹曲量值	79
圖 4.1 10 Model I, II, III：靠近澆口處之纖維排向分析：(a) Node A, (b) Node B, (c) Node C.	82
圖 4.1 11 Model I, II, III：試片幾何中心處之纖維排向分析：(a) Node D, (b) Node E, (c) Node F	83
圖 4.1 12 Model I, II, III：試片流動末端之纖維排向分析：(a) Node G, (b) Node H, (c) Node I	84
圖 4.1 13 三種進澆方式Model在進澆口處及流動末端處纖維排向比較	86
圖 4.1 14 單一參數設變之翹曲比較圖	89
圖 4.1 15 Model I 入口效應 (紅圈處)	90
圖 4.1 16 Model I之纖維排向模擬結果 (核心層處)	91
圖 4.1 17 Model I,II在試片中心處之纖維排向	92
圖 4.2 1 15%充填之模擬波前流動及螺桿行程15 mm之實射圖	94
圖 4.2 2 22%充填之模擬波前流動及螺桿行程25 mm之實射圖	94
圖 4.2 3 50%充填之模擬波前流動及螺桿行程60 mm之實射圖	94
圖 4.2 4 65%充填之模擬波前流動及螺桿行程75 mm之實射圖	95
圖 4.2 5 75%充填之模擬波前流動及螺桿行程90 mm之實射圖	95
圖 4.2 6 100%充填之模擬波前流動及螺桿行程117 mm之實射圖	95
圖 4.2 7 實射產品幾何與尺寸(單位：mm)	96
圖 4.2 8 實際完成之產品(正視圖)	97
圖 4.2 9 實際完成之產品(上視圖)	97
圖 4.2 10 實際完成之產品(右視圖)	97
圖 4.2 11 整體試片四邊之翹曲量值模擬分析與實驗結果比較	100
圖 4.2 12 PP+30% GF (3 mm)實驗產品之整體翹曲	101
圖 4.2 13 實際射出產品之Model II上翹示意圖	101
圖 4.2 14 四種材料之四邊翹曲變形實驗結果比較	103
圖 4.2 15 纖維長度為25 mm裁切後試片	104
圖 4.2 16 燒結前試片秤重量	104
圖 4.2 17 燒結後之試片	105
圖 4.2 18 量測纖維長度之方法	105
圖 4.2 19 四種材料在原始設計之下不同Model的拉伸強度比較	109
圖 4.2 20 四種材料之拉伸應力與變形量之關係圖	109
圖 4.2 21 拉伸後試片示意圖	110
圖 4.2 22 Model I試片前端中心處撒點位置(右邊紅圈處)	111
圖 4.2 23 Model II試片前端中心處撒點位置(紅圈處)	111
圖 4.2 24 12 mm及25 mm在不同model 之拉伸破斷面	114
圖 4.2 25 雙邊進澆之熔膠匯集示意圖	115
圖 4.2 26 試片拉伸後之破斷面	116
圖 4.2 27 SEM量測之觀察位置定義	116
圖 4.2 28 Model I之破斷面於位置A之SEM圖 (50X)	117
圖 4.2 29 Model II之破斷面於位置B和C之SEM圖 (50X)	117
圖 4.2 30 Model I之破斷面於位置C2之SEM圖 (130X)	118
圖 4.2 31 Model II之破斷面於位置C3之SEM圖 (75X)	118
圖 4.2 32 Model III之破斷面於位置A之SEM圖 (50X)	119
圖 4.2 33 Model III之破斷面於位置B之SEM圖 (50X)	119
圖 8.1 1 Model I, II, III：靠近澆口處之纖維排向分析：(a) Node A, (b) Node B, (c) Node C	129
圖 8.1 2 Model I, II, III：試片幾何中心處之纖維排向分析：(a) Node D, (b) Node E, (c) Node F	130
圖 8.1 3 Model I, II, III：試片流動末端之纖維排向分析：(a) Node G, (b) Node H, (c) Node I	131
圖 8.2 1 Model I改變射速之纖維排向比較	132
圖 8.2 2 Model I改變料溫之纖維排向比較	133
圖 8.2 3 Model I改變模溫之纖維排向比較	133
圖 8.2 4 Model II改變射速之纖維排向比較	134
圖 8.2 5 Model II改變料溫纖維排向比較	134
圖 8.2 6 Model II改變模溫纖維排向比較	135
圖 8.2 7 Model III改變射速纖維排向比較	135
圖 8.2 8 Model III改變料溫纖維排向比較	136
圖 8.2 9 Model III改變模溫纖維排向比較	136
圖 8.2 10 射速30%在不同model下之纖維排向比較	137
圖 8.2 11 射速60%在不同model下之纖維排向比較	137
圖 8.2 12 射速90%在不同model下之纖維排向比較	138
圖 8.2 13 料溫240°C在不同model下之纖維排向比較	138
圖 8.2 14 料溫260°C在不同model下之纖維排向比較	139
圖 8.2 15 料溫280°C在不同model下之纖維排向比較	139
圖 8.2 16 模溫25°C在不同model下之纖維排向比較	140
圖 8.2 17 模溫50°C在不同model下之纖維排向比較	140
圖 8.2 18 模溫80°C在不同model下之纖維排向比較	141

表目錄
表 3-1不同網格尺寸進行模擬分析相關資訊	42
表 3-2本研究模擬分析專案使用之網格相關資訊	43
表 3-3本研究所使用四種材料之型號與纖維含量	46
表 3-4原始設計之加工參數設定表(模擬)	50
表 3-5單一參數設定表(模擬)	50
表 3-6 固定不變之參數表	51
表 3-7射出成型機規格	60
表 3-8實驗研究執行之射出成型操作參數設定	61
表 3-9掃描式電子顯微鏡規格	69
表 4-1單個拉伸試片之模擬量測結果(單位：mm)	74
表 4-2兩點參考點之翹曲差值比較(單位：mm)	89
表 4-3種不同拉伸試片之收縮翹曲模擬與實驗比對(單位：mm)	99
表 4-4 PP+50% GF (12 mm)之纖維含量實驗	106
表 4-5 PP+50% GF (25 mm)之纖維含量實驗	106
表 4-6 PP+50% GF (12 mm)之纖維長度實驗	107
表 4-7 PP+50% GF (25 mm)之纖維長度實驗	107
表 4-8熔膠一進入Model I表皮層及核心層處之速度	112
表 4-9熔膠一進入Model II表皮層及核心層處之速度	112
表 4-10 Model I及Model II在核心層之速度與纖維排向比較	112

[1].	M. Rohde, A. Ebel, F. F. Wolff, V. Altstädt, (2011) "Influence of processing parameters on the fiber length and impact properties of injection molded long glass fiber reinforced polypropylene", International Polymer Processing, Vol. 26(3), pp.292-303. 
[2].	S. Goris, U. Gandhi, Y. Y. Song, T Osswald, (2011) "Analysis of the process-induced microstructure in injection molding of long glass fiber-reinforced thermoplastics", Annual Technical Conference-Antec, Conference Proceedings, pp.318-326.
[3].	P. H. Foss, H. C. Tseng, J. Snawerdt, Y J Chang, W H Yang, (2014) "Prediction of fiber orientation distribution in injection molded parts using Moldex3D simulation", Polymer Composites, Vol. 35(4), 671-680. 
[4].	楊淑麗，(2006)，“長纖維增強熱塑性塑料加工和應用的新發展. 工程塑料應用”，pp.12-14。
[5].	劉維亞、張曉明，(2007). “纖維增強熱塑性復合材料及其應用”，北京化學工業出版社。
[6].	T. D. Papathanasiou, (1996) "Flow and rheology in polymer composites manufacturing", Vol.10, Journal of Chemical Technology & Biotechnology, Vol. 67(1), pp.106-106.
[7].	H. C. Tseng, R. Y. Chang, C. H. Hsu, (2013) "Phenomenological improvements to predictive models of fiber orientation in concentrated suspensions. Journal of Rheology", Vol. 57(6), pp.1597-1631. 
[8].	H. C. Tseng, R. Y. Chang, C. H. Hsu, (2013) "Method and computer readable media for determining orientation of fibers in a fluid", U.S. Patent No. 8571828.
[9].	P. H. Foss, H. C. Tseng, J. Snawerdt, Y. J. Chang, W. H. Yang, (2013) "Prediction of fiber orientation distribution in injection molded parts using Moldex3D simulation", Polymer Composites, Vol. 35(4), pp.671-680. 
[10].	K. Serope, S. Steven  Chaper 9, "Manufacturing, Engineering & Technology", Fifth Edition. 
[11].	C. T. Huang, H. C. Tseng; R. Y. Chang; S. J. Hwang, (2016) SPE Technical Papers, ANTEC2016, pp.1120-1125. 
[12].	H. Mavridis, A. N. Hrymak, J. Vlachopoulos, (1986) "Finite element simulation of fountain flow in injection molding", Polymer Engineering & Science, Vol. 26(7), pp.449-454. 
[13].	M. R. Kamal, S. K. Goyal, E. Chu, (1988) "Simulation of injection mold filling of viscoelastic polymer with fountain flow", ,Journal, Vol. 34(1), pp.94-106.
[14].	S. Kenig, (1986) "Fiber orientation development in molding of polymer composites", Polymer Composites, Vol. 7(1), pp.50-55. 
[15].	M. Gupta, K. K. Wang, (1993) "Fiber orientation and mechanical properties of short‐fiber‐reinforced injection‐molded composites: Simulated and experimental results", Polymer Composites, Vol.14(5), pp.367-382. 
[16].	K. C. Ho, M. C. Jeng, (1996) "Fiber Orientation of Short Glass Fiber Reinforced Polycarbonate Composites under Various Injection Molding Conditions", Plastics, Rubber and Composites Processing and Applications, Vol. 25, pp.469-476. 
[17].	H. Phelps J, (2009) "Processing-microstructure models for short- and long-fiber thermoplastic composites", Vol. 70-06, Section B, pp.3752.
[18].	D. Zeng, L. Lu, J. Zhou, Y. Li, Z. C. Xia, (2015) "Modeling of Long Fiber Reinforced Plastics", SAE Technical Paper 2015-01-0698.
[19].	J. F. Morris, (2009). "A review of microstructure in concentrated suspensions and its implications for rheology and bulk flow. Rheologica", Vol. 48(8), pp.909-923. 
[20].	T. Matsuoka, J. I. Takabatake, A. Koiwai, Y. Inore, (1991) "Integrated simulation to predict warpage of injection molded parts", Polymer Engineering & Science, Vol. 31(14), pp.1043-1050. 
[21].	F. Folgar, C. L. Tucker, (1984) "Orientation Behavior of Fibers in Concentrated Suspensions", Journal of Reinforced Plastics and Composites, Vol. 3(2), pp.98-119. 
[22].	H. C. Tseng., R. Y. Chang, C. H. Hsu, (2017) "Improved fiber orientation predictions for injection molded fiber composites", Composites Part A: Applied Science and Manufacturing, pp.65-75. 
[23].	S. G. Advani, C. L. Tucker, (1990) "Closure approximations for three‐dimensional structure tensors", Journal of Rheology, Vol. 34(3), pp.367-386. 
[24].	J. H. Phelps, A. E. Rahman, V. Kunc, C. L. Tucker, (2013) "A model for fiber length attrition in injection-molded long-fiber composites", Composites Part A: Applied Science and Manufacturing, Vol. 51, pp.11-21. 
[25].	S. M. Dinh, & R. C. Armstrong, (1984) A Rheological Equation of State for Semiconcentrated Fiber Suspensions. Journal of Rheology, Vol. 28(3), pp.207-227.
[26].	S. H. Maron, P. E. Pierce, (1956) "Application of ree-eyring generalized flow theory to suspensions of spherical particles", Journal of Colloid Science, Vol. 11(1), pp.80-95. 
[27].	T. Kitano, T. Kataoka, (1980) "The effect of the mixing methods on viscous properties of polyethylene melts filled with fibers", Rheologica acta, Vol. 19(6), pp.753-763. 
[28].	J. F. Morris, F. Boulay, (1999) "Curvilinear flows of noncolloidal suspensions: The role of normal stresses", Journal of Rheology, Vol. 43(5), pp.1213-1237. 
[29].	R. Zheng, N. McCaffrey, K. Winch, H. Yu, P. Kennedy, (1996) "Predicting Warpage of Injection Molded Fiber-Reinforced Plastics", Journal of  Thermoplastic Cmposite MㄋㄠVol. 9.
[30].	R. Thomas, N. McCaffrey, (1989) "The Prediction of Real Product Shrinkages from a Simulation of the Injection Molding Process", SPE Technical Papers, Vol. 35: pp.371-375. 
[31].	M. Fujiyama, H. Awaya, S. Kimura, (1977) "Mechanical anisotropy in injection‐molded polypropylene", Journal of Applied Polymer Science, Vol. 21(12), pp.3291-3309. 
[32].	H. W Cox, C. C. Mentzer, (1986) "Injection molding: The effect of fill time on properties", Polymer Engineering & Science, Vol. 26(7), pp.488-498. 
[33].	P. J. .Hogg, (1987) "The influence of flow-induced anisotropy on the impact behaviour of injection-moulded short-fibre composites", Composites Science and Technology, Vol. 29(2), pp.89-102. 
[34].	C.-I. Hsu, (2010) "Optimization of Injection Molding Process and Mechanical Properties of POM/GF Composites", 機械工程學類, p.27. 
[35].	羅壬成，(2006)，“模流分析與射出成型控制參數的優化”，工學院精密與自動化工程學程。
[36].	陳劉旺、丁金超，(2013)，“高分子加工”，高立圖書有限公司。
[37].	張榮語，(1995)，“射出成型模具設計─操作實務”，高立圖書有限公司。
[38].	塑膠射出流程，塑膠射出模具網。
[39].	魏綸群，(2008) ，“塑膠射出成型製程時間最佳化設計”，機械工程學系，元智大學。
[40].	張榮語，(1998) ，“射出成型模具設計－材料特性”，高立圖書有限公司。
[41].	聚丙烯的簡介及基本用途，(2014)，鉅亨網新聞中心。
[42].	劉佳宙，(1997) ，“不銹鋼纖維/碳纖維/ABS射出成型之製作與性質”，碩士論文，逢甲大學紡織工程研究所。
[43].	胡雲淇，(2007)，“射出成型寶典”，塑膠工業技術發展中心。
[44].	劉大佼，(1997)，“高分子加工原理與應用”，揚智文化事業股份有限公司。
[45].	T. Osswald, P. J Gramann, (2001) "Injection Molding Handbook".
[46].	謝淵清，(1990)，“工程塑膠之特性及其加工”，徐氏基金會。
[47].	S. Advani, E. Sozer, L. Mishnaevsky, (2003) "Process Modeling in Composites Manufacturing", Applied Mechanics Reviews, 56, B69. 
[48].	J. L. Thomason, M. A. Vlug, G Schipper, (1996) "Influence of fibre length and concentration on the properties of glass fibre-reinforced polypropylene: Part 3. Strength and strain at failure", omposites Part A: Applied Science and Manufacturing. Vol. 27. pp.1075-1084
[49].	L. A. Mondy, H. Brenner, S. A. Altobelli, J. R. Abbott, A. L. Graham, (1994) "Shear‐induced particle migration in suspensions of rods", Journal of Rheology, Vol. 38(2), pp.444-452.
[50].	M. Cross, (1979) "Relation between viscoelasticity and shear-thinning behaviour in liquids", Rheologica acta, Vol. 18(5), pp.609-614. 
[51].	S. G. Advani., C. L. Tucker, (1987) "The Use of Tensors to Describe and Predict Fiber Orientation in Short Fiber Composites", Journal of Rheology, Vol. 31(8), pp.751-784. 
[52].	Moldex3D. R14 Online help.
[53].	Mesh Database. CoreTech System, Hsinchu, Taiwan. 
[54].	P. C. Chang, S. J. Hwang, H. H. Lee, D. Y Huang, (2006) "Design and verification of a clamping system for micro-injection molding machine", Vol. 30, pp.413-428.
[55].	P. C. Chang, S. J. Hwang, H. H. Lee, D. Y Huang, (2007) "Development of an external-type microinjection molding module for thermoplastic polymer", Journal of Materials Processing Technology, Vol. 184(1), pp.163-172. 
[56].	張沛頎，黃登淵，李輝煌，黃聖杰，(2006)，“外掛式閥口型熱澆道微射出模組，先進精密模具、成型與相關設計/分析、設備及檢測研習會”。
[57].	楊岡錦，黃聖杰，李輝煌，(2006)，“外掛式微射出模組的控制，第五屆精密製造學術研討會”，Kaohsiung, Taiwan, pp.233-238。
[58].	S. Hamid, G. Masood, (2007) "Effects of Some Injection Molding Process Parameters on Fiber Orientation Tensor of Short Glass Fiber Polystyrene Composite (SGF/PS) ", Journal of Reinforced Plastics and Composites, Vol. 26, Issue 17, pp.1729-1741.