张增明_高速旋转轮盘模态分析

1、高速旋转轮盘模态分析在进行高速旋转机械的转子系统动力设计时，需要对转动部件进行模态分析，求解出其固有频率和相应的模态振型。通过合理的设计使其工作转速尽量远离转子系统的固有频率。而对于高速部件，工作时由于受到离心力的影响，其固有频率跟静止时相比会有一定的变化。为此，在进行模态分析时需要考虑离心力的影响。了解如何用ANSYS进行有预应力的结构的模态分析。一问题描述对某高速旋转轮盘进行考虑离心载荷引起的预应力的模态分析，求解出该轮盘的前10阶固有频率及其对应的模态振型。轮盘截面形状如图所示，该轮盘安装在某转轴上以12000转/分的速度高速旋转。相关参数为：弹性模量EX=2.1E5Mpa,泊松比PRX

2、Y=0.3,密度DENS=7.8E-9Tn/mm3。关键点坐标:1(-10,150,0)2(-10,140,0)3(-3,140,0)4(-4,55,0)5(-15,40,0)6(-15,10,0)RS=5二分析具体步骤1. 定义工作名、工作标题、过滤参数 定义工作名：Utilitymenu>File>Jobname：yuanpan 工作标题：Utilitymenu>File>ChangeTitle：yuanpan2. 选择单元类型本实验将选用六面体结构实体单元来分析，但在建模过程中需要使用四边形平面单元，所有需要定义两种单元类型：PLANE42和SOLID45，具体操

3、作如下：MainMenu'Preprocessor>ElementType>Add/Edit/Delete “StructuralSolid”“Quad4node42”Apply（添加PLANE42为1号单元） “StructuralSolid”“brick8node45”ok（添加六面体单元SOLID45为2号单元）3. 设置材料属性由于要进行的是考虑离心力引起的预应力作用下的轮盘的模态分析，材料的弹性模量EX和密度DENS必须定义。 定义材料的弹性模量EXMainMenu'Preprocessor>MaterialProps>MaterialMode

4、ls'Structural>Linear>Elastic'Isotropic弹性模量EX=2.1E5泊松比PRXY=0.3定义材料的密度DENSMainMenu'Preprocessor>MaterialProps>MaterialModels>densityDENS=7.8E-9ANStSMunh)rwc«UtilityM*nu(yuanpn)n-Lilt£lotPlot£tri«JorkPlanePytorsMacroMefluCtrlsyelp刮W勺空丄tlnnrftflypnDefinedEl

5、««entTvcesHlCLo“|SoGeneralPostprocTiaoHistPoatpraTopologicalOptROMToolDesignOptProbDesignRadiationOpt：Em,5»OI5io«4«：ia|円到®iajqq到Igj勻lAJadlAJQQl»!B他他到IJ£J"iCJ:九mAKSrsCmid(HEF7)AAN5YSMulliohvMuUUHtyMmu<yu«VMn-BD1.eforkPian*Pax<uwt®rrJjacroXr

6、QuCtrliUlp型凹£|創£回|X3闵ANSYSToolbar：H'v-':'iJ：*i4'HilAXSYSMa:nManuPraferancasPreprocessorsEleaentType* RelConstantsnMatarialProps MaterialLibrary TeaperatureUnits ElectrcmaHUnits FPWirErr1 ConvertALPx ChangoMatNub FailuroCriteria WritetoFite RnadfroaFtle«Section*$Modelin

