ANSYS的非线性收敛准则

ANSYS的非线性收敛准则
ANSYS的非线性收敛准则

ANSYS的非线性收敛准则

ANSYS, 准则, 非线性ANSYS, 准则, 非线性

ANSYS的非线性收敛准则描述:计算收敛过程图

图片:

CNVTOL, Lab, V ALUE, TOLER, NORM, MINREF

ANSYS中,非线性收敛准则主要有力的收敛,位移的收敛,弯矩的收敛和转角的收敛。一般用力的控制加载时,可以使用残余力的2-范数控制收敛;而位移控制加载时,最好用位移的范数控制收敛。

When SOLCONTROL,ON, TOLER Defaults to 0.005 (0.5%) for force and moment, and 0.05 (5%) for displacement when rotational DOFs are not present.

When SOLCONTROL,OFF, defaults to 0.001 (0.1%) for force and moment.

收敛精度一般可放宽至5%,以提高收敛速度。

加快收敛的方法有一下几种:

1可以增大荷载子步数,nsubst,nsbstp,nsbmn,carry

2修改收敛准则,cnvtol,lab,value,toler,norm,minref

3 打开优化的非线性默认求解设置和某些强化的内部求解算法,

solcontrol,key1,key2,key3,vtol (一般情况下,默认是打开的)

4重新划分网格,网格的单元不宜太大或太小, 一般在5~10厘米左右

5 检查模型的正确性

下面计算收敛过程图中的各个曲线的具体含义是什么?

非线性计算是一个迭代计算的过程,曲线表示两次迭代之间的误差,图中分别表示力和位移在迭代过程中的每次迭代之间的误差

关于ansys中收敛准则(cnvtol)理解

ansys中依据缺省的收敛准则,程序将对不平衡力SRSS与V ALUE*TOLER的值进行比较;而V ALUE的缺省值是在SRSS和MINREF中取较大值。现假如TOLER 的缺省值是0.1的话,这个准则是不是可以理解成后一次的SRSS是前一次的SRSS的01倍就收敛啦?

请指点

我是这样理解的例如下面的命令流:

cnvtol,f,5000,0.0005,0

cnvtol,u,10,0.001,2

如果不平衡力(独立的检查每一个自由度)小于等于5000*0.0005(也就是2.5),并且如果位移的变化小于等于10*0.001时,认为子步是收敛的。

ANSYS中收敛准则,程序默认力与位移共同控制,并且收敛的控制系数好像是0.001。这样的收敛精度一般很难使塑性分析收敛,对于一般的塑性分析收敛问题,前几个荷载步(弹性阶段)用力与位移共同控制,进入塑性后用力控制或位移控制,也可以先用力后用位移控制(位移控制比较容易收敛),至于控制系数取多少,自己根据需要逐步放大直至收敛!也有人建议最后用能量来控制收敛,

convergence value 是收敛值,convergence norm是收敛准则。ansys可以用cnvtol 命令,如:cnvtol,f,10000,0.00001,2,,其中f是指采用力结果,10000是收敛绝对值,0.00001是收敛系数,2是收敛2范数。

收敛准则应该是指选取那种结果进行收敛判定,通常有三种选择,分别是力(f),位移(u)、和能量。当然这三种形式可以单独使用也可以联合使用。收敛准则的

另一层意思应该是选取什么范数形式(1、2、3范数)。一般结构通常都选取2范数格式。

而收敛值只是收敛准则中的一部分,如cnvtol命令中的收敛绝对值与收敛系数的乘积就应该是你所指的收敛值(convergence value)。

ansys 使用收敛准则有L1,L2,L~~(无穷大)三个收敛准则。

在工程中,一般使用收敛容差(0.05)就可以拉。

建议使用位移收敛准则( cnvtol,u,0.05,,, )与力收敛准则( cnvtol,f,0.05,,, )。因为仅仅只使用一个收敛准则,会存在较大的误差。

假如你只能是使用一个收敛准则,建议你提高收敛容差(0.01以下)。

ansys计算非线性时会绘出收敛图,其中横坐标是cumulative iteration number 纵坐标是absolute convergence norm。他们分别是累积迭代次数和绝对收敛范数,用来判断非线性分析是否收敛。

ansys在每荷载步的迭代中计算非线性的收敛判别准则和计算残差。其中计算残差是所有单元内力的范数,只有当残差小于准则时,非线性叠代才算收敛。ansys的位移收敛是基于力的收敛的,以力为基础的收敛提供了收敛量的绝对值,而以位移为基础的收敛仅提供表现收敛的相对量度。一般不单独使用位移收敛准则,否则会产生一定偏差,有些情况会造成假收敛.(ansys非线性分析指南--基本过程Page.6) 。因此ansys官方建议用户尽量以力为基础(或力矩)的收敛误差,如果需要也可以增加以位移为基础的收敛检查。ANSYS缺省是用L2范数控制收敛。其它还有L1范数和L0范数,可用CNVTOL命令设置。在计算中L2值不断变化,若L2

由于ANSYS缺省的criterion计算是你全部变量的平方和开平方(SRSS)*valuse(你设置的值),所以crition也有小小变化。如有需要,也可自己指定crition 为某一常数,CNVTOL,F,10000,0.0001,0

就指定力的收敛控制值为10000*0.0001=1。

