| 2 | 1/1 | 返回列表 |
| 查看: 1058 | 回復(fù): 1 | ||
創(chuàng)意笑笑鐵蟲 (初入文壇)
|
[求助]
使用read_data后,如何建立模型group?
|
|
小弟最近在學(xué)習(xí)lammps軟件進(jìn)行分子模擬,自己建了一個(gè)polycrystal模型,想對模型進(jìn)行g(shù)roup,但是運(yùn)行后報(bào)錯(cuò),應(yīng)該是region語句的格式問題,但是剛開始學(xué)習(xí)不知道該怎么改,請蟲友們指點(diǎn),附件上是in文件(in.4txt),data文件(fetwin.lmp)和勢函數(shù)(用的是鋁的勢函數(shù),隨便找了一個(gè),應(yīng)該不是勢函數(shù)的問題) 這是in文件 # ------------------------ INITIALIZATION ---------------------------- units metal dimension 3 boundary p p p atom_style atomic # variable latparam equal 4.05 # ----------------------- ATOM DEFINITION ---------------------------- # lattice fcc ${latparam} # region whole block 0 10 0 10 0 10 # create_box 1 whole # lattice fcc ${latparam} orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 # create_atoms 1 region whole read_data fetwin.data region up block -INF INF -INF INF 80 INF region low block -INF INF -INF INF -INF 20 group up region up group low region low # ------------------------ FORCE FIELDS ------------------------------ pair_style eam/alloy pair_coeff * * Al99.eam.alloy Al # ------------------------- SETTINGS --------------------------------- compute csym all centro/atom fcc compute peratom all pe/atom compute stress all stress/atom NULL ###################################### reset_timestep 0 timestep 0.001 velocity all create 300 12345 mom yes rot no fix 1 all npt temp 300 300 1 iso 0 0 1 drag 1 # Set thermo output thermo 10 thermo_style custom step lx ly lz press pxx pyy pzz pe temp run 500 unfix 1 # Store final cell length for strain calculations variable tmp equal "lx" variable L0 equal ${tmp} print "Initial Length, L0: ${L0}" 報(bào)錯(cuò)的內(nèi)容:: PS C:\Users\xxpc-office\Desktop\tensile-twin> lmp_mpi -in in.4txt LAMMPS (11 Aug 2017) OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (../comm.cpp:90) using 1 OpenMP thread(s) per MPI task Reading data file ... orthogonal box = (0 0 0) to (100 100 100) 1 by 1 by 1 MPI processor grid reading atoms ... 84355 atoms ERROR: Expected floating point parameter in input script or data file (../region_block.cpp:38) Last command: region up block -INF INF -INF INF 80 INF data文件開頭內(nèi)容:: 84355 atoms 1 atom types 0.00000000 100.00000000 xlo xhi 0.00000000 100.00000000 ylo yhi 0.00000000 100.00000000 zlo zhi Masses 1 55.84500000 # Fe Atoms # atomic 1 1 25.11560200 83.72774500 76.47610300 2 1 28.27335500 85.32570400 75.71484200 3 1 22.03587500 80.44872200 74.03231600 4 1 21.11270900 82.67894400 74.88441700 5 1 17.95495500 81.08098500 75.64567900 6 1 23.57573800 82.08823300 75.25421000 7 1 24.27046200 84.27690300 74.12315500 8 1 28.96807900 87.51437300 74.58378700 9 1 30.58597000 87.47282000 72.60063100 10 1 27.42821600 85.87486200 73.36189300 11 1 25.81032600 85.91641400 75.34504900 12 1 32.12583300 89.11233200 73.82252500 13 1 32.82055700 91.30100100 72.69147000 |
鐵蟲 (初入文壇)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
孅華 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料292調(diào)劑 +5 | 橘頌思美人 2026-03-23 | 5/250 |
|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +5 | kotoko_ik 2026-03-23 | 6/300 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +4 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 4/200 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 求材料,環(huán)境專業(yè)調(diào)劑 +3 | 18567500178 2026-03-18 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 3/150 |
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|