CN109830267B - 一种螺旋状位错线原子结构的建模方法 - Google Patents
一种螺旋状位错线原子结构的建模方法 Download PDFInfo
- Publication number
- CN109830267B CN109830267B CN201910118388.8A CN201910118388A CN109830267B CN 109830267 B CN109830267 B CN 109830267B CN 201910118388 A CN201910118388 A CN 201910118388A CN 109830267 B CN109830267 B CN 109830267B
- Authority
- CN
- China
- Prior art keywords
- dislocation line
- atomic structure
- spiral
- file
- dislocation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Crystals, And After-Treatments Of Crystals (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种螺旋状位错线原子结构的建模方法。该方法的主要内容包括,在给定包含晶体模型原子结构信息的文件的前提下,根据拟创建螺旋状位错线原子结构的Burgers矢量、位错线位置、滑移面和位错线形态的要求,使用编程语言提取文件中晶体模型的原子结构信息,自动计算出包含符合要求的螺旋状位错线原子结构的晶体模型的原子坐标,然后按计算机仿真技术能识别的文件格式输出数据到文件。本发明可方便快捷地在晶体内部指定位置直接创建指定方位、组态和形态的螺旋状位错线原子结构,为计算机仿真技术对螺旋状位错线的形态及行为的精准研究创造了有利条件。
Description
技术领域
本发明涉及计算机仿真技术领域,尤其涉及一种螺旋状位错线原子结构的建模方法。
背景技术
晶体的塑性变形、晶体生长、形变强化、滞弹性、断裂、相变、晶体的电磁性能、晶体的光学性质、超导性以及其他许多物理、化学性质都与位错有重要的关联。因此,位错的研究无论对于科学研究还是实际应用都有着重要意义。实验上研究位错的方法有浸蚀法、缀饰法、透射电子显微分析法、X射线衍射分析法、场离子显微分析法等。这些实验技术被广泛地应用于分析研究位错的密度、分布和组态以及它们的运动和交互作用等。但是在原子尺度上的研究,计算机仿真技术(如分子动力学模拟)具有重要的地位。各种位错原子结构的直接创建有利于计算机仿真技术对位错行为更为精准的研究。本发明公开了一种螺旋状位错线原子结构的建模方法,解决了计算机仿真研究中直接创建螺旋状位错线原子结构的问题,本发明公开的方法可方便快捷地在晶体模型内部指定位置直接创建指定方位、组态和形态的螺旋状位错线原子结构。
发明内容
本发明所要解决的技术问题,就是提供一种方便快捷地创建螺旋状位错线原子结构的方法,本方法在给定包含晶体模型的原子结构信息的文件的前提下,根据拟创建螺旋状位错线原子结构的Burgers矢量、位错线位置、滑移面和位错线形态的要求,使用编程语言提取文件中晶体模型的原子结构信息,自动计算出包含符合要求的螺旋状位错线原子结构的晶体模型的原子坐标,然后按计算机仿真技术能识别的文件格式输出数据到文件。
本发明采用的技术解决方案如下:
设拟创建的螺旋状位错线原子结构的螺旋轴线经过点P(xp yp zp),位错的Burgers矢量为[uvw]a,a为晶格常数,滑移面为(hkl),螺距为h,旋转半径为r。
步骤一:准备包含晶体模型原子结构信息的文件。
步骤二:使用编程语言提取上述文件中的原子结构信息,将坐标系原点移动到点P处。旋转坐标系,使x轴的正向与[uvw]方向一致,y轴的正向与[hkl]方向一致,然后计算晶体模型内所有原子在新坐标系中的坐标值。
步骤三:设位错周围晶格扭曲较明显的区域在x和y方向的范围为以位错中心为中心的2b×2c的矩形。为创建螺旋状位错线原子结构,晶体模型内的原子要进行相应的位移,本方法根据螺旋状位错线周围原子分布的特点,提出以下计算位移的公式,设原子在x方向的位移为q,y方向和z方向不发生位移,计算公式如下:
设:d=a(u2+v2+w2)1/2,x’=x-r cos(2πz/h),y’=y-r sin(2πz/h),
当-b≤x’≤b且0≤y’≤c时,q=-d x’(1-y’/c)/(4b),
当x’>b且0≤y’≤c时,q=-d(1-y’/c)/4,
当x’<-b且0≤y’≤c时,q=d(1-y’/c)/4,
当y’>c时,q=0,
当-b≤x’≤b且-c≤y’<0时,q=d/2+d x’(1+y’/c)/(4b),
当x’>b且-c≤y’<0时,q=d/2+d(1+y’/c)/4,
当x’<-b且-c≤y’<0时,q=d/2-d(1+y’/c)/4,
当y’<-c时,q=d/2。
步骤四:根据上面计算得到的每个原子的位移值q,计算晶体模型内所有原子位移后的坐标值,由此在晶体内部指定位置直接构建出位向、组态和形态符合指定要求的螺旋状位错线原子结构。
步骤五:按照步骤二的逆向移动坐标系,使坐标系恢复成原来的位向。
步骤六:按计算机仿真技术能识别的格式输出数据到文件。
以上内容即为本发明公开的螺旋状位错线原子结构的建模方法的主要内容。
本发明所公开的螺旋状位错线原子结构的建模方法,可方便快捷地在晶体内部指定位置直接创建指定方位、组态和形态的螺旋状位错线原子结构,有利于计算机仿真技术更为精准地对螺旋状位错线的形态及行为进行研究。
附图说明
图1是本发明实施实例中创建的不含螺旋状位错线的B2型NiAl金属间化合物40×40×40超胞的原子结构图。
图2是本发明实施实例中超胞内已创建螺旋状位错线原子结构的Ovito软件显示的位错识别图,箭头指示出位错的Burgers矢量的方向。
具体实施方式
下面结合实施例及附图对本发明作进一步详细说明,但本发明的实施方式不限于此,在不脱离本发明上述思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的范围内。
实施例:
本实例公开了一种螺旋状位错线原子结构的建模方法。本实例在一个B2型NiAl金属间化合物40×40×40超胞内创建螺旋轴线经过超胞的中心、位错的Burgers矢量为1/2[1-11]a、位错的滑移面为(011)、螺旋半径为2a、螺距为6a的螺旋状位错线原子结构,a为晶格常数。
步骤一:采用Materials Studio创建40×40×40超胞,如图1所示,然后以car的格式输出数据文件。
步骤二:使用C/C++语言提取上述文件中的原子结构信息,将坐标原点移到超胞的中心,旋转坐标系,使x轴正向沿着晶向[1-11],y轴正向沿着晶向[011],计算晶体模型内所有原子在新坐标系中的坐标值。
步骤三:设定螺旋状位错线周围晶格扭曲较明显的区域在x和y方向的范围为以位错中心为中心的矩形8a×8a,a为晶格常数,计算晶体超胞内所有原子在x方向的位移q,y方向和z方向不发生位移,主要程序代码如下:
for (i=0; i< total_no_atoms; i++) { // total_no_atoms为晶体模型内原子的总数
a=2.882; pi=3.1415926; //晶格常数和常数π
aa=4a; bb=4a; //晶格扭曲较明显的范围参数
r=2a; h=6a; //螺旋半径和螺距
u=1; v=-1; w=1; h=0; k=1; l=1; //晶向和晶面指数
d=a*sqrt(u*u+v*v+w*w)/2; //计算位错的Burgers矢量的长度
x1=atoms[i].x[0]; y1=atoms[i].x[1]; z1=atoms[i].x[2]; //提取第i个原子的坐标
x2=x1-r*cos(2*pi/h*z1); y2=y1-r*sin(2*pi/h*z1);
if ((x2>=-aa) &&(x2<=aa) && (y2>=0) && (y2<=bb)) {q=-d/4*x2/aa*(1-y2/bb);}
if ((x2>aa) && (y2<=bb) && (y2>=0)){q=-d/4*(1-y2/bb);}
if ((x2<-aa) && (y2<=bb) && (y2>=0)){q=d/4*(1-y2/bb);}
if ((y2>bb)) {q=0;}
if ((x2>=-aa) &&(x2<=aa) && (y2<0) && (y2>=-bb)) {q=d/2+d/4*x2/aa*(1+y2/bb);}
if ((x2>aa) && (y2>=-bb) && (y2<0)){q=d/2+d/4*(1+y2/bb);}
if ((x2<-aa) && (y2>=-bb) && (y2<0)){q =d/2-d/4*(1+y2/bb);}
if ((y2<-bb)) {q =d/2;}}。
步骤四:根据上面计算得到的每个原子的位移值q,计算晶体超胞内所有原子位移后的坐标值,由此构建出螺旋轴线经过超胞的中心、位错的Burgers矢量为1/2[1-11]a、位错的滑移面为(011)、螺旋半径为2a、螺距为6a的螺旋状位错线原子结构,程序代码如下:
for(i=0; i<total_no_atoms; i++){
atoms[i].x[0]+=q;}。
步骤五:按照步骤二的逆向移动坐标系,使坐标系恢复成原来的位向。
步骤六:按计算机仿真软件能识别的格式输出数据到文件。
由此,完成了符合指定要求的螺旋状位错线原子结构的建模。图2为使用Ovito软件的位错识别工具显示的上述过程所构建的在一个B2型NiAl金属间化合物40×40×40超胞内螺旋轴线经过超胞的中心、位错的Burgers矢量为1/2[1-11]a、位错的滑移面为(011)、螺旋半径为2a、螺距为6a的螺旋状位错线原子结构。
Claims (1)
1.一种螺旋状位错线原子结构的建模方法,其特征在于该方法在给定包含晶体模型的原子结构信息的文件的前提下,根据拟创建螺旋状位错线原子结构的Burgers矢量、位错线位置、滑移面和位错线形态的要求,使用编程语言提取文件中晶体模型的原子结构信息,自动计算出包含符合要求的螺旋状位错线原子结构的晶体模型的原子坐标,然后按计算机仿真技术能识别的文件格式输出数据到文件,主要步骤如下:
步骤一:准备包含晶体模型原子结构信息的文件;
步骤二:使用编程语言提取上述文件中的原子结构信息,设拟创建的螺旋状位错线原子结构的螺旋轴线经过点P(xp yp zp),位错的Burgers矢量为[uvw]a,a为晶格常数,滑移面为(hkl),螺距为h,旋转半径为r;将坐标系原点移动到点P处,旋转坐标系,使x轴的正向与[uvw]方向一致,y轴的正向与[hkl]方向一致,然后计算晶体模型内所有原子在新坐标系中的坐标值;
步骤三:设位错周围晶格扭曲较明显的区域在x和y方向的范围为以位错中心为中心的2b×2c的矩形;为创建螺旋状位错线原子结构,晶体模型内的原子要进行相应的位移,本方法根据螺旋状位错线周围原子分布的特点,提出以下计算位移的公式,设原子在x方向的位移为q,y方向和z方向不发生位移,计算公式如下:
设:d=a(u2+v2+w2)1/2,x’=x-r cos(2πz/h),y’=y-r sin(2πz/h),
当-b≤x’≤b且0≤y’≤c时,q=-d x’(1-y’/c)/(4b),
当x’>b且0≤y’≤c时,q=-d(1-y’/c)/4,
当x’<-b且0≤y’≤c时,q=d(1-y’/c)/4,
当y’>c时,q=0,
当-b≤x’≤b且-c≤y’<0时,q=d/2+d x’(1+y’/c)/(4b),
当x’>b且-c≤y’<0时,q=d/2+d(1+y’/c)/4,
当x’<-b且-c≤y’<0时,q=d/2-d(1+y’/c)/4,
当y’<-c时,q=d/2;
步骤四:根据上面计算得到的每个原子的位移值q,计算晶体模型内所有原子位移后的坐标值,由此在晶体内部指定位置直接构建出位向、组态和形态符合指定要求的螺旋状位错线原子结构;
步骤五:按照步骤二的逆向移动坐标系,使坐标系恢复成原来的位向;
步骤六:按计算机仿真技术能识别的格式输出数据到文件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910118388.8A CN109830267B (zh) | 2019-02-16 | 2019-02-16 | 一种螺旋状位错线原子结构的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910118388.8A CN109830267B (zh) | 2019-02-16 | 2019-02-16 | 一种螺旋状位错线原子结构的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109830267A CN109830267A (zh) | 2019-05-31 |
CN109830267B true CN109830267B (zh) | 2021-08-24 |
Family
ID=66863624
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910118388.8A Active CN109830267B (zh) | 2019-02-16 | 2019-02-16 | 一种螺旋状位错线原子结构的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109830267B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109192249A (zh) * | 2018-10-24 | 2019-01-11 | 辽宁石油化工大学 | 一种以刃型位错为轴的正弦波形位错原子结构的建模方法 |
CN109215744A (zh) * | 2018-10-24 | 2019-01-15 | 辽宁石油化工大学 | 一种正弦波形位错攀移原子结构的建模方法 |
CN109243542A (zh) * | 2018-10-24 | 2019-01-18 | 辽宁石油化工大学 | 一种以螺型位错为轴的正弦波形位错原子结构的建模方法 |
CN109271740A (zh) * | 2018-10-21 | 2019-01-25 | 辽宁石油化工大学 | 一种螺型位错分叉原子结构的建模方法 |
CN109300509A (zh) * | 2018-10-11 | 2019-02-01 | 辽宁石油化工大学 | 一种构建刃型位错原子结构的方法 |
CN109299786A (zh) * | 2018-10-11 | 2019-02-01 | 辽宁石油化工大学 | 一种构建混合型位错原子结构的方法 |
CN109326331A (zh) * | 2018-10-11 | 2019-02-12 | 辽宁石油化工大学 | 一种位错攀移原子结构的建模方法 |
CN109346132A (zh) * | 2018-10-17 | 2019-02-15 | 辽宁石油化工大学 | 一种构建螺型位错原子结构的方法 |
CN109344509A (zh) * | 2018-10-04 | 2019-02-15 | 辽宁石油化工大学 | 一种螺型位错原子结构构建的方法 |
-
2019
- 2019-02-16 CN CN201910118388.8A patent/CN109830267B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109344509A (zh) * | 2018-10-04 | 2019-02-15 | 辽宁石油化工大学 | 一种螺型位错原子结构构建的方法 |
CN109300509A (zh) * | 2018-10-11 | 2019-02-01 | 辽宁石油化工大学 | 一种构建刃型位错原子结构的方法 |
CN109299786A (zh) * | 2018-10-11 | 2019-02-01 | 辽宁石油化工大学 | 一种构建混合型位错原子结构的方法 |
CN109326331A (zh) * | 2018-10-11 | 2019-02-12 | 辽宁石油化工大学 | 一种位错攀移原子结构的建模方法 |
CN109346132A (zh) * | 2018-10-17 | 2019-02-15 | 辽宁石油化工大学 | 一种构建螺型位错原子结构的方法 |
CN109271740A (zh) * | 2018-10-21 | 2019-01-25 | 辽宁石油化工大学 | 一种螺型位错分叉原子结构的建模方法 |
CN109192249A (zh) * | 2018-10-24 | 2019-01-11 | 辽宁石油化工大学 | 一种以刃型位错为轴的正弦波形位错原子结构的建模方法 |
CN109215744A (zh) * | 2018-10-24 | 2019-01-15 | 辽宁石油化工大学 | 一种正弦波形位错攀移原子结构的建模方法 |
CN109243542A (zh) * | 2018-10-24 | 2019-01-18 | 辽宁石油化工大学 | 一种以螺型位错为轴的正弦波形位错原子结构的建模方法 |
Non-Patent Citations (2)
Title |
---|
分子动力学中晶体位错的三维可视化系统设计与实现;宁媛等;《电脑知识与技术》;20140630;第10卷(第16期);第3896-3899页 * |
石墨烯中位错的结构与滑移;姚茵;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20180615(第6期);B014-3 * |
Also Published As
Publication number | Publication date |
---|---|
CN109830267A (zh) | 2019-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109192249A (zh) | 一种以刃型位错为轴的正弦波形位错原子结构的建模方法 | |
CN109215744A (zh) | 一种正弦波形位错攀移原子结构的建模方法 | |
CN109299786A (zh) | 一种构建混合型位错原子结构的方法 | |
CN109271740A (zh) | 一种螺型位错分叉原子结构的建模方法 | |
CN109243542A (zh) | 一种以螺型位错为轴的正弦波形位错原子结构的建模方法 | |
CN109346132A (zh) | 一种构建螺型位错原子结构的方法 | |
CN109326331A (zh) | 一种位错攀移原子结构的建模方法 | |
CN109300509A (zh) | 一种构建刃型位错原子结构的方法 | |
CN109344509A (zh) | 一种螺型位错原子结构构建的方法 | |
Gaiotto et al. | Aspects of $\Omega $-deformed M-theory | |
WO1982003705A1 (en) | Method of forming curved surface | |
CN109448112B (zh) | 三维模型库的建立方法及装置 | |
CN113988449A (zh) | 基于Transformer模型的风电功率预测方法 | |
CN109766667B (zh) | 刃位错和螺位错垂直相接的位错线原子结构的建模方法 | |
CN109830268B (zh) | 一种倾扭混合型晶界原子结构的建模方法 | |
CN109560849A (zh) | 一种基于二阶锥规划的波束赋形的鲁棒性自适应算法 | |
CN109830267B (zh) | 一种螺旋状位错线原子结构的建模方法 | |
CN109411024B (zh) | 一种位错环原子结构的建模方法 | |
CN109840381B (zh) | 一种蜷线位错原子结构的建模方法 | |
ODA | On a technique to obtain an optimum strength shape by the finite element method | |
CN109411025B (zh) | 一种混合型位错分叉原子结构的建模方法 | |
CN109817286B (zh) | 一种以刃型位错为轴的方波形位错线原子结构的建模方法 | |
Qi et al. | Folding dynamics and its intermittency in turbulence | |
Tao et al. | Chaotic time series prediction based on radial basis function network | |
CN109920488B (zh) | 一种以螺型位错为轴的方波形位错线原子结构的建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20211020 Address after: 110169 door 1, No. 99 A-12, Chuangxin first road, Hunnan District, Shenyang City, Liaoning Province Patentee after: Liaoning zhongkeliler Testing Technology Service Co.,Ltd. Address before: Wanghua District West Dandong Road 113001 Liaoning city of Fushun province No. 1 Patentee before: Liaoming Petrochemical University |
|
TR01 | Transfer of patent right |