24小时热门版块排行榜    

查看: 2800  |  回复: 7

排骨饭oo

铁虫 (初入文坛)

[求助] 弹性模量的提取 已有1人参与

下图是一种用LAMMPS模拟的聚合物拉伸曲线,曲线振荡比较大,请问用什么方法提取弹性模量,请大神指教!!

弹性模量的提取
聚合物拉伸应力应变曲线.png
回复此楼

» 收录本帖的淘帖专辑推荐

分子动力学模拟

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
排骨饭oo: 金币+10, ★★★很有帮助 2016-06-08 10:29:40
你应该拉伸一点,弛豫一段时间,求出一个压力平均值和误差,再拉伸一点,弛豫一段时间,求出一个压力平均值和误差。每次弛豫时间要足够长,比如十万步,才能得到精确的统计值。lammps我不会用,但这个原理我是懂得。你的涨落大,是因为没有做统计平均。

[ 发自手机版 http://muchong.com/3g ]
2楼2016-06-07 02:03:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

建议每次拉伸千分之一,每次弛豫十万步。不用拉很长,一共拉百分之一或者百分之二就可以求弹性模量了。

[ 发自手机版 http://muchong.com/3g ]
3楼2016-06-07 02:05:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

排骨饭oo

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by brucefan at 2016-06-07 02:05:54
建议每次拉伸千分之一,每次弛豫十万步。不用拉很长,一共拉百分之一或者百分之二就可以求弹性模量了。

非常感谢,你的建议很有参考性!
4楼2016-06-08 10:26:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

biancheng159

铜虫 (正式写手)

俺也发现msi2lmp后不易继续用lammps模拟
例如,不知是应该自定义pair_coeff等呢(转换后的data中已含各类_coeff),还是从外置力场中读入。
楼主若能搞清合适的in的写法,也贴来学习下。

PE的拉伸模拟input(使用了dreiding,但不是由msi2lmp产生的data)
https://icme.hpc.msstate.edu/med ... rphous_Polyethylene
so,hu r u?
5楼2016-07-14 10:49:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青丘白浅

新虫 (初入文坛)

您好,楼主,请问这个问题最后怎么解决的,能请教一下吗?
6楼2017-12-04 11:22:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青丘白浅

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by brucefan at 2016-06-07 02:03:11
你应该拉伸一点,弛豫一段时间,求出一个压力平均值和误差,再拉伸一点,弛豫一段时间,求出一个压力平均值和误差。每次弛豫时间要足够长,比如十万步,才能得到精确的统计值。lammps我不会用,但这个原理我是懂得。 ...

您好老师,想请教一下压力平均值怎么获取,是在驰豫阶段获取吗?如果我用npt系综做驰豫保持一个大气压,那么驰豫时做压力的统计平均好像没有意义啊。
7楼2018-06-19 22:09:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

引用回帖:
7楼: Originally posted by 青丘白浅 at 2018-06-19 22:09:41
您好老师,想请教一下压力平均值怎么获取,是在驰豫阶段获取吗?如果我用npt系综做驰豫保持一个大气压,那么驰豫时做压力的统计平均好像没有意义啊。...

平均就是时间平均,一般取每一个阶段的后一部分的数据平均。拉伸时可以在垂直方向控压。你统计的应该是拉伸方向的压强。

[ 发自手机版 http://muchong.com/3g ]
8楼2018-06-19 23:40:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 排骨饭oo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 13/650 2026-11-09 06:50 by l48hFd8Nwae2
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 21:29 by l48hFd8Nwae2
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见