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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 冲击结果异常
51/1返回列表
查看: 1778  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lchunmei001

木虫 (正式写手)

[求助] 冲击结果异常

请教大家，我沿x方向给粒子2KM/s的速度冲击加载，具体命令如下，但得到的结果很奇怪，压强很低的情况下，温度就已达到几万K，显然很不合理，但是在不知道是哪里出问题了，测试了几种物质，均出现这种情况，应该是命令哪里没对，还请高手指点一二，感激不尽！
velocity        all set 20  NULL NULL  sum yes units box
fix             1 all nve
fix             2 all wall/reflect xlo xhi
回复此楼

» 收录本帖的淘帖专辑推荐

condensed matter physics

» 猜你喜欢
1楼 2013-01-16 10:46:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lchunmei001: 金币+5 2013-01-16 16:23:13
zh1987hs: 金币+3, 谢谢 2013-01-19 14:04:12
楼主，我有点儿不懂，请教你哈：
1、你是要模拟一个什么样的冲击呢？我把你的in文件运行了一下，4万多个原子都一起向x方向运动吗？安装我的理解，所谓的“冲击”应该是一个物体以较高的速度撞击另一个物体，比如炮弹攻击掩体的时候。我看过使用lammps模拟纳米颗粒冲击石墨烯的。等等这些。我没弄明白你的模拟要干啥。
2、在你的in文件中，你能告诉我哪句命令是thermostat的？你设置的温度是多少？


PS 把dump命令改成  dump                id all atom 10000   dump.lammpstrj  吧，这样方便模拟完了使用VMD观看轨迹动画。
赞一下 回复此楼
5楼2013-01-16 15:20:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2013-01-16 14:29:01
把整个in文件，和相关输出的log贴出来，否则无法判断。
赞一下 回复此楼
2楼2013-01-16 11:11:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lchunmei001

木虫 (正式写手)
输入和结果如下，请予以指点，谢谢！
输入：
units                 metal
boundary        s p p
lattice         fcc 4.07  
atom_style        atomic
region                box block 0 100 0 10 0 10
create_box        1 box
create_atoms        1 box

pair_style eam
pair_coeff * * Pt_u3.eam

timestep        1.0e-3
neighbor        1.0 bin
neigh_modify    every 1 delay 10 check yes

thermo          100
thermo_style custom step press vol temp pxx pyy pzz lx ly lz

velocity        all set 20  NULL NULL  sum yes units box
fix             1 all nve
fix             2 all wall/reflect xlo xhi

dump                id all atom 10000   dump.melt
run             50000

结果log文件：

Step Press Volume Temp Pxx Pyy Pzz Lx Ly Lz
       0    19146.265    674326.27    31285.986    535234.08   -238897.64   -238897.64     407.0814         40.7         40.7
     100    23424.423    673363.58    31311.738     540451.5   -235089.11   -235089.11    406.50024         40.7         40.7
     200    27653.106    671960.25    31331.714    546177.73   -231609.21   -231609.21    405.65307         40.7         40.7
     300    31059.835    671084.06    31347.923    550838.79   -228829.64   -228829.64    405.12413         40.7         40.7
     400    36428.574    669579.59    31362.412    557393.91   -224054.09   -224054.09    404.21589         40.7         40.7
     500    39253.628    668750.35    31377.061    561471.39   -221855.25   -221855.25     403.7153         40.7         40.7
     600    44925.431     667277.7    31394.038    568443.51   -216833.61   -216833.61    402.82628         40.7         40.7
     700    47685.276    666310.63    31413.841     572717.7   -214830.93   -214830.93    402.24247         40.7         40.7
。。。。。。
赞一下 回复此楼
3楼2013-01-16 11:36:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
jiaoyixiong: 金币+3, 赞！鼓励交流 2013-01-16 15:10:23
lchunmei001: 金币+5 2013-01-16 16:23:08
引用回帖:
3楼: Originally posted by lchunmei001 at 2013-01-16 11:36:01
输入和结果如下，请予以指点，谢谢！
输入：
units                 metal
boundary        s p p
lattice         fcc 4.07  
atom_style        atomic
region                box block 0 100 0 10 0 10
create_box        1 box
create_atoms        1 box

p ...

你的原子初速度很高，原子从边界反弹会与后面的原子频繁撞击。 X方向的压力很高， 有可能产生瞬间高温。 看能量，你的过程还算平稳。 有几点可以尝试：
1）你用dump命令输出各个原子的速度值，看看其最大值是多少？
2） 盒子的X尺寸在变小，是你所需要的吗？ 边界用p p p试试？
赞一下 回复此楼
4楼2013-01-16 14:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +10 枫桥ZL 2026-03-18 12/600 2026-03-21 22:02 by peike
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 一志愿深大，0703化学，总分302，求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 277材料科学与工程080500求调剂 +6 自由煎饼果子 2026-03-16 6/300 2026-03-21 17:21 by 学员8dgXkO
[考研] 0856材料专硕353求调剂 +3 NIFFFfff 2026-03-20 3/150 2026-03-21 10:23 by luoyongfeng
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 15/750 2026-03-21 10:16 by 脱颖而出
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 一志愿华中科技大学，080502，354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕，总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 295材料求调剂，一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
信息提示
关闭
请填处理意见
关闭