另外,非线性计算中用到的一个开关是SOLCONTROL

如关闭SOLCONTROL 选项,那么软件默认收敛准则:力或弯矩的收敛容差是0.001,而不考虑位移的收敛容差;如果打开SOLCONTROL 选项,同

样的默认收敛准则:力或弯矩的收敛容差是0.005,而位移收敛容差是0.05。

非线性收敛非常麻烦,与网格精度、边界条件、荷载步等一系列因素有关,单元的特点对收敛的影响很大,单元的性态不好收敛则困难些;合理的步长可以使求解在真解周围不至于振荡,步长过小,计算量太大,步长过大,会由于过大的荷载步造成不收敛。网格密度适当有助于收敛,网格太密计算量太大,当然太稀计算结果会有较大的误差。究竟多少往往要针对问题进行多次试算。

如果不收敛,可以考虑一下方法改进

1.放松非线性收敛准则。

(CNVTOL #Sets convergence values for nonlinear analyses).

2.增加荷载步数。

(NSUBST #Specifies the number of substeps to be taken this load step)

3.增加每次计算的迭代次数(默认的25次)

(NEQIT #Maximum number of equilibrium iterations allowed each substep) 4 重新划分单元试试,后续会得到不同的答案。

ansys非线性分析指南

ANSYS 非线性分析指南(1) 基本过程 第一章结构静力分析 1. 1 结构分析概述 结构分析的定义: 结构分析是有限元分析方法最常用的一个应用领域。结构这个术语是一个广义的概念,它包括土木工程结构,如桥梁和建筑物;汽车结构,如车身、骨架;海洋结构,如船舶结构;航空结构,如飞机机身、机翼等,同时还包括机械零部件,如活塞传动轴等等。 在ANSYS 产品家族中有七种结构分析的类型,结构分析中计算得出的基 本未知量- 节点自由度,是位移;其他的一些未知量,如应变、应力和反力, 可通过节点位移导出。 七种结构分析的类型分别是: a. 静力分析- 用于求解静力载荷作用下结构的位移和应力等。静力分析 包括线性和非线性分析。而非线性分析涉及塑性、应力刚化、大变形、大应变、超弹性、接触面和蠕变,等。 b. 模态分析- 用于计算结构的固有频率和模态。 c. 谐波分析- 用于确定结构在随时间正弦变化的载荷作用下的响应。 d. 瞬态动力分析- 用于计算结构在随时间任意变化的载荷作用下的响应,并且可计及上述提到的静力分析中所有的非线性性质。 e. 谱分析- 是模态分析的应用拓广,用于计算由于响应谱或PSD 输入 随机振动引起的应力和应变。 f. 屈曲分析- 用于计算屈曲载荷和确定屈曲模态,ANSYS 可进行线性特征值和非线性屈曲分析。 g. 显式动力分析- ANSYS/LS-DYNA可用于计算高度非线性动力学和复 杂的接触问题。 除了前面提到的七种分析类型,还有如下特殊的分析应用: ? 断裂力学 ? 复合材料 ? 疲劳分析

? p-Method 结构分析所用的单元:绝大多数的ANSYS 单元类型可用于结构分析。单元类型从简单的杆单元和梁单元一直到较为复杂的层合壳单元和大应变实体单元 1.2 结构线性静力分析 静力分析的定义: 静力分析计算在固定不变的载荷作用下结构的响应。它不考虑惯性和阻尼的影响,如结构受随时间变化载荷的情况。可是静力分析可以计算那些固定不变的惯性载荷对结构的影响,如重力和离心力;以及那些可以近似为等价静力作用的随时间变化载荷,如通常在许多建筑规范中所定义的等价静力风载和地震载荷。 静力分析中的载荷: 静力分析用于计算由那些不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位移、应力、应变和力。固定不变的载荷和响应是一种假定,即假定载荷和结构的响应随时间的变化非常缓慢,静力分析所施加的载荷包括: ? - 外部施加的作用力和压力 ? - 稳态的惯性力如中力和离心力 ? - 位移载荷 ? - 温度载荷 线性静力分析和非线性静力分析 静力分析既可以是线性的也可以是非线性的。非线性静力分析包括所有的非线性类型:大变形、塑性、蠕变、应力刚化、接触、间隙单元、超弹性单元等,本节主要讨论线性静力分析,非线性静力分析在下一节中介绍。 线性静力分析的求解步骤 1 建模 2 施加载荷和边界条件求解 3 结果评价和分析

ansys前后处理的一些技巧

收集的一些ansys前后处理技巧 1.ANSYS后处理时如何按灰度输出云图? 1)你可以到utilitymenu-plotctrls-style-colors-window colors试试 2)直接utilitymenu-plotctrls-redirect plots 2 将云图输出为JPG 菜单->PlotCtrls->Redirect Plots->To JPEG Files 3.怎么在计算结果实体云图中切面? 命令流 /cplane /type 图形界面操作 <1.设置工作面为切面 <2.PlotCtrls-->Style-->Hidden line Options 将[/TYPE]选项选为section 将[/CPLANE]选项选为working plane 4.非线性计算过程中收敛曲线实时显示 solution>load step opts>output ctrls>grph solu track>on 5.运用命令流进行计算时,一个良好的习惯是: 使用SELECT COMMEND后.........其后再加上ALLSEL......... 6.应力图中左侧的文字中,SMX与SMN分别代表最大值和最小值 如你plnsolv,s,eqv 则 SMX与SMN分别代表最大值等效应力和最小值等效应力 如你要看的是plnsolv,u 则SMX与SMN分别代表位移最大值和位移最小值 不要被S迷惑 mx(max) mn(min) 7.在非线性分析中,如何根据ansys的跟踪显示来判断收敛? 在ansys output windows 有 force convergenge valu 值和 criterion 值当前者小于后者时,就完成一次收敛 你自己可以查看 两条线的意思分别是: F L2:不平衡力的2范数 F CRIT:不平衡力的收敛容差, 如果前者大于后者说明没有收敛,要继续计算 当然如果你以弯矩M为收敛准则那么就对应 M L2 和 M CRIT 希望你现在能明白 8.两个单元建成公共节点,就成了刚性连接,不是接触问题了。做为接触问题,两个互相接触的单元的节点必须是不同的。 9.接触单元 主要分为有厚度和无厚度的,有厚度主要以desai 为代表,无厚度的则以goodman 为代表。尽管古得曼也提出了相应的本构关系,但是如今goodman 单元成了无厚度接触单元的代名词,相应的本构关系现在也作了较大的改进。

ANSYS求解非线性问题

ANSYS求解非线性问题 牛顿一拉森方法 ANSYS程序的方程求解器计算一系列的联立线性方程来预测工程系统的响应。然而,非线性结构的行为不能直接用这样一系列的线性方程表示。需要一系列的带校正的线性近似来求解非线性问题。 逐步递增载荷和平衡迭代 一种近似的非线性救求解是将载荷分成一系列的载荷增量。可以在几个载荷步内或者在一个载步的几个子步内施加载荷增量。在每一个增量的求解完成后,继续进行下一个载荷增量之前程序调整刚度矩阵以反映结构刚度的非线性变化。遗憾的是,纯粹的增量近似不可避免地随着每一个载荷增量积累误差,导种结果最终失去平衡,如图1所示所示。 (a)纯粹增量式解(b)全牛顿-拉普森迭代求解 图1 纯粹增量近似与牛顿-拉普森近似的关系 ANSYS程序通过使用牛顿-拉普森平衡迭代克服了这种困难,它迫使在每一个载荷增量的末端解达到平衡收敛(在某个容限范围内)。图1(b)描述了在单自由度非线性分析中牛顿-拉普森平衡迭代的使用。在每次求解前,NR方法估算出残差矢量,这个矢量是回复力(对应于单元应力的载荷)和所加载荷的差值。程序然后使用非平衡载荷进行线性求解,且核查收敛性。如果不满足收敛准则,重新估算非平衡载荷,修改刚度矩阵,获得新解。持续这种迭代过程直到问题收敛。 ANSYS程序提供了一系列命令来增强问题的收敛性,如自适应下降,线性搜索,自动载荷步,及二分等,可被激活来加强问题的收敛性,如果不能得到收敛,那么程序或者继续计算下一个载荷前或者终止(依据你的指示)。 对某些物理意义上不稳定系统的非线性静态分析,如果你仅仅使用NR方法,正切刚度矩阵可能变为降秩短阵,导致严重的收敛问题。这样的情况包括独立实体从固定表面分离的静态接触分析,结构或者完全崩溃或者“突然变成”另一个稳定形状的非线性弯曲问题。对这样的情况,你可以激活另外一种迭代方法,弧长方法,来帮助稳定求解。弧长方法导致NR平衡迭代沿一段弧收敛,从而即使当正切刚度矩阵的倾斜为零或负值时,也往往阻止发散。这种迭代方法以

ansys学习-非线性静态分析实例

ansys学习-非线性静态分析实例 问题描述 一个子弹以给定的速度射向壁面。壁面假定是刚性的和无摩擦的。将研究子弹和壁面接触后达80微秒长的现象。目的是确定子弹的整个变形,速度历程,以及最大等效Von Mises应变。求解使用SI单位。 用轴对称单元模拟棒。求解最好能通过单一载荷步实现。在这个载荷步中,将同时施加初始速度和约束。将圆柱体末端的节点Y方向约束住以模拟一固壁面。打开自动时间分步来允许ANSYS 确定时间步长。定义分析结束的时间为8E-5秒,以确保有足够长的时间来扑捉整个变形过程。 问题详细说明 下列材料性质应用于这个问题: EX=117.0E09 (杨氏模量) DENS=8930.0 (密度) NUXY=0.35(泊松比) Yield Strength=400.0OE06(屈服强度) Tangent Modulus (剪切模量) 下列尺寸应用于这个问题: 长=32.4E-3m 直径=6.4E-3m 对于这个问题的初始速度是227.0。 图1铜圆柱体图解 求解步骤: 步骤一:设置分析标题 1、选择菜单路径:Utility Menn>File>ChangeTitle。

2、键入文字“Coppery Cylinder Impacting a Rigid Wall” 3、单击OK。 步骤二:定义单元类型 1、选择菜单路径Mail Menu>Preprocessor>Element Type>All/Edit/Delete。 2、单击Add。Library of Element Types(单元类型库)对话框出现。 3、在靠近左边的列表中,单击“Visio Solid”仅一次。 4、选靠近右边的列表中,单击“4node Plas 106”仅一次。 5、单击OK。Library of Element Types 对话框关闭。 6、单击Options (选项)。VISCO106 element type Options(visco106单元类型选项)对话框出现。 7、在关于element behavior(单元特性)的卷动柜中,卷动到“Axisymmetric” 且选中它。 8、单击OK。 9、单击Element Types (单元类型)对话框中的Close。 步骤三:定义材料性质 1、选择菜单路径Main Menu>Preprocessor>Material Props>-Constant-Isotropic. Isotropic Matersal Properties (各向同性材料性质)对话框出现。 2、单击OK来指定材料号为1。另一个I sotropic Material Properties对话框出现。 3、对杨氏模量(EX)键入117.0E09 4、对密度(DENS)键入8930。 5、对泊松比(NUXY)键入0.35。 6、单击OK。 步骤四:定义双线性各向同性强化数据表(BISO) 1、选择菜单路径Main Menu>Preprocessor>Matersal Props>Data Tables> Define/Activate . Define/Activate Data Table(定义数据表)对话柜出现。 2、在关于type of data table(数据表类型)的卷动框中,卷动到“Bilin isotr BISO”且选中它。 3、对material reference number(材料参考号)健入1。 4、对number of temperatures(温度数)键入1和单击OK。 5、选择菜单路径Main Menu>Preprocessor>Material Props>Data Tables>Edit Active. Data Table BISO对话框出现。 6、对YLD Strs(屈服应力)键入400.0e06。 7、对 Tang Mod(剪切模量)键入100.0e06。 8、选择File>Apply & Quit。 9、选择菜单路径Main Menu>Preprosessor>Material Porps>Data Tables>Graph. Graph Data Tables(图形表示数据表)对话框出现。 10、单击OK接受绘制BISO表的缺省。一个BISO表的标绘图出现在ANSYS图形窗口中。 11、在ANSYS TooLbar上单击SAVE_DB。 步骤五、产生矩形 在这一步中,你产生一个代表柱体半横截面积的矩形。

ansys收敛问题

1.一般首先是改变初值,尝试不同的初始化,事实上好像初始化很关键,对于收敛~ 2.Ansys的收敛最基础的是网格的质量,这个靠经验 3.首先查找网格问题,如果问题复杂比如与模型、边界、初始条件都有关系。 4.边界条件、网格质量 5.有时初始条件和边界条件严重影响收敛性,曾经作过一个计算反反复复,通过修改网格,重新定义初始条件,包括具体的选择的模型,还有老师经常用的方法就是看看哪个因素不收敛,然后寻找和它有关的条件,改变相应参数。就收敛了 6.A.检查是否哪里设定有误.B.从算至发散前几步,看presure分布,看不出来的话,再算几步, 看看问题大概出在那个区域,连地方都知道的话,应该不难想出问题所在. C.网格,配合第二点作修正, 或是认命点,就重建个更漂亮的,或是更粗略的来除错... D.再找不出来的话,我会换个solver... 7.我解决的办法是设几个监测点,比如参数变化较大的地方,若这些地方的参数变化很小,就可以认为是收敛了,尽管此时残值曲线还没有降下来。 8.记得好像调节松弛因子也能影响收敛,不过代价是收敛速度。 9.网格有一定的影响,最主要的还是初始和边界条件 载荷步、载荷子步均是对所施加荷载的一种描述方式。在施加荷载的时候需要对载荷步、载荷子步进行定义。 载荷步仅仅是为了获得解答的载荷配置,它的作用是在给定时间间隔内的一组荷载。在线性静态或稳态分析中,可以使用不同的载荷步、施加不同的载荷组合。在瞬态分析中,多个载荷步载荷历程曲线的不同区段来描述荷载随时间的变化情况。 在有一些分析中需要用到载荷子步。载荷子步是正在求解的载荷步中的时间点,是对载荷步描述的进一步细化。 在所有的静态和瞬态分析中,ANSYS通过指定分析中载荷步结束的时间来定义载荷步。这样,在瞬态分析或其他有关速率的静态分析中,时间具有实际意义;在于速率无关的分析中,时间是作为识别载荷步以及载荷子步的“计数器”来跟踪载荷步,并无实际意义。 时间步则一般在非线性问题或瞬态动力学问题求解中使用,是每一次迭代求解的步长,设置过大则容易不收敛。非线性求解被分成三个操作级别:载荷步、子步、平衡迭代。“顶层”级别由在一定“时间”范围内你明确定义的载荷步组成。假定载荷在载荷步内是线性地变化的。在每一个载荷是步内,为了逐步加载可以控制程序来执行多次求解(子步或时间步)。在每一个子步内,程序将进行一系列的平衡迭代以获得收敛的解。小的时间步通常导致较好的精度,但这是以增多的运行时间为代价的。 ANSYS提供两种方法来控制子步数:直接设定子步数或时间步长或自动时间步长。采用自动时间步长时,ANSYS程序,基于结构的特性和系统的响应,来自动调整每一子步的时间步长。如果结构的行为从线性变化到非线性,或者想要在系统响应的非线性部分期间变化时间步长,可以激活自动时间步长以便随需要调整时间步长,获得精度和计算代价之间的良好平衡。如果不能确保问题成功地收敛,则可以使用自动时间分布来激活ANSYS程序的二分法。二分法是对收敛失败自动矫正的方法。应用二分法时,只要平衡迭代不收敛,就把时间步长减半,然后从最后收敛的子步自动重新计算。如果已二分的时间步再次收敛失败,程序将再次分割时间步长然后重新计算,持续过程一直到获得收敛或者到达设置的最小时间步长

Ansys非线性接触分析和设置

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下非线性分析不收敛

解决非线性分析不收敛的技巧 影响非线性收敛稳定性及其速度的因素很多: 1、模型——主要是结构刚度的大小。对于某些结构,从概念的角度看,可以认为它是几何不变的稳定体系。但如果结构相近的几个主要构件刚度相差悬殊,在数值计算中就可能导致数值计算的较大误差,严重的可能会导致结构的几何可变性——忽略小刚度构件的刚度贡献。 如出现上述的结构,要分析它,就得降低刚度很大的构件单元的刚度,可以加细网格划分,或着改用高阶单元(BEAM->SHELL,SHELL->SOLID)。构件的连接形式(刚接或铰接)等也可能影响到结构的刚度。 2、线性算法(求解器)。ANSYS中的非线性算法主要有:稀疏矩阵法(SPARSE DIRECT SOLVER)、预共轭梯度法(PCG SOLVER)和波前法(FRONT DIRECT SLOVER)。稀疏矩阵法是性能很强大的算法,一般默认即为稀疏矩阵法(除了子结构计算默认波前法外)。预共轭梯度法对于3-D实体结构而言是最优的算法,但当结构刚度呈现病态时,迭代不易收敛。为此推荐以下算法: 1)、BEAM单元结构,SHELL单元结构,或以此为主的含3-D SOLID的结构,用稀疏矩阵法; 2)、3-D SOLID的结构,用预共轭梯度法; 3)、当你的结构可能出现病态时,用稀疏矩阵法; 4)、当你不知道用什么时,可用稀疏矩阵法。 3、非线性逼近技术。在ANSYS里还是牛顿-拉普森法和弧长法。牛顿-拉普森法是常用的方法,收敛速度较快,但也和结构特点和步长有关。弧长法常被某些人推崇备至,它能算出力加载和位移加载下的响应峰值和下降响应曲线。但也发现:在峰值点,弧长法仍可能失效,甚至在非线性计算的线性阶段,它也可能会无法收敛。 为此,尽量不要从开始即激活弧长法,还是让程序自己激活为好(否则出现莫名其妙的问题)。子步(时间步)的步长还是应适当,自动时间步长也是很有必要的。 4、加快计算速度 在大规模结构计算中,计算速度是一个非常重要的问题。下面就如何提高计算速度作一些建议: 充分利用ANSYS MAP分网和SWEEP分网技术,尽可能获得六面体网格,这一方面减小解题规模,另一方面提高计算精度。 在生成四面体网格时,用四面体单元而不要用退化的四面体单元。比如95号单元有20节点,可以退化为10节点四面体单元,而92号单元为10节点单元,在此情况下用92号单元将优于95号单元。 选择正确的求解器。对大规模问题,建议采用PCG法。此法比波前法计算速度要快10倍以上(前提是您的计算机内存较大)。对于工程问题,可将ANSYS缺省的求解精度从1E-8改为1E-4或1E-5即可。 5、荷载步的设置直接影响到收敛。应该注意以下几点: 1、设置足够大的荷载步(将MAXMIUM SUBSTEP=1000000),可以更容易收敛,避免发散的出现(nsub,nsbstp,nsbmx,nsbmn);

