• 回答数

    2

  • 浏览数

    1038

  • 收藏数

    0

作者:侧耳倾听 发表于 2018-11-19 17:38:39
跳转到指定楼层
这是我的IN文件,只有最后一点不一样
units metal
boundary s s s
atom_style atomic
lattice fcc 4.05
region box1 block 0 300 5 0 5
create_box 1 box1        # 创建盒子,名为box1
region box2 block 0 30 1.5 3.5 1.5 3.5 side out #创建空心体,名box2
create_atoms 1 region box2
pair_style eam/alloy
pair_coeff * * Al99.eam.alloy Al
neighbor 2.5 bin
neigh_modify delay 5
region   left block 29 INF INF INF INF INF
region   right blockINF 1 INF INF INF INF  #加载组定义
group  left region left
group  rightregion right
group  boundary unionleft right
group   mobile subtract all boundary    #group定义
#computenew1boundarystress/atom
compute new2all centro/atom  fcc     #原子适配度
compute new3 mobile temp          #统计温度
variable Pstress equal -pxx/10000       #XX应力分量,单位为Gpa
variable Pvelocity equal  1
variable Pdisplacement equal  step*v_Pvelocity/1000  #位移,单位为0.1埃(10的-11次方)
minimize 10e-8 10e-8 5000 5000         
dump 123 all cfg 500 aa*.cfg id type xs ys zs c_new2
minimize 10e-8 10e-8 5000 5000
      
fix 1 all nve
fix 2 left setforce0.0 NULL NULL
fix 3 rightsetforce0.0 NULL NULL
fix 4 mobile temp/rescale 10 30.0 30.0 5.0 1.0
fix_modify 4 temp new3
velocity mobile create 30 887723 temp new3
timestep 0.001  
thermo 1000      
run 1000          #弛豫1000步
reset_timestep 0       #时间步重置   
fix   stress boundary ave/time1 100 100 v_Pstress file Stress.txt
fix   strain boundary  ave/time  1 1 100 v_PdisplacementfileDisplacement.txt
velocityleft set v_Pvelocity 0 0   #位移加载
(常量情况就是把v_Pvelocity改成1)
run 30000
用后处理软件看的时候明显发现用常量的时候拉伸快很多,输出的应力分量也完全不一样

下面图中,拉伸明显较快的是常量1加载的
分享:
回复

使用道具

新手上路

Rank: 1

积分
29
极客币
53
主题
8
帖子
27
注册时间
2018-11-5
在线时间
1 小时
性别
保密
沙发
 楼主| 发表于 2018-11-19 17:39:12 | 只看该作者
这两张图是输出的盒子长度结果,发现用variable的时候,加载的位移是规律的1000步移动1的数值,而用常量1加载时则显得不规则,这是个什么情况?
回复

使用道具 举报

新手上路

Rank: 1

积分
29
极客币
47
主题
8
帖子
21
注册时间
2018-11-5
在线时间
1 小时
性别
保密
板凳
发表于 2018-11-19 17:39:48 | 只看该作者
我的car文件怎不能用表格打开呢?而且我的电脑目前没有能打开它的软件,请问你是用什么打开的?
还有,我把导出来的模型保存成.pdb的形式可以么?
回复

使用道具 举报

高级模式 评论
您需要登录后才可以回帖 登录 | 立即注册