24小时热门版块排行榜    

查看: 520  |  回复: 1

527751273

银虫 (正式写手)

[求助] 用LJ模拟裂纹扩展,可是结果出来都是从加力的地方断。。。。求解决啊。。。。

我想用LJ模拟裂纹扩展,可是结果出来都是从加力的地方断,而试验是在约束的地方开始断的,实在不懂啊。。

回复此楼

» 猜你喜欢

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

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

527751273

银虫 (正式写手)

in文件如下:
# 2d LJ crack simulation
dimension 2
boundary s s p
atom_style atomic
neighbor 0.3 bin
neigh_modify delay 3
# create geometry
lattice hex 0.93
region box block 0 10 0 40 -5 5
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
# LJ potentials
pair_style  lj/cut 2.5
pair_coeff  * * 1.0 1.0 2.5
# define groups
region 1 block 9.99 10 0 0.01 INF INF
group load1 region 1
region 2 block 9.99 10 23.99 24 INF INF
group load2 region 2
group load union load1 load2
group mobile subtract all load
region 3 block 0 0.01 15.99 16 INF INF
region 4 block 0 0.01 39.99 40 INF INF
group lower1 region 3
group lower2 region 4
group boundary union lower1 lower2
region  5 block INF 2 19.99 20.01 INF INF
region  6 block 8 INF 19.99 20.01 INF INF
group  lowernotch region 5
group  uppernotch region 6
group  notch union lowernotch uppernotch
delete_atoms  group notch
set  group load1 type 2
set  group load2 type 3
set  group lower1 type 4
set  group lower2 type 5
# initial velocities
velocity   load set -2.5e-3 0.0 0.0
minimize   1.0e-5 1.0e6 1000 1000
# fixes
fix  1 all nve
fix  2 load1 setforce 1000 NULL 0.0
fix  3 load2 setforce 5000 NULL 0.0
fix  4 boundary setforce 0.0 0.0 0.0
# run
timestep  0.00003
thermo  200
dump 1  all atom 500 dump8.lammpstrj
run  6000
2楼2011-08-24 21:24:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 527751273 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 081700 调剂 267分 +4 迷人的哈哈 2026-03-23 4/200 2026-03-23 13:02 by ztnimte
[考研] 336化工调剂 +3 王大坦1 2026-03-23 4/200 2026-03-23 11:07 by laoshidan
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 307求调剂 +11 冷笙123 2026-03-17 11/550 2026-03-22 20:16 by edmund7
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +3 jiajunX 2026-03-22 3/150 2026-03-22 19:32 by brblmd
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 17/850 2026-03-22 15:18 by 脱颖而出
[考研] 生物学一志愿985,分数349求调剂 +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
请填处理意见