文档库 最新最全的文档下载
当前位置:文档库 › ANSYS中的接触分析教程

ANSYS中的接触分析教程

ANSYS中的接触分析教程
ANSYS中的接触分析教程

一般的接触分类 (2)

ANSYS接触能力 (2)

点─点接触单元 2

点─面接触单元 2

面─面的接触单元 3

执行接触分析 (3)

面─面的接触分析 4

接触分析的步骤: 4

步骤1:建立模型,并划分网格 (4)

步骤2:识别接触对 (4)

步骤2:指定接触面和目标面 4

步骤4:定义刚性目标面 (5)

步骤5:定义柔性体的接触面 (8)

步骤6:设置实常数和单元关键字 (9)

步骤7:控制刚体目标的运动 (19)

步骤8:给变形体单元加必要的边界条件 (20)

步骤9:定义求解和载荷步选项

20

第十步:检查结果 (21)

点─面接触分析 (23)

点─面接触分析的步骤 (24)

点-点的接触 (32)

接触分析实例(GUI方法) (34)

非线性静态实例分析(命令流方式) (37)

接触分析

接触问题是一种高度非线性行为,需要较大的计算资源,为了进

行实为有效的计算,理解问题的特性和建立合理的模型是很重要的。

接触问题存在两个较大的难点:其一,在你求解问题之前,你不知道接触区域,表面之间是接触或分开是未知的,突然变化的,这随载荷、材料、边界条件和其它因素而定;其二,大多的接触问题需要计算摩擦,有几种摩擦和模型供你挑选,它们都是非线性的,摩擦使问题的收敛性变得困难。

一般的接触分类

接触问题分为两种基本类型:刚体─柔体的接触,半柔体─柔体的接触,在刚体─柔体的接触问题中,接触面的一个或多个被当作刚体,(与它接触的变形体相比,有大得多的刚度),一般情况下,一种软材料和一种硬材料接触时,问题可以被假定为刚体─柔体的接触,许多金属成形问题归为此类接触,另一类,柔体─柔体的接触,是一种更普遍的类型,在这种情况下,两个接触体都是变形体(有近似的刚度)。

ANSYS接触能力

ANSYS支持三种接触方式:点─点,点─面,平面─面,每种接触方式使用的接触单元适用于某类问题。

为了给接触问题建模,首先必须认识到模型中的哪些部分可能会相互接触,如果相互作用的其中之一是一点,模型的对立应组元是一个结点。如果相互作用的其中之一是一个面,模型的对应组元是单元,例如梁单元,壳单元或实体单元,有限元模型通过指定的接触单元来识别可能的接触匹对,接触单元是覆盖在分析模型接触面之上的一层单元,至于ANSTS使用的接触单元和使用它们的过程,下面分类详述。

点─点接触单元

点─点接触单元主要用于模拟点─点的接触行为,为了使用点─点的接触单元,你需要预先知道接触位置,这类接触问题只能适用于接触面之间有较小相对滑动的情况(即使在几何非线性情况下)如果两个面上的结点一一对应,相对滑动又以忽略不计,两个面挠度(转动)保持小量,那么可以用点─点的接触单元来求解面─面的接触问题,过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。

点─面接触单元

点─面接触单元主要用于给点─面的接触行为建模,例如两根梁的相互接触。

如果通过一组结点来定义接触面,生成多个单元,那么可以通过点─面的接触单元来模拟面─面的接触问题,面即可以是刚性体也可

以是柔性体,这类接触问题的一个典型例子是插头到插座里。

使用这类接触单元,不需要预先知道确切的接触位置,接触面之间也不需要保持一致的网格,并且允许有大的变形和大的相对滑动。

Contact48和Contact49都是点─面的接触单元,Contact26用来模拟柔性点─刚性面的接触,对有不连续的刚性面的问题,不推荐采用Contact26因为可能导致接触的丢失,在这种情况下,Contact48通过使用伪单元算法能提供较好的建模能力。

面─面的接触单元

ANSYS支持刚体─柔体的面─面的接触单元,刚性面被当作“目标”面,分别用Targe169和Targe170来模拟2─D和3—D的“目标”面,柔性体的表面被当作“接触”面,用Conta171,Conta172,Conta173,Conta174来模拟。一个目标单元和一个接单元叫作一个“接触对”程序通过一个共享的实常号来识别“接触对”,为了建立一个“接触对”给目标单元和接触单元指定相同的实常的号。

与点─面接触单元相比,面─面接触单元有好几项优点,

·支持低阶和高阶单元

·支持有大滑动和摩擦的大变形,协调刚度阵计算,单元提法不对称刚度阵的选项。

·提供工程目的采用的更好的接触结果,例如法向压力和摩擦应力。

·没有刚体表面形状的限制,刚体表面的光滑性不是必须允许有自然的或网格离散引起的表面不连续。

·与点─面接触单元比,需要较多的接触单元,因而造成需要较小的磁盘空间和CPU时间。

·允许多种建模控制,例如:

·绑定接触

·渐变初始渗透

·目标面自动移动到补始接触

·平移接触面(老虎梁和单元的厚度)

·支持死活单元

使用这些单元,能模拟直线(面)和曲线(面),通常用简单的几何形状例如圆、抛物线、球、圆锥、圆柱采模拟曲面,更复杂的刚体形状能使用特殊的前处理技巧来建模。

执行接触分析

不同的接触分析类型有不同的过程,下面分别讨论

面─面的接触分析

在涉及到两个边界的接触问题中,很自然把一个边界作为“目标”面而把另一个作为“接触”面,对刚体─柔体的接触,“目标”面总是刚性的,“接触”面总是柔性面,这两个面合起来叫作“接触对”使用Targe169和Conta171或Conta172来定义2-D接触对,使用Targe170和Conta173或Conta174来定义3-D接触对,程序通过相同的实常收号来识别“接触对”。

接触分析的步骤:

执行一个典型的面─面接触分析的基本步骤列示如下:

1.建立模型,并划分网格

2.识别接触对

3.定义刚性目标面

4.定义柔性接触面

5.设置单元关键字和实常的

6.定义/控制刚性目标面的运动

7.给定必须的边界条件

8.定义求解选项和载荷步

9.求解接触问题

10.查看结果

步骤1:建立模型,并划分网格

在这一步中,你需要建立代表接触体几何形状的实体模型。与其它分析过程一样,设置单元类型,实常的,材料特性。用恰当的单元类型给接触体划分网格。

命令:AMESH

VMESH

GUI:Main Menu>Preprocessor>mesh>Mapped>3 or4 Sided

Main Menu>Pneprocessor>mesh>mapped>4 or 6 sided

步骤二:识别接触对

你必须认识到,模型在变形期间哪些地方可能发生接触,一是你已经识别出潜在的接触面,你应该通过目标单元和接触单元来定义它们,目标和接触单元跟踪变形阶段的运动,构成一个接触对的目标单元和接触单元通过共享的实常号联系起来。

接触环(区域)可以任意定义,然而为了更有效的进行计算(主要指CPU时间)你可能想定义更小的局部化的接触环,但能保证它足以描述所需要的接触行为,不同的接触对必须通过不同的实常数号来定义(即使实常数号没有变化)。

由于几何模型和潜在变形的多样形,有时候一个接触面的同一区

域可能和多个目标面产生接触关系。在这种情况下,应该定义多个接触对(使用多组覆盖层接触单元)。每个接触对有不同的实常数号。

步骤三:定义刚性目标面