7、sMeshinc CheckiryjClrlt(NumborfegCtrls ArchiveModel Coupling/CoqnzMultifioldSotUpaLoads* Phy»ics* P«<hOperationsSolutionGeneralPostprocTiaoHistPostproTopologicalOptioi««iiaADKWwMatwulMo<MIMaterialModelsDefinedMaterialMod»l>Available.AUmvImWOBIc*9*MMffWNumImy1LinearI

8、sotropicMaterialProperforMaterialNuab«r1T1TMperaturBfl|AddT»o(>»ratur*IDeleteT««p»ratur»j3誇7|«何刮®IIAJQIIJLIIIJlJ®Qa97l®a®llEH土创土.tlZJ4.实体建模对于本实例的有限元模型，首先需要建立轮盘的截面几何模型，然后对其进行网格划分,最后通过截面的有限元网格扫描出整个轮盘的有限元模型。具体的操作过程如下。创建关键点操作：MainMenu>Pre

9、processor>Modeling>Create>Keypoints>InActiveCSPOINTSTYPENUMANJUL5201510:58:37112231011 由关键点生成线的操作：MainMenu>Preprocessor>Modeling>Create>Lines>Lines>InActiveCoord 建立圆角：Preprocessor>Modeling>Create>Lines>Lines>LinesFilletPOINTSTYPENUMANJUL5201510:58:3711223

10、101149压空卩 生成面：MainMenu>Preprocessor>Modeling>Create>Areas'Arbitrary>ByLines5划分网络MainMenu'Preprocessor>Meshing>MeshTool 对全局进行设置。单击SizeControls（尺寸控制区）全局设置项（Global）的Set按钮，将弹出GlobalElementSizes单元尺寸全局设置对话框在对话框中输入ElementedgeLength（单元边长度）为6。 单击SizeControls（尺寸控制区）Lines（线设置项）的按钮，

11、将弹出ElementSizeonPickedLines（在所选线上定义单元尺寸）的拾取对话框。用鼠标左键在图形输出窗口中拾取圆角对应的线。单击ok按钮，将弹出ElementSizeonPickedLines（在所选线上定义单元尺寸）对话框,在对话框中输入No.ofelementdivisions（每条线将要分成的单元数）为1,设定圆角处对应的线只分一个单元。（由于是模态分析，只要能反应出需要知道的前几阶模态就行，而不需要知道具体的应力值，所以不需要对此处进行单元细化。） 对分网进行控制。在分网控制区的Mesh下拉框中选定分网类型为Area（面），Shape（网格形状）设置为Quad（四边形），

12、分网方式设置为Free（自由分网）。 对面进行分网。在MeshTool（分网工具）对话框中单击Mesh按钮，将弹出MeshAeras（对面划分网格）拾取对话框。从图形输出窗口中拾取创建的面，单击ok按钮。完成网格划分。1ELEMENTSANJUL5201511:09:216.出整个轮盘的有限元模型通过将面绕轴旋转成有限元实体模型的功能，将前面建立的轮盘截面有限元网格，围绕定义的旋转轴扫 定义旋转轴。可以通过定义旋转轴所在轴线上的两个关键点来，指定旋转轴的位置。MainMenu'Preprocessor'Modeling>Create>Keypoints>InA

13、ctiveCS生成两个关键点20、21。关键点20：X,Y,Z位置分别为-10,0，0关键点21：X,Y,Z位置分别为10，0，0 设置单元生成选项MainMenu'Preprocessor'Modeling'Operate>Extrude>ElemExtOpts，弹出ElementExtrusionOptions（单元挤出选项）对话框，在对话框中的Elementtypenumber（单元类型序号）下拉框中选择2号单元SOLID45。单元尺寸选项中的分割单元数（VAL1NO.ElemDivs）设置为18,即在挤压出的每个体上将沿周向被分成18份。拉伸比例为0

14、,保持等截面拉伸。将Cleararea（s）afterext（删除原始面）设置为Yes，在挤压的单元完成之后将删除原来的面以及其上的单元。单击ok按钮，完成对单元选项的设置。 绕轴旋转截面MainMenu'Preprocessor'Modeling>Operate>Extrude>Areas>AboutAxis,将弹出SweepAreasaboutAxis（绕轴扫描面）的拾取对话框。从图形输出窗口中选择创建好的平面网格，单击拾取对话框中的按钮。然后从图形窗口中选取定义旋转轴的关键点20,21,单击ok按钮,将弹出SweepAreasaboutAxis对话

15、框，在对话框中输入旋转角度为（Arclengthindegrees）360.No.ofvolumesegments（一周创建体的数目）为4,单击ok按钮。创建如图所示的整个盘的有限元模型。 观察创建的网格形式。ANSYSMultiphysicsUtilityMenu(yuanpan)£ileelectJUist£lotPlottrlsWorkPlanePax；ainetersMacroMenuCtirlsHelpANSYSToolgr|savedb|re$umdbqd工t|ANSYSMainMenu<refofencesreprocessor*aElemeirtTyp

16、e3Rea.1Constanfcs日alPropsaSectionsmModeling田Create曰Operate曰Ext.rudeE3ElemExtOptsejAreasAlongNorma1UtilityMenu>Plot>Element，图形窗口中将会显示出由平面网格扫掠而成的实体单元网格情况。存盘，SAVE_DB。至此，完成了创建轮盘有限元模型的所有工作。A1ongL»ines田L.lne田KeypointsKExtendL.ine田BooleansOScaleCoilcGeomI-terns田Move/Modify田Copy田Reflect田CheckGeom

17、田De1etemCyclieSectosr田CMS厅“c7.节点的坐标变换根据轮盘的工作情况其约束条件为盘心轴向和周向约束，这种约束条件在直角坐标系下无法定义，而柱坐标下可以非常方便地定义。根据ANSYS程序中坐标系的定义规则，需要将柱坐标系的Z轴和旋转轴重合，Y轴表示转角，X轴表示径向。ANSYS程序提供的全局柱坐标系不满足要求。通常可以有两种办法来解决这个问题：a. 将所建有限元模型进行旋转使其轴向和柱坐标Z轴方向一致。b. 重新建立一个柱坐标系使其的Z向和旋转轴一致。本实例米用第二种方法。具体操作过程如下：UtilityMenu>WorkPlane>OffsetWPbyInc

18、rements，弹出OffsetWP（工作平面偏移）菜单,拖动Degrees滑动条，将Degrees（旋转角度）值设置为90。单击更刃按钮，使工作平面绕Y轴正向旋转90度，单击ok按钮，将工作平面的WZ轴和总体坐标系的X轴方向重合。ANSYSMultiphysksUtilityMenu(yuanpan)3siMflOffsetWP2dJJX.Y.ZOthrtt3EJ£y-I£耐iri90AlW.2.Z<AnjtMufeJfeJfeJblJ&J»l&J0§JT171釘tej釦fejfcjANSYSToolbarSAVE.DBRESfM

19、DBQin|POWRGRPHAfile'electList£lotPlot£trlsWorkPlaneParametersMacroMefluCtrlsHelpANSYSMainMenureferencesreprocessorJElementTypeiRealConstantsiMaterialPropsJSectionsJModeling®CreateOperate曰ExtrudeElenExtOpts日AreasAlongNormalByXYZOffsetAlongLines田Lines田KeypointsExtendLine田Booleans田Sc

20、ale田CalcGeo<nIterns田Move/ModifyaCopy田Reflect田CheckGeoffl田Delete田CyclicSector田CMSz"二.GloMS0Y0Z>0厂DynameMode 在工作平面原点创建柱坐标系。UtilityMenu>WorkPlane>LocalCoordinateSystems>CreateLocalCS>AtWPOrigin，将弹出CreateLocalCSatWPOrigin（在工作平面原点创建本地坐标系）对话框，Refnumberofnewcoordsys（新坐标系的参考序号）缺省值为11，

21、一般就使用缺省值，也可自己设定。在Typeofcoordinatesystem（坐标系类型）下拉框中选取Cylindrical1（柱坐标系）,其它设置为缺省值。单击ok按钮。将完成要求的柱坐标系的创建，并且将新建坐标系定义为当前激活坐标系。 将所有节点移到当前柱坐标系中。MainMenu>Preprocessor'Modeling>Create>Nodes>RotateNodeCS>ToActiveCS将弹出RotateNodesintoCS菜单。在菜单中单击PickA11按钮，将所有的节点都移到当前激活柱坐标系下。8.进行静力分析由于对轮盘模态的分析需要

22、考虑离心力引起的应力对模态的影响，所以需要先对其进行静力分析，求解出离心力产生的应力，及其对刚度阵的影响，将结果写入数据库文件。有预应力的轮盘静力分析具体过程如下。 指定分析类型及分析选项a. MainMenu>Solution>NewAnalysis"Static"b. MainMenu>Solution>Sol'nControls，将弹出求解控制（SolutionControls）对话框。单击标签“Basic”，在CalculatePrestresseffects选项前打“V"。打开预应力选项。单击OK 定义边界条件和转速对于本

23、实例分析的轮盘，由安装条件知道其边界条件应该是，在轮盘盘心的节点轴向和周向固定，而径向自由。其离心载荷是由于高速旋转产生的，因此需要在分析时指定轮盘的旋转速度，具体操作过程如下：a. 将要加载模型放大。在图形窗口中选取盘心部分，对其进行放大，以便在加载时能够准确地选择所要加载的模型元素（本实例为节点）。b. 定义约束MainMenu'Solution>DefineLoads>Apply'Structural'Displacement>OnNodes，将弹出ApplyU,ROTonNodes（给节点施加约束）对话框，单击对话框中的Circle单选按钮，选

24、用圆形区域选择方式（因为需要约束的为盘心节点）。将光标移至盘心，按住鼠标左键并拖动光标，在图形窗口中将会出现一个选择用的圆形选择框随光标移动，当圆形框刚好将盘心所有的第一排节点选中，而没有选择第二排节点时，放开按钮。将会将所有盘心的节点选中,后面的约束也将施加在这些节点上。单击ok按钮，将弹出ApplyU,ROTonNodes（给节点施加约束）的对话框，在对话框的DOFstobeconstrained（被约束自由度）列表框中选择UY,UZ两个自由度，分别表示对周向和轴向施加约束。在Displacementvalue（位移约束值）文本框中输入0（缺省值为：0）。C.定义转速MainMenu

25、9;Solution>DefineLoads>Apply'Structural>Inertia>AngularVelocity>Global，将弹出Applyangularvelocity（角速度定义）对话框，在对话框中的GlobalCartesianX-comp（绕总体坐标系中X轴旋转角速度）文本框中输入1256.64,指定轮盘的旋转速度为1256.64弧度/秒。单击ok按钮，完成对旋转角速度的定义。存盘SAVE_DB。 进行求解MainMenu|Solution|CurrentLS9. 进行模态分析设定由于模态分析时位移约束条件和进行静力分析时一样，而

26、静力分析时所定义的转速在模态分析时将被忽略，因此进行模态分析时再不需要重新定义边条，可以直接进入求解菜单进行求解。 指定分析选项MainMenu>Solution'AnalysisType>NewAnalysisf“Modal"（模态分析） 分析选型设定MainMenu>Solution>AnalysisType>AnalysisOptions将弹出模态分析选项对话框，在对话框中，指定Modeextractionmethod（模态提取方法）为BlockLanczos（分块兰索斯法），并指定No.ofmodesextract（提取模态的阶数）为10

27、,将Expandmodeshapes（模态扩展）设置为“YES”，在No.ofmodestoexpand（模态扩展阶数）文本框中输入10。将Inclprestresseffects（预应力效应）设置为“YES”，这样在进行模态分析时ANSYS程序将会把前面静力分析中求解得到的离心力产生的应力对刚度的影响考虑进去。单击ok按钮,将会弹出BlockLanczosMethod（兰索斯法模态分析选项）对话框，在对话框中对模态提取的范围进行定义，一般定义的范围要求将所关心的所有频率都包含在指定的范围之内，这里指定的值为：StartFreq（开始频率）是0,EndFrequency（结束频率）是99999

28、999。单击ok按钮，完成对分析选项的设置。 进行求解MainMenuISolutionICurrentLS10.后处理列出固有频率一一MainMenu>GeneralPostproc>ResultsSummary，将列出轮盘的所有求解的固有频率，在文本框里列出了轮盘的前10阶固有频率，可以看出有些频率值相同，这是由于轮盘结构和边条都是对称的，会出现振型和频率相同但相位不同的情况。观察解得的模态本实例中由于设置了对模态进行扩展，所以对于求得的每一阶固有频率，程序同时都求解了其对应的模态振型，来反映在该固有频率时，轮盘的各节点的位移情况。可以利用通用后处理器方便地对其进行观察和分析，并且可以对各阶模态振型进行动画显示。具体操作过程如下：a. 选取菜单路径MainMenu>GeneralPostproc>ReadResults>FirstSet,选择轮盘第一阶模态。b. 选取菜单路径MainMenu>Gener