ANSYS中文翻译官方手册_接触分析

一般的接触分类 (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使用的接触单元和使用它们的过程,下面分类详述。 点─点接触单元 点─点接触单元主要用于模拟点─点的接触行为,为了使用点─点的接触单元,你需要预先知道接触位置,这类接触问题只能适用于接触面之间有较小相对滑动的情况(即使在几何非线性情况下) 如果两个面上的结点一一对应,相对滑动又以忽略不计,两个面挠度(转动)保持小量,那么可以用点─点的接触单元来求解面─面的接触问题,过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。 点─面接触单元 点─面接触单元主要用于给点─面的接触行为建模,例如两根梁的相互接触。 如果通过一组结点来定义接触面,生成多个单元,那么可以通过点─面的接触单元来模拟面─面的接触问题,面即可以是刚性体也可以是柔性体,这类接触问题的一个典型例子是插头到插座里。

ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析

!ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析 !学习重点: !1、强化非线性屈曲知识 首先了解屈曲问题。在理想化情况下,当F < Fcr时, 结构处于稳定平衡状态,若引入一个小的侧向扰动力,然后卸载, 结构将返回到它的初始位置。当F > Fcr时, 结构处于不稳定平衡状态, 任何扰动力将引起坍塌。当F = Fcr时,结构处于中性平衡状态,把这个力定义为临界载荷。在实际结构中, 几何缺陷的存在或力的扰动将决定载荷路径的方向。在实际结构中, 很难达到临界载荷,因为扰动和非线性行为, 低于临界载荷时结构通常变得不稳定。 要理解非线性屈曲分析,首先要了解特征值屈曲。特征值屈曲分析预测一个理想线弹性结构的理论屈曲强度,缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲强度,特征值屈曲一般产生非保守解, 使用时应谨慎。 !理论解,根据Euler公式。其中μ取决于固定方式。 !有限元方法, 已知在特征值屈曲问题: 求解,即可得到临界载荷 而非线性屈曲问题: 其中为结构初始刚度,为有缺陷的结构刚度,为位移矩阵,为载荷矩阵。 非线性屈曲分析时考虑结构平衡受扰动(初始缺陷、载荷扰动)的非线性静力分析,该分析时一直加载到结构极限承载状态的全过程分析,分析中可以综合考虑材料塑性、几何非线性、接触、大变形。非线性屈曲比特征值屈曲更精确,因此推荐用于设计或结构的评价。 !2、熟悉WB中非线性屈曲分析流程 (1) 前处理,施加单元载荷,进行预应力静力分析。 (2) 基于预应力静力分析,指定分析类型为特征值屈曲分析,完成特征值屈曲分析。 (3) 在APDL模块将一阶特征屈曲模态位移乘以适当系数,将此变形后的形状当做非线性分析的初始模型。

Ansys使用技巧-非线性收敛准则

ansys计算非线性时会绘出收敛图,其中横坐标是cumulative iteration number 纵坐标是absolute convergence norm。他们分别是累积迭代次数和绝对收敛范数,用来判断非线性分析是否收敛。 ansys在每荷载步的迭代中计算非线性的收敛判别准则和计算残差。其中计算残差是所有单元内力的范数,只有当残差小于准则时,非线性叠代才算收敛。ansys的位移收敛是基于力的收敛的,以力为基础的收敛提供了收敛量的绝对值,而以位移为基础的收敛仅提供表现收敛的相对量度。一般不单独使用位移收敛准则,否则会产生一定偏差,有些情况会造成假收敛.(ansys非线性分析指南--基本过程Page.6) 。因此ansys官方建议用户尽量以力为基础(或力矩)的收敛误差,如果需要也可以增加以位移为基础的收敛检查。ANSYS缺省是用L2范数控制收敛。其它还有L1范数和L0范数,可用CNVTOL命令设置。在计算中L2值不断变化,若L2

非线性-接触分析

接触分析 接触问题是一种高度非线性行为,需要较大的计算资源,为了进行实为有效的计算,理解问题的特性和建立合理的模型是很重要的。 接触问题存在两个较大的难点:其一,在你求解问题之前,你不知道接触区域,表面之间是接触或分开是未知的,突然变化的,这随载荷、材料、边界条件和其它因素而定;其二,大多的接触问题需要计算摩擦,有几种摩擦和模型供你挑选,它们都是非线性的,摩擦使问题的收敛性变得困难。 一般的接触分类 接触问题分为两种基本类型:刚体─柔体的接触,半柔体─柔体的接触,在刚体─柔体的接触问题中,接触面的一个或多个被当作刚体,(与它接触的变形体相比,有大得多的刚度),一般情况下,一种软材料和一种硬材料接触时,问题可以被假定为刚体─柔体的接触,许多金属成形问题归为此类接触,另一类,柔体─柔体的接触,是一种更普遍的类型,在这种情况下,两个接触体都是变形体(有近似的刚度)。 ANSYS接触能力 ANSYS支持三种接触方式:点─点,点─面,平面─面,每种接触方式使用的接触单元适用于某类问题。 为了给接触问题建模,首先必须认识到模型中的哪些部分可能会相互接触,如果相互作用的其中之一是一点,模型的对立应组元是一个结点。如果相互作用的其中之一是一个面,模型的对应组元是单元,例如梁单元,壳单元或实体单元,有限元模型通过指定的接触单元来识别可能的接触匹对,接触单元是覆盖在分析模型接触面之上的一层单元,至于ANSTS 使用的接触单元和使用它们的过程,下面分类详述。 点─点接触单元 点─点接触单元主要用于模拟点─点的接触行为,为了使用点─点的接触单元,你需要预先知道接触位置,这类接触问题只能适用于接触面之间有较小相对滑动的情况(即使在几何非线性情况下) 如果两个面上的结点一一对应,相对滑动又以忽略不计,两个面挠度(转动)保持小量,那么可以用点─点的接触单元来求解面─面的接触问题,过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。 点─面接触单元 点─面接触单元主要用于给点─面的接触行为建模,例如两根梁的相互接触。 如果通过一组结点来定义接触面,生成多个单元,那么可以通过点─面的接触单元来模拟面─面的接触问题,面即可以是刚性体也可以是柔性体,这类接触问题的一个典型例子是插头到插座里。 使用这类接触单元,不需要预先知道确切的接触位置,接触面之间也不需要保持一致的网格,并且允许有大的变形和大的相对滑动。 Contact48和Contact49都是点─面的接触单元,Contact26用来模拟柔性点─刚性面的接触,对有不连续的刚性面的问题,不推荐采用Contact26因为可能导致接触的丢失,在这种情况下,Contact48通过使用伪单元算法能提供较好的建模能力。 面─面的接触单元 ANSYS支持刚体─柔体的面─面的接触单元,刚性面被当作“目标”面,分别用Targe169和Targe170来模拟2─D和3—D的“目标”面,柔性体的表面被当作“接触”面,用Conta171,Conta172,Conta173,Conta174来模拟。一个目标单元和一个接单元叫作一个“接触对”程序通过一个共享的实常号来识别“接触对”,为了建立一个“接触对”给目标单元和接触单元指定相同的实常的号。 与点─面接触单元相比,面─面接触单元有好几项优点, *支持低阶和高阶单元 *支持有大滑动和摩擦的大变形,协调刚度阵计算,单元提法不对称刚度阵的选项。 *提供工程目的采用的更好的接触结果,例如法向压力和摩擦应力。 *没有刚体表面形状的限制,刚体表面的光滑性不是必须允许有自然的或网格离散引起的表面不连续。

ANSYS 非线性_结构分析

目录 非线性结构分析的定义 (1) 非线性行为的原因 (1) 非线性分析的重要信息 (3) 非线性分析中使用的命令 (8) 非线性分析步骤综述 (8) 第一步:建模 (9) 第二步:加载且得到解 (9) 第三步:考察结果 (16) 非线性分析例题(GUI方法) (20) 第一步:设置分析标题 (21) 第二步:定义单元类型 (21) 第三步:定义材料性质 (22) 第四步:定义双线性各向同性强化数据表 (22) 第五步:产生矩形 (22) 1

第六步:设置单元尺寸 (23) 第七步:划分网格 (23) 第八步:定义分析类型和选项 (23) 第九步:定义初始速度 (24) 第十步:施加约束 (24) 第十一步:设置载荷步选项 (24) 第十二步:求解 (25) 第十三步:确定柱体的应变 (25) 第十四步:画等值线 (26) 第十五步:用Post26定义变量 (26) 第十六步:计算随时间变化的速度 (26) 非线性分析例题(命令流方法) (27) 非线性结构分析 非线性结构的定义 在日常生活中,会经常遇到结构非线性。例如,无论何时用钉书针钉书,金 2

属钉书钉将永久地弯曲成一个不同的形状。(看图1─1(a))如果你在一个木架上放置重物,随着时间的迁移它将越来越下垂。(看图1─1(b))。当在 汽车或卡车上装货时,它的轮胎和下面路面间接触将随货物重量的啬而变化。(看图1─1(c))如果将上面例子所载荷变形曲线画出来,你将发现它们都显示了非线性结构的基本特征--变化的结构刚性. 图1─1 非线性结构行为的普通例子 3

非线性行为的原因 引起结构非线性的原因很多,它可以被分成三种主要类型: 状态变化(包括接触) 许多普通结构的表现出一种与状态相关的非线性行为,例如,一根只能拉伸的电缆可能是松散的,也可能是绷紧的。轴承套可能是接触的,也可能是不接触的, 冻土可能是冻结的,也可能是融化的。这些系统的刚度由于系统状态的改变在不同的值之间突然变化。状态改变也许和载荷直接有关(如在电缆情况中),也可能由某种外部原因引起(如在冻土中的紊乱热力学条件)。ANSYS程序中单元的激活与杀死选项用来给这种状态的变化建模。 接触是一种很普遍的非线性行为,接触是状态变化非线性类型形中一个特殊而重要的子集。 几何非线性 如果结构经受大变形,它变化的几何形状可能会引起结构的非线性地响应。一个例的垂向刚性)。随着垂向载荷的增加,杆不断弯曲以致于动力臂明显地减少,导致杆端显示出在较高载荷下不断增长的刚性。 4