刚性目标面可能是2—D的或3─D的。在2—D情况下,刚性目标面的形状可以通过一系列直线、圆弧和抛物线来描述,所有这些都可以用TAPGE169来表示。另外,可以使用它们的任意组合来描述复杂的目标面。

在3—D情况下,目标面的形状可以通过三角面,圆柱面,圆锥面和球面来推述,所有这些都可以用TAPGE170来表示,对于一个复杂的,任意形状的目标面,应该使用三角面来给它建模。

控制结点(Pilot)

刚性目标面可能会和“pilot结点“联系起来,它实际上是一个只有一个结点的单元,通过这个结点的运动可以控制整个目标面的运动,因此可以把pilot结点作为刚性目标的控制器。整个目标面的受力和转动情况可以通过pilot结点表示出来,“pilot结点”可能是目标单元中的一个结点,也可能是一个任意位置的结点,只有当需要转动或力矩载荷时,“pilot结点”的位置才是重要的,如果你定义了“pilot 结点”ANSYS程序只在“pilot结点”上检查边界条件,而忽略其它结点上的任何约束。

对于圆、圆柱、圆锥、和球的基本图段,ANSYS总是使用条一个结点作为“pilot结点”

基本原型

你能够使用基本几形状来模拟目标面,例如:“圆、圆柱、圆锥、球。直线、抛物线、弧线、和三角形不被允许、虽然你不能把这些基本原型彼此合在一起,或者是把它们和其它的目标形状合在一起以便形成一个同一实常数号的复杂目标面。但你可以给每个基本原型指定它自己的实常的号。

单元类型和实常数

在生成目标单元之前,首先必须定义单元类型(TARG169或TARG170)。

命令:ET

GUI:main menu>preprocessor>Element Type> Add/Edit/Delete

随后必须设置目标单元的实常数。

命令:Real

GUI:main menn>preprocessor>real constants

对TARGE169和TARGE170仅需设置实常数R1和R2,而只有在使用直接生成法建立目标单元时,才需要从为指定实常数R1、R2,另外除了直接生成法,你也可以使用ANSYS网格划分工具生成目标

单元,下面解释这两种方法。

使用直接生成法建立刚性目标单元

为了直接生成目标单元,使用下面的命令和菜单路径。

命令:TSHAP

GUI:main menu>preprocessor>modeling-create>Elements>Elem Attributes

随后指定单元形状,可能的形状有:

·straight line (2D)

·parabola (2-D)

·clockwise arc(2-D)

·counterclokwise arc (2-D)

·circle(2-D)

·Triangle (3-D)

·Cylinder (3-D)

·Cone (3-D)

·Sphere (3-D)

·Pilot node (2-D和3-D)

一旦你指定目标单元形状,所有以后生成的单元都将保持这个形状,除非你指定另外一种形状。

然后你就可以使用标准的ANSYS直接生成技术生成结点和单元。

命令:N

E

GUI:main menu>pnoprocessor> modeling- create> nodes

main menu>pnoprocessor> modeling- create>Elements 在建立单元之后,你可以通过列示单元来验证单元形状

命令:ELIST

GUI:utility menu>list>Elements>Nodes+Attributes

使用ANSYS网格划分工具生成刚性目标单元

你也可以使用标准的ANSYS网格划分功能让程序自动地生成目标单元,ANSYS程序将会以实体模型为基础生成合适的目标单元形状而忽略TSHAP命令的选项。

为了生成一个“PILOT结点”使用下面的命令或GUI路径:

命令:Kmesh

GUI:main menu>proprocessor>meshing-mesh>keypoints

注意:KMESH总是生成“PILOT结点”

为了生成一个2─D目标单元,使用下面的命令和GUI路径:ANSYS在每条直线上生成一条单一的线,在样条曲线上生成抛物线部分,在每条圆弧和倒角上生成圆弧部分,如果所有的圆弧形成

一个封闭的圆,ANSYS生成一个单一的圆段。

命令:LMESH

GUI:main menu>pneprocessor>mesling-mesh>lines

为了生成3─D的目标单元,使用下面的命令或GUI路径。

如果实体模型的表面部分形成了一个完整的球,圆柱或圆锥,那么ANSYS程序自动生成一个基本的3─D目标单元,因为生成较少的单元,从而使你分析计算更有效率,对任意形状的表面,应该使用Amesh命令来生成目标单元,在这种情况下,网格形状的质量不是重要的,而目标单元的形状是否能完成好的模拟刚性面的表面几何形状显得更重要。

命令:AMESH

GUI:main menu>preprocessor>-meshing-mesh>Area

ANSYS在所有可能的面上推荐使用三角形的映射网格划分,如果在表面的边界上没有曲率,则在网格划分时,指定那条边界分为一分,下面的命令或GUI路径将尽可能的生成一个映射网格(如果不能进行映射,它将生成自由网格)

命令:MSHKFY,2

GUI:main menu>preprocessor>-meshling-mesh>-Ares-Target Surf

建模和网格划分的注意点:

一个目标面可能由两个或多个面断的区域组成,你应该尽可能地通过定义多个目标面来使接触区域局部比(每个目标面有一个不同的实常数号)刚性目标面上由的离散能足够指述出目标面的形状,过粗的网格离散可能导致收敛问题。如果刚性面有一个实的凸角,求解大的滑动问题时很难获得收敛结果,为了避免这些建模问题,在实体模型上,使用线或面的倒角来使尖角光滑比,或者在曲率突然变化的区域使用更细的网格。

注意:不能使用镜面对称技术(ARSYSM,LSYMM)来映射圆、圆柱、圆锥或球面到对称平面的另一边,因为每个实常数的设置不能同时赋给多个基本原型段。

检验目标面的接触方向。

目标面的结点号顺序是重要的,因为它定义了接触主向,对2─D接触问题,当沿着目标线从第一个结点移向第二个结点时,变形体的接触单元必须位于

目标面的右边。

对3─D接触问题,目标三角形单元号应该使刚性面的外法线方向指向接触面,外法线通过右手原则来定义

为了检查法线方向,显示单元坐标系

命令:/PSYMS,ESYS,1

GUI:Utility menu>plotctrls>symbols

如果单元法向不指向接触面,选择单元反转表面的法向的方向。

命令:ESURF,,REVE

GUI:main menu>preprocossor>create>Element>on free surf

步骤4:定义柔性体的接触面

为了定义柔性体的接触面,必须使用接触单元CONFA171或CONFA172(对2─D)或CONTA173或CONTA174(对3─D)来定义表面

程序通过组成变形体表面的接触单元来定义接触表面,接触单元与下面覆盖的变形体单元有同样的几何特性,接触单元与下面覆盖的变形体单元必须处于同一阶次(低阶或高阶)下面的变形体单元可能是实体单元、壳单元、梁单元或超单元,接触面可能壳或梁单元任何一边。

与目标面单元一样,你必须定义接触面的单元类型,然后选择正确的实常数号(实常数号必须与它对应目标的实常数号相同)最后生成接触单元。

单元类型:

下面简单描述四种类型的接触单元

CONTA171:这是一种2─D,2个结点的低附线单元,可能位于2─D实体,壳或梁单元的表面

CONTA172:这是一个2─D的,3结点的高阶抛物线形单元,可能位于有中结点的2─D实体或梁单元的表面

CONTA173:这是一个3─D的,4结点的低阶四边形单元可能位于3─D实体或壳单元的表面,它可能褪化成一个结点的三角形单元。

CONTA174:这是一个3─D,8结点的高阶四边形单元,可能位于有中结点的3─D实体或壳单元的表面,它可能褪化成6结点的三角形单元。

