流体算例空气阻力系数计算dragforce

空气阻力系数计算本题是为了演示空气阻力系数的计算方法。首先导入一个汽车模型，如下图所示，此汽车模型是经过简化的。点击菜单ADINA-M>ImportParasolidModel，导入car_simple.x_t。导入汽车模型建立流场空间点击菜单ADINA-M>DefineBody，如下图所示建立一个立方体。建立流场空间点击菜单ADINA-M>BooleanOperator，如下图所示用第二个body减去第一个body，剪完之后剩下的部分就是真正的流场空间。注意，目前只有parasolid体才可以做布尔运算。在功能选择模块做如下设置，进行流场的稳态计算。进入流体模块，进行设置点击菜单Model>FlowAssumptions，在打开的窗口中做如下设置，表示三维模型、不考虑热、采用SA湍流模型。定义材料点击菜单Model>Materials>ManageMaterials，在打开的窗口中点击Spalart-AllmarasModel，定义一个SA湍流模型的材料。仅输入粘度和密度就可以，其它参数均采用默认值。定义约束点击菜单Model>UsualBoundaryConditions/Loads>ZeroValues，点击Add，定义两个约束条件，分别为YF和ZF施加约束点击菜单Model>UsualBoundaryConditions/Loads>ApplyZeroValues，给body2的1、2、4、5四个face施加约束条件：定义并施加荷载点击菜单Model>UsualBoundaryConditions/Loads>Apply，定义一个27.78m/s的速度载荷，约相当于100KM/h的车速。定义并施加特殊边界条件点击菜单Model>SpecialBoundaryConditions，如下图所示给body2的第7至65个face施加wall边界条件。可采用Auto的方式输入face编号。定义单元组点击菜单Meshing>ElementGroup，定义单元组指定网格大小点击菜单Meshing>MeshDensity>Face，给body2的face1至face6指定网格大小为1.5，face7至face65指定网格大小为0.08。可以采用Auto的方式输入面号。注意，外围的网格不重要，所以可以大些。划分网格点击菜单Meshing>CreateMesh>Body，给body2划分网格。定义求解控制点击菜单Control>SolutionProcess，选择使用FCBI-C单元，点击OuterIetration按钮，在打开的窗口中再点击AdvancedSettings按钮，进行如下设置。这里采用了FCBI-C算法，也可以采用FCBI，有关区别可参考手册。显示图形点击按钮，再点击按钮：保存命令流文件单击Save，保存为car.in文件。生成求解文件、并进行计算菜单：Solution>DataFile/Run后处理--切片进入后处理模块，打开car.por文件，点击图标，做一个切片，以便看到模型内部。画压力云图点击图标，显示模型外边界，点击，画出压力云图。定义一个ModelPoint菜单：Definition>ModelPoint（Special）>ReactionSum，定义一个ModelPoint菜单：List>ValueList>ModelPoint，列出x方向的反力，也就是空气阻力。求出阻力计算空气阻力系数可以得到空气阻力为388N，根据公式D=1/2ρV^2×S×CD，可计算出空气阻力系数。公式中：D：空气