24小时热门版块排行榜    

查看: 2560  |  回复: 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的回帖

chinaboy798

银虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
!请问你有什么问题?你拿着一个in文件来问有什么问题,难道你不知道有什么问题?
8楼2014-11-12 14:07:34
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 一志愿东华大学化学070300,求调剂 +6 2117205181 2026-03-21 7/350 2026-03-22 14:59 by allen-yin
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +10 枫桥ZL 2026-03-18 12/600 2026-03-21 22:02 by peike
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 310求调剂 +3 baibai1314 2026-03-16 3/150 2026-03-21 03:56 by JourneyLucky
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 材料工程(专)一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
请填处理意见