应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【求助】求一模拟程序（fortran）
141/2返回列表12下一页
查看: 963  |  回复: 13
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

jianchaoyv2651

铁虫 (初入文坛)

[交流] 【求助】求一模拟程序（fortran）

本人初次编程解决实际问题，感觉力不从心！！请求各位高手帮忙！
问题如下：用fortran语言编程即时计算和展现十个氢原子在模拟盒中的运动情况。初速度满足麦克斯韦分布，采取硬球完全弹性碰撞

[ Last edited by lei0736 on 2009-11-25 at 11:58 ]
回复此楼

» 猜你喜欢
1楼2009-04-12 15:12:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyy816

至尊木虫 (著名写手)
★ ★ ★
jianchaoyv2651(金币+1,VIP+0): 4-12 15:21
lei0736(金币+2,VIP+0):谢谢 4-20 17:32
这个就是相当于编写一个单原子分子动力学模拟的源程序啊。。。
建议楼主去参考computer simulation of liquids里面的例子
赞一下(13人) 回复此楼
2楼2009-04-12 15:18:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianchaoyv2651

铁虫 (初入文坛)
能否给说一下具体参考《computer simulation of liquids》中那一部分？？
赞一下 回复此楼
3楼2009-04-12 15:26:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianchaoyv2651

铁虫 (初入文坛)
能否给说一下具体参考《computer simulation of liquids》中那一部分？？
赞一下 回复此楼
4楼2009-04-12 15:32:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州
★ ★ ★ ★ ★
jianchaoyv2651(金币+3,VIP+0): 4-12 16:02
lei0736(金币+2,VIP+0):谢谢 4-20 17:33
http://www.mdbbs.org/viewthread.php?tid=6423
这里有个模拟256个Ar原子，LJ势的供学习用的frotran代码以及一些解释。

H原子因为强的量子力学效应不适合MD
赞一下(13人) 回复此楼
5楼2009-04-12 15:38:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianchaoyv

金虫 (小有名气)
jianchaoyv2651(金币+11,VIP+0): 4-13 11:39
lei0736(金币+0,VIP+0):马甲行为 自己给自己转金币 请解释！！！ 5-7 11:58
lei0736(金币-11,VIP+0):请解释 5-9 13:32
问题已解决！！
回复此楼
6楼2009-04-13 11:36:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianchaoyv

金虫 (小有名气)
lei0736(金币+0,VIP+0):请履行承诺奖励金币给解决问题的人 5-7 11:59
引用回帖:
Originally posted by bay__gulf at 2009-4-12 15:38:
http://www.mdbbs.org/viewthread.php?tid=6423
这里有个模拟256个Ar原子，LJ势的供学习用的frotran代码以及一些解释。

H原子因为强的量子力学效应不适合MD

请问所给的模拟256个Ar原子程序中函数RNDF的作用是什么？
调用子程序MXWELL的作用是什么？

[ Last edited by jianchaoyv on 2009-5-7 at 10:33 ]
赞一下(10人) 回复此楼
7楼2009-05-07 10:28:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianchaoyv

金虫 (小有名气)
请看清楚不是一个人的，只不过我们认识吧了
赞一下 回复此楼
8楼2009-05-07 14:02:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州
★ ★
lei0736(金币+2,VIP+0):谢谢 5-9 13:30
引用回帖:
Originally posted by jianchaoyv at 2009-5-7 10:28:
请问所给的模拟256个Ar原子程序中函数RNDF的作用是什么？
调用子程序MXWELL的作用是什么？
[ Last edited by jianchaoyv on 2009-5-7 at 10:33 ]

RNDF生成随机数函数
MAXWELL用于给原子赋速度初值
赞一下(10人) 回复此楼
9楼2009-05-07 14:11:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianchaoyv

金虫 (小有名气)
请问我是按所给程序运行的，所得结果怎样转化为我在其他文献上看到此程序所得的总能量随步数演化的图形、模拟中动能的分布图以及其他一些物理量的值等，请问应该怎样去操作？
下面是我运行结果的前面的一部分，后面是从100到3000部分我未写出
MOLECULAR DYNAMICS SIMULATION PROGRAM

