Bonded: This is the default configuration and applies to all contact regions (surfaces, solids, lines, faces, edges). If contact regions are bonded, then no sliding or separation between faces or edges is allowed.
No Separation: This contact setting is similar to the bonded case. It only applies to regions of faces (for 3-D solids) or edges (for 2-D plates). Separation of faces in contact is not allowed, but small amounts of frictionless sliding can occur along contact faces. [Not supported for Explicit Dynamics analyses.]
Frictionless: This setting models standard unilateral contact; that is, normal pressure equals zero if separation occurs. It only applies to regions of faces (for 3-D solids) or edges (for 2-D plates).
Rough: Similar to the frictionless setting, this setting models perfectly rough frictional contact where there is no sliding. It only applies to regions of faces (for 3-D solids) or edges (for 2-D plates). By default, no automatic closing of gaps is performed. This case corresponds to an infinite friction coefficient between the contacting bodies. [Not supported for Explicit Dynamics analyses.]
Frictional: In this setting, two contacting faces can carry shear stresses up to a certain magnitude across their interface before they start sliding relative to each other. It only applies to regions of faces.
Workbench中提供了5种接触类型，单从字面上很难理解这几种接触的区别，下面将帮助中关于这几个接触类型的描述翻译出来，供参考:
Bonded（绑定）:这是AWE中关于接触的默认设置。如果接触区域被设置为绑定，不允许面或线间有相对滑动或分离。可以将此区域看做被连接在一起。因为接触长度/面积是保
持不变的，所以这种接触可以用作线性求解。如果接触是从数学模型中设定的，程序将填充所有的间隙，忽略所有的初始渗透。
No Separation（不分离）:这种接触方式和绑定类似。它只适用于面。不允许接触区域的面分离，但是沿着接触面可以有小的无摩擦滑动。
Frictionless（无摩擦）:这种接触类型代表单边接触，即，如果出现分离则法向压力为零。只适用于面接触。因此，根据不同的载荷，模型间可以出现间隙。它是非线性求解，因为在载荷施加过程中接触面积可能会发生改变。假设摩擦系数为零，因此允许自由滑动。使用这种接触方式时，需注意模型约束的定义，防止出现欠约束。程序会给装配体加上弱弹簧，帮助固定模型，以得到合理的解。
Rough（粗糙的）:这种接触方式和无摩擦类似。但表现为完全的摩擦接触，即没有相对滑动。只适用于面接触。默认情况下，不自动消除间隙。这种情况相当于接触体间的摩擦系数为无穷大。
Frictional（有摩擦）:这种情况下，在发生相对滑动前，两接触面可以通过接触区域传递一定数量的剪应力。有点像胶水。模型在滑动发生前定义一个等效的剪应力，作为接触压力的一部分。一旦剪应力超过此值，两面将发生相对滑动。只适用于面接触。摩擦系数可以是任意非负值。
压电变换器的自振频率分析及详细过程 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).查看结果 模态分析的结果包括结构的频率、振型、相对应力和力等
1 引言 在ANSYS平台上，所谓子结构技术就是将一组单元用矩阵凝聚为一个单元过程的技术，切吧这个单一的矩阵单元称为超单元。在ANSYS分析中，超单元可以象其他单元类型一样使用。唯一的区别就是必须先进行结构生成分析以生成能够利用的超单元。但子结构并非在所有ANSYS模块中都能利用，目前ANSYS子结构技术可以在ANSYS/Mutiphysics，ANSYS/Mechanical和ANSYS/Structural中使用。 在ANSYS平台上，使用子结构的目的主要是为了节省机时，并且允许在比较有限的计算机设备资源的基础上求解超大规模的问题。比如进行非线性分析和带有大量重复几何结构的分析。在非线性分析中，可以将模型线性部分作成子结构，这部分的单元矩阵就不用在非线性迭代过程中重复计算。而在有重复几何结构的模型中（如有四条腿的桌子），可以对于重复的部分生成超单元，然后将它拷贝到不同的位置，这样做可以节省大量的计算时间和计算机资源。 子结构还用于模型有大转动的情况下。对于这些模型，ANSYS假定每个结构都是围绕其质心转动的。在三维情况下，子结构有三个转动自由度和三个平动自由度。在大转动模型中，用户在使用部分之前无须对子结构施加约束，因为每个子结构都是作为一个单元进行处理，是允许刚体位移的。 对于大型三维问题的分析而言，需用磁盘空间相对于一个普通计算机系统来说太庞大了，在这种情况下，用户可以通过子结构将问题分块进行分析，从而使得每一块对于计算机系统来说都是可以计算和承受的。 2 ANSYS子结构使用步骤 ANSYS子结构使用过程分为以下三个步骤： 1）ANSYS子结构生成部分 生成部分就是将普通的有限元单元凝聚为一个超单元。凝聚是通过定义一组主自由度来实现的。主自由度用于定义超单元与模型中其他单元的边界，提取模型的动力学特性。图1是一个板状构件用接触单元分析的示意。由于接触单元需要迭代计算，将板状构件形成子结构将显著地节省机时。本例中，主自由度是板与接触单元相连的自由度。 图1 子结构使用示例 2）ANSYS子结构使用部分 用部分就是将超单元与模型整体相连进行分析的部分。整个模型可以是一个超单元，也可以象上例一样是超单元与非超单元相连的。使用部分的计算只是超单元的凝聚（自由度计算仅限于主自由度）和非超单元的全部计算。
前面一篇基于Anｓys经典界面得接触分析例子做完以后,不少朋友希望了解该例子在Wｏrkbencｈ中就是如何完成得。我做了一下，与大家共享,不一定正确。毕竟这种东西,教科书上也没有,我只就是按照自己得理解在做,有错误得地方,恳请指正。 1．问题描述 一个钢销插在一个钢块中得光滑销孔中。已知钢销得半径就是０、5 ｕｎits, 长就是２、５ｕnitｓ,而钢块得宽就是4 Units，长4 Ｕnits,高为１Units,方块中得销孔半径为0、４９uｎｉtｓ,就是一个通孔。钢块与钢销得弹性模量均为3６e６，泊松比为0、3、由于钢销得直径比销孔得直径要大,所以它们之间就是过盈配合。现在要对该问题进行两个载荷步得仿真。 (1)要得到过盈配合得应力。 (2)要求当把钢销从方块中拔出时,应力，接触压力及约束力。 2．问题分析 由于该问题关于两个坐标面对称,因此只需要取出四分之一进行分析即可。 进行该分析,需要两个载荷步:
第一个载荷步,过盈配合。求解没有附加位移约束得问题,钢销由于它得几何尺寸被销孔所约束,由于有过盈配合,因而产生了应力。 第二个载荷步,拔出分析。往外拉动钢销1、7 ｕｎｉts,对于耦合节点上使用位移条件。打开自动时间步长以保证求解收敛。在后处理中每10个载荷子步读一个结果。 本篇只谈第一个载荷步得计算。 3．生成几何体 上述问题就是ＡＮSYS自带得一个例子。对于几何体,它已经编制了生成几何体得命令流文件。所以,我们首先用经典界面打开该命令流文件，运行之以生成四分之一几何体;然后导出为一个IGS文件,再退出经典界面,接着再到WORKＢＥＮＣH中，打开该IGS文件进行操作。 (3、1)首先打开ＡNSYＳAPＤL1４、5、 (３、2）然后读入已经做好得几何体。从【工具菜单】-->【＞【Ｒead Iｎｐut Frｏm】打开导入文件对话框
手把手教你用ANSYS workbench 本文的目的主要是帮助那些没有接触过ansys workbench的人快速上手使用这个软件。在本文里将展示ansys workbe nch如何从一片空白起步，建立几何模型、划分网格、设置约束和边界条件、进行求解计算，以及在后处理中运行疲劳分析模块，得到估计寿命的全过程。 一、建立算例 打开ansys workbench这时还是一片空白。 ■A Un■$曲虑日Project - Wor^L-bemdi FI E Vievi Took Units EKlhenMrs Hep 口百］牙.匾1丿狂存*■::_____________ 4J Import-■■ ?b RBConn^dt | 半］Project Lbd盘B Project g pp^iijT 咗nifint 世Eiqen/alue Ekxkfing Q Elqenwlue Bucktig ［samcef) 醪Flwtnc 闵E^pict Cynannics ? Fluid F I M -M UN Mud凶『山山理］ ◎Hud Ftaw - Estrusoi (PdyflEMiJ ? Fluid Flow (CFX) 也rlud Flow :FkirflL) Q Hud How (Pdvftouf) I朗Hermoinic IResporiSB 営H>d,qdyr>amic DiFFractlon I岂？H^drcclj/riarw Resrwnw 讐 JCEnjina = 逝MocW 爲Moda (阳AQU5) fjy Muds 口■ ii』) 肚| H^ndorn wbracior 迦| Spedtium Riyid D/ruriL^ 国StStIC ^truchjral 冒Static Structural 卜对Static■Strucbj-cl (5aTiccF) 1 5Zac\-5taZ Wrnml D Ihemnal 0 5tcady-5Uts Ihcmal (Sancd7) 密Thnrra^-FlPirrrir 电j Tlroughlkw ◎Il i oughfki^ ^DiaJcGcrO innsflnr strudturAi 回7rans?n: Structural (ABiQUS) 褪Tr slismL 5trudtural (Stfncsf) A 怕Ment rhenr^l 首先我们要清楚自己要计算的算例的分析类型，一般对于结构力学领域，有 静态分析(Static Structural)、动态分析(Rigid Dynamics)、模态分析(Modal)。
①下面对非对称行为接触表面的正确选择给出选择指导: –如果一凸的表面要和一平面或凹面接触，应该选取平面或凹面为目标面. –如果一个表面有粗糙的网格而另一个表面网格细密，则应选择粗糙网格表面为目标面. –如果一个表面比另一个表面硬，则硬表面应为目标面. –如果一个表面为高阶而另一个为低阶，则低阶表面应为目标面. –如果一个表面大于另一个表面，则大的表面应为目标面. ②法向刚度WB-Mechanical系统默认自动设定。 –用户可以输入“法向刚度因子Normal Stiffness Factor” (FKN) 它是计算刚度代码的乘子.因子越小，接触刚度就越小。 ?默认 FKN =10 (对于绑定和不分离的接触) ?默认 FKN=10(其他形式接触) 默认 FKN1.0 (其他形式接触) ?接触问题法向刚度选择一般准则: –体积为主的问题: 用“Program Controlled”或手动输入“Normal Stiffness Factor”为“1” –弯曲为主的问题: 手动输入“Normal Stiffness Factor”为“0.01”到“0.1”之 间的数值。 -在大变形问题的无摩擦或摩擦接触中建议使用“Augmented Lagrange” 法向接触刚度 knormal是影响精度和收敛行为最重要的参数. –刚度越大，结果越精确，收敛变得越困难. –如果接触刚度太大，模型会振动，接触面会相互弹开。 - 其中update stifness 设置可以控制计算收敛与否。
③ -刚度增加, 渗透减少，而最大压力增加. 并且通常会有更多的迭代和更长运行时间 ④ 不管使用了何种接触行为 (对称或反对称), 模型的变形和等效应力本质 是相同的. 对称行为可以提高收敛. 但对称接触结果不容易解释，为接触面与目标面结果的平均值。 0.0032902 0.0033033 0.0033052 0.0033055 0.0033053 565.05Mp a 774.12Mp a 811.34Mp a 816.26Mp a 812.78Mp a 0.011864 0.0016253 0.0017035 0.000017138 0.000019984 17 17 20 24 57
1、 ANSYS12.1 Workbench界面相关分析系统和组件说明 【Analysis Systems】分析系统【Component Systems】组件系统【CustomSystems】自定义系统【Design Exploration】设计优化 分析类型说明 Electric (ANSYS) ANSYS电场分析 Explicit Dynamics (ANSYS) ANSYS显式动力学分析 Fluid Flow (CFX) CFX流体分析 Fluid Flow (Fluent) FLUENT流体分析 Hamonic Response (ANSYS) ANSYS谐响应分析 Linear Buckling (ANSYS) ANSYS线性屈曲Magnetostatic (ANSYS) ANSYS静磁场分析 Modal (ANSYS) ANSYS模态分析 Random Vibration (ANSYS) ANSYS随机振动分析Response Spectrum (ANSYS) ANSYS响应谱分析 Shape Optimization (ANSYS) ANSYS形状优化分析 Static Structural (ANSYS) ANSYS结构静力分析 Steady-State Thermal (ANSYS) ANSYS稳态热分析 Thermal-Electric (ANSYS) ANSYS热电耦合分析Transient Structural(ANSYS) ANSYS结构瞬态分析Transient Structural(MBD) MBD 多体结构动力分析Transient Thermal(ANSYS) ANSYS瞬态热分析 组件类型说明 AUTODYN AUTODYN非线性显式动力分析BladeGen 涡轮机械叶片设计工具 CFX CFX高端流体分析工具
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)
1概述 接触问题是一种高度非线性行为，需要较大的计算资源，为了进行实为有效的计算，理解问题的特性和建立合理的模型是很重要的。 接触问题存在两个较大的难点：其一，在你求解问题之前，你不知道接触区域，表面之间是接触或分开是未知的，突然变化的，这随载荷、材料、边界条件和其它因素而定；其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。 一般的接触分类 接触问题分为两种基本类型：刚体─柔体的接触，半柔体─柔体的接触。 （1）刚-柔接触 在刚体─柔体的接触问题中，接触面的一个或多个被当作刚体，（与它接触的变形体相比，有大得多的刚度），一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体─柔体的接触，许多金属成形问题归为此类接触。 （2）柔-柔接触 柔体─柔体的接触，是一种更普遍的类型，在这种情况下，两个接触体都是变形体（有近似的刚度）。 2ANSYS接触能力 ANSYS支持三种接触方式：点─点，点─面，平面─面，每种接触方式使用的接触单元适用于某类问题。为了给接触问题建模，首先必须认识到模型中的哪些部分可能会相互接触，如果相互作用的其中之一是一点，模型的对立应组元是一个结点。如果相互作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元，有限元模型通过指定的接触单元来识别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元，至于ANSYS使用的接触单元和使用它们的过程，下面分类详述。 2.1点─点接触单元 点─点接触单元主要用于模拟点─点的接触行为，为了使用点─点的接触单元，你需要预先知道接触位置，这类接触问题只能适用于接触面之间有较小相对滑动的情况（即使在几何非线性情况下） 如果两个面上的结点一一对应，相对滑动又以忽略不计，两个面挠度（转动）保持小量，那么可以用点─点的接触单元来求解面─面的接触问题，过盈装配问题是一个用点─点的接触单元来模拟面─面的接触问题的典型例子。
参考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使用的接触单元和使用它们的过程，下面分类详述。 点─点接触单元 点─点接触单元主要用于模拟点─点的接触行为，为了使用点─点的接触单元，你需要预先知道接触位置，这类接触问题只能适用于接触面之间有较小相对滑动的情况（即使在几何非线性情况下） 如果两个面上的结点一一对应，相对滑动又以忽略不计，两个面挠度（转动）保持小量，那么可以用点─点的接触单元来求解面─面的接触问题，过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。 点─面接触单元 点─面接触单元主要用于给点─面的接触行为建模，例如两根梁的相互接触。 如果通过一组结点来定义接触面，生成多个单元，那么可以通过点─面的接触单元来模拟面─面的接触问题，面即可以是刚性体也可以是柔性体，这类接触问题的一个典型例子是插头到插座里。
面与面接触实例：插销拨拉问题分析 定义单元类型 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
◆前提： ◇有限元模型。 ◇已识别接触面及目标面。（*可应用自由度耦合来替代接触。） 选择目标面和接触面的准则： 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 附注： +值作为比例因子，-值作为绝对值； 带*号的实常数比较重要，关乎接触分析的收敛； 一般实常数可为缺省值。
）间的球形界面的压力形貌。
上机实验报告: 软件版本：ANSYS workbench 19.2 1.主要分析过程及注意事项 分析过程： ●打开workbench，从左侧的“analysis system”中拖入“static structural”到中间空白区域 ●由于材料已经是默认的结构钢，所以我们不用修改，但是单位和它的显示模式我们要改 成像下图中的(Tonne,mm,…)和“display values in project units”。 ●在geometry中导入“ball-socket.x_t ”之前，先在右边的属性栏里，找到analysis type， 将3D改为2D，改完之后再导入“ball-socket.x_t ”。
●双击model进去“mechanical”，选中Geometry，在Definition中把2D Behavior改为 Axisymmetric。同时检查工作单位制是否是Metric (mm,kg,N,s,mV,mA) ●选中“contacts”，插入“Frictional” Frictional Coefficient设为0.4，behavior改为auto asymmetric（自动非对称），formulation改为augmented lagrange（后面的试验结果表明，formulation设为program controlled，结果都一样）
●在analysis setting里把Large Deflection改为ON ●鼠标选中mesh，我们可以在下面的element size 改变网格大小，本上机实验中会分别试验 1.0mm和0.5mm，修改完后右键generate mesh可观看效果 ●选中static structural，插入fixed support ，选中socket的上边线，并apply，然后在插入loads 里的force，这时选择ball的下边线，并apply，在define by里选择component，并在y方向上输入-1000。
ANSYS 中如何使用接触向导定义接触对 在ANSYS 中定义接触通常有两种方法： 1. 用户自己手工创建接触单元和目标单元。这种方法，在定义接触和目标单元时还比较简单，但是在设置或修改单元属性和定义实常数时却比较复杂。需要用户对接触有较深刻的理解和通过实践积累丰富的经验。 2. 使用接触管理器中的接触向导定义接触对：使用接触管理器 (接触向导) 定义接触对(即接触单元和目标单元) 时，可以定义除了点-点接触以外的各种接触类型；它可以自动生成接触单元和目标单元，并提供了一组默认的单元属性和实常数值。使用这些默认的设置，加上适当的求解设置，对于多数接触问题都能够获得收敛的结果。而且，如果使用默认设置时，计算不收敛或对结果不太满意，也可以通过接触管理器(接触向导) 对单元属性和实常数方便的进行修改和调整。 因此，我们推荐，在可能的情况下，尽量使用接触管理器(接触向导) 来定义接触。本文将通过一个实例介绍接触管理器的基本使用方法。 所使用的例子如下： 两块平板，中间夹一个圆球。上面平板的上表面承受压力，分析模型的变形和应力随压力的变化。 两块平板，尺寸都是(100*100*20)，相距100。中间夹一个半径50 的圆球。两个平板分别与圆球的上下边缘接触。尺寸单位为mm。几何模型如图1。
图 1 中，为了能够划分映射网格，分别对体积进行了切割材料属性为：两块平板： E = 201000 Mpa；μ= 0.3 圆球： E = 70100 Mpa；μ= 0.33 接下来对各个Volumes 划分网格，单元类型采用solid186 (20 节点六面体)，单元边长统一取 6 mm。网格划分结果如图 2 所示：
前言 WokBench 是众所周知的好东西，以下是自己琢磨的一个小应用，肯定有不对的地方， 欢迎指出，便于大家共同提高。 问题描述 这是一个塑料小卡扣的例子，主要想使用WorkBench 了解在使用中，塑料件的变形是否足够。模型是用ProE 制作的，为了简化，只切取了关于变形的部分，如下图： 其中蓝色的部分是活动的，只有一个方向的运动，红色的部分是固定的。 大体的尺寸如下，单位是毫米：
注意：在模型中，蓝色和红色部件的距离要控制好（这是由ProE 中，模型装配关系 决定的），如果太近，软件将自动计算出一个接触区域，但对于这个例子，还需要手 动扩大接触区域。如果距离太远，在手动设置Pinball 类型的接触区域时，Pinball 的 半径要设得很大，可能导致无法计算。请参考上面的尺寸图纸调节两个部件之间的距 离。 之后，设置接触面（2、3）：需要将两个部件在运动过程中，会接触的地方一一标出， 千万不要加无用的面。 将Pinball Region 设置为Radius 方式（4），并将Radius 设置一个合适的值（5），本例设置了3 毫米（如图，会形成一个蓝色的大圆球），求解的时候软件会使用这个PinBall 自动探测接触。 还需要将接触方式设置为无摩擦的（6）。 最后将接触面计算方式设置为Adjust To Touch（7）。也可以尝试其他的方式，不过对 于这个仅研究红色部件变形的例子就无所谓了。
关于单元格 WorkBench 中可以不自行划分单元格（在解算的时候，如果没有手动的设置，软件就会先自动划分），软件帮你自动产生。如果你的其他设置正确，即便是这个自动的值也能很精确了。 添加分析 这个分析用静力学就可以了（1）。 之后要设置Analysis Setting（2）。将Nuber Of Step 设置为2（3）。 注意： 1）蓝色部件在运动的过程中，先压迫红色部件，再逐渐松开，因此必须将这个过 程至少分解为至少两个阶段（阶段指“Step”）。 2）对于一个阶段而言，Ansys 求解时，会先考察它的开始和结束两个点的状态。
ANSYS WORKBENCH 疲劳分析指南 第一章简介 1.1 疲劳概述 结构失效的一个常见原因是疲劳，其造成破坏与重复加载有关。疲劳通常分为两类：高周疲劳是当载荷的循环（重复）次数高(如1e4 -1e9)的情况下产生的。因此，应力通常比材料的极限强度低，应力疲劳（Stress-based）用于高周疲劳；低周疲劳是在循环次数相对较低时发生的。塑性变形常常伴随低周疲劳，其阐明了短疲劳寿命。一般认为应变疲劳（strain-based）应该用于低周疲劳计算。 在设计仿真中，疲劳模块拓展程序（Fatigue Module add-on）采用的是基于应力疲劳（stress-based）理论，它适用于高周疲劳。接下来，我们将对基于应力疲 劳理论的处理方法进行讨论。 1.2 恒定振幅载荷 在前面曾提到，疲劳是由于重复加载引起： 当最大和最小的应力水平恒定时，称为恒定振幅载荷，我们将针对这种最简 单的形式，首先进行讨论。 否则，则称为变化振幅或非恒定振幅载荷。 1.3 成比例载荷 载荷可以是比例载荷，也可以非比例载荷： 比例载荷，是指主应力的比例是恒定的，并且主应力的削减不随时间变化， 这实质意味着由于载荷的增加或反作用的造成的响应很容易得到计算。 相反，非比例载荷没有隐含各应力之间相互的关系，典型情况包括： σ1/σ2=constant 在两个不同载荷工况间的交替变化； 交变载荷叠加在静载荷上； 非线性边界条件。 1.4 应力定义 考虑在最大最小应力值σmin和σmax作用下的比例载荷、恒定振幅的情况： 应力范围Δσ定义为(σmax-σmin) 平均应力σm定义为(σmax+σmin)/2 应力幅或交变应力σａ是Δσ/2 应力比R是σmin/σmax 当施加的是大小相等且方向相反的载荷时，发生的是对称循环载荷。这就是 σm=0，R=-1的情况。 当施加载荷后又撤除该载荷，将发生脉动循环载荷。这就是σm=σmax/2，R=0的情况。 1.5 应力-寿命曲线 载荷与疲劳失效的关系，采用的是应力-寿命曲线或S-N曲线来表示： （1）若某一部件在承受循环载荷, 经过一定的循环次数后,该部件裂纹或破坏将会发展，而且有可能导致失效； （2）如果同个部件作用在更高的载荷下,导致失效的载荷循环次数将减少； （3）应力-寿命曲线或S-N曲线,展示出应力幅与失效循环次数的关系。
Workbench -Mechanical Introduction Introduction 作业3.1 31 接触控制
作业3.1 –目标 Workshop Supplement ?作业3.1调查了一个简单组件的接触行为。目的是为了说明由于不适当接触导致的刚体运动是怎么产生的。 ?问题描述： 问题描述 –模型从一个简单Parasolid组件文件获得 –我们的目标是在组件的各部件中建立接触，查看非对称加载对结果有何影响 我们的目标是在组件的各部件中建接触，查看非对称加载对结果有何影响
作业3.1 –假设 Workshop Supplement ?假设arm shaft 和side plate上的孔间的摩擦忽略不计，同样arm shaft 和stop shaft 之间的接触也忽略不计。最后假设stop shaft固定在两个side plate之间。 之间 Arm Shaft Side Plate Side Plate p Stop Shaft
作业3.1 –Project Schematic Workshop Supplement ?打开Project page（项目页） ?通过“Units” 菜单确定： –Project单位设置为“US Customary (lbm, in, s, F, A, lbf, V). –选择“Display Values in Project Units”
. . .作业3.1 –Project Schematic Workshop Supplement 1.在Toolbox（工具箱）中双击 Static Structural建立新的分析系 统 1. 2.Geometry上点击鼠标右键选择 2在 Import Geometry导入 2. Contact_Arm.x_t文件
ANSYS高级接触问题 第一章接触问题概述 在工程中会遇到大量的接触问题，如齿轮的啮合、法兰联接、机电轴承接触、卡头与卡座、密封、板成形、冲击等等。接触是典型的状态非线性问题，它是一种高度非线性行为。接触例子如图1： 分析中常常需要确定两个或多个相互接触物体的位移、接触区域的大小和接触面上的应力分布。 接触分析存在两大难点： 在求解之前，你不知道接触区域、表面之间是接触或分开是未知的，表面之间突然接触或突然不接触会导致系统刚度的突然变化。 大多数接触问题需要计算摩擦。摩擦是与路径有关的现象，摩擦响应还可能是杂乱的，使问题求解难以收敛。 1.1 接触分类 1.1.1 刚－柔 一个表面是完全刚性的—除刚体运动外无应变、应力和变形，另一表面为软材料构成是可变形的。 只在一个表面特别刚硬并且不关心刚硬物体的应力时有效。 1.1.2 柔－柔
两个接触体都可以变形。 1.2 接触单元 ANSYS采用接触单元来模拟接触问题：跟踪接触位置；保证接触协调性（防止接触表面相互穿透）；在接触表面之间传递接触应力（正压力和摩擦）。 接触单元就是覆盖在分析模型接触面上的一层单元。在ANSYS中可以采用三种不同的单元来模拟接触：面—面接触单元；点—面接触单元；点—点接触单元。 不同的单元类型具有完全不同的单元特性和分析过程。 1.2.1 面—面接触单元 用于任意形状的两个表面接触，不必事先知道接触的准确位置；两个面可以具有不同的网格；支持大的相对滑动；支持大应变和大转动。例如：面一面接触可以模拟金属成型，如轧制过程。 1.2.2 点—面接触单元 用于某一点和任意形状的面的接触，可使用多个点－面接触单元模拟棱边和面的接触；不必事先知道接触的准确位置；两个面可以具有不同的网格；支持大的相对滑动；支持大应变和大转动。例：点面接触可以模拟棱边和面之间的接触。 1.2.3 点－点接触单元 用于模拟单点和另一个确定点之间的接触。建立模型时必须事先知道确切的接触位置；多个点－点接触单元可以模拟两个具有多个单元表面间的接触；每个