均匀直杆的子空间法模态分析1.模态分析的定义及其应用模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型),即结构的固有频率和振型,它们是承受动态载荷结构设计中的重要参数。同时,也可以作为其它动力学分析问题的起点,例如瞬态动力学分析、谐响应分析和谱分析,其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析,后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。ANSYS提供的模态提取方法有:子空间法(subspace)、分块法(blocklancets),缩减法(reduced/householder)、动态提取法(powerdynamics)、非对称法(unsymmetric),阻尼法(damped),QR阻尼法(QRdamped)等,大多数分析都可使用子空间法、分块法、缩减法。ANSYS的模态分析是线形分析,任何非线性特性,例如塑性、接触单元等,即使被定义了也将被忽略。2.模态分析操作过程一个典型的模态分析过程主要包括建模、模态求解、扩展模态以及观察结果四个步骤。(1).建模模态分析的建模过程与其他分析类型的建模过程是类似的,主要包括定义单元类型、单元实常数、材料性质、建立几何模型以及划分有限元网格等基本步骤。(2).施加载荷和求解包括指定分析类型、指定分析选项、施加约束、设置载荷选项,并进行固有频率的求解等。指定分析类型,MainMenu-Solution-AnalysisType-NewAnalysis,选择Modal。指定分析选项,MainMenu-Solution-AnalysisType-AnalysisOptions,选择MODOPT(模态提取方法〕,设置模态提取数量MXPAND.定义主自由度,仅缩减法使用。施加约束,MainMenu-Solution-DefineLoads-Apply-Structural-Displacement。求解,MainMenu-Solution-Solve-CurrentLS。(3).扩展模态如果要在POSTI中观察结果,必须先扩展模态,即将振型写入结果文件。过程包括重新进入求解器、激话扩展处理及其选项、指定载荷步选项、扩展处理等。激活扩展处理及其选项,MainMenu-Solution-LoadStepOpts-Expansionpass-SingleExpand-Expandmodes。指定载荷步选项。扩展处理,MainMenu-solution-Solve-CurrentLS。注意:扩展模态可以如前述办法单独进行,也可以在施加载荷和求解阶段同时进行。本例即采用了后面的方法(4).查看结果模态分析的结果包括结构的频率、振型、相对应力和力等实例1均匀直杆的模态分析3.问题描述及解析解图8-1所示为一根长度为L的等截面直杆,一端固定,一端自由。己知杆材料的弹性模量E=2×1011N/m2护,密度p=7800kg/m3,杆长L=0.1m.要求计算直杆纵向振动的固有频率。根据振动学理论,假设直杆均匀伸缩,如图8-1所示等截面直杆纵向振动第i阶固有频率为i=Li2)12(Erad/s(i=1,2…)将角频率i转化为周频率f,并将已知参数代入,可得fi=2i=Li412E=1.0412i780010211=12659(2i-1)Hz(8-1)按式(8-1)计算出直杆的前5阶频率,列表如表8-1所示。阶次12345频率(Hz)126593797863296886151139334.分析步骤4.1改变工作名拾取菜单UtilityMenu-File-ChangeJobname弹出如图8-2所示的对话框,在[/FILNAM]文本框中输入EXAMPLE8,单击“OK按钮。4.2创建单元类型拾取菜单MainM--Preprocessor-ElementType-Add/Edit/Delete,弹出对话框,单击Add按钮:弹出对话框,在左侧列表中选择StructuralSolid,在右侧列表中选择Brick20node186,单击OK按钮:单击对话框的“Close按钮。4.3定义材料属性拾取菜单MainMenu-Preprocessor--MaterialProps-MaterialModels.弹出对话框,在右侧列表中依次双击Structural,Linear,Elastic,Isotropic,弹出对话框,在EX文本框中输入2e11(弹性模量).在PRXY文本框中输入0.3(泊松比).单击OK按钮:再双击右侧列表中Structural下Density,弹出对话框,在DENS文本框中输入7800(密度),单击OK按钮。然后对话框。4.4创建块拾取菜单MainMenu-Preprocessor–Modeling-Create-Volumes-Block–ByDimension。弹出对话框,在X1,X2文本框中输入0,0.01,在Y1,Y2文本框中输入0,0.01,在Z1,Z2文本框中输入0,0,1,单击OK按钮。4.5划分单元拾取菜单MainMenu-Preprocessor-Meshing-MeshTool。弹出对话框,单击''Size.Controls'',区域中Lines后Set,按钮,弹出拾取窗口,任意拾取块x轴和y轴方向的边各一条(短边),单击OK按钮,弹出对话框,在NDIV文本框中输入3,单击Apply按钮:再次弹出拾取窗口,拾取块z轴方向的边(长边),单击OK按钮。在NDIV文木框中输入15,单击OK按钮。在Mesh区域,选择单元形状为Hex(六面体),选择划分单元的方法为Mapped(映射)单击Mesh按钮,弹出拾取窗口,单击OK按钮。图1单元划分4.6施加约束拾取菜单MainMenu--Solution-DefineLoads-Apply-Structural-Displacement-OnAreas弹出拾取窗口,拾取z=0的平面,单击OK按钮.弹出对话框,在列表中选择UZ,单击Apply按钮;再次弹出拾取窗口,拾取y=0的平面,单击OK按钮,弹出对话框,在列表中选择UY,单击Apply按钮再次弹出拾取窗口,拾取x=0的平面,单击OK按钮,弹出对话框,在列表中选择UX,单击OK按钮。所加约束与图8-1不同.主要是为了与推导式(8-1)所作的轴向振动假设一致。约束施加的正确与否,对结构模态分析的影响十分显著,因此对于该问题应十分注意,保证对模型施加的约束与实际情况尽量符合。4.7指定分析类型拾取菜单MainMenu-Solution-AnalysisType-NewAnalysis。弹出对话框,选择TypeofAnalysis为Modal,单击OK按钮。4.8指定分析选项拾取菜单MainMenu-Solution-AnalysisType-AnalysisOptions。弹出对话框,在No.ofmodestoextract文本框中输入5,单击“OK按钮:弹出BlockLanczosMethod,单击OK按钮。4.9指定要扩展的模态数拾取菜单MainMenu-Solution-LoadStepOpts-Expansionpass-SingleExpand-Expandmodes。弹出对话框,在NMODE文本框中输入5,单击OK按钮。4.10求解拾取菜单MainMenu-Solution-Solve-CurrentLS。单击“SolveCurrentLoadStep”对话框的“OK”按钮。出现“Solutionisdone!”提示时,求解结束,即可查看结果了。图2求解结果4.11列表固有频率拾取菜单MainMe-GeneralPostproc-ResultsSummary。弹出窗口,列表中显示了模型的前5阶频率,与表8-1相对照,可以看出结果虽然存在一定的误差,但与解析解是基本符合的。查看完毕后,关闭该窗口。4.12从结果文件读结果拾取菜单MainMenu-GeneralPostproc-ReadResults-FirstSet.4.13改变视点为便于更好地观察模型的模态。拾取菜单UtilityMenu-PlotCtrls-PanZoomRotate.在弹出的对话框中,单击“Left”按钮.4.14用动画观察模型的一阶模态拾取菜单UtilityMenu-PlotCtrls-Animate-ModeShape。弹出对话框,单击“OK”按钮。观察完毕,单击“AnimationController”对话框的close按钮。4.15观察其余各阶模态拾取菜单MainMenu-GeneralPostproc-ReadResults-NextSet.依次将其余各阶模态的结果读入,然后重复步骤15观察完模型的各阶模态后,请读者自行分析频率结果产生误差的原因,并改进以上分析过程。5命令流/CLEAR/FILNAME,EXAMPLE8/PREP7ET,1,SOLID186MP,EX,1,2E11MP,PRXY1,0.3MP,DENS,1,7800BLOCK,0,0.01,0,0.01,0,0.1LESIZE,1,,,3LESIZE.2,,,3LESIZE,9,,,15MSHAPE,OMSHKEY,1VMESH,1FINISH/SOLUANTYPE,MODALMODOPT,LANB,5MXPAND,5DA,1,UZDA,3,UYDA,5,UXSOLVESAVEFINISH/POST1SET,LISTSET,FIRST/VIEW,1,-1/REPLOTPLDIANMODE,10,0.5,,0SET,NEXTPLDIANMODE,10,0.5,,0FINISH