 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 7337 | 回復(fù): 14 | |||||||
JackieXiao銅蟲 (初入文壇)
|
[求助]
lammps做拉伸試驗，求大神進來指導(dǎo)下（烈日炎炎，有勞了`(*∩_∩*)′） 已有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，而我用不同的系綜得到的結(jié)果差距有點大。 2、周期性邊界條件如何確定，是三維周期性邊界條件還是自由邊界。 自由邊界條件，nve系綜下應(yīng)力應(yīng)變曲線 三維周期邊界條件，npt系綜下應(yīng)力應(yīng)變曲線 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 量化,第一性,MD筆記 | 分子動力學(xué) | 學(xué)習(xí)與研究 | 分子動力學(xué)模擬 |
| 拉抻 |
» 猜你喜歡
一志愿華中科技大學(xué)071000，求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)一志愿985，分數(shù)349求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有4人回復(fù)
-
085600材料與化工306
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有4人回復(fù)
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有8人回復(fù)
-
資源與環(huán)境 調(diào)劑申請(333分)
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- LAMMPS進行拉伸模擬的一些問題
已經(jīng)有4人回復(fù)
- 哪位大神能告訴我LAMMPS可以模擬聚合物的聚合過程嗎？
已經(jīng)有4人回復(fù)
- lammps 中 做拉伸時 sss和ssp邊界條件不同，結(jié)果為什么不一樣？？？？
已經(jīng)有5人回復(fù)
- 硅納米管用lammps做拉伸
已經(jīng)有3人回復(fù)
- 請教大神如何在64bit的win7系統(tǒng)下正確安裝Lammps？計算機是4核的處理器
已經(jīng)有14人回復(fù)
- 求大神幫忙分析以下拉伸斷口掃描圖片
已經(jīng)有6人回復(fù)
- Origin處理拉伸曲線！求大神解決
已經(jīng)有8人回復(fù)
- lammps實驗拉伸Al納米線的in文件的錯誤提示好久找不出？
已經(jīng)有8人回復(fù)
- lammps 拉伸應(yīng)力應(yīng)變曲線
已經(jīng)有4人回復(fù)
chinaboy798
銀蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 299.3
- 散金: 20
- 紅花: 8
- 帖子: 177
- 在線: 69.6小時
- 蟲號: 2007979
- 注冊: 2012-09-17
- 性別: GG
- 專業(yè): 金屬材料的力學(xué)行為
【答案】應(yīng)助回帖
|
這個問題我已經(jīng)困擾很久了，在這里我說下我的一些已經(jīng)測試過的結(jié)果，也希望大家能和我討論下： 1.用fix npt系綜，就只能在x y z都是P P P的周期性邊界條件下做拉伸，那么對于速度拉伸法就不能用了，所以只能用fix deform命令，而fix deform命令的拉伸原理是變形盒子，再讓所有原子的coords有一個remap來實現(xiàn)（remap x），我的理解是原子跟著盒子重新分布位置？不管怎么說，這種方法下原子截面是不會變化的，而且體積是逐漸增大的，不符合物理實際，對于很多文獻上用的此法，我一直保持懷疑態(tài)度，但不知作者們作何合理解釋。 2.采用速度方法來實現(xiàn)拉伸，這種情況下就必須設(shè)置拉伸方向的邊界為非周期條件（其他兩個方向不做要求），然后固定兩端邊界一定層厚的原子（上下兩層）不動，再施加某一層原子有一個移動速度來實現(xiàn)拉伸。這里又有一個問題，是先讓兩層原子不動以后再做平衡，還是做了平衡以后再選定邊界兩層原子以實現(xiàn)拉伸？后者顯然不太合理，而前者操作起來就不容易了，因為邊界層原子setfoece 0 0 0 以后做平衡拉伸方向的應(yīng)力無法消除，在我的測試中，隨后的拉伸初始應(yīng)力很大. 希望大神指導(dǎo)！ |
7樓2014-11-11 15:05:02
JackieXiao
銅蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1.8
- 散金: 50
- 帖子: 26
- 在線: 29.6小時
- 蟲號: 3036740
- 注冊: 2014-03-10
- 專業(yè): 金屬材料的凝固與結(jié)晶學(xué)
5樓2014-08-10 08:27:07
Ergon
金蟲 (著名寫手)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 1332.5
- 散金: 220
- 紅花: 5
- 沙發(fā): 2
- 帖子: 1553
- 在線: 105.5小時
- 蟲號: 3295777
- 注冊: 2014-06-27
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
10樓2015-06-23 21:36:58
★ ★
感謝參與，應(yīng)助指數(shù) +1
JackieXiao(月只藍代發(fā)): 金幣+2, 鼓勵交流！ 2014-07-26 16:50:39
感謝參與，應(yīng)助指數(shù) +1
JackieXiao(月只藍代發(fā)): 金幣+2, 鼓勵交流！ 2014-07-26 16:50:39
|
本帖內(nèi)容被屏蔽 |
2樓2014-07-23 16:45:05
JackieXiao
銅蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1.8
- 散金: 50
- 帖子: 26
- 在線: 29.6小時
- 蟲號: 3036740
- 注冊: 2014-03-10
- 專業(yè): 金屬材料的凝固與結(jié)晶學(xué)
3樓2014-07-23 23:13:20
Saraing
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 65.5
- 帖子: 13
- 在線: 6.6小時
- 蟲號: 3058337
- 注冊: 2014-03-16
- 性別: MM
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
|
你好，我跟你做的有些相似，也是做拉伸模擬，不知道你的應(yīng)力應(yīng)變曲線是怎么做出來的？ # 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文件的類型，接下來我應(yīng)該是怎么算呢？我的程序也是nve的自由周期邊界！ |
4樓2014-08-09 10:16:22
fcwu
木蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 3323.9
- 散金: 30
- 紅花: 2
- 帖子: 234
- 在線: 142.8小時
- 蟲號: 3366411
- 注冊: 2014-08-14
- 性別: GG
- 專業(yè): 金屬材料的力學(xué)行為
6樓2014-08-14 19:35:40
Saraing
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 65.5
- 帖子: 13
- 在線: 6.6小時
- 蟲號: 3058337
- 注冊: 2014-03-16
- 性別: MM
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
8樓2014-12-23 10:49:25
9樓2015-03-04 16:54:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 JackieXiao 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 311求調(diào)劑 +13 | 冬十三 2026-03-15 | 14/700 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
|
[考研] 281求調(diào)劑（0805） +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 機械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|