2025《基于SolidWorks的面条机挤面部分的流固耦合仿真分析案例》2500字.docxVIP

PAGE51

PAGE51

基于SolidWorks的面条机挤面部分的流固耦合仿真分析案例

目录

TOC\o1-3\h\u9402基于SolidWorks的面条机挤面部分的流固耦合仿真分析案例 1

48621.1流体区域与固体区域的模型构建 1

267001.2面团流体仿真 2

233611.2.1网格划分 2

287261.2.2物性条件与计算条件设置 3

216441.2.3流体计算结果 4

142981.3挤面轴固体力学仿真 6

120811.3.1网格划分与材料属性 6

119861.3.2固体计算结果 7

252651.4小结 8

1.1流体区域与固体区域的模型构建

将面条机SolidWorks三维模型保存为x_t格式，导入ANSYSWorkbench软件库中的Spaceclaim软件，对其进行切分，只保留其挤面部分固体域，如图1.1所示。

图1.1面条机挤面部分固体域模型

易知方孔为面团入口，面模上的圆孔为面条出口，面仓内表面和螺旋绞龙挤面轴外表面为面团边界。

基于面条机挤面部分固体域模型，采用布尔操作将流体域抽出，流体域几何模型如图1.2所示。

图1.2流体域几何模型

流体域几何模型可分为两部分，面团入口处为静止流体域，螺旋挤面处为旋转流体域。划分静止流体域和旋转流体域的目的是为了方便后续进行流体仿真工作。

1.2面团流体仿真

1.2.1网格划分

使用ANSYSWorkbench软件库中的Fluent软件对面条机挤面部分的面团流动进行仿真，使用Fluent软件中的Mesh模块对流体域进行网格划分。使用多面体网格（polyhedramesh）对流体域进行网格划分，多面体网格具有网格质量好、收敛性好、计算速度快、集合适应性好等优点。流体域网格划分结果如图1.3所示。

图1.3流体域网格划分结果

对螺旋绞龙挤面轴外表面进行网格加密，在外表面划分边界层网格共5层，网格数量为74万，网格最大扭斜度（skewness）为0.69。

1.2.2物性条件与计算条件设置

面粉与水经混合、搅拌之后形成面团。面团组成成分较为复杂，其主要成分为水、淀粉及蛋白质，可将其近似看成高分子聚合物的浓溶液，其流体性质属于粘流态非牛顿流体REF_Re\r\h[26]，是一种幂律流体，且是幂律流体中的假塑性流体，其流变特性常用如下方程式描述REF_Re\r\h[27]。

τ=ky

式中，τ——剪切应力，Pa；

y——剪切速率，s?1

k——稠度系数，或称为幂律系数，Pa·sn；

n——流性指数，或称为幂律指数。

从式（4-1）中我们可以看出，面团的黏度随剪切速率的增加而降低。这是因为面团中含有的高分子胶体粒子在静止时，粒子中的链状高分子相互缠结，黏度较大；而当流速较高时，流层之间的剪应力使得散乱的链状高分子滚动旋转并收缩成团，相互之间的缠结减少，黏度降低。使用Fluent进行流体仿真时，对于面团这种假塑性流体，大多采用carreau黏性模型REF_Re\r\h[28]计算流体黏度，该模型方程式如下：

η=η

式中，n——幂指数；

λ——时间常数，s；

η0

η∞

查阅相关文献REF_Re\r\h[29]可知面团参数如表4-1所示。

表4-1面团参数

名称

符号

单位

数值

幂律指数

n

0.4

时间常数

λ

s

2.12

零剪切黏度

η

Pa·s

24350

无限剪切黏度

η

Pa·s

0

密度

ρ

kg/

850

在Fluent软件的Setup模块中设置物性条件与计算条件。该流体域几何模型由静止流体域和旋转流体域两部分组成，形状也较为复杂，为达到较精确的计算结果，使用双精度稳态求解器。面团流动速度为低速，使用层流模型求解面团流动情况。使用单参考系法（SRF）定义旋转流体域的转动，转动速度为48r/min。压力与速度的耦合方法使用couple方法，计算格式均使用二阶迎风格式。对旋转绞龙表面的力和进出口流量进行监测，计算至力和出入口流量稳定视为计算完成。

1.2.3流体计算结果

在Fluent软件中的Solution进行求解，并最终在Results中得到求解结果。

螺旋绞龙挤面轴表面和垂直于入口平面截面的面团流场压力分布如图1.4所示。

图1.4面团流场压力分布云图

可见流体压力从入口至出口呈递增规律，到达螺旋绞龙挤面轴叶片底部时压力最大，压力最大值为39.74Mpa。螺旋绞龙挤面轴叶片的前表面压力大于后表面，这