| 查看: 6306 | 回复: 28 | |||||||||
| 本帖产生 2 个 模拟EPI ，点击这里进行查看 | |||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||
[交流]
【讨论】想学lammps是否应该先学习C++ 已有16人参与
|
|||||||||
| 最近打算学习lammps，是不是应该先把C++学好呢？还是只要懂点基础就行了？ |
回复此楼» 收录本帖的淘帖专辑推荐
 MD分子动力学 | Lammps 相关 | 分子动力学 | 计算 |
| 正能量 | LAMMPS |
» 猜你喜欢
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在Lammps中安装xmovie工具时，Makefile文件怎么修改X库路径
已经有3人回复
- lammps tutorial
已经有7人回复
- lammps中固定键长
已经有4人回复
- 【求助】LAMMPS-运行报错
已经有9人回复
- 【分享】LAMMPS的语法中文解释
已经有137人回复
- 【求助】fedora中LAMMPS编译问题
已经有4人回复
- 【分享】尝试lammps中, 分享中...
已经有15人回复
- 【求助】LAMMPS、NAMD和CHARMM谁最适合生物大分子模拟？
已经有8人回复
- 【其他】免费软件学习资料大征集（abinit, pwscf, lammps, CPMD, XMD 等等免费软件）
已经有21人回复
jiangtao9571
新虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 238.2
- 散金: 1046
- 红花: 1
- 帖子: 454
- 在线: 236.2小时
- 虫号: 706996
- 注册: 2009-02-23
- 专业: 高分子材料结构与性能
3楼2011-04-13 16:27:23
2楼2011-01-22 17:54:54
jiangtao9571
新虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 238.2
- 散金: 1046
- 红花: 1
- 帖子: 454
- 在线: 236.2小时
- 虫号: 706996
- 注册: 2009-02-23
- 专业: 高分子材料结构与性能
4楼2011-04-13 16:27:44
御剑江湖
荣誉版主 (著名写手)
天池冶海
- 模拟EPI: 15
- 应助: 7 (幼儿园)
- 贵宾: 6.831
- 金币: 9082.9
- 散金: 8688
- 红花: 129
- 沙发: 4
- 帖子: 2917
- 在线: 1251.7小时
- 虫号: 1229656
- 注册: 2011-03-11
- 专业: 电机与电器
- 管辖: 第一性原理
★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+6, 模拟EPI+1): 谢谢，精彩的答疑 2011-04-14 09:56:37
小木虫: 金币+2, 帖子真精彩 2014-03-19 12:38:20
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+6, 模拟EPI+1): 谢谢，精彩的答疑 2011-04-14 09:56:37
小木虫: 金币+2, 帖子真精彩 2014-03-19 12:38:20
|
其实我觉得你刚学lammps，先看看网上关于in文件的详细解释步骤，再学习Fortran语言，最后才C++，因为那时就是自己编译开发源程序，楼上“懂一点基础比较好，尤其是建立模型和数据处理的时候用途很大”也是正确的，不过建立模型，现在一般用MS，然后read data比较方便，数据处理用origin也不错，呵呵，这只是我一点建议而已。 给你看看关于in文件的详细解释步骤 boundary p s s # 边界条件，拉伸方向是周期性，其余是自由边界；如果是薄膜拉伸则是两个周期性，块体则是三个周期性 units metal #单位制定义为metal atom_style atomic #原子类型自动 neighbor 2.0 bin #截断半径相关的东西 neigh_modifydelay 1 check yes # 邻近原子列表更新速度 #create geometry lattice fcc 3.61 # 定义晶胞为fcc，晶格常数3.61A region box block 0 30 0 3 0 3 #定义一个长方体区域叫box，长30，宽和高是3 create_box 1 box #创建了这样一个box create_atoms 1 box #在box里创建了一种原子 mass 1 63.546 #定义这种原子的质量是63.546 # potentials pair_style eam # 定义势函数是EAM pair_coef * * Cu_u3.eam #势所需要的参数在此文件里 #define groups region 1 block INF 1 INF INF INF INF # 定义了一个叫1的区域 group left region 1 #定义此区域里的原子叫left region 2 block 29 INF INF INF INF INF #定义了一个叫2的区域 group right region 2 #定义此区域里的原子叫right group boundary union left right #定义left+ right = boundary group mobile subtract all left # 定义mobile= all - left # initialvelocities velocity left set 0.0 0.0 0.0 # 设置原子初速度为0 computep all pressure thermo_temp # 计算应力，计算结果记为p variable pressx equal c_p[1] #定义变量pressx=c_p[1],c_p[1]的意思是p里第一个值 variable pressy equal c_p[2] #定义变量pressy=c_p[2],c_p[2]的意思是p里第二个值 variable pressz equal c_p[3] #定义变量pressz=c_p[3],c_p[3]的意思是p里第三个值 thermo_style custom step temp etotal press v_pressx v_pressy v_pressz vol dump 1 all atom 1000 dump.tensile # 输出结果到dump.tensile minimize 1.0e-6 1.0e-6 1000 1000 # 驰豫 fix 1 left setforce 0.0 NULL NULL #固定左边的原子 fix 2 all npt 1.0 1.0 1.0 aniso 0.0 0.0 NULL NULL NULL NULL 10.0 drag 1.0 thermo 1000 # 每1000步输出一次结果 timestep 0.002 #时间步长0.002ps run 10000 #运行10000步 unfix2 fix 3 all nve fix 4 all temp/rescale 100 1.0 1.0 5 1.0 fix 5 all deform 1000 x erate 0.002 units box #设置沿x方向拉伸，应变率为0.002(1/ps) compute s all stress/atom #计算每个原子上的应力，计算结果记为s compute strea all reduce sum c_s[1] #将s里第一个值求和，计算结果记为strea variable stressa equal c_strea/vol # 定义变量stressa= strea/vol, vol是体积 compute streb all reduce sum c_s[2] variable stressb equalc_streb/vol compute strec all reduce sum c_s[3] variable stressc equal c_strec/vol thermo_style custom step tempetotal press v_stressa v_stressb v_stressc vol run 300000 |
5楼2011-04-13 23:45:03
40