關鍵詞：lammps模擬，裂紋擴展，拉伸，單晶鋁，ovito

隨著納米技術的發展，人們的注意力逐漸從宏觀物體轉向微觀物體。由于納米晶體金屬及合金材料具有優越的物理、化學、力學特性，越來越受到人們的重 視，但是材料的缺陷嚴重影響著人們的安全，所以研究裂紋的擴展機制成為一項重要的研究課題。 由于裂紋擴展在原子尺度上進行，目前傳統的宏觀連續介質力學已經無法滿足材料微觀尺度變形機理的研究。近幾十年來，分子動力學方法作為一種計算機模擬技術，解決了由大量原子組成的系統動力學問題，它能夠揭示在微觀尺度下材料的變形和斷裂的實質過程。尤其近幾年來，計算機的飛速發展也為研究裂紋擴展提供了可能。

為了深入研究單晶鋁在裂紋存在時的行為，建立了兩種三維單晶鋁帶有預制初始裂紋的模型。這兩種模型是基于單晶鋁的嚴格面心立方晶格結構設計的，其中晶格常數a被設定為0.405 nm，這是鋁在室溫下的典型晶格尺寸。采用可視化分析處理軟件ovito對編程得到的原子坐標數據，具體模型如圖a、b所示：

圖(a)和圖(b)分別為帶有不同裂紋的單晶鋁初始模型，使用顏色將模型簡單分區，在黃色區域加載Z方向正向載荷拉伸，考慮拉伸過程中的裂紋擴展情況。兩種模型的大小、尺寸相同，使用相同的EAM勢函數進行單向載荷加載，得到的應力應變曲線、楊氏模量及屈服應力如圖所示：

首先，圖(c)和圖(D)分別展示了基于圖(a)和圖(b)模型的應力-應變曲線。這兩條曲線直觀地反映了材料在受到外力作用下的力學響應。從應力-應變曲線中我們可以看出，盡管兩種模型具有不同的初始裂紋形態，但它們對單晶鋁的屈服應力影響并不顯著。這意味著在裂紋擴展之前，材料的彈性變形階段和屈服點附近的力學行為是相似的，裂紋形態并不是決定屈服應力的主要因素。然而，當目光轉向楊氏模量這一參數時，情況有所不同。楊氏模量是描述材料在彈性階段對力的響應程度的物理量，它反映了材料抵抗彈性變形的能力。從圖(D)中我們可以觀察到，對于圖(b)所示的模型，其缺陷形態（即特定的初始裂紋形態）導致了楊氏模量的相對增大。這表明在這種裂紋形態下，材料在彈性階段對力的響應更為敏感，需要更大的力才能使材料發生相同的彈性變形。

Ovito可視化圖：

在單晶鋁材料的裂紋擴展研究中，不同的裂紋形態在相同的勢函數和加載速率下會展現出顯著的差異。這些差異不僅影響裂紋擴展的速率，還直接關系到裂紋擴展的寬度，這在材料的力學性能和壽命評估中具有重要意義。在拉伸過程中，當單晶鋁受到外力作用時，裂紋會開始擴展。在這個過程中，裂紋的擴展寬度是衡量裂紋擴展程度的一個重要指標。在相同的應變條件下，即材料受到相同程度的拉伸變形時，不同的裂紋形態會導致裂紋擴展寬度的明顯不同。圖(e)中的裂紋在拉伸過程中展現出了較大的擴展寬度d1，這意味著裂紋的擴展更為迅速和顯著。這可能是由于圖(e)中的裂紋形態更容易在拉伸過程中形成應力集中，從而促進裂紋的擴展。相比之下，圖(f)中的裂紋在相同應變下的擴展寬度d2較小，這表明裂紋的擴展相對較慢，或者裂紋的擴展受到了某種程度的抑制。

圖(g)和圖(h)分別為圖(a)和(b)在拉伸載荷作用下的原子應力集中示意圖，值得注意的是，盡管預裂紋在圖(a)和圖(b)中的建立方式有所不同，但兩者在拉伸載荷下的原子應力集中現象都呈現出相似的模式。具體來說，無論是哪種預裂紋模型，應力集中的區域都主要集中在裂紋的四周。這一發現揭示了裂紋擴展過程中應力分布的一個普遍規律，也表明兩種建立方式并不影響模擬的準確性且值得探索其微觀變形機理。

最后，有相關分子動力學需求，歡迎通過公眾號“320科技工作室”與我們聯絡