模态提取方法之PowerDynamics法
一.模态分析的定义和目的
模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型)，即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。
ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析ANSYS提供的模态提取方法有:子空间法(subspace)、分块法(block lancets),缩减法(reduced/householder)、动态提取法(power dynamics)、非对称法(unsymmetric),阻尼法(damped), QR 阻尼法(QR damped)等，大多数分析都可使用子空间法、分块法、缩减法。
模态分析可以使结构设计避免共振或以特定频率进行振动（例如扬声器），使工程师可以认识到结构对于不同类型的动力载荷是如何响应的；有助于在其它动力分析中估算求解控制参数（如时间步长）。
二.模态分析概念
1.通用运动方程
[]{}[]{}[]{}()
{}t
M=
+
+
u
u
F
K
u
C
2.假定为自由振动并忽略阻尼：
[]{}[]{}{}0=+u K u
M
3. 假定为谐运动：
[][](){}{}
02
=-u M K ω
这个方程的根是i ，即特征值，i 的范围从1到自由度的数目，相应的向量是 {u}I ，即特征向量。模态分析假定结构是线性的(如, [M]和[K]保持为常数)，简谐运动方程u = u0cos(wt), 其中 w 为自振圆周频率(弧度/秒）。特征值的平方根是 wi ， 它是结构的自然圆周频率（弧度/秒），并可得出自然频率 fi = wi /2p 特征向量 {u}i 表示振型， 即假定结构以频率 fi 振动时的形状模态提取 是用来描述特征值和特征向量计算的术语。 三.模态提取方法- Power Dynamics 法
PowerDynamics 法适用于提取很大的模型（100.000个自由度以上）的较少振型（< 20）。这种方法明显比 Block Lanczos 法或子空间法快，但是：需要很大的内存；当单元形状不好或出现病态矩阵时，用这种方法可能不收敛；建议只将这种方法作为对大模型的一种备用方法。子空间技术使用Power 求解器(PCG)和 一直质量矩阵；不执行Sturm 序列检查(对于遗漏模态); 它可能影响多个重复频率的模型；一个包含刚体模态的模型, 如果你使用PowerDynamics 方法,必须执行RIGID 命令(或者在分析设置对话框中指定RIGID 设置)。
Power Dynamics法的四个步骤
1.模型建立：模态分析是线性分析，如果在分析中指定了非线性单元，程序在计算过程中将忽略其非线性行为，故模态分析尽可能选用线性单元。在材料特性中密度DENS一定要定义，以构建质量矩阵；另外必须指定弹性模量EX。材料的性质可以是线性的、非线性的、恒定的或与温度相关的，但非线性性质将被忽略。
2.选择分析类型和分析选项：进入/SOLU中定义模态分析，声明模态分析方法，结构外力负载（通常指结构约束条件，如果有结构外力，则是预应力问题），主自由度的选择（如选用降阶法）。求解，退出/SOLU。
3.施加边界条件并求解：进入/SOLU，将所得结果扩展至全结构，求解，并保存至结果文件以便在后处理器中检查结果。
4.评价结果进入通用后处理器POST1，列出各自然频率，观察振型，观察模态应力。
四.综合实例分析
(一)现有方钢，截面面积25*25mm，长度为450mm，弹性模量为3*105 MPa，泊松比为0.3，密度为7.9*10-6kg/mm3.
模型
（二）GUI操作
1.定义工作文件名和工作标题
（1）定义工作文件名。
GUI：【Utility Menu】/【File】/【Change Jobname】。
弹出下图对话框，在此出现的对话框输入“liang”，并将“New long and error file”复选框选为“yes”，单击【ok】按钮。
图1
（2）定义工作标题。
GUI：【Utility Menu】/【File】/【Change Title】命令。
在出现的对话框中输入“liangmotai”，单击【ok】按钮，如下图。
图2
（3）重新显示。
GUI：【Utility Menu】/【Plot】/【Replot】。
2.过滤界面
GUI：【Main Menu】/【Preferences】。
弹出如下对话框，选中“structural”和“h-Method”项，单击【ok】。
图3
3.创建单元类型
GUI：【Main Menu】/【Preprocessor】/【Elements Type】/
【Add/Edit/Delete】。
在弹出下图左侧对话框中，单击【Add】按钮。弹出的对话框，在左侧列表中选择“Structural Beam”，在右侧列表中选择“2 node 188”单击【ok】按钮，单击右侧框图的【Close】按钮。
图4
图5
4.定义材料特性
GUI：【Main Menu】/【Preprocessor】/【Materials Props】/【Material Models】。
弹出的如图6所示对话框，在右表中依次选择【Structural】/
【Linear】/【Elastic】/【Isotropic】命令，弹出图7所示对话框，在“EX”文本框输入3e5（弹性模量），在“PRXY”文本框输入0.3（泊松比），单击【ok】按钮，单击【Density】按钮，弹出如图8所示对话框，在“DENS”文本框输入7.9e-6（密度），然后关闭图8所示对话框。
图6定义材料特性
图7定义杨氏模量与泊松比
图8定义密度
5.设置截面
GUI：【Main Menu】/【Preprocessor】/【Section】/【Beam】/【Common Section】。
弹出如图9所示对话框，在“B”“H”中分别输入25和25，单击【ok】按钮。
图9 定义截面
6.建模和单元划分
（1）创建关键点
GUI：【Main Menu】/【Preprocessor】/【Modeling】/【Create】/【Keypoints】/【In Active CS】。
弹出如图10所示对话框，在“Keypoints Number”文本框中输入1，在“X，Y，Z”分别输入0，0，0，单击【Apply】按钮，在“Keypoints Number”输入2，在“X，Y，Z”分别输入450，0，0单击【ok】按钮。
图10
（2）创建直线
GUI:【Main Menu】/【Preprocessor】/【Modeling】/【Create】/【Lines】/【Straight Line】。
弹出窗口，分别选取关键点1和2，创建直线，单击【ok】。
图11 创建直线
（3）划分单元
GUI:【Main Menu】/【Preprocessor】/【Meshing】/【MeshTool】。弹出图12所示对话框，单击“Size Controls”区域中“Lines”后【Set】按钮，弹出拾取对话框，拾取直线，单击【ok】弹出图13对话框，在“NDIV”中输入100，单击【ok】。在“Mesh”区采取默认格式，单击【Mesh】按钮，弹出拾取对话框，拾取直线，单击【ok】。
图12 单元划分编辑框
图13 编辑单元格数
（4）显示单元格
GUI:【Utility Menu】/【PlotCtrls】/【Lines】/【Syle】/【Size and Shape】。
弹出14图所示对话框，将“Display of element”选中单击【ok】。
图14 单元显示
（5）改变视角
GUI：【Utility Menu】/【PlotCtrls】/【Pan Zoom Rotate】在弹出的拾取框中依次单击【Iso】【Fit】，单击【Close】，如图15所示，最终图形如图16
图15 改变视角
图16 建立模型
7.定义分析类型
GUI:【Main Menu】/【Solution】/【Analysis Type】/【New Analysis】. 在弹出的对话框中选取“Modal”单击【ok】按钮。
图17
8.定义模态分析方法
GUI：【Main Menu】/【Solution】/【Analysis Type】/【Analysis Options】。选择“Block Lanczos”。在“NO. Of modes extract ”输入20 ，在“NO. Of modes to expand”中输入20，其他选项接受默认，单击【ok】，弹出19图，全部默认，单击【ok】。
图18 定义模态分析方法
图19参数设置
9.施加约束
GUI：【Main Menu】/【Solution】/【Define Loads】/【Apply】/【Structural】【Displacement】/【On Keypoints】。
在弹出图20对话框。拾取两端点，单击【ok】，弹出图21对话框。在列表中选择“All DOF”，单击【ok】按钮。
图20 施加约束
图21 约束面自由度
图22得到结果
10.模态扩展设置
GUI:【Main Menu】/【Solution】/【Load Step Opts】/【ExpansionPass】/【Single Expand】/【Expand Modes】。
弹出图23所示，勾选“Calculate elem results”，其他默认。
图23模态扩展设置
11.求解
GUI:【Main Menu】/【Solution】/【Solve】/【Current LS】。
如图24所示，单击“Solve Current Load Step”对话框的【ok】按钮，出现“Solution is done”提示时，如图25所示，求解完成。
图25开始求解
ANSYS模态分析实例 5.2ANSYS建模 该课题研究的弹性联轴器造型如下图5.2： 在ANSYS中建立模型，先通过建立如5.2所式二分之一的剖面图，通过绕中轴线旋转建立模拟模型如下图5.3
5.3单元选择和网格划分 由于模型是三给实体模型，故考虑选择三维单元，模型中没有圆弧结构，用六面体单元划分网格不会产生不规则或者畸变的单元，使分析不能进行下去，所以采用六面体单元。经比较分析，决定采用六面体八结点单元SOLID185，用自由划分的方式划分模型实体。课题主要研究对象是联轴器中橡胶元件，在自由划分的时候，中间件2网格选择最小的网格，smart size设置为1，两端铁圈的smart size设置为6，网格划分后模型如图5.4。 5.4边界约束 建立柱坐标系R-θ-Z，如5-5所示，R为径间，Z为轴向
选择联轴器两个铁圈的端面，对其面上的节点进行坐标变换，变换到如图5.5所示的柱坐标系，约束节点R，Z方向的自由度，即节点只能绕Z轴线转 5.5联轴器模态分析 模态分析用于确定设计中的结构或者机器部件振动特性（固有频率和振型），也是瞬态变动力学分析和谐响应分析和谱分析的起点。 在模态分析中要注意：ANSYS模态分析是线性分析，任何非线性因素都会被忽略。因此在设置中间件2的材料属性时，选用elastic材料。 5.5.1联轴器材料的设置 材料参数设置如下表5-1： 表5.1材料参数设置 表5.1材料参数设置 铁圈1 中间件2 铁圈3 泊松比0.3 0.4997 0.3 弹性模量Mpa 2E5 1.274E3 2E5 密度kg/m 7900 1000 7900 5.5.2联轴器振动特性的有限元计算结果及说明 求解方法选择Damped方法，频率计算结果如表5-2，振型结果为图5.6： 表5.2固有频率 SET TEME/FREQ LOAO STEP SUBSTEP CUMULATIVE 1 40.199 1 1 1 1 73.63 2 1 2 2 3 132.42 1 3 3 4 197.34 1 4 4
有限元分析概念 有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型：它是真实系统理想化的数学抽象。由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。 有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。如果采用高效的代数方程组求解方法，也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别： 1）非线性问题的方程是非线性的，一般需要迭代求解； 2）非线性问题不能采用叠加原理； 3）非线性问题不总有一致解，有时甚至没有解。 有限元求解非线性问题可分为以下三类：
1）材料非线性问题 材料的应力和应变是非线性的，但应力与应变却很微小，此时应变与位移呈线性关系，这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。 2）几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。 3）非线性边界问题 在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。 平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。
有限单元法与有限元分析 1.有限单元法 在数学中，有限元法（FEM，Finite Element Method）是一种为求解偏微分方程边值问题近似解的数值技术。求解时对整个问题区域进行分解，每个子区域都成为简单的部分，这种简单部分就称作有限元。它通过变分方法，使得误差函数达到最小值并产生稳定解。类比于连接多段微小直线逼近圆的思想，有限元法包含了一切可能的方法，这些方法将许多被称为有限元的小区域上的简单方程联系起来，并用其去估计更大区域上的复杂方程。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。 随着电子计算机的发展，有限单元法是迅速发展成一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 1.1.有限元法分析本质 有限元法分析计算的本质是将物体离散化。即将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算精度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。 1.2.特性分析 1）选择位移模式： 在有限单元法中，选择节点位移作为基本未知量时称为位移法；选择节点力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。 当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量如
压电变换器的自振频率分析及详细过程 1.模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型)，即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS提供的模态提取方法有:子空间法(subspace)、分块法(block lancets),缩减法(reduced/householder)、动态提取法(power dynamics)、非对称法(unsymmetric),阻尼法(damped), QR阻尼法(QR damped)等，大多数分析都可使用子空间法、分块法、缩减法。 ANSYS的模态分析是线形分析，任何非线性特性，例如塑性、接触单元等，即使被定义了也将被忽略。 2.模态分析操作过程 一个典型的模态分析过程主要包括建模、模态求解、扩展模态以及观察结果四个步骤。 (1).建模 模态分析的建模过程与其他分析类型的建模过程是类似的，主要包括定义单元类型、单元实常数、材料性质、建立几何模型以及划分有限元网格等基本步骤。 (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 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).查看结果 模态分析的结果包括结构的频率、振型、相对应力和力等
模态分析中的几个基本概念 物体按照某一阶固有频率振动时，物体上各个点偏离平衡位置的位移是满足一定的比例关系的，可以用一个向量表示，这个就称之为模态。模态这个概念一般是在振动领域所用，你可以初步的理解为振动状态，我们都知道每个物体都具有自己的固有频率，在外力的激励作用下，物体会表现出不同的振动特性。一阶模态是外力的激励频率与物体固有频率相等的时候出现的，此时物体的振动形态叫做一阶振型或主振型；二阶模态是外力的激励频率是物体固有频率的两倍时候出现，此时的振动外形叫做二阶振型，以依次类推。一般来讲，外界激励的频率非常复杂，物体在这种复杂的外界激励下的振动反应是各阶振型的复合。模态是结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。有限元中模态分析的本质是求矩阵的特征值问题，所以“阶数”就是指特征值的个数。将特征值从小到大排列就是阶次。实际的分析对象是无限维的，所以其模态具有无穷阶。但是对于运动起主导作用的只是前面的几阶模态，所以计算时根据需要计算前几阶的。一个物体有很多个固有振动频率（理论上无穷多个），按照从小到大顺序，第一个就叫第一阶固有频率，依次类推。所以模态的阶数就是对应的固有频率的阶数。振型是指体系的一种固有的特性。它与固有频率相对应，即为对应固有频率体系自身振动的形态。每一阶固有频率都对应一种振型。振型与体系实际的振动形态不一定相同。振型对应于频率而言，一个固有频率对应于一个振型。按照频率从低到高的排列，来说第一振型，第二振型等等。此处的振型就是指在该固有频率下结构的振动形态，频率越高则振动周期越小。在实验中，我们就是通过用一定的频率对结构进行激振，观测相应点的位移状况，当观测点的位移达到最大时，此时频率即为固有频率。实际结构的振动形态并不是一个规则的形状，而是各阶振型相叠加的结果。 固有频率也称为自然频率( natural frequency)。物体做自由振动时，其位移随时间按正弦或余弦规律变化，振动的频率与初始条件无关，而仅与系统的固有特性有关（如质量、形状、材质等），称为固有频率，其对应周期称为固有周期。 物体做自由振动时，其位移随时间按正弦规律变化，又称为简谐振动。简谐振动的振幅及初相位与振动的初始条件有关，振动的周期或频率与初始条件无关，而与系统的固有特性有关，称为固有频率或者固有周期。 物体的频率与它的硬度、质量、外形尺寸有关，当其发生形变时，弹力使其恢复。弹力主要与尺寸和硬度有关，质量影响其加速度。同样外形时，硬度高的频率高，质量大的频率低。一个系统的质量分布，内部的弹性以及其他的力学性质决定 模态扩展是为了是结果在后处理器中观察而设置的，原因如下： 求解器的输出内容主要是固有频率，固有频率被写到输出文件Jobname.OUT 及振型文件Jobnmae.MODE 中，输出内容中也可以包含缩减的振型和参与因子表，这取决于对分析选项和输出控制的设置，由于振型现在还没有被写到数据库或结果文件中，因此不能对结果进行后处理，要进行后处理，必须对模态进行扩展。在模态分析中，我们用“扩展”这个词指将振型写入结果文件。也就是说，扩展模态不仅适用于Reduced 模态提取方法得到的缩减振型，而且也适用与其他模态提取方法得到的完整振型。因此，如果想在后处理器中观察振型，必须先扩展模态。谱分析中的模态合并是因为激励谱是其实是由一系列的激励组合成的一个谱，里面的频率不会是只有一个，而不同的激励频率对于结构产生的结果是不一样的，对于结果的贡献也是不一样的，所以要选择模态组合法对模态进行组合，得到最终的响应结果。
ANSYS 模态分析实例 5.2ANSYS 建模 该课题研究的弹性联轴器造型如下图 5.2 : 图勺2弹性联轴器 1-联接柴油机大铁圈；茁橡胶膜片；3-联接电动机小铁圈 在ANSYS 中建立模型，先通过建立如 5.2所式二分之一的剖面图，通过绕中轴线 旋转建立模拟模型如下图 5.3资料个人收集整理，勿做商业用途 _.:q: 4 1(. 片三 _」」_止
5.3单元选择和网格划分 由于模型是三给实体模型，故考虑选择三维单元，模型中没有圆弧结构，用六面体单元划分网格不会产生不规则或者畸变的单元，使分析不能进行下去，所以采用六面体单元。经比较分析，决定采用六面体八结点单元SOLID185，用自由划分的方式划 分模型实体。课题主要研究对象是联轴器中橡胶元件，在自由划分的时候，中间件2 网格选择最小的网格，smart size设置为1,两端铁圈的smart size设置为6,网格划分 后模型如图5.4。资料个人收集整理，勿做商业用途 5.4边界约束 建立柱坐标系R- &Z,如5-5所示，R为径间，Z为轴向
选择联轴器两个铁圈的端面，对其面上的节点进行坐标变换，变换到如图5.5所示的柱坐标系，约束节点R，Z方向的自由度，即节点只能绕Z轴线转资料个人收集整理，勿做商业用途 5.5联轴器模态分析 模态分析用于确定设计中的结构或者机器部件振动特性（固有频率和振型），也是瞬态变动力学分析和谐响应分析和谱分析的起点。资料个人收集整理，勿做商业用途在模态分析中要注意：ANSYS模态分析是线性分析，任何非线性因素都会被忽略。因此在设置中间件2的材料属性时，选用elastic材料。资料个人收集整理，勿做商业用途 5.5.1联轴器材料的设置 材料参数设置如下表5-1 : 表5.1材料参数设置 表5.1材料参数设置 5.5.2联轴器振动特性的有限元计算结果及说明 求解方法选择Damped方法，频率计算结果如表5-2，振型结果为图5.6： 表5.2固有频率
基于ANSYS 的齿轮模态分析 齿轮传动是机械传动中最重要的传动部件，被广泛的应用在各个生产领域中，经常用在重要的场合；传动齿轮在工作过程中受到周期性载荷力的作用，有可能在标定转速内发生强烈的共振，动应力急剧增加，致使齿轮过早出现扭转疲劳和弯曲疲劳。静力学计算不能完全满足设计要求，因此有必要对齿轮进行模态分析，研究其振动特性，得到固有频率和主振型(自由振动特性)。同时，模态分析也是其它动力学分析如谐响应分析、瞬态动力学分析和谱分析的基础。 本文运用UG 对齿轮建模并用有限元软件ANSYS 对齿轮进行模态分析，为齿轮动态设计提供了有效的方法。 1.模态分析简介 由弹性力学有限元法，可得齿轮系统的运动微分方程为： []{}[]{}[]{}{()}M X C X K X F t ++= （1） 式中，[]M ，[]C ，[]K 分别为齿轮质量矩阵、阻尼矩阵和刚度矩阵；分别为齿轮振动加速度向量、速度向量和位移向量，{}X 、{}X 、{}X 分别为齿轮振动加速度向量、速度向量和位移向量，12{}{,, ,}T n X x x x =；{()}F t 为齿轮所受外界激振力向量，{}12{()},,T n F t f f f =。若无外力作用，即{}{()}0F t =，则得 到系统的自由振动方程。在求齿轮自由振动的频率和振型即求齿轮的固有频率和固有振型时，阻尼对它们影响不大，因此，可以作为无阻尼自由振动问题来处理 [2]。无阻尼项自由振动的运动方程为： []{}[]{}0M X K X += （2） 如果令 {}{}sin()X t φωφ=+ 则有 2{}{}sin()X t ωφωφ=+ 代入运动方程，可得 2([][]){}0i i K M ωφ-= （3） 式中i ω为第I 阶模态的固有频率，i φ为第I 阶振型，1,2, ,i n =。 2.齿轮建模 在ANSYS 中直接建模有一定的难度，考虑到其与多数绘图软件具有良好的数据接口，可以方便的转化，而UG 软件以其参数化、全相关的特点在零件造型方面表现突出，可以通过参数控制模型尺寸的变化，因此本文采用通过UG 软件对齿轮进行参数化建模，保存为IGES 格式，然后将模型导入到ANSYS 软件中的方法。设有模数m=2.5mm ，齿数z=20，压力角β=20°，齿宽b=14mm ，孔径为￠20mm 的标准齿轮模型。如图1
模态分析中的几个基本概念 一、模态定义：物体按照某一阶固有频率振动时，物体上各个点偏离平衡位置的位移是满足一定的比例关系的，可以用一个向量表示。 模态分析一般是在振动领域应用，每个物体都具有自己的固有频率，在外力的激励作用下，物体会表现出不同的振动特性： 一阶模态是外力的激励频率与物体固有频率相等的时候出现的，此时物体的振动形态叫做一阶振型或主振型； 二阶模态是外力的激励频率是物体固有频率的两倍时候出现，此时的振动外形叫做二阶振型，以依次类推。 一般来讲，外界激励的频率非常复杂，物体在这种复杂的外界激励下的振动反应是各阶振型的复合。 二、模态分析：模态是结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。 有限元中模态分析的本质是求矩阵的特征值问题，所以“阶数”就是指特征值的个数。将特征值从小到大排列就是阶次。 实际的分析对象是无限维的，所以其模态具有无穷阶。但是对于运动起主导作用的只是前面的几阶模态，所以计算时根据需要计算前几阶的。 一个物体有很多个固有振动频率（理论上无穷多个），按照从小到大顺序，第一个就叫第一阶固有频率，依次类推。所以模态的阶数就是对应的固有频率的阶数。 三、振型是指体系的一种固有的特性。它与固有频率相对应，即为对应固有频率体系自身振动的形态。每一阶固有频率都对应一种振型。振型与体系实际的振动形态不一定相同。振型对应于频率而言，一个固有频率对应于一个振型。按照频率从低到高的排列，来说第一振型，第二振型等等。此处的振型就是指在该固有频率下结构的振动形态，频率越高则振动周期越小。在实验中，我们就是通过用一定的频率对结构进行激振，观测相应点的位移状况，当观测点的位移达到最大时，此时频率即为固有频率。实际结构的振动形态并不是一个规则的形状，而是各阶振型相叠加的结果。 四、模态扩展是为了是结果在后处理器中观察而设置的，原因如下： 求解器的输出内容主要是固有频率，固有频率被写到输出文件Jobname.OUT及振型文件Jobnmae.MODE中，输出内容中也可以包含缩减
有限元仿真分析学习心得 1 有限元分析方法原理 有限元分析（FEA，Finite Element Analysis）利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。还利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 有限元法是随着电子计算机发展而迅速发展起来的一种工程力学问题的数值求解方法。20世纪50年代初，它首先应用于连续体力学领域—飞机结构静、动态特性分析之中，用以求得结构的变形、应力、固有频率以及阵型。由于其方法的有效性，迅速被推广应用于机械结构分析中。随着电子计算机的发展，有限元法从固体力学领域扩展到流体力学、传热学、电磁学、生物工程学、声学等。 随着计算机科学与应用技术的发展，有限元理论日益完善，随之涌现了一大批通用和专业的有限元计算软件。其中，通用有限元软件以ANSYS，MSC公司旗下系列软件为杰出代表，专业软件以ABAQUS、LS-DYNA、Fluent、ADAMS 为代表。 ANSYS作为最著名通用和有效的商用有限元软件之一，集机构、传热、流体、电磁、碰撞爆破分析于一体，具有强大的前后处理及计算分析能力，能够进行多场耦合，结构-热、流体-结构、电-磁场的耦合处理求解等。 有限元分析一般由以下基本步骤组成： ①建立求解域，并将之离散化成有限个单元，即将问题分解成单元和节点； ②假定描述单元物理属性的形(shape)函数，即用一个近似的连续函数描述每个单元的解； ③建立单元刚度方程； ④组装单元，构造总刚度矩阵； ⑤应用边界条件和初值条件，施加载荷； ⑥求解线性或者非线性微分方程组得到节点值，如不同节点的位移； ⑦通过后处理获得最大应力、应变等信息。 结构的离散化是有限元的基础。所谓离散化就是将分析的结构分割成为有限
高速旋转轮盘模态分析 在进行高速旋转机械的转子系统动力设计时，需要对转动部件进行模态分析，求解出其固有频率和相应的模态振型。通过合理的设计使其工作转速尽量远离转子系统的固有频率。而对于高速部件，工作时由于受到离心力的影响，其固有频率跟静止时相比会有一定的变化。为此，在进行模态分析时需要考虑离心力的影响。通过该实验掌握如何用ANSYS进行有预应力的结构的模态分析。 一．问题描述 本实验是对某高速旋转轮盘进行考虑离心载荷引起的预应力的模态分析，求解出该轮盘的前5阶固有频率及其对应的模态振型。轮盘截面形状如图所示，该轮盘安装在某转轴上以12000转/分的速度高速旋转。相关参数为：弹性模量EX＝2.1E5Mpa，泊松比PRXY＝0.3， 密度DENS＝7.8E-9Tn/mm 3。 1-5关键点坐标： 1(-10, 150, 0) 2(-10, 140, 0) 3(-3, 140, 0) 4(-4, 55, 0) 5(-15, 40, 0) L=10+（学号×0.1） RS=5 二．分析具体步骤 1.定义工作名、工作标题、过滤参数 ①定义工作名：Utility menu > File > Jobname ②工作标题：Utility menu > File > Change Title（个人学号） 2.选择单元类型 本实验将选用六面体结构实体单元来分析，但在建模过程中需要使用四边形平面单元，所有需要定义两种单元类型：PLANE42和SOLID45，具体操作如下： Main Menu >Preprocessor > Element Type > Add/Edit/Delete
①“ Structural Solid”→“ Quad 4node 42” →Apply（添加PLANE42为1号单元） ②“ Structural Solid”→“ Quad 8node 45” →ok（添加六面体单元SOLID45为2号单元） 在Element Types (单元类型定义)对话框的列表框中将会列出刚定义的两种单元类型：PLANE42、SOLID45，关闭Element Types (单元类型定义)对话框，完成单元类型的定义。 3.设置材料属性 由于要进行的是考虑离心力引起的预应力作用下的轮盘的模态分析，材料的弹性模量EX 和密度DENS必须定义。 ①定义材料的弹性模量EX Main Menu >Preprocessor > Material Props > Material Models> Structural > Linear > Elastic >Isotropic 弹性模量EX＝2.1E5 泊松比PRXY＝0.3 ②定义材料的密度DENS Main Menu >Preprocessor > Material Props > Material Models>density DENS =7.8E-9 4.实体建模 对于本实例的有限元模型，首先需要建立轮盘的截面几何模型，然后对其进行网格划分，最后通过截面的有限元网格扫描出整个轮盘的有限元模型。具体的操作过程如下。 ①创建关键点操作：Main Menu > Preprocessor > Modeling > Create > Keypoints > In Active CS 列出各点坐标值Utility menu >List > Keypoints >Coordinate only
模态提取方法之PowerDynamics法 一.模态分析的定义和目的 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型)，即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析ANSYS提供的模态提取方法有:子空间法(subspace)、分块法(block lancets),缩减法(reduced/householder)、动态提取法(power dynamics)、非对称法(unsymmetric),阻尼法(damped), QR 阻尼法(QR damped)等，大多数分析都可使用子空间法、分块法、缩减法。 模态分析可以使结构设计避免共振或以特定频率进行振动（例如扬声器），使工程师可以认识到结构对于不同类型的动力载荷是如何响应的；有助于在其它动力分析中估算求解控制参数（如时间步长）。 二.模态分析概念 1.通用运动方程 []{}[]{}[]{}() {}t M= + + u u F K u C
2.假定为自由振动并忽略阻尼： []{}[]{}{}0=+u K u M 3. 假定为谐运动： [][](){}{} 02 =-u M K ω 这个方程的根是i ，即特征值，i 的范围从1到自由度的数目，相应的向量是 {u}I ，即特征向量。模态分析假定结构是线性的(如, [M]和[K]保持为常数)，简谐运动方程u = u0cos(wt), 其中 w 为自振圆周频率(弧度/秒）。特征值的平方根是 wi ， 它是结构的自然圆周频率（弧度/秒），并可得出自然频率 fi = wi /2p 特征向量 {u}i 表示振型， 即假定结构以频率 fi 振动时的形状模态提取 是用来描述特征值和特征向量计算的术语。 三.模态提取方法- Power Dynamics 法 PowerDynamics 法适用于提取很大的模型（100.000个自由度以上）的较少振型（< 20）。这种方法明显比 Block Lanczos 法或子空间法快，但是：需要很大的内存；当单元形状不好或出现病态矩阵时，用这种方法可能不收敛；建议只将这种方法作为对大模型的一种备用方法。子空间技术使用Power 求解器(PCG)和 一直质量矩阵；不执行Sturm 序列检查(对于遗漏模态); 它可能影响多个重复频率的模型；一个包含刚体模态的模型, 如果你使用PowerDynamics 方法,必须执行RIGID 命令(或者在分析设置对话框中指定RIGID 设置)。
Ansys模态分析详细论述 1、有限元概述 将求解域分解成若干小域，有限元模型由单元组成，单元之间通过节点连接，并承受载荷，节点自由度是随着连接该点单元类型变化的。 1.1分析前准备 （1）研读相关理论基础； （2）参考别人的分析方法和思路； （3）考虑时间和设备，做适当的简化假设，设定条件、材料并决定分析方式；（4）了解力学现象、分析关键位置并预先评估。 1.2 Von Mises 应力 Von Mises 应力是非负值，应力表达式可表示为： 1.3结果的分析 （1）建立疏密不同的三至五种网络，选择适中密度，不能以存在应力集中点处的结果做对比； （2）检验网格，分析结果的合理性，选择安全系数，并且要分析应力集中的真实性与危险性。 （3）接触收敛速度的提高：在不影响结构的前提下，控制或减少接触单元生成数目，并采用线性搜索，与打开自适应开关来提高收敛速度。 2、模态分析中的几个基本概念 物体按照某一阶固有频率振动时，物体上各个点偏离平衡位置的位移是满足一定的比例关系的，可以用一个向量表示，这个就称之为模态。模态这个概念一般是在振动领域所用，可以初步的理解为振动状态，我们都知道每个物体都具有自己的固有频率，在外力的激励作用下，物体会表现出不同的振动特性。2.1主要模态 一阶模态是外力的激励频率与物体固有频率相等的时候出现的，此时物体的振动形态叫做一阶振型或主振型；二阶模态是外力的激励频率是物体固有频率
的两倍时候出现，此时的振动外形叫做二阶振型，以依次类推。一般来讲，外界激励的频率非常复杂，物体在这种复杂的外界激励下的振动反应是各阶振型的复合。模态是结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。有限元中模态分析的本质是求矩阵的特征值问题，所以“阶数”就是指特征值的个数。将特征值从小到大排列就是阶次。 实际的分析对象是无限维的，所以其模态具有无穷阶。但是对于运动起主导作用的只是前面的几阶模态，所以计算时根据需要计算前几阶的。一个物体有很多个固有振动频率（理论上无穷多个），按照从小到大顺序，第一个就叫第一阶固有频率，依次类推。 所以模态的阶数就是对应的固有频率的阶数。振型是指体系的一种固有的特性。它与固有频率相对应，即为对应固有频率体系自身振动的形态。每一阶固有频率都对应一种振型。振型与体系实际的振动形态不一定相同。振型对应于频率而言，一个固有频率对应于一个振型。按照频率从低到高的排列来说，有第一振型，第二振型等等，此处的振型就是指在该固有频率下结构的振动形态，频率越高则，振动周期越小。在实验中，我们就是通过用一定的频率对结构进行激振，观测相应点的位移状况，当观测点的位移达到最大时，此时频率即为固有频率。实际结构的振动形态并不是一个规则的形状，而是各阶振型相叠加的结果。 固有频率也称为自然频率(natural frequency)。物体做自由振动时，其位移随时间按正弦或余弦规律变化，振动的频率与初始条件无关，而仅与系统的固有特性有关（如质量、形状、材质等），称为固有频率，其对应周期称为固有周期。物体做自由振动时，其位移随时间按正弦规律变化，又称为简谐振动。简谐振动的振幅及初相位与振动的初始条件有关，振动的周期或频率与初始条件无关，而与系统的固有特性有关，称为固有频率或者固有周期。 物体的频率与它的硬度、质量、外形尺寸有关，当其发生形变时，弹力使其恢复。弹力主要与尺寸和硬度有关，质量影响其加速度。同样外形时，硬度高的频率高，质量大的频率低。一个系统的质量分布，内部的弹性以及其他的力学性质决定。 2.2模态扩展
第11章SAMCEF Field模态分析实例 11.1 连杆模态分析 11.1.1 概述 此例题是在上例的基础上，对连杆结构进行模态分析，求出连杆结构的前十阶固有频率和相应的模态振型。为初次接触Samcef Field的用户介绍如何利用应力分析模型进行模态分析，并给出模态分析的过程，以便用户能够容易地跟随操作。 在本书光盘中提供有包含此例题所有建模、分析和结果确认过程的最终数据库文件和动画文件。通过动画先对整个分析过程获得一定了解的话，可以进一步提高跟随操作的效果。 此例题所介绍的各阶段的分析步骤与一般实际工作中的分析过程基本相同。其具体内容如下。 11.1.2 分析数据 材料特性如下： 弹性模量 E = 2.1e11 N/m2； 泊松比ν = 0.3； 密度ρ = 7.8e3 Kg/m3。 右侧大圆孔固定。 11.1.3 模态分析过程 1. 读取应力分析模型 首先读取上例中生成的分析模型。 (1)将显示在桌面上的Samcef Field图标 或相应目录中的Samcef Field连 击两次以打开程序。 (2)在主菜单选择File\Open，浏览并选择模型文件crank-stress.sfield。 (3)单击
266 图11-1 转换分析类型 (2) 如图11-2所示，从弹出菜单中
均匀直杆的子空间法模态分析 1.模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型)，即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS提供的模态提取方法有:子空间法(subspace)、分块法(block lancets),缩减法(reduced/householder)、动态提取法(power dynamics)、非对称法(unsymmetric),阻尼法(damped), QR阻尼法(QR damped)等，大多数分析都可使用子空间法、分块法、缩减法。 ANSYS的模态分析是线形分析，任何非线性特性，例如塑性、接触单元等，即使被定义了也将被忽略。 2.模态分析操作过程 一个典型的模态分析过程主要包括建模、模态求解、扩展模态以及观察结果四个步骤。 (1).建模 模态分析的建模过程与其他分析类型的建模过程是类似的，主要包括定义单元类型、单元实常数、材料性质、建立几何模型以及划分有限元网格等基本步骤。 (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 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).查看结果 模态分析的结果包括结构的频率、振型、相对应力和力等
一个典型的ANSYS分析过程可分为以下6个步骤： 1定义参数 2创建几何模型 3划分网格 4加载数据 5求解 6结果分析 1定义参数 1.1指定工程名和分析标题 启动ANSYS软件，选择File→Change Jobname命令 选择File→Change Title菜单命令 1.2定义单位 (2) 设置计算类型 ANSYS Main Menu: Preference→Material Props →Material Models →Structural →OK (3) 定义分析类型 ANSYS Main Menu: Preprocessor →Loads →Analysis Type →New Analysis→STATIC →OK 1.3定义单元类型 选择Main Menu→Preprocessor→Element Type→Add/Edit/Delete命令 单击[Options]按钮，在[Element behavior]下拉列表中选择[Plane strs w/thk]选项，单击确定 1.4定义单元常数 在ANSYS程序主界面中选择Main Menu→Preprocessor→Real Constants→Add/Edit/Delete命令 单击[Add]按钮，进行下一个[Choose Element Type]对话框 1.5定义材料参数 在ANSYS程序主界面，选择Main Menu→Preprocessor→Material Props→Material Models命令 (1)选择对话框右侧Structural→Linear→Elastic→Isotropic命令，并单击[Isotropic]选项，接着弹出如下所示[Linear Isotropic Properties for Material Number 1]对话框。 在[EX]文本框中输入弹性模量“200000”，在[PRXY]文本框中输入泊松比“0.3”，单击OK 2创建几何模型 在ANSYS程序主界面，选择Main Menu→Preprocessor→Modeling→Creat→Areas→Rectangle →By 2Corners命令 选择Main Menu→Preprocessor→Modeling→Creat→Areas→Circle→Solid Circle命令 3网格划分(之前一定要进行材料的定义和分配) 选择Main Menu→Preprocessor→Modeling→Operate→Booleans→Subtract→Arears Circle命令 选择Main Menu→Preprocessor→Meshing→Mesh→Areas→Free命令，弹出实体选择对话框，单击[Pick All]按钮，得到如下所示网格 4加载数据 (1)选择Main Menu→Preprocessor→Loads→Define Loads→Apply→Structural→Displacement→On Lines命令， 出现如下所示对话框，选择约束[ALL DOF]选项，并设置[Displacement value]为0，单击OK。
有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC四个比较知名比较大的公司。 常见软件 有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC四个比较知名比较大的公司，其中ADINA、ABAQUS在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件，ANSYS、MSC进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析，但是除ADINA以外其它三个必须与别的软件搭配进行迭代分析，唯一能做到真正流固耦合的软件只有ADINA。 软件对比 ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。 结构分析能力排名：1、ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名：1、ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名：1、ADINA、ANSYS、MSC、ABAQUS 性价比排名：最好的是ADINA，其次ABAQUS、再次ANSYS、最后MSC ABAQUS软件与ANSYS软件的对比分析 1．在世界范围内的知名度 两种软件同为国际知名的有限元分析软件，在世界范围内具有各自广泛的用户群。ANSYS软件在致力于线性分析的用户中具有很好的声誉，它在计算机资源的利用，用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题，其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。 由于ANSYS产品进入中国市场早于ABAQUS，并且在五年前ANSYS的界面是当时最好的界面之一，所以在中国，ANSYS软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北
有限元分析方法在工程中的应用 Application of finite element analysis method in Engineering 一、引言 从20世纪50年代诞生到现在，有限元方法和技术经历了60年的发展历程，已经成为当今科学与工程领域中分析和求解微分方程的系统化数值计算方法。由于有限元分析方法适用性强、形式简单、理论可靠等众多优点，近年来已被推广应用到航空航天、土木建筑、机械等相关科学领域。本文以ANSYS软件为例，介绍其功能和应用，包括几何建模技术、网格划分与有限元建模技术、施加载荷与求解过程、结果后处理技术等。图1是用有限元方法分析工程问题时的具体步骤[1]。 本文以车轮钢的疲劳性能研究为例，介绍有限元分析方法在其中的应用。 图1. 有限元方法进行计算机辅助工程分析的步骤 二、ANSYS操作步骤 ANSYS的基本操作步骤包括建模、划分网格、加载求解和后处理等步骤。进入ANSYS系统后有六个系统，提供使用者和软件之间的交流凭借这六个窗口可以实现输入命令、检查模型的建立、观察分析结果及图形输出与打印。ANSYS
各窗口及工具条如图2所示。 图2. ANSYS的窗口及工具条 1、建立模型 首先必须指定作业名和分析标题，接着使用PREP7前处理器定义单元类型、单元实常数、材料特性，然后建立几何模型。需要注意的是，ANSYS的GUI界面下没有类似WORD中的后退操作按钮，所以就出现了一个常见问题：做错一步操作如何后退？这里可以采用三种方法：（1）建模阶段可以使用Delete(删除)图元命令，划分网格阶段可以使用Clear（清除）单元命令。（2）每完成一个模块的操作，都用SA VE AS保存数据到不同名的数据库文件中，出错后点击Resum Form恢复。（3）使用命令：UNDO，ON以便激活ANSYS内部的返回命令。 本文以车轮钢为例，建立好的模型与图2类似，只是未划分网格。 2、单元网格划分 一个实体模型进行网格划分（meshing）之前必须指定所产生的单元属性（element attribute）。ANSYS有限元网格划分是进行数值模拟分析至关重要的一步，它直接影响着后续数值计算分析结果的精确性。ANSYS软件平台提供了映射网格划分和自由网格划分的策略。映射网格划分用于曲线、曲面、实体的网格划分方法，自由网格划分方法用于空间自由曲面和复杂实体。