24小时热门版块排行榜    

查看: 2579  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 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
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见