24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4080  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dyc_2008

捐助贵宾 (著名写手)

[求助] lammps得出的应力张量如何可视化?已有1人参与

lammps是可以通过 compute stress/atom 来计算各个原子的应力的。通过dunp,可以把这个应力信息输出到文件。按道理说,应该可以对应力场进行可视化。但是应力张量,有6个分量,可视化好难啊。

哪位大侠有做个这种图啊

感激不尽
回复此楼
一布福田,二修菩提
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

凉yeah

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by werkstoff at 2014-02-27 18:11:38
不难的,你可视化应力无非是hydrostatical 或者是 von Mises
正应力 hydro. 就是 -(++)/3, 这个应力引起体积变化的
von Miese 的稍微长一点 , 看这个
http://kbs.cnki.net/forums/25971/ShowThread.aspx,等效应 ...

你好,请问 我计算应力用的是compute stress/atom,我想做一个应力的图,lammps怎么做输出?是输出一个TXT文件吗?
11楼2016-03-19 22:05:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

werkstoff

银虫 (初入文坛)

★ ★ ★
月只蓝: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

werkstoff

银虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dyc_2008

捐助贵宾 (著名写手)

引用回帖:
3楼: Originally posted by werkstoff at 2014-02-27 18:27:37
也就是说,你最后要输出的就一个值
要么你的正应力 (hys.)
要么等效应力 (von Miese).
取决于你要说明什么,和体积变化相关的用HYS,考虑切应力作用的用等效应力

我是力学盲。关于这个正应力,= -([1]+[2]+[3])/3不。
用这个公式前
要将张量【也就是矩阵】对角化不?
一布福田,二修菩提
4楼2014-02-28 07:27:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见