鋼鐵製程中所使用的設備，如螺絲、冷卻輥、齒輪、旋轉軸、補強肋等零部件，常有段差或凹槽設計以滿足功能性需求或避免干涉，這些設計特徵容易因楔型缺口而產生應力奇異性。當此類設備因斷裂而損壞時，斷裂面或延伸中的裂縫起始點通常出現在這些容易引發應力奇異性的部位。因此，必須計算破壞參數以評估裂縫可能出現的位置或裂縫延伸的趨勢。目前，師門已能使用異向性彈性力學自創軟體 AEPH 中的 H積分(s543(2-D)&s544(3-D))模組計算破壞參數：應力的奇異階次及其相關的應力強度因子。相比之下，商用軟體 Ansys 的破壞參數計算工具需要在預處理階段先定義裂縫後才能求解破壞參數，同時無法計算如楔型缺口或界面角或界面裂縫等奇異點之破壞參數，這限制了其應用範圍。 AEPH 只需獲得特定位置周圍的應力、應變、位移值，即可用以計算破壞參數，因此其適用範圍更為廣泛。然而，使用 AEPH 計算時需提供從商用軟體匯出並整理成 AEPH 可接受的輸入格式的相關輸入檔，這在評估多個應力奇異點的破壞參數時較為耗時，不符合實際應用需求。
本文主要探討如何利用 Ansys 巨集語言與撰寫 AEPH 的所採用之 Matlab 相互連結，最終僅需在 Ansys 完成結構計算後，便能直接調用 AEPH 計算使用者指定位置的破壞參數，並在原圖面上直接標示形成裂縫圖譜，供設計者進行後續評估，並比對標準案例的結果，以進一步驗證此計算方法的正確性與便利性。

In steel manufacturing, components such as screws, cooling rollers, and gears often feature wedge-shaped notches that create stress singularities and become common sites for fracture initiation. Accurate fracture parameter calculations are crucial for predicting crack propagation. Our research utilizes AEPH, an anisotropic elasticity software that includes H-integral modules (s543 for 2-D and s544 for 3-D) to compute fracture parameters such as stress intensity factors and related singular orders. Unlike commercial software such as ANSYS, AEPH does not require explicit crack geometry and can handle singularities at wedge-shaped notches, interface corners, and interfacial cracks. However, AEPH requires input files in a specific format, which complicates its use in practical engineering workflows. To address this, our study integrates the ANSYS Parametric Design Language (APDL) with MATLAB scripts to automate AEPH invocation following ANSYS analysis. This enables fracture parameters to be computed at specified locations and automatically annotated on the original model. The method is validated through benchmark cases to ensure accuracy and efficiency, providing a reliable tool for fracture analysis in engineering design.

[1] A. Öchsner, Continuum damage mechanics: Springer, 2016.
[2] C. Hwu, Anisotropic elasticity with Matlab: Springer, 2021.
[3] C. Hwu, Anisotropic elastic plates: Springer Science & Business Media, 2010.
[4] W.-T. Zhang, “Evaluationof Stress Intensity Factors at Points Displaying Interface Singularity.,” 2024.
[5] 吳偉特, “破壞準則(Failure Criteria),” 地工技術, vol. 30, pp. 106-112, 1986.
[6] H. Yu, and M. Kuna, “Interaction integral method for computation of crack parameters K–T–A review,” Engineering Fracture Mechanics, vol. 249, pp. 107722, 2021.
[7] D. Weichert, and M. Schulz, “J-integral concept for multi-phase materials,” Computational materials science, vol. 1, no. 3, pp. 241-248, 1993.
[8] 龔曙光, 謝桂蘭, and 黃雲清, ANSYS 參數化編程與命令手冊: 機械工業出版社, 2009.
[9] H.-H. Lee, Finite element simulations with ANSYS workbench 2023: Theory, applications, case studies: SDC publications, 2023.