《ANSYS模态分析实例和详细过程21461.pdf》由会员分享,可在线阅读,更多相关《ANSYS模态分析实例和详细过程21461.pdf(31页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、均匀直杆的子空间法模态分析 1.模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型),即结构的固有频率和振型,它们是承受动态载荷结构设计中的重要参数。同时,也可以作为其它动力学分析问题的起点,例如瞬态动力学分析、谐响应分析和谱分析,其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。ANSYS 的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析,后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。ANSYS 提供的模态提取方法有:子空间法(subspace)、分块
2、法(block lancets),缩减法(reduced/householder)、动态提取法(power dynamics)、非对称法(unsymmetric),阻尼法(damped),QR 阻尼法(QR damped)等,大多数分析都可使用子空间法、分块法、缩减法。ANSYS 的模态分析是线形分析,任何非线性特性,例如塑性、接触单元等,即使被定义了也将被忽略。2.模态分析操作过程 一个典型的模态分析过程主要包括建模、模态求解、扩展模态以及观察结果四个步骤。(1).建模 模态分析的建模过程与其他分析类型的建模过程是类似的,主要包括定义单元类型、单元实常数、材料性质、建立几何模型以及划分有限元
3、网格等基本步骤。(2).施加载荷和求解 包括指定分析类型、指定分析选项、施加约束、设置载荷选项,并进行固有频率的求解等。指定分析类型,Main Menu-Solution-Analysis Type-New Analysis,选择 Modal。指定分析选项,Main Menu-Solution-Analysis Type-Analysis Options,选择 MODOPT(模态提取方法,设置模态提取数量 MXPAND.定义主自由度,仅缩减法使用。施加约束,Main Menu-Solution-Define Loads-Apply-Structural-Displacement。求解,Main
4、 Menu-Solution-Solve-Current LS。(3).扩展模态 如果要在 POSTI 中观察结果,必须先扩展模态,即将振型写入结果文件。过程包括 重新进入求解器、激话扩展处理及其选项、指定载荷步选项、扩展处理等。激活扩展处理及其选项,Main Menu-Solution-Load Step Opts-Expansionpass-Single Expand-Expand modes。指定载荷步选项。扩展处理,Main Menu-solution-Solve-Current LS。注意:扩展模态可以如前述办法单独进行,也可以在施加载荷和求解阶段同时进行。本例即采用了后面的方法(4
5、).查看结果 模态分析的结果包括结构的频率、振型、相对应力和力等 实例 1 均匀直杆的模态分析 3.问题描述及解析解 图 8-1 所示为一根长度为 L 的等截面直杆,一端固定,一端自由。己知杆材料的弹性模量 E=2 1011 N/m2护,密度 p=7800 kg/m3,杆长 L=0.1 m.要求计算直杆纵向振动的固有频率。根据振动学理论,假设直杆均匀伸缩,如图 8-1 所示等截面直杆纵向振动第 i 阶固有频率为 i=Li2)12(E rad/s (i=1,2)将角频率i转化为周频率 f,并将已知参数代入,可得 fi=2i=Li412 E=1.0412i780010211=12659(2i-1)
6、Hz (8-1)按式(8-1)计算出直杆的前 5 阶频率,列表如表 8-1 所示。阶次 1 2 3 4 5 频率(Hz)12659 37978 63296 88615 113933 4.分析步骤 改变工作名 拾取菜单 Utility Menu-File-Change Jobname 弹出如图 8-2 所示的对话框,在/FILNAM 文本框中输入 EXAMPLE8,单击“OK按钮。创建单元类型 拾取菜单 Main M-Preprocessor-Element Type-Add/Edit/Delete,弹出对话框,单击Add按钮:弹出对话框,在左侧列表中选择Structural Solid,在右侧
7、列表中选择Brick 20node 186,单击OK按钮:单击对话框的“Close按钮。定义材料属性 拾取菜单 Main Menu-Preprocessor-Material Props-Material Models.弹出对话框,在右侧列表中依次双击Structural,Linear,Elastic,Isotropic,弹出对话框,在EX文本框中输入2e11(弹性模量).在PRXY文本框中输入(泊松比).单击OK按钮:再双击右侧列表中Structural下Density,弹出对话框,在DENS文本框中输入 7800(密度),单击OK按钮。然后对话框。创建块 拾取菜单 Main Menu-Pr
8、eprocessor Modeling-Create-Volumes-Block By Dimension。弹出对话框,在X1,X2文本框中输入 0,,在Y1,Y2文本框中输入 0,在Z1,Z2文本框中输入 0,0,1,单击OK按钮。划分单元 拾取菜单 Main Menu-Preprocessor-Meshing-MeshTool。弹出对话框,单击Size.Controls,区域中Lines后Set,按钮,弹出拾取窗口,任意拾取块 x 轴和 y 轴方向的边各一条(短边),单击 OK按钮,弹出对话框,在NDIV文本框中输入 3,单击Apply按钮:再次弹出拾取窗口,拾取块 z 轴方向的边(长边)
9、,单击OK按钮。在NDIV文木框中输入 15,单击OK按钮。在 Mesh 区域,选择单元形状为Hex(六面体),选择划分单元的方法为Mapped(映射)单击 Mesh 按钮,弹出拾取窗口,单击OK按钮。图 1 单元划分 施加约束 拾取菜单 Main Menu-Solution-Define Loads-Apply-Structural-Displacement-On Areas 弹出拾取窗口,拾取 z=0 的平面,单击OK按钮.弹出对话框,在列表中选择UZ,单击Apply按钮;再次弹出拾取窗口,拾取 y=0 的平面,单击OK按钮,弹出对话框,在列表中选择UY,单击Apply按钮再次弹出拾取窗口
10、,拾取 x=0 的平面,单击OK按钮,弹出对话框,在列表中选择UX,单击OK按钮。所加约束与图 8-1 不同.主要是为了与推导式(8-1)所作的轴向振动假设一致。约束施加的正确与否,对结构模态分析的影响十分显著,因此对于该问题应十分注意,保证对模型施加的约束与实际情况尽量符合。指定分析类型 拾取菜单 Main Menu-Solution-Analysis Type-New Analysis。弹出对话框,选择Type of Analysis为Modal,单击OK按钮。指定分析选项 拾取菜单 Main Menu-Solution-Analysis Type-Analysis Options。弹出对
11、话框,在No.of modes to extract文本框中输入 5,单击“OK按钮:弹出Block Lanczos Method,单击OK按钮。指定要扩展的模态数 拾取菜单 Main Menu-Solution-Load Step Opts-Expansionpass-Single Expand-Expand modes。弹出对话框,在NMODE文本框中输入 5,单击OK按钮。求解 拾取菜单 Main Menu-Solution-Solve-Current LS。单击“Solve Current Load Step”对话框的“OK”按钮。出现“Solution is done!”提示时,求解
12、结束,即可查看结果了。图 2 求解结果 列表固有频率 拾取菜单 Main Me-General Postproc-Results Summary。弹出窗口,列表中显示了模型的前 5 阶频率,与表 8-1 相对照,可以看出结果虽然存在一定的误差,但与解析解是基本符合的。查看完毕后,关闭该窗口。从结果文件读结果 拾取菜单 Main Menu-General Postproc-Read Results-First Set.改变视点 为便于更好地观察模型的模态。拾取菜单 Utility Menu-PlotCtrls-Pan Zoom Rotate.在弹出的对话框中,单击“Left”按钮.用动画观察模型
13、的一阶模态 拾取菜单 Utility Menu-PlotCtrls-Animate-Mode Shape。弹出对话框,单击“OK”按钮。观察完毕,单击“Animation Controller”对话框的close按钮。观察其余各阶模态 拾取菜单 Main Menu-General Postproc-Read Results-Next Set.依次将其余各阶模态的结果读入,然后重复步骤 15 观察完模型的各阶模态后,请读者自行分析频率结果产生误差的原因,并改进以上分 析过程。5 命令流/CLEAR/FILNAME,EXAMPLE8/PREP7 ET,1,SOLID186 MP,EX,1,2E11 MP,PRXY 1,MP,DENS,1,7800 BLOCK,0,0,0,LESIZE,1,3,3 LESIZE,9,15 MSHAPE,O MSHKEY,1 VMESH,1 FINISH/SOLU ANTYPE,MODAL MODOPT,LANB,5 MXPAND,5 DA,1,UZ DA,3,UY DA,5,UX SOLVE SAVE FINISH/POST1 SET,LIST SET,FIRST/VIEW,1,-1 /REPLOT PLDI ANMODE,10,0 SET,NEXT PLDI ANMODE,10,0 FINISH
限制150内