| 4 | 1/1 | 返回列表 |
| 查看: 3479 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問題,作者lupei努力將贈送您 5 個金幣 | |||
lupei努力新蟲 (初入文壇)
|
[求助]
lammps使用velocity方法進(jìn)行石墨烯拉伸
|
||
|
最近使用lammps做石墨烯之字形方向的拉伸模擬,按照網(wǎng)上各種資源設(shè)計了自己的in文件,一直出現(xiàn)拉不斷的情況,邊界原子鍵角已經(jīng)很平了,可是就是不斷裂。 應(yīng)力已經(jīng)到了170gpa,應(yīng)變到了0.3還沒斷。我看很多文獻(xiàn)之字形的理論抗拉強(qiáng)度是120,GPA,130GPA等等。不知道什么地方出現(xiàn)問題。下面是我的in文件 boundary s s s units metal atom_style full read_data perfect.data region left block INF 5 INF INF INF INF region right block 93 INF INF INF INF INF group left region left group right region right group boundary union left right group mobile subtract all boundary # potentials pair_style airebo 1.92 1 1 pair_coeff * * CH.airebo C neighbor 2.0 nsq neigh_modify every 1 delay 0 check yes # define computes compute myTemp all temp compute disp all displace/atom compute pot all pe/atom compute ke all ke/atom compute s1 all stress/atom NULL compute press all pressure thermo_temp compute lx all reduce sum c_s1[1] variable tmp equal "lx" variable L0 equal ${tmp} variable strain equal "(lx-v_L0)/v_L0" thermo_style custom step temp c_myTemp etotal ke pe lz ly lx vol press pxx pyy pzz v_tmp v_L0 v_strain thermo 100 thermo_modify lost ignore timestep 0.001 min_style sd minimize 1.0e-6 1.0e-6 1000 1000 #relax for equilibrium velocity all create 300 4928459 mom yes rot no fix 1 all nve fix 2 all langevin 300.0 300.0 0.08 48279 run 40000 dump 5 all custom 1000 atom.lammpstrj id type x y z vx vy vz xs ys zs c_pot c_ke c_disp[1] c_disp[2] c_disp[3] c_s1[1] c_s1[2] c_s1[3] c_s1[4] c_s1[5] c_s1[6] #set the deformation process velocity boundary set 0.0 0.0 0.0 fix 3 boundary setforce 0.0 0.0 0.0 velocity right set 0.0995 0.0 0.0 sum yes velocity mobile ramp vx 0.0 0.0995 x 5 93 sum yes run 2000000 |
新蟲 (初入文壇)
銅蟲 (小有名氣)
新蟲 (初入文壇)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +9 | 王曉陽- 2026-03-09 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生 +24 | another12 2026-03-04 | 32/1600 |
|
|
[考研] 0817化學(xué)工程319求調(diào)劑 +5 | lv945 2026-03-08 | 6/300 |
|
|
[考研] 313分生物學(xué)求調(diào)劑 +4 | Yyt楊1 2026-03-09 | 5/250 |
|
|
[考研] 290 材料與化工求調(diào)劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
[考研] 0701-322 求調(diào)劑 +3 | jiliuxian 2026-03-06 | 8/400 |
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +7 | T123 tt 2026-03-04 | 7/350 |
|
|
[考研] 誠求調(diào)劑,323分 有專利/科研/實習(xí)經(jīng)歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
[考研] 267化工調(diào)劑求助 +7 | 聰少OZ 2026-03-04 | 7/350 |
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
[考研] 276求調(diào)劑 本科天津工業(yè)大學(xué) 一志愿北京化工大學(xué)0855 +4 | 路lyh123 2026-03-06 | 6/300 |
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
[考研] 材料與化工,291,求調(diào)劑 +7 | 咕嚕咕嚕123123 2026-03-05 | 8/400 |
|
|
[考研] 一志愿西交化工專碩288專業(yè)課93求調(diào)劑 +6 | 好運(yùn)好運(yùn)接接 2026-03-04 | 6/300 |
|
|
[論文投稿]
100+4
|
Stray2021 2026-03-03 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑,總分320分,一志愿中南民族大學(xué) +3 | 19286712440 2026-03-04 | 3/150 |
|
|
[考研] 一志愿中科大080500總分324求調(diào)劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|