lammps該怎么入門？我自己都快忘了我研究生畢設(shè)用的就是lammps。如今都工作了小十年了，在研究所從事結(jié)構(gòu)動(dòng)力學(xué)試驗(yàn)工作，lammps忘的差不多了，翻看自己的論文還有之前的幾十G的文獻(xiàn)資料和計(jì)算數(shù)

lammps該怎么入門？

我自己都快忘了我研究生畢設(shè)用的就是lammps。如今都工作了小十年了，在研究所從事結(jié)構(gòu)動(dòng)力學(xué)試驗(yàn)工作，lammps忘的差不多了，翻看自己的論文還有之前的幾十G的文獻(xiàn)資料和計(jì)算數(shù)據(jù)，都不敢相信是自己做的。

這里我只能提供一個(gè)思路，細(xì)節(jié)上的東西實(shí)在是忘了。用lammps首先得熟悉分子動(dòng)力學(xué)，我們是兩年半研究生學(xué)制，大四研一都在搞有限元計(jì)算，研二我用了差不多一年時(shí)間讀分子動(dòng)力學(xué)相關(guān)文獻(xiàn)，啥都沒干，研三開始計(jì)算和寫論文，半年搞定。個(gè)人感覺碩士課題雖淺點(diǎn)，但我的基礎(chǔ)還是比較牢固。

用lammps還得過linux這一關(guān)，大學(xué)我喜歡折騰linux發(fā)行版，筆記本格式化了N次，實(shí)際上省略了學(xué)習(xí)linux的過程，直接上手編譯lammps。

程序就選從官網(wǎng)下載最新版，勢(shì)函數(shù)記不清從哪找的，基本根據(jù)文獻(xiàn)中來(lái)，也嘗試過自己寫勢(shì)函數(shù)，結(jié)合atom以及VMD可視化，用matlab做分析，origin畫圖，基本就夠了。

學(xué)習(xí)路徑基本是從易到難，先做最簡(jiǎn)單的算例，然后需要復(fù)現(xiàn)文獻(xiàn)中的結(jié)果，根據(jù)自己的課題挑選合適的勢(shì)函數(shù)，需要做大量驗(yàn)證，包括計(jì)算材料參數(shù)，計(jì)算結(jié)晶過程等等。當(dāng)然也得根據(jù)計(jì)算機(jī)的算力確定計(jì)算規(guī)模。

遇到不懂的就查文獻(xiàn)，查文獻(xiàn)，查文獻(xiàn)，我當(dāng)初也在小木蟲論壇混過，里面大神很多。總之多嘗試多做算例，肯定能上手。

我現(xiàn)在其實(shí)基本是門外漢，說(shuō)的不對(duì)的地方多多包涵。目前還可以探討一下結(jié)構(gòu)動(dòng)力學(xué)試驗(yàn)的相關(guān)問題。

lammps中怎么固定一族原子？

可以用 dump 命令， 里面加上 fx fy fz 這三項(xiàng)就是每個(gè)原子受到的力在三個(gè)方向上的分量。 例如：dump 2 subgroup custom 1000 dump.atom id fx fy fz LAMMPS會(huì)每隔1000個(gè)時(shí)間步，把叫做“subgroup”這個(gè)組里面的所有原子的序號(hào)和三個(gè)方向上所受