《飞轮有限元,坎贝尔图》由会员分享，可在线阅读，更多相关《飞轮有限元,坎贝尔图（21页珍藏版）》请在金锄头文库上搜索。
1、报 告项目背景：某公司计划设计一个飞轮，它的工作转速是0-12000rpm。飞轮由不锈钢制成，其中杨氏模量为140GPa，泊松比为0.28，密度为4400kg/m3，屈服应力为850MPa。飞轮在滚动轴承支座上，底部的垂直位移约束使飞轮绕中心轴旋转。问题：1、飞轮在惯性载荷作用下的最大应力水平。 2、是否在工作转速下发生或接近共振。建模：有限元分析首先需要建立模型。在此,首先将“fiywheeliges2016.igs“导入ansys，然后建立有限元模型。具体步骤如下：a. 导入IGES文件。用file importIGES导入后的文件如下图：b. 用PreprocessorModelingD
2、eleteLine and Below删除模型外多余的线。c. Use PreprocessorModelingOperateBooleansDivideLine by Line to divide lines and create additional keypoints needed for creating the flywheel radial cross-section.d. 模型被分开，用PreprocessorModelingDeleteLine and Below删除图形中多余的线。e. 用PreprocessorModelingcreateAreasArbitraryBy L
3、ine创建模型。分析：1、 用二维模型进行静力分析。a. 定义单元属性用PreprocessorElement TypeAdd/Edit/Delete定义单元类型，选择PLANE182单元。并在选项里选择axisymmetric。定义材料属性PreprocessorMaterial PropsMaterial Models。材料为前面要求的不锈钢，注意：此处注意单位换算！因为模型的单位是MM，需要将问题中的材料属性单位均换算成毫米制。b. 划分网格用PreprocessorMeshingMesh tool设定网格尺寸，完成二维模型的网格划分。此处不考虑网格收敛性问题，属于初次划分分析。划分之后
4、的模型如下图：c. 施加载荷按照设计要求，用SolutionApplyStructuralDefine LoadsDisplacement对飞轮两端轴承处施加固定约束。将最大转速12000rpm定义为载荷施加在模型上。在施加转速约束时，对于只建立了一个截面的模型，首先要用SolutionApplyStructuralDefineLoadsDisplacementSymmetry B.C.选择转轴作为位移的对称约束。用SolutionApplyStructuralDefine LoadsInertiaAngular VelocGlobal施加绕Y轴的转速d. 求解用SolutionAnalysi
5、s TypeNew Analysis定义静力分析的基本类型。并求解。应力结果如下图：即模型的最大应力为295.014MPa，在材料的许用应力范围之内。即在最大转速下，结构不会被破坏。2、 用三维模型进行模态分析a. 在导入二维模型的基础上，用PreprocessorModelingOperateExtrudeAreaAbout Axis选择截面，OK，然后选择轴，点ok。通过选择截面绕Y轴旋转得到三维飞轮模型。生成模型如下图：b. 定义单元类型和材料属性与静力分析一样选择单元类型为SOLID186，此处不需要选择axisymmetric。然后定义材料属性，同静力分析。c. 划分网格对三维实体进
6、行网格划分。同样是没有考虑网格收敛性的初划分，得到划分好网格之后的模型如下图：d. 施加约束和载荷按照设计要求，对飞轮与轴承接触的面施加固定约束，并将最大转速12000rpm定义为载荷施加在模型上。除了约束是施加在面上以外，其余方法与静力分析类似。e. 求解用SolutionAnalysis TypeNew Analysis选择模态分析。并用SolutionAnalysis Type Analysis Options设置分析6阶模态。各阶模态图如下：3、 转速分析由于陀螺效应，旋转结构的特征频率与其旋转速度相关。计算不痛旋转速度时的频率，可以得到各个模态频率随转动速度的变化曲线。在分析时，选择
7、0.955rpm和11997.736rpm转速计算特征频率，随载荷步的增加，转速增加，生成坎贝尔图需要在各个不同转速下的载荷步频率结果。运行后这两个转速下的频率结果如下图：输出坎贝尔图如下：可以看出频率与转速的结果如下：4、 网格收敛性分析网格的划分对结果会产生影响，这个影响的大小需要经过网格收敛性分析来看。a. 二维静力分析的网格收敛性在分析2中模型的节点和单元参数如下，定义为case1：Largest Number Number Number Defined Selected Nodes . . . . . . . . . . . 3575 1122 1122 Elements. . .
8、. . . . . . . 549 549 549在此情况下的最大应力为：295.014MPa接下来，更改网格划分时节点和单元数量如下，定义为case2：Largest Number Number Number Defined Selected Nodes . . . . . . . . . . . 2879 977 977 Elements. . . . . . . . . . 426 426 426重新进行应力分析，最大应力值为：298.301MPaCase1与Case2收敛性水平为：Case2-Case1Case1100%=298.301-295.014295.014100%=1.11%
9、再进行一组比对，定义为case3：Largest Number Number Number Defined Selected Nodes . . . . . . . . . . . 3920 1467 1467 Elements. . . . . . . . . . 861 861 861最大应力值为：308.218MPaCase3与Case2收敛性水平为： Case3-Case2Case2100%=308.218-298.301298.301100%=3.32%两个均小与5%，那么可以认为网格大小对应力分析的结果没有造成影响，结果可靠。b. 三维模态分析的网格收敛性在分析3中模型的节点和单元
10、参数如下，定义为case4： Largest Number Number Number Defined Selected Nodes . . . . . . . . . . . 45358 45358 45358 Elements. . . . . . . . . . 28122 28122 28122在此情况下的基频和振幅分别为：27.8411Hz和8.60145mm。接下来，更改网格划分时节点和单元数量如下，定义为case5：Largest Number Number Number Defined Selected Nodes . . . . . . . . . . . 60811 608
11、11 60811 Elements. . . . . . . . . . 38711 38711 38711重新进行模态分析，基频和振幅分别为：27.8162Hz和8.60156mm。Case4与Case5收敛性水平为：基频：Case5-Case4Case4100%=27.8162-27.841127.8411100%=0.1%振幅：Case5-Case4Case4100%=8.60156-8.601458.60145100%=0.001%再进行一组比对，定义为case6： Largest Number Number Number Defined Selected Nodes . . . .
12、. . . . . . . 30248 30248 30248 Elements. . . . . . . . . . 18132 18132 18132基频和振幅分别为：27.8998Hz和8.60155mm。Case5与Case6收敛性水平为：基频：Case6-Case5Case5100%=27.8998-27.816227.8162100%=0.3%振幅：Case6-Case5Case5100%=8.60155-8.601568.60156100%0两个均小与5%，那么可以认为网格大小对结果没有造成影响，结果可靠。建议：a. 从应力分析结果可以看出，飞轮的最大应力位于从内向外的第一个薄壁处。因此，在不影响其他部件运动和安装的前提下，适当的增加此处的厚度可以提高飞轮的强度。
