[發明專利]一種構建刃型位錯原子結構的方法

在審

• 申請號: 201811183635.4
• 申請日: 2018-10-11
• 公開（公告）號: CN109300509A
• 公開（公告）日: 2019-02-01
• 發明（設計）人: 呂柏林;張藝程;蘇輝;高學朋;吳敏
• 申請（專利權）人: 遼寧石油化工大學
• 主分類號: G16C10/00
• 分類號: G16C10/00;G16C20/30
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種構建刃型位錯原子結構的方法，其特征在于該方法在給定包含晶體的原子結構信息的文件的前提下，使用編程語言提取晶體的原子結構信息，自動計算出在晶體內部包含擬構建的刃型位錯的原子坐標，然后按分子動力學軟件能識別的文件格式輸出文件，包含以下主要步驟：

步驟一：準備包含晶體模型原子結構信息的文件；

步驟二：使用編程語言提取上述文件中的原子結構信息，將坐標系原點移動到擬構建的刃型位錯線上；旋轉坐標系，使x軸在位錯線的滑移方向上，y軸垂直于位錯滑移面，然后計算晶體模型內所有原子在新坐標系中的坐標值；

步驟三：設定刃型位錯周圍晶格扭曲明顯的區域在x和y方向的范圍為2a×2b，根據下面的計算公式計算晶體模型內所有原子在x方向的位移q，y方向和z方向不發生位移，計算公式如下：

當-a≤x≤a且0≤y≤b時，q=-d/4×x/a×(1-y/b)，

當x>a且0≤y≤b時，q=-d/4×(1-y/b)，

當x<-a且0≤y≤b時，q=d/4×(1-y/b)，

當y>b時，q=0，

當-a≤x≤a且-b≤y<0時，q=d/2+d/4×x/a×(1+y/b)，

當x>a且-b≤y<0時，q=d/2+d/4×(1+y/b)，

當x<-a且-b≤y<0時，q=d/2-d/4×(1+y/b)，

當y<-b時，q=d/2，

其中，d為擬創建位錯的Burgers矢量的長度；

步驟四：根據上面計算得到的每個原子的位移值，計算晶體模型內所有原子位移后的坐標值，由此構建出相應位向的刃型位錯；

步驟五：按照步驟二的逆向移動坐標系，使坐標系恢復成原來的位向，為新的操作做準備；

步驟六：按分子動力學軟件能識別的格式輸出數據到文件。