不能在高阶柔性体单元的表面上分成低阶接触单元,反之也不行,不能在高阶接触单元上消去中结点。

命令:ET

GUI:main menu>preprocessor>Element type>Add/Edit/Delete 实常数和材料特性

在定义了单元类型之后,需要选择正确的实常数的设置,每个接触对的接触面和目标面必须有相同的实常数号,而每个接触对必须有它自己不同的实常数号。

ANSYS使用下面柔性体单元的材料特性来计算一个合适的接触

(或罚)刚度,如果下面的单元是一个超单元。接触单元的材料的设置必须与超单元形成时的原始结构单元相同,生成接触单元。

我们既可以通过直接生成法生成接触单元,也可以在柔性体单元的外表面上自动生成接触单元,我们推荐采用自动生成法,这种方法更为简单和可靠。

可以通过下面三个步骤来自动生成接触单元

1、选择结点

选择已划分网格的柔性体表面的结果,如果你确定某一部分结点永远不会接触到目标面,你可以忽略它以便减少计算时间,然而,你必须保证设有漏掉可能会接触到目标面的结点。

命令:NSEL

GUI:main menu>preprocessor>create>Element>on>free surf

2、生成接触单元

命令:ESURF

GUI:main menu>preprocessor>create>Element>on free surf

如果接触单元是附在已用实体单元划分网格的面或体上,程序会自动决定接触计算所需的外法向,如果下面的单元是梁或壳单元,则必须指明哪个表面(上表面或下表面)是接触面

命令:ESURF,TOP OR BOTIOM

GUI:main menu>preprocessor>create>Element>on free surf

使用上表面生成接触单元,则它们的外法向与梁或壳单元的法向相同,使用下表面生成接触单元,则它们的外法向与梁或壳单元的法向相反,如果下面的单元是实体单元,则TOP或BOTTOM选项不起作用

3、检查接触单元外法线的方向,当程序进行是否接触的检查时,接触面的外法线方向是重要的,对3─D单元,按结点程序号以右手定则来决定单元的外法向,接面的外法向应该指向目标面,否则,在开始分析计算时,程序可能会认为有面的过度渗透而很难找到初始解。在此情况下,程序一般会立即停止执行,你可以检查单元外法线方向是否正确。

命令:/PSYMB

GUI:Utility menu>plotctrls>symbols

当发现单元的外法线方向不正确时,必须通过倒不正确单元的结点号来改变它们。

命令:ESURF,REVE

GUI:main menu>preprocossor>Create>Elements on free surf

步骤5:设置实常数和单元关键字

程序使用九个实常数和好几个单元关键字来控制面─面接触单

元的接触行为。

实常数

9个实常数中,两个(R1和R2)用采定义目标面单元的几何形状,乘下的7个用来控制接触行为。

R1和R2定义目标单元几何形状

FKN 定义法向接触刚度因子

FTOLN 定义最大的渗透范围

ICONT 定义初始靠近因子

PINB 定义“Pinball"区域

PMIN和PMAX 定义初始渗透的容许范围

TAUMAR 指定最大的接触摩擦

命令:R

GUI:main menu> preprocessor>real constant

对实常数FKN,FTOLN,ICONT,PINB,PMAX,和PMIN,你既可以定义一个正值也可以定义一个负值,程序将正值作为比例因子,将负值作为真实值,程序将下面覆盖原单元的厚度作为ICON,FTOLN,PINB,PMAX和PMIN的参考值,例如对ICON,0.1表明初始间隙因子是0.1*下面覆盖层单元的厚度。然而,-0.1表明真实缝隙是0.1,如果下面覆盖层单元是超单元,则将接触单元的最小长度作为厚度。

单元关键字

每种接触单元都饭知好几个关键字,对大多的接触问题缺省的关键字是合适的,而在某些情况下,可能需要改变缺省值,来控制接触行为。

接触算法(罚函数+拉格郎日或罚函数)(KEYOPT(2))

出现超单元时的应力状态(DEYOPT(3))

接触方位点的位置(KEYOPI(4))

刚度矩阵的选择(KEYOPT(6))

时间步长控制(KEYOPT(7))

初始渗透影响(KEYOPT(9))

接触表面情况(KEYOPT(12))

命令:KEYOPT

ET

GUI:main menu>preprocessor>Elemant Type>Add/Edit/Delete

选择接触算法:

对面─面的接触单元,程序可以使用扩增的拉格朗日算法或罚函数方法,通过使用单元关键字KETOPT(2)来指定。

扩张的拉格朗日算法是为了找到精确的拉格朗日乘子而对罚函

数修正项进行反复迭代,与罚函数的方法相比,拉格朗日方法不易引起病态条件,对接触刚度的灵敏度较小,然而,在有些分析中,扩增的拉格朗日方法可能需要更多的迭代,特别是在变形后网格变得太扭曲时。

使用拉格朗日算法的同时应使用实常数FTOLN

FTOLN为搠格朗日算法指定容许的最大渗艉,如果程序发现渗透大于此值时,即使不平衡力和位移增量已经满足了收敛准则,总的求解仍被当作不收敛处理,FTLON的缺省值为0.1,你可以改变这个值,但要注意如果此值太小可能会造成太多的迭代次数或者不收敛。

决定接触刚度

所有的接触问题都需要定义接触刚度,两个表面之间渗量的大小取决了接触刚度,过大的接触刚度可能会引起总刚矩阵的病态,而造成收敛困难,一般来谘,应该选取足够大的接触刚度以保证接触渗透小到可以接受,但同时又应该让接触刚度足够小以使不会引起总刚矩阵的病态问题而保证收敛性。

程序会根据变形体单元的材料特性来估计一个缺省的接触刚度值,你能够用实常数FKN来为接触刚度指定一个比例因子或指定一个真正的值,比例因子一般在0.01和10之间,当避免过多的迭代次数时,应该尽量使渗透到达极小值。

为了取得一个较好的接触刚度值,又可需要一些经验,你可以按下面的步骤过行。

1、开始时取一个较低的值,低估些值要比高估些值好因为由一个

较低的接触刚度导致的渗透问题要比过高的接触刚度导致

的收敛性困难,要容易解决。

2、对前几个子步进行计算

3、检查渗透量和每一子步中的平衡迭代次数,如果总体收敛困难

是由过大的渗透引起的(而不是由不平衡力和位移增量引起

的),那么可能低估了FKN的值或者是将FTOLN的值取得大

小,如果总体的收敛困难是由于不平衡力和位移增量达到收敛

值需要过多的迭代次数,而不是由于过大的渗透量,那么FKN

的值可能被高估。

4、按需要调查FKN或FTOLN的值,重新分析。

选择摩擦类型。

在基本的库仑摩擦模型中,两个接触面在开始相互滑动之前,在它们的界面上会有达到某一大小的剪应力产生,这种状态则作粘合状态(stick)库仑摩擦模型定义了一个等效剪应力。)一旦剪应力超过此值后,两个表面之间将开始相互滑动,这种状态,叫作滑动状态(Sliding)粘合\滑动计算决定什么时候一个点从粘合状态到滑动状

