应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » fix deform 指令拉伸求教
111/2返回列表12下一页
查看: 4314  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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图像
回复此楼

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

力学性能

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
我是菜鸟，所以要早起
1楼 2012-05-28 11:35:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

天泊乐

金虫 (正式写手)

【答案】应助回帖

★
chaizhm: 金币+1, 谢谢~ 2013-08-08 09:29:37
Any box dimension varied by this fix must be periodic.fix deform命令应该用于周期性边界
赞一下(1人) 回复此楼
音乐you我！
5楼2013-08-07 11:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

韩欢欢

金虫 (初入文坛)
引用回帖:
7楼: Originally posted by ldx626 at 2014-03-09 22:35:29
请问，c_strec/(vol*10000)是什么意思呢？10000是怎么来的？

单位换算，1GPa=10000bar
赞一下(1人) 回复此楼
Demi
8楼2014-03-10 09:32:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

天泊乐

金虫 (正式写手)

【答案】应助回帖

★
chaizhm: 金币+1, 谢谢~ 2013-08-08 09:29:47
非周期边界条件使用displace atom命令
赞一下 回复此楼
音乐you我！
6楼2013-08-07 11:49:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldx626

木虫 (初入文坛)
请问，c_strec/(vol*10000)是什么意思呢？10000是怎么来的？
赞一下 回复此楼
7楼2014-03-09 22:35:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dyc_2008

捐助贵宾 (著名写手)

【答案】应助回帖

1、应力很大是正常的，这个模拟出来的，是理论强度，动辄以Gpa为单位。
2、deform，是均匀变形，所以延伸率会非常的大。有次我运行一个测试，3个方向都是周期性，变形到了300%还没断。
赞一下 回复此楼
一布福田，二修菩提
9楼2014-03-10 09:57:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldx626

木虫 (初入文坛)
引用回帖:
8楼: Originally posted by 韩欢欢 at 2014-03-10 09:32:12
单位换算，1GPa=10000bar...

那62.480000003-1是怎么来的？
赞一下 回复此楼
10楼2014-03-10 10:22:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuzhihua1988 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +4 vv迷 2026-03-22 4/200 2026-03-22 23:29 by king123！
[考研] 352求调剂 +3 大米饭！ 2026-03-22 3/150 2026-03-22 23:28 by king123！
[考研] 0854电子信息求调剂 +3 α____ 2026-03-22 3/150 2026-03-22 21:28 by zhq0425
[考研] 北科281学硕材料求调剂 +6 tcxiaoxx 2026-03-20 6/300 2026-03-22 20:23 by edmund7
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 一志愿华中科技大学，080502，354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
关闭
请填处理意见
关闭