24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4391  |  回复: 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的回帖

ldx626

木虫 (初入文坛)

请问,c_strec/(vol*10000)是什么意思呢?10000是怎么来的?
7楼2014-03-09 22:35:29
已阅   回复此楼   关注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的回帖

yinhanqing

金虫 (著名写手)


chaizhm: 金币+1, 谢谢~ 2012-06-17 20:35:45
fix deform命令不能应用于周期性边界条件
4楼2012-06-17 20:32:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料相关专业344求调剂双非工科学校或课题组 +24 hualkop 2026-04-12 26/1300 2026-04-18 23:33 by 路病情
[考研] 291求调剂 +10 关忆北. 2026-04-14 10/500 2026-04-18 23:32 by 路病情
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 化工学硕294分,求导师收留 +33 yzyzx 2026-04-12 37/1850 2026-04-17 23:00 by wunaiy88
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 335求调剂 +20 想上岸呀!! 2026-04-12 23/1150 2026-04-17 10:50 by cuisz
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 5/250 2026-04-17 10:02 by bobvan
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 297,工科调剂? +10 河南农业大学-能 2026-04-14 10/500 2026-04-15 21:50 by noqvsozv
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
信息提示
请填处理意见