态或从滑动状态变到粘合状态,摩擦系数可以是任一非负值。程序缺省值为表面之间无摩擦,对rough 或bonded 接触(KEYOPT (2)=1(或3),程序将不管给定的MV 值而认为摩擦阻力无限大。

程序提供了一个不管接触压力的故而人为指定最大等效剪应力的选项,如果等效剪应力达到此值时,滑动发生。看图,4─1,为了指定接触界面上最大许可剪应力,设置常数 TAUMAX (缺省为1.0E20),这种限制剪应力的情况一般用于接触压力非常大的时候,以至于用库仑理论计算出的界面剪应力超过了材料的屈服极限。一对TAUMAX 的一个合理高估为y

σ3(y

σ3是材料的mises 屈服应力)。

图4—1 摩擦模式

对无摩擦rough 和bonded 接触,接触单元刚度矩阵是对称的,而涉及到摩擦的接触问题产生一个不对称的刚度,而在每次迭代使用不对称的求解器比对称的求解器需要更多的计算时间,因此ANSYS 程序采用对称化算法。通过采用这种算法大多的摩擦接触问题能够使用对称系统的求解器来求解。如果摩擦应力在整个位移范围内有相当大的影响,并且摩擦应力的大小高度依赖于求解过程。对刚度阵的任何对称近似都可能导致收敛性的降低,在这种情况下,选择不对称求解选项(KEYOPT (6)=1)来改善收敛性。

选择检查接触与否的位置

接触检查点位于接触单元的积分点上,在积分点上,接触单元不渗透进入目标面,然而,目标面能渗透进入接触面,看图4—2。

图4—2 接触检查点位于高斯积分点上

ANSYS面─接触单元使用GAUSS积分点作为缺省值,GAUSS 积分点通常会比Newton-Cotes/robatto结点积分项产生更精确的结果,Newton-cotes/lobatto使用结点本身作为积分点,通过KEYOPT(4)来选择,你想使用的方法,然而,使用结点本身作为积分点仅应该用于角接触问题(看图4─3)。

图4—3 接触检查点位于高斯结点上

然而,使用结点作为接触发现点,可能会导致其它的收敛性问题,例如“滑脱”(结点滑下目标面的边界)看图4─4,对大多的点─面的接触问题,我们推荐使用其它的点─面的接触单元,例如CONTA26、CONTA48和CONTA49。

图4—4 结点滑脱

调整初始接触条件

在动态分析中,刚体运动一般不会引起问题,然而在静力分析中,当物体没有足够的约束时会产生刚体运动,有可能引起错误而终止计算。

在仅仅通过接触的出现来约束刚体运动时,必须保证在初始几何体中,接触对是接触的,换句话说,你要建立模型以便接触对是“刚好接触”的,然而这样作可能会遇到以下问题:

·刚体外形常常是复杂的,很难决定第一个接触点发生在哪儿

·既使实体模型是在初始接触状态,在网格划分后余于数值舍入误差;两个面的单元网格之间也可能会产生小的缝隙。

·接触单元的积分点和目标单元之间可能有小的缝隙。

同理,在目标面和接触面之间可能发生过大的初始渗透,在这种情况下,接触单元可能会高估接触力,导致不收敛或民接触面之间脱离开接触关系。定义初始接触也许是建立接触分析模型时最重要的方面,因此,程序提供了几种方法来调整接触对的初始接触条件。

注意:下面的技巧可以在开始分析时独立执行成几个联合起来执行,它们是为了消除由于生成网格造成的数值舍入误差而引起的小缝隙或渗透,而不是为了改正网格或几何数据的错误。

1、使用实常数ICONT来指定一个好的初始接触环,初始接触环

是指沿着目标面的“调整环”的深度,如果没有人为指定

ICONT的值,程序会根据几何尺寸来给ICONT提供一个小

值,同时输出一个表时什么值被指定的警告信息,对ICONT

一个正值表示相对于下面变形体单元厚度的比例因子,一个

负值表示接触环的真正值,任何落在“调整环”敬域内的接

触检查点被自动移到目标面上,(看图4—5(a))建议使用一个

小的ICONT值否则,可能会发生大的不连续(看图4─5(b))

图4—5 用ICON进行接触面的调整

(a) 调整前

(b) 调整后

2、使用实常数PMIN和PMAX来指定初始容许的渗透范围,当

指定PMAX或PMIN后,在开始分析时,程序会将目标面移

到初始接触状态,如果初始渗透大于PMAX,程序会调整目标

面的减少渗透,接触状态的初始调节仅仅通过平移来实现。

对给定载或给定位移的刚性目标面将会执行初始接触状态的初始调节。同样,对没有指定边界条件的目标面也可以进行初始接触的调整。

当目标面上的节点,有给的零位移值时,使用PMAX和PMIN 的初始调节将不会被执行。

注意:ANSYS程序独立地处理目标面上节点的自由度,例如:如果你指定自中度UX值为“0”,那么,沿着X方向就没有初始调查,然而,在Y和Z方向仍然会激活PMAX和PMIN选项。

初始状态调整是一个迭代过程,程序最多进行20次迭代,如果目标面不能进入可接受的渗透范围,程序会给出一个警告信息,你可能需要调整你的初始几何模型。

图4─6给出了一个初始接触调整迭代失败的例子。目标面的UY 被约束住。因此,初始接触唯一容许的调整是在X方向,然而,在这个问题中,刚性目标面在X方向的任何运动都不会引起初始接触。

图4—6 一个初始调整失败的例子

3、设置KEYOPI(9)=1来消除初始渗透,看图4─7。

图4—7 消除初始始渗透

在某些情况下,例如过盈装配问题,期望有过度的渗,为了缓解收敛性困难,在第一个载荷步中设置KEYOTI(9)=2来使过度渗透渐进到0,看图4─8。当使用这种方法时,在第一个载步中不要给定其它任何载荷,也就是说要保证载荷是渐进的(KBC,0)

图4—8 渐进初始渗透

在开始分析时,程序会给出每个目标面的初始接触状态的输出信息,(在输出窗口或输出文件中),这个信息有助于决定每个目标面的最大渗透成最小间隙。

对于给定的目标面如果没有发现接触,可能是目标面离接触面太远(超出了Piaball区域或者是接触/目标单元已经被杀死。

决定接触状态和Pinball区域。

接触单元相对于目标面的运动和位置决定了接触单元的状态;程序检测每个接触单元并给出一种状态

·STAT=0 未合的远区接触

·STAT=1 未合的近区接触

·STAT=2 滑动接触

·STAT=3 粘合接触

当目标面进入pinball区域后,接触单元就被当作未合上的近区域接触,pinball区域是以接触单元的积分点为中心的。使用实常数PINB来为pinball指定一个比例因子(正值)或其实值(负值),缺省时,程序将pinball区域定义为一个以4*变形体单元厚度为半径的圆(对2-D问题)或球(对3-D问题)。

检查接触的计算时间依赖于pinball区域的大小,远区接触单元的计算是简单的且计算时间较少,近区接触计算将要接触的接触单元是较慢的,并是较复杂,当单元已经接触时,计算最为复杂。

如果刚性面有好几个凸形区域,为了克服伪接触定义,设置一个合适的pinball区域是有用的,名而对大多数问题,缺省值是合适的选择表面作用模式

通过设置kcyopt(12)来选择下面的某种作用模式

·法问单边接触(KEYOPT(12)=0)

·.粗糙接触,用来模拟无滑动的,表面相当粗糙的摩擦接触

问题,这种设置对应于摩擦系数无限大(MU),因此用户

定义的摩擦系的(MU)被忽略KEYOPT(12)=1)

·.不分开的接触,用来模拟那种一是接触就再不分开的问题,

这种不分开是指对法方接触而言,允许有相对滑动。

(KEYOPT(12)=2)

·.绑定接触用来模拟那种接触一是发生表面在所省方向都

被绑定的问题。一旦接触就再也不能脱开也不允许有相对

滑动(KEYOPT(12)=3)

用超单元建立接触模型

面一面的接触单元能模拟刚体和另一个有的运动的线单性体的接触,而线单性体又以体用超单元来建模,这大大降低了进行接触代的自由度数,记住任荷接触结定都必须是超单元的主自由度。

既然超单元仅仅由一组保留的结点自由度组成,它没有用来定义接触的表面几何形状,因此,必须在形成超单元之前在单元表面上成接触单元,来自超单元的信息包括结点连结和组合刚度,但是没有材料特性和应力状态,(是否轴支称,平面应力或平面应变),一个限制是接触单元的材料特性设置必须与形成超单元之前的原始单元的材料特性相同。

使用KEYOPT(3)来提供接触分析的信息,对2─D单元(CONTA171 CONTA172)关键字选项如下所示:

·.不使用超单元(KETOPT(3)=0)

·.轴对称(KEYOPT(3)=1)

·.平面应变或单位厚度的平面应力(KEYOPT(3)─2)

·.需要厚度输X的平面应力(KEYOPT(3)=3),对这种情

况使用实常数的R2来指定指定厚度

对3─D单元(CONTA173,CONTA174)关键字选项如下示:·使用H单元(KETOPI(3)=0)

·使用超单元(KEYOPI(3)=1)

考虑厚度影响

程序够用KEYOPI(11)来考虑壳(2-D和3-D)和梁(2-D)的厚度缺省时,程序不考虑单元厚度,用或中面来表示它。当设置KFTOPI(11)=1时则考虑梁或壳的厚度,从底面或顶面来计算接触距离,建模时要考虑到厚度,记住刚性目标面会向任一边移动,半个梁或壳单元的厚度,当使用壳单元181号时,在变形期间厚度的变化也将被考虑。

使用时间步长控制

时间步长控制是一个自动时间步长特征,这个特征预测什么时间接触单元的状态,将发生变化或者需要二分当前的时间步长,使用KEYOPT(7)来选择下列四种行为之一来控制时间步长。KEYOPT (7)=0时不提供控制,KEYOPT(7)=3提供最多的控制。

·KEYOPI(7)=0,设有控制,时间步的大小不受预测影响,当自动时间步长被激活且允许一个很小的时间步长时,这个设

置是合适的。

·KETOPI(7)=1如果一次迭代期间有太大的渗透发生或者接触状态突然变化,则进行时间步长二分。

·KEYOPI(7)=2对下一个子步预测一个合理的时间增量

·KETOPI(7)=3对下一个子步,预测一个最小的时间增量

使用死活单元选项

面─面的接触单元允许激活或杀死单元,能够在分析的某一阶段中杀死这个单元而在以后的阶段再重新激活它,这个特征对于模拟复杂的金属戍形过程是有用的、在此过程的不同分析阶段有多个目标需要和接触面相互作用,回弹模拟常常需要在成形过程的后期移走刚性工具。

步骤六:

控制刚性目标的运动。

按照物体的原始外形来建立的且整个表刚性目标面是面的运动是通过“pilot”结点上的给定来定义的,(如果没有定义“pilot”结点,则通过刚性目标面上的不同结点。)

为了控制整个目标面的运动,在下面的任何情况下都必须使用"pilot"结点。

·目标面上作用着给定的外力

·目标面发生旋转

·目标面和其它单元相连(例结构质量单元)

"pilot"结点的厚度代表着整个刚性面的运动,你可以在"pilot"结点上给定边界条件(位移、初速度)集中载转动等等,为了考虑刚体的质量,在"pilot"结点上定义一个质量单元。

当使用"pilot"结点时,记住下面的几点局限性

·每个目标面只能有一个“Pilot"的结点

·.圆、圆锥、圆柱、球的第一个结点(结点工)是”pilot“结

点,你不能另外定义或改变"pilot"结点

·.程序忽略不是"lilot"结点的所有其它结点上的边条件。

·.只有“pilot”结点能与其它单元相连

·.当定义了“pilot”结点后,不能使用约束方程(CF)或结

点来耦合(CP)来控制目标面的自由度,如果你在刚性面

上给定任意载荷或者约束,你必须定义“pilot”结点,是在

"pilot"结点上加载,如果没有使用“pilot”结点,则只能有

刚体运动。

在每个载步的开始,程序检查每个目标面的边界条件,如果下面的条件都满足,那么程序将目标面作为固定处理:

·在目标面结点上没有明确定义边界条件或给定力

·.目标面结点没有和其它单元相连

·.没有目标面结上使用约束方程或结点来合

在每个载体步的末尾,程序将会放松被内部设置的约束条件

步骤7:给变形体单元加必要的边界条件

现在可以按需要加上任你边界条件。加载过程与其它的分析类型相同

步骤8:定义求解和载步选项

接触问题的收敛性随问题不同而不同,下面列式了一些典型的在大多数面—面的接触分析中推荐使用的选项

·时间步长必须足够以描述适当的接触。如果时间步太大,则接触力的光滑传递会被破坏,设置精确时间步长的可信赖的方法是打开自动时间步长。

命令:Autots,on

GUI:Main Menu>Solution>-load step opts-Time/Frequence>Time&Time step

/Time& substeps

·.如果在迭代期间接触状态变化,可能发生不连续,为了避免收敛太慢,使用修改的刚度阵,将牛顿一拉普森选项设置成FULL 命令:NROPT,FULL,,OFF

GUI:Main Menu>Solution>Analysis options

不要使用自下降因子,对面一面的问题,自适应下降因子通常不会提供任何帮助,因此我们建议关掉它。

·设置合理的平衡迭代次数,一个合理的平衡迭代次数通常在25和50之间

有限元ansys静力分析的一个小例子

有限元 学院:机电学院 专业: 姓名: 学号:

一、问题描述 如图所示的平面,板厚为0.01m,左端固定,右端作用50kg的均布载荷,对其进行静力分析。弹性模量为210GPa,泊松比为0.25. 二、分析步骤 1.启动ansys,进入ansys界面。 2.定义工作文件名 进入ANSYS/Multiphsics的的程序界面后,单击Utility Menu菜单下File中Change Jobname的按钮,会弹出Change Jobname对话框,输入gangban为工作文件名,点击ok。 3.定义分析标题 选择菜单File-Change Title在弹出的对话框中,输入Plane Model作为分析标题,单击ok。 4.重新显示 选择菜单Plot-Replot单击该按钮后,所命令的分析标题工作文件名出现在ANSYS 中。 5.选择分析类型 在弹出的对话框中,选择分析类型,由于此例属于结构分析,选择菜单Main Menu:Preferences,故选择Structural这一项,单击ok。 6.定义单元类型 选择菜单Main Menu-Preprocessor-Element Type-Add/Edit/Delete单击弹出对话框中的Add按钮,弹出单元库对话框,在材料的单元库中选Plane82单元。即在左侧的窗口中选取Solid单元,在右侧选择8节点的82单元。然后单击ok。 7.选择分析类型 定义完单元类型后,Element Type对话框中的Option按钮被激活,单击后弹出一个对话框,在Elenment behavior中选择Plane strs w/ thk,在Extra Element output 中,选择Nodal stress,单击close,关闭单元类型对话框。 8.定义实常数 选择菜单Main Menu-Preprocessor-Real Constants Add/Edit/Delete执行该命令后,在弹出Real Constants对话框中单击Add按钮,确认单元无误后,单击ok,弹出Real Constants Set Number 1,for Plane 82对话框,在thickness后面输入板的厚度0.01单击ok,单击close。 9.定义力学参数 选择菜单Main Menu-Preprocessor-Material Props-Material Model 在弹出的对

ANSYS动力学分析

第5章动力学分析 结构动力学研究的是结构在随时间变化载荷下的响应问题,它与静力分析的主要区别是动力分析需要考虑惯性力以及运动阻力的影响。动力分析主要包括以下5个部分:模态分析:用于计算结构的固有频率和模态。 谐波分析(谐响应分析):用于确定结构在随时间正弦变化的载荷作用下的响应。 瞬态动力分析:用于计算结构在随时间任意变化的载荷作用下的响应,并且可涉及上述提到的静力分析中所有的非线性性质。 谱分析:是模态分析的应用拓广,用于计算由于响应谱或PSD输入(随机振动)引起的应力和应变。 显式动力分析:ANSYS/LS-DYNA可用于计算高度非线性动力学和复杂的接触问题。 本章重点介绍前三种。 【本章重点】 ?区分各种动力学问题; ?各种动力学问题ANSYS分析步骤与特点。 5.1 动力学分析的过程与步骤 模态分析与谐波分析两者密切相关,求解简谐力作用下的响应时要用到结构的模态和振型。瞬态动力分析可以通过施加载荷步模拟各种何载,进而求解结构响应。三者具体分析过程与步骤有明显区别。 5.1.1 模态分析 1.模态分析应用 用模态分析可以确定一个结构的固有频率利振型,固有频率和振型是承受动态载荷结构设计中的重要参数。如果要进行模态叠加法谐响应分析或瞬态动力学分析,固有频率和振型也是必要的。可以对有预应力的结构进行模态分析,例如旋转的涡轮叶片。另一个有用的分析功能是循环对称结构模态分析,该功能允许通过仅对循环对称结构的一部分进行建模,而分析产生整个结构的振型。 ANSYS产品家族的模态分析是线性分析,任何非线性特性,如塑性和接触(间隙)单元,即使定义也将被忽略。可选的模态提取方法有6种,即Block Lanczos(默认)、Subspace、Power Dynamics、Reduced、Unsymmetric、Damped及QR Damped,后两种方法允许结构中包含阻尼。 2.模态分析的步骤

ansys动力学分析全套讲解

第一章模态分析 §模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型),即结构的固有频率和振型,它们是承受动态载荷结构设计中的重要参数。同时,也可以作为其它动力学分析问题的起点,例如瞬态动力学分析、谐响应分析和谱分析,其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析,后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性,如塑性和接触(间隙)单元,即使定义了也将被忽略。ANSYS提供了七种模态提取方法,它们分别是子空间法、分块Lanczos法、PowerDynamics法、缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存在阻尼。后面将详细介绍模态提取方法。 §模态分析中用到的命令 模态分析使用所有其它分析类型相同的命令来建模和进行分析。同样,无论进行何种类型的分析,均可从用户图形界面(GUI)上选择等效于命令的菜单选项来建模和求解问题。 后面的“模态分析实例(命令流或批处理方式)”将给出进行该实例模态分析时要输入的命令(手工或以批处理方式运行ANSYS时)。而“模态分析实例(GUI方式)” 则给出了以从ANSYS GUI中选择菜单选项方式进行同一实例分析的步骤。(要想了解如何使用命令和GUI选项建模,请参阅<>)。<>中有更详细的按字母顺序列出的ANSYS命令说明。 §模态提取方法 典型的无阻尼模态分析求解的基本方程是经典的特征值问题: 其中: =刚度矩阵, =第阶模态的振型向量(特征向量), =第阶模态的固有频率(是特征值), =质量矩阵。 有许多数值方法可用于求解上面的方程。ANSYS提供了7种方法模态提取方法,下面分别进行讨论。 1.分块Lanczos法 2.子空间(Subspace)法 Dynamics法

