| 4 | 1/1 | 返回列表 |
| 查看: 3465 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者lupei努力將贈(zèng)送您 5 個(gè)金幣 | ||
lupei努力新蟲(chóng) (初入文壇)
|
[求助]
lammps使用velocity方法進(jìn)行石墨烯拉伸
|
|
|
最近使用lammps做石墨烯之字形方向的拉伸模擬,按照網(wǎng)上各種資源設(shè)計(jì)了自己的in文件,一直出現(xiàn)拉不斷的情況,邊界原子鍵角已經(jīng)很平了,可是就是不斷裂。 應(yīng)力已經(jīng)到了170gpa,應(yīng)變到了0.3還沒(méi)斷。我看很多文獻(xiàn)之字形的理論抗拉強(qiáng)度是120,GPA,130GPA等等。不知道什么地方出現(xiàn)問(wè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 |
新蟲(chóng) (初入文壇)
銅蟲(chóng) (小有名氣)
新蟲(chóng) (初入文壇)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 320材料與化工,求調(diào)劑 +4 | 鶴遨予卿 2026-03-04 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿314求調(diào)劑 +6 | 202111120625 2026-03-03 | 6/300 |
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
[考研] 264求調(diào)劑 +3 | thext 2026-03-03 | 3/150 |
|
|
[考研] 264求調(diào)劑 +7 | 26調(diào)劑 2026-03-03 | 7/350 |
|
|
[考博] 26申博 求博導(dǎo) +3 | 愛(ài)讀書(shū)的小帥 2026-02-28 | 5/250 |
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
[考研] 306求調(diào)劑 +7 | 張張張張oo 2026-03-03 | 7/350 |
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
[考研] 278求調(diào)劑 +3 | 滿天星11_22 2026-03-02 | 3/150 |
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級(jí)已過(guò) 求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
[考研] 考研復(fù)試調(diào)劑,過(guò)國(guó)家線的同學(xué)都可報(bào)名 +7 | 黑!在干嘛 2026-02-28 | 8/400 |
|
|
[考研] 材料復(fù)試調(diào)劑 +5 | 學(xué)材料的點(diǎn) 2026-03-01 | 6/300 |
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 熬夜的貓頭鷹 2026-03-02 | 3/150 |
|
|
[考研] 一志愿鄭大材料學(xué)碩298分,求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|