lammps学习（一）单晶硅纳米磨削

学习lammps写的第一个简单的例子。参考文献《单晶硅纳米级磨削过程的理论研究》。

一、in文件部分

# 金刚石半径4.5nm，晶格常数3.57
# 单晶硅 x：15nm,y：15nm ,z：15nm#echo screen #输出调试、仿真到屏幕###设定基本参数units metal #金属类型
dimension 3 #三维
atom_style atomic #无键
boundary f f f #非边界条件###生成原子/区域
region box block -200 200 -200 200 -200 200 units box #box区域，长宽高1000的盒子。
create_box 2 box #创建包含两种原子的盒子lattice diamond 3.57 #金刚石模型，晶格常数3.57
region s1 sphere 60 -35 63 30 units box #区域名称s1 球形区域 球心坐标xyz 球半径4.5 单位与box相同
region s2 block 0 100 -100 0 0 63 units box #区域名s2 长方体区域 xyz的最高最低坐标 单位同box
region se intersect 2 s1 s2 #区域名 联合区域 联合数 联合s1 联合s2#create_box 2 se #原子类型数目 区域create_atoms 1 region se units box #原子类型 区域 单位###创建多层Si
lattice diamond 5.4
region b1 block  0 120 0 88.23 39.2 58.82 units box
region b2 block  0 120 0 88.23 19.60 39.2 units box
region b3 block  0 120 0 88.23 0 19.60 units box
region b4 union 3 b1 b2 b3create_atoms 2 region b4 units box###设置质量
mass 1 12.0 #C质量12.0g/mol
mass 2 28.0 #Si质量28g/mol###分组
group JGS type 1 #将第一类原子归为JGS类
neighbor 2.0 bin #skin距离
neigh_modify every 2 delay 10 check yes exclude group JGS JGS
#在上一次构建list1后20个步长内不进行构建，延时后两个步长才开始构建，并且要满足当某个原子超过skin距离的一般。排除JGS group的相互作用。group Si type 2 #将2类原子归为JGS类
group layer1 region b1 #将b1归为layer1类（牛顿层）
group layer2 region b2 #恒温层
group layer3 region b3 #固定边界层###设置作用势
pair_style tersoff #按tersoff计算相互作用势
pair_coeff * * SiC.tersoff C Si  #两种原子类型，1C，2Si
###设置计算方法
compute mytemp layer2 temp/com #按temp/com计算恒温层温度#nve下驰豫
velocity JGS set NULL 0.1667 NULL #JGS组vx设置为1，其余不设置
velocity layer2 create 293 12345 dist gaussian #在293K温度下以12345随机种子赋予layer2，且原子速度按高斯分布
fix 10 JGS nve/noforce #nve：保证系统中原子数量n、体系总体积v和体系总能量e保持不变。noforce：忽略原子力，不更新速度，更新位置。fix 20 layer1 nve
fix 30 layer2 nve
fix 31 layer2 temp/rescale 1 293 293 1.0 1.0 #基本原理是当温度超出设定温度时，通过缩放的方式调节温度使其达到设定值。
fix_modify 31 temp mytemp
fix 40 layer3 nve/noforce  ###输出体系参数
thermo 200 #每200步输出结果
thermo_style custom step temp c_mytemp ke pe etotal enthalpy #自定义输出格式
thermo_modify lost ignore###输出原子参数
dump 3 JGS custom 500 velocity.txt vx vy vz
dump 2 Si custom 1000 Force.txt id x y z fx fy fz
dump 1 all custom 500 C.lammpsstrj id type x y z #dumpID 所有组 自定义输出形式 每一步输出一次 输出坐标xyztimestep 0.001
run 60000 #run一个步长write_data C-Si.data #输出data文件
write_restart C-Si.restart

二、ovito磨削过程可视化

三、数据处理

脚本自动化处理还没怎么学，自己手动拉excel将所有Si中所有原子的力求的和输出的受力—时间曲线。

lammps学习（一）单晶硅纳米磨削相关推荐

圆锥形SiC纳米磨削单晶硅源代码
本文以圆锥形SiC磨粒磨削单晶Si为例,介绍纳米磨削模拟过程以及力的输出. 模拟流程: (1)建模:单晶硅分为三个区域:边界区.恒温区.牛顿区.磨粒为圆锥形,材料为SiC. (2)模拟初始化,固定边界 ...
lammps学习（二）联合原子模型聚乙烯拉伸
相较于之前做的单晶硅纳米磨削,聚合物的建模要复杂的多,用lammps自带密令通常无法解决.常用的建模方法有用python/C/MATLAB代码生成聚合物data文件,不过就算建模个简单的聚合物往往也需 ...
matlab模拟砂轮形貌,基于多颗磨粒随机分布的虚拟砂轮建模及磨削力预测
[1] Syoji K. Grinding technology[M]. Guo Y B, Wang Z Z, translated. Beijing: China Machine Press, 20 ...
用ANSYS画矩形_钢轨打磨用复合砂轮磨削温度场的研究
转载自金刚石与磨料磨具工程闫薪霖, 肖博, 吴恒恒, 肖冰, 张益权(南京航空航天大学机电学院, 南京 210016)摘要为研究钢轨打磨用钎焊金刚石插片复合砂轮磨削时的温度场,用复合砂轮和树脂 ...
青春痘疤的治疗,微晶磨削 TCA FREXEL (转)
我长痘痘已经5年了,现在情况好多了可留下了不少痘坑从长痘痘起就被灌输了错误的思想不可以挤痘所以我从来不敢碰痘痘的等它自己平了之后可怕的是留了一个硬块在皮肤里过几天脱落后就形成了一个小坑 ...
LAMMPS学习系列（1）
LAMMPS学习系列(1) LAMMPS模拟软件的功能有很多: 微观结构演化结构优化过渡态搜索热导率计算等等(太多了,有兴趣的可自行百度)..... 其中它最常用的功能就是模拟微观结构的演化, ...
磨削工具类毕业论文文献有哪些？
本文是为大家整理的磨削工具主题相关的10篇毕业论文文献,包括5篇期刊论文和5篇学位论文,为磨削工具选题相关人员撰写毕业论文提供参考. 1.[期刊论文]磨削工具抑尘装置的研究与应用期刊:<中国环 ...
德国jager电主轴在牙科义齿加工，PCB加工，机器人雕铣磨削应用简介！
德国Jager电主轴应用行业德国Alfred Jager电主轴目前主要应用于宝石加工.牙科.机械表壳加工.玻璃磨边倒角.轮胎模具刻字刻花.金属切削.大幅面雕刻.毛刺修整.金属高光.内圆磨削.PCB板 ...
用计算机计算磨削,磨削参数计算软件
磨削参数计算软件是用于计算磨削线速度的工具,输入砂轮直径的参数后,就能计算出砂轮线速度和砂轮转速. 相关软件软件大小版本说明下载地址磨削参数计算软件是用于计算磨削线速度的工具,输入砂轮直径的参数后, ...

lammps学习（一）单晶硅纳米磨削

一、in文件部分

二、ovito磨削过程可视化

三、数据处理

lammps学习（一）单晶硅纳米磨削相关推荐

最新文章

热门文章