24小时热门版块排行榜    

查看: 2893  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

云中梦华

禁虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

13楼2018-10-22 12:51:03
已阅   回复此楼   关注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的回帖
★ ★ ★
小木虫: 金币+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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +10 Kittylucky 2026-08-27 10/500 2026-08-28 07:10 by sincosx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
请填处理意见