anays渗流.docVIP

水工流体力学问题的ANSYS模拟 ANSYS软件具有专门针对流体的计算模块，可以模拟理想流体的流动、稳定和非稳定渗流。本章主要介绍利用ANSYS进行简单的渗流分析和流体分析。 水工渗流场模拟 需要指出的是，在ANSYS中并不存在专门的渗流分析模块，但由于渗流场与温度场的有限元计算公式相同，因此可以利用ANSYS中的热分析模块进行分析，只需相应参数对应采取即可。 鉴于问题的复杂性，这里仅用一个重力坝的例子进行分析。 数值模拟对象 考虑混凝土大坝下水的渗流。假定坝下土壤的渗流系数K约为每天15米，即K＝15m/天，试确定该土壤的渗流速度分布。 本问题为一个稳态渗流问题，可以利用ANSYS的稳态热传导进行分析。 有限元模型 本问题中，考虑混凝土坝体不透水，渗流主要发生在坝基部分，因此，模型仅涉及坝基，具体建模过程如下： 1．模型的定义 启动ANSYS，设置好文件夹及文件名； 2．设置分析类型 点击Preferences，在弹出的对话框中选择Thermal，选择热分析，点击OK； 随后进入前处理模块：Main MenuPreProcessor； 3．设置单元类型 对于二维渗流（热）分析，采用Plane55单元。 Element TypeAdd/Edit/Delete…，弹出对话框： 4．定义材料参数 稳态渗流分析中，主要设置材料的渗透系数（本问题中，坝基材料的渗透系数为K＝15m/天，其对应热分析中的热传导系数）。设置方法如下： Material PropsMaterial Models…，弹出对话框中： 5．建立数学模型 本问题中，通过创建关键点，再由关键点直接生成二维（坝基）面。 在命令行中，分别输入关键点命令，生成关键点： K,1,0,0；K,2,5,0；K,3,9,0；K,4,16,0；K,5,16,4；K,6,16,5；K,7,9,5；K,8,9,4；K,9,5,4；K,10,5,5；K,11,0,5；K,12,0,4 连接关键点生成面 ModellingCreatAreasArbitaryThrouth KPs…，弹出点选择对话框： 顺时针选择关键点1，2，9，12，点击OK后生成一个面，同理连接（2，3，8，9）、（3，4，5，8）、（8，5，6，7）、（12、9，10，11），生成其它四个面。 6．划分网格 指定面的属性 这里，仅有一种材料和一种单元，采用默认值。 MeshingMesh AttributeAll Areas…，弹出对话框中，默认设置，点击OK。 指定线划分密度 选定所有的线，指定单元长度为1。 MeshingMesh Tool 划分网格 MeshingMesh Tool，在弹出的对话框中： 稳定渗流计算 首先进入求解模块：Main MenuSolution。 1．指定渗流边界条件 对于渗流而言，存在多种边界条件，包括流量、水头、流速等。本问题中为水头边界条件，及上游水头为10m，下游水头为0.5m（ANSYS中对应为温度边界条件），对于不指定的边界默认为不透水边界（本问题中假定模型左右两侧、底部及坝底为不透水）。具体做法为： 点击Define LoadsApplyThermalTemperatureOn Nodes，弹出点选择对话框： 指定上游边界 如上图，框选上游河床底节点，点击OK后，弹出下面对话框，Lab2框中选择TEMP，在VALUE中输入10，点击Apply，完成上游边界指定。 指定下游边界 如下图，框选下游河床底节点，点击OK后，弹出下面对话框，Lab2框中选择TEMP，在VALUE中输入0.5，点击Apply，完成上游边界指定。 2．求解 点击SolveCurrent Ls，完成求解。 3．计算结果查看 进入后处理模块：Main MenuGeneral PostProc。 查看水头等值线分布 Plot ResultsContour PlotNodal Solu…，弹出对话框： 查看渗流梯度等值线分布 Plot ResultsContour PlotNodal Solu…，弹出对话框： 查看渗流速度矢量分布 Plot ResultsVector PlotPredefined…，弹出对话框： 查看渗流速度沿路径A-B的分布 Path OperationsDefine PathBy Nodes…，弹出节点选择框，依次点选A到B及AB之间的节点，从而创建名为AB的路径PATH。 映射数据到路径AB上：Path OperationsMap Onto Path… 绘制AB路径上的渗流速度变化：Path OperationsPlot Path ItemOn Graph… 点击或保存文件。 理想流体的流动 在管道（隧洞或其它类型管道）内流动的液体可以分成两种不同的运动形式：层流和紊流。根据著