 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 2/1 | 返回列表 |
| 查看: 7329 | 回復: 14 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
JackieXiao銅蟲 (初入文壇)
|
[求助]
lammps做拉伸試驗，求大神進來指導下（烈日炎炎，有勞了`(*∩_∩*)′） 已有2人參與
|
|
|
lammps輸入文件如下： #input file for uniaxial tensile loading of single crystal Mg # ------------------------ INITIALIZATION ---------------------------- units metal dimension 3 boundary p p p atom_style atomic # ----------------------- ATOM DEFINITION ---------------------------- read_data data.300K # ------------------------ FORCE FIELDS ------------------------------ pair_style meam pair_coeff * * ./library.meam Mg ./Mg.meam Mg neighbor 1.0 bin neigh_modify every 1 delay 0 check yes # ------------------------- SETTINGS --------------------------------- compute csym all centro/atom 12 compute peratom all pe/atom ###################################### # EQUILIBRATION reset_timestep 0 timestep 0.003 fix 1 all npt temp 300 300 1 iso 0 0 1 drag 1 # Set thermo output thermo 100 thermo_style custom step lx ly lz press pxx pyy pzz pe temp # Run for at least 10 picosecond (assuming 1 fs timestep) run 10000 write_data data.300K_relax unfix 1 # Store final cell length for strain calculations variable tmp equal "lx" variable L0 equal ${tmp} print "Initial Length, L0: ${L0}" ###################################### # DEFORMATION reset_timestep 0 fix 1 all npt temp 300 300 1 y 0 0 1 z 0 0 1 drag 1 //此處系綜控制 variable srate equal 1.0e10 variable srate1 equal "v_srate / 1.0e12" fix 2 all deform 1 x erate ${srate1} units box remap x //此處拉伸 # Output strain and stress info to file # for units metal, pressure is in [bars] = 100 [kPa] = 1/10000 [GPa] # p2, p3, p4 are in GPa variable strain equal "(lx - v_L0)/v_L0" variable p1 equal "v_strain" variable p2 equal "-pxx/10000" variable p1 equal "v_strain" variable p2 equal "-pxx/10000" variable p3 equal "-pyy/10000" variable p4 equal "-pzz/10000" fix def1 all print 100 "${p1} ${p2} ${p3} ${p4}" file Al_SC_100.def1.txt screen no # Use cfg for AtomEye dump 1 all custom 5000 dump.tensile.atom id type xs ys zs c_csym c_peratom fx fy fz dump_modify 1 element Mg # Display thermo thermo 100 thermo_style custom step v_strain temp v_p2 v_p3 v_p4 ke pe press run 100000 以上就是輸入文件，存在下面幾點疑問 1、拉伸過程中要采用什么系綜？上面的腳本是用的NPT系綜，但是我看到有的文獻用NVT，有的用NVE，而我用不同的系綜得到的結果差距有點大。 2、周期性邊界條件如何確定，是三維周期性邊界條件還是自由邊界。 自由邊界條件，nve系綜下應力應變曲線 三維周期邊界條件，npt系綜下應力應變曲線 |
回復此樓» 猜你喜歡
【考研調劑】化學專業(yè) 281分，一志愿四川大學，誠心求調劑
已經(jīng)有4人回復
-
317求調劑
已經(jīng)有5人回復
-
0817 化學工程 299分求調劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有10人回復
-
一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調劑。
已經(jīng)有3人回復
-
一志愿福大288有機化學，求調劑
已經(jīng)有3人回復
-
085600材料與化工求調劑
已經(jīng)有6人回復
-
一志愿西安交通大學材料工程專業(yè) 282分求調劑
已經(jīng)有6人回復
-
085600材料與化工調劑 324分
已經(jīng)有10人回復
-
294求調劑材料與化工專碩
已經(jīng)有11人回復
-
一志愿南昌大學，327分，材料與化工085600
已經(jīng)有3人回復
89807dlydly
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 205.8
- 帖子: 11
- 在線: 36.2小時
- 蟲號: 4971048
- 注冊: 2016-09-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
13樓2017-03-23 15:14:55
★ ★
感謝參與，應助指數(shù) +1
JackieXiao(月只藍代發(fā)): 金幣+2, 鼓勵交流！ 2014-07-26 16:50:39
感謝參與，應助指數(shù) +1
JackieXiao(月只藍代發(fā)): 金幣+2, 鼓勵交流！ 2014-07-26 16:50:39
|
本帖內容被屏蔽 |
2樓2014-07-23 16:45:05
JackieXiao
銅蟲 (初入文壇)
- 應助: 1 (幼兒園)
- 金幣: 1.8
- 散金: 50
- 帖子: 26
- 在線: 29.6小時
- 蟲號: 3036740
- 注冊: 2014-03-10
- 專業(yè): 金屬材料的凝固與結晶學
3樓2014-07-23 23:13:20
Saraing
鐵蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 65.5
- 帖子: 13
- 在線: 6.6小時
- 蟲號: 3058337
- 注冊: 2014-03-16
- 性別: MM
- 專業(yè): 金屬材料的微觀結構
|
你好，我跟你做的有些相似，也是做拉伸模擬，不知道你的應力應變曲線是怎么做出來的？ # Output strain and stress info to file # for units metal, pressure is in [bars] = 100 [kPa] = 1/10000 [GPa] # p2, p3, p4 are in GPa variable p1 equal "v_strain" variable p2 equal "-pxx/10000" variable p3 equal "-pyy/10000" variable p4 equal "-pzz/10000" variable p5 equal "-pxy/10000" variable p6 equal "-pxz/10000" variable p7 equal "-pyz/10000" variable p8 equal "pe" variable p9 equal "v_strain_x" fix equil1 all print 1 "${p1} ${p2} ${p3} ${p4} ${p5} ${p6} ${p7} ${p8} ${p9}" file data.1.txt screen no fix equil2 all print 1 "${a}" file data.2.txt 這個是輸出data文件的類型，接下來我應該是怎么算呢？我的程序也是nve的自由周期邊界！ |
4樓2014-08-09 10:16:22
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調劑 324分 +10 | llllkkkhh 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿985，本科211，0817化學工程與技術319求調劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 求材料調劑 +10 | 隔壁陳先生 2026-03-12 | 10/500 |
|
|
|
[考研] 330求調劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
|
[考研] 331求調劑（0703有機化學 +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 304求調劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 考研求調劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 材料與化工求調劑 +6 | 為學666 2026-03-16 | 6/300 |
|
|
|
[考研] 材料工程專碩274一志愿211求調劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 283求調劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 機械專碩325，尋找調劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 318求調劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 321求調劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 070303 總分349求調劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 070305求調劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 復試調劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 328求調劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 311求調劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|