PFC单向流固耦合——模拟颗粒落入流动的水中

PFC中流固耦合有三种方式:

1、单向流固耦合(one_way):也就是颗粒受流体作用,但是流体不受影响。

2、利用达西定律实现双向耦合

3、和第三方的算法或者流体软件进行耦合(比如OpenFOAM)

这里做一个单向耦合的小例子——模拟颗粒落入流动的水中。

由于当水比较多的时候,流速不太容易受到下落的颗粒影响,这里简化为单向耦合是合理的。

首先生成cfd网格和颗粒。这里的网格使用我之前帖子中生成方形网格的小程序生成节点和单元文件。


new
domain extent -3 3
wall generate box -2 2 -1 1 -0.5 2
wall delete walls range id 2[x_pos=0.5]
[height=1]
[box_chicun=0.5]
[rdMin=0.01]
[rdMax=0.03]ball generate radius [rdMin] [rdMax] number 1000 tries 2000000 range x [x_pos+rdMin] [x_pos+box_chicun-rdMin] ...
            y [-box_chicun*0.5+rdMin] [box_chicun*0.5-rdMin] z [height+rdMin] [height+box_chicun*2-rdMin]
           cmat default model linear method deform emod 100e6 kratio 1.5 property fric 0.5
ball attribute density 2.7e3 damp 0.5
set gravity 0 0 -9.8
configure cfd
cfd read nodes Node.dat
cfd read elements Elem.dat
element cfd attribute density 1000.0
element cfd attribute viscosity 1.5
define set_fluid_velocity
  loop foreach local ele element.cfd.list
    element.cfd.vel.x(ele) = -0.5
    element.cfd.vel.y(ele) = 0.0
    element.cfd.vel.z(ele) = 0.0
  end_loop
end
@set_fluid_velocitysave sample


之后进行一个下落:

keliweiyi.gif

这里颗粒的颜色代表粒径,可以看到粒径大的颗粒先落入河底,粒径小的颗粒随水流飘得比较远,这里也是符合常理的。

这里显示一下水流对颗粒的拖拽力:

cfdForce.gif

可以看到大颗粒所受的力比较大,而且拖拽力也是在变化的。

这里显示颗粒粒径在x向的分布:

lijingfenbu.gif

可以比较数值化的看出不同粒径在水流中的分离。

PFC2D/3DCFD

PFC单向流固耦合——模拟颗粒落入流动的水中的评论3条

  • 仿真小白_7656
    0
    大神,我用压力边界分别在cfd网格左右两侧施加不同压力形成压力差,为啥小球不动?请指教。@lobby
  • Abel_7172
    0
    老师这个网格的两个文件的参数是多少?

PFC单向流固耦合——模拟颗粒落入流动的水中的相关案例教程

0 引言 对于含集料混凝土,集料强度相对于水泥强度高很多,所以一般不会考虑到集料的可破坏性。所以大部分情况下,使用clump来模拟集料混凝土也就够了。但是在一些特别的受集中力的工况下,集料的破坏也不少见,这种情况就需要去考虑集料的破坏了。 本文基于已有的研究成果,提出一种以cluster为集料颗粒的混凝土成样方式,并用单轴实验来进行测试。 1 成样 对于混凝土,集料一般不是很多,这里在空间中生成松
0 引言 最近刚好刷到《战狼1》的短视频,里面有一个经典的子*弹穿墙镜头,想着或许可以使用离散元的方法进行模拟。对于弹道分析、穿墙速度损失等方面应该会有一定的指导意义。案例中使用到了很多技术方法,比如子*弹的形状是用圆锥二次化后生成的,墙体的边界处理等。 1 成样 破碎方面一般只在被击中的砖周边,所以我们只建了三行三列的砖墙,其余的边界用有限元的zone模拟。砖头采用240*115*90mm尺寸,
前言 最近熟悉6.0最大的感受就是,之前的连续体和离散体混合模型都可以使用6.0了,这里给大家介绍一下在PFC5.0版本实现起来特别费力并且效果一般的柔性三轴实验。 在6.0中实现柔性三轴方便了不是一星半点,而且解决了很多颗粒膜存在的问题。有点实话就是: 耦合的方法使得之前的颗粒膜方法成为了笑话。 我个人用颗粒膜方法,存在的一个最大的问题就是,颗粒膜会出现大的变形,因为膜的弹性模量比较小,导致在加
1 引言 离散元作为一款基于颗粒力学的软件,得益于胶结模型的开发应用,使其在模拟胶结材料上展现出了令人惊讶的效果。对于散体材料,使用ball模拟可以很好的模拟颗粒在变形过程中的滑动、跃迁。但是美中不足的是,ball作为刚性体,不能体现出颗粒的破碎效果。 目前主流的是有两种方法,一种是监测ball的应力状态,当其超过强度阈值时,用若干小颗粒来代替大颗粒完成破碎效果。这种方法比较好的是计算效率高,且计
前言: 柔性双轴和三轴的代码一年半以前就已经写好框架了,但有一些Bug,一直没有时间去整理。最近计算力有点空闲,刚好今天有点酒意,将柔性双轴的代码先整理了一下,基本上土的力学特性是可以反应了,剩下的就靠各位去继续完善了。 首先我们得搞清楚为什么要做柔性双轴或者三轴。真三轴其实是更加符合土单元概念的力学试验,但是在现实中真三轴的难度可以说是假三轴的百倍以上。这就导致了目前很多土力学实验都是假三轴。而
硕士
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    3 3