 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1582 | 回復(fù): 3 | |||
charliecaili鐵蟲 (初入文壇)
|
[求助]
石墨烯拉伸模擬 已有1人參與
|
|
本人做石墨烯拉伸模擬，想固定石墨烯板Y軸兩端各2個(gè)單位長(zhǎng)度，但是模擬結(jié)果出來卻沒有對(duì)兩端固定，請(qǐng)大神幫我看看問題出在哪里，跪謝��！ ## INITIALIZATION ###################### units metal dimension 3 boundary f p f atom_style atomic newton on ## ATOM DEFINITION ###################### lattice custom 2.4595 a1 1 0 0 a2 0 1.73203 0 a3 0 0 1 basis 0 0.666666666 0 basis 0.5 0.5 0 basis 0.5 0.16666666 0 & basis 0 0 0 region box block -20 20 -10 10 -50 50 create_box 1 box region mid1 block -8 8 -8 8 -0.5 0.5 create_atoms 1 region mid1 mass 1 12.010000 ## SETTINGS ###################### pair_style airebo 3.0 pair_coeff * * potentials/CH.airebo C region lower block -8 8 -8 -6 -0.5 0.5 region upper block -8 8 6 8 -0.5 0.5 group upper region upper group lower region lower group boundary union upper lower group mobile subtract all boundary # initial velocities fix 1 all nvt temp 300 300 0.1 velocity boundary set 0.0 0.0 0.0 fix 2 boundary setforce 0.0 0.0 0.0 min_style cg minimize 1e-25 1e-25 5000 10000 unfix 1 unfix 2 timestep 0.001 fix 1 all npt temp 300 300 0.05 y 0 0 0.5 thermo 20 compute myTemp all temp compute 1 all stress/atom myTemp compute 2 all pe/atom pair bond compute 3 all reduce sum c_1[1] c_1[2] c_1[3] shell mkdir 10.9duxg3 shell cd 10.9duxg3 thermo_style custom step temp pe ke etotal lx ly press pxx pyy c_3[1] c_3[2] c_3[3] dump 1 all custom 1000 dump.chiyu.* id type x y z ## RUNING-RELAXATION ###################### run 30000 ## RUNNING-DEFORMATION ###################### unfix 1 undump 1 fix 1 all nvt temp 300 300 0.05 fix 2 all ave/atom 1 1000 1000 c_1[1] c_1[2] c_1[3] c_2 fx fy fz dump 1 all custom 1000 dump.new.* id type x y z vx vy vz c_1[1] c_1[2] c_1[3] & c_2 f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7] variable srate equal 1.0e10 variable srate1 equal "v_srate / 1.0e12" fix 3 all deform 1 y erate ${srate1} units box remap x ################### variable tmp equal "ly" variable ly0 equal ${tmp} print "Initial Length, ly0: ${ly0}" variable strain equal "(ly- v_ly0)/v_ly0" # Output strain and stress info to file # for units metal, pressure is in [bars] = 100 [kPa] = 1/10000 [GPa] # p2 are in GPa variable p1 equal "v_strain" variable p2 equal "-pyy/10000" variable p3 equal "pe" fix equill all print 500 "${p1} ${p2} ${p3}" file data.1.txt screen no run 40000 ########### ###################################### # SIMULATION DONE print "All done" 結(jié)果圖.png |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : CH.airebo
2015-09-04 14:47:54, 941.56 K
» 收錄本帖的淘帖專輯推薦
 分子動(dòng)力學(xué)模擬 |
» 猜你喜歡
0703化學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有10人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- lammps中的錯(cuò)誤
已經(jīng)有1人回復(fù)
- 科學(xué)家模擬合成新型石墨烯（轉(zhuǎn)載）
已經(jīng)有0人回復(fù)
- 最薄的壓電材料——氮化石墨烯
已經(jīng)有4人回復(fù)
- 石墨烯z方向褶皺
已經(jīng)有1人回復(fù)
- 石墨烯拉伸的動(dòng)力學(xué)模擬，用MS建好模型之后接著該干嘛。
已經(jīng)有0人回復(fù)
- 兩層石墨pull mdp編寫
已經(jīng)有7人回復(fù)
- 請(qǐng)問做原子模擬或者分子模擬時(shí)勢(shì)函數(shù)的選取
已經(jīng)有6人回復(fù)
- 聚合物拉伸
已經(jīng)有0人回復(fù)
- 石墨烯單軸拉伸模擬
已經(jīng)有2人回復(fù)
- 石墨烯拉伸力學(xué)性質(zhì)模擬
已經(jīng)有1人回復(fù)
panbin
木蟲 (正式寫手)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 4751.9
- 散金: 120
- 紅花: 4
- 帖子: 517
- 在線: 173.4小時(shí)
- 蟲號(hào): 1806466
- 注冊(cè): 2012-05-09
- 性別: GG
- 專業(yè): 金屬基復(fù)合材料
2樓2015-09-04 16:06:34
charliecaili
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 231
- 帖子: 12
- 在線: 11.7小時(shí)
- 蟲號(hào): 2920736
- 注冊(cè): 2014-01-08
- 性別: GG
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
3樓2015-09-05 08:34:20
panbin
木蟲 (正式寫手)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 4751.9
- 散金: 120
- 紅花: 4
- 帖子: 517
- 在線: 173.4小時(shí)
- 蟲號(hào): 1806466
- 注冊(cè): 2012-05-09
- 性別: GG
- 專業(yè): 金屬基復(fù)合材料
4樓2015-09-05 20:11:42
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 charliecaili 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +6 | 煙汐憶海 2026-03-16 | 16/800 |
|
|---|---|---|---|---|
|
|
[考研] 求材料調(diào)劑 +8 | 隔壁陳先生 2026-03-12 | 8/400 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 材料專碩英一數(shù)二306 +4 | z1z2z3879 2026-03-18 | 4/200 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +6 | d如愿上岸 2026-03-12 | 9/450 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +4 | 芬達(dá)46 2026-03-12 | 4/200 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|