24小时热门版块排行榜

返回列表

我在东洞庭边

金虫 (正式写手)

应助: 3 (幼儿园)
金币: 470.7
散金: 138
红花: 4
沙发: 4
帖子: 412
在线: 110小时
虫号: 4304430
注册: 2015-12-22
性别: GG
专业: 机械摩擦学与表面技术

[求助] 求助lammps tools 中的chain 调用？

lammps 额外提供了一些tools 文件，其中包含了我需要的可以构建高聚物分子链条单体的chain工具，可是现在遇到了一些问题。
我之前已经成功运行examples中的例子，已经有了一点思路，简单来说就是改变in文件中的一些参数，或者改变语句（删除%使之有效）再直接在终端运行in 文件，可是现在关于chain的文件不包括in文件。有一些疑惑：
1.现在我认为可能与chain工具有关的有三个 def.chain,data.file,chain.f.文件，而且其中的内容研读了好几遍，但是具体怎么使用呢？
2.打开chain .f 发现其中有一些代码感觉是不是要在别的地方运行啊，不在lammps中运行，如下所示。
   program chain

   integer swaptype
   integer, allocatable :: nchain(

,nmonomer(

integer, allocatable :: ntype(

,nbondtype(

integer, allocatable :: type(

,molecule(

integer, allocatable :: imagex(

,imagey(

,imagez(

real*8, allocatable :: x(

,y(

,z(

real*8, allocatable :: bondlength(

,restrict(

   common xprd,yprd,zprd,xboundlo,xboundhi,
   $    yboundlo,yboundhi,zboundlo,zboundhi
   real*8 random
900  format(a)
901  format(2f15.6,a)
902  format(i3,f5.1)
903  format(i10,i8,i8,3f10.4,3i4)
904  format(i9,i3,2i9)

3.想建立高聚合物模型，改变分子链分布（chain数量和monomer数量），在施加力的情况研究他的摩擦力，是不是不需要用这个chain工具呢？我会认真仔细听取别的思路。

真心向各位前辈求教，谢谢～！

回复此楼