| 查看: 6799 | 回复: 4 | |||||
[求助]
lammps计算弹性模量应力应变 已有1人参与
|
|
大家好，我现在在用lammps做弹性模量，在计算的过程中我发现有两个计算应力的，一个是通过 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学模拟 | 拉抻 | 精华 |
» 猜你喜欢
293求调剂
已经有5人回复
-
290求调剂
已经有10人回复
-
285化工学硕求调剂（081700）
已经有8人回复
-
211本，11408一志愿中科院277分，曾在中科院自动化所实习
已经有5人回复
-
材料与化工304求B区调剂
已经有8人回复
-
有没有大佬发小论文能带我个二作
已经有4人回复
-
考研调剂
已经有5人回复
-
267一志愿南京工业大学0817化工求调剂
已经有6人回复
-
293求调剂
已经有4人回复
-
302求调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 拉伸曲线下的面积计算
已经有14人回复
- 求助固溶处理前后应力应变变化
已经有4人回复
- 如何分析真应力应变曲线？
已经有6人回复
- lammps得出的应力张量如何可视化？
已经有12人回复
- Origin 拟合应力应变曲线获得应变硬化指数n
已经有6人回复
- 在受力不知道的情况下且为弹塑性材料，如何通过应变曲线得到应力？？？
已经有8人回复
- deform中应力应变
已经有5人回复
- lammps 拉伸应力应变曲线
已经有4人回复
- 请问为什么拉伸比压缩更容易破坏材料？
已经有10人回复
- 如何得出应力-应变曲线
已经有14人回复
- castep计算类金刚石结构的应力应变曲线
已经有18人回复
- 站里关于MS Analysis计算材料的应力应变与Cubic以及Trigonal cell弹性常数的讨论
已经有34人回复
|
本帖内容被屏蔽 |
2楼2015-06-30 16:32:57
3楼2015-07-01 10:43:23
|
本帖内容被屏蔽 |
4楼2015-07-01 16:46:20
honest199174 
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 77.5
- 散金: 5
- 红花: 1
- 帖子: 24
- 在线: 27.4小时
- 虫号: 3995677
- 注册: 2015-07-29
- 性别: GG
- 专业: 固体力学
5楼2017-05-10 14:39:39
