 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1057 | 回復(fù): 1 | ||
創(chuàng)意笑笑鐵蟲 (初入文壇)
|
[求助]
使用read_data后，如何建立模型group？
|
|
小弟最近在學(xué)習(xí)lammps軟件進(jìn)行分子模擬，自己建了一個polycrystal模型，想對模型進(jìn)行g(shù)roup，但是運行后報錯，應(yīng)該是region語句的格式問題，但是剛開始學(xué)習(xí)不知道該怎么改，請蟲友們指點，附件上是in文件（in.4txt），data文件（fetwin.lmp）和勢函數(shù)（用的是鋁的勢函數(shù)，隨便找了一個，應(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}" 報錯的內(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 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : Al99.eam.alloy
- 附件 2 : fetwin.lmp
2018-07-03 10:24:26, 762.16 K
2018-07-03 10:24:28, 5.55 M
» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
生物學(xué)學(xué)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
-
求材料，環(huán)境專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
335求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
創(chuàng)意笑笑
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 8.5
- 帖子: 11
- 在線: 10.8小時
- 蟲號: 3127829
- 注冊: 2014-04-11
- 專業(yè): 固體力學(xué)
2樓2018-07-03 10:30:58
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 創(chuàng)意笑笑 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +5 | 某某某某位 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 335分 | 材料與化工專碩 | GPA 4.07 | 有科研經(jīng)歷 +4 | cccchenso 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +7 | 生物工程調(diào)劑 2026-03-17 | 12/600 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
|
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|