| 9 | 1/1 | 返回列表 |
| 查看: 2628 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +4 | 吃吃吃才有意義 2026-03-19 | 4/200 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 6/300 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個(gè)專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|