 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 2200 | 回復: 2 | ||
奈落之夜木蟲 (正式寫手)
|
[求助]
關于lammps 質(zhì)心速度vcm（）如何輸出問題
|
| 新手剛學，要輸出質(zhì)心速度 想用vcm這個函數(shù) 根據(jù)說明書 我先定義了變量 variable v_m equal vcm(motion,v_m) motion是一個group 然后輸出 dump v_data motion cfg 10 dump.snap.*.cfg $v_m 可能很多語法錯誤 manual里面vcm（group，dim）前面應該是定義的組，后面的自變量不知道是什么定了，也不知道該怎么輸出 |
回復此樓» 猜你喜歡
北科281學碩材料求調(diào)劑
已經(jīng)有13人回復
-
289求調(diào)劑
已經(jīng)有11人回復
-
網(wǎng)絡空間安全0839招調(diào)劑
已經(jīng)有5人回復
-
材料277求調(diào)劑
已經(jīng)有4人回復
-
274求調(diào)劑
已經(jīng)有7人回復
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有3人回復
-
302求調(diào)劑
已經(jīng)有3人回復
-
312求調(diào)劑
已經(jīng)有5人回復
-
329求調(diào)劑
已經(jīng)有3人回復
-
【2026考研調(diào)劑】制藥工程 284分 求相關專業(yè)調(diào)劑名額
已經(jīng)有8人回復
bluewhale
鐵桿木蟲 (正式寫手)
- 應助: 58 (初中生)
- 金幣: 8438.3
- 散金: 50
- 紅花: 10
- 帖子: 890
- 在線: 362.2小時
- 蟲號: 278289
- 注冊: 2006-09-10
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
【答案】應助回帖
★ ★ ★ ★
感謝參與，應助指數(shù) +1
月只藍: 金幣+2, 感謝指導！ 2013-11-01 08:35:56
奈落之夜: 金幣+2, ★有幫助 2013-11-01 09:10:33
感謝參與，應助指數(shù) +1
月只藍: 金幣+2, 感謝指導！ 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
奈落之夜
木蟲 (正式寫手)
- 應助: 6 (幼兒園)
- 金幣: 1890.6
- 紅花: 1
- 帖子: 393
- 在線: 73.9小時
- 蟲號: 2550025
- 注冊: 2013-07-17
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
3樓2013-11-01 09:11:14
相關版塊跳轉(zhuǎn)
我要訂閱樓主 奈落之夜 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +17 | 陌の森林 2026-03-18 | 17/850 |
|
|
|
[考研] 287求調(diào)劑 +10 | 晨昏線與星海 2026-03-19 | 11/550 |
|
|
[考研] 0703化學求調(diào)劑 +6 | 奶油草莓. 2026-03-22 | 7/350 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分 +6 | @taotao 2026-03-19 | 6/300 |
|
|
|
[考研] 305分求調(diào)劑（食品工程） +5 | Sxy112 2026-03-21 | 7/350 |
|
|
|
[考研] 一志愿華東理工大學081700，初試分數(shù)271 +5 | kotoko_ik 2026-03-23 | 6/300 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實習經(jīng)歷 +9 | 材料學257求調(diào)劑 2026-03-23 | 10/500 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|