例子的来源是Abaqus CLE的官方教程,可是写的太粗线条,我还是搞了两天才做出了这个例子。其实就是个滚筒洗衣机带着洗衣机里的水一起转的问题。1. 分别为Eulerian domain和Lagrangian domain建立两个part建立Lagrangian domain的Part,类型设置为Discrete rigid,并设置Reference Point。
建立Eulerian domain的Part,类型设置为Eulerian,要注意Eulerian domain和Lagrangian domain要保证有重叠的部分,这是一种弱耦合,数据在两个区域间抛来抛去,所以网格要有重叠部分。这导致在Eulerian domain里有的部分是有材料的,有的地方是没有材料的。为了之后设置材料分布时候方便,要把part实现划出几个辅助的partition。黄色虚线是在划分partition时,为了指明Extrude/Sweep方向用到的辅助坐标轴。
2. 定义水的材料属性选择状态方程模型EOS中Us-Up,设置声速c0=1483m/s;密度为1000kg/m3;粘度为0.001kg/ms。并把截面属性赋给Eulerian domain。
3. 把两个Part组装起来
4. 新建一个Step-15. 为Eulerian domain和Lagrangian domain划分网格
6. 设置接触新建一个Contact Property,因为不是普通的面和面的接触,水中的任何的一个部分可能在流动区域里的任何一个地方和Lagrangian domain接触,设置Tangential Behavior为Rough,赋给水和洗衣机之间的关系。新建一个Interaction,把刚才的Contact Property赋给它。
更重要的是设置接触的两个Surface。其中一个Surface是Lagrangian domain部分的内侧面,为Geometry类型,另一个Surface是Eulerian domain的全部网格,为Mesh类型。
7. 加载Load为全模型加载重力场。
为Eulerian domain施加一个Boundary Conidition,为了在前后两个方向限制住水,限制要覆盖所有水可能会流到的区域。水运动到这个区域就让它的z方向速度变为0,这样水就不会流出洗衣机。
下面为Lagrangian domain施加Boundary Condition,让洗衣机转动,在1s内转动3个弧度,分别在Initial和Step-1设置BC-2。
最后设置在Eulerian domain的材料填充,虽然已经把水的材料的截面属性赋给了整个Eulerian domain区域,但是默认情况下所有Eulerian domain都是Void,即没有材料(下图中材料为0,Void为1),这也是在新建Eulerian截面属性的时候,出现这两个材料的原因。
为Eulerian domain的Part中初始存在水的部分设置Predefined Field材料。为初始存在水的partition新建一个Set,在Part的这个set中设置材料为1,Void为0,其他部分不用设置,依然是默认的材料为0,Void为1,即没有材料填充。
8. 设置Field Output Request尤其要选中Volume/Thickness/Coordinates中的EVF,element volume fraction,为了观察水的流动。
9. 提交Job计算10. 后处理在Result-Field Output设置显示EVF_VOID看液体的流动。
红色部分void为1,没有水,蓝色部分void为0,有水。
为了让显示的更加清晰,设置不显示没有水的Eulerian domain区域。
将阈值改为100%,Spectrums改为Blue to red。
设置ODB Display Options和Options-Common,让结果显示为理想状态。最后将动画结果存储为avi格式。
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删