| 查看: 4911 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
清风明月P6铜虫 (初入文坛)
|
[求助]
求一个LAMMPS里面聚合物适用的拉伸in文件
|
|
|
最近在做拉伸模拟。自己尝试着写了一个in文件，但是运行之后确得不到想要的结果。。。 我用聚乙烯做的测试，现在MS里面建好的模型，后来到处转化为lammps文件 在论坛找了好久基本上都是金属的，求大神给个高分子可以用的吧。。谢谢了 下面是我的in文件，也不知道问题出在哪里。。附件是得到的输出文件 #------模拟环境初始化------ boundary p p p #边界周期性设定，周期、周期、自由非周期 units real #模型单元模式设定为real dimension 3 #几何维度设定为3D atom_style full #原子模式设定为full neighbor 2.0 bin #设定neighbor截断半径-力场截断半径（类似MS中Buffer width） neigh_modify every 1 delay 0 check yes #------微观模型的建立------ # 定义原子结构 pair_style lj/class2/coul/cut 10.0 9.5 #airebo 3.0 1 0 lj/cut/coul/cut 10.0 #设定PCFF势函数的对势（类似MS中的Nonbond potential）信息：对势包含相互作用为范德华力和库仑力，范德华力和库仑力截断半径（类似MS中Cutoff distance） bond_style class2 #设定Compass势函数的键势 angle_style class2 #设定Compass势函数的角势 dihedral_style class2 #设定Compass势函数的两面角势 improper_style class2 #设定Compass势函数的两面角扭转势 read_data ../bench/20140924/PeDeform.lammps05 #------设置原子相关参数------ velocity all create 300 102486 rot yes dist & gaussian units box #设定所有原子的速率 # 能量优化 minimize 1.0e-7 0 1000000 1000000 #设置minimize，指定能量收敛值为1.0e-7、力矢量在坐标轴上的分量长度收敛值为0、最大迭代步数为1000000、计算力长度/能量的最大次数为1000000 #------动力学设置------ fix NPT all npt temp 300.0 300.0 100 iso 1.0 1.0 1000 #对整体设定npt系综，控温方法为Nose/Hoover，设定起始、终止温度和压力，衰减常数（温度一般取100倍timesteps，这里取100.0，压力一般取1000倍timesteps，这里取1000.0） #拉伸模拟 reset_timestep 0 region 1 block INF INF INF INF INF 3 units box group low region 1 region 2 block INF INF INF INF 72 INF units box group up region 2 group boundary union low up group mobile subtract all boundary velocity up set 0 0 0.004 velocity low set 0 0 -0.004 # 动力学输出设置 thermo 1000 #每1000步输出一次热动力学信息（温度、能量、压力） thermo_style custom step lx ly lz press pxx pyy pzz pe temp thermo_modify lost warn #检查丢失原子，显示警报 #fix myrelax all ave/time 1 1000 1000 c_lz c_pxx ../bench/20140924/data/relax_gpc.dat1 timestep 1 #设定时间步长为1fs dump train1 all xyz 1000 ../bench/20140924/g-PEc.xyz #每1000步输出所有区间的dcd坐标，写入到g_PE_C.xyz文件中 run 100000 #运行10ps |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : myout.out
2014-09-28 17:19:56, 11.55 K
» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有4人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
什么时候开奖？
已经有10人回复
-
2026国自然函评费到账
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- MS粗粒化模型如何导入lammps
已经有10人回复
- LAMMPS进行拉伸模拟的一些问题
已经有4人回复
- 如何将amber的轨迹文件mdcrd格式转换成gromacs的轨迹文件格式 trr呢
已经有5人回复
- 新手求助一个quantum-espresso的使用问题！
已经有13人回复
- lammps 中 做拉伸时 sss和ssp边界条件不同，结果为什么不一样？？？？
已经有5人回复
- 硅纳米管用lammps做拉伸
已经有3人回复
- 请教一个C++读文件的问题，
已经有12人回复
- lammps做拉伸和压缩的应力问题
已经有9人回复
- 跑完Gromacs的动力学模拟之后的轨迹文件.trr里面没有内容，为什么？
已经有4人回复
- 求回答聚合物被拉伸时，在屈服点和应变软化以及应变硬化时内部分子链是怎么运动的的？
已经有9人回复
- lammps实验拉伸Al纳米线的in文件的错误提示好久找不出？
已经有8人回复
- 关于lammps fix npt/deform做拉伸时温度的变化
已经有6人回复
- lammps 拉伸应力应变曲线
已经有4人回复
- windows下lammps编译
已经有5人回复
- 【求助】DL_POLY的config文件如何生成??
已经有98人回复
南gege
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 540.7
- 散金: 10
- 帖子: 45
- 在线: 50.3小时
- 虫号: 3807463
- 注册: 2015-04-13
- 性别: GG
- 专业: 高分子材料结构与性能
3楼2020-10-10 10:12:23
tang_wei
金虫 (正式写手)
- 应助: 24 (小学生)
- 金币: 1354.3
- 红花: 28
- 帖子: 428
- 在线: 81.7小时
- 虫号: 4067192
- 注册: 2015-09-12
- 专业: 高分子物理与高分子物理化
2楼2017-12-08 14:30:25
40