24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 4366  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wuzhihua1988

新虫 (初入文坛)

[求助] fix deform 指令拉伸求教 已有1人参与

面是写的IN文件,但是模拟得出的结果,一直到应变为0.5时才将近拉断,而且感觉其中应力一度还出现了特别大的值,所以有哪位用过fix deform指令拉伸的还望指教一下?
#lammps graphene  test
units           metal
boundary        p s s  
atom_style      atomic
neighbor 2.0 nsq
neigh_modify delay 1 check yes
read_data       data.graphene
pair_style      tersoff                    
pair_coeff      * * SiCGe.tersoff C
region          1  block  INF 3  INF INF INF INF units box      
group           left region 1
region          2  block  61 INF  INF INF INF INF units box
group           right  region  2
group           boundary union left right
group           mobile subtract all boundary
velocity        mobile create 0.01 4928459 dist gaussian units box
velocity        boundary set 0.0 0.0 0.0 units box
#fix
fix    1 all nvt temp 0.01 0.01 0.001
fix    2 boundary setforce 0.0 0.0 0.0
min_style    sd
minimize    1.0e-12 1.0e-12 1000 1000
#run
timestep    0.001
compute    2 all pe/atom
compute    3 all reduce sum c_2
variable    avepo equal c_3/atoms
compute    4 all stress/atom
compute    strec all reduce sum c_4[1]
variable    s1 equal c_strec/(vol*10000)    #X方向应力
thermo_style custom step temp etotal vol v_avepo v_s1 lx
thermo          1000
run   100000
dump    1 all custom 1000 dump.out id type x y z
#for loop
label    forloop
variable    i loop 400
fix 5 all deform 1000 x erate 0.0015 units box
run    1000
variable    strainx equal lx/62.480000003-1    #X应变
print "@ @ @ @ : ${strainx}  ${s1} "
next    i
jump    in.graphene forloop

拉伸后的VMD图像
回复此楼

» 收录本帖的淘帖专辑推荐

力学性能

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

我是菜鸟,所以要早起
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinhanqing

金虫 (著名写手)


chaizhm: 金币+1, 谢谢~ 2012-06-17 20:35:45
fix deform命令不能应用于周期性边界条件
4楼2012-06-17 20:32:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

shrekren

禁虫 (小有名气)

★ ★
感谢参与,应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-06-01 10:25:37
本帖内容被屏蔽

2楼2012-05-31 00:49:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuzhihua1988

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by shrekren at 2012-05-31 00:49:08
先顶一下,是不是势函数或者模型的原因

谢谢了,势函数改了个CH.airebo的但感觉也不是很理想
我是菜鸟,所以要早起
3楼2012-06-01 09:56:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天泊乐

金虫 (正式写手)

【答案】应助回帖


chaizhm: 金币+1, 谢谢~ 2013-08-08 09:29:37
Any box dimension varied by this fix must be periodic.fix deform命令应该用于周期性边界
音乐you我!
5楼2013-08-07 11:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307中医考研调剂 +3 于以采蘩 2026-04-14 3/150 2026-04-15 11:16 by xtt..
[考研] 271求调剂 +31 2261744733 2026-04-11 31/1550 2026-04-15 10:50 by shenshen207
[考研] 302分求调剂 +10 凡语祈愿 2026-04-08 11/550 2026-04-14 16:50 by jiangguiquan11
[考研] 0854调剂 +11 长弓傲 2026-04-12 14/700 2026-04-14 16:07 by 逆水乘风
[考研] 297,工科调剂? +6 河南农业大学-能 2026-04-14 6/300 2026-04-14 16:04 by Art1977
[考研] 297,工科调剂?河南农业大学本科 +9 河南农业大学-能 2026-04-14 9/450 2026-04-14 16:03 by Art1977
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 机械工程313分找工科调剂 +4 双一流本科机械 2026-04-08 4/200 2026-04-14 07:32 by Abskk
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 材料考研调剂 +29 云木达达 2026-04-11 31/1550 2026-04-13 13:32 by lyh鲁老师
[考研] 一志愿华工085600 331分 +7 天下ww 2026-04-09 7/350 2026-04-13 09:01 by lhj2009
[考研] 22408 352分求调剂 +5 努力的夏末 2026-04-09 5/250 2026-04-12 19:17 by wj165256
[考研] 322求调剂 +6 123安康 2026-04-12 13/650 2026-04-12 15:51 by 123安康
[考研] 280求调剂 +7 兮兮夜夜 2026-04-09 10/500 2026-04-12 00:33 by 蓝云思雨
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
信息提示
请填处理意见