| 5 | 1/1 | 返回列表 |
| 查看: 6801 | 回復(fù): 4 | |||||
changmao新蟲 (初入文壇)
|
[求助]
lammps計算彈性模量應(yīng)力應(yīng)變 已有1人參與
|
|
大家好,我現(xiàn)在在用lammps做彈性模量,在計算的過程中我發(fā)現(xiàn)有兩個計算應(yīng)力的,一個是通過 compute stress 另一個是press/v 我做的是y方向的拉伸,那我得應(yīng)力是直接取stress或者press/v,還是用他們Y方向的分量?我只有10個金幣,把我全部的金幣都拿出來了,求大家指點,非常感激!我in文件如下 units metal boundary p s p atom_style atomic lattice bcc 3.165 orient x 1 0 0 orient y 0 1 1 orient z 0 -1 1 region box block 0 10 0 10 0 5 create_box 3 box create_atoms 1 box pair_style eam/alloy pair_coeff * * W_zhou.eam.alloy W neighbor 0.3 bin neigh_modify delay 0 every 1 check yes region lower block INF INF INF 0.9 INF INF region top block INF INF INF INF 8.9 INF group lower region lower group top region top group boundary union lower upper group mobile subtract all boundary set group lower type 2 set group upper type 3 velocity lower set 0.0 0.0 0.0 compute new3d mobile temp compute new2d mobile temp/partial 1 0 1 compute cna all cna/atom 3.0 compute stress all stress/atom compute stressx all reduce sum c_stress[1] compute stressy all reduce sum c_stress[2] compute stressz all reduce sum c_stress[3] # equilibrate velocity mobile create 300.0 85277 temp new3d fix 1 all nve fix 2 upper setforce 0.0 0.0 0.0 fix 4 lower setforce 0.0 0.0 0.0 fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0 fix_modify 3 temp new3d thermo 50 thermo_modify temp new3d timestep 0.001 run 100 variable tmp equal "ly" variable L0 equal ${tmp} variable strain equal "(ly-v_L0)/ly" # shear velocity upper set 0 0.1 0 velocity mobile ramp vy 0 0.1 y 0.9 8.9 sum yes unfix 3 fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0 fix_modify 3 temp new2d variable p1 equal "v_strain" variable p2 equal –{pzz/(3.165^3/2)}/10000 variable p3 equal –{pzz/(3.165^3/2)}/10000 variable p4 equal –{pzz/(3.165^3/2)}/10000 variable pm equal "" –{(pxx+pyy+pzz0/(3.165^3/2)}/10000 variable st equal "step" fix def1 all print 100 "${st} ${p1} ${pm} ${p2} ${p3} ${p4}" file tensile_d.txt screen no dump 1 all cfg 100 tensile_*.cfg mass type xs ys zs c_cna c_stress[1] c_stress[2] c_stress[3] thermo 100 thermo_modify temp new2d reset_timestep 0 run 100000 |
分子動力學(xué)模擬 | 拉抻 | 精華 |
|
本帖內(nèi)容被屏蔽 |
新蟲 (初入文壇)
|
本帖內(nèi)容被屏蔽 |
銅蟲 (初入文壇)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研]
|
.6lL 2026-03-18 | 6/300 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|