NUMBER OF PARTICLE IS             256
SIDE LENGTH OF THE BOX IS   6.752840   
CUT OFF IS                  2.500000   
REDUCED TEMPERATURE IS     0.7220000   
BASIC TIME STEP IS         6.4000003E-02

TOTAL LINEAR MOMENTUM IN X DIRECTION IS  6.7786448E-02
TOTAL LINEAR MOMENTUM IN Y DIRECTION IS -6.6548891E-02
TOTAL LINEAR MOMENTUM IN Z DIRECTION IS  0.1201641   

VELOCITY ADJUSTMENT
TEMPERATURE BEFORE SCALING IS   45.41020   
SCALE FACTOR IS  0.1260933   
        2     270.799927   -1699.370728   -1428.570801
        4     252.771790   -1681.374146   -1428.602295
        6     218.253845   -1646.222656   -1427.968750
        8     171.139252   -1598.697998   -1427.558716
       10     132.214249   -1559.054932   -1426.840698
       12     116.545746   -1543.839722   -1427.293945
       14     119.052315   -1546.832153   -1427.779785
       16     127.827530   -1555.922119   -1428.094604
       18     134.569000   -1562.949707   -1428.380737
       20     137.285446   -1565.422729   -1428.137329
       22     137.336853   -1565.384644   -1428.047852
       24     136.563614   -1564.932617   -1428.369019
       26     134.083801   -1562.593628   -1428.509766
       28     129.720490   -1558.355103   -1428.634644
       30     126.987198   -1555.609863   -1428.622681
       32     127.910660   -1556.748169   -1428.837524
       34     131.486847   -1560.345459   -1428.858643
       36     135.715729   -1564.671997   -1428.956299
       38     138.538544   -1567.548828   -1429.010254
       40     139.642487   -1568.478882   -1428.836426
       42     140.040283   -1568.885742   -1428.845459
       44     140.190887   -1569.073975   -1428.883057
       46     139.712982   -1568.837280   -1429.124268
       48     138.952988   -1568.239868   -1429.286865
VELOCITY ADJUSTMENT
TEMPERATURE BEFORE SCALING IS  0.3634446   
SCALE FACTOR IS   1.409449   
       50     276.164978   -1568.357178   -1292.192139
       52     271.862701   -1564.260986   -1292.398315
       54     259.602783   -1552.116089   -1292.513306
       56     242.570480   -1535.268799   -1292.698364
       58     225.266388   -1518.148560   -1292.882202
       60     211.759369   -1504.485107   -1292.725708
       62     205.603088   -1498.357300   -1292.754150
       64     209.019333   -1501.927979   -1292.908691
       66     213.457825   -1506.205566   -1292.747803
       68     215.538666   -1508.412109   -1292.873413
       70     217.871109   -1510.706299   -1292.835205
       72     223.228973   -1516.197266   -1292.968262
       74     229.282425   -1522.213135   -1292.930664
       76     230.661316   -1523.383911   -1292.722656
       78     227.155609   -1519.937866   -1292.782227
       80     221.844513   -1514.623657   -1292.779175
       82     215.889771   -1508.402832   -1292.513062
       84     211.333252   -1503.536133   -1292.202881
       86     211.333344   -1503.780518   -1292.447144
       88     216.985687   -1509.352051   -1292.366333
       90     227.358459   -1519.550293   -1292.191895
       92     235.516296   -1527.578247   -1292.062012
       94     236.259933   -1528.119873   -1291.859985
       96     232.384796   -1524.102417   -1291.717651
       98     227.971985   -1519.898438   -1291.926514
VELOCITY ADJUSTMENT
TEMPERATURE BEFORE SCALING IS  0.5878155   
SCALE FACTOR IS   1.108276   
      100     276.164978   -1516.690308   -1240.525391
赞一下 回复此楼
10楼2009-05-07 16:42:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jianchaoyv2651 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