24小时热门版块排行榜    

查看: 2578  |  回复: 1

zzhwise

金虫 (小有名气)

[求助] 关于使用setforce施加压力相关问题 已有1人参与

我使用以下命令,拉动上平板沿着x方向运动,同时给其施加一个沿着z轴负方向的力,为什么它会穿过另一个平板
units            real                                                                              
dimension        3                                                                                   
boundary         p p p                                                                              
newton           on                                                                                 
atom_style       atomic                                                                                
                                                                             
lattice                       fcc 4.08   
region                       box block 0 11 0 11 -10  20
create_box             4 box   
region           upboundary  block 0  11    0  11  12 INF  units lattice
region           lowboundary  block 0  11    0  11  INF  -3   units lattice
region                       force block 0  11   0  11   9  12   units lattice                        
region                       lower block 0  11   0  11   -3  0   units lattice                                         
region                       upper block 0  11   0  11    7  9     units lattice

create_atoms     1  region  lower
create_atoms     1  region  force
mass                         1   196.096655
group                         lower region lower   
group                         upper region upper  
group                         force region force     
                 
group            Rigid   union force         #拉动组                                          
group            FREEZ   union lower                                                         
neigh_modify     every 1 delay  10  exclude group Rigid Rigid  #刚体内相互作用不考虑         
                                                                                             

pair_style       eam                 
                                                                                                               
pair_coeff       * *   Au_u3_real.eam  
                                                        
neighbor         0.2  bin    # rcutoff+0.2                                                           
neigh_modify     every 1 delay  10                                                                     
                                                                                                   
timestep               2                                                                                         
                                                                                                     
thermo                 200                                                                                      
                                                                                                     
#结构优化   
                                                                                   
min_style        sd                                                                                    
minimize          1.0e-4  1.0e-4   10000   1000                                                            
timestep                 2   

#  fix    80  all nve
velocity         upper create 300 12345 mom yes rot no
fix              10 upper nvt temp 300.0 300.0 100.0  #加温
fix              20 Rigid rigid group 1 Rigid force 1 off off on  #刚体,xy方向无力,z方向不取消力   
fix              21 Rigid move linear 0.1 0.0  NULL               #右拉
fix              22 Rigid setforce 0.0 0.0 -1.0                   #下压
fix              30 FREEZ  setforce 0.0 0.0 0.0                         #固定不动下版

dump 1 all xyz 500 model1.xyz   

run 10000

关于使用setforce施加压力相关问题
1.png


关于使用setforce施加压力相关问题-1
2.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iSimuLy

捐助贵宾 (正式写手)

资深专家顾问


【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
zzhwise: 金币+5, ★★★很有帮助 2017-04-27 19:24:38
查查你的相互作用吧
如果相互作用弱,重叠很正常
2楼2017-04-27 18:03:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zzhwise 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +8 ziyangfang 2026-08-25 11/550 2026-09-01 09:51 by zzuzxg
[基金申请] 面上函评意见出来了,像什么等级? 20+3 Tsingking1 2026-08-27 14/700 2026-09-01 08:29 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见