南宁seo网站建设,山东企业网站建设推荐,建站推荐网站,湖南优化推广本期给大家带来lammps模拟带缺陷镍板剪切变形的in文件及其详细注释。初始模型如图一所示#xff1a;图1 生成的初始模型 in文件及注释如下#xff1a;#利用eam势函数模拟带缺陷镍板的剪切#模型构成——上下镍板夹可动镍块#xff0c;镍块中有圆柱形缺陷#xff0c;移动上镍… 本期给大家带来lammps模拟带缺陷镍板剪切变形的in文件及其详细注释。初始模型如图一所示 图1 生成的初始模型 in文件及注释如下#利用eam势函数模拟带缺陷镍板的剪切#模型构成——上下镍板夹可动镍块镍块中有圆柱形缺陷移动上镍板使其发生剪切变形# 3d metal shear simulationunits metal #单位制为metalBoundary s s p #剪切和厚度方向为非周期性浮动边界宽度方向为周期性边界条件atom_style atomiclattice fcc 3.52 #面心立方晶格晶格常数3.52region box block 0 16.0 0 10.0 0 2.828427 #盒子总体大小长16埃米宽10厚2.828427create_box 3 boxlattice fcc 3.52 orient x 1 0 0 orient y 0 1 1 orient z 0 -1 1 #从原点开始沿上述的晶体取向生长create_atoms 1 box#在盒子内填满原子pair_style eampair_coeff * * Ni_u3.eam#调用eam势函数neighbor 0.3 binneigh_modify delay 5region lower block INF INF INF 0.9 INF INF #定义底板区域region upper block INF INF 6.1 INF INF INF#定义顶板区域group lower region lowergroup upper region uppergroup boundary union lower upper#顶板顶板设为边界群group mobile subtract all boundary#除去边界的地方为可动群set group lower type 2#设置底板原子类型为2set group upper type 3 #设置顶板原子类型为3# voidregion void cylinder z 8 3.535534 2.5 INF INF#以(83.535534)为轴线2.5埃米为半径定义出圆柱区域delete_atoms region void#删除圆柱区域的原子# temp controllerscompute new3d mobile temp#定义温度的计算new3d(mobile区域统计平均)compute new2d mobile temp/partial 0 1 1#定义温度的计算new2d,忽略x方向速度# equilibratevelocity 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.001psrun 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 3fix 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 100thermo_modify temp new2dreset_timestep 0 #重新设当前为第0步run 3000 #总共运行3000步 最后发生剪切后模型如下所示 图2 发生剪切后的模型最后如果您有分子动力学相关需求可以通过微信公众号联系我们呦~微信公众号320科技工作室
