| 5 | 1/1 | 返回列表 |
| 查看: 6802 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 329求調(diào)劑 +5 | 想上學(xué)吖吖 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
[考研] 331求調(diào)劑(0703有機化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|