剪切文件_lammps模拟带缺陷镍板剪切变形(in文件及注释)
本期给大家带来lammps模拟带缺陷镍板剪切变形的in文件及其详细注释。初始模型如图一所示:
图1 生成的初始模型
in文件及注释如下:
#利用eam势函数模拟带缺陷镍板的剪切
#模型构成——上下镍板夹可动镍块,镍块中有圆柱形缺陷,移动上镍板使其发生剪切变形
# 3d metal shear simulation
units metal
#单位制为metal
Boundary s s p
#剪切和厚度方向为非周期性浮动边界,宽度方向为周期性边界条件
atom_style atomic
lattice fcc 3.52
#面心立方晶格,晶格常数3.52
region box block 0 16.0 0 10.0 0 2.828427
#盒子总体大小,长16埃米,宽10,厚2.828427
create_box 3 box
lattice fcc 3.52 orient x 1 0 0 orient y 0 1 1 orient z 0 -1 1
#从原点开始,沿上述的晶体取向生长
create_atoms 1 box
#在盒子内填满原子
pair_style eam
pair_coeff * * Ni_u3.eam
#调用eam势函数
neighbor 0.3 bin
neigh_modify delay 5
region lower block INF INF INF 0.9 INF INF
#定义底板区域
region upper block INF INF 6.1 INF INF INF
#定义顶板区域
group lower region lower
group upper region upper
group boundary union lower upper
#顶板+顶板设为边界群
group mobile subtract all boundary
#除去边界的地方为可动群
set group lower type 2
#设置底板原子类型为2
set group upper type 3
#设置顶板原子类型为3
# void
region void cylinder z 8 3.535534 2.5 INF INF
#以(8,3.535534)为轴线,2.5埃米为半径定义出圆柱区域
delete_atoms region void
#删除圆柱区域的原子
# temp controllers
compute new3d mobile temp
#定义温度的计算new3d(mobile区域统计平均)
compute new2d mobile temp/partial 0 1 1
#定义温度的计算new2d,忽略x方向速度
# equilibrate
velocity mobile create 300.0 5812775 temp new3d
#设置原子初始速度
fix 1 all nve
#nve系综
fix 2 boundary setforce 0.0 0.0 0.0
#固定边界原子,钢化原子,便于加载
fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0
#直接调温法,温度控制为300K+-10K之内
fix_modify 3 temp new3d
#fix 3 调温时采用new3d的计算方案
thermo 25
#每25步输出一次热力学统计量
thermo_modify temp new3d
#每次输出的温度按照new3d方案来
timestep 0.001
#时间步长0.001ps
run 100
#弛豫100步
# shear 剪切
velocity upper set 1.0 0 0
#使顶板原子以1埃/ps的速度沿x正向运动
velocity mobile ramp vx 0.0 1.0 y 1.4 8.6 sum yes
#mobile原子的初始速度从0到1线性变化
unfix 3
#结束弛豫时的温度约束fix 3
fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0
#开始新的控温,直接调温法,温度控制为300K+-10K之内
fix_modify 3 temp new2d
#fix 3 调温时采用new2d的计算方案
dump 1 all atom 100 dump.shear.void
#每隔100步输出一次原子信息
thermo 100
thermo_modify temp new2d
reset_timestep 0 #重新设当前为第0步
run 3000 #总共运行3000步
最后,发生剪切后模型如下所示:
图2 发生剪切后的模型
最后,如果您有分子动力学相关需求,可以通过微信公众号联系我们呦~
微信公众号:320科技工作室
剪切文件_lammps模拟带缺陷镍板剪切变形(in文件及注释)相关推荐
- 亲测好用的caj文件转换成带目录的文字版PDF文件工具
参考GitHub上caj2pdf/caj2pdf: Convert CAJ (China Academic Journals) files to PDF. 转换中国知网 CAJ 格式文献为 PDF.佛 ...
- zynq7020安装linux图形界面,ZYNQ开发板的NFS文件系统安装 - ZYNQ7020学习
ZYNQ开发板的NFS文件系统安装 - ZYNQ7020学习 ZYNQ开发板的NFS文件系统安装 - ZYNQ7020学习 这篇博客主要讲解NFS的安装和开发板与Linux虚拟机如何mount成功.重 ...
- 选择的串口 _ 不存在或开发板没有连接_PC 和开发板之间传输文件
8.1 传输文件的多种方法 在嵌入式开发中,涉及交叉编译:在 PC 上编译,在开发板上运行.所以就涉及 PC 和开发板之间的文件 传输,方法有多种:网络传输.串口传输.USB 传输,当然,还有最笨的方 ...
- tftpd32、虚拟机、PC机、开发板之间的文件传输
问题:如何在电脑虚拟机上的Linux C程序下载到开发板呢? 首先,虚拟机上安装Linux系统 Linux上的操作有: 编辑好Hello.c源程序 接着编译,因为最后是要在开发板执行程序,所以使用编译 ...
- [RK3288] AP6335WIFI模块移植(板级dts文件代码编写)
1.经过一周的努力的今天终于解决了RK3288安卓系统中打不开WIFI的问题,特意记录一下. 2. 一般拿到板子厂家会提供RK3288的SDK框架包(我的是android8.1),在这里我附上移植WI ...
- atitit. 文件上传带进度条 atiUP 设计 java c# php
atitit. 文件上传带进度条atiUP设计java c# php 1. 设计要求 1 2. 原理and架构 1 3. ui 2 4. spring mvc 2 5. springMVC.xml 3 ...
- i.MX6 u-boot 怎么确定板级头文件
/*********************************************************************** i.MX6 u-boot 怎么确定板级头文件* 说明: ...
- java web 默认页面配置文件_Tomcat中配置全局的错误页面(如404)+删除Tomcat中webapps目录下的自带项目,防止Tomcat默认文件泄露...
进入tomcat目录中conf文件编辑web.xml,将下面代码复制到文件末尾,如下 400 /error.html 404 /404/404.html 500 /404/500.html 2.添加4 ...
- Linux(shell)遍历目录删除指定文件,解决文件夹名称带空格问题
最近打算把一些学习资源解压上传到阿里网盘上,可是在解压的时候出现了问题,无故产生很多隐藏文件(猪哥mac电脑,用的keka解压软件). 且文件夹层级很多,不可能实现手动一个一个文件夹删除,所以就考虑写 ...
最新文章
- 数据科学之基石:数据科学家必须掌握的10个统计学概念
- 线程中使用SaveFileDialog不能弹出窗体
- 【职场】面试中最难回答频率最高的70个问题
- java王子救公主的游_计蒜客 王子救公主(DFS)
- 小孩都看得懂的推荐系统
- namefilter 前台反斜杠格式_001获取小猴子的信息并格式化输出
- forms、forms_toolbar例子理解
- 书法是什么?书法的美从何说起?
- 【ROS学习笔记】(五)话题消息的定义与使用
- 69. x 的平方根
- Unity中的存档与读档
- 月报总结|Moonbeam 11月份大事一览
- 提升领导力这7个法则,你不得不知道
- 经纬财富:乐山稳健投资之道,积小胜为大胜
- 射频原理图设计checklist
- WPF TabControl 美化
- 廖雪峰webApp部署
- 关于电脑WLAN突然消失
- 算法刷题-钥匙与房间
- 都说程序员的工资高,程序员工资真的很高吗?小编今天带你解析一下程序员薪资的真实情况