M-積分是研究具有裂縫之物體破壞行為之重要參數。本論文結合有限元素法研究針對具有多裂縫之三維線彈性材料物體，受混合載重作用下，計算其M-積分的數值分析。首先依序針對具有任意形狀之二維裂縫問題以及具有任意形狀之三維裂縫問題，進行M-積分式的理論介紹，其次証明M-積分具有與積分曲面無關的性質。
在三維問題，對單裂縫的M-積分計算結果顯示與積分曲面無關和具有與原點無關的特性；關於多裂縫問題的M-積分計算，考慮一特殊狀況，在一個無限域的線彈性結構體內包含多裂縫，受到遠端均勻載重作用，在此特殊的邊界情況下，選取任意積分原點且任意包含所有裂縫的積分區域，計算出的M-積分值皆會相同。
接著利用不同裂縫間距去探討裂縫間距對於M-積分的影響，而結果呈現出，隨裂縫間距縮小的情況下所得到的M積分值也越小。
二維情況下M-積分的物理意義為裂縫面形成時能量變化的兩倍，三維中M-積分的物理意義則為裂縫面形成時能量變化的三倍。

M-integral is an important parameter for fracture analysis. In this research, a numerical procedure, incorporated with the finite element method, that is developed for calculation of the 3D linear elastic solid subjected to mixed-mode load with cracks. First, we verify M-integral for arbitrary-shaped cracks in 2D problems and the arbitrary-shaped cracks in 3D problems. Second,we verify the property of surface independence.
In 3D single crack problems, the M-integral is verified to be surface independent and origin independent. In 3D multiple cracks problems, the result of M-integral also has the property of surface independence and origin independence on specific boundary conditions.
We also discuss the relation between cracking space and the results of the M-integral.
In 2-D situation, M-integral is equal to twice the surface energy required for the formation of the whole cracks. In 3-D situation, due to the different geometry of the cracks, the M-integral appears to be equal to triple the surface energy required for the formation of the whole cracks.

[1] J.H. Chang , A.J. Chien , “Evaluation of M-integral for anisotropic elastic media with multiple defects”, International Journal of Fracture , Volume 114, Number 3, 2002 , Pages 267-289.
[2] J.H. Chang, D.J. Wu, “Stress intensity factor computation along a non-planar curved crack in three dimensions”, International Journal of Solids and Structures, Volume 44, Issue 2, 15 January 2007, Pages 371-386.
[3] Chao-Shi Chen, Chia-Hau Chen, Ernian Pan , “Three-dimensional stress intensity factors of a central square crack in a transversely isotropic cuboid with arbitrary material orientations”, Engineering Analysis with Boundary Elements, Volume 33, Issue 2, February 2009, Pages 128-136.
[4] G.Z. Wang, X.L. Liu, F.Z. Xuan, S.T. Tu, “Effect of constraint induced by crack depth on creep crack-tip stress field in CT specimens”, International Journal of Solids and Structures, Volume 47, Issue 2, February 2010, Pages 51-57.
[5] V.I. Kushch, A.S. Sangani , “Stress intensity factor and effective stiffness of a solid containing aligned penny-shaped cracks” ,International Journal of Solids and Structures, Volume 37, Issue 44, November 2000, Pages 6555-6570.
[6] Yi-Heng Chen,“M-integral analysis for two-dimensional solids with strongly interacting microcracks. Part I: in an infinite brittle solid” ,International Journal of Solids and Structures, Volume 38, Issue 18, May 2001, Pages 3193-3212.
[7] Han J.J., Dhanasekar M., “Modelling cracks in arbitrarily shaped
finite bodies by distribution of dislocation”, International Journal of Solid and Structures, Volume 41, Issue 2, January 2004, Pages 399-411.
[8] Jan Sladek , Vladimir Sladek, “Evaluation of the Elastic T-stress in Three-dimensional Crack Problems Using an Integral Formula”, Volume 101, Number 4, 2000 , Pages 47-52.
[9] P.Isaksson , R. Hagglund ,“Strain energy distribution in a crack-tip region in random fiber networks”, International Journal of Fracture , Volume 156, Number 1, 2009 , Pages 1-9.
[10] M. Lorentzon, K. Eriksson , “A path independent integral for the crack extension force of the circular arc crack”, Engineering Fracture Mechanics , Volume 66, Issue 5, July 2000 , Pages 423-439.
[11] Y.Z. Chen , X.Y. Lin , Z.X. Wang , “Evaluation of the stress intensity factors and the T-stress in periodic crack problem”, International Journal of Fracture , Volume 156, Number 2, 2009 , Pages 203-216.
[12] Xin Wang , “Two-parameter characterization of elastic–plastic crack front fields: Surface cracked plates under tensile loading”, Engineering Fracture Mechanics , Volume 76, Issue 7, May 2009 ,Pages 827-982.
[13] X.L. Fu, G.F. Wang, X.Q. Feng , “Surface effects on mode-I crack tip fields: A numerical study”, Engineering Fracture Mechanics , Volume 77 Issue 5 , May 2010 , Pages 1031-1202.
[14] P.V. Jogdand, K.S.R.K. Murthy, “A finite element based interior collocation method for the computation of stress intensity factors and T-stresses”, Engineering Fracture Mechanics Volume 77 Issue 5 , May 2010 , Pages 1031-1202.
[15] Jui-Hung Chang and Deh-Jun Peng, (2004).
[16] 葉俊彬，「應用Jk積分於均質與非均質材料在裂縫延伸時之能量釋放計算」，碩士論文，國立中央大學土木工程研究所，中壢(1996)。
[17] 陳哲彬，「複合材料垂直於介面上裂縫之Jk積分計算」，碩士論文，國立中央大學土木工程研究所，中壢(1997)。
[18] 鄔德傳，「三維裂縫之Jk積分與應力強度因子之數值計算」，博士論文，國立中央大學土木工程研究所，中壢(2005)。
[19] 康宇權，「混合載重下的三維多裂縫問題之M-積分」，碩士論文，國立中央大學土木工程研究所，中壢(2010)。
[20] 鐘靈光，「三維多裂縫問題之M-積分」，碩士論文，國立中央大學土木工程研究所，中壢(2011)。