24小时热门版块排行榜    

查看: 2526  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +5 Zhao anqi 2026-03-22 5/250 2026-03-22 17:38 by barlinike
[考研] 一志愿东华大学化学070300,求调剂 +6 2117205181 2026-03-21 7/350 2026-03-22 14:59 by allen-yin
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 384求调剂 +3 子系博 2026-03-22 4/200 2026-03-22 11:04 by 搏击518
[考研] 0856材料专硕353求调剂 +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
信息提示
请填处理意见