解决ansys下非线性分析不收敛
ANSYS结构非线性分析指南
ANSYS结构非线性分析指南ANSYS是一个强大的工程仿真软件，能够对各种复杂的结构进行分析。
其中，结构非线性分析是其中一种重要的分析方法，它能够模拟结构在非线性载荷和变形条件下的行为。
本文将为您提供一个ANSYS结构非线性分析的指南，帮助您更好地理解和应用这个方法。
首先，我们需要明确结构非线性分析的目标。
一般来说，结构非线性分析主要用于研究结构在大变形、材料非线性、接触或摩擦等复杂条件下的响应。
例如，当结构受到极大的外力作用时，其产生的变形可能会导致材料的非线性行为，这时我们就需要进行非线性分析。
在进行非线性分析之前，我们需要进行准备工作。
首先，我们需要准备一个几何模型，可以通过CAD软件导入或者直接在ANSYS中绘制。
然后，我们需要选择合适的材料模型，这将直接影响分析结果的准确性。
ANSYS提供了多种材料模型，例如线弹性模型、塑性模型和粘弹性模型等。
接下来，我们需要定义边界条件和载荷。
边界条件指明了结构的固定边界和自由边界，这决定了结构的位移约束。
载荷是作用在结构上的外力或者外界约束，例如压力、点载荷或者摩擦力等。
在非线性分析中，载荷的大小和施加方式可能会导致结构的非线性响应，因此需要仔细选择。
接下来，我们需要选择适当的非线性分析方法。
ANSYS提供了多种非线性分析方法，例如几何非线性分析、材料非线性分析和接触非线性分析等。
几何非线性分析适用于大变形情况下的分析，材料非线性分析适用于材料的弹塑性行为分析，而接触非线性分析适用于多个结构之间的接触行为分析。
在进行非线性分析之前，我们需要对模型进行预处理，包括网格划分和解算控制参数的设置。
网格划分的精度会直接影响分析结果的准确性，因此需要进行适当的剖分。
解算控制参数的设置涉及到收敛性和稳定性的问题，需要进行合理的调整。
然后，我们可以进行非线性分析了。
ANSYS提供了多种求解器，例如Newton-Raphson方法和弧长法等。
这些求解器可以通过迭代算法来求解非线性方程组，得到结构的响应结果。
关于ansys中收敛的介绍
5、荷载步的设置直接影响到收敛。应该注意以下几点：
1、设置足够大的荷载步（将MAXMIUM SUBSTEP=1000000），可以更容易收敛，避免发散的出现(nsub,nsbstp,nsbmx,nsbmn)；
如出现上述的结构，要分析它，就得降低刚度很大的构件单元的刚度，可以加细网格划分，或着改用高阶单元（BEAM->SHELL,SHELL->SOLID)。构件的连接形式（刚接或铰接）等也可能影响到结构的刚度。
2、线性算法（求解器）。ANSYS中的非线性算法主要有：稀疏矩阵法(SPARSE DIRECT SOLVER)、预共轭梯度法(PCG SOLVER)和波前法(FRONT DIRECT SLOVER)。稀疏矩阵法是性能很强大的算法，一般默认即为稀疏矩阵法（除了子结构计算默认波前法外）。预共轭梯度法对于3-D实体结构而言是最优的算法，但当结构刚度呈现病态时，迭代不易收敛。为此推荐以下算法：
4重新划分网格,网格的单元不宜太大或太小, 一般在5～10厘米左右
5 检查模型的正确性
五.计算收敛过程图中的各个曲线的具体含义是什么？
非线性计算是一个迭代计算的过程，曲线表示两次迭代之间的误差，图中分别表示力和位移在迭代过程中的每次迭代之间的误差
关于ansys中收敛准则(cnvtol)理解
4、加快计算速度
在大规模结构计算中，计算速度是一个非常重要的问题。下面就如何提高计算速度作一些建议：
充分利用ANSYS MAP分网和SWEEP分网技术，尽可能获得六面体网格，这一方面减小解题规模，另一方面提高计算精度。
ANSYS的非线性收敛准则
ANSYS的非线性收敛准则ansys非线性问题ANSYS的非线性收敛准则--转自中华钢结构论坛CNVTOL, Lab, VALUE, TOLER, NORM, MINREFANSYS中,非线性收敛准则主要有力的收敛，位移的收敛，弯矩的收敛和转角的收敛。
一般用力的控制加载时，可以使用残余力的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,carry2修改收敛准则,cnvtol,lab,value,toler,norm,minref3 打开优化的非线性默认求解设置和某些强化的内部求解算法，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,0cnvtol,u,10,0.001,2如果不平衡力（独立的检查每一个自由度）小于等于5000*0.0005（也就是2.5），并且如果位移的变化小于等于10*0.001时，认为子步是收敛的。
ANSYS材料非线性分析
【分享】ANSYS7.0超弹材料的定义－新的曲线拟合功能--摘自ansys用户专区几何非线性几何非线性不受敛主要原因1.网格质量，特别是warpage2.约束方程，少用刚性连接3.收敛准则，可适当加大容差4.荷载步设置，可适当加大步数最近碰到一个对我来说很意外的问题：如果确实如此希望大家以后小心大家知道定义接触后会自动生成一组实常数，前几天我碰到一个问题，需定义超过10组实常数，接触对很多，好像有20多处，按照常规步骤划分完所有网格，当时因为有一个实常数参数没确定，便预留了最后一组（第10组）实常数里面的参数为空，接下来就定义了所有的接触对，由于所有接触对里的设置一样，ANSYS在我保存db完重新打开后便把我所有的接触对综合成一个了！接下来我就把第十组实常数里面的参数补上了，但在求解时却提示我该实常数同时被两种单元（包括CNTACT单元）同时占用，出现错误！！检查了半天才发现自动生成的接触对实常数把第10组实常数也占用了！我实在没找到什么好的解决办法，只得把接触对删除了重新定义，那可是上百多个面的选取过程，痛苦不堪简直！ANSYS里接触对面的选取时还不能针对Component操作！ANSYS7.0超弹材料的定义－新的曲线拟合功能ANSYS7.0中的超弹材料模拟能力得到了很大的加强，在ANSYS6.1的超弹材料模型的基础上又增加了Gent, Yeoh, Blatz-Ko, and Ogden (Foam)四种超弹性材料模型，使得其超弹模拟能力得到了进一步扩展。
ANSYS7.0中对超弹能力最吸引人的增强还不在于此，而是在于其曲线拟合能力的大幅度扩展，不再像ANSYS6.1以前的版本一样曲线拟合仅仅局限于Mooney-Rivlin模型，而是将其扩展到所有的超弹模型，这样，用户可以利用实验得到的应力应变数据直接让程序自己拟合出任意一种超弹材料模型的参数，大大方便了用户的使用。
以下就ANSYS7.0的超弹拟合功能做一简单介绍。
ANSYS讲义非线性分析
t1
t2
时间 t
XJTU
自动时间步(续)
• 自动时间步算法是 非线性求解控制 中包含的多种算法的一种。
（在以后的非线性求解控制中有进一步的讨论。） • 基于前一步的求解历史与问题的本质，自动时间步算法或者增加
或者减小子步的时间步大小。
XJTU
5) 输出文件的信息
在非线性求解过程中，输出窗口显示许多关于收敛的信息。输出 窗口包括：
子步
时间 ”相关联。
“时间
两个载荷步的求解 ”
XJTU
在非线性求解中的 “ 时间 ”
• 每个载荷步与子步都与 “ 时间 ”相关联。 子步 也叫时间步。
• 在率相关分析（蠕变，粘塑性）与瞬态分析中，“ 时间 ”代表真实 的时间。
• 对于率无关的静态分析，“ 时间 ” 表示加载次序。在静态分析中， “ 时间 ” 可设置为任何适当的值。
最终结果偏离平衡。
u 位移
XJTU
1) Newton-Raphson 法
ANSYS 使用Newton-Raphson平衡迭代法 克服了增量
求解的问题。 在每个载荷增量步结束时，平衡迭代驱 使解回到平衡状态。
载荷
F
4 3 2
1
u 位移
一个载荷增量中全 Newton-Raphson 迭代 求解。（四个迭代步如 图所示）
XJTU
非线性分析的应用(续)
宽翼悬臂梁的侧边扭转失 稳
一个由于几何非线性造 成的结构稳定性问题
XJTU
非线性分析的应用(续)
橡胶底密封 一个包含几何非线 性（大应变与大变 形），材料非线性 （橡胶），及状态 非线性（接触的例 子。
XJTU
非线性分析的应用(续)
关于ansys中收敛的介绍
关于ansys中收敛的介绍默认分类2010-09-26 16:25:12 阅读277 评论1 字号：大中小订阅一、何为收敛？在这里我引用一个会员的提问来解释这个问题：Q:结构非线性静力分析经常出现收敛这个词，如：收敛容限，收敛准则，收敛的解，位移收敛检验等，请解释，thanks!A: 个人是这样理解的：谈到收敛总会和稳定性联系在一起，简单的说，就是在进行求解过程中的一些中间值的误差对于结果的影响的大小，当中间量的误差对于你的数值积分的结果没有产生影响，就说明你的积分方法是稳定的，最终你的数值积分的结果就会收敛于精确解；当中间量的误差导致数值积分结果与精确解有很大的差别时，就说明你的方法稳定性不好，你的数值积分结果不会收敛于精确解。
我想当你对于稳定性和收敛的概念真正理解后，那些名词对于你来说，并不是问题，力学的问题最终都会和数学联系在一起，建议你看看数值积分方面的教程，学好了数学，力学对于你来说就是a piece of cake。
Q:那么说收不收敛，最终都是因为采用的计算方法和计算参数选取的问题了？A:就本人所学的专业来说，很大程度上取决于所采用的算法，我学的是结构工程，举个例子吧:当在进行结构动力时程分析时，采用的几种方法有线性加速度法，威尔逊－theta法，对于线性加速度法，当时间步长大于周期的0.5倍时，计算结果很可能出现不收敛，而当时间步长小于0.1倍的周期时，才有可能获得稳定的计算结果；而威尔逊－theta法，实质上就是线性加速度法的修正形式，很多实例表明当theta 值大于1.37时，这种算法是无条件稳定的。
当然影响计算结果是否收敛的原因有很多，比如初始条件，我所指的仅仅是我所学专业的一个问题的很小的一个方面。
A: 说白了，就是数学。
牵涉到实际的计算问题时，才发现数学实在是太有用了，不过可惜数学实在学得不好。
A: 收敛的问题，就好像你往水里扔一块石头激起的波浪，慢慢会平息下来，这就收敛了。
ANSYS非线性接触问题分析汇总
ANSYS 分线性接触问题分析汇总接触非线性是一门复杂的学科，ANSYS 关于计算非线性接触的设置选项多只又多，很多人摸不到头脑，本文就基于ANSYS 模拟过的几个接触实例，研究了相关设置选项对接触结果的影响。
实例1：橡胶密封圈配合接触研究—非线性求解设置对结果的影响密封圈配合模型简图见图1，左右两端为刚体，中间圆部分为橡胶密封圈，将刚体2沿刚体1方面移动，从而实现橡胶圈密封作用，采用plane182单元，设置轴对称行为，建立橡胶密封圈与刚体接触模型，见图2。
图1 密封圈配合模型简图 图2 密封圈配合有限元模型图接触对采用默认设置，摩擦系数取0.10，研究非线性求解器设置对收敛方面的影响，大变形静态（Large Displacement Static ）效应打开，自动时间步长（Automatic time stepping ）打开，子步数（Number of substeps ）设置为50，线性搜索（Line search ）打开。
1 收敛准则对结果的影响此实例收敛准则默认采用力收敛结合力矩收敛准则（基于L2范数），收敛容差（Tolerance ）默认为0.001，工程上认为0.05的收敛容差足够满足要求。
表 1 收敛容差对计算结果的影响收敛容差 最大应力/ MPa报错与否？ 0.001 4.12364报错 0.05 4.12785 报错 0.14.12996报错查看报错信息，见图3，表示单元过于扭曲，建议提高子步数或降低时间步长，需要提高网格质量，也要考虑材料属性，接触对及约束方程的合理性，若在第一步迭代就如此，需要预先执行单元形状检查。
图3 报错信息刚体1刚体2密封圈橡胶密封圈配合Von Mises应力云图见图4。
图4 橡胶密封圈配合Von Mises应力2 子步数对结果的影响此实例子步数设置为50、100、200、500，收敛容差（Tolerance）默认为0.001，研究子步数对收敛的影响。
Ansys非线性接触分析和设置
Ansys非线性接触分析和设置设置实常数和单元关键选项程序使用20个实常数和数个单元关键选项，来控制面─面接触单元的接触。
参见《ANSYS Elements Reference》中对接触单元的描述。
实常数在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 静、动摩擦衰减系数。
命令：RGUI：main menu> preprocessor>real constant对实常数FKN, FTOLN, ICONT, PINB, PMAX, PMIN, FKOP 和FKT，用户既可以定义一个正值，也可以定义一个负值。
程序将正值作为比例因子，将负值作为绝对值。
程序将下伏单元的厚度作为ICON，FTOLN，PINB，PMAX 和PMIN 的参考值。
例如ICON = 表明初始闭合因子是“*下层单元的厚度”。
然而，ICON = 则表示真实调整带是单位。
如果下伏单元是超单元，则将接触单元的最小长度作为厚度。
参见图5-8。
图5-8 下层单元的厚度在模型中，如果单元尺寸变化很大，而且在实常数如ICONT, FTOLN, PINB, PMAX, PMIN 中应用比例系数，则可能会出现问题。
ansys电磁场分析的一些问题
1、进入求解器之后出现如下警告：Using 1 iteration per substep may result in unconvergent solutions fornonlinear analysis and the program may not indicate divergence in thiscase. Check your results.不要担心，这是ANSYS程序设置中的一个友好提示，并不说明你的程序就是错误的（当然出现这个提示也不能证明你的程序就是无误的），你可以接着运行程序。
2、进行非线性求解的时候，如果你输入了材料的磁化曲线，但是很奇怪的是在进行二维分析的时候是成功的，但是三维分析的时候被告知得到的解是不收敛。
原因一：你需要检查一下程序的问题，是不是求解的步长设置的不合理。
一般情况下，ANSYS 默认的是25，这个值可以满足大部分计算的需要，不过有的时候你需要把这个步长变大。
原因二、你输入的磁化曲线是正确的吗？是否满足ANSYS的要求？很多参考书上的磁化曲线数据并不满足ANSYS程序的要求，你需要用命令TBPLOT绘制出材料的BH、NB、MH三条曲线，看这三条曲线是不是都是光滑的单调的。
3、为什么绘制出的磁力线只有寥寥几根，如何加密呢？ANSYS中有很多中显示方式，你不妨去看看命令/SHOW中的说明.4、边界条件在电磁场分析中至关重要，那么何时采用通量平行条件何时采用通量垂直条件呢？在我看来，在ANSYS中，电磁场分析的边界条件是相当直观的。
帮助文档中是这么定义的：Flux-normal: Set the normal component of A to zeroFlux-parallel: Set in-plane components of A to zero.Far-field: Use element INFIN111.Far-field zer Use AX = AY = AZ = 0.Periodic: Use ANSYS' cyclic symmetry capability.Imposed external field: A(X,Y,Z) does not equal zero.一般情况下，通量平行条件是可以不加的，但是通量垂直条件一定要加。
ansys常见问题解答
ANSYS常见问题及应用技巧本篇开始讲述ANSYS在使用过程中常见的问题和在使用时一些常用的使用技巧，对与初学者来说，理解和弄清楚这些问题的根源和掌握这些使用技巧，能够更好的理解ANSYS这个软件本身。
1.ANSYS中的等效应力是什么物理含义？1）ANSYS中等效应力最大应力s1有什么区别，平常讨论应力分布，应该用等效应力还是最大应力s1呢？2）计算等效应力时是否需要输入等效泊松比呢？3）在实际的应用中，例如在讨论平板上的圆孔应力集中的应力分布问题时，应该用等效应力来描述应力集中的现象，还是采用主应力s1来反应集中的程度呢？还是采用一个单方向的sx来说明问题呢？答：1）这个等效应力应该就是弹塑性力学里的VonMises应力，他主要考察的是材料在各个方向上的应力差值，因为在实验室里获得材料强度都是单向载荷作用下的强度（当然现在也有三轴应力实验仪），所以有时候材料所受的单向载荷可能很大，但并没有造成破坏，这是就是看他的等效应力，具体计算公式是: σ等效=sqrt{0.5[(σ1-σ2)^2+(σ2-σ3)^2+(σ3-σ1)^2]}2）等效应力是三项主应力的组合如s，int即为max(si-sj),si,sj为三项主向应力。
i，j=1，2，3 i≠j即tresca型s,eqv为sqrt(0.5*∑(si-sj)**2)，i，j=1，2，3 i≠j即mises型3）个人认为应该采用等小应力来描述应力集中的现象，因为在实际中很难找到真正的单轴拉压的情况，一般结构的受力都没有这么简单，所以在分析的时候需要用等效应力来将各主应力进行转化，因此应该用等效应力来描述应力集中的现象。
4）等效泊松比就是泊松比，等效应力计算时不会用到泊松比，不过在计算mises 等效应变时会用到。
对于泊松比的取值原则应遵循以下两条：a:对于elastic & thermal strains 泊松比取为材料的泊松比；b:对于plastic creep hyperelastic strains 泊松比取为0.5。
ANSYS接触非线性应用技巧
Workbench 默认算法，计算速度快
Augmented Lagrange
ANSYS 经典界面默认算法，计算精度较高，收敛性好
Pure Lagrange
点对点的接触（较少使用）
Pure Lagrange on Normal; Penalty on Tangent
常用于模拟接触面之间零穿透并且存在摩擦的情况，例 如 模拟橡胶－金属之间的接触
Single contact pair
Contact pair-1
2.接触刚度不合适
缺省的接触刚度不再适用
ANSYS求解接触问题时，接触刚度对收敛性影响显著。大多数情况下缺省设置 有效，但不适于弯曲为主的接触问题。 不好收敛的信号是，收敛曲线平行于收敛准则，如下图。选择小的罚刚度值（ FKN），可以使接触容易收敛。
Error in element formulation
案例：非线性培训教材例子
项目背景
– 使用简化的2D 平面应力模型同样收敛失败
案例：非线性培训教材例子
解决方案
– 回到 DM 界面,增加倒圆角，右下角可看到小半径
案例：非线性培训教材例子
结果后处理.
案例：某机车结构件
工程背景
– 计算橡胶材料在 挤压作用下变形 和应力应变
Multipoint Constraint (MPC)
经常和Bonded联合使用，用于大模型的装配定义，可大大 降低计算时间，还或用于不同类型单元之间的连接
ANSYS接触不收敛原因及对策
遗漏接触对 接触刚度大 载荷步少 奇异 网格大变形 刚体位移 振荡 其他...
1.遗漏接触对
• 如果有些区域接触了（如自接触），但是没有定义接触对，单元会扭曲。
Ansys常见问题
Ansys常见问题ansys problems you often meet计算中途停止计算：假如觉得计算时间太长或感觉某些方面设置不对要求重新计算或停止计算，提前查看已经计算的结果（直接关闭ANSYS方法显然不可取），可以在计算的时候按ctrl+c,这样计算就停止了，然后在output 窗口中输入quit 就可以退出计算。
绘制等值线期刊上大都不用彩色，所以打出的云图一片模糊，无法识别，这时候可以选择出等值线图，但是等值线图也是彩色的，如何把它转成黑白的呢？开始是抓图后用Photoshop处理，太麻烦，ansys自己行不行呢？方法如下：1 用命令jpgprf,500,100,1将背景变为白色；2 plotctrls>device option中，把vector mode改为on,画出等值线图；3 plotctrls>style>contour>contour labeling, 将key vector mode contour labels设为on every Nth ele,对N输入一个数值，值越大，图中的label越少；4 plotctrls>style>colors>contour colors,将所有的系列都改为黑色;5 如果不喜欢ansys给出的MX,MN标志，可以用plotctrls>window controls>window options把它们去掉，将MINM 后的Mix-Min Symbols改为off就可以了。
这时候，一幅清晰的黑白等值线图就出来了。
1：ansys中的等效应力是什么物理含义？ (2)2．ansys后处理中负值的应力是压应力还是拉应力？ (3)3解决非线性分析不收敛的技巧！ (3)4非线性计算完的收敛图线，如何看他的收敛性呢，每条颜色的线代表什么意思呢？ (4)5求教accat及lccat命令 (5)6有关分块后的merge问题。
关于ansys中收敛的介绍
关于ansys中收敛的介绍1、何为收敛？在这里我引用一个会员的提问来解释这个问题：Q:结构非线性静力分析经常出现收敛这个词，如：收敛容限，收敛准则，收敛的解，位移收敛检验等，请解释，thanks!A: 个人是这样理解的：谈到收敛总会和稳定性联系在一起，简单的说，就是在进行求解过程中的一些中间值的误差对于结果的影响的大小，当中间量的误差对于你的数值积分的结果没有产生影响，就说明你的积分方法是稳定的，最终你的数值积分的结果就会收敛于精确解；当中间量的误差导致数值积分结果与精确解有很大的差别时，就说明你的方法稳定性不好，你的数值积分结果不会收敛于精确解。
我想当你对于稳定性和收敛的概念真正理解后，那些名词对于你来说，并不是问题，力学的问题最终都会和数学联系在一起，建议你看看数值积分方面的教程，学好了数学，力学对于你来说就是 a piece of cake。
Q:那么说收不收敛，最终都是因为采用的计算方法和计算参数选取的问题了？A:就本人所学的专业来说，很大程度上取决于所采用的算法，我学的是结构工程，举个例子吧 :当在进行结构动力时程分析时，采用的几种方法有线性加速度法，威尔逊－theta法，对于线性加速度法，当时间步长大于周期的0.5倍时，计算结果很可能出现不收敛，而当时间步长小于0.1倍的周期时，才有可能获得稳定的计算结果；而威尔逊－theta法，实质上就是线性加速度法的修正形式，很多实例表明当theta值大于1.37时，这种算法是无条件稳定的。
当然影响计算结果是否收敛的原因有很多，比如初始条件，我所指的仅仅是我所学专业的一个问题的很小的一个方面。
A: 说白了，就是数学。
牵涉到实际的计算问题时，才发现数学实在是太有用了，不过可惜数学实在学得不好。
A: 收敛的问题，就好像你往水里扔一块石头激起的波浪，慢慢会平息下来，这就收敛了。
计算的时候就是这样，数据在每次迭代的时候在精确解的周围震荡，最后无限趋向于精确解。
ANSYS技巧与常见问题解决
ANSYS 查询函数（Inquiry Function）在ANSYS操作过程或条件语句中，常常需要知道有关模型的许多参数值，如选择集中的单元数、节点数，最大节点号等。
此时，一般可通过*GET命令来获得这些参数。
现在，对于此类问题，我们有了一个更为方便的选择，那就是查询函数— Inquiry Function。
Inquiry Function类似于ANSYS的 *GET 命令，它访问ANSYS数据库并返回要查询的数值，方便后续使用。
ANSYS每执行一次查询函数，便查询一次数据库，并用查询值替代该查询函数。
假如你想获得当前所选择的单元数，并把它作为*DO循环的上界。
传统的方法是使用*GET命令来获得所选择的单元数并把它赋给一个变量，则此变量可以作为*DO循环的上界来确定循环的次数*get, ELMAX,elem,,count*do, I, 1, ELMAX……*enddo现在你可以使用查询函数来完成这件事，把查询函数直接放在*DO循环内，它就可以提供所选择的单元数*do, I, ELMIQR(0,13)……*enddo这里的ELMIQR并不是一个数组，而是一个查询函数，它返回的是现在所选择的单元数。
括弧内的数是用来确定查询函数的返回值的。
第一个数是用来标识你所想查询的特定实体（如单元、节点、线、面号等等），括弧内的第二个数是用来确定查询函数返回值的类型的（如选择状态、实体数量等）。
同本例一样，通常查询函数有两个变量，但也有一些查询函数只有一个变量，而有的却有三个变量。
查询函数的种类和数量很多，下面是一些常用、方便而快速快捷的查询函数1 AREA—arinqr(areaid,key)areaid—查询的面，对于key=12,13,14可取为0；key—标识关于areaidr的返回信息=1，选择状态=12，定义的数目=13，选择的数目=14，定义的最大数=-1，材料号=-2，单元类型=-3，实常数=-4，节点数=-6，单元数…arinqr(areaid,key)的返回值对于key=1=0, areaid未定义=-1，areaid未被选择=1， areaid被选择…2 KEYPOINTS—kpinqr(kpid,key)kpid—查询的关键点，对于key=12,13,14为0 key —标识关于kpid的返回信息=1，选择状态=12，定义的数目=13，选择的数目=14，定义的最大数目=-1，数料号=-2，单元类型=-3，实常数=-4，节点数，如果已分网=-7，单元数，如果已分网kpinqr(kpid,key)的返回值对于key=1=-1,未选择=0，未定义=1，选择3 LINE—lsinqr(lsid,key)lsid—查询的线段，对于key=12,13,14为0 key—标识关于lsid的返回信息=1，选择状态=2, 长度=12，定义的数目=13，选择的数目=14，定义的最大数=-1，材料号=-2，单元类型=-3，实常数=-4，节点数=-6，单元数…4 NODE—ndinqr(node,key)node—节点号，对于key=12,13,14为0 key—标识关于node的返回信息=1，选择状态=12，定义的数目=13，选择的数目=14，定义的最大数=-2，超单元标记=-3，主自由度=-4，激活的自由度=-5，附着的实体模型ndinqr(node,key)的返回值对于key=1=-1,未选择=0，未定义=1，选择5 VOLUMES—vlinqr(vnmi,key)vnmi—查询的体，对于key=12,13,14为0key—标识关于vnmi的返回信息=1，选择状态=12，定义的数目=13，选择的数目=14，定义的最大数目=-1，数料号=-2，单元类型=-3，实常数=-4，节点数=-6，单元数=-8，单元形状=-9，中节点单元=-10，单元坐标系vlinqr(vnmi,key)的返回值对于key=1=-1,未选择=0，未定义=1，选择ANSYS能实现直接流-固耦合分析吗?ANSYS流固耦合分析有三种形式，可以实现全直接或半直接耦合分析：一： ANSYS/Mechanical模块或含该模块的软件包中的流固耦合分析功能，但此处的流体是非流动的流体，而是静流体，它计算流体由于重力、惯性力、波动压力等引起的分布压力载荷与结构的相互作用。
ANSYS非线性_几何非线性分析
几何非线性分析随着位移增长，一个有限单元已移动的坐标可以以多种方式改变结构的刚度。
一般来说这类问题总是是非线性的，需要进行迭代获得一个有效的解。
大应变效应一个结构的总刚度依赖于它的组成部件（单元）的方向和单刚。
当一个单元的结点经历位移后，那个单元对总体结构刚度的贡献可以以两种方式改变变。
首先，如果这个单元的形状改变，它的单元刚度将改变。
（看图2─1(a)）。
其次，如果这个单元的取向改变，它的局部刚度转化到全局部件的变换也将改变。
（看图2─1（b)）。
小的变形和小的应变分析假定位移小到足够使所得到的刚度改变无足轻重。
这种刚度不变假定意味着使用基于最初几何形状的结构刚度的一次迭代足以计算出小变形分析中的位移。
（什么时候使用“小”变形和应变依赖于特定分析中要求的精度等级。
相反，大应变分析说明由单元的形状和取向改变导致的刚度改变。
因为刚度受位移影响，且反之亦然，所以在大应变分析中需要迭代求解来得到正确的位移。
通过发出NLGEOM，ON（GUI路径Main Menu>Solution>Analysis Options)，来激活大应变效应。
这效应改变单元的形状和取向，且还随单元转动表面载荷。
（集中载荷和惯性载荷保持它们最初的方向。
）在大多数实体单元（包括所有的大应变和超弹性单元），以及部分的壳单元中大应变特性是可用的。
在ANSYS/Linear Plus程序中大应变效应是不可用的。
图1─11 大应变和大转动大应变处理对一个单元经历的总旋度或应变没有理论限制。
（某些ANSYS单元类型将受到总应变的实际限制──参看下面。
）然而，应限制应变增量以保持精度。
因此，总载荷应当被分成几个较小的步，这可以〔NSUBST，DELTIM，AUTOTS〕，通过GUI路径Main Menu>Solution>Time/Prequent)。
无论何时当系统是非保守系统，来自动实现如在模型中有塑性或摩擦，或者有多个大位移解存在，如具有突然转换现象，使用小的载荷增量具有双重重要性。
关于ANSYS中的收敛问题共14页word资料
关于ansys中收敛准则(cnvtol)理解ansys中依据缺省的收敛准则，程序将对不平衡力SRSS与VALUE*TOLER的值进行比较；而VALUE的缺省值是在SRSS和MINREF中取较大值。
现假如TOLER的缺省值是0.1的话，这个准则是不是可以理解成后一次的SRSS是前一次的SRSS的0.1倍就收敛啦？请指点我是这样理解的例如下面的命令流:cnvtol,f,5000,0.0005,0cnvtol,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~~（无穷大）三个收敛准则。
解决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的收敛问题
ansys中收敛准则(cnvtol)理解ansys中依据缺省的收敛准则，程序将对不平衡力SRSS与VALUE*TOLER的值进行比较；而VALUE的缺省值是在SRSS和MINREF中取较大值。
现假如TOLER的缺省值是0.1的话，这个准则是不是可以理解成后一次的SRSS是前一次的SRSS的0.1倍就收敛啦？请指点我是这样理解的例如下面的命令流:cnvtol,f,5000,0.0005,0cnvtol,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~~（无穷大）三个收敛准则。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
解决非线性分析不收敛的技巧
影响非线性收敛稳定性及其速度的因素很多：
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)；
2、设置足够大的平衡迭代步数，默认为25，可以放大到很大(100)(eqit,eqit)；
3、将收敛准则调整，以位移控制时调整为0.05，以力控制为0.01(CNVTOL,lab,value,toler,norm,minref)。
4、对于线性单元和无中间节点的单元（SOLID65和SOLID45），关闭EXTRA DISPLACEMENTS OPTIONS（在OPTIONS中）。
5、对于CONCRETE材料，可以关闭压碎功能，将CONCRETE中的单轴抗压强度设置为-1(tadata,mat,shrcf-op,shrcf-cl,UntensSt,UnCompSt(-1))。
你的模型我看了一下，基本上有以下想法：
1. 模型很粗糙，做焊接的网格划那么粗，真汗
2. 绝对零度、玻尔兹曼常数未设置
3. 做焊接居然用温度直接达到1300度来模拟，太不准确
4. 瞬时热传导是直接将1300度加到单元上，当然不容易收敛，可以尝试用amplitude改进一下
5. 单元生死做的有问题，第一步中其实你已经kill掉了所有节点，却还在T-1节点上加温度。
6. 这种焊接模拟我真没有做过，其实按你的模型，再近似点，把焊缝同时升温降温也可以。
7. 我本来想改模型，不过改着改着，发现改的太慢了，心灰意冷，直接将你的前36步删除，把T-1换成weld-kill（我已经将weld-kill改成焊缝的所有节点和单元），加了个1300度，用的amp，没有出现不收敛。
1.用solid185单元
2。
CUTCONTROL,PLSLIMIT, 0.15, ! Maximum equivalent plastic strain allowed within a time-step 可以大一些
3。
substep 多一些
4。
检查约束条件，和材料属性。
我做焊接瞬态热应力分析有段时间了，你这个问题，我也遇到过，成功方法为：主要把你的子步设置的更大一些，如手工改动：NSUBST,25, , ,1，或者改更大一些NSUBST,100, , ,1，当然计算时间必成倍增长！这种方法基本上解决楼主以上的错误信息，计算正确！
另外，也可能是求解器默认的等值塑性应变比计算中的小，如设置为CUTCONTROL,PLSLIMIT,2，系统一般默认0.25左右，所以对大变形是不能适应的！！这种方法可以解决你计算中的可能的大变形问题！！
另外，再改动默认设置如下：CUTCONTROL,CRPLIMITexp,0.2,0
CUTCONTROL,CRPLIMITimp,0.25,1
CUTCONTROL,NPOINT,15
以上的三种方法完全可以使你的计算正确，这都是我在计算时花很长时间摸索出并加以解决的
实际上，这三种方法解决的就是你给的错误信息的方法！！
当然，最不建议你的划分全规则的六面体网格，保证较小的边长比！！
分享给大家，祝你成功，不能顺利计算尽管找我！！。