| 查看: 1116 | 回复: 6 | ||
Myprayert木虫 (正式写手)
|
[求助]
各位大神看看这个脚本哪里错了，十分感谢 已有2人参与
|
|
###################################### # Input file for dislocation nucleation in single crystal Al # Mark Tschopp, July 2010 # lmp_exe < Al_SC_100.in ###################################### ###################################### # INITIALIZATION units metal dimension 3 boundary p p p atom_style atomic ###################################### # ATOM BUILD lattice fcc 4.05 region whole block 0 10 0 10 0 10 create_box 1 whole region sc block INF INF INF INF INF INF units box lattice fcc 4.05 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 create_atoms 1 region sc variable natoms equal "count(all)" print "# of atoms are: ${natoms}" ###################################### # INTERATOMIC POTENTIAL pair_style eam/alloy pair_coeff * * Al99.eam.alloy Al # Define computes required compute csym all centro/atom ###################################### # EQUILIBRATION reset_timestep 0 timestep 0.001 velocity all create 300 12345 mom yes rot no fix 1 all npt 300 300 1 xyz 0 0 1 drag 1 # Set thermo output thermo 100 thermo_style custom step lx ly lz press pxx pyy pzz pe temp # Run for at least 10 picosecond (assuming 1 fs timestep) run 10000 unfix 1 # Store final cell length for strain calculations variable tmp equal "lx" variable L0 equal ${tmp} print "Initial Length, L0: ${L0}" ###################################### # DEFORMATION reset_timestep 0 fix 1 all npt 300 300 1 aniso NULL NULL 0 0 0 0 1 drag 1 variable srate equal 1.0e9 variable srate1 equal "v_srate / 1.0e12" fix 2 all deform 1 x erate ${srate1} units box remap x # Output strain and stress info to file # for units metal, pressure is in [bars] = 100 [kPa] = 1/10000 [GPa] # p2, p3, p4 are in GPa variable strain equal "(lx - v_L0)/v_L0" variable p1 equal "v_strain" variable p2 equal "-pxx/10000" variable p3 equal "-pyy/10000" variable p4 equal "-pzz/10000" fix def1 all print 100 "${p1} ${p2} ${p3} ${p4}" file Al_SC_100.def1.txt screen no # Dump coordinates to file (for void size calculations) dump 1 all custom 2500 dump.* x y z c_csym # Display thermo thermo 100 thermo_style custom step v_strain temp v_p2 v_p3 v_p4 ke pe press run 150000 ###################################### # SIMULATION DONE print "All done" |
回复此楼» 猜你喜欢
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
今年审到国自然15份，谈谈感受
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
-
半夜喝咖啡
已经有5人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
2楼2017-11-29 17:32:00
Myprayert
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2671.4
- 散金: 760
- 红花: 1
- 帖子: 537
- 在线: 158.2小时
- 虫号: 3701611
- 注册: 2015-02-28
- 性别: GG
- 专业: 固体力学
3楼2017-11-29 19:26:32
jluchemlib
铁杆木虫 (正式写手)
- 应助: 119 (高中生)
- 金币: 11607.4
- 散金: 200
- 红花: 19
- 帖子: 521
- 在线: 356.1小时
- 虫号: 589656
- 注册: 2008-08-30
- 性别: GG
- 专业: 高分子物理与高分子物理化
4楼2017-11-30 08:45:13
5楼2017-11-30 09:07:18
Myprayert
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2671.4
- 散金: 760
- 红花: 1
- 帖子: 537
- 在线: 158.2小时
- 虫号: 3701611
- 注册: 2015-02-28
- 性别: GG
- 专业: 固体力学
6楼2017-11-30 19:49:22
Myprayert
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2671.4
- 散金: 760
- 红花: 1
- 帖子: 537
- 在线: 158.2小时
- 虫号: 3701611
- 注册: 2015-02-28
- 性别: GG
- 专业: 固体力学
7楼2017-11-30 19:51:43