ANSYS的非线性收敛准则

ANSYS的非线性收敛准则 ansys非线性问题 ANSYS的非线性收敛准则--转自中华钢结构论坛 CNVTOL, Lab, VALUE, TOLER, NORM, MINREF ANSYS中,非线性收敛准则主要有力的收敛,位移的收敛,弯矩的收敛和转角的收敛。一般用力的控制加载时,可以使用残余力的2-范数控制收敛;而位移控制加载时,最好用位移的范数控制收敛。 When SOLCONTROL,ON, TOLER Defaults to 0.005 (0.5%) for force and moment, and 0.05 (5%) for displacement when rotational DOFs are not present. When SOLCONTROL,OFF, defaults to 0.001 (0.1%) for force and moment. 收敛精度一般可放宽至 5%,以提高收敛速度。 加快收敛的方法有一下几种: 1可以增大荷载子步数,nsubst,nsbstp,nsbmn,carry 2修改收敛准则,cnvtol,lab,value,toler,norm,minref 3 打开优化的非线性默认求解设置和某些强化的内部求解算法,solcontrol,key1,key2,key3,vtol(一般情况下,默认是打开的) 4重新划分网格,网格的单元不宜太大或太小, 一般在5~10厘米左右 5 检查模型的正确性 下面计算收敛过程图中的各个曲线的具体含义是什么? 非线性计算是一个迭代计算的过程,曲线表示两次迭代之间的误差,图中分别表示力和位移在迭代过程中的每次迭代之间的误差 关于ansys中收敛准则(cnvtol)理解 ansys中依据缺省的收敛准则,程序将对不平衡力SRSS与VALUE*TOLER的值进行比较;而VALUE的缺省值是在SRSS和MINREF中取较大值。现假如TOLER的缺省值是0.1的话,这个准则是不是可以理解成后一次的SRSS是前一次的SRSS的01倍就收敛啦? 请指点 我是这样理解的例如下面的命令流: cnvtol,f,5000,0.0005,0 cnvtol,u,10,0.001,2

