 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1112 | 回復(fù): 0 | ||
| 【懸賞金幣】回答本帖問題，作者Augus537將贈(zèng)送您 20 個(gè)金幣 | ||
Augus537新蟲 (初入文壇)
|
[求助]
lammps設(shè)置應(yīng)力拉伸的時(shí)候，應(yīng)變卻一直是負(fù)數(shù)？
|
|
|
如題所示，我在用lammps進(jìn)行y軸單軸加卸載的時(shí)候，明明施加的是拉伸力（4GPa），但是應(yīng)變卻一直負(fù)向增大到了-0.05是什么原因？ # 讀取結(jié)構(gòu)+勢(shì)函數(shù)+能量最小化 clear boundary p p p units metal atom_style atomic dimension 3 timestep 0.002 read_data NiTi_100K.data pair_style meam pair_coeff * * library.meam Ni Ti NiTi.meam Ni Ti neighbor 2.0 bin neigh_modify every 10 delay 0 check yes # ------------ 能量最小化 --------------- fix relax all box/relax iso 0.0 vmax 0.001 min_style cg minimize 1e-25 1e-25 10000 10000 unfix relax # 弛豫過程 fix 1 all npt temp 100 100 0.2 iso 0.0 0.0 2.0 drag 2.0 thermo 1000 run 40000 unfix 1 # 設(shè)置加卸載應(yīng)力 compute atompe all pe/atom compute atomke all ke/atom compute atomcentro all centro/atom bcc compute cna all cna/atom 3.64 #加卸載過程應(yīng)變 variable tmp equal "ly" variable L0 equal ${tmp} print "Initial Length, L0: ${L0}" variable p1 equal "(ly-v_L0)/v_L0" variable p2 equal "-pxx/10000" variable p3 equal "-pyy/10000" variable p4 equal "-pzz/10000" variable pe equal "pe" variable temp equal "temp" #輸出 reset_timestep 0 fix stress all ave/time 1 100 200 v_p3 file Stress5de9cf.txt fix strain all ave/time 1 100 200 v_p1 file Strain5de9cf.txt fix def1 all print 200 "${p1} ${p3} ${p2} ${p4}" file NiTi_def1.txt screen no # Display thermo thermo 500 thermo_style custom step v_p1 temp lx ly lz vol density ke pe etotal enthalpy epair press v_p2 v_p3 v_p4 #----------------------------Loading----------------------------------- variable s_xx equal -3.0 print "set on command line s_xx=${s_xx}" #x-force on rigid body variable stmp equal ${s_xx}*10000 fix load all npt temp 100 100 0.2 y 0.0 ${stmp} 2.0 x 0.0 0.0 2.0 z 0.0 0.0 2.0 nreset 10 dump 1 all custom 500 dump_loading_110K.lammpstrj id type x y z fx fy fz c_atompe c_atomke write_data loading_110K.data run 40000 unfix load undump 1 #----------------------------Unloading----------------------------------- print "stress at the end of loading, p3: ${p3}" print "unloading" fix unload all npt temp 100 100 0.2 y ${stmp} 0.0 2.0 x 0.0 0.0 2.0 z 0.0 0.0 2.0 nreset 10 dump 2 all custom 500 dump_unloading_110K.lammpstrj id type x y z fx fy fz c_atompe c_atomke write_data unloading_110K.data run 40000 unfix unload undump 2 unfix def1 |
回復(fù)此樓» 猜你喜歡
085701環(huán)境工程 求調(diào)劑
已經(jīng)有3人回復(fù)
-
264求調(diào)劑
已經(jīng)有7人回復(fù)
-
085600 材料與化工 298
已經(jīng)有8人回復(fù)
-
266求調(diào)劑
已經(jīng)有5人回復(fù)
-
0854總分272
已經(jīng)有6人回復(fù)
-
EST拒稿重投
已經(jīng)有3人回復(fù)
-
歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考
已經(jīng)有11人回復(fù)
-
0703化學(xué)學(xué)碩 288分
已經(jīng)有4人回復(fù)
-
接收調(diào)劑
已經(jīng)有23人回復(fù)
-
化學(xué) 0703求調(diào)劑 總分293 一志愿211
已經(jīng)有3人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Augus537 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +8 | pin8023 2026-02-28 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +5 | 1一11 2026-03-03 | 5/250 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
|
[考研] 中國林科院林化所（南京）2026年招收化學(xué)/材料/環(huán)境工程等背景碩士研究生3名 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
|
[考研] 085600材料與化工 298 調(diào)劑 +3 | 小西笑嘻嘻 2026-03-03 | 3/150 |
|
|
|
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
|
[考研] 材料工程269求調(diào)劑 +5 | 白刺玫 2026-03-02 | 5/250 |
|
|
|
[考研] 化工專碩348，一志愿985求調(diào)劑 +8 | 弗格個(gè) 2026-02-28 | 11/550 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +15 | February_Feb 2026-03-01 | 17/850 |
|
|
|
[考研] 材料類考研調(diào)劑 +6 | gemmgemm 2026-03-01 | 7/350 |
|
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
|
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5 | awaystay 2026-03-02 | 6/300 |
|
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
|
[考研]
 高分子化學(xué)與物理調(diào)劑
+6
|
好好好1233 2026-02-28 | 16/800 |
|
|
|
[考研] 材料284求調(diào)劑，一志愿鄭州大學(xué)英一數(shù)二專碩 +15 | 想上岸的土撥鼠 2026-02-28 | 15/750 |
|
|
|
[考研] 材料與化工328求調(diào)劑 +3 | 。，。，。，。i 2026-03-02 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 熬夜的貓頭鷹 2026-03-02 | 3/150 |
|
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|