 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1124 | 回復(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ù)此樓» 猜你喜歡
售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】
已經(jīng)有5人回復(fù)
-
售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】
已經(jīng)有3人回復(fù)
-
接受26屆調(diào)劑生
已經(jīng)有10人回復(fù)
-
售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】
已經(jīng)有8人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有46人回復(fù)
-
化學(xué)工程與技術(shù)316求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料調(diào)劑
已經(jīng)有13人回復(fù)
-
085600材料與化工 292分 求調(diào)劑
已經(jīng)有3人回復(fù)
-
273求調(diào)劑
已經(jīng)有7人回復(fù)
-
334求調(diào)劑
已經(jīng)有3人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Augus537 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 哈爾濱理工大學(xué)2026年研究生調(diào)劑，材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 我嘗試在小木蟲上招計(jì)算機(jī)專業(yè)的調(diào)劑學(xué)生 +3 | littlehu66 2026-03-04 | 5/250 |
|
|
|
[考研] 材料專碩290求調(diào)劑 +8 | 杰尼龜aaa 2026-03-04 | 8/400 |
|
|
|
[考研] 【求調(diào)劑】293分環(huán)境工程求調(diào)劑材料/化工，服從調(diào)劑，抗壓能力強(qiáng)！ +8 | xiiiia 2026-03-04 | 8/400 |
|
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +10 | pin8023 2026-02-28 | 13/650 |
|
|
[考研]
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
 5+3
|
蘇堯幺幺 2026-03-03 | 8/400 |
|
|
|
[論文投稿]
 主編終審卡者不動(dòng)
100+4
|
Stray2021 2026-03-03 | 4/200 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 267化工調(diào)劑求助 +5 | 聰少OZ 2026-03-04 | 5/250 |
|
|
|
[考研] 一志愿中石油，化學(xué)工程，301求調(diào)劑 +3 | 豚饅枳枳 2026-03-04 | 3/150 |
|
|
|
[考研] 282求調(diào)劑 +5 | 2103240126 2026-03-02 | 8/400 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|
|
|
[考研] 化工270求調(diào)劑 +10 | 什么名字qwq 2026-03-02 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|