全面讲解ANSYS高速旋转轮盘模态分析讲解 高速旋转轮盘模态分析 在进行高速旋转机械的转子系统动力设计时，需要对转动部件进行模态分析，求解出其固有频率和相应的模态振型。通过合理的设计使其工作转速尽量远离转子系统的固有频率。而对于高速部件，工作时由于受到离心力的影响，其固有频率跟静止时相比会有一定的变化。为此 ，在进行模态分析时需要考虑离心力的影响。我通过该例子学习到了如何用ANSYS进行有预应力的结构的模态分析。 一．例子描述 本例子是对某高速旋转轮盘进行考虑离心载荷引起的预应力的模态分析，求解出该轮盘的 前10阶固有频率及其对应的模态振型。轮盘截面形状如图1所示，该轮盘安装在某转轴上以120 00转/分的速度高速旋转。相关参数为：弹性模量EX＝2.1E5Mpa，泊松比PRXY＝0.3，密度DE NS＝7.8E-9T/mm3。 图1、轮盘截面图 1-5关键点坐标： 1(-10, 150, 0) 2(-10, 140, 0) 3(-3, 140, 0) 4(-4, 55, 0)
5(-15, 40, 0) L=15 RS=5 二．A nsys求解的具体步骤 1.启动ansys，定义工作名、工作标题 ①定义工作名：Example of dynamic ②工作标题：dynamic analysis of a disc 2、选择单元类型 本例将选用六面体结构实体单元来分析，但在建模过程中需要使用四边形平面单元，所有需要定义两种单元类型：PLANE42和SOLID45，设置完成后，如图2，在Element Types (单元类型定义)对话框的列表框中将会列出刚定义的两种单元类型：PLANE42、 SOLID45， 图2、定义单元类型 3、设置材料属性 由于要进行的是考虑离心力引起的预应力作用下的轮盘的模态分析，材料的弹性模量EX 和密度DENS必须定义。 ①定义材料的弹性模量EX 弹性模量 EX＝2.1E5 泊松比 PRXY＝0.3 ②定义材料的密度DENS DENS =7.8E-9 4、建立实体模型 对于本实例的有限元模型，首先需要建立轮盘的截面几何模型，然后对其进行网格划分，最后通过截面的有限元网格扫描出整个轮盘的有限元模型。具体的操作过程如下。
/PREP7 !进入前处理 ANTYPE,STATIC !设置分析类型为静力结构分析 PSTRES,ON !用于后面的模态分析中考虑预应力(该开关不影响静力分析) ET,1,LINK10 !选取单元类型1(单向杆单元) KEYOPT,1,3,0 !设置仅承受拉应力，KEYOPT(3)=0 R,1,306796E-8,543248E-8 !设置实常数，包括绳索截面积(306796E-8)，初始应变(543248E-8) MP,EX,1,30E6 !定义材料的弹性模量(1号材料) MP,DENS,1,73E-5 !定义材料的密度(1号材料) N,1 ! 定义第1号节点 N,14,100 ! 定义第14号节点 FILL ! 均分填满第2号至第13号节点 E,1,2 !由节点1及节点2生成单元 EGEN,13,1,1 !依序复制生成13个单元 D,ALL,ALL ! 对所有节点施加固定约束 FINISH ! 前处理结束 /SOLU ! 进入求解模块，求解预应力引起的应力状态 SOLVE ! 求解 FINISH ! 退出求解模块 /POST1 ! 进入一般的后处理 ETABLE,STRS,LS,1 !针对LINK10单元,建立单元列表STRS,通过LS及特征号1来获得单元的轴向应力 *GET,STRSS,ELEM,13,ETAB,STRS !针对单元列表STRS, 提取13号单元的应力 FINISH ! 后处理结束 /POST26 ! 进入时间历程后处理，处理支反力 RFORCE,2,1,F,X !将1号节点上的x方向支反力提取，并存储到2号变量中 STORE ! 存储 *GET,FORCE,V ARI,2,EXTREM,VMAX !将2号变量的最大值赋给参数FORCE /SOLU ! 再次进入求解模块，模态分析 ANTYPE,MODAL ! 模态分析 MODOPT,SUBSP,3 ! 选择子空间迭代法，求3阶模态 MXPAND,3 ! 设定3阶模态扩展 PSTRES,ON ! 用于在模态分析中考虑预应力(还需在前面的静力分析中也同时打开) DDELE,2,UX,13 ! 删除从2号节点到13号节点上的UX约束 DDELE,2,UY,13 !删除从2号节点到13号节点上的UY约束 SOLVE !求解 *GET,FREQ1,MODE,1,FREQ ! 提取第1阶模态共振频率，并赋值给参数FREQ1 *GET,FREQ2,MODE,2,FREQ ! 提取第2阶模态共振频率，并赋值给参数FREQ2 *GET,FREQ3,MODE,3,FREQ ! 提取第3阶模态共振频率，并赋值给参数FREQ3 *STATUS !列出所有参数的实际内容
压电变换器的自振频率分析及详细过程 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).查看结果 模态分析的结果包括结构的频率、振型、相对应力和力等
/FILNAME,Elastic-Plasitc,1 /TITLE, Elastic-Plasitc Analysis !前处理。 /PREP7 !**定义梁单元189。 ET,1,BEAM189 !定义单元。 !**梁截面1。 SECTYPE, 1, BEAM, HREC, , 0 !定义梁截面。SECOFFSET, CENT SECDATA,50,100,6,6,6,6,0,0,0,0 !定义梁截面完成。 !**定义材料。 MPTEMP,,,,,,,, !定义弹塑性材料模型。MPTEMP,1,0 MPDATA,EX,1,,2.05e5 MPDATA,PRXY,1,,0.3 TB,BISO,1,1,2, TBTEMP,0 TBDATA,,150,18600,,,, !定义弹塑性材料模型。!**建立几何模型。 K,1, , , , K,2 ,900, K,3 ,,50 LSTR, 1, 2 !**网格划分。 FLST,5,1,4,ORDE,1 !定义网格密度。FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义网格密度完成。CM,_Y,LINE !网格划分。 LSEL, , , , 1 CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , 3, ,1 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 LMESH, 1 !网格划分完成。 !施加载荷及求解。 FINISH /SOL
!**施加约束。 FLST,2,1,3,ORDE,1 !施加约束。FITEM,2,1 /GO DK,P51X, , , ,0,UX,UY,UZ,ROTX, , , FLST,2,1,3,ORDE,1 FITEM,2,2 /GO DK,P51X, , , ,0,UY,UZ,ROTX, , , , !施加约束完成。 !**加载。 FLST,2,50,2,ORDE,2 FITEM,2,1 FITEM,2,-50 SFBEAM,P51X,1,PRES,100, , , , , , LSWRITE,1, !定义载荷步1完成。FLST,2,50,2,ORDE,2 !定义载荷步2。FITEM,2,1 FITEM,2,-50 SFEDELE,P51X,1,PRES LSWRITE,2, !定义载荷步2完成。!设定求解步并求解。 LSSOLVE,1,2,1,
Ansys非线性接触分析和设置 5.4.9 设置实常数和单元关键选项 程序使用20个实常数和数个单元关键选项，来控制面─面接触单元的接触。参见《ANSYS Elements Reference》中对接触单元的描述。 5.4.9.1 实常数 在20个实常数中，两个(R1和R2)用来定义目标面单元的几何形状。剩下的用来控制接触面单元。 R1和R2 定义目标单元几何形状。 FKN 定义法向接触刚度因子。 FTOLN 是基于单元厚度的一个系数，用于计算允许的穿透。 ICONT 定义初始闭合因子。 PINB 定义“Pinball"区域。 PMIN和PMAX 定义初始穿透的容许范围。 TAUMAR 指定最大的接触摩擦。 CNOF 指定施加于接触面的正或负的偏移值。 FKOP 指定在接触分开时施加的刚度系数。 FKT 指定切向接触刚度。 COHE 制定滑动抗力粘聚力。 TCC 指定热接触传导系数。 FHTG 指定摩擦耗散能量的热转换率。 SBCT 指定 Stefan-Boltzman 常数。 RDVF 指定辐射观察系数。 FWGT 指定在接触面和目标面之间热分布的权重系数。
FACT 静摩擦系数和动摩擦系数的比率。 DC 静、动摩擦衰减系数。 命令： R GUI：main menu> preprocessor>real constant 对实常数 FKN, FTOLN, ICONT, PINB, PMAX, PMIN, FKOP 和 FKT，用户既可以定义一个正值，也可以定义一个负值。程序将正值作为比例因子，将负值作为绝对值。程序将下伏单元的厚度作为ICON，FTOLN，PINB，PMAX 和 PMIN 的参考值。例如 ICON = 0.1 表明初始闭合因子是“0.1*下层单元的厚度”。然而，ICON = -0.1 则表示真实调整带是 0.1 单位。如果下伏单元是超单元，则将接触单元的最小长度作为厚度。参见图5-8。 图5-8 下层单元的厚度 在模型中，如果单元尺寸变化很大，而且在实常数如 ICONT, FTOLN, PINB, PMAX, PMIN 中应用比例系数，则可能会出现问题。因为从比例系数得到的实际结果，取决于下层单元的厚度，这就可能引起大、小单元之间的重大变化。如果出现这一问题，请用绝对值代替比例系数。 TCC, FHTG, SBCT, RDVF 和 FWGT 仅用于热接触分析[KEYOPT(1)=1]。 5.4.9.2 单元关键选项 每种接触单元都包括数个关键选项。对大多的接触问题，缺省的关键选项是合适的。而在某些情况下，可能需要改变缺省值。下面是可以控制接触行为的一些关键选项： 自由 度 KEYOPT(1) 接触算法(罚函数+拉格朗日乘子或罚函数) KEYOPT(2) 存在超单元时的应力状态(仅2D) KEYOPT(3)
ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立 的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT
1概述 接触问题是一种高度非线性行为，需要较大的计算资源，为了进行实为有效的计算，理解问题的特性和建立合理的模型是很重要的。 接触问题存在两个较大的难点：其一，在你求解问题之前，你不知道接触区域，表面之间是接触或分开是未知的，突然变化的，这随载荷、材料、边界条件和其它因素而定；其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。 一般的接触分类 接触问题分为两种基本类型：刚体─柔体的接触，半柔体─柔体的接触。 （1）刚-柔接触 在刚体─柔体的接触问题中，接触面的一个或多个被当作刚体，（与它接触的变形体相比，有大得多的刚度），一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体─柔体的接触，许多金属成形问题归为此类接触。 （2）柔-柔接触 柔体─柔体的接触，是一种更普遍的类型，在这种情况下，两个接触体都是变形体（有近似的刚度）。 2ANSYS接触能力 ANSYS支持三种接触方式：点─点，点─面，平面─面，每种接触方式使用的接触单元适用于某类问题。为了给接触问题建模，首先必须认识到模型中的哪些部分可能会相互接触，如果相互作用的其中之一是一点，模型的对立应组元是一个结点。如果相互作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元，有限元模型通过指定的接触单元来识别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元，至于ANSYS使用的接触单元和使用它们的过程，下面分类详述。 2.1点─点接触单元 点─点接触单元主要用于模拟点─点的接触行为，为了使用点─点的接触单元，你需要预先知道接触位置，这类接触问题只能适用于接触面之间有较小相对滑动的情况（即使在几何非线性情况下） 如果两个面上的结点一一对应，相对滑动又以忽略不计，两个面挠度（转动）保持小量，那么可以用点─点的接触单元来求解面─面的接触问题，过盈装配问题是一个用点─点的接触单元来模拟面─面的接触问题的典型例子。
模态分析步骤 第1步： 载入模型Plot>Volumes 第2步： 指定分析标题并设置分析范畴 1设置标题等Utility Menu>File>Change Title Utility Menu>File> Change Jobname Utility Menu>File>Change Directory 2选取菜单途径MainMenu>Preference ,单击Structure,单击OK第3步： 定义单元类型 MainMenu>Preprocessor>ElementType>Add/Edit/Delete,出现Element Types 对话框,单击Add出现Library of Element Types对话框,选择Structural Solid,再右滚动栏选择Brick 20node 95,然后单击OK，单击Element Types对话框中的Close 按钮就完成这项设置了。 第4步： 指定材料性能 选取菜单途径MainMenu>Preprocessor>MaterialProps>MaterialModels。出现DefineMaterialModelBehavior对话框,在右侧Structural>Linear>Elastic>Isotropic,指定材料的弹性模量和泊松系数，Structural>Density指定材料的密度，完成后退出即可。 第5步： 划分网格
选取菜单途径Main Menu>Preprocessor>Meshing>MeshTool,出现MeshTool 对话框，一般采用只能划分网格，点击SmartSize,下面可选择网格的相对大小（太小的计算比较复杂，不一定能产生好的效果，一般做两三组进行比较），保留其他选项，单击Mesh出现Mesh Volumes对话框，其他保持不变单击Pick All，完成网格划分。 第6步： 进入求解器并指定分析类型和选项 选取菜单途径Main Menu>Solution>Analysis Type>New Analysis，将出现New Analysis对话框,选择Modal单击OK。 选取Main Menu>Solution> Analysis Type>Analysis Options，将出现Modal Analysis对话框，选中Subspace模态提取法，在Number ofmodes to extract处输入相应的值（一般为5或10，如果想要看更多的可以选择相应的数字），单击OK，出现Subspace Model Analysis对话框，选择频率的起始值，其他保持不变，单击OK。 第7步： 施加边界条件.选取 MainMenu>Solution>Defineloads>Apply>Structural>Displacement,出现 ApplyU,ROTonKPS对话框，选择在点、线或面上施加位移约束,单击OK会打开约束种类对话框，选择（AllDOF,UX,UY,UZ）相应的约束,单击apply或OK即可。 第8步： 指定要扩展的模态数。选取菜单途径 MainMenu>Solution>LoadStepOpts>ExpansionPass>ExpandModes,出现Expand Modes对话框,在number of modes to expand处输入第6步相应的数字，单击OK 即可。（当选取MainMenu>Solution>AnalysisType>AnalysisOptions，将出现ModalAnalysis对话框，选中Subspace模态提取法，在Number of modes to extract处输入相应的值（一般为5或10，如果想要看更多的可以选择相应的数字），同时选择number of modes to expand输入相应值时，这步可以省略）
/FILNAME, Beam,1 ！定义工作文件名。 /TITLE, Beam Analysis ！定义工作标题。/PREP7 !定义单元。 ET,1,BEAM188 !定义材料属性。 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.1e5 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7.9e-6 ! 定义杆件截面■200。 SECTYPE, 1, BEAM, RECT, , 0 SECOFFSET, CENT SECDATA,10,10,0,0,0,0,0,0,0,0 !建立几何模型。 K,1, ,, , K,2,350,, , !生成立柱。 LSTR, 1, 2 !以上完成几何模型。 !以下进行网格划分。 FLST,5,1,4,ORDE,1 FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义单元大小。!分配、划分平板结构。 LMESH, 1 !分析类型施加载荷并求解。 FINISH /SOLUTION ANTYPE,2 ！定义分析类型及求解设置。MSAVE,0 ！模态提取方法。
MODOPT,LANB,10 EQSLV,SPAR MXPAND,10, , ,0 ！模态扩展设置。 LUMPM,0 PSTRES,0 MODOPT,LANB,10,0,0, ,OFF MXPAND,10,0,0,1,0.001, !施加约束。 FLST,2,2,3,ORDE,2 FITEM,2,1 FITEM,2,-2 /GO DK,P51X, , , ,0,ALL, , , , , , !求解。 FINISH /SOL /STATUS,SOLU SOLVE !定义谱分析。 FINISH /SOLUTION ANTYPE,8 SPOPT,PSD,10,1 PSDUNIT,1,DISP,386.4, PSDFRQ,1, ,13.8,40,50.6,73,120 !定义谱—频率表。PSDFRQ,1, ,134,178,233, , PSDV AL,1,1,4,0.6,3,5 PSDV AL,1,6,2,6, , FLST,2,2,1,ORDE,2 !施加谱。 FITEM,2,1 FITEM,2,-2 D,P51X,UX,1.0 PFACT,1,BASE, !计算PSD激励参与系数。PSDRES,DISP,REL !设置输出选项。PSDRES,VELO,OFF PSDRES,ACEL,OFF
参考ANSYS的中文帮助文件 接触问题（参考ANSYS的中文帮助文件） 当两个分离的表面互相碰触并共切时，就称它们牌接触状态。在一般的物理意义中，牌接触状态的表面有下列特点： 1、不互相渗透； 2、能够互相传递法向压力和切向摩擦力； 3、通常不传递法向拉力。 接触分类：刚性体－柔性体、柔性体－柔性体 实际接触体相互不穿透，因此，程序必须在这两个面间建立一种关系，防止它们在有限元分析中相互穿过。 ――罚函数法。接触刚度 ――lagrange乘子法，增加一个附加自由度（接触压力），来满足不穿透条件――将罚函数法和lagrange乘子法结合起来，称之为增广lagrange法。 三种接触单元：节点对节点、节点对面、面对面。 接触单元的实常数和单元选项设臵： FKN：法向接触刚度。这个值应该足够大，使接触穿透量小；同时也应该足够小，使问题没有病态矩阵。FKN值通常在0.1~10之间，对于体积变形问题，用值1.0（默认），对弯曲问题，用值0.1。 FTOLN：最大穿透容差。穿透超过此值将尝试新的迭代。这是一个与接触单元下面的实体单元深度（h）相乘的比例系数，缺省为0.1。此值太小，会引起收敛困难。 ICONT：初始接触调整带。它能用于围绕目标面给出一个“调整带”，调整带内任何接触点都被移到目标面上；如果不给出ICONT值，ANSYS根据模型的大小提供一个较小的默认值（＜0.03＝ PINB：指定近区域接触范围（球形区）。当目标单元进入pinball区时，认为它处于近区域接触，pinball区是围绕接触单元接触检测点的圆（二维）或球（三维）。可以用实常数PINB调整球形区（此方法用于初始穿透大的问题是必要的）PMIN和PMAX：初始容许穿透容差。这两个参数指定初始穿透范围，ANSYS把整个目标面（连同变形体）移到到由PMIN和PMAX指定的穿透范围内，而使其成为闭合接触的初始状态。初始调整是一个迭代过程，ANSYS最多使用20个迭代步把目标面调整到PMIN和PMAX范围内，如果无法完成，给出警告，可能需要修改几何模型。 TAUMAX：接触面的最大等效剪应力。给出这个参数在于，不管接触压力值多大，只要等效剪应力达到最大值TAUMAX，就会发生滑动。该剪应力极限值通常用于接触压力会变得非常大的情况。 CNOF：指定接触面偏移。+CNOF增加过盈、-CNOF减少过盈或产生间隙、CNOF能与几何穿透组合应用。 FKOP：接触张开弹簧刚度。针对不分离或绑定接触模型，需要设臵实常数FKOP，该常数为张开接触提供了一个刚度值。FKOP阻止接触面的分离；FKOP默认为1.0，用于建立粘结模型，用一个较小值（1e-5）去建立软弹簧模型。 FKT：切向接触刚度。作为初值，可以采用－FKT＝0.01*FKN，这是大多数ANSYS 接触单元的缺省值。 COHE：粘滞力。即没有法向压力时开始滑动的摩擦应力值。 FACT，DC：定义摩擦系数变化规律
一般的接触分类 (2) ANSYS接触能力 (2) 点─点接触单元 (2) 点─面接触单元 (2) 面─面的接触单元 (3) 执行接触分析 (4) 面─面的接触分析 (4) 接触分析的步骤： (4) 步骤1：建立模型，并划分网格 (4) 步骤二：识别接触对 (4) 步骤三：定义刚性目标面 (5) 步骤4：定义柔性体的接触面 (8) 步骤5：设置实常数和单元关键字 (10) 步骤六： (21) 步骤7：给变形体单元加必要的边界条件 (21) 步骤8：定义求解和载步选项 (22) 第十步：检查结果 (23) 点─面接触分析 (25) 点─面接触分析的步骤 (26) 点－点的接触 (35) 接触分析实例（GUI方法） (38) 非线性静态实例分析（命令流方式） (42) 接触分析 接触问题是一种高度非线性行为，需要较大的计算资源，为了进行实为有效的计算，理解问题的特性和建立合理的模型是很重要的。 接触问题存在两个较大的难点：其一，在你求解问题之前，你不知道接触区域，表面之间是接触或分开是未知的，突然变化的，这随载荷、材料、边界条件和其它因素而定；其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。
一般的接触分类 接触问题分为两种基本类型：刚体─柔体的接触，半柔体─柔体的接触，在刚体─柔体的接触问题中，接触面的一个或多个被当作刚体，（与它接触的变形体相比，有大得多的刚度），一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体─柔体的接触，许多金属成形问题归为此类接触，另一类，柔体─柔体的接触，是一种更普遍的类型，在这种情况下，两个接触体都是变形体（有近似的刚度）。 ANSYS接触能力 ANSYS支持三种接触方式：点─点，点─面，平面─面，每种接触方式使用的接触单元适用于某类问题。 为了给接触问题建模，首先必须认识到模型中的哪些部分可能会相互接触，如果相互作用的其中之一是一点，模型的对立应组元是一个结点。如果相互作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元，有限元模型通过指定的接触单元来识别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元，至于ANSTS使用的接触单元和使用它们的过程，下面分类详述。 点─点接触单元 点─点接触单元主要用于模拟点─点的接触行为，为了使用点─点的接触单元，你需要预先知道接触位置，这类接触问题只能适用于接触面之间有较小相对滑动的情况（即使在几何非线性情况下） 如果两个面上的结点一一对应，相对滑动又以忽略不计，两个面挠度（转动）保持小量，那么可以用点─点的接触单元来求解面─面的接触问题，过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。 点─面接触单元 点─面接触单元主要用于给点─面的接触行为建模，例如两根梁的相互接触。 如果通过一组结点来定义接触面，生成多个单元，那么可以通过点─面的接触单元来模拟面─面的接触问题，面即可以是刚性体也可以是柔性体，这类接触问题的一个典型例子是插头到插座里。
平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
面与面接触实例：插销拨拉问题分析 定义单元类型 Element/add/edit/delete 定义材料属性 Material Props/Material Models Structural/Linear/Elastic/Isotropic 定义材料的摩擦系数 … 建立几何模型 Modeling/Create/Volumes/Block/By Dimensions X1=Y1=0，X2=Y2=2，Z1=，Z2=
Modeling/Create/Volumes/Cylinder/By Dimensions Modeling/Operate/Booleans/Subtract/Volumes 先拾取长方体，再拾取圆柱体。 Modeling/Create/Volumes/Cylinder/By Dimensions 、 划分掠扫网格 Meshing/Size Cntrls/ManualSize/Lines/Picked Lines 拾取插销前端的水平和垂直直线，输入NDIV=3再拾取插座前端的曲线，输入NDIV=4
PlotCtrls/Style/Size and Shape，在Facets/element edge列表中选择2 facets/edge 建立接触单元 : Modeling/Create/Contact pair，弹出Contact Manager对话框，如图所示。 单击最左边的按钮，启动Contact Wizard（接触向导），如图所示。
单击Pick Target，选择目标面。 选择接触面 定义位移约束 施加对称约束，Define Loads/Apply/Structural/Displacement/Symmetric On Areas，选择对称面。 再固定插座的左侧面。 ） 设置求解选项 Analysis Type/Sol’s Control
大桥全长2996.8m，其中主桥采用跨度为101.5+188.5+580+217.5+159.5+116m的钢桁梁斜拉桥；非通航孔正桥采用6孔跨径64m预应力混凝土简支箱梁；东引桥采用16孔梁长32.6m预应力混凝土简支箱梁；跨大堤桥采用48.9+86+48.8m预应力混凝土连续箱梁；西引桥采用15孔梁长32.6m预应力混凝土简支梁及2孔梁长24.6m预应力混凝土简支梁，其中宁安线采用箱梁，阜景线采用T梁。 主桥采用103+188.5+580+217.5+159.5+117.5m两塔钢桁斜拉桥方案，全长1366m。主梁为三片主桁钢桁梁，桁间距2x14m，节间长14.5m，桁高15m。主塔为钢筋混凝土结构，塔顶高程+204.00m，塔底高程-6.00m，斜拉索为空间三索面，立面上每塔两侧共18对索，全桥216根斜拉索。所有桥墩上均设竖向和横向约束，4#塔与主梁之间设纵向水平约束，3#塔与梁间使用带限位功能的粘滞阻尼器。主梁为”N”字型桁式，横向采用三片桁结构，主桁的横向中心距各为14m，桁高15m，节间距14.5m[2]。 结构构造 主桥采用两塔钢桁斜拉桥方案，主梁为三片主桁钢桁梁，主桁上下弦杆均为箱型截面，上弦杆内高1000mm，内宽1200mm，板厚20～48mm。下弦杆内高1400mm，宽1200mm，板厚20～56mm。下弦杆顶板向桁内侧加宽700mm与整体桥面板焊接。腹杆主要采用H型截面。H型杆件宽1200mm，高720和760mm，板厚20～48mm。根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段，腹杆采用箱型截面杆件。主桁采用焊接杆件，整体节点。在节点外以高强度螺栓拼接的结构形式，上下弦杆四面等强对接拼装。H型腹杆采用插入式连接。箱型腹杆采用四面与主桁节点对拼的连接形式。主桁拼接采用M30高强螺栓。
第一章模态分析 §模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性（固有频率和振型），即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性，如塑性和接触（间隙）单元，即使定义了也将被忽略。ANSYS提供了七种模态提取方法，它们分别是子空间法、分块Lanczos法、PowerDynamics法、缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存在阻尼。后面将详细介绍模态提取方法。 §模态分析中用到的命令 模态分析使用所有其它分析类型相同的命令来建模和进行分析。同样，无论进行何种类型的分析，均可从用户图形界面（GUI）上选择等效于命令的菜单选项来建模和求解问题。 后面的“模态分析实例（命令流或批处理方式）”将给出进行该实例模态分析时要输入的命令（手工或以批处理方式运行ANSYS时）。而“模态分析实例（GUI方式）” 则给出了以从ANSYS GUI中选择菜单选项方式进行同一实例分析的步骤。（要想了解如何使用命令和GUI选项建模，请参阅<
◆前提： ◇有限元模型。 ◇已识别接触面及目标面。（*可应用自由度耦合来替代接触。） 选择目标面和接触面的准则： 1.凸面和凹面或平面接触是，选平面或凹面为目标面。2、接触的两个面网格划分有粗细的话，选粗网格所在面为目标面。3两个面刚度不同时，选择刚度大的面为目标面4如果两个面为一个高阶单元，一个为低阶单元，选低阶单元为目标面 5.如果一个面比另一个面大选大的面为目标面。 2. ◆定义接触单元及实常数
◇（刚性）目标单元—— TARGE169 TARGE170 ； ◇（柔性）接触单元—— CONTA171~CONTA172。 ***Commands*** ET,K,169 !K - 指定的单元编号 ET,K+1,172 *** **** ◇实常数——一个接触对对应同一个实常数号。 TARGE单元的实常数包括：R1、R2 —定义目标单元几何形状 CONTA单元的实常数包括： No. Name Description 1 R1 Target circle radius（刚性环半径） 2 R2 Superelement thickness（单元厚度） *3 FKN Normal penalty stiffness factor（法向接触刚度因子） *4 FTOLN Penetration tolerance factor（最大允许的穿透） *5 ICONT Initial contact closure（初始闭合因子） 6 PINB Pinball region（“Pinball”区域） *7 PMAX Upper limit of initial allowable penetration（初始穿透的最大值）*8 PMIN Lower limit of initial allowable penetration（初始穿透的最小值）*9 TAUMAX Maximum friction stress（最大的接触摩擦） *10 CNOF Contact surface offset（施加于接触面的正或负的偏移值） 11 FKOP Contact opening stiffness or contact damping *12 FKT Tangent penalty stiffness factor（切向接触刚度） 13 COHE Contact cohesion（滑动抗力粘聚力） 14 TCC Thermal contact conductance（热接触传导系数） 15 FHTG Frictional heating factor（摩擦耗散能量的热转换率） 16 SBCT Stefan-Boltzmann constant 17 RDVF Radiation view factor 18 FWGT Heat distribution weighing factor 19 ECC Electric contact conductance 20 FHEG Joule dissipation weight factor 21 FACT Static/dynamic ratio（静摩擦系数和动摩擦系数的比率） 22 DC Exponential decay coefficient（摩擦衰减系数） 23 SLTO Allowable elastic slip 24 TNOP Maximum allowable tensile contact pressure 25 TOLS Target edge extension factor 附注： +值作为比例因子，-值作为绝对值； 带*号的实常数比较重要，关乎接触分析的收敛； 一般实常数可为缺省值。