 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2202 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
奈落之夜木蟲 (正式寫手)
|
[求助]
關(guān)于lammps 質(zhì)心速度vcm（）如何輸出問題
|
|
| 新手剛學(xué)，要輸出質(zhì)心速度 想用vcm這個(gè)函數(shù) 根據(jù)說明書 我先定義了變量 variable v_m equal vcm(motion,v_m) motion是一個(gè)group 然后輸出 dump v_data motion cfg 10 dump.snap.*.cfg $v_m 可能很多語法錯(cuò)誤 manual里面vcm（group，dim）前面應(yīng)該是定義的組，后面的自變量不知道是什么定了，也不知道該怎么輸出 |
回復(fù)此樓» 猜你喜歡
307求調(diào)劑
已經(jīng)有13人回復(fù)
-
071010 323 分求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料調(diào)劑
已經(jīng)有9人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有8人回復(fù)
-
284求調(diào)劑
已經(jīng)有15人回復(fù)
-
085602化工求調(diào)劑（331分）
已經(jīng)有3人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有5人回復(fù)
奈落之夜
木蟲 (正式寫手)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1890.6
- 紅花: 1
- 帖子: 393
- 在線: 73.9小時(shí)
- 蟲號(hào): 2550025
- 注冊: 2013-07-17
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
3樓2013-11-01 09:11:14
bluewhale
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 58 (初中生)
- 金幣: 8438.3
- 散金: 50
- 紅花: 10
- 帖子: 890
- 在線: 362.2小時(shí)
- 蟲號(hào): 278289
- 注冊: 2006-09-10
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
【答案】應(yīng)助回帖
★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣+2, 感謝指導(dǎo)！ 2013-11-01 08:35:56
奈落之夜: 金幣+2, ★有幫助 2013-11-01 09:10:33
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣+2, 感謝指導(dǎo)！ 2013-11-01 08:35:56
奈落之夜: 金幣+2, ★有幫助 2013-11-01 09:10:33
|
參照下面的例子： heated Brownian particle # typical input script used to investigate the Brownian motion of a heated nanoparticle. The diffusion coefficient is computed as the plateau value of the integral of the velocity autocorrelation function. Details in: # L. Joly, S. Merabia and J.-L. Barrat, EPL 94, 50007 (2011) # http://iopscience.iop.org/0295-5075/94/5/50007/ # the data file is generated by the perl script 'data.pl'. # fix ave/sph is similar to fix ave/spatial but using spherical bins (enclosed sources might be broken with latest lammps versions). # to apply a rescaling thermostat only to liquid atoms lying beyond 15 sigma from the center of the particle (see article for details), fix recenter is used to constrain the center-of-mass of the particle to remain at the origin, with keyword 'shift all' so that it does not alter the dynamics of the system. atom_style bond read_data data.brown group liq type 1 group sol type 2 region far sphere 0.0 0.0 0.0 15.9 side out units box region int1 sphere 0.0 0.0 0.0 14.4 side out units box region int2 sphere 0.0 0.0 0.0 15.9 units box region int intersect 2 int1 int2 units box region ins sphere 0.0 0.0 0.0 14.4 units box compute Tsol sol temp/com compute Tfar liq temp/region far compute Tint liq temp/region int group far region far group int region int group ins region ins group insliq intersect ins liq neigh_modify delay 0 every 1 check yes pair_style lj/cut 2.5 pair_modify shift yes bond_style fene special_bonds fene mass * 1.0 pair_coeff * * 1.0 1.0 bond_coeff * 30.0 1.5 1.0 1.0 variable ts equal step variable vx equal vcm(sol,x)-vcm(liq,x) variable vy equal vcm(sol,y)-vcm(liq,y) variable vz equal vcm(sol,z)-vcm(liq,z) fix 3 all print 5 "${ts} ${vx} ${vy} ${vz}" file print.brown screen no title "# ts vx vy vz" fix 4vx all ave/correlate 5 2401 1200000 v_vx type full ave running start 200000 file cvx.brown fix 4vy all ave/correlate 5 2401 1200000 v_vy type full ave running start 200000 file cvy.brown fix 4vz all ave/correlate 5 2401 1200000 v_vz type full ave running start 200000 file cvz.brown thermo_style custom step cpu c_Tsol c_Tint c_Tfar press vol thermo 1000 compute KEatom all ke/atom variable Tatom atom (2.0/3.0)*c_KEatom fix 5 all ave/sph 100 2000 200000 x origin 0.1 density/number v_Tatom units box file sph.brown fix 1 liq nph aniso 0.015 0.015 5.0 dilate partial fix 1far liq temp/rescale 1 0.75 0.75 0.0 1.0 fix_modify 1far temp Tfar fix 1int liq temp/rescale 1 0.75 0.75 0.0 1.0 fix_modify 1int temp Tint fix 1sol sol nvt temp 3.5 3.5 0.5 fix_modify 1sol temp Tsol variable vol equal vol fix avevol all ave/time 10 8000 160000 v_vol variable l equal 0.5*((f_avevol)^(1./3.)) fix 2 sol recenter 0.0 0.0 0.0 shift all units box velocity sol create 0.75 36582 dist gaussian mom yes rot yes velocity int create 0.75 36582 dist gaussian mom yes rot yes velocity far create 0.75 36582 dist gaussian mom yes rot yes velocity insliq create 0.75 36582 dist gaussian mom yes rot yes run 160000 unfix 1 fix 1 liq nve fix dxyz liq deform 1 x final -$l $l y final -$l $l z final -$l $l units box unfix 2 fix 2 sol recenter 0.0 0.0 0.0 shift all units box run 10000 unfix dxyz run 30000 #variable F equal 10.0 #variable Fs equal $F/count(sol) #fix 9s sol addforce ${Fs} 0.0 0.0 #variable Fl equal -$F/count(liq) #fix 9l liq addforce ${Fl} 0.0 0.0 #unfix 2 #fix 2 sol recenter 0.0 0.0 0.0 shift all units box run 12000000 |
2樓2013-10-31 20:45:20
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 283求調(diào)劑（080500） +14 | A child 2026-03-27 | 14/700 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10 | 景隨風(fēng) 2026-03-29 | 15/750 |
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +6 | minht 2026-03-29 | 6/300 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研]
 304材料求調(diào)劑
+4
|
鐘llll 2026-03-26 | 4/200 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|