| 查看: 1545 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
LAMMPS模拟均聚物拉伸，采用dpd、srp、msrp方法得到的应力应变图都是趋近0 已有1人参与
|
|||
in文件中，第一个是dpd，第二个是msrp。求指点 链长40-3000原子-dpd.jpg  40-msrp.jpg  文献图.gif |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : data.chain
- 附件 2 : in.dpd
- 附件 3 : in.msrp
2017-02-28 15:21:42, 179.07 K
2017-02-28 15:26:44, 1.5 K
2017-02-28 15:27:01, 1.8 K
» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
3楼2017-04-12 08:09:08
|
in.dpd如下 variable simname index nc10_c11000 units lj atom_style full boundary p p p special_bonds lj/coul 1 1 1 newton on on # save an extra atom type for bond particles read_data data.chain neighbor 2.0 bin neigh_modify every 10 check yes bond_style harmonic bond_coeff * 225.0 0.85 comm_modify vel yes # must use pair hybrid, since srp bond particles # do not interact with other atoms types pair_style dpd 1.0 1.0 373692 pair_coeff 1 1 60.0 4.5 1.0 thermo 1000 thermo_style custom step temp pe press atoms lx ly lz dump 1 all custom 10000 dump1.lammpstrj id type x y z fix 1 all nve timestep 0.01 run 200000 unfix 1 ##################################################### # Uniaxial Tensile Deformation run 0 variable tmp equal "lx" variable L0 equal ${tmp} variable strain equal "(lx - v_L0)/v_L0" variable p1 equal "v_strain" variable p2 equal "-pxx/10000*1.01325" variable p3 equal "-pyy/10000*1.01325" variable p4 equal "-pzz/10000*1.01325" variable p5 equal "lx" variable p6 equal "ly" variable p7 equal "lz" variable p8 equal "temp" fix 1 all npt temp 1 1 1 y 0 0 1000 z 0 0 1000 drag 2 fix 2 all deform 1 x erate 1e-5 units box remap x fix def1 all print 1000 "${p1} ${p2} ${p3} ${p4} ${p5} ${p6} ${p7} ${p8}" file ${simname}.def1.txt screen no dump 2 all custom 10000 dump2.lammpstrj id type x y z thermo_style custom step temp pxx pyy pzz lx ly lz thermo 1000 timestep 0.01 reset_timestep 0 run 500000 unfix 2 unfix def1 |
2楼2017-02-28 15:24:56
投票:
10