| 查看: 1427 | 回复: 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
» 猜你喜欢
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
国自然申请面上模板最新2026版出了吗？
已经有11人回复
-
推荐一本书
已经有12人回复
-
基金申报
已经有4人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有4人回复
-
纳米粒子粒径的测量
已经有7人回复
-
常年博士招收(双一流，工科)
已经有4人回复
-
参与限项
已经有5人回复
-
有没有人能给点建议
已经有5人回复
|
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
3楼2017-04-12 08:09:08