关于ansys中收敛的介绍

一.何为收敛?在这里我引用一个会员的提问来解释这个问题: Q:结构非线性静力分析经常出现收敛这个词,如:收敛容限,收敛准则,收敛的解,位移收敛检验等,请解释,thanks! A: 个人是这样理解的: 谈到收敛总会和稳定性联系在一起, 简单的说,就是在进行求解过程中的一些中间值的误差对于结果的影响的大小,当中间量的误差对于你的数值积分的结果没有产生影响,就说明你的积分方法是稳定的,最终你的 数值积分的结果就会收敛于精确解;当中间量的误差导致数值积分结果与精确解有很大的差别时,就说明你的方法稳定性不好,你的数值积分结果不会收敛于精确解。 我想当你对于稳定性和收敛的概念真正理解后,那些名词对于你来说,并不是问题,力学的问题最终都会和数学联系在一起,建议你看看数值积分方面的教程,学好了数学,力学对于你来说就是a piece of cake。 Q:那么说收不收敛,最终都是因为采用的计算方法和计算参数选取的问题了? A: 就本人所学的专业来说,很大程度上取决于所采用的算法,我学的是结构工程,举个例子吧 :当在进行结构动力时程分析时,采用的几分方法有线性加速度法,威尔逊-theta法,对于线性加速度法,当时间步长大于周期的0.5倍时,计算结果很可能出现不收敛,而当时间步长小于0.1倍的周期时,才有可能获得稳定的计算结果;而威尔逊-theta法,实质上就是线性加速度法的修正形式,很多实例表明当theta值大于1.37时,这种算法是无条件稳定的。当然影响计算结果是否收敛的原因有很多,比如初始条件,我所指的仅仅是我所学专业的一个问题的很小的一个方面。 A: 说白了,就是数学。 牵涉到实际的计算问题时,才发现数学实在是太有用了,不过可惜数学实在学得不好。

