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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 有人做裂纹模拟的吗？看看我这个裂纹的in文件有什么问题
汕头大学海洋科学接受调剂
131/1返回列表
查看: 2628  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kite175

银虫 (初入文坛)

[交流] 有人做裂纹模拟的吗？看看我这个裂纹的in文件有什么问题

# 3d Morse crack simulation


units          metal
boundary        s s p
atom_style        atomic


# create geometry

lattice                fcc 4.05  orient x 1 0 0  orient y 0 1 0 orient z 0 0 1
region                box block 0 60 0 40 -0.25 0.25
create_box        5 box
create_atoms        1 box

mass                1 1.0
mass                2 1.0
mass                3 1.0
mass                4 1.0
mass                5 1.0


# Morse potentials

pair_style        morse 2.5
pair_coeff        * * 1.0 1.0 2.5
neighbor        0.3 bin
neigh_modify        every 5 delay 0 check yes

# define groups

region                1 block INF INF INF 2 INF INF
group                lower region 1
region                2 block INF INF 38 INF INF INF
group                upper region 2
group                boundary union lower upper
group                mobile subtract all boundary


region                leftupper block INF 20 20 INF INF INF
region                leftlower block INF 20 INF 20 INF INF

region          bonds   block  INF 20  19.75 20.25 -0.25 0.25
delete_atoms    region  bonds

group                leftupper region leftupper
group                leftlower region leftlower


group                halfupper1 region halfupper1
group                halfupper2 region halfupper2

set                group leftupper type 3
set                group leftlower type 4


set                group lower type 1
set                group upper type 2

neigh_modify        exclude type 3 4
neigh_modify        exclude type 1 2


# 能量最小化
min_style sd
minimize       1.0e-12 1.0e-12  10000 10000

# initial velocities 弛豫

compute                  new mobile temp
velocity        mobile create 0.001 887723 temp new
timestep        0.001
thermo                200
run 10000

velocity        upper  set 0.0 0.01 0.0
velocity        lower  set 0.0 -0.01 0.0
velocity        mobile ramp vy 0.01 0.01 y 20 38 sum yes

# fixes
fix            1 all nve  
fix            2 lower setforce  0.0 NULL 0.0
fix            3 upper setforce  0.0 NULL 0.0
# run
thermo                200
thermo_modify        temp new

dump                1 all atom 1000 dump.crack.lammpstrj
run                300000
回复此楼

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-08-27 09:17:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kite175

银虫 (初入文坛)
怎么没有人回答呀？纠结呀
赞一下 回复此楼
2楼2012-09-05 10:11:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1779064387

新虫 (初入文坛)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
jiaoyixiong: 金币+2, 鼓励交流 2012-09-27 19:49:16
建议把你的data文件贴一下~
赞一下 回复此楼
3楼2012-09-27 15:19:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
chaizhm: 金币+2, 谢谢~ 2012-09-28 09:38:30
引用回帖:
2楼: Originally posted by kite175 at 2012-09-05 10:11:52
怎么没有人回答呀？纠结呀

不知为何纠结？

想知道你的in文件有何问题，最简单的办法就是你运行一下啊，
有什么问题，立刻就知道了。
如果顺利运行，那就没什么问题，多好啊。
赞一下 回复此楼
4楼2012-09-27 19:51:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kite175

银虫 (初入文坛)
引用回帖:
4楼: Originally posted by jiaoyixiong at 2012-09-27 19:51:35
不知为何纠结？

想知道你的in文件有何问题，最简单的办法就是你运行一下啊，
有什么问题，立刻就知道了。
如果顺利运行，那就没什么问题，多好啊。...

谢谢，谢谢你最美妙的回答，
赞一下 回复此楼
5楼2012-09-28 08:49:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tfl03

木虫 (职业作家)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
zh1987hs: 金币+3, 谢谢 2012-10-01 11:02:11
# 3d Morse crack simulation


units          metal
boundary        s s p
atom_style        atomic


# create geometry

lattice                fcc 4.05  orient x 1 0 0  orient y 0 1 0 orient z 0 0 1
region                box block 0 60 0 40 -0.25 0.25
create_box        5 box
create_atoms        1 box

mass                1 1.0
mass                2 1.0
mass                3 1.0
mass                4 1.0
mass                5 1.0


# Morse potentials