ANSYS静力分析的简单步骤

第一步,启动工作台软件,然后选择与启动DS模块弹出得界面。 第二步,导入三维模型。根据操作步骤进行。首先,单击“几何体”,选择“文件”,然后选择弹出窗口中的3D模型文件,如果当时catia文件格式不符,可以把三维图先转换为“.stp”的格式,即可导入。 第三步,选择零件材料:文件导入软件后,在这个时候,依次选择“几何”下的“零件”,并且在左下角的“Details of ‘Part’”中以调整零件材料属性,本次钟形壳的材料是刚。 第四步,划分网格:选择“Project”树中的“Mesh”,右键选择“Generate Mesh”即可在这一点上,你可以在左下角的“网格”对话框的细节调整网格的大小(体积元)。 第五步,添加类型分析:第一选择顶部工具栏上的“分析”按钮,添加需要的类型分析,因为我们需要做的是在这种情况下的静态分析。所以选择结构静力。 第六步,添加固定约束:首先选择“Project”树中的“Static Structural”按钮,右键点击支持插入固定树。这时候在左下角的“Details of ‘Fixed Support’”对话框中“Geometry”会被选中,会要求输入固定的支撑面。在这种情况下,固定支架的类型是表面支持,确定六凹面(此时也可点击“Edge”来确定“边”)。然后一直的按住“CTRL”键,连续选择其它几个弧面为支撑面,在点击“Apply”进行确认, 第七步,添加载荷:选择“Project”树中的“结构静力”,右键选择“Insert”中的“Force”,然后在选择载荷的作用面,再次点击“Apply”按钮进行确定。 第八步,添加变形:右键点击选择“Project”树中的“Solution”,随后依次选择插入,变形,Total”,添加变形。 第九步,添加等效应变:右键单击“项目”的树,“>插入应变->解决方案->添加等效,等效应变。 第十步,添加等效应力:首先右键点击“Project”树中的“Solution—>Insert —> Stress—>Equivalent”,添加等效应力。 第十一步,求解:右键点击选择“Project”树中的“Solution”,随后选择“Solve”求解