ansys非线性收敛总结(转载)

ansys非线性收敛总结 智创仿真 2016年8月6日1750 文章来源于网络,讲解很系统,可以经典收藏,由于无法查证出处,无意冒犯,如有不 妥,请联系我 ansys非线性收敛总结 ansys计算非线性时会绘出收敛图,其中横坐标是cumulative iteration number 纵坐标是absolute convergence norm。他们分别是累积迭代次数和绝对收敛范数,用来判断非线性分析是否收敛。 ansys在每荷载步的迭代中计算非线性的收敛判别准则和计算残差。其中计算残差是所有单元内力的范数,只有当残差小于准则时,非线性叠代才算收敛。 ansys的收敛是基于力的收敛的,以力为基础的收敛提供了收敛量的绝对值,而以位移为基础的收敛仅提供表现收敛的相对量度。一般不单独使用位移收敛准则,否则会产生一定偏差,有些情况会造成假收敛.(ansys非线性分析指南--基本过程Page.6) 。因此ansys官方建议用户尽量以力为基础(或力矩)的收敛误差,如果需要也可以增加以位移为基础的收敛检查。 ANSYS缺省是用L2范数控制收敛。其它还有L1范数和L0范数,可用CNVTOL命令设置。在计算中L2值不断变化,若L2

