| 查看: 3368 | 回复: 6 | ||
leoluode银虫 (小有名气)
|
[求助]
Lammps拉伸应变计算 已有1人参与
|
|
请问大家，如果输入文件中timestep=0.001(ps)，拉伸变形命令 fix 1 all deform 100 x erate 0.0001 units box remap x , run 40000步，大家是怎么计算最终的应变的啊？ 按目前我的看法是，工程应变速率为0.0001/ps,总共运行40000步，每步的时间为0.001ps,那么总应变为0.001ps*40000*0.0001/ps=0.0004=0.4%！但是感觉核deform后面的100没有啥关系啊，换成1000，总应变依然是0.4%! 所以我想问问大家对deform后面的100是如何理解的？ |
回复此楼» 猜你喜欢
国自然申请面上模板最新2026版出了吗？
已经有10人回复
-
推荐一本书
已经有12人回复
-
基金申报
已经有4人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有4人回复
-
溴的反应液脱色
已经有6人回复
-
纳米粒子粒径的测量
已经有7人回复
-
常年博士招收(双一流，工科)
已经有4人回复
-
参与限项
已经有5人回复
-
有没有人能给点建议
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps碳纳米管拉伸in文件报错求助 急
已经有1人回复
- lammps计算弹性模量应力应变
已经有4人回复
- 为什么lammps中的应变率这么大？
已经有4人回复
- lammps做拉伸试验，求大神进来指导下（烈日炎炎，有劳了`(*∩_∩*)′）
已经有14人回复
- LAMMPS进行拉伸模拟的一些问题
已经有4人回复
- 关于lammps计算应力——应变的几个问题
已经有2人回复
- lammps做拉伸和压缩的应力问题
已经有9人回复
- lammps里，关于应力应变的问题
已经有0人回复
- LAMMPS模拟宏观应力可以实现吗？
已经有2人回复
- lammps实验拉伸Al纳米线的in文件的错误提示好久找不出？
已经有8人回复
- 关于lammps fix npt/deform做拉伸时温度的变化
已经有6人回复
- lammps 拉伸应力应变曲线
已经有4人回复
- fix deform 指令拉伸求教
已经有10人回复
2楼2016-07-25 11:42:49
leoluode
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 346.5
- 帖子: 62
- 在线: 53.4小时
- 虫号: 3487208
- 注册: 2014-10-20
- 性别: GG
- 专业: 金属结构材料
|
哈哈，多谢一直的回答啊！ 关于您说的，我还是有点疑问！ （1）deform后面的N是每N 步变形一次，那么这个N设为1，就是每1步变形一次、总共运行40000次，那么模型就变形4万次，总变形时间为40000*0.001ps=40ps，应变速率为0.0001/ps，则总变形为0.004； 如果N后面设为100，就是每100步变形一次，总共运行40000次，实际上变形只有400次，则总共变形时间为400*0.001ps=0.4ps，应变速率为0.0001/ps，则总变形为0.00004。 所以N是会影响总变形的； （2）deform后面的N为1或100的区别在于，N为1时，是每运行1步，就变形1次；而N为100时，就是第一步变形了，第2-100步没有变形，第101步变形，第102-200步没有变形。第201步变形。。。。。。。 N为1或100区别就是变形后的弛豫时间步相同 不知道上面的理解您认为有问题没？ |
3楼2016-07-25 15:16:16
leoluode
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 346.5
- 帖子: 62
- 在线: 53.4小时
- 虫号: 3487208
- 注册: 2014-10-20
- 性别: GG
- 专业: 金属结构材料
4楼2016-07-25 15:18:41
5楼2016-07-25 16:20:21
leoluode
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 346.5
- 帖子: 62
- 在线: 53.4小时
- 虫号: 3487208
- 注册: 2014-10-20
- 性别: GG
- 专业: 金属结构材料
|
追梦2016 你好， 通过近一段时间的实际模拟测试，我发现关于拉伸应变计算结果与你的解释是一致的！ 对于timestep=0.001(ps)，拉伸变形命令 fix 1 all deform N x erate 0.0001 units box remap x , run 300000步， 总应变是由erate=0.0001/ps、timestep=0.001ps、run 300000步三者决定，总应变为0.03，与deform后的N值无关！ （1）N=1时，意味着每运行一步，就变形一次，总共运行300000步要达到变形量0.03，也就是说每次运行的应变为10^-7，即每运行一步，应变增加10^-7，运行到10000步时，总应变为10^-3;则在step=1~10000的过程中，总应变依次为10^-7,2*10^-7,3*10^-7.........,10^-3 (2)N =10000时，意味着每运行10000步，就变形一次，总共运行300000步要达到变形量0.03，亦即在每个10000th变形一次，且应变为10^-3;;则在step=1~10000的过程中，总应变依次为0,0 0,...，0,10^-3，step达到10000时，应变一下子就由0增加大10^-3。 上述分析说明，实际上N=1和N=10000时，模型的实际变形速率确实是不一致的，N=1是每一步变形，变形速率为10^-4/ps; N=10000时在10000步之前都保持0应变，而在每个10000 th时，直接变形达到规定值,变形速率为1/ps。N=10000时的实际变形速率是N=1的10000倍。 无论如何，erate后规定的应变速率是一样的（10^-4/ps），但是N改变时实际变形速率也改变了，只有当N=1时，erate规定的变形速率才与实际变形速率是一样的！那么变形速率我们在写论文时到底说采用哪一个呢？ 还望“追梦2016”同学能继续交流啊，能够指出我理解上的错误 |
6楼2016-07-27 15:30:22
|
好的，首先感谢对我的信任。 我的理解是，咱们所说的变形速率（或者说论文时体现的变形速率）应该是erate所规定的那个值。可能5楼解释的不太准确，这里N=1或=10000，只是告诉计算机，你多长时间（或者多少步）计算一次，但计算的依据是根据你说的erate、timestep和run共同决定的。只是说，如果N越小，我们就能得到越多的值，规律可能就越明显。举个例子，我们要找变形的中间规律，如果我们设定N=10000的话，那么用30000步得到0.03变形量我们得到30个数，可以画出变形规律，但如果N=1的话，那将得到30000个值，这样我们得到的规律就更准确。 所以，我认为，N值只是给计算机看的，他对我们研究的结论、属性没有什么关系。 当然，仅为个人见解，可能会有错误，希望多加讨论。 |
7楼2016-07-27 16:07:00