计算化学公社

 找回密码 Forget password
 注册 Register

deform计算聚合物应力应变出现问题

查看数: 1884 | 评论数: 15 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-10-27 21:14

正文摘要:

使用MD模拟计算应力应变有两种方式,一种是采用velocity,一种是采用deform方式因此,我想使用deform的命令计算聚合物的应力应变之间的关系 应变是通过(L-L0)/L0计算,而应力则是Pres-zz(bar转换成MPa应除10) 但在 ...

回复 Reply

liuxdhs 发表于 Post on 3 day ago
tangwenke 发表于 2024-11-20 10:04
如果设置负值  不就是变成压缩了吗

压缩xy方向,相当于计算z方向的应力应变
tangwenke 发表于 Post on 3 day ago
liuxdhs 发表于 2024-11-18 09:34
deform的设置可能出现问题,方向上。
可参考以下mdp设置。
Pcoupl     = Berendsen

如果设置负值  不就是变成压缩了吗
liuxdhs 发表于 Post on 5 day ago
tangwenke 发表于 2024-11-17 21:58
你好  为什么我算出来的应力是负值呢

deform的设置可能出现问题,方向上。
可参考以下mdp设置。
Pcoupl     = Berendsen
pcoupltype = anisotropic
tau_p = 2.0
ref_p = 100.0 100.0 0 0 0 0
compressibility = 0 0 4.5e-5 0 0 0

deform = -0.005 -0.005 0 0 0 0
deform-init-flow = yes
tangwenke 发表于 Post on 6 day ago
你好  为什么我算出来的应力是负值呢

202411172158171911..png (16.42 KB, 下载次数 Times of downloads: 0)

202411172158171911..png
liuxdhs 发表于 Post on 2024-3-29 22:44:31
foging 发表于 2024-3-10 19:39
你好,我也遇到这个情况。请问后来是怎么做出这个pxx和pyy曲线的呀,我输出三个方向的pressure都是不规则 ...

gmx energy,把所有数据导出来,画一下,然后你就知道哪个是你的曲线了。
liuxdhs 发表于 Post on 2024-3-29 22:43:33
herbrex 发表于 2024-3-6 13:09
您好 我现在正在寻找如何分析deform过程中应力的方法 能请教一下吗 有相关文献吗?

可以直接google,deform,molecular dynamics simulation, your materials。会有很多
liuxdhs 发表于 Post on 2024-3-29 22:42:29
牧生 发表于 2023-1-15 10:13
我有一个疑问,纵坐标达到了上百MPa,如此大的应力才发生形变,这个如何理解啊。

我觉得这是跟力场跟控压有关的。这是全原子模型,力场占据了非常重要的因素。
我最近做bead-spring model,改变bond coeff的k值,stree也会随之变化。
liuxdhs 发表于 Post on 2024-3-29 22:40:38
Sui1999 发表于 2023-8-9 10:02
您好,可以分享一下修改后的mdp文件吗,主要是compressibility和defrom那两行。

compressibility = 0 0 6.5e-5 6.5e-5 6.5e-5 6.5e-5
deform = -0.003 -0.003 0 0 0 0
foging 发表于 Post on 2024-3-10 19:39:01
liuxdhs 发表于 2022-10-29 15:17
x,y方向计算出来的Pxx,Pyy如下图。

你好,我也遇到这个情况。请问后来是怎么做出这个pxx和pyy曲线的呀,我输出三个方向的pressure都是不规则抖动
herbrex 发表于 Post on 2024-3-6 13:09:40
您好 我现在正在寻找如何分析deform过程中应力的方法 能请教一下吗 有相关文献吗?
Sui1999 发表于 Post on 2023-8-9 10:02:52
您好,可以分享一下修改后的mdp文件吗,主要是compressibility和defrom那两行。
牧生 发表于 Post on 2023-1-15 10:13:04
liuxdhs 发表于 2022-10-29 15:17
x,y方向计算出来的Pxx,Pyy如下图。

我有一个疑问,纵坐标达到了上百MPa,如此大的应力才发生形变,这个如何理解啊。
liuxdhs 发表于 Post on 2022-10-29 15:17:03
本帖最后由 liuxdhs 于 2022-10-29 18:10 编辑
DwyaneWan 发表于 2022-10-28 02:34
换成Berendsen试一下。 另外你给的是x和y 方向上的deformation, z方向不会有规整的curve的吧。

你好,感谢回复。这是压缩的过程,已经理解,非常感谢!

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-23 01:03 , Processed in 0.629756 second(s), 32 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list