24小时热门版块排行榜    

查看: 1019  |  回复: 2

jianying8996

木虫 (著名写手)


[交流] LAMMPS模拟宏观应力可以实现吗?

想模拟层状金属的弯曲或拉伸应力应变,从中得到一些数值变化。不知LAMMPS可否实现?
没用过LAMMPS,希望大家提供一些意见,多谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tudragon

捐助贵宾 (小有名气)


★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
jianying8996: 金币+6, 感谢指教! 2013-09-28 22:42:43
可以做到,常用的有两种方法,第一种是把每层的原子分别选为一个group,然后计算每层的各压强张量,取负值即得应力,这种方法的缺点好像是计算时每层的体积被考虑为初始时的体积,因此有一定误差;第二种是输出每个原子的原子应力张量(都是能量量纲)和伏龙诺依原子体积,然后后处理时对每层的相同取向的原子应力张量进行统计,然后除以该层所有原子的伏龙诺依原子体积之和,这种方法的缺点有二,其一是是否应该考虑位力项存在争议,其二是如果某层存在缺陷或空位,会导致缺陷和空位部分的伏龙诺依原子体积与真实情况相偏离,如果某层缺陷较多,则此法计算结果可能存在较大误差。
2楼2013-09-28 15:24:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianying8996

木虫 (著名写手)


引用回帖:
2楼: Originally posted by tudragon at 2013-09-28 15:24:01
可以做到,常用的有两种方法,第一种是把每层的原子分别选为一个group,然后计算每层的各压强张量,取负值即得应力,这种方法的缺点好像是计算时每层的体积被考虑为初始时的体积,因此有一定误差;第二种是输出每个原 ...

每层原子组为一个group的话就不必考虑每个原子的受力了,体积误差是原子的间距的导致的?
第二个方法我可否加入修正项改过来,应该可以的,谢谢了!
3楼2013-09-28 23:00:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jianying8996 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 uKOSCocCb6Fm 2026-11-12 8/400 2026-11-15 02:26 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 uKOSCocCb6Fm 2026-11-12 10/500 2026-11-15 02:15 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 9/450 2026-11-15 02:04 by xSrQNSSjcgoX
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 16:55 by zzy933268
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 13:02 by ECLTXnJnRSXf
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:22 by HajghKOMynmz
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 9/450 2026-11-12 11:22 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:03 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:00 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:49 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:38 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:30 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 FChKX8meu1ds 2026-11-09 12/600 2026-11-12 09:43 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 5/250 2026-11-12 08:15 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 06:25 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
信息提示
请填处理意见