关于ansys非线性分析的几点忠告

关于非线性分析的几点忠告 了解程序的运作方式和结构的表现行为 如果你以前没有使用过某一种特别的非线性特性,在将它用于大的,复杂的模型前,构造一个非常简单的 模型(也就是,仅包含少量单元),以及确保你理解了如何处理这种特性。 通过首先分析一个简化模型,以便使你对结构的特性有一个初步了解。对于非线性静态模型,一个初步的 线性静态分析可以使你知道模型的哪一个区域将首先经历非线性响应,以及在什么载荷范围这些非线性将 开始起作用。对于非线性瞬态分析,一个对梁,质量块及弹簧的初步模拟可以使你用最小的代价对结构的 动态有一个深入了解。在你着手最终的非线性瞬时动态分析前,初步非线性静态,线性瞬时动态,和/或模 态分析同样地可以有助于你理解你结构的非线性动态响应的不同的方面。 阅读和理解程序的输出信息和警告。至少,在你尝试后处理你的结果前,确保你的问题收敛。对于与路程 相关的问题,打印输出的平衡迭代记录在帮助你确定你的结果是有效还是无效方面是特别重的。 简化 尽可能简化最终模型。如果可以将3─D结构表示为2─D平面应力,平面应变或轴对称模型,那么这样做, 如果可以通过对称或反对称表面的使用缩减你的模型尺寸,那么这样做。(然而,如果你的模型非对称加 载,通常你不可以利用反对称来缩减非线性模型的大小。由于大位移,反对称变成不可用的。)如果你可 以忽略某个非线性细节而不影响你模型的关键区域的结果,那么这样做。 只要有可能就依照静态等效载荷模拟瞬时动态加载。 考虑对模型的线性部分建立子结构以降低中间载荷或时间增量及平衡迭代所需要的计算时间。 采用足够的网格密度 考虑到经受塑性变形的区域要求一个合理的积分点密度。每个低阶单元将提供和高阶单元所能提供的一样

相关文档
最新文档