| 5 | 1/1 | 返回列表 |
| 查看: 2078 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問題,作者歐陽肥腸將贈(zèng)送您 15 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
歐陽肥腸新蟲 (初入文壇)
|
[求助]
瀝青/二氧化硅復(fù)合材料界面拉伸的分子動(dòng)力學(xué)模擬
|
|
|
最近在做瀝青和二氧化硅的界面拉拔模擬中,遇到了一個(gè)問題,就是固定底下部分瀝青,拉拔二氧化硅整體。結(jié)果當(dāng)二氧化硅發(fā)生二埃左右的位移時(shí),又會(huì)被瀝青的力給“扯回去”。請問一下這要怎么解決? 這個(gè)是全部in文件: # Initialization units real boundary p p s atom_style full neighbor 0.3 bin #Dreiding potential information bond_style harmonic angle_style harmonic dihedral_style opls pair_style lj/cut/coul/cut 9.5 9.5 pair_modify mix geometric improper_style cvff special_bonds lj/coul 0.0 0.0 0.5 read_data Layer.data #Area division1 region sio2 block INF INF INF INF 62.0 INF units box region base block INF INF INF INF 70 INF units box group sio2 region sio2 group base region base group asphalt subtract all sio2 #minimization min_style cg minimize 1e-10 1e-10 5000 5000 #Equilibration Stage1 (Langevin dynamics at 500 K) velocity all create 500.0 1231 fix 1 all nve/limit 0.05 fix 2 all langevin 500.0 500.0 10.0 904297 thermo_style custom step temp press thermo 100 timestep 1 run 100000 unfix 1 unfix 2 write_data after.data1 #Equilibration Stage2 fix 1 all nvt temp 500.0 500.0 50 fix 2 all momentum 1 linear 1 1 1 reset_timestep 0 thermo_style custom step temp press thermo 100 timestep 1 run 100000 unfix 1 unfix 2 write_data after.data2 #Equilibration Stage3 fix 1 all npt temp 500.0 300.0 50 x 1 1 1000 y 1 1 1000 fix 2 all momentum 1 linear 1 1 1 reset_timestep 0 thermo_style custom step temp press thermo 100 timestep 1 run 100000 unfix 1 unfix 2 write_data after.data3 #Equilibration Stage4 fix 1 all npt temp 300.0 300.0 50 x 1 1 1000 y 1 1 1000 fix 2 all momentum 1 linear 1 1 1 reset_timestep 0 thermo_style custom step temp press epair ebond eangle edihed pxx pyy pzz lx ly lz thermo 100 timestep 1 run 100000 unfix 1 unfix 2 write_data after.data4 write_restart restart.equil #Area division2 reset_timestep 0 region lower block INF INF INF INF INF 30 units box group lower region lower group mobile subtract all sio2 lower #changebox to larger change_box all z delta 0 100 units box #Compute settings compute 1 mobile stress/atom NULL compute fzz mobile reduce sum c_1[3] compute v mobile voronoi/atom only_group compute vv mobile reduce sum c_v[1] variable sigma equal c_fzz/(c_vv*10) compute d sio2 displace/atom compute dd sio2 reduce ave c_d[3] thermo_style custom step c_dd v_sigma thermo 10 #Pull-off test velocity sio2 set 0.0 0.0 1e-4 fix 1 lower setforce 0.0 0.0 0.0 velocity lower set 0.0 0.0 0.0 fix 2 all nvt temp 300 300 50 dump 1 all atom 100 tension.xyz run 20000 |
新蟲 (小有名氣)
新蟲 (初入文壇)
新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 271求調(diào)劑 +6 | 月色c 2026-03-05 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 復(fù)試調(diào)劑 +4 | 呼呼?~+123456 2026-03-05 | 6/300 |
|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線,總分 281 求調(diào)劑 +4 | 香香 Lu 2026-03-03 | 8/400 |
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
[考研] 紡織、生物、化學(xué)、材料等專業(yè) +3 | Eember. 2026-03-05 | 7/350 |
|
|
[考研] 304分材料專碩求調(diào)劑 +11 | qiuzhigril 2026-03-03 | 14/700 |
|
|
[考研] 0856材料與化工,270求調(diào)劑 +17 | YXCT 2026-03-01 | 20/1000 |
|
|
[考研] 材料與化工,291分,求調(diào)劑 +9 | 咕嚕咕嚕123123 2026-03-04 | 11/550 |
|
|
[考研] 材料學(xué)碩080500復(fù)試調(diào)劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
[考研] 0703化學(xué) 一志愿東北大學(xué) 279分求調(diào)劑 +4 | 聶樂平 2026-03-03 | 4/200 |
|
|
[考研] 能動(dòng)297求調(diào)劑,本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
[考研] 264求調(diào)劑 +3 | thext 2026-03-03 | 3/150 |
|
|
[考研] 環(huán)境工程學(xué)碩288求助調(diào)劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
[考研] 26考研報(bào)考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
[考研] 324求調(diào)劑 +4 | wxz2 2026-03-03 | 5/250 |
|
|
[論文投稿]
通訊作者寫誰,問題是你意想不到的問題
15+3
|
阿爾法啊 2026-03-01 | 3/150 |
|
|
[考研] 085600材料工程一志愿中科大總分312求調(diào)劑 +9 | 吃宵夜1 2026-02-28 | 11/550 |
|