 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 分子模擬 » MS » 石墨烯拉伸模擬
| 5 | 1/1 | 返回列表 |
| 查看: 1581 | 回復: 3 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
charliecaili鐵蟲 (初入文壇)
|
[求助]
石墨烯拉伸模擬 已有1人參與
|
||
|
本人做石墨烯拉伸模擬，想固定石墨烯板Y軸兩端各2個單位長度，但是模擬結果出來卻沒有對兩端固定，請大神幫我看看問題出在哪里，跪謝��！ ## 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" 結果圖.png |
回復此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產權問題，其責任在于用戶本人，如對版權有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : CH.airebo
2015-09-04 14:47:54, 941.56 K
» 收錄本帖的淘帖專輯推薦
 分子動力學模擬 |
» 猜你喜歡
281求調劑（0805）
已經有16人回復
-
304求調劑
已經有6人回復
-
材料工程專碩調劑
已經有6人回復
-
一志愿天大材料與化工（085600）總分338
已經有4人回復
-
085700資源與環(huán)境308求調劑
已經有3人回復
-
求材料調劑
已經有8人回復
-
294求調劑材料與化工專碩
已經有5人回復
-
一志愿華中科技大學，080502，354分求調劑
已經有4人回復
-
一志愿吉林大學材料學碩321求調劑
已經有6人回復
-
085410人工智能專碩317求調劑（0854都可以）
已經有3人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- lammps中的錯誤
已經有1人回復
- 科學家模擬合成新型石墨烯（轉載）
已經有0人回復
- 最薄的壓電材料——氮化石墨烯
已經有4人回復
- 石墨烯z方向褶皺
已經有1人回復
- 石墨烯拉伸的動力學模擬，用MS建好模型之后接著該干嘛。
已經有0人回復
- 兩層石墨pull mdp編寫
已經有7人回復
- 請問做原子模擬或者分子模擬時勢函數的選取
已經有6人回復
- 聚合物拉伸
已經有0人回復
- 石墨烯單軸拉伸模擬
已經有2人回復
- 石墨烯拉伸力學性質模擬
已經有1人回復
charliecaili
鐵蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 231
- 帖子: 12
- 在線: 11.7小時
- 蟲號: 2920736
- 注冊: 2014-01-08
- 性別: GG
- 專業(yè): 金屬材料的微觀結構
3樓2015-09-05 08:34:20
panbin
木蟲 (正式寫手)
- 應助: 6 (幼兒園)
- 金幣: 4751.9
- 散金: 120
- 紅花: 4
- 帖子: 517
- 在線: 173.4小時
- 蟲號: 1806466
- 注冊: 2012-05-09
- 性別: GG
- 專業(yè): 金屬基復合材料
2樓2015-09-04 16:06:34
panbin
木蟲 (正式寫手)
- 應助: 6 (幼兒園)
- 金幣: 4751.9
- 散金: 120
- 紅花: 4
- 帖子: 517
- 在線: 173.4小時
- 蟲號: 1806466
- 注冊: 2012-05-09
- 性別: GG
- 專業(yè): 金屬基復合材料
4樓2015-09-05 20:11:42
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專碩330分求調劑 +5 | 小材化本科 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 311求調劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 0703化學調劑，求各位老師收留 +9 | 秋有木北 2026-03-14 | 9/450 |
|
|
|
[考研] 環(huán)境工程調劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 308求調劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 274求調劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 一志愿211 0703方向310分求調劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 085600調劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研] 285求調劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 288求調劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[基金申請] 現在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 復試調劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 材料080500調劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] 求調劑 +5 | 一定有學上- 2026-03-12 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|
|
|
[論文投稿]
投稿問題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 070303一志愿西北大學學碩310找調劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|