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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 有人做裂纹模拟的吗？看看我这个裂纹的in文件有什么问题
131/1返回列表
查看: 2556  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 318求调剂 +3 plum李子 2026-03-21 6/300 2026-03-22 11:35 by plum李子
[考研] 305分求调剂（食品工程） +3 Sxy112 2026-03-21 5/250 2026-03-22 11:35 by 杨杨杨紫
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 16/800 2026-03-22 11:20 by 脱颖而出
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 有没有道铁/土木的想调剂南林，给自己招师弟中～ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
信息提示
关闭
请填处理意见
关闭