pair_style        morse 2.5
pair_coeff        * * 1.0 1.0 2.5
neighbor        0.3 bin
neigh_modify        every 5 delay 0 check yes

# define groups

region                1 block INF INF INF 2 INF INF
group                lower region 1
region                2 block INF INF 38 INF INF INF
group                upper region 2
group                boundary union lower upper
group                mobile subtract all boundary


region                leftupper block INF 20 20 INF INF INF
region                leftlower block INF 20 INF 20 INF INF

region          bonds   block  INF 20  19.75 20.25 -0.25 0.25
delete_atoms    region  bonds

group                leftupper region leftupper
group                leftlower region leftlower


group                halfupper1 region halfupper1
group                halfupper2 region halfupper2

set                group leftupper type 3
set                group leftlower type 4


set                group lower type 1
set                group upper type 2

neigh_modify        exclude type 3 4
neigh_modify        exclude type 1 2


# 能量最小化
min_style sd
minimize       1.0e-12 1.0e-12  10000 10000

# initial velocities 弛豫

compute                  new mobile temp
velocity        mobile create 0.001 887723 temp new
timestep        0.001
thermo                200
run 10000

velocity        upper  set 0.0 0.01 0.0
velocity        lower  set 0.0 -0.01 0.0
velocity        mobile ramp vy 0.01 0.01 y 20 38 sum yes

# fixes
fix            1 all nve  
fix            2 lower setforce  0.0 NULL 0.0
fix            3 upper setforce  0.0 NULL 0.0
# run
thermo                200
thermo_modify        temp new

dump                1 all atom 1000 dump.crack.lammpstrj
run                300000
赞一下 回复此楼
6楼2012-09-29 11:51:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kite175

银虫 (初入文坛)
引用回帖:
6楼: Originally posted by tfl03 at 2012-09-29 11:51:41
# 3d Morse crack simulation


units          metal
boundary        s s p
atom_style        atomic


# create geometry

lattice                fcc 4.05  orient x 1 0 0  orient y 0 1 0 orient ...

你这个in文件是和我那个一样的。你想告诉我什么呢？
赞一下 回复此楼
7楼2012-10-05 10:37:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chinaboy798

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
!请问你有什么问题？你拿着一个in文件来问有什么问题，难道你不知道有什么问题？
赞一下 回复此楼
8楼2014-11-12 14:07:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lnyuge1989

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
velocity        mobile ramp vy 0.01 0.01 y 20 38 sum yes
这句有错误
应该是
velocity        mobile ramp vy -0.01 0.01 y 20 38 sum yes
赞一下 回复此楼
9楼2014-11-12 15:17:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

闵磊子

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by 1779064387 at 2012-09-27 15:19:20
建议把你的data文件贴一下~

data文件是什么，我只知道in文件
赞一下 回复此楼
10楼2015-04-01 23:46:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Stoudemire29

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
group                halfupper1 region halfupper1
group                halfupper2 region halfupper2
将这两行注释掉就可以了吧
赞一下 回复此楼
11楼2018-10-19 16:21:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Stoudemire29

金虫 (小有名气)
楼主之前是研究裂纹的吗？新手一个，求指导.
赞一下 回复此楼
12楼2018-10-19 16:22:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

云中梦华

禁虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
13楼2018-10-22 12:51:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kite175 的主题更新
131/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291分调剂 +8 上岸小莹加油 2026-04-09 9/450 2026-04-15 00:41 by zzzggc
[考研] 通信工程求调剂！！！ +4 zlb770521 2026-04-14 4/200 2026-04-14 18:19 by lbsjt
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 297，工科调剂?河南农业大学本科 +9 河南农业大学-能 2026-04-14 9/450 2026-04-14 16:03 by Art1977
[考研] 食品与营养（0955）271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 求调剂288 +7 ioodiiij 2026-04-10 9/450 2026-04-13 08:33 by Hayaay
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 一志愿华中农微生物，288分，三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 一志愿郑州大学 22408 305分求调剂 +5 安小满zzz 2026-04-08 5/250 2026-04-12 00:41 by 蓝云思雨
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 求调剂 +6 电气300求调剂不 2026-04-08 6/300 2026-04-11 20:14 by 逆水乘风
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 085402通信工程调剂，有4项学科竞赛国奖（电赛国二），硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 0858求调剂 5+5 Gky09300550， 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 314求调剂 +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
信息提示
关闭
请填处理意见
关闭