本文基於破壞力學理論，模擬一具有邊緣初始裂縫之試片，承受垂直於裂縫方向之疲勞負載，探討在各種不確定因子的影響下所造成裂縫延伸的情形與變化。藉由既有實驗結果驗證電腦模擬的實用性，並由最小平方誤差之觀念，建立計算模擬所需材料常數之數值方法。將多次模擬所得之破裂週次繪製分佈函數，結果發現以對數常態分佈函數最為適合。比較各因子對破裂週次的影響，發現除了初始裂縫長度外，試片寬亦有顯著的影響。若能將影響較顯著的變異納入模擬考量中，將可使模擬更趨近於實驗結果，可供應用於疲勞裂縫成長壽命的評估。

參考文獻
[1]姚枚.陳鴻賓合著,"金屬機械性能",全華科技圖書股份有限公司出版,1996.
[2]P.C. Paris, M. P. Gomez and W. P. Anderson, "A rational analytic theory of fatigue",The Trend in Engineering, Vol.13, pp.9-14, 1961.
[3]P. C. Paris and F. Erdogan, "A critical analysis of crack propagation law," Transaction of the ASME, Journal of Basic Engineering, Vol. 85, No. 4, pp.528-534, 1963.
[4],J. Schijve and D. Broke, "The influence of the mean stress on the propagation of fatigue cracks in aluminum alloy sheets, " Nat. Aerospace Inst. Amsterdam TR-M-2111, 1963.
[5] E.K. Walker, "Effects of environments and complex load history on fatigue life, "ASME STP 462, pp.1-14, 1970.
[6] R. G. Forman, and V. E. Kearney, and R. M. Engle, "Numerical analysis of crack propagation in cyclic-loaded structures," Transaction ASME, Journal of Basic Engineering, Vol. 89, No. 3, pp. 459-464, 1976.
[7]R. J. Donahue, H. M. Clark, P. Atanmo, R. Kumble, andA. J. McEvily, "Crack opening displacement and the rate of fatigue crack growth, " International Journal of Fracture Mechanics, Vol. 8, pp.209-219, 1972.
[8]A. J. McEvily, "On closure in fatigue crack growth, " ASTM STP 982, American Society for Testing and Materials, Philadelphia, pp.35-43, 1988.
[9]C.R. Chiang, "A unified theory of fatigue and crack growth: a statistical approach," International Journal of Fracture, Vol. 53, pp.337-342, 1992.
[10]C.R. Chiang, "Numerical simulation on fracture behavior of materials, "國科會專題研究計劃成果報告NSC86-2212-E007-006, 1997.
[11]黃郁青,"裂縫的疲勞成長電腦模擬",碩士論文,國立清華大學,1999.
[12]王琮瑋,"疊層複合材料中疲勞裂縫成長之電腦模擬",碩士論文,國立清華大學, 2000.
[13]C.R. Chiang, "Kinked crack in an anisotropic material, "
Engineering Fracture Mechanics, Vol. 39, No. 5, pp. 927-930, 1991.
[14]陳兆勛 譯著,"破裂力學之實際應用,"國立編譯館出版,1999.
[15]A. Ohta and E Sasaki, "Influence of stress ratio on the threshold level for fatigue crack propagation in high strength steels, "Engineering Fracture Mechanics, Vol 9,no.2, pp. 307-315, 1977.
[16]D. N. Lal and T. K. G. Namboodhiri, " Model for the effect of the mean stress on the threshold condition for fatigue crack propagation , " Materials Science and Engineering , A. Vol. A130, no. 1, pp. 37-49. 1990.
[17]C. R. Chiang, "Threshold stress intensity factor of fatigue cracks, "Engineering Fracture Mechanics, Vo.49, No. 1, pp. 29-33, 1994.
[18]D. Broek,陳文華、張士欽 譯,"基本工程破裂力學",國立編譯館出版,1995.
[19]E. Campbell, W. Gerberich and H. Underwood, " Application of fracture mechanics for selection of metallic structural materials", American Society for Metals,1982
[20]M. H. Kalos and P. A. Whitlock, " Monte Carlo Methods", Wiley & Sons, New York, 1986.
[21]P.A.W. Lewis, A.S. Goodmen, and J. M. Miller, "A pseudo random number for the system/360,"IBM Syst. J., Vol. 8, pp. 136-146, 1969.
[22] A.M. Law and W.D. Kelton. "Simulation Modeling and Analysis" McGraw-Hill, New York, 2nd edition, 1991.
[23]M. B. Buczek and C. T. Herakovich, "A normal stress criterion for crack extension direction in orthotropic composite materials, "Journal of Composite Materials, Vol. 19, pp.544-553, 1985.
[24]陳志忠,"鋁合金裂縫成長行為之可靠度研究",碩士論文,
國立台灣大學,2001.
[25]張智星 編著,"MATLAB程式設計與應用",清蔚科技出版,2001
[26]周啟,蔡秀芬,"普通物理實驗數據之讀取、統計、處理與分析"，全國第三屆大學校院普通物理學及實驗改進研討會,1990。
[27]許淵賓, "張開型裂縫成長方向的統計預測",碩士論文,國立清華大學, 2008.