24小时热门版块排行榜    

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

043114076

木虫 (正式写手)

[交流] 【求助】dl_poly 如何固定键角 已有3人参与

比如A-B-C三个原子,我想固定键角,尝试将AB键长和AC键长固定来达到目的,结果报错:
DL_POLY terminated due to error   105
error - shake algorithm failed to converge
不知有没有其他方法来有效的固定键角?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

★ ★
ghcacj(金币+2):谢谢你把解决方案放上来,如果自己解决了,发上来给大家分享,这种方式很好。 2010-05-19 19:58:01
基本上已经解决了。我使用了两种方法
1.定义A-B-C角度,设置力常数足够的大。这样角度变化就可以忽略。但必须把时间步设短,否则也会报错。原因是力常数大了,振动频率也会增大。
2.其实控制AB和AC键长(AC之间距离)也是可以的,把时间步设定足够小就行。然后逐渐改大。
3楼2010-05-19 18:02:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

老虎大王

木虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+2):谢谢 2010-05-19 11:36:51
也许,需要把三个原子弄成一个rigid body,但又怕和别的设置相矛盾。
2楼2010-05-19 11:19:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongyesies

金虫 (正式写手)

红叶

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+1):谢谢 2010-06-01 16:04:12
dlpoly FILED里可以用RIGID和constraints来达到这个目的。
4楼2010-06-01 15:30:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

★ ★
ghcacj(金币+2):谢谢 2010-06-01 20:26:10
引用回帖:
Originally posted by hongyesies at 2010-06-01 15:30:05:
dlpoly FILED里可以用RIGID和constraints来达到这个目的。

rigid不行,因为我还有二面角
constraints只能固定键长吧?
5楼2010-06-01 18:36:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见