| 查看: 4825 | 回复: 12 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
lammps得出的应力张量如何可视化？ 已有1人参与
|
||||
|
lammps是可以通过 compute stress/atom 来计算各个原子的应力的。通过dunp，可以把这个应力信息输出到文件。按道理说，应该可以对应力场进行可视化。但是应力张量，有6个分量，可视化好难啊。 哪位大侠有做个这种图啊 感激不尽 |
回复此楼» 收录本帖的淘帖专辑推荐
 资源收集 | 材料 |
» 猜你喜欢
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
申请科研项目纯靠实力能成功吗？
已经有7人回复
-
极限拉扯，如何抉择
已经有14人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助张量如何表示第二应力不变量
已经有3人回复
★ ★
月只蓝: 金币+2, 鼓励交流！ 2014-02-28 18:22:12
月只蓝: 金币+2, 鼓励交流！ 2014-02-28 18:22:12
|
不需要的，就是你每个原子的PXX PYY PZZ的输出值除以原子体积。 我想你做力学单位应该是METAL，所以体积默认单位就是A3， 你PXX的输出值的单位是 BarA3 我举个例子，好比PXX=1000， 单位是BAR， 然后原子体积假定BCC的铁，晶格常数2.87，原子体积就是 2.87^3/2，最后你需要的值是 1000/（2.87^3/2)/10000， 最后那个10000从Bar写到GPa, PYY, PZZ都照做以后， 加起来取负除以3，就可以PLOT了。 如果你算VON MIESE，上面也适用PXY， PYZ和PXZ，就是公式复杂点 够清楚了吧，希望可以帮到你  |
5楼2014-02-28 17:11:01
★ ★ ★
月只蓝: 金币+3, 鼓励交流！ 2014-02-28 08:16:53
月只蓝: 金币+3, 鼓励交流！ 2014-02-28 08:16:53
|
不难的，你可视化应力无非是hydrostatical 或者是 von Mises 正应力 hydro. 就是 -([1]+[2]+[3])/3, 这个应力引起体积变化的 von Miese 的稍微长一点 , 看这个 http://kbs.cnki.net/forums/25971/ShowThread.aspx，等效应力 但是我要提醒你的是，LAMMPS算出的STRESS/ATOM 不是p, 而是pV. 对于每个原子来说V就是原子体积 a^3/(nr of the atom in unit lattice). a是晶格常数。 如果你晶格变化不大还OK，但是如果出现缺陷，那一块的V是无法定义的，只能近似地使用上面的公式。 |
2楼2014-02-27 18:11:38
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
dyc_2008: 金币+200, ★★★★★最佳答案, 非常感谢 2014-02-28 07:27:03
月只蓝: 金币+1, 鼓励交流！ 2014-02-28 18:21:59
感谢参与，应助指数 +1
dyc_2008: 金币+200, ★★★★★最佳答案, 非常感谢 2014-02-28 07:27:03
月只蓝: 金币+1, 鼓励交流！ 2014-02-28 18:21:59
|
也就是说，你最后要输出的就一个值 要么你的正应力 （hys.) 要么等效应力 （von Miese). 取决于你要说明什么，和体积变化相关的用HYS，考虑切应力作用的用等效应力 |
3楼2014-02-27 18:27:37
4楼2014-02-28 07:27:04
10
