| 5 | 1/1 | 返回列表 |
| 查看: 6816 | 回復: 4 | |||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||
changmao新蟲 (初入文壇)
|
[求助]
lammps計算彈性模量應力應變 已有1人參與
|
||||
|
大家好,我現(xiàn)在在用lammps做彈性模量,在計算的過程中我發(fā)現(xiàn)有兩個計算應力的,一個是通過 compute stress 另一個是press/v 我做的是y方向的拉伸,那我得應力是直接取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 |
分子動力學模擬 | 拉抻 | 精華 |
新蟲 (初入文壇)
|
本帖內容被屏蔽 |
|
本帖內容被屏蔽 |
銅蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調劑 +9 | 小熊~~ 2026-03-31 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 311求調劑 +7 | 李芷新1 2026-03-31 | 7/350 |
|
|
[考研] 08開頭看過來。! +3 | wwwwffffff 2026-03-31 | 5/250 |
|
|
[考研] 本科211生物醫(yī)學工程085409求調劑339分 +7 | 里子木yy 2026-03-29 | 7/350 |
|
|
[考研] 0703化學 +20 | 妮妮ninicgb 2026-03-27 | 20/1000 |
|
|
[考研] 272求調劑,接受跨專業(yè)調劑! +3 | 閑魚盧 2026-03-31 | 3/150 |
|
|
[考研] 313求調劑 +6 | 賣個關子吧 2026-03-31 | 6/300 |
|
|
[考研] 335求調劑 +3 | 321* 2026-03-31 | 3/150 |
|
|
[考研] 一志愿中海洋320化學工程與技術學碩求調劑 +8 | 披星河 2026-03-30 | 8/400 |
|
|
[考研] 本科211總分289,08工學真心求調劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[考研] 332求調劑 +14 | 032500 2026-03-25 | 14/700 |
|
|
[考研] 289求調劑 +16 | 新時代材料 2026-03-27 | 16/800 |
|
|
[考研] 292求調劑 +13 | 是妍子也是研子 2026-03-30 | 13/650 |
|
|
[考研] 318一志愿吉林大學生物與醫(yī)藥 求調劑 +5 | 篤行致遠. 2026-03-28 | 5/250 |
|
|
[考研] 290求調劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
[考研] 085602 化工專碩 338分 求調劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 330一志愿中國海洋大學 化學工程 085602 有讀博意愿 求調劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 一志愿上海理工能源動力(085800)310分求調劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 0703化學338求調劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 一志愿吉林大學材料與化工303分求調劑 +4 | 為學666 2026-03-24 | 4/200 |
|