用Abaqus来玩玩slinky下楼梯和“诡异”的自由落体。 案例1. slinky下楼梯仿真 它的轴向刚度 =0.002829095N/mm 一般的弹簧轴向刚度是多大呢?1~200N/mm !!! 相比之下,slinky几乎没有轴向刚度可言! 正因为这样,slinky第一阶梯走完之后,才有可能由收尾惯性段抵抗很小的轴向刚度,甩向下一阶梯,并重复上一状态。 这是一个常规的动力学分析,材料属性是普通的塑料线弹性本构,模型的重点还是在边界条件的处理上,建模时以下两点需要注意: 案例2. “诡异”的自由落体仿真 像上图那样丢过slinky的人可能会发现,它的底端,会神奇地表现出“反重力”的特点,是不是和你头脑中的自由落体的感觉不太一样?
为什么会这样呢? 让我们用Abaqus来一探究竟! 首先,要明白一个概念:初始应力状态。一般情况下我们丢的东西初始状态都很放松,所以都表现正常。这个slinky呢,在吊起来被释放之前,由于重力,体内储存了大量的应变能,也就是通常说的弹性势能,那么,约束解除之后,应变能去哪了? 先看看力学过程,这个仿真可以分两个过程,第一个过程是求解弹簧在重力作用下变形的静力学过程;第二个过程是求解变形后的弹簧在重力作用下自由落体的动力学过程。而仿真建模的关键是在静力学分析与动力学分析之间传递数据,共享应力状态:Standard to Explicit(用Abaqus做过冲压成型-回弹分析的肯定对这个再熟悉不过了,过程和这里刚好相反:Explicit to Standard)。 建模要点: “诡异”的自由落体仿真结果 上图中,最左侧是从地面(惯性系)来看slinky的效果;最右侧是经过坐标变换,跟随一同下落的刚体物块来(非惯性系)看slinky的效果;中间曲线则是slinky的质心Y向坐标在惯性系内随时间变化关系。 通过仿真分析,我们会发现,slinky的质心的确是在做自由落体运动,而前面提到的应变能,则转化为弹簧从两端趋向质心加速的动能,从非惯性系看这个过程,和弹簧拉伸后自由释放的过程并没什么两样,而在惯性系的我们来看,应变能释放完之前,底端是停滞的。 老规矩,inp文件收点辛苦费~~~只有下楼梯这个,那两个简单的就不贴了。
|