 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1328 | 回復(fù): 0 | ||
xdr金蟲 (正式寫手)
|
[求助]
聚合物拉伸
|
|
之前在lammps網(wǎng)站上看到一個模擬聚合物拉伸的例子，但是運行之后發(fā)現(xiàn)應(yīng)力值有問題，試著修改了很多參數(shù)還是沒有找到原因，請各位指點： in文件 # Initialization units real boundary p p p atom_style molecular read_data data-lammps.txt # Dreiding potential information neighbor 0.4 bin neigh_modify every 10 one 10000 bond_style harmonic bond_coeff 1 350 1.53 angle_style harmonic angle_coeff 1 60 109.5 dihedral_style multi/harmonic dihedral_coeff 1 1.73 -4.49 0.776 6.99 0.0 pair_style lj/cut 10.5 pair_coeff 1 1 0.112 4.01 10.5 # Equilibration Stage 4 (NPT dynamics at 100 K) fix 1 all npt temp 100 100 50 iso 0 0 1000 drag 2 fix 2 all momentum 1 linear 1 1 1 thermo_style custom step temp press epair ebond eangle edihed pxx pyy pzz lx ly lz thermo 1000 timestep 0.5 reset_timestep 0 run 50000 unfix 1 unfix 2 # Uniaxial Tensile Deformation run 0 fix 1 all npt temp 100 100 50 y 0 0 1000 z 0 0 1000 drag 2 fix 2 all deform 1 x erate 1e-5 units box remap x thermo_style custom step temp pxx pyy pzz lx ly lz epair ebond eangle edihed thermo 100 timestep 1 reset_timestep 0 dump 1 all atom 1000 dump1.lammpstrj run 171800 data文件 100 atoms 99 bonds 98 angles 97 dihedrals 1 atom types 1 bond types 1 angle types 1 dihedral types 0.0000 158.5000 xlo xhi 0.0000 158.5000 ylo yhi 0.0000 100.0000 zlo zhi Masses 1 14.02 Atoms 1 1 1 5.6240 5.3279 51.6059 ……（省去） Bonds 1 1 1 2 ……（省去） Angles 1 1 1 2 3 ……（省去） Dihedrals 1 1 1 2 3 4 ……（省去） log文件（只列出拉伸的部分數(shù)據(jù)） Step Temp Pxx Lx 0 99.789705 0.38065973 275.47105 100 101.88813 0.5177075 275.74652 200 107.26163 0.31814294 276.02199 300 91.681759 -0.054211006 276.29746 400 104.34186 -0.30922012 276.57293 500 93.074487 0.80893489 276.8484 600 93.14631 0.67074702 277.12387 700 95.795549 0.79213256 277.39934 800 85.602899 0.08175127 277.67482 900 101.17604 0.31384049 277.95029 1000 109.7724 0.1937424 278.22576 1100 98.72743 0.11628565 278.50123 1200 112.30888 -0.3747989 278.7767 1300 98.345275 -0.73043229 279.05217 1400 97.97989 0.54237926 279.32764 看到一些例子，都是金屬或者納米管、石墨烯類似，atom_style atomic，有人試過molecular的嗎？ |
回復(fù)此樓» 猜你喜歡
283求調(diào)劑（080500）
已經(jīng)有3人回復(fù)
-
286求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分
已經(jīng)有4人回復(fù)
-
308求調(diào)劑
已經(jīng)有6人回復(fù)
-
275求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
287求調(diào)劑
已經(jīng)有9人回復(fù)
-
274求調(diào)劑
已經(jīng)有23人回復(fù)
-
085602 化工專碩 338分 求調(diào)劑
已經(jīng)有6人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有16人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xdr 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +4 | 祁祺祺 2026-03-25 | 4/200 |
|
|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|
|
[考研] 085600材料與化工306 +10 | z1z2z3879 2026-03-21 | 11/550 |
|
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +10 | 77chaselx 2026-03-20 | 10/500 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +13 | Faune 2026-03-21 | 13/650 |
|
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|