ANSYS静力学分析APDL建模实例-应力集中

计算分析模型如图所示, 习题文件名: scf 材料参数:E=205GPa, v = 0.3 力载:4500N 注意单位的一致性:使用N, mm, MPa单位制 建模教程 在ANSYS工作文件夹内新建“stress concentration factor”目录,以存放模型文件。 注意定期保存文件,注意不可误操作,一旦误操作,不可撤销。 1.1 进入ANSYS 开始→程序→ANSYS 14.5→Mechanical APDL Product Launcher14.5→然后在弹出的启动界面输入相应的working directory及文件名scf 如通过Mechanical APDL 14.5进入,则进入预设的working directory working directory必须设置在电脑最后一个分区(因为教学用电脑只有最后一个分区不受系统保护) 至此ANSYS静力学分析模块启动,ANSYS在“stress concentration factor”目录下自动创建了.log、.err等必要的文件。 2.2设置计算类型 ANSYS Main Menu: Preferences →select Struc tural → OK 2.3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4 node 182 →OK (back to Element Types window)→ Options… →select K3: Plane Strs w/thk →OK→Close (the Element Type window)

ANSYS-结构稳态(静力)分析之经典实例-命令流格式

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

实验四-五:结构静力分析与ANSYS模态分析

注:3月20号,周二课程内容主要是完成下面实验四 特别注意:本周六没课,本五周23号,8:00--12:00有课------------------------------------------------------------------------------------- 实验四MEMS薄膜压力传感器静力学分析 一、实验目的 1、掌握静力学分析 2、验证理论分析结果 3、对不同形状膜的分析结果进行对比 二、实验器材 能够安装ANSYS软件,内存在512MHz以上,硬盘有5G空间的计算机 三、实验说明 (一)基本思路 1、建模与网格化 2、静力学分析 3、对结果进行分析和比较 (二)问题描述: 由于许多压力传感器的工作原理是将受压力作用而变形的薄膜硅片中的应变转换成所需形式的电输出信号,所以我们要研究比较一下用什么样形状的膜来作为压力传感器的受力面比较好。我们比较的膜形状有三种,分别是圆形. 正方形. 长方形。在比较的过程中,三种形状膜的面积.,厚度和承受的压力是都是相等的。设置参数具体为:F=0.1MPa, EX=1.9e11,PRXY=0.3,DENS=2.33e3.单元尺寸为5e-006。为了选择合适的网格化类型,首先我们拿圆的结构进行一下比较,最后选择比较接近理论计算的网格化类型,通

过比较,我们知道映射网格化类型比较优越,所以后面的两种类型膜结构选择了映射网格化。 四、实验内容和步骤 圆形薄膜1 1.先建立一个圆形薄膜:Main Menu>Preprocessor>modeling>Create>volumes>solid cylinder.弹出以个对话 框如图,输入数据如图4-1,单击OK. 图4-1 2.设置单元类型:Main Menu>Preprocessor>element type>add/edit/delete,弹出一个对话框,点击add,显示library of element type对话框如图:在library of element type下拉列表框中选择structural solide 项,在其右侧下拉表框中选择brick 8node 45选项,单击OK. 在点击close.如图4-2. 图4-2 3.设置材料属性:Main Menu>Preprocessor>material props> material models,弹出一个对话框,在material

ANSYS WORKBENCH 11.0静力结构分析

ANSYS WORKBENCH 11.0培训教程(DS)

第四章 静力结构分析

序言 ?在DS中关于线性静力结构分析的内容包括以下几个方面: –几何模型和单元 –接触以及装配类型 –环境(包括载荷及其支撑) –求解类型 –结果和后处理 ?本章当中所讲到的功能同样适用与ANSYS DesignSpace Entra及其以上版本. –本章当中的一些选项可能需要高级的licenses,但是这些都没有提到。 –模态,瞬态和非线性静力结构分析在这里没有讨论,但是在相关的章节当中将会有所阐述。

线性静力分析基础 ?在线性静力结构分析当中,位移矢量{x} 通过下面的矩阵方程得到: 在分析当中涉及到以下假设条件: –[K] 必须是连续的 ?假设为线弹性材料?小变形理论 ?可以包括部分非线性边界条件–{F} 为静力载荷 ?不考虑随时间变化的载荷 ?不考虑惯性(如质量,阻尼等等)影响 ?在线性静力分析中,记住这些假设是很重要的。非线性分析和动力学分析将在随后的章节中给予讨论。 []{}{} F x K =

A. 几何结构 ?在结构分析当中,可以使用所有DS 支持的几何结构类型. ?对于壳体,在几何菜单下厚度选项是必须要指定的。 ?梁的截面形状和方向在DM已经指定并且可以自动的传到DS模型当中。 –对于线性体,仅仅可以得到位移结果. ANSYS License Availability DesignSpace Entra x DesignSpace x Professional x Structural x Mechanical/Multiphysics x

Ansys静力分析

Ansys静力分析实例: 1 问题描述: 如图所示支架简图,支架材料为结构钢,厚度10mm,支架左侧的两个通孔为固定孔,顶面的开槽处受均布载荷,载荷大小为500N/mm。 2 启动Ansys Workbench,在界面中选择Simulation启动DS模块。

3 导入三维模型,操作步骤按下图进行,单击“Geometry”,选择“From File”。 从弹出窗口中选择三维模型文件,如果文件格式不符,可以把三维图转换为“.stp”格式文件,即可导入,如下图所示。 4 选择零件材料:文件导入后界面如下图所示,这时,选择“Geometry”下的“Part”,在左下角的“Details of ‘Part’”中可以调整零件材料属性。

5 划分网格:如下图,选择“Project”树中的“Mesh”,右键选择“Generate Mesh”即可。【此时也可以在左下角的“Details of ‘Mesh’”对话框中调整划分网格的大小(“Element size”项)】。

生成网格后的图形如下图所示:

6 添加分析类型:选择上方工具条中的“New Analysis”,添加所需做的分析类型,此例中要做的是静力分析,因此选择“Static Structural”,如下图所示。 7 添加固定约束:如下图所示,选择“Project”树中的“Static Structural”,右键选择“Insert”中的“Fixed Support”。

这时左下角的“Details of ‘Fixed Support’”对话框中“Geometry”被选中,提示输入固定支撑面。本例中固定支撑类型是面支撑,因此要确定图示6位置为“Face”,【此处也可选择“Edge”来选择“边”】 然后按住“CTRL”键,连续选择两个孔面为支撑面,按“Apply”确认,如下图所示。

ANSYS动力学分析的几个入门例子

ANSYS动力学分析的几个入门例子 问题一:悬臂梁受重力作用发生大变形,求其固有频率。图片附件: 1.jpg ( 4.85 K ) 基本过程: 1、建模 2、静力分析 NLGEOM,ON STRES,ON 3、求静力解 4、开始新的求解:modal STRES,ON UPCOORD,1,ON 修正坐标 SOLVE... 5、扩展模态解 6、察看结果

/PREP7 ET,1,BEAM189 !使用beam189梁单元MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,210e9 MPDATA,PRXY,1,,0.3 MPDATA,DENS,1,,7850 SECTYPE, 1, BEAM, RECT, secA, 0 !定义梁截面secA SECOFFSET, CENT SECDATA,0.005,0.01,0,0,0,0,0,0,0,0 K, ,,,, !建模与分网 K, ,2,,, K, ,2,1,, LSTR, 1, 2 LATT,1, ,1, , 3, ,1 LESIZE,1, , ,20, , , , ,1 LMESH, 1 FINISH /SOL !静力大变形求解 ANTYPE,0 NLGEOM,1 PSTRES,ON !计及预应力效果 DK,1, , , ,0,ALL, , , , , , ACEL,0,9.8,0, !只考虑重力作用 TIME,1 AUTOTS,1 NSUBST,20, , ,1 KBC,0 SOLVE FINISH /SOLUTION ANTYPE,2 !进行模态求解 MSA VE,0 MODOPT,LANB,10 MXPAND,10, , ,0 !取前十阶模态 PSTRES,1 !打开预应力效应MODOPT,LANB,10,0,0, ,OFF UPCOORD,1,ON !修正坐标以得到正确的应力PSOLVE,TRIANG !三角化矩阵 PSOLVE,EIGLANB !提取特征值和特征向量FINISH /SOLU

ansys静力分析实例1

结构分析实验指导书 1.问题描述: 这是一个关于角支架的单载荷步的结构静力分析。如图所示,左上角的销孔由于焊接而被固定死。右下角的销孔上作用一分布力。本问题的目标是熟悉ANSYS分析的基本过程。使用的是美国的单位体系。 材料的氏模量为30E6 psi,泊松比0.27。 2.几何建模: 第一步:定义矩形 1.Main Menu> Preprocessor> Modeling> Create> Areas> Rectangle> By Dimensions 2.Enter the following: X1 = 0 ,X2 = 6,Y1 = -1,Y2 = 1 3.Apply to create the first rectangle. 4.Enter the following: X1 = 4,X2 = 6,Y1 = -1,Y2 = -3 5.OK to create the second rectangle and close the dialog box.

第二步:更改绘图属性和重绘。 1.Utility Menu> Plot Ctrls> Numbering 2.Turn on area numbers. 3.OK to change controls, close the dialog box, and replot. 4.Toolbar: SAVE_DB. 第三步:更改工作平面为极坐标系并创建第一个圆 1.Utility Menu> WorkPlane> Display Working Plane (toggle on) 2.Utility Menu> WorkPlane> WP Settings 3.Click on Polar. 4.Click on Grid and Triad. 5.Enter 0.1 for snap increment.

ansys动力学分析全套讲解

a n s y s动力学分析全套讲 解 This model paper was revised by the Standardization Office on December 10, 2020

第一章模态分析 §模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型),即结构的固有频率和振型,它们是承受动态载荷结构设计中的重要参数。同时,也可以作为其它动力学分析问题的起点,例如瞬态动力学分析、谐响应分析和谱分析,其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析,后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性,如塑性和接触(间隙)单元,即使定义了也将被忽略。ANSYS提供了七种模态提取方法,它们分别是子空间法、分块Lanczos法、PowerDynamics法、缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存在阻尼。后面将详细介绍模态提取方法。 §模态分析中用到的命令 模态分析使用所有其它分析类型相同的命令来建模和进行分析。同样,无论进行何种类型的分析,均可从用户图形界面(GUI)上选择等效于命令的菜单选项来建模和求解问题。 后面的“模态分析实例(命令流或批处理方式)”将给出进行该实例模态分析时要输入的命令(手工或以批处理方式运行ANSYS时)。而“模态分析实例(GUI方式)” 则给出了以从ANSYS GUI中选择菜单选项方式进行同一实例分析的步骤。(要想了解如何使用命令和GUI选项建模,请参阅<>)。<>中有更详细的按字母顺序列出的ANSYS命令说明。 §模态提取方法 典型的无阻尼模态分析求解的基本方程是经典的特征值问题: 其中:

ansys静力分析基本步骤

定义单元类型: Main menu | preprocessor | Element Type | Add/Edit/Delete命令出现Element Types 对话框,单击Add按钮,出现Library of Element Types 对话框。在列表框中,选择单元类型,然后单击OK。 定义材料性能参数: 选择Main menu | preprocessor | Material props | Material Models 命令,出现Define Material Model Behavior 对话框。在Material Models Available 一栏中依次双击Structural、Linear、Elastic、isotropic 选项。在EX输入栏中输入(弹性模量)2.2E11,(2.2E11=220Gpa) 在PRXY中输入(泊松比)0.3单击OK关闭对话框。 创建几何模型,划分网格: 选择Utility Menu | file | import 导入有限元模型。 通过Main menu | preprocessor | Modeling | create | point/line.... 给有限元模型添加相关的重要的点和线。 选择Main menu | preprocessor | Meshing | Mesh tool | 勾选Smart Size ,然后在下拉菜单Mesh 中选择V olumes,然后单击凸起按钮Mesh | 弹出选择对话框,用鼠标选取有限元模型,使得整个模型变色后单击OK开始网格划分,划分结束后关掉弹出的多余对话框。 运行计算: Main menu | Solution | analysis Type | New Analysis | 选择static 单击OK Main menu | Solution | Define Loads | Apply |Structural | Displacement | On Nodes/point/line/area.... 出现拾取菜单,用鼠标选择需要添加约束的点,线或者面,单击OK, 出现另一对话框,在列表对话框选择需要约束的方向。UX/UY/UZ或者全部。在Displacement Value中输入O 这一步是定义位移约束。如有多个约束按上步骤再做几次。 Main menu | Solution | Define Loads | Apply | Structural | pressure | On Nodes/point/line/area.... 。出现拾取菜单,用鼠标选择选择要加载力的点、线、或者面。单击OK出现加载力对话框

ANSYS静力学分析

实验报告(四) 一、实验目的:掌握静力学分析;验证理论分析结果;对不同形状膜的分析结果进行对比。 二、实验器材:能够安装ANSYS^件,内存在512MHz以上,硬盘有5G空间的计算机。 三、实验说明 1 、基本思路:建模与网格化;静力学分析;对结果进行分析和比较。 2、问题描述:许多压力传感器的工作原理是将受压力作用而变形的薄膜硅片中的应变转换成所需形式的电输出信号,比较的膜形状有圆形. 正方形. 长方形。 四、实验内容及步骤(一)、圆形薄膜(自由式网格划分) 1 、先建立一个圆形薄膜:Main Menu —>Preprocessor —> M odeling —>Create —>Volumes—>Solid Cylinder. 弹出 一个对话框,输入WPX=0 ,WPY=0 ,Radius=300e-6,Depth=13.887e-6 ,单击[OK] 。设置单元类型:Main Menu —>Preprocessor —>Element Type —>Add/Edit/Delete, 弹出一个对话框,点击[Add], 在Library Of Element Type 下拉列表框中选择Structural Solide 项,在其右侧下选择Brick 8node 45 选项,单击[OK] 。再点击

Close.

2、设置材料属性: Main Menu — 'Preprocessor — 'Material Props — > Material Models, 弹出一个对话框,在 Material Models Available 下双击打开 Structural — 'Linear —'Elastic — 'Isotropic ,又弹出一对话框,在 EX 后面输 入1.9E11,在PRXY 后面输入栏中输入 0.3,在双击Density, 在DENS 后面输入2.33e3,单击[0K ]。关闭对话框。设置单元 尺寸:Main Menu — 'Preprocessor — 'Meshing — >MeshTool, 弹出一个对话框,单击 Global 中的Set 按钮,弹出Global Element Sizes 对话框,输入 size=5e-006,然后单击[0K ] < 09060242-44-sunguoliang 3、 采用自由式网格化生成单元:Main Menu'Preprocessor'meshing'meshtool 弹出一个对话框,选 1 VOLUMES TYPE NUM AN MAR 20 2012

ansys动力学分析全套讲解 (2)

第一章模态分析 §1.1模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型),即结构的固有频率和振型,它们是承受动态载荷结构设计中的重要参数。同时,也可以作为其它动力学分析问题的起点,例如瞬态动力学分析、谐响应分析和谱分析,其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析,后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性,如塑性和接触(间隙)单元,即使定义了也将被忽略。ANSYS 提供了七种模态提取方法,它们分别是子空间法、分块Lanczos法、PowerDynamics法、缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存在阻尼。后面将详细介绍模态提取方法。 §1.2模态分析中用到的命令 模态分析使用所有其它分析类型相同的命令来建模和进行分析。同样,无论进行何种类型的分析,均可从用户图形界面(GUI)上选择等效于命令的菜单选项来建模和求解问题。 后面的“模态分析实例(命令流或批处理方式)”将给出进行该实例模态分析时要输入的命令(手工或以批处理方式运行ANSYS时)。而“模态分析实例(GUI方式)”则给出了以从ANSYS GUI中选择菜单选项方式进行同一实例分析的步骤。(要想了解如何使用命令和GUI选项建模,请参阅<>)。<>中有更详细的按字母顺序列出的ANSYS命令说明。 §1.3模态提取方法 典型的无阻尼模态分析求解的基本方程是经典的特征值问题: 其中: =刚度矩阵, =第阶模态的振型向量(特征向量), =第阶模态的固有频率(是特征值), =质量矩阵。 有许多数值方法可用于求解上面的方程。ANSYS提供了7种方法模态提取方法,下面分别进行讨论。 1.分块Lanczos法 2.子空间(Subspace)法 3.Power Dynamics法 4.缩减(Reduced /Householder)法 5.非对称(Unsymmetric)法 6.阻尼(Damp)法(阻尼法求解的是另一个方程,参见<>中关于此法的详细信息) 7.QR阻尼法(QR阻尼法求解的是另一个方程,参见<>中关于此法的详细信息) 注意—阻尼法和非对称法在ANSYS/Professional中不可用。

相关文档