| 查看: 1535 | 回复: 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
» 猜你喜欢
能否退出参与的面上项目解除限项
已经有18人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有3人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
|
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
