 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1811 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
劍雪封侯新蟲 (小有名氣)
|
[求助]
cannot zero langevin force of 0 atoms
|
|
|
你看這兩個(gè)group里明明有原子�。� 為什么會(huì)出現(xiàn)這種情況？ 未命名.jpg 未命名1.jpg |
回復(fù)此樓» 猜你喜歡
【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額
已經(jīng)有7人回復(fù)
-
334分 一志愿武理-080500 材料求調(diào)劑
已經(jīng)有3人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有9人回復(fù)
-
347求調(diào)劑
已經(jīng)有4人回復(fù)
-
0854人工智能方向招收調(diào)劑
已經(jīng)有4人回復(fù)
-
求b區(qū)院校調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有4人回復(fù)
-
284求調(diào)劑
已經(jīng)有15人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【分享】嘗試lammps中, 分享中...
已經(jīng)有15人回復(fù)
劍雪封侯
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 220.7
- 散金: 3
- 帖子: 75
- 在線: 24.3小時(shí)
- 蟲號(hào): 1717920
- 注冊(cè): 2012-03-26
- 專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
5樓2013-05-13 11:10:32
ddll_slsl
木蟲 (正式寫手)
- 應(yīng)助: 101 (高中生)
- 金幣: 2939.2
- 紅花: 6
- 帖子: 670
- 在線: 255.3小時(shí)
- 蟲號(hào): 2121784
- 注冊(cè): 2012-11-12
- 專業(yè): 動(dòng)力學(xué)與控制
2樓2013-05-13 09:44:29
劍雪封侯
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 220.7
- 散金: 3
- 帖子: 75
- 在線: 24.3小時(shí)
- 蟲號(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
木蟲 (正式寫手)
- 應(yīng)助: 101 (高中生)
- 金幣: 2939.2
- 紅花: 6
- 帖子: 670
- 在線: 255.3小時(shí)
- 蟲號(hào): 2121784
- 注冊(cè): 2012-11-12
- 專業(yè): 動(dòng)力學(xué)與控制
4樓2013-05-13 10:56:46
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|---|---|---|---|---|
|
[考研]
 材料與化工考研調(diào)劑
+9
|
孅華 2026-03-22 | 9/450 |
|
|
|
[考研] 286求調(diào)劑 +11 | Faune 2026-03-21 | 11/550 |
|
|
|
[考研] 303求調(diào)劑 +3 | 藍(lán)山月 2026-03-25 | 3/150 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 食品專碩 一志愿雙一流 328 +3 | xiaom99 2026-03-21 | 4/200 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +4 | 曼111 2026-03-24 | 5/250 |
|
|
|
[考研]
 085404電子信息284分求調(diào)劑
+4
|
13659058978 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研]
求調(diào)劑材料學(xué)碩080500，總分289分
 5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
|
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|