[发明专利]一种混合型位错分叉原子结构的建模方法

摘要

本发明公开了一种混合型位错分叉原子结构的建模方法。该方法的主要内容包括，在给定包含晶体模型原子结构信息的文件的前提下，根据拟创建的混合型位错分叉原子结构各位错的Burgers矢量、位错线的方向、滑移面和分叉位置的要求，利用编程语言提取文件中的晶体模型原子结构信息，自动计算出包含符合要求的混合型位错分叉原子结构的晶体模型的原子坐标，然后按计算机仿真技术能识别的文件格式输出数据到文件。本发明可方便快捷地在晶体内部指定位置直接创建指定方位和组态的混合型位错分叉原子结构，为计算机仿真技术对混合型位错分叉原子结构的形态及位错间反应行为的精准研究创造了有利条件。

主权项

1.一种混合型位错分叉原子结构的建模方法，其特征在于该方法在给定包含晶体模型的原子结构信息的文件的前提下，根据拟创建混合型位错分叉原子结构各位错的Burgers矢量、位错线的方向、滑移面及分叉位置的要求，使用编程语言提取文件中晶体模型的原子结构信息，自动计算出包含符合要求的混合型位错分叉原子结构的晶体模型的原子坐标，然后按计算机仿真技术能识别的文件格式输出文件，主要步骤如下：步骤一：准备包含晶体模型原子结构信息的文件；步骤二：使用编程语言提取上述文件中的原子结构信息，设拟创建的混合型位错分叉原子结构的分叉点P的坐标为(xp yp zp)，三条位错的Burgers矢量分别为[u1v1w1]a、[u2v2w2]a和[u3v3w3]a，[u1v1w1]a + [u2v2w2]a = [u3v3w3]a，a为晶格常数，位错线的方向分别为[e1f1g1]、[e2f2w2]和[e3f3g3]，滑移面为(hkl)；将坐标系原点移动到点P处，旋转坐标系，使x轴的正向与[u1v1w1]方向一致，y轴的正向与[hkl]方向一致，然后计算晶体模型内所有原子在新坐标系中的坐标值；步骤三：设位错周围晶格扭曲较明显的区域在x和y方向的范围为以位错中心为中心的2b×2c的矩形；为创建混合型位错分叉原子结构，晶体模型内的原子要进行相应的位移，本方法根据混合型位错分叉结构周围原子分布的特点，提出两步位移法，首先计算第一步位移，设原子在x方向的位移为q1，y方向和z方向不发生位移，计算公式如下：d1=a (u12+v12+w12)1/2，z1=v1 l ‑ w1 k，z2=w1 h ‑ u1 l，z3=u1 k ‑ v1 h，b1= (e1 u1+f1 v1+g1 w1) (z12+z22+z32)1/2/((e1 z1+f1 z2+g1 z3) (u12+v12+w12)1/2)，b3= (e3 u1+f3 v1+g3 w1) (z12+z22+z32)1/2/((e3 z1+f3 z2+g3 z3) (u12+v12+w12)1/2)，当z≤0时，x’= x ‑ b1z，当z>0时，x’= x ‑ b3z，当‑b≤x’≤b且0≤y≤c时，q1=‑d1x’(1‑y/c)/(4b)，当x’>b且0≤y≤c时，q1=‑d1(1‑y/c)/4，当x’<‑b且0≤y≤c时，q1=d1(1‑y/c)/4，当y>c时，q1=0，当‑b≤x’≤b且‑c≤y<0时，q1=d1/2+d1 x’(1+y/c)/(4b)，当x’>b且‑c≤y<0时，q1=d1/2+d1(1+y/c)/4，当x’<‑b且‑c≤y<0时，q1=d1/2‑d1(1+y/c)/4，当y<‑c时，q1=d1/2；步骤四：根据上面计算得到的每个原子的位移值q1，计算晶体模型内所有原子位移后的坐标值；步骤五：旋转坐标系，使x轴的正向与[u2v2w2]方向一致，y轴的正向与[hkl]方向一致，计算晶体模型内所有原子在新坐标系中的坐标值；步骤六：计算第二步位移，设在新坐标系中原子在x方向的位移为q2，y方向和z方向不发生位移，计算公式如下：d2=a (u22+v22+w22)1/2，z1=v2 l2 ‑ w2 k2，z2=w2 h2 ‑ u2 l2，z3=u2 k2 ‑ v2 h2，b2=(e2 u2+f2 v2+g2 w2) (z12+z22+z32)1/2/((e2 z1+f2 z2+g2 z3) (u22+v22+w22)1/2)，b3=(e3 u2+f3 v2+g3 w2) (z12+z22+z32)1/2/( (e3 z1+f3 z2+g3 z3) (u22+v22+w22)1/2)，当z≤0时，x’= x ‑ b2z，当z>0时，x’= x ‑ b3z，当‑b≤x’≤b且0≤y≤c时，q2=‑d2x’(1‑y/c)/(4b)，当x’>b且0≤y≤c时，q2=‑d2(1‑y/c)/4，当x’<‑b且0≤y≤c时，q2=d2(1‑y/c)/4，当y>c时，q2=0，当‑b≤x’≤b且‑c≤y<0时，q2=d2/2+d2x’(1+y/c)/(4b)，当x’>b且‑c≤y<0时，q2=d2/2+d2(1+y/c)/4，当x’<‑b且‑c≤y<0时，q2=d2/2‑d2(1+y/c)/4，当y<‑c时，q2=d2/2；步骤七：根据上面计算得到的每个原子的位移值q2，计算晶体模型内所有原子位移后的坐标值，由此创建出在指定位置分叉的位向和组态符合要求的混合型位错分叉原子结构；步骤八：按照步骤二和步骤五的逆向移动坐标系，使坐标系恢复成原来的位向；步骤九：按计算机仿真技术能识别的格式输出数据到文件。