| 9 | 1/1 | 返回列表 |
| 查看: 2627 | 回復(fù): 8 | ||
zhoudeyuan銀蟲 (小有名氣)
|
[求助]
納米壓痕過程中程序出錯(cuò),哪位大神幫我解答下!
|
|
我的納米壓痕程序如下,運(yùn)行結(jié)束后,牛頓原子層很不穩(wěn)定,且log文件提示錯(cuò)誤,不知道哪位大神能夠幫我解答下,小弟萬分感謝! # 3d contact simulation(Sphere_Slid) units metal dimension 3 boundary p p p atom_style atomic neighbor 3.0 bin neigh_modify delay 5 newton on # create geometry region box block 0.0 40.0 0.0 40.0 0.0 50.0 units box create_box 4 box mass 1 29.0 mass 2 29.0 mass 3 29.0 mass 4 12.0 # atom regions region reg block 0.0 40.0 0.0 40.0 0.0 25.0 units box region reg1 block 3.0 37.0 3.0 37.0 3.0 25.0 units box region newtonian_layer block 8.0 32.0 8.0 32.0 8.0 25.0 units box region asperity sphere 20 20 30 5 units box lattice diamond 5.431 create_atoms 1 region reg units box lattice diamond 3.5667 create_atoms 2 region asperity units box # Tersoff Potentials pair_style tersoff pair_coeff * * SiC.tersoff Si Si Si C # define groups group reg region reg group reg1 region reg1 group newtonian_layer region newtonian_layer group thermostat_layer subtract reg1 newtonian_layer group boundary_layer_down subtract reg reg1 group asperity region asperity group integrater union reg asperity set group newtonian_layer type 1 set group thermostat_layer type 2 set group boundary_layer_down type 3 set group asperity type 4 # initial velocities compute newTemp thermostat_layer temp velocity thermostat_layer create 293.0 5812775 temp newTemp units box # fixed fix 1 integrater nve fix 2 boundary_layer_down setforce 0.0 0.0 0.0 velocity newtonian_layer set NULL NULL NULL units box velocity thermostat_layer create 293.0 5812775 temp newTemp units box fix 3 thermostat_layer temp/rescale 10 293.0 293.0 10 1 fix_modify 3 temp newTemp fix 4 asperity setforce 0.0 0.0 0.0 compute forcez newtonian_layer property/atom fz compute force newtonian_layer reduce sum c_forcez compute forcex newtonian_layer property/atom fx compute lateralforce newtonian_layer reduce sum c_forcex timestep 0.001 thermo 5 thermo_style custom step temp press pe ke etotal c_force c_lateralforce vol thermo_modify temp newTemp dump 1 all atom 50 dump.contact_block_slid fix 6 asperity move linear 0.0 0.0 -1.0 units box run 2000 fix 6 asperity move linear 0.0 0.0 1.0 units box run 2000 log文件出錯(cuò): WARNING: Temperature for thermo pressure is not for group all (thermo.cpp:432) WARNING: One or more atoms are time integrated more than once (modify.cpp:216), 我不知道問題出在哪兒了,希望高手給予指點(diǎn),謝謝! |

木蟲 (正式寫手)
銀蟲 (小有名氣)

木蟲 (正式寫手)
銅蟲 (小有名氣)
銀蟲 (小有名氣)

銀蟲 (小有名氣)

銀蟲 (小有名氣)
木蟲 (初入文壇)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招收調(diào)劑碩士 +3 | lidianxing 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 274求調(diào)劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|