ansys非线性分析指南
ANSYS教程,非线性结构分析过程

ANSYS教程,非线性结构分析过程尽管非线性分析比线性分析变得更加复杂,但处理基本相同。
只是在非线形分析的适当过程中,添加了需要的非线形特性。
非线性结构分析的基本分析过程也主要由建模、加载并求解和观察结果组成。
下面来讲解其主要步骤和各个选项的处理方法。
建模这一步对线性和非线性分析都是必需的,尽管非线性分析在这一步中可能包括特殊的单元或非线性材料性质,如果模型中包含大应变效应,应力─应变数据必须依据真实应力和真实(或对数)应变表示。
加载求解在建立好有限元模型之后,将进入ANSYS求解器(GUI:Main Menu | Solution),并根据分析的问题指定新的分析类型(ANTYPE)。
求解问题的非线性特性在ANSYS中是通过指定不同的分析选项和控制选项来定义的。
非线性分析不同于线性分析之处在于,它通常要求执行多荷载步增量和平衡迭代。
下面就详细讲解一下进行非线性结构分析需要定义的各个求解选项、分析选项和控制选项是如何设置的,以及他们的意义是什么。
求解控制对于一些基本的非线性问题的分析选项,可以通过ANSYS提供的求解控制对话框中的选项设置来完成。
选择菜单路径:Main Menu | Solution | Analysis Type | Sol’n Controls,将弹出求解控制(Solution Controls)对话框,如下图所示。
从图中可以看出该对话框主要包括5个选项卡:基本选项(Basic)、瞬态选项(Transient)、求解选项(Sol’n Options)、非线性选项(Nonlinear)和高级非线性选项(Advanced NL)。
如果开始一项新的分析,在设置分析类型和非线性选项时,选择“Large Displacement Static”选项(不是所有的非线性分析都支持大变形)。
如果想要重新启动一个失败的非线性分析,则选择“Restart Current Analysis”选项。
选中下面的“Calculate prestress effects”单选按钮用于有预应力的模态分析时的预应力计算,具体内容见模态分析部分。
ANSYS结构非线性分析指南

ANSYS结构非线性分析指南ANSYS是一个强大的工程仿真软件,能够对各种复杂的结构进行分析。
其中,结构非线性分析是其中一种重要的分析方法,它能够模拟结构在非线性载荷和变形条件下的行为。
本文将为您提供一个ANSYS结构非线性分析的指南,帮助您更好地理解和应用这个方法。
首先,我们需要明确结构非线性分析的目标。
一般来说,结构非线性分析主要用于研究结构在大变形、材料非线性、接触或摩擦等复杂条件下的响应。
例如,当结构受到极大的外力作用时,其产生的变形可能会导致材料的非线性行为,这时我们就需要进行非线性分析。
在进行非线性分析之前,我们需要进行准备工作。
首先,我们需要准备一个几何模型,可以通过CAD软件导入或者直接在ANSYS中绘制。
然后,我们需要选择合适的材料模型,这将直接影响分析结果的准确性。
ANSYS提供了多种材料模型,例如线弹性模型、塑性模型和粘弹性模型等。
接下来,我们需要定义边界条件和载荷。
边界条件指明了结构的固定边界和自由边界,这决定了结构的位移约束。
载荷是作用在结构上的外力或者外界约束,例如压力、点载荷或者摩擦力等。
在非线性分析中,载荷的大小和施加方式可能会导致结构的非线性响应,因此需要仔细选择。
接下来,我们需要选择适当的非线性分析方法。
ANSYS提供了多种非线性分析方法,例如几何非线性分析、材料非线性分析和接触非线性分析等。
几何非线性分析适用于大变形情况下的分析,材料非线性分析适用于材料的弹塑性行为分析,而接触非线性分析适用于多个结构之间的接触行为分析。
在进行非线性分析之前,我们需要对模型进行预处理,包括网格划分和解算控制参数的设置。
网格划分的精度会直接影响分析结果的准确性,因此需要进行适当的剖分。
解算控制参数的设置涉及到收敛性和稳定性的问题,需要进行合理的调整。
然后,我们可以进行非线性分析了。
ANSYS提供了多种求解器,例如Newton-Raphson方法和弧长法等。
这些求解器可以通过迭代算法来求解非线性方程组,得到结构的响应结果。
ansys高级非线性分析-第九章几何不稳定性

失稳准则、米泽斯失稳准则和霍夫失稳准则等。
这些判据可以帮助我们确定结构的临界载荷和失稳模式,从而
03
采取相应的措施来提高结构的稳定性。
几何不稳定性的影响因素
材料性质
材料的弹性模量、泊松比、屈 服强度等都会影响结构的稳定
性。
结构形状和尺寸
结构的形状、尺寸、支承条件 等都会影响其稳定性。
外部载荷
外部载荷的大小、方向和分布 也会影响结构的稳定性。
案例二:高层建筑的几何不稳定性分析
总结词
高层建筑的几何不稳定性分析是确保高层建筑结构安全的重要环节。
详细描述
利用ANSYS的高级非线性分析功能,可以对高层建筑在不同风载、地震等载荷作 用下的结构响应进行模拟,评估其稳定性和安全性,为设计提供依据。
案例三:重型机械的几何不稳定性分析
总结词
重型机械的几何不稳定性分析是确保 重型机械在各种工况下安全运行的关 键。
02
几何不稳定性分析在复杂边界条件、多物理场耦合等方面的研究尚不够深入, 需要进一步拓展研究范围,完善分析方法。
03
随着计算机技术和数值计算方法的不断发展,几何不稳定性分析的计算效率和 精度将得到进一步提高,为工程实际提供更加准确和可靠的理论支持。
THANKS
感谢观看
现象。
在非线性分析中,需要考虑 结构在变形过程中形状和尺 寸的变化,以及由此引起的
力和位移的重新分布。
几何非线性行为通常出现在大 变形、应力刚化、旋转软化和
塑性流动等情况下。
几何不稳定性判据
01
几何不稳定性是指结构在某些条件下失去稳定性,发生屈曲或 失稳的现象。
02
判据是用来判断结构是否稳定的准则,常用的判据包括:欧拉
ANSYS非线性分析:1-非线性分析概述

第一章钢筋混凝土结构非线性分析概述1.1 钢筋混凝土结构的特性1.钢筋混凝土结构由两种材料组成,两者的抗拉强度差异较大,在正常使用阶段,结构或构件就处在非线性工作阶段,用弹性分析方法分析的结构内力和变形无法反映结构的真实受力状况;2.混凝土的拉、压应力-应变关系具有较强的非线性特征;3.钢筋与混凝土间的黏结关系非常复杂,特别是在反复荷载作用下,钢筋与混凝土间会产生相对滑移,用弹性理论分析的结果不能反映实际情况;4.混凝土的变形与时间有关:徐变、收缩;5.应力-应变关系莸软化段:混凝土达到强度峰值后有应力下降段;6.产生裂缝以后成为各向异形体。
混凝土结构在荷载作用下的受力-变形过程十分复杂,是一个变化的非线性过程。
11.2 混凝土结构分析的目的和主要内容《混凝土结构设计规范》(GB50010-2002)中新增的主要内容:(1)混凝土的本构关系和多轴强度:给出了单轴受压、受拉非线性应力-应变(本构)关系,混凝土二轴强度包络图、三轴抗压强度图和三轴应力状态破坏准则;(2)结构分析:规范概括了用于混凝土结构分析的5类方法,列入了结构非线性分析方法。
一、结构分析的基本目的:计算在各类荷载作用下的结构效应——内力、位移、应力、应变根据设计的结构方案确定合理的计算简图,选择不利荷载组合,计算结构内力,以便进行截面配筋计算和采取构造措施。
二、结构分析的主要内容:(1)确定结构计算简图:考虑以下因素:(a)能代表实际结构的体形和尺寸;(b)边界条件和连接方式能反映结构的实际受力状态,并有可靠的构造措施;(c)材料性能符合结构的实际情况;(d)荷载的大小、位置及组合应与结构的实际受力吻合;(e)应考虑施工偏差、初始应力及变形位移状况对计算简图进行适当修正;(f)根据结构受力特点,可对计算简图作适当简化,但应有理论或试验依据,或有可靠的工程经验;(g)结构分析结果应满足工程设计的精度要求。
(2)结构作用效应分析:根据结构施工和使用阶段的多种工况,分别进行结构分析,确定最不利荷载效应组合。
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结构非线性分析指南(一到三章)屈服准则概念:1.理想弹性材料物体发生弹性变形时,应力与应变完全成线性关系,并可假定它从弹性变形过渡到塑性变形是突然的。
2.理想塑性材料(又称全塑性材料)材料发生塑性变形时不产生硬化的材料,这种材料在进入塑性状态之后,应力不再增加,也即在中性载荷时即可连续产生塑性变形。
3.弹塑性材料在研究材料塑性变形时,需要考虑塑性变形之前的弹性变形的材料这里可分两种情况:Ⅰ.理想弹塑性材料在塑性变形时,需要考虑塑性变形之前的弹性变形,而不考虑硬化的材料,也即材料进入塑性状态后,应力不再增加可连续产生塑性变形。
Ⅱ.弹塑性硬化材料在塑性变形时,既要考虑塑性变形之前的弹性变形,又要考虑加工硬化的材料,这种材料在进入塑性状态后,如应力保持不变,则不能进一步变形。
只有在应力不断增加,也即在加载条件下才能连续产生塑性变形。
4.刚塑性材料在研究塑性变形时不考虑塑性变形之前的弹性变形。
这又可分两种情况:Ⅰ.理想刚塑性材料在研究塑性变形时,既不考虑弹性变形,又不考虑变形过程中的加工硬化的材料。
Ⅱ.刚塑性硬化材料在研究塑性变形时,不考虑塑性变形之前的弹性变形,但需要考虑变形过程中的加工硬化材料。
屈服准则的条件:1.受力物体内质点处于单向应力状态时,只要单向应力大到材料的屈服点时,则该质点开始由弹性状态进入塑性状态,即处于屈服。
2.受力物体内质点处于多向应力状态时,必须同时考虑所有的应力分量。
在一定的变形条件(变形温度、变形速度等)下,只有当各应力分量之间符合一定关系时,质点才开始进入塑性状态,这种关系称为屈服准则,也称塑性条件。
它是描述受力物体中不同应力状态下的质点进入塑性状态并使塑性变形继续进行所必须遵守的力学条件,这种力学条件一般可表示为)=Cf(σij又称为屈服函数,式中C是与材料性质有关而与应力状态无关的常数,可通过试验求得。
屈服准则是求解塑性成形问题必要的补充方程。
1.1 什么是结构非线性在日常生活中,经常会遇到结构非线性。
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材料非线性分析
弹簧单元模拟土壤
杆为平面应变分析
土壤单元选项及实常数
材料性能:d-p模型(土壤)
铁桩材料
布尔运算:叠合。
并压缩实体编号(方便划网)
选择所有
all
相加分离的线
delte 划网注意划网属性与面编号对应
定义映射的边再划分面
转换为柱坐标系
创建节点并复制
通过节点再创建单元,注意选设置好单元属性
选择r=180(圆边界节点)(from all)及外节点180*(2*1.732)/3(also select)进行全约束
转换为直角坐标
对桩进行压力加载500(Mpa)先选择节点(通过坐标值)
创建静态分析
并设置非线性收敛标准
求解保存
结果(应力、位移、剪力)。
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非线形分析指南:弹塑性
ANSYS⾮线形分析指南:弹塑性⽬录什么是塑性 (1)路径相关性 (1)率相关性 (1)⼯程应⼒、应变与真实应⼒、应变 (1)什么是激活塑性 (2)塑性理论介绍 (2)屈服准则 (2)流动准则 (3)强化准则 (3)塑性选项 (5)怎样使⽤塑性 (6)ANSYS输⼊ (7)输出量 (7)程序使⽤中的⼀些基本原则 (8)加强收敛性的⽅法 (8)查看结果 (9)塑性分析实例(GUI⽅法) (9)塑性分析实例(命令流⽅法) (14)弹塑性分析在这⼀册中,我们将详细地介绍由于塑性变性引起的⾮线性问题--弹塑性分析,我们的介绍⼈为以下⼏个⽅⾯:什么是塑性塑性理论简介ANSYS程序中所⽤的性选项怎样使⽤塑性塑性分析练习题什么是塑性塑性是⼀种在某种给定载荷下,材料产⽣永久变形的材料特性,对⼤多的⼯程材料来说,当其应⼒低于⽐例极限时,应⼒⼀应变关系是线性的。
另外,⼤多数材料在其应⼒低于屈服点时,表现为弹性⾏为,也就是说,当移⾛载荷时,其应变也完全消失。
由于屈服点和⽐例极限相差很⼩,因此在ANSYS 程序中,假定它们相同。
在应⼒⼀应变的曲线中,低于屈服点的叫作弹性部分,超过屈服点的叫作塑性部分,也叫作应变强化部分。
塑性分析中考虑了塑性区域的材料特性。
路径相关性:即然塑性是不可恢复的,那么这种问题的就与加载历史有关,这类⾮线性问题叫作与路径相关的或⾮保守的⾮线性。
路径相关性是指对⼀种给定的边界条件,可能有多个正确的解—内部的应⼒,应变分布—存在,为了得到真正正确的结果,我们必须按照系统真正经历的加载过程加载。
率相关性:塑性应变的⼤⼩可能是加载速度快慢的函数,如果塑性应变的⼤⼩与时间有关,这种塑性叫作率⽆关性塑性,相反,与应变率有关的性叫作率相关的塑性。
⼤多的材料都有某种程度上的率相关性,但在⼤多数静⼒分析所经历的应变率范围,两者的应⼒-应变曲线差别不⼤,所以在⼀般的分析中,我们变为是与率⽆关的。
⼯程应⼒,应变与真实的应⼒、应变:塑性材料的数据⼀般以拉伸的应⼒—应变曲线形式给出。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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 结果评价和分析1.3 非线性结构静力分析非线性结构的定义:在日常生活中,经常会遇到结构非线性,例如:无论何时用钉书针钉书,金属钉书钉将永久地弯曲成一个不同的形状,如图1─1a;如果你在一个木架上放置重物,随着时间的迁移它将越来越下垂,如图1─1b;当在汽车或卡车上装货时它的轮胎和下面路面间接触将随货物的重量而变化,如图1─1c;如果将上面例子的载荷- 变形曲线画出来,你将发现它们都显示了非线性结构的基本特征- 变化的结构刚性,即结构变形与载荷间不是线性关系。
非线性行为的原因引起结构非线性的原因很多,它可以被分成三种主要类型:1. 状态变化,包括接触许多普通的结构表现出一种与状态相关的非线性行为,例如一根只能拉伸的电缆可能是松散下垂的,也可能是绷紧的。
轴承套可能是接触的,也可能是不接触的;冻土可能是冻结的,也可能是融化的。
这些系统的刚度由于系统状态的改变在不同的值之间突然变化,状态改变也许和载荷直接有关,如在电缆情况中;也可能由某种外部原因引起,如在冻土中的紊乱热力学条件。
ANSYS 程序中单元的激活与杀死选项用来给这种状态的变化建模。
接触是一种很普遍的非线性行为,接触是状态变化非线性类型形中一个特殊而重要的子集。
2. 几何非线性如果结构经受大变形,它变化的几何形状可能会引起结构的非线性响应。
一个例子是钓鱼竿的垂向刚性,随着垂向载荷的增加,钓鱼杆不断弯曲,以致于动力臂明显地减少,导致杆端显示出在较高载荷下不断增长的刚性。
3. 材料非线性非线性的应力- 应变关系是结构非线性的常见原因,许多因素可以影响材料的应力- 应变性质,包括加载历史(如在弹- 塑性响应状况下);环境状况(如温度),加载的时间总量(如在蠕变响应状况下)。
4. 牛顿一拉森方法ANSYS 程序的方程求解器,计算一系列的联立线性方程来预测工程系统的响应。
然而非线性结构的行为不能直接用这样一系列的线性方程表示,需要一系列的带校正的线性近似来求解非线性问题。
5. 逐步递增载荷和平衡迭代一种近似的非线性求解是将载荷分成一系列的载荷增量,可以在几个载荷步内,或者在一个载荷步的几个子步内施加载荷增量。
在每一个增量的求解完成后,继续进行下一个载荷增量之前,程序调整刚度矩阵以反映结构刚度的非线性变化。
遗憾的是纯粹的增量近似不可避免地随着每一个载荷增量积累误差,可能导致结果最终失去平衡,如图1─3a所示:ANSYS 程序通过使用牛顿拉普森平衡迭代克服了这种困难。
它迫使在每一个载荷增量的末端,解达到平衡收敛在某个容限范围内。
图1─3b描述了在单自由度非线性分析中牛顿拉普森平衡迭代的使用,在每次求解前NR 方法估算出残差矢量,这个矢量是回复力对应于单元应力的载荷和所加载荷的差值,然后程序使用非平衡载荷进行线性求解,且核查收敛性,如果不满足收敛准则重新估算非平衡载荷,修改刚度矩阵获得新解;持续这种迭代过程,直到问题收敛。
ANSYS 程序提供了一系列命令来增强问题的收敛性,如自适应下降、线性搜索、自动载荷步及二分等,可被激活来加强问题的收敛性。
如果不能得到收敛,那么程序或者继续计算下一个载荷步或者终止(依据你的指示)。
5. 弧长法对某些物理意义上不稳定系统的非线性静态分析,如果你仅仅使用NR 方法,正切刚度矩阵可能变为降秩短阵,导致严重的收敛问题。
这样的情况包括独立实体从固定表面分离的静态接触分析,结构或者完全崩溃或者突然变成另一个稳定形状的非线性弯曲问题。
对这样的情况你可以激活另外一种迭代方法- 弧长方法来帮助稳定求解。
弧长方法导致NR 平衡迭代沿一段弧收敛,从而即使当正切刚度矩阵的倾斜为零或负值时也往往阻止发散。
这种迭代方法以图形表示在图1─4中:6. 非线性求解的组织级别非线性求解被分成三个操作级别:载荷步、子步、平衡迭代。
- 顶层级别由在一定时间范围内你明确定义的载荷步组成。
假定载荷在载荷步内是线性地变化的。
- 在每一个载荷步内,为了逐步加载,可以控制程序来执行多次求解(子步或时间步)。
- 在每一个子步内,程序将进行一系列的平衡迭代以获得收敛的解。
图1─5说明了一段用于非线性分析的典型的载荷历史。
7. 收敛容限当你对平衡迭代确定收敛容限时,你必须回答这些问题:- 你想基于载荷变形还是联立二者来确定收敛容限?- 既然径向偏移(以弧度度量) 比对应的平移小,你是不是想对这些不同的条目建立不同的收敛准则?当你确定收敛准则时,ANSYS 程序会给你一系列的选择:你可以将收敛检查建立在力、力矩、位移、转动或这些项目的任意组合上,另外每一个项目可以有不同的收敛容限值,对多自由度问题你同样也有收敛准则的选择问题。
当你确定你的收敛准则时,记住以力为基础的收敛提供了收敛的绝对量度,而以位移为基础的收敛仅提供了表观收敛的相对量度。
因此你应当如果需要总是使用以力为基础(或以力矩为基础) 的收敛容限,如果需要可以增加以位移为基础(或以转动为基础) 的收敛检查,但是通常不单独使用位移或转动为基础。
图1─6说明了一种单独使用位移收敛检查导致出错的情况。
在第二次迭代后计算出的位移很小可能被认为是收敛的解,尽管问题仍旧远离真正的解,要防止这样的错误应当使用力收敛检查。
8. 保守行为与非保守行为过程依赖性如果通过外载输入系统的总能量当载荷移去时复原,我们说这个系统是保守的;如果能量被系统消耗(如由于塑性应变或滑动摩擦),我们说系统是非保守的。
一个非守恒系统的例子显示在图1─7。
一个保守系统的分析是与过程无关的,通常可以任何顺序和以任何数目的增量加载而不影响最终结果;相反地,一个非保守系统的分析是过程相关的,必须紧紧跟随系统的实际加载历史以获得精确的结果。
如果对于给定的载荷范围可以有多于一个的解是有效的(如在突然转变分析中)。
这样的分析也可能是过程相关的。
过程相关问题通常要求缓慢加载(也就是使用许多子步),逐步达到最终的载荷值。
9. 子步当使用多个子步时,你需要考虑精度和代价之间的平衡。
更多的子步(也就是小的时间步) 通常导致较好的精度,但以增多的运行时间为代价。
ANSYS 提供两种方法来控制子步数:子步数或时间步长我们既可以通过指定实际的子步数也可以通过指定时间步长控制子步数。
自动时间步长ANSYS 程序基于结构的特性和系统的响应来调整时间步长。
子步数如果你的结构在它的整个加载历史期间显示出高度的非线性特点,而且你对结构的行为理解足够好,可以确保得到收敛的解,那么你也许能够自己确定多小的时间步长是必需的,且对所有的载荷步使用这同一时间步务必允许足够大的平衡迭代数。
自动时间分步如果你预料你的结构的行为将从线性到非线性变化,你也许想要在系统响应的非线性部分采用变化的时间步长。
在这样一种情况,你可以激活自动时间分步,以便根据需要调整时间步长,获得精度和代价之间的良好平衡。
同样地,如果你不确信你的问题将成功地收敛,你也许想要使用自动时间分步来激活ANSYS 程序的二分特点。
二分法提供了一种对收敛失败自动矫正的方法。
无论何时,只要平衡迭代收敛失败,二分法将把时间步长分成两半,然后从最后收敛的子步自动重启动,如果已二分的时间步再次收敛失败,二分法将再次分割时间步长,然后重启动。
持续这一过程直到获得收敛或到达最小时间步长(由你指定)。
载荷和位移方向当结构经历大变形时,应该考虑到载荷将发生了什么变化。
在许多情况中,无论结构如何变形,施加在系统中的载荷都保持恒定的方向;而在另一些情况中,力将改变方向随着单元方向的改变而变化。
ANSYS 程序对这两种情况都可以建模,依赖于所施加的载荷类型。
加速度和集中力将不管单元方向的改变而保持它们最初的方向,表面载荷作用在变形单元表面的法向且可被用来模拟跟随力。
图1─8说明了恒力和跟随力。
注意- 在大变形分析中不修正结点坐标系方向,因此计算出的位移在最初的方向上输出。
非线性瞬态过程的分析用于分析非线性瞬态行为的过程与对线性静态行为的处理相似,以步进增量加载程序在每一步中进行平衡迭代。