24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2539  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工329分求调剂 +19 叶zilin 2026-04-13 20/1000 2026-04-14 00:23 by 福农资环_环境基
[考研] 332求调剂 +15 蕉蕉123 2026-04-10 15/750 2026-04-13 23:12 by pies112
[考研] 考研调剂 +11 长弓傲 2026-04-13 12/600 2026-04-13 22:48 by pies112
[考研] 复试调剂 +17 积极向上; 2026-04-10 19/950 2026-04-13 19:53 by 闲人终南山
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 266调剂 +10 daya sun 2026-04-07 11/550 2026-04-13 10:12 by fenglj492
[考研] 药学求调剂 +3 RussHu 2026-04-12 4/200 2026-04-12 17:49 by 陈皮皮
[硕博家园] 新一代电子信息294求调剂 不挑学校 +7 Ytyt11 2026-04-09 8/400 2026-04-12 16:57 by ajpv风雷
[考研] 322求调剂 +6 123安康 2026-04-12 13/650 2026-04-12 15:51 by 123安康
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 农学0904 312求调剂 +3 Say Never 2026-04-11 3/150 2026-04-11 17:22 by daydayup2005
[考研] 269求调剂 +11 啊啊我我 2026-04-07 11/550 2026-04-11 16:45 by vgtyfty
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 一志愿华东师范生物学326分,求调剂 +8 刘墨墨 2026-04-09 8/400 2026-04-10 12:00 by pengliang8036
[考研] 一志愿中科院105500专业总分315求调剂 +6 lallalh 2026-04-09 7/350 2026-04-09 17:51 by lallalh
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 招收有机化学、化工,药学,食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
信息提示
请填处理意见