 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (594)
- >考研 (63)
- >蟲(chóng)友互識(shí) (41)
- >導(dǎo)師招生 (22)
- >基金申請(qǐng) (11)
- >考博 (9)
- >教師之家 (8)
- >碩博家園 (8)
- >論文投稿 (8)
- >文獻(xiàn)求助 (5)
- >無(wú)機(jī)/物化 (3)
- >論文道賀祈福 (2)
- >文獻(xiàn)求助完結(jié)子版 (2)
- >找工作 (2)
- >公派出國(guó) (2)
- >復(fù)合材料 (1)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 4746 | 回復(fù): 2 | ||
清風(fēng)明月P6銅蟲(chóng) (初入文壇)
|
[求助]
求一個(gè)LAMMPS里面聚合物適用的拉伸in文件
|
|
最近在做拉伸模擬。自己嘗試著寫(xiě)了一個(gè)in文件，但是運(yùn)行之后確得不到想要的結(jié)果。。。 我用聚乙烯做的測(cè)試，現(xiàn)在MS里面建好的模型，后來(lái)到處轉(zhuǎn)化為lammps文件 在論壇找了好久基本上都是金屬的，求大神給個(gè)高分子可以用的吧。。謝謝了 下面是我的in文件，也不知道問(wèn)題出在哪里。。附件是得到的輸出文件 #------模擬環(huán)境初始化------ boundary p p p #邊界周期性設(shè)定，周期、周期、自由非周期 units real #模型單元模式設(shè)定為real dimension 3 #幾何維度設(shè)定為3D atom_style full #原子模式設(shè)定為full neighbor 2.0 bin #設(shè)定neighbor截?cái)喟霃?力場(chǎng)截?cái)喟霃剑愃芃S中Buffer width） neigh_modify every 1 delay 0 check yes #------微觀模型的建立------ # 定義原子結(jié)構(gòu) pair_style lj/class2/coul/cut 10.0 9.5 #airebo 3.0 1 0 lj/cut/coul/cut 10.0 #設(shè)定PCFF勢(shì)函數(shù)的對(duì)勢(shì)（類似MS中的Nonbond potential）信息：對(duì)勢(shì)包含相互作用為范德華力和庫(kù)侖力，范德華力和庫(kù)侖力截?cái)喟霃剑愃芃S中Cutoff distance） bond_style class2 #設(shè)定Compass勢(shì)函數(shù)的鍵勢(shì) angle_style class2 #設(shè)定Compass勢(shì)函數(shù)的角勢(shì) dihedral_style class2 #設(shè)定Compass勢(shì)函數(shù)的兩面角勢(shì) improper_style class2 #設(shè)定Compass勢(shì)函數(shù)的兩面角扭轉(zhuǎn)勢(shì) read_data ../bench/20140924/PeDeform.lammps05 #------設(shè)置原子相關(guān)參數(shù)------ velocity all create 300 102486 rot yes dist & gaussian units box #設(shè)定所有原子的速率 # 能量?jī)?yōu)化 minimize 1.0e-7 0 1000000 1000000 #設(shè)置minimize，指定能量收斂值為1.0e-7、力矢量在坐標(biāo)軸上的分量長(zhǎng)度收斂值為0、最大迭代步數(shù)為1000000、計(jì)算力長(zhǎng)度/能量的最大次數(shù)為1000000 #------動(dòng)力學(xué)設(shè)置------ fix NPT all npt temp 300.0 300.0 100 iso 1.0 1.0 1000 #對(duì)整體設(shè)定npt系綜，控溫方法為Nose/Hoover，設(shè)定起始、終止溫度和壓力，衰減常數(shù)（溫度一般取100倍timesteps，這里取100.0，壓力一般取1000倍timesteps，這里取1000.0） #拉伸模擬 reset_timestep 0 region 1 block INF INF INF INF INF 3 units box group low region 1 region 2 block INF INF INF INF 72 INF units box group up region 2 group boundary union low up group mobile subtract all boundary velocity up set 0 0 0.004 velocity low set 0 0 -0.004 # 動(dòng)力學(xué)輸出設(shè)置 thermo 1000 #每1000步輸出一次熱動(dòng)力學(xué)信息（溫度、能量、壓力） thermo_style custom step lx ly lz press pxx pyy pzz pe temp thermo_modify lost warn #檢查丟失原子，顯示警報(bào) #fix myrelax all ave/time 1 1000 1000 c_lz c_pxx ../bench/20140924/data/relax_gpc.dat1 timestep 1 #設(shè)定時(shí)間步長(zhǎng)為1fs dump train1 all xyz 1000 ../bench/20140924/g-PEc.xyz #每1000步輸出所有區(qū)間的dcd坐標(biāo)，寫(xiě)入到g_PE_C.xyz文件中 run 100000 #運(yùn)行10ps |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : myout.out
2014-09-28 17:19:56, 11.55 K
» 猜你喜歡
309求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- MS粗�；Ｐ腿绾螌�(dǎo)入lammps
已經(jīng)有10人回復(fù)
- LAMMPS進(jìn)行拉伸模擬的一些問(wèn)題
已經(jīng)有4人回復(fù)
- 如何將amber的軌跡文件mdcrd格式轉(zhuǎn)換成gromacs的軌跡文件格式 trr呢
已經(jīng)有5人回復(fù)
- 新手求助一個(gè)quantum-espresso的使用問(wèn)題！
已經(jīng)有13人回復(fù)
- lammps 中 做拉伸時(shí) sss和ssp邊界條件不同，結(jié)果為什么不一樣？？？？
已經(jīng)有5人回復(fù)
- 硅納米管用lammps做拉伸
已經(jīng)有3人回復(fù)
- 請(qǐng)教一個(gè)C++讀文件的問(wèn)題，
已經(jīng)有12人回復(fù)
- lammps做拉伸和壓縮的應(yīng)力問(wèn)題
已經(jīng)有9人回復(fù)
- 跑完Gromacs的動(dòng)力學(xué)模擬之后的軌跡文件.trr里面沒(méi)有內(nèi)容，為什么？
已經(jīng)有4人回復(fù)
- 求回答聚合物被拉伸時(shí)，在屈服點(diǎn)和應(yīng)變軟化以及應(yīng)變硬化時(shí)內(nèi)部分子鏈?zhǔn)窃趺催\(yùn)動(dòng)的的？
已經(jīng)有9人回復(fù)
- lammps實(shí)驗(yàn)拉伸Al納米線的in文件的錯(cuò)誤提示好久找不出？
已經(jīng)有8人回復(fù)
- 關(guān)于lammps fix npt/deform做拉伸時(shí)溫度的變化
已經(jīng)有6人回復(fù)
- lammps 拉伸應(yīng)力應(yīng)變曲線
已經(jīng)有4人回復(fù)
- windows下lammps編譯
已經(jīng)有5人回復(fù)
- 【求助】DL_POLY的config文件如何生成??
已經(jīng)有98人回復(fù)
tang_wei
金蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 1353.5
- 紅花: 28
- 帖子: 428
- 在線: 81.7小時(shí)
- 蟲(chóng)號(hào): 4067192
- 注冊(cè): 2015-09-12
- 專業(yè): 高分子物理與高分子物理化
2樓2017-12-08 14:30:25
南gege
銀蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 540.7
- 散金: 10
- 帖子: 45
- 在線: 50.3小時(shí)
- 蟲(chóng)號(hào): 3807463
- 注冊(cè): 2015-04-13
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
3樓2020-10-10 10:12:23
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 清風(fēng)明月P6 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 286分 材料求調(diào)劑 +10 | 麻辣魷魚(yú) 2026-03-27 | 11/550 |
|
|---|---|---|---|---|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +9 | 景隨風(fēng) 2026-03-29 | 14/700 |
|
|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +5 | 真名有冰 2026-03-29 | 6/300 |
|
|
[考研] 343求調(diào)劑 +6 | 愛(ài)羈絆 2026-03-29 | 6/300 |
|
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書(shū) 2026-03-23 | 4/200 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 一志愿國(guó)科過(guò)程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|