刚体双摆分析无法收敛

wgm096350 · 发表于 2018-1-15 21:30:12

本帖最后由 wgm096350 于 2018-1-16 15:34 编辑

复摆受重力载荷作用，使用B31单元进行模拟，梁1的刚体约束参考点为rep1，梁2的刚体约束参考点为rep2，具体设置如图
分析监控里面有如下错误

这样简单的模型如法收敛，肯定是哪里设置有问题，求各位朋友指导一二！
相同的模型，只要将连接单元类型改为Beam就能正常计算，非常奇怪！
inp文件见附件

tommmy · 发表于 2018-1-29 21:34:53

本帖最后由 tommmy 于 2018-1-29 21:41 编辑

可以试试explicit，结果和全柔性体implicit的计算差不多（中间是末端位移总量-时间曲线）：

个人觉得在MBD方面，explicit表现要比implicit强一些，当然并不能说明ABAQUS在MBD方面很差，相反，理解了连接器的思路，你会发现ABAQUS在多体方面十分强大，就求解精度方面之于ADAMS完全有过之而无不及，它有各种各样的连接器供你使用，而且处理刚柔耦合不需要MNF，十分方便的。下面是完全刚体约束下implicit的计算结果:

下面是刚体条件下explicit与implicit的表现，几乎完全一致。

你的建模没问题，也就说不是ABAQUS MBD方法的问题，而是收敛性问题，多说一句哈，遇到问题多琢磨琢磨对寄几对大家都好的，加油

wgm096350 · 发表于 2018-1-16 12:42:45

麻烦大家帮忙看看！！！！！

wgm096350 · 发表于 2018-1-16 13:11:52

如果将两个梁都改为柔性体，也能成功计算收敛！！

wgm096350 · 发表于 2018-1-16 14:27:39

将梁1保持为刚体，梁2设置为柔性体，也能收敛！！！

断之桥 · 发表于 2018-1-16 14:30:50

给顶一下

wgm096350 · 发表于 2018-1-17 14:37:55

这个东西放在ADAMS里面秒算，我只能说ABAQUS多体分析真的cuo

断之桥 · 发表于 2018-1-18 19:11:18

帮忙顶一下

:lol:lol

popmuck2005 · 发表于 2018-1-19 09:29:19

自己不会用赖软件哈哈哈

wgm096350 · 发表于 2018-1-24 16:47:24

popmuck2005 发表于 2018-1-19 09:29
自己不会用赖软件哈哈哈

大佬，请指点一下呗！这个问题一直没解决，我现在都只能用ADAMS做项目了。

wgm096350 · 发表于 2018-1-27 15:42:15

有酬回复！能帮我解决问题的朋友必有重谢！

wgm096350 · 发表于 2018-1-28 15:54:38

popmuck2005 发表于 2018-1-19 09:29
自己不会用赖软件哈哈哈

朋友，可否指出我的问题？

popmuck2005 · 发表于 2018-1-30 10:40:50

既然是刚体动力学，那为啥不用adams和recurdyn呢？

wgm096350 · 发表于 2018-1-31 16:29:44

tommmy 发表于 2018-1-29 21:34
可以试试explicit，结果和全柔性体implicit的计算差不多（中间是末端位移总量-时间曲线）：

个人觉得在MBD ...

谢谢朋友！这两天仔细研究了你的答案，并自己建模研究了显式算法，有如下疑问：
1、请问你使用的模型是从我的inp里的导入的吗？我自己使用这个 inp 文件中的模型则无法计算，请问你是改了什么设置才能计算的呢？这个 inp 文件中单个摆的长度为 300mm，截面为 10mm 半径的圆截面，是无法计算的，但是我保持其他参数不变，把摆的长度改到 20mm，截面改为 0.1mm 半径的圆截面则能够计算。这是很让人费解的事情，请问该如何解释？
2、我使用 implicit 计算上述长度 20mm，截面 0.1mm 半径圆截面的摆时，采用自动步长，从监控器里边可以看到其采用的最小步长是 1e-5s，1s 时刻自由端节点合位移大小为 21.48mm。使用 explicit 算法计算时，我测试三种步长（由于全是刚体，explicit无法使用自动步长）：1e-5、1e-6、1e-7s（全部开双精度计算）。1s 时刻自由端节点合位移分别为 34.49、34.50、33.99mm，可见显式算法不同步长下结果差别不大，但与前面隐式算法的结果相差甚远。请问朋友你的显式算法步长是如何决定的呢？

wgm096350 · 发表于 2018-1-31 19:37:06

tommmy 发表于 2018-1-29 21:34
可以试试explicit，结果和全柔性体implicit的计算差不多（中间是末端位移总量-时间曲线）：

个人觉得在MBD ...

朋友，请问可否方便把你模型对应的 inp 文件发给我一下？

YHQ881892 · 发表于 2018-1-31 19:46:58

2797871439加我，我可以搞定　　隐式算得慢，显示算得快，结果基本一样，图为X向位移量

tommmy · 发表于 2018-1-31 20:33:00

wgm096350 发表于 2018-1-31 19:37
朋友，请问可否方便把你模型对应的 inp 文件发给我一下？

用的就是你的模型啊，你把求解参数改为固定时间步长，increment size调小点，就行了。注意允许的最大步数一定要够哦（≥t/increment size）,只是对于求解时间，做好思想准备。
相比之下，同样精度，explicit要快多了。

wgm096350 · 发表于 2018-2-1 08:44:51

popmuck2005 发表于 2018-1-30 10:40
既然是刚体动力学，那为啥不用adams和recurdyn呢？

进行刚体动力学分析只是我现在在做课题的一部分，我最终想要用脚本全部实现自动化做机构的优化分析。其中还会涉及到柔性体分析，还要与外部数据交换。adamas 与 recurdyn 参数化应该没有 abaqus 方便。
-----------------------其实主要是没那么多精力去学。

-----------------------------

wgm096350 · 发表于 2018-2-1 12:41:35

YHQ881892 发表于 2018-1-31 19:46
2797871439加我，我可以搞定　　隐式算得慢，显示算得快，结果基本一样，图为X向位移量 ...

谢谢！为什么用你的模型算得非常非常快呢？太神奇了！MPC 约束比 connector 好用一些吗？

wgm096350 · 发表于 2018-2-1 14:56:33

tommmy 发表于 2018-1-31 20:33
用的就是你的模型啊，你把求解参数改为固定时间步长，increment size调小点，就行了。注意允许的最大步数 ...

我把隐式分析步长设为1e-3、1e-4、1e-5计算，都在计算了几步之后终止，提示步长过大。请问你使用的固定步长是多大？可否上传inp文件？

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[动态分析] 刚体双摆分析无法收敛

本帖子中包含更多资源

本帖子中包含更多资源

本帖子中包含更多资源