欢迎参与abaqus模拟牛顿摆(会有奖励)

北鹰南飞

刚看到http://forum.simwe.com/thread-993739-1-2.html出的这个题，感觉有点意思，比较适合大众口味，我们版也借此模型

让大家练练手，主要练习abaqus中的接触算法（算法任选），可以采用static,general;dynamic implicit ;dynamic explicit中的

任意算法，希望大家踊跃参与，好的计算模型，重赏！！

************************************************************************************************************

******模型尺寸和材质：球直径20mm，钢球，悬线长500mm，初始钢球偏置角度为30度.********

************************************************************************************************************

￥￥￥奖励规则：

所展示的模型必须与已经贴出的模型有所差异，不能雷同,否则不予奖励。
   
模型背景介绍【转载】：

牛顿摆是由法国物理学家伊丹·马略特（Edme Mariotte）最早于1676年提出的。当摆动最右侧的球并在回摆时碰撞紧密排列

的另外四个球，最左边的球将被弹出，并仅有最左边的球被弹出。当然此过程也是可逆的，当摆动最左侧的球撞击其它球时，

最右侧的球会被弹出。当最右侧的两个球同时摆动并撞击其他球时，最左侧的两个球会被弹出。同理相反方向同样可行，并

适用于更多的球，三个，四个，五个……

實驗原理 要瞭解由五顆大小、質量相同的鋼珠所組成的牛頓擺的擺動原理之前，先瞭解較簡單的兩顆大小、質量相同的鋼珠

的碰撞情形。當理想的剛體在碰撞時，必定遵守力學能守恆、動量守恆，我們先討論一維碰撞的情形：

zsq-w

非常有意思的话题，欢迎大家踊跃参与，待会儿我也来试一试。

王兆喜

肿么没有人啊？？？

sqy001

很有意思的题目，有空做做试试

响马0726

建议abaqus 出个 20.经典案例共析 模块，也不一定要怎么经典，平时制作或练习所用到的即可，供大家模仿、学习之用，（简单的说明+cae文件+inp文件）

北鹰南飞

5# 响马0726
好建议！

jinzhuce

5# 响马0726
这个建议不错！！希望班里的能手们，整理一下！！我们也好学习学习！！

北鹰南飞

7# jinzhuce
新手也可以试着做做，为减小难度，可以采用一些方法处理，如：
悬绳可以先用联接代替，也可以考虑用索单元或者非线性弹簧，为减少计算量和难度，也可以先做个二维的，可以先做两个球的，再三个球的，4个，5个。。。。尽量采用刚体球减少计算量。。。。等等。

zsq-w

7# jinzhuce  
新手也可以试着做做，为减小难度，可以采用一些方法处理，如：
悬绳可以先用联接代替，也可以考虑用索单元或者非线性弹簧，为减少计算量和难度，也可以先做个二维的，可以先做两个球的，再三个球的 ...
北鹰南飞 发表于 2011-7-17 21:36

赞同将军。  我刚才用显式动力算了下，计算量不小。光是第一个球的个体运动就有0.36 s左右。 明天再好好做做。


PS: 这是新手学习的一个好机会，新手做出来加分从重哦~

北鹰南飞

9# zsq-w
"光是第一个球的个体运动就有0.36 s左右。"
也可以用这种思路来做：对于没有接触的运动过程，用隐式比较快，如果小球运动到了接触状态，再切换为显式来算，这样可能会快点，不过这样做需要调试确定切换时间点。

北鹰南飞

9# zsq-w
突然发现，原来就是在aba版，老大们的回帖，我也不能编辑，，这个嘛~~~~ ~~~~

我本想把“PS: 这是新手学习的一个好机会，新手做出来加分从重哦~”突出的。

zsq-w

9# zsq-w  
"光是第一个球的个体运动就有0.36 s左右。"
也可以用这种思路来做：对于没有接触的运动过程，用隐式比较快，如果小球运动到了接触状态，再切换为显式来算，这样可能会快点，不过这样做需要调试确定切换 ...
北鹰南飞 发表于 2011-7-17 22:07

我也是这么想的，用个import吧

hoop247

新手比如我对于这种问题，果然很无助！不知道用有限元做这个，会不会小题大做，如果纯考验ABAQUS使用的话，还是有意义的。

tonnyw

Once this is done. Some capable guy can simulate this model.
http://sciencedemonstrations.fas ... b_pagecontent341734

zsq-w

Once this is done. Some capable guy can simulate this model.
http://sciencedemonstrations.fas.harvard.edu/icb/icb.do?keyword=k16940&pageid=icb.page80863&pageContentId=icb.pagecontent341734&state=maxi ...
tonnyw 发表于 2011-7-18 21:11

欢迎拿其他软件试一试啦，发上来大家一样可以探讨。

jackie.lee

我的做法是用transient dynamics来模拟：
~~~
1）如果用显式来算，很耗时；
2）球被处理成了刚体；
3）表面压印，为接触做准备，这样，网格可以划得很细；
4）计算单摆周期，确定大致的接触时间点；
5）未接触的行程用大的时间步，接触行程用很小的时间步来捕捉接触；
6）悬线用remote point来模拟；
7）接触里面设置impact constraint；
8）计算一次（半个周期），耗时仅5~10分钟；
9）这都是在ANSYS Workbench里面完成的。
~~~

敦诚

做了一个觉得挺好玩的，但是结果还需要改正。以前做过多体动力学的东西，现在都快忘光了，正好复习一下。

北鹰南飞

17# 敦诚
敦兄能否将模型模拟出一边的球下落，撞击，另一边的球弹起，中间球几乎不动的效果，如果做出那种效果就更完美了，这个也是本模型仿真的技术点。

seawideyp

等上班了找时间玩玩
蛮有意思

jackie.lee

我的做法是用transient dynamics来模拟：
~~~
1）如果用显式来算，很耗时；
2）球被处理成了刚体；
3）表面压印，为接触做准备，这样，网格可以划得很细；
4）计算单摆周期，确定大致的接触时间点；
5）未接触 ...
jackie.lee 发表于 2011-7-18 22:40

球越多，越难达到理想结果~

这是3个球的碰撞，不好意思，只能做到这样了 ~~~
总体感觉，用ANSYS做这种碰撞，尤其是刚体之间的碰撞，有点儿勉为其难了，以后有时间用LS-DYNA试试，应该比ANSYS好处理一点~~~