| 4 | 1/1 | 返回列表 |
| 查看: 3473 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問題,作者lupei努力將贈送您 5 個金幣 | ||
lupei努力新蟲 (初入文壇)
|
[求助]
lammps使用velocity方法進行石墨烯拉伸
|
|
|
最近使用lammps做石墨烯之字形方向的拉伸模擬,按照網(wǎng)上各種資源設(shè)計了自己的in文件,一直出現(xiàn)拉不斷的情況,邊界原子鍵角已經(jīng)很平了,可是就是不斷裂。 應(yīng)力已經(jīng)到了170gpa,應(yīng)變到了0.3還沒斷。我看很多文獻之字形的理論抗拉強度是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ā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +9 | yang婷 2026-03-02 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 中科大材料299求調(diào)劑 +7 | DAIjiayo 2026-03-05 | 12/600 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)0856材料與化工,前三科206,總分283,求調(diào)劑 +7 | 26考研求調(diào)劑 2026-03-06 | 7/350 |
|
|
[考研] 材料調(diào)劑 +5 | ounce. 2026-03-05 | 13/650 |
|
|
[考研] 273求調(diào)劑 +5 | 星星111222 2026-03-02 | 7/350 |
|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線,總分 281 求調(diào)劑 +6 | 香香 Lu 2026-03-03 | 11/550 |
|
|
[考研] 求調(diào)劑 +5 | danyyyy 2026-03-04 | 5/250 |
|
|
[考研] 274求調(diào)劑 +9 | 一個學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
[考研] 274環(huán)境工程求調(diào)劑 +6 | 扶柳盈江 2026-03-05 | 6/300 |
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +7 | 1一11 2026-03-03 | 7/350 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財 2026-03-04 | 3/150 |
|
|
[考研] 能動297求調(diào)劑,本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
[考研]
材料325求調(diào)劑
30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研]
|
glwshine 2026-03-02 | 5/250 |
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研] 材料調(diào)劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|
|
[考研] 一志愿鄭大材料學(xué)碩298分,求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|