本文藉由擴張史磋公式維度，設計相對應之邊界元素程式，解決二維磁電彈材料的孔洞、裂紋與界面等問題。並利用商用軟體Ansys與解析解佐以模擬與驗證。
另外，過去由於矩陣運算的便利性，因而採用以矩陣運算見長的MATLAB軟體來撰寫此邊界元素法，但又因邊界元素法需使用大量迴圈運算，使運算速度大幅拖慢。本研究以MATLAB軟體與FORTRAN軟體撰寫相同架構之邊界元素程式，客觀且公正的表現兩者於撰寫難度與迴圈運算速度的優劣比較。

By extending the dimension of Stroh formalism and using the associated boundary element, we can successfully solve two dimensional linear magneto-electro-elasticity defect problems. We also use commercial software and explicit solution to prove the correctness of our boundary element design.
In the past, because the powerful matrix operators of Matlab, it has been used to design the boundary element code. But, the massive number of loop slows down the computation process. In this paper, we use the computer language "Fortran" to rewrite the boundary element code under the same design structure and compare the computational efficiency of these two different computer languages.

[1] Suchtelen, J.V., “Product Properties: A New Application of Composite Materials,” Philips Research Reports, Vol. 27, No. 1, pp. 28-37, 1972.

[2] Run, A. M. J. G. V., Terrell, D. R., Scholing, J. H., “An in situ grown eutextic magnetoelectric composite material,” Journal of Materials Science, Volume 9, Issue 10, pp. 1710-1714, 1974.

[3] 林睿哲，“尖晶石CoFe2O4與鈣鈦礦Pb(Zr, Ti)O3積層複合薄膜的多重鐵性與磁電耦合特性研究”，博士論文，國立清華大學材料科學工程學系，2008.

[4] Ryu, J., Alfredo, V. C., Kenji, U., Kim, H. E.,“Magnetoelectric Properties in Piezoelectric and Magnetostrictive Laminate Composite,” Jpn. J. Appl. Phys. 40, pp. 4948-4951, 2001.

[5] Marco, A., “Magnetoelectric Effect in Piezoelectric /Magnetostrictive Multilayer (2-2) Composites,” Journal of Intelligent Material Systems and Structures, 1994 pp. 501-513, 1994.

[6] Nan, C. W., “Magnetoelectric effect in composites of piezoelectric and piezomagnetic phases,” Phys. Rev. B 50, pp. 6082–6088, 1994.

[7] Benveniste, Y., “Magnetoelectric effect in fibrous composites with piezoelectric and piezomagnetic phases,” Phys. Rev. B 51, pp16424–16427, 1995.

[8] Daga, A., Ganesan, N., Shankar, K., “Harmonic Response of Three-Phase Magneto-electro-elastic Beam Under Mechnical,Electrical and Magnetic Environment,” Journal of Intelligent Material Systems and Structures, Volume 20, issue 10 , pp. 1203-1220, 2009.

[9] Eerenstein, W., Mathur, N. D., Scott, J. F. , “Multiferroic and magnetoelectric　materials,” Nature 442, pp. 759-765, 2006.

[10] Liang, Y. C. , Hwu, C., “Electromechanical analysis of defects in piezoelectric, ” Smart Mater. Struct 5, pp. 314-320, 1996.

[11] Chen, Y. C., Hwu, C., “Green’s Functions for Anistropic/Piezoelectric Bimaterials and Their Applications to Boundary Element Analysis,” CMES, vol.57, no.1, pp.31-50, 2010.

[12] Haojiang, D., Aimin, J., “A boundary integral formulation and solution for 2D problems in magneto-electro-elastic media,” Computers & Structures Volume 82, Issues 20-21, pp. 1599-1607, 2004.

[13] 李政錡，“多孔洞異向性板之邊界元素設計”，碩士論文，成功大學航空太空工程研究所，2012.

[14] 詹育儒，“含任一方向裂紋之磁壓電材料破壞分析”，碩士論文，成功大學機械工程研究所，2006.

[15] 劉晉奇，“含界面裂紋之雙磁電彈楔形結構反平面破壞分析”，博士論文，成功大學航空太空工程研究所，2005.

[16] Hsieh, M. C., “Extended Stroh-Like Formalism for The Magneto-Electro-Elastic Composite Laminates, ”

[17] Hwu, C., “Anisotropic Elastic Plates,” Springer, 2010.

[18] 李仲喬，“結合有限元素法及邊界元素法分析多孔洞之問題”，碩士論文，成功大學航空太空工程研究所，2013.