 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (4942)
- >考研 (795)
- >導(dǎo)師招生 (767)
- >文獻(xiàn)求助 (293)
- >蟲(chóng)友互識(shí) (205)
- >考博 (131)
- >碩博家園 (96)
- >博后之家 (80)
- >論文投稿 (74)
- >基金申請(qǐng) (69)
- >第一性原理 (59)
- >招聘信息布告欄 (58)
- >教師之家 (55)
- >休閑灌水 (43)
- >公派出國(guó) (32)
- >外文書(shū)籍求助 (22)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 1806 | 回復(fù): 9 | ||
劍雪封侯新蟲(chóng) (小有名氣)
|
[求助]
cannot zero langevin force of 0 atoms
|
|
你看這兩個(gè)group里明明有原子��？ 為什么會(huì)出現(xiàn)這種情況？ 未命名.jpg 未命名1.jpg |
回復(fù)此樓» 猜你喜歡
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有7人回復(fù)
-
一志愿C9材料與化工專業(yè)總分300求調(diào)劑
已經(jīng)有4人回復(fù)
-
321求調(diào)劑
已經(jīng)有4人回復(fù)
-
307求調(diào)劑
已經(jīng)有4人回復(fù)
-
293求調(diào)劑
已經(jīng)有6人回復(fù)
-
085602 289分求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北京化工大學(xué)材料與化工 264分各科過(guò)A區(qū)國(guó)家線
已經(jīng)有3人回復(fù)
-
277分求調(diào)劑，跨調(diào)材料
已經(jīng)有3人回復(fù)
-
材料292調(diào)劑
已經(jīng)有6人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【分享】嘗試lammps中, 分享中...
已經(jīng)有15人回復(fù)
ddll_slsl
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 101 (高中生)
- 金幣: 2939.2
- 紅花: 6
- 帖子: 670
- 在線: 255.3小時(shí)
- 蟲(chóng)號(hào): 2121784
- 注冊(cè): 2012-11-12
- 專業(yè): 動(dòng)力學(xué)與控制
2樓2013-05-13 09:44:29
劍雪封侯
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 220.7
- 散金: 3
- 帖子: 75
- 在線: 24.3小時(shí)
- 蟲(chóng)號(hào): 1717920
- 注冊(cè): 2012-03-26
- 專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
|
dimension 3 units real boundary p p f atom_style full neighbor 3.0 bin neigh_modify delay 5 every 1 check yes read_data gufinal2-oh1.data #lattice diamond 3.5667 #region box1 cylinder z 9 16.5 3 3 18 units box #create_atoms 1 region box1 #group tip region box1 # force field pair_style hybrid reax 6.0 1 0 1.0e-6 airebo 2.5 0 0 lj/cut 10.0 pair_coeff * * reax ffield.reax 1 1 1 1 3 2 pair_coeff * * airebo CH.airebo C C C C NULL NULL pair_coeff 1 2*4 lj/cut 0.0654904 3.40 pair_coeff 2 3*4 lj/cut 0.0654904 3.40 pair_coeff 3 4 lj/cut 0.0654904 3.40 pair_coeff 2*4 5 lj/cut 0.101 3.283 pair_coeff 2*4 6 lj/cut 0.0475 3.025 pair_coeff 5 5 lj/cut 0.155402 3.166 pair_coeff 6 6 lj/cut 0.03459 2.65 # define region before define group region rfix block INF INF INF INF INF -8 units box region r1 block INF INF INF INF -8 -4.5 units box region r2 block INF INF INF INF -4.5 -1 units box region r3 block INF INF INF INF -1 INF units box region th1 block INF INF INF 3.5 -8 INF units box region th2 block INF INF 29.5 INF -8 INF units box group gfix region rfix group g1 region r1 group g2 region r2 group g3 region r3 group gth1 region th1 group gth2 region th2 group gth union gth1 gth2 group gthg12 subtract gth g3 group gthg3 subtract gth gthg12 group g12 union g1 g2 group gf12 union gfix g12 # compute and velocity variable fg3x equal fcm(g3,x)#x component of the total force on group g3 variable fg3y equal fcm(g3,y) variable fg3z equal fcm(g3,z) variable vg3x equal vcm(g3,x)#x component of velocity of the center of mass of group g3 variable vg3y equal vcm(g3,y)#y component of velocity of the center of mass of group g3 variable pg3x equal xcm(g3,x)#x component of position of the center of mass of group g3 variable pg3y equal xcm(g3,y)#y component of position of the center of mass of group g3 variable pg3z equal xcm(g3,z)#z component of position of the center of mass of group g3 variable pg2x equal xcm(g2,x)#x component of position of the center of mass of group g2 variable pg2y equal xcm(g2,y)#y component of position of the center of mass of group g2 variable pg2z equal xcm(g2,z)#z component of position of the center of mass of group g2 variable vg2x equal vcm(g2,x)#x component of velocity of the center of mass of group g2 variable vg2y equal vcm(g2,y)#y component of velocity of the center of mass of group g2 #displace_atoms g3 rotate 21.27 25.45 8.08 0.0 0.0 1.0 10.0 units box #displace_atoms g3 move 0.0 0.0 0.5 units box compute temp1 g3 temp compute temp2 g3 temp/com compute temp3 g12 temp compute temp4 gthg12 temp compute temp5 gthg3 temp compute temp6 gthg3 temp/com compute cf1 g3 group/group gf12 compute ke1 g3 ke compute pe1 g3 pe/atom compute peg3 g3 reduce sum c_pe1 compute pe2 g12 pe/atom compute peg12 g12 reduce sum c_pe2 velocity gfix set 0.0 0.0 0.0 units box #velocity tip set 0.0 0.0 0.0 units box velocity g12 create 10 4928459 rot yes dist gaussian units box velocity g3 create 10 4928459 rot yes dist gaussian units box # fixes fix 1 all nve fix 2 gthg12 langevin 10 10 300 699483 zero yes fix 3 gthg3 langevin 10 10 100 699483 zero yes fix_modify 3 temp temp6 #fix fnvt1 g12 nvt temp 10 10 0.1 drag 2.0 #fix fnvt2 g3 nvt temp 10 10 0.1 drag 2.0 #fix_modify fnvt2 temp temp2 #fix 2 g12 temp/berendsen 10 10 0.3 #fix 3 g3 temp/berendsen 10 10 0.3 #fix_modify 3 temp temp2 fix 4 gfix setforce 0.0 0.0 0.0 #fix ftip tip setforce 0.0 0.0 0.0 #fix ftip1 tip rigid single #relax #min_modify dmax 0.01 #minimize 1.0e-6 2.306e-7 1000 100000 #neigh_modify exclude group gfix gfix #neigh_modify exclude group tip tip timestep 1 thermo 100 thermo_style custom step temp pe ke etotal press pzz c_ke1 c_peg3 c_peg12 c_temp1 c_temp3 c_temp4 c_temp6 v_fg3x v_fg3y v_fg3z & v_pg3x v_pg3y v_pg3z v_vg3x v_vg3y c_cf1 c_cf1[1] c_cf1[2] c_cf1[3] v_pg2x v_pg2y v_pg2z v_vg2x v_vg2y thermo_modify temp temp2 dump 1 all xyz 200 gufinal2-oh.xyz dump 2 all atom 200 gufinal2-oh.lammpstrj run 50000 #Run #variable dz1 equal ramp(17,15.6) #fix fw1 all wall/harmonic zhi v_dz1 0.1 1.0 2.5 units box #fix_modify fw1 energy yes #run 30000 #fix fw1 all wall/harmonic zhi 15.6 0.1 1.0 2.5 units box #fix_modify fw1 energy yes #min_modify dmax 0.01 #minimize 1.0e-6 1.0e-8 1000 100000 #run 40000 velocity g3 set 0.001 NULL NULL sum yes units box run 100000 |
3樓2013-05-13 10:45:13
ddll_slsl
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 101 (高中生)
- 金幣: 2939.2
- 紅花: 6
- 帖子: 670
- 在線: 255.3小時(shí)
- 蟲(chóng)號(hào): 2121784
- 注冊(cè): 2012-11-12
- 專業(yè): 動(dòng)力學(xué)與控制
4樓2013-05-13 10:56:46
劍雪封侯
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 220.7
- 散金: 3
- 帖子: 75
- 在線: 24.3小時(shí)
- 蟲(chóng)號(hào): 1717920
- 注冊(cè): 2012-03-26
- 專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
5樓2013-05-13 11:10:32
ddll_slsl
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 101 (高中生)
- 金幣: 2939.2
- 紅花: 6
- 帖子: 670
- 在線: 255.3小時(shí)
- 蟲(chóng)號(hào): 2121784
- 注冊(cè): 2012-11-12
- 專業(yè): 動(dòng)力學(xué)與控制
6樓2013-05-13 11:19:38
劍雪封侯
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 220.7
- 散金: 3
- 帖子: 75
- 在線: 24.3小時(shí)
- 蟲(chóng)號(hào): 1717920
- 注冊(cè): 2012-03-26
- 專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
7樓2013-05-13 16:57:54
ddll_slsl
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 101 (高中生)
- 金幣: 2939.2
- 紅花: 6
- 帖子: 670
- 在線: 255.3小時(shí)
- 蟲(chóng)號(hào): 2121784
- 注冊(cè): 2012-11-12
- 專業(yè): 動(dòng)力學(xué)與控制
8樓2013-05-13 17:02:33
劍雪封侯
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 220.7
- 散金: 3
- 帖子: 75
- 在線: 24.3小時(shí)
- 蟲(chóng)號(hào): 1717920
- 注冊(cè): 2012-03-26
- 專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
9樓2013-05-13 18:58:56
gswylq
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 10 (幼兒園)
- 金幣: 2546.6
- 散金: 120
- 紅花: 1
- 帖子: 782
- 在線: 730.1小時(shí)
- 蟲(chóng)號(hào): 659377
- 注冊(cè): 2008-11-21
- 專業(yè): 金屬材料的凝固與結(jié)晶學(xué)
10樓2013-10-01 22:15:41
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 劍雪封侯 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +6 | 加一一九 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研]
 085404電子信息284分求調(diào)劑
+4
|
13659058978 2026-03-24 | 4/200 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 7/350 |
|
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|