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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 有人做裂纹模拟的吗？看看我这个裂纹的in文件有什么问题
51/1返回列表
查看: 2542  |  回复: 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

银虫 (初入文坛)
引用回帖:
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的回帖
查看全部 13 个回答

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的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854复试调剂 276 +4 wmm9 2026-03-01 6/300 2026-03-02 09:28 by 热情沙漠
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[考研] 材料复试调剂 +4 学材料的点 2026-03-01 5/250 2026-03-02 08:26 by houyaoxu
[考研] 哈工大计算机刘劼团队招生 +3 hit_aiot 2026-03-01 5/250 2026-03-02 07:48 by 得鹿梦鱼111
[考研] 欢迎采矿、地质、岩土、计算机、人工智能等专业的同学报考 +4 pin8023 2026-02-28 6/300 2026-03-02 06:35 by 汪！？！
[考研] 0857调剂 +4 一ll半 2026-02-28 5/250 2026-03-02 02:33 by 908055542
[考研] 279求调剂 +3 dua1 2026-03-01 4/200 2026-03-02 00:23 by 大脸蛋子
[基金申请] 本子写完了，给DS兄弟看了，得了92分 +3 Doma 2026-03-01 7/350 2026-03-02 00:00 by jnzsy
[考研] 材料化工调剂 +12 今夏不夏 2026-03-01 13/650 2026-03-01 23:32 by L135790
[考研] 0856调剂 +5 刘梦微 2026-02-28 5/250 2026-03-01 22:30 by wang_dand
[考研] 274求调剂 +3 cgyzqwn 2026-03-01 6/300 2026-03-01 21:24 by cgyzqwn
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 7/350 2026-03-01 17:20 by 2026计算机_诚心
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 303求调剂 +4 今夏不夏 2026-03-01 4/200 2026-03-01 14:46 by 嘟嘟小浣熊
[考研] 求调剂 +6 repeatt?t 2026-02-28 6/300 2026-03-01 14:37 by Sakura绘
[考研] 材料284求调剂，一志愿郑州大学英一数二专硕 +10 想上岸的土拨鼠 2026-02-28 10/500 2026-03-01 14:12 by yc258
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