| 4 | 1/1 | 返回列表 |
| 查看: 3466 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問題,作者lupei努力將贈(zèng)送您 5 個(gè)金幣 | ||
lupei努力新蟲 (初入文壇)
|
[求助]
lammps使用velocity方法進(jìn)行石墨烯拉伸
|
|
|
最近使用lammps做石墨烯之字形方向的拉伸模擬,按照網(wǎng)上各種資源設(shè)計(jì)了自己的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ā)表 | |
|---|---|---|---|---|
|
[考研] 本科太原理工采礦工程,求調(diào)劑 +3 | onlx 2026-03-01 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境調(diào)劑 +6 | 柒槿levana 2026-03-01 | 6/300 |
|
|
[考研] 一志愿東北大學(xué)080500,英一數(shù)二343分求調(diào)劑 +3 | Seele_v5 2026-03-04 | 3/150 |
|
|
[考研] 材料化工調(diào)劑 +15 | 今夏不夏 2026-03-01 | 18/900 |
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
[考研] 環(huán)境工程學(xué)碩288求助調(diào)劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
[考研] 成績(jī)276,專業(yè)代碼0856求調(diào)劑 +7 | 小陳朵 2026-03-03 | 7/350 |
|
|
[考研]
材料270求調(diào)劑
6+6
|
Eiiiio 2026-03-01 | 11/550 |
|
|
[考研]
|
glwshine 2026-03-02 | 5/250 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 調(diào)劑 +5 | 13853210211 2026-03-02 | 7/350 |
|
|
[考研] 課題組接收材料類調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
[考研] 一志愿中石油(華東)本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-02 | 3/150 |
|
|
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5 | 半截的詩(shī)0927 2026-03-02 | 5/250 |
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-02 | 3/150 |
|
|
[基金申請(qǐng)] 成果系統(tǒng)訪問量大,請(qǐng)一小時(shí)后再嘗試。---NSFC啥時(shí)候好哦,已經(jīng)兩天這樣了 +4 | NSFC2026我來(lái)了 2026-02-28 | 4/200 |
|