| 3 | 1/1 | 返回列表 |
| 查看: 4747 | 回復: 2 | ||
清風明月P6銅蟲 (初入文壇)
|
[求助]
求一個LAMMPS里面聚合物適用的拉伸in文件
|
|
最近在做拉伸模擬。自己嘗試著寫了一個in文件,但是運行之后確得不到想要的結果。。。 我用聚乙烯做的測試,現(xiàn)在MS里面建好的模型,后來到處轉化為lammps文件 在論壇找了好久基本上都是金屬的,求大神給個高分子可以用的吧。。謝謝了 下面是我的in文件,也不知道問題出在哪里。。附件是得到的輸出文件 #------模擬環(huán)境初始化------ boundary p p p #邊界周期性設定,周期、周期、自由非周期 units real #模型單元模式設定為real dimension 3 #幾何維度設定為3D atom_style full #原子模式設定為full neighbor 2.0 bin #設定neighbor截斷半徑-力場截斷半徑(類似MS中Buffer width) neigh_modify every 1 delay 0 check yes #------微觀模型的建立------ # 定義原子結構 pair_style lj/class2/coul/cut 10.0 9.5 #airebo 3.0 1 0 lj/cut/coul/cut 10.0 #設定PCFF勢函數(shù)的對勢(類似MS中的Nonbond potential)信息:對勢包含相互作用為范德華力和庫侖力,范德華力和庫侖力截斷半徑(類似MS中Cutoff distance) bond_style class2 #設定Compass勢函數(shù)的鍵勢 angle_style class2 #設定Compass勢函數(shù)的角勢 dihedral_style class2 #設定Compass勢函數(shù)的兩面角勢 improper_style class2 #設定Compass勢函數(shù)的兩面角扭轉勢 read_data ../bench/20140924/PeDeform.lammps05 #------設置原子相關參數(shù)------ velocity all create 300 102486 rot yes dist & gaussian units box #設定所有原子的速率 # 能量優(yōu)化 minimize 1.0e-7 0 1000000 1000000 #設置minimize,指定能量收斂值為1.0e-7、力矢量在坐標軸上的分量長度收斂值為0、最大迭代步數(shù)為1000000、計算力長度/能量的最大次數(shù)為1000000 #------動力學設置------ fix NPT all npt temp 300.0 300.0 100 iso 1.0 1.0 1000 #對整體設定npt系綜,控溫方法為Nose/Hoover,設定起始、終止溫度和壓力,衰減常數(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 # 動力學輸出設置 thermo 1000 #每1000步輸出一次熱動力學信息(溫度、能量、壓力) thermo_style custom step lx ly lz press pxx pyy pzz pe temp thermo_modify lost warn #檢查丟失原子,顯示警報 #fix myrelax all ave/time 1 1000 1000 c_lz c_pxx ../bench/20140924/data/relax_gpc.dat1 timestep 1 #設定時間步長為1fs dump train1 all xyz 1000 ../bench/20140924/g-PEc.xyz #每1000步輸出所有區(qū)間的dcd坐標,寫入到g_PE_C.xyz文件中 run 100000 #運行10ps |
金蟲 (正式寫手)
銀蟲 (初入文壇)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 286分 材料求調劑 +10 | 麻辣魷魚 2026-03-27 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 一志愿華東師范大學有機化學專業(yè),初試351分,復試被刷求調劑! +5 | 真名有冰 2026-03-29 | 6/300 |
|
|
[考研] 材料專碩調劑 +5 | 椰椰。 2026-03-29 | 5/250 |
|
|
[考研] 0856求調劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
[考研] 0856求調劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
[考研] 289求調劑 +13 | 新時代材料 2026-03-27 | 13/650 |
|
|
[考研] 11408軟件工程求調劑 +3 | Qiu學ing 2026-03-28 | 3/150 |
|
|
[考研] 343求調劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研] 壓國家一區(qū)線,求導師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 一志愿華北電力大學能動專碩,293,求調劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
[考研] 312,生物學求調劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
[考研] 346求調劑 一志愿070303有機化學 +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
[考研] 275求調劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 085600材料與化工調劑 +10 | A-哆啦Z夢 2026-03-23 | 16/800 |
|
|
[考研] 336材料求調劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 總分322求生物學/生化與分子/生物信息學相關調劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 材料專碩 335 分求調劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 上海電力大學材料防護與新材料重點實驗室招收調劑研究生(材料、化學、電化學,環(huán)境) +4 | 我愛學電池 2026-03-23 | 4/200 |
|
|
[考研] 一志愿南航材料專317分求調劑 +5 | 炸呀炸呀炸薯條 2026-03-23 | 5/250 |
|