 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1899 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
s455761311鐵蟲 (初入文壇)
|
[交流]
哪位大佬能幫著看看這個(gè) in文件有哪些問題 已有7人參與
|
||
|
units metal # 單位為lammps 中的metel類型 boundary p p p # 邊界條件 atom_style atomic # 原子模式 #special_bonds fene region workpiece1 block 0 47.385 0 50 0 45 units box region workpiece2 block 47.385 47.91 0 50 0 45 units box region workpiece3 block 47.91 150 0 50 0 45 units box region workpiece union 3 workpiece1 workpiece2 workpiece3 units box region grain cone y 177 22.5 2.5 34.50067 40 60 units box region total union 4 workpiece1 workpiece2 workpiece3 grain units box #建鐵素體 create_box 4 total mass 1 55.84 #Fe原子 mass 2 12.011 # C 原子 mass 3 10.811 #B原子 mass 4 14.007 #N原子 lattice custom 2.86 a1 1.0 0.0 0.0 & a2 0.0 1.0 0.0 & a3 0.0 0.0 1.0 & basis 0.0 0.0 0.0 & basis 0.5 0.5 0.5 & basis 0.0 0.0 0.5 & basis 0.0 0.5 0.5 & basis 0.0 0.5 0.0 & basis 0.5 0.0 0.0 & basis 0.5 0.5 0.0 & basis 0.5 0.0 0.5 create_atoms 1 region workpiece1 basis 1 1 basis 2 1 basis 3 2 basis 4 2 basis 5 2 basis 6 2 basis 7 2 basis 8 2 units box #建奧氏體 lattice custom 3.58 a1 1.0 0.0 0.0 & a2 0.0 1.0 0.0 & a3 0.0 0.0 1.0 & basis 0.0 0.0 0.0 & basis 0.5 0.5 0.5 & basis 0.0 0.0 0.5 & basis 0.0 0.5 0.5 & basis 0.0 0.5 0.0 & basis 0.5 0.0 0.0 & basis 0.5 0.5 0.0 & basis 0.5 0.0 0.5 create_atoms 1 region workpiece2 basis 1 1 basis 2 2 basis 3 2 basis 4 1 basis 5 2 basis 6 2 basis 7 1 basis 8 1 units box #建馬氏體 lattice custom 3.018 a1 0.943 0.0 0.0 & a2 0.0 0.943 0.0 & a3 0.0 0.0 1 & basis 0.0 0.0 0.0 & basis 0.5 0.5 0.5 & basis 0.0 0.0 0.5 & basis 0.0 0.5 0.5 & basis 0.0 0.5 0.0 & basis 0.5 0.0 0.0 & basis 0.5 0.5 0.0 & basis 0.5 0.0 0.5 create_atoms 1 region workpiece3 basis 1 1 basis 2 1 basis 3 2 basis 4 2 basis 5 2 basis 6 2 basis 7 2 basis 8 2 lattice diamond 3.615 #CBN砂輪 create_atoms 3 region grain basis 1 4 basis 2 4 basis 3 4 basis 4 4 basis 5 3 basis 6 3 basis 7 3 basis 8 3 units box group grain region grain pair_style hybrid eam/fs tersoff lj/cut 10.0 lcbop pair_coeff * * eam/fs Fe-C.eam.fs Fe C NULL NULL #Fe-C pair_coeff * * lcbop C.lcbop NULL C NULL NULL #C-C #pair_coeff 1 2 lj/cut 0.00759200 3.5 9.7625 #Fe-C pair_coeff * * tersoff BNC.tersoff NULL NULL B N #B-N #pair_coeff 1 1 lj/cut 0.105 3.85 10 #Fe-Fe pair_coeff 2 4 lj/cut 0.004068 3.367 #N-C pair_coeff 2 3 lj/cut 0.003311 3.411 #B-C pair_coeff 1 3 lj/cut 0.00284 3.4 #Fe-B pair_coeff 1 4 lj/cut 0.00759200 3.5 9.7625 #Fe-N neighbor 2 bin neigh_modify every 1 delay 5 check yes #bond_style harmonic #bond_coeff 1 80.0 2.0 region 1 block 5 150 5 50 2.5 42.5 units box region 2 block 10 150 10 50 5 40 units box group 1 region 1 group 2 region 2 group workpiece region workpiece group boundary subtract workpiece 1 group hengwen subtract 1 2 group newton subtract workpiece boundary hengwen group workpiece union newton hengwen #min_style sd # 能量最小化模式，sd #minimize 1.0e-20 1.0e-30 20000 20000 # 能量最小化參數(shù)，指數(shù)越大最小化程度越深 #write_restart restart.txt compute 3 all pe/atom # 計(jì)算每個(gè)原子的勢能 compute 4 all ke/atom # 計(jì)算每個(gè)原子的動(dòng)能 compute temperature2 hengwen temp compute temperature3 newton temp/region 2 compute p all pressure thermo_temp # 計(jì)算應(yīng)力，計(jì)算結(jié)果記為p variable pressx equal c_p[1] #定義變量pressx=c_p[1],c_p[1]的意思是p里第一個(gè)值 variable pressy equal c_p[2] #定義變量pressy=c_p[2],c_p[2]的意思是p里第二個(gè)值 compute m workpiece stress/atom # 計(jì)算工件應(yīng)力，計(jì)算結(jié)果記為p compute n workpiece reduce sum c_m[1] c_m[2] c_m[3] variable px equal -(c_n[1])/(3*vol) variable py equal -(c_n[2])/(3*vol) compute k grain stress/atom # 計(jì)算磨粒應(yīng)力，計(jì)算結(jié)果記為p compute g grain reduce sum c_k[1] c_k[2] c_k[3] variable pgx equal -(c_g[1])/(3*vol) variable pgy equal -(c_g[2])/(3*vol) velocity boundary set 0.0 0.0 0.0 velocity hengwen create 293.0 5812775 temp temperature2 units box velocity newton create 293.0 5812775 units box fix 3 all nve fix 4 boundary setforce 0.0 0.0 0.0 fix 5 hengwen temp/rescale 10 293.0 293.0 10.0 1.0 fix 6 grain setforce 0.0 0.0 0.0 compute force grain group/group newton fix 7 grain move linear -1.4 0 0.0 units box # tan=0.025Vw/Vs dump result all custom 100 dump.*.lammpstrj id type mass x y z xs ys zs c_3 c_4 c_5 #dump result workpiece cfg 100 *.cfg id type xs ys zs x y z thermo 100 thermo_style custom step c_temperature3 c_force[1] c_force[2] c_force[3] etotal press v_pressx v_pressy v_px v_py v_pgx v_pgy timestep 0.001 # 步長 0.001fs run 21429 # 運(yùn)行 100 步 |
回復(fù)此樓» 猜你喜歡
歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有24人回復(fù)
-
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生
已經(jīng)有10人回復(fù)
-
招調(diào)劑
已經(jīng)有11人回復(fù)
-
0856材料與化工求調(diào)劑！
已經(jīng)有11人回復(fù)
-
材料與化工304求調(diào)劑
已經(jīng)有10人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有31人回復(fù)
-
中科大材料299求調(diào)劑
已經(jīng)有7人回復(fù)
-
301求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料考研339求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料專碩290求調(diào)劑
已經(jīng)有5人回復(fù)
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
8樓2020-01-10 15:07:23
2樓2019-12-19 07:48:07
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
3樓2019-12-31 14:40:39
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
4樓2020-01-02 15:44:20
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +4 | Bahati 2026-03-05 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 復(fù)試調(diào)劑 +4 | 呼呼？~+123456 2026-03-05 | 7/350 |
|
|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線，總分 281 求調(diào)劑 +4 | 香香 Lu 2026-03-03 | 8/400 |
|
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +6 | 青科11 2026-03-02 | 6/300 |
|
|
|
[考研] 0856材料求調(diào)劑 +17 | hyf hyf hyf 2026-02-28 | 18/900 |
|
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +7 | 1一11 2026-03-03 | 7/350 |
|
|
|
[考研] 環(huán)境調(diào)劑 +7 | 柒槿levana 2026-03-01 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +7 | 哇塞王帥 2026-03-03 | 7/350 |
|
|
|
[考研] 一志愿中石油，化學(xué)工程，301求調(diào)劑 +3 | 豚饅枳枳 2026-03-04 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +4 | 熬夜的貓頭鷹 2026-03-02 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +6 | 等春來， 2026-03-04 | 6/300 |
|
|
|
[考研] 292求調(diào)劑 +9 | yhk_819 2026-02-28 | 9/450 |
|
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
|
[考研] 一志愿中石油（華東）本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-02 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|