0
  • 聊天消息
  • 系统消息
  • 评论与回复
登录后你可以
  • 下载海量资料
  • 学习在线课程
  • 观看技术视频
  • 写文章/发帖/加入社区
会员中心
创作中心

完善资料让更多小伙伴认识你,还能领取20积分哦,立即完善>

3天内不再提示

浅析基于Fluent的二自由度圆柱涡激振动数值模拟与应用

8XCt_sim_ol 来源:仿真秀App 作者:姜工 2022-11-29 16:02 次阅读

一、案例描述

本案例源自美国康奈尔大学经典圆柱水槽试验,在实际的工程应用中,管柱被认为是一个有弹性支撑的刚性结构,结合切片理论可将计算模型简化至二维,可以看作是一个质量弹簧阻尼系统,如图中所示

f67d0fd4-6cb0-11ed-8abf-dac502259ad0.png

圆柱涡激振动示意图

计算参数如下表所示

f6984cea-6cb0-11ed-8abf-dac502259ad0.png  

二、采用Designmodeler进行几何模型建立

Workbench中导入Fluid Flow

右击Geometry,选择designmodeler

详细几何建模过程不进行描述

f6b6569a-6cb0-11ed-8abf-dac502259ad0.png

流体计算域模型如图所示,模型使用的计算域为30D×20D矩形区域。入口边界距离圆柱中心10D处,出口边界距离圆柱中心20D处,上部边界和下部位于距圆柱体中心10D处,并被定义为对称边界。

其中圆柱位于入口下游10D处,位于横流方向的中心线(y=0)。阻塞率为圆柱直径与流场域宽度的比值,本文阻塞率为0.05,不超过0.05,故满足要求。采用“刚性运动区域+动网格区域+静止网格区域”多区域划分的方法来建立圆截面涡激振动流场区域,整体流场分为外流域部分和内流域部分。

注:为保证流场在内外流域交界地方一致性,选中内外流域,鼠标右击from a part,可保证后续画网格时具有网格一致性。

f6e05634-6cb0-11ed-8abf-dac502259ad0.png

f708cab0-6cb0-11ed-8abf-dac502259ad0.png  

三、采用Meshing进行网格划分

Workbench中导入Fluid Flow

双击Meshing,详细网格划分步骤在这里不介绍了

f72f0086-6cb0-11ed-8abf-dac502259ad0.png

整体流场在进行网格划分时,采用“刚性运动区域+动网格区域+静止网格区域”的 方法来建立圆截面涡激振动网格划分区域,如图中所示,具体划分区域如下:

f7429c18-6cb0-11ed-8abf-dac502259ad0.png  

(1) 动网格区域采用非结构化三角形网格进行网格划分,Fluent动态网格模型算法采用Smoothing+Remeshing的方式更新动网格,此区域为变形区,为圆柱的振动提供了一个运动的空间,其网格会因圆柱的运动而进行重新划分,因此网格的初始质量好但随着重新划分可能有所下降,但可通过参数设置使重新划分后的网格质量尽可能好。

此区域的大小由圆柱的振幅大小决定, 考虑到文献实验中的最大振幅为1.5D左右,因此动网格区域尺寸取为6D×6D,留有充足的空间允许圆柱进行大范围的的振动,能够很好地满足数值计算的需要(避免振幅超过容许空间从而挤压外边界形成负网格);

(2)刚性运动区域主要包括圆柱流固耦合边界与边界层区域,圆柱近壁面边界层网格采用膨胀层网格处理,使在主要流速下预估的Y+=1,边界层第一层高度可按经验公式计算,近壁面划分15层四边形结构化边界层网格,网格增加率设置为1.2,向外过渡至与动网格的交界处。

边界层区域网格尺寸较小、网格划分质量较好并具有与圆截面相同的运动速度,因此将随着圆截面的振动而振动,由此来保证圆截面附近的网格一直处于初始的高质量状态从而确保各时间步都具有足够的计算精度。

(3) 静止网格区域采用四边形结构化网格,整体网格尺寸网格尺寸保持不变,用以为模拟尾流中的涡街提供更好的计算精度。

四、Fluent 基本设置

步骤一:启动fluent

Workbench中导入Fluid Flow

双击Setup,不考虑运行计算的话,默认即可

f772c820-6cb0-11ed-8abf-dac502259ad0.png

步骤二:General 设置

点击左侧模型树节点 General ,出现 General 面板,如图所示。

f7925d5c-6cb0-11ed-8abf-dac502259ad0.png

对面板进行参数设置 Scale…

f7bcdbea-6cb0-11ed-8abf-dac502259ad0.png

点击scale…

将单位修改为与几何建模时尺寸单位保持一致

点击close

Check

点击 Check 按钮,屏幕输出网格信息如图所示,确保网格最小体积为正值

f7d2fdc6-6cb0-11ed-8abf-dac502259ad0.png

Time

f7ff34ea-6cb0-11ed-8abf-dac502259ad0.png

由于涡激振动属于周期性振动,Time选项设置为transient,进行瞬态研究 其余保持默认即可

步骤三 Models设置

选择模型树节点 Models ,在右侧 Models 列表项中双击,弹出湍流模型设置对话框。设置参数:

选择 Model 为 k-omega ( 2 eqn )

选择k-omega Model 为 SST

选择OK按钮关闭对话框

其他保持默认即可

f81e880e-6cb0-11ed-8abf-dac502259ad0.png

步骤四:Materials设置

选择模型树Materials,分为流体材料fluid和固体材料solid,双击材料空气air进行设置,弹出材料设置窗口

f845b3b6-6cb0-11ed-8abf-dac502259ad0.png

f85f407e-6cb0-11ed-8abf-dac502259ad0.png

点击Fluent Datebase…

下拉菜单至最后,找到材料水(water)

点击Copy

其余选择默认,点击close关闭

f879a770-6cb0-11ed-8abf-dac502259ad0.png

步骤五 Cell Zone Conditions

选择模型树Cell Zone Conditions,双击弹出区域设置窗口,如图所示

分别选择part-dongquyu和solid-part-jignzhiyu,将属性设置为流体

将material name 选择为 water-liquid

其余保持默认,点击close关闭

f8a31826-6cb0-11ed-8abf-dac502259ad0.png

步骤六:Boundary Conditions

点击模型树节点 Boundary Conditions ,参数面板如图所示。

f8c6a584-6cb0-11ed-8abf-dac502259ad0.png

1、Inlet 边界设置

鼠标双击 inlet 列表项。设置参数:

设置 Velocity Magnitude 为 0.09m/s

设置 Specification Method 为 Intensity and Hydraulic Diameter

设置 Turbulent Intensity 为5%

设置Hydraulic Diameter 为 0.0381

选择OK按钮关闭对话框

f8eab41a-6cb0-11ed-8abf-dac502259ad0.png

2 、 Outlet 设置

鼠标选中 outlet 项,设置 Type 为 Pressure-outlet,Gauge pressure设置为0,其他参数保持默认。

步骤七 :Dynamic Mesh

点击模型树节点 Dynamic Mesh ,参数面板如图所示。

f90408fc-6cb0-11ed-8abf-dac502259ad0.png

勾选Dynamic Mesh

勾选Smoothing和Remeshing

点击settings…,选择 spring

选择remeshing,点击default,点击ok

f926a52e-6cb0-11ed-8abf-dac502259ad0.png

f94ccf10-6cb0-11ed-8abf-dac502259ad0.png

导入外部UDF文件

f97ef6fc-6cb0-11ed-8abf-dac502259ad0.pngf99b8b96-6cb0-11ed-8abf-dac502259ad0.png

点击菜单栏User-Defined

点击Functions,选择第二个Compiled

选择需要的UDF文件

点击Build,点击加载Load

边界层网格分离

f9babc96-6cb0-11ed-8abf-dac502259ad0.png

f9e2fdfa-6cb0-11ed-8abf-dac502259ad0.png

fa29d2c0-6cb0-11ed-8abf-dac502259ad0.png

点击菜单栏Domain

选择Adapt,点击Manage…

选择Cell Registers,点击New

点击boundary

fa3f0096-6cb0-11ed-8abf-dac502259ad0.png

fa5bf520-6cb0-11ed-8abf-dac502259ad0.png

选择Boundary zone 中的fsi

设置number of cells 为相应边界层网格的层数

点击save/display

点击close关闭

fa7d157a-6cb0-11ed-8abf-dac502259ad0.png

点击Domain 中的 separate,选择cell…

选择boundary_0和part-dongquyu

点击separate

点击close关闭

边界运动或变形的指定

fa94ee20-6cb0-11ed-8abf-dac502259ad0.png

faa21c30-6cb0-11ed-8abf-dac502259ad0.png

fac12210-6cb0-11ed-8abf-dac502259ad0.png

点击动网格部分的dynamic mesh zones的create

Zone name 选择为fsi

Type 选择为刚体运动 Rigid body

UDF选择导入的fluid::libudf

将分离出的边界层网格part-dongquyu:013进行相同的操作

其余默认,点击close关闭

步骤八 Reference Values

选择模型树节点 Reference Values,如图所示

fae3f34e-6cb0-11ed-8abf-dac502259ad0.png

Compute from 选项选择为 inlet

Area 设置为0.0381

Density设置为998

其余保持默认

步骤九 Solution Methods

此模型保持默认参数

Step 11 :Solution Controls

保持默认参数

步骤十 Monitors

选择模型树节点 Monitors,如图所示,选择Report Files,鼠标右键,点击edit‘’

fafc87c4-6cb0-11ed-8abf-dac502259ad0.png

点击report files 中new…

点击new report file中的new,选择force report,点击lift…

fb216c42-6cb0-11ed-8abf-dac502259ad0.png

Name 可自己命名名字

点击lift force,选择fsi

勾选report file和report plot

点击OK,关闭窗口即可

步骤十一 Initialization

点击模型树中的 initialization,出现如图所示窗口

fb7cb75a-6cb0-11ed-8abf-dac502259ad0.png

选择 Standard Initialization

Compute from 选择 inlet

点击 initialize

步骤十二 Calculation activities

双击 solution animations,如图所示

fba19f98-6cb0-11ed-8abf-dac502259ad0.png

点击new object

选择contours…,如图所示

fbd8658c-6cb0-11ed-8abf-dac502259ad0.png

fbf03086-6cb0-11ed-8abf-dac502259ad0.png

fc16f1c6-6cb0-11ed-8abf-dac502259ad0.png

云图Contours of 选择 为velocity

点击save,出现初始速度云图

鼠标滚轮调整到合适窗口,选择contour2,点击use active

点击OK

步骤十三 Run calculation

点击左侧模型树节点Run calculation ,出现 Run calculation 面板,如图所示。

fc39536a-6cb0-11ed-8abf-dac502259ad0.png

将number of time steps 设置为需要时间步数

每一个时间步时长 time step size设置为0.01

点击calculate,进行计算

步骤十四 查看计算结果

点击模型树Results,可查看残差曲线、各种云图(速度、压力等),速度矢量图、动画等

fc653750-6cb0-11ed-8abf-dac502259ad0.png
fc810b10-6cb0-11ed-8abf-dac502259ad0.pngfc9a14d4-6cb0-11ed-8abf-dac502259ad0.png






审核编辑:刘清

声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
  • UDF
    UDF
    +关注

    关注

    0

    文章

    4

    浏览量

    6464
  • mesh模块
    +关注

    关注

    1

    文章

    5

    浏览量

    1399

原文标题:基于Fluent二自由度圆柱涡激振动数值模拟与应用研究

文章出处:【微信号:sim_ol,微信公众号:模拟在线】欢迎添加关注!文章转载请注明出处。

收藏 人收藏

    评论

    相关推荐

    DSP在六自由度电磁跟踪系统中的应用

    自由度电磁跟踪系统根据电磁感应的原理来计算出目标的六个自由度参数,从而唯一确定目标的位置和姿态。怎样由接收到的模拟信号计算得到目标的位置和姿态参数是影响整个系统性能的关键所在。在一些早期的产品中对感应天线接收信号的处理主要采用
    发表于 01-25 00:11 1296次阅读
    DSP在六<b class='flag-5'>自由度</b>电磁跟踪系统中的应用

    汽车驾驶模拟自由度运动系统

    汽车驾驶模拟自由度运动系统
    发表于 08-06 13:06

    请问如何模拟弹簧阻尼振动系统?

    我想模拟自由度阻尼振动,并采集振动速度信号,请问应该如何实现?
    发表于 04-29 00:25

    自由度两弹簧系统振动实验设计与分析

    设计了自由度两弹簧振子系统,并进行了系统理论模态分析与ANSYS仿真,通过激振实验对系统进行了频率响应特性和模态分析,研究了系统振动动态特性测试。研究结果表明,实
    发表于 02-28 17:07 21次下载

    基于单亲遗传算法自由度PID控制器设计

    针对自由度PID 控制器参数整定困难的问题,提出了一种基于单亲遗传算法的自由度参数整定方法。运用该算法分别对单自由度PID 控制器和
    发表于 12-26 13:56 23次下载

    SIMTEC HEXADRIVE 6自由度运动模拟平台

    SimTec仿真技术公司总部位于德国北部Braunschweig,是一家拥有超过20年历史的国际领先的运动仿真系统制造公司。 SimTec专门电液伺服液多自由度运动模拟系统及运动数据采集系统的的软硬
    发表于 04-07 16:34 22次下载

    自由度内模控制方法

    针对典型的工业过程控制对象,提出了一种自由度内模控制方法,此方法可以使系统同时具有良好的目标值跟踪特性、干扰抑制特性和鲁棒性,改进了常规内模控制的不足,而且控制算法简
    发表于 05-23 15:12 27次下载
    <b class='flag-5'>二</b><b class='flag-5'>自由度</b>内模控制方法

    自由度转台的设计与实现

    基于MSP430F149单片机设计了一种三自由度转台系统。在该系统中,利用旋转编码器的信息反馈,通过对电机的闭环控制,使三自由度转台能够模拟飞机的飞行姿态。
    发表于 09-26 13:58 2537次阅读
    三<b class='flag-5'>自由度</b>转台的设计与实现

    自由度行波型超声波电机定子的优化研究

    自由度行波型超声波电机定子的优化研究_傅平
    发表于 01-04 16:45 1次下载

    基于自由度模型驱动PID的CFB锅炉床温控制_吴丹丹

    基于自由度模型驱动PID的CFB锅炉床温控制_吴丹丹
    发表于 01-13 21:36 0次下载

    自由度汽车油气悬架系统模型

    建立了自由度汽车油气悬架系统模型,分析了油气悬架和车身模型的非线性特性。运用AMEsim仿真软件,通过模型对不同的振频率的时域和频谱响应分析,验证了油气悬架非线性振动模型的正确性。
    发表于 01-16 16:49 2次下载
    <b class='flag-5'>二</b><b class='flag-5'>自由度</b>汽车油气悬架系统模型

    风电机组如何避免振动

    假若塔筒的自振频率与漩涡的发放频率相接近,就会使塔筒发生共振破坏,因此这种振动是极其有害的,而风机设备高昂的造价使其无法承受一次这种类型的破坏,因此,最大限度地避免被
    的头像 发表于 06-09 11:20 1.3w次阅读

    浅谈工业机器人自由度

    工业机器人自由度一般分为4个自由度,或者6个自由度,关节机器人主要是模仿人的手臂进行设计的。
    发表于 03-16 12:20 2828次阅读

    如何制作一个2自由度并联仿生腿?

    本文实例将实现2自由度并联仿生腿模组运动起来,模拟实现狗腿行走的动作。
    的头像 发表于 06-18 10:35 561次阅读
    如何制作一个2<b class='flag-5'>自由度</b>并联仿生腿?

    基于Simulink的汽车2自由度模型

    本期采用Simulink建立汽车2自由度动力学模型。为了便于掌握汽车稳定性的基本特性,对汽车简化为线性自由度的汽车模型,忽略转向系统的影响,直接一前轮转角作为输入;忽略悬架的作用,认为汽车车厢只作用于地面的平面运动。
    发表于 06-28 16:39 2380次阅读
    基于Simulink的汽车2<b class='flag-5'>自由度</b>模型