ANSYS几何非线性概述
- 格式:pdf
- 大小:315.10 KB
- 文档页数:5
ANSYS基础教程,非线性分析
由荷载-变形曲线将会发现非线性结构的基本特征:变化的结构刚度。
引起非线性的原因
引起非线性行为的原因很多,这里介绍三种主要原因:
几何非线性
如果结构经受大变形,它变化的几何形状可能会引起结构的非线性响应,例如:随着钓鱼竿钓到鱼,竖向荷载就增加,杆不断弯曲以至于动力臂明显减少,导致杆端显示出较高的荷载下不断增长的刚性。
材料非线性
非线性的应力-应变关系是造成结构的非线性的常见的原因。
许多因素可以影响材料的应力应变性质,包括加载历史(如在弹塑性响应状况下)、环境状况(如温度)、加载的时间总量(如在蠕变响应情况下)。
状态非线性
许多普通结构表现出一种与状态相关的非线性行为,例如,一根只能拉伸的电缆可能是松散的,也可能是绷紧的;轴承套可能是接触的,也可能是不接触的;冻土可能是冻结的,也可能是融化的。
这些系统的刚度由于系统状态的改变在不同的值之间变化。
状态改变也许
和荷载直接有关(如在电缆情况下),也可能由某种外部原因引起(如冻土中的紊乱力学条件)。
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的超弹拟合功能做一简单介绍。
第一章钢筋混凝土结构非线性分析概述1.1 钢筋混凝土结构的特性1.钢筋混凝土结构由两种材料组成,两者的抗拉强度差异较大,在正常使用阶段,结构或构件就处在非线性工作阶段,用弹性分析方法分析的结构内力和变形无法反映结构的真实受力状况;2.混凝土的拉、压应力-应变关系具有较强的非线性特征;3.钢筋与混凝土间的黏结关系非常复杂,特别是在反复荷载作用下,钢筋与混凝土间会产生相对滑移,用弹性理论分析的结果不能反映实际情况;4.混凝土的变形与时间有关:徐变、收缩;5.应力-应变关系莸软化段:混凝土达到强度峰值后有应力下降段;6.产生裂缝以后成为各向异形体。
混凝土结构在荷载作用下的受力-变形过程十分复杂,是一个变化的非线性过程。
11.2 混凝土结构分析的目的和主要内容《混凝土结构设计规范》(GB50010-2002)中新增的主要内容:(1)混凝土的本构关系和多轴强度:给出了单轴受压、受拉非线性应力-应变(本构)关系,混凝土二轴强度包络图、三轴抗压强度图和三轴应力状态破坏准则;(2)结构分析:规范概括了用于混凝土结构分析的5类方法,列入了结构非线性分析方法。
一、结构分析的基本目的:计算在各类荷载作用下的结构效应——内力、位移、应力、应变根据设计的结构方案确定合理的计算简图,选择不利荷载组合,计算结构内力,以便进行截面配筋计算和采取构造措施。
二、结构分析的主要内容:(1)确定结构计算简图:考虑以下因素:(a)能代表实际结构的体形和尺寸;(b)边界条件和连接方式能反映结构的实际受力状态,并有可靠的构造措施;(c)材料性能符合结构的实际情况;(d)荷载的大小、位置及组合应与结构的实际受力吻合;(e)应考虑施工偏差、初始应力及变形位移状况对计算简图进行适当修正;(f)根据结构受力特点,可对计算简图作适当简化,但应有理论或试验依据,或有可靠的工程经验;(g)结构分析结果应满足工程设计的精度要求。
(2)结构作用效应分析:根据结构施工和使用阶段的多种工况,分别进行结构分析,确定最不利荷载效应组合。
第一章钢筋混凝土结构非线性分析概述1.1 钢筋混凝土结构的特性1.钢筋混凝土结构由两种材料组成,两者的抗拉强度差异较大,在正常使用阶段,结构或构件就处在非线性工作阶段,用弹性分析方法分析的结构内力和变形无法反映结构的真实受力状况;2.混凝土的拉、压应力-应变关系具有较强的非线性特征;3.钢筋与混凝土间的黏结关系非常复杂,特别是在反复荷载作用下,钢筋与混凝土间会产生相对滑移,用弹性理论分析的结果不能反映实际情况;4.混凝土的变形与时间有关:徐变、收缩;5.应力-应变关系莸软化段:混凝土达到强度峰值后有应力下降段;6.产生裂缝以后成为各向异形体。
混凝土结构在荷载作用下的受力-变形过程十分复杂,是一个变化的非线性过程。
1.2 混凝土结构分析的目的和主要内容《混凝土结构设计规范》(GB50010-2002)中新增的主要内容:(1)混凝土的本构关系和多轴强度:给出了单轴受压、受拉非线性应力-应变(本构)关系,混凝土二轴强度包络图、三轴抗压强度图和三轴应力状态破坏准则;(2)结构分析:规范概括了用于混凝土结构分析的5类方法,列入了结构非线性分析方法。
一、结构分析的基本目的:计算在各类荷载作用下的结构效应——内力、位移、应力、应变根据设计的结构方案确定合理的计算简图,选择不利荷载组合,计算结构内力,以便进行截面配筋计算和采取构造措施。
二、结构分析的主要内容:(1)确定结构计算简图:考虑以下因素:(a)能代表实际结构的体形和尺寸;(b)边界条件和连接方式能反映结构的实际受力状态,并有可靠的构造措施;(c)材料性能符合结构的实际情况;(d)荷载的大小、位置及组合应与结构的实际受力吻合;(e)应考虑施工偏差、初始应力及变形位移状况对计算简图进行适当修正;(f)根据结构受力特点,可对计算简图作适当简化,但应有理论或试验依据,或有可靠的工程经验;(g)结构分析结果应满足工程设计的精度要求。
(2)结构作用效应分析:根据结构施工和使用阶段的多种工况,分别进行结构分析,确定最不利荷载效应组合。
第二章材料本构关系§2.1本构关系的概念本构关系:应力与应变关系或内力与变形关系结构的力学分析,必须满足三类基本方程:(1)力学平衡方程:结构的整体或局部、静力荷载或动力荷载作用下的分析、精确分析或近似分析都必须满足;(2)变形协调方程:根据结构的变形特点、边界条件和计算精度等,可精确地或近似地满足;(3)本构关系:是连接平衡方程和变形协调方程的纽带,具体表达形式有:材料的应力-应变关系,截面的弯矩-曲率关系,轴力-变形(伸长、缩短)关系,扭矩-转角关系,等等。
所有结构(不同材料、不同结构形式和体系)的力学平衡方程和变形协调方程原则上相同、数学形式相近,但本构关系差别很大。
有弹性、弹塑性、与时间相关的粘弹性、粘塑性,与温度相关的热弹性、热塑性,考虑材料损伤的本构关系,考虑环境对材料耐久性影响的本构关系,等等。
正确、合理的本构关系是可靠的分析结果的必要条件。
混凝土结构非线性分析的复杂性在于:钢筋混凝土---复杂的本构关系:有限元法---结构非线性分析的工具:非线性全过程分析---解决目前结构分析与结构设计理论矛盾的途径:§2.2 一般材料本构关系分类1.线弹性(a) 线性本构关系; (b) 非线性弹性本构关系图2-1 线弹性与非线性弹性本构关系比较在加载、卸载中,应力与应变呈线性关系:}]{[}{εσD = (图2-1a ) 适用于混凝土开裂前的应力-应变关系。
2. 非线性弹性在加载、卸载中,应力与应变呈非线性弹性关系。
即应力与应变有一一对应关系,卸载沿加载路径返回,没有残余变形(图2-1b )。
}{)]([}{εεσD = 或 }{)]([}{εσσD =适用于单调加载情况结构力学性能的模拟分析。
3. 弹塑性图2 – 2 弹塑性本构关系(a)典型弹塑性;(b)理想弹塑性;(c)线性强化;(d)刚塑性典型的钢筋拉伸应力、应变曲线 (图2-2(a ))包含弹性阶段(OA )、流动阶段(AB )及硬化阶段(BC )。
ANSYS结构非线性分析指南ANSYS是一款非常强大的有限元分析软件,广泛应用于各种工程领域的结构分析。
在常规的结构分析中,通常会涉及到线性分析,但一些情况下,结构出现了非线性行为,这时就需要进行非线性分析。
非线性分析可以更准确地模拟结构的真实行为,包括材料的非线性、几何的非线性和接触非线性等。
在进行ANSYS结构非线性分析时,需要考虑以下几个方面:1.材料的非线性:在材料的应力-应变关系中,材料的性质可能会发生变化,如塑性变形、损伤、软化等。
因此在非线性分析中,需要考虑材料的非线性特性,并正确选取材料模型。
2.几何的非线性:在一些情况下,结构本身的几何形态可能会发生较大变化,如大变形、屈曲等。
这需要考虑结构的几何非线性,并在分析中充分考虑结构的形变情况。
3.接触非线性:当结构中存在接触面时,接触面之间的接触力可能是非线性的,如摩擦力、法向压力等。
在进行非线性分析时,需要考虑接触面上的非线性行为,确保接触的可靠性。
在进行ANSYS结构非线性分析时,可以按照以下步骤进行:1.建立模型:首先需要根据实际情况建立结构的有限元模型,包括几何形状、边界条件和加载条件等。
在建立模型时,需要考虑到结构的材料、几何和接触情况,并进行合理的网格划分。
2.设置分析类型:在ANSYS中,可以选择静力分析、动力分析等不同的分析类型。
在进行非线性分析时,需要选择适合的非线性分析模块,并设置相应的参数。
3.设置材料模型:根据结构的材料特性,选择合适的材料模型,如弹塑性模型、本构模型等。
根据实际情况,设置材料的材料参数,确保材料的非线性行为能够得到准确的描述。
4.设置几何非线性:考虑结构的几何非线性时,需要选择合适的几何非线性选项,并设置合适的几何参数。
在进行大变形分析时,需要选择几何非线性选项,确保结构的形变情况能够得到准确的描述。
5.设置接触非线性:当结构存在接触面时,需要考虑接触面上的非线性行为。
在ANSYS中,可以设置接触类型、摩擦系数等参数,确保接触的可靠性。
ansys 非线性分析原理ANSYS中的非线性分析是指通过考虑材料的非线性行为、几何非线性和边界条件的非线性等因素,对结构进行分析和计算。
非线性分析的原理主要包括以下几个方面。
1. 材料的非线性行为:考虑到材料在受载作用下的非线性行为,一般采用弹塑性分析方法。
弹塑性材料在受力时会出现应力-应变曲线的非线性特征,这需要使用合适的本构模型来描述。
ANSYS中常用的本构模型有弹塑性模型、弹性模型等,根据问题的实际情况选择适当的本构模型进行分析。
2. 几何的非线性效应:当结构在受载作用下出现较大的变形时,就需要考虑几何非线性效应。
一般情况下,当结构的变形较小时可以忽略几何非线性,反之则需要进行几何非线性分析。
几何非线性的分析可通过使用大变形理论来描述结构的非线性变形,并进行相应的计算。
3. 边界条件的非线性效应:非线性分析还需要考虑边界条件的非线性效应。
在实际工程中,边界条件往往是随着结构的变形而变化的,如约束条件的变化、边界载荷的变化等。
这些非线性边界条件会对结构的响应产生影响,因此需要将其考虑在内进行非线性分析。
在ANSYS中进行非线性分析时,通常需要进行以下步骤:1. 定义材料的本构模型:选择合适的弹塑性模型或弹性模型,并设置相应的参数。
2. 构建几何模型:根据实际工程要求,构建结构的几何模型,并对其进行离散化,即将结构分割成有限元网格。
3. 施加边界条件和载荷:根据实际工况,为结构施加边界条件和载荷。
4. 求解非线性方程组:通过非线性方程的迭代求解方法,求解得到结构的非线性响应。
5. 分析结果的后处理:对求解得到的结果进行分析和后处理,获取所需的工程参数和信息。
总之,非线性分析在ANSYS中是通过考虑材料的非线性行为、几何的非线性效应和边界条件的非线性效应等因素,对结构进行全面分析和计算的方法。
引用小健哥的ANSYS 非线性分析命令解析ANSYS应用基于问题物理特性的自动求解控制方法,把各种非线性分析控制参数设置到合适的值。
如果用户对这些设置不满意,还可以手工设置。
下列命令的缺省设置已进行了优化处理:AUTOTS PRED MONITORDELTIM NROPT NEQITNSUBST TINTP SSTIFCNVTOL CUTCONTROL KBCLNSRCH OPNCONTROL EQSLVARCLEN CDWRITE LSWRITE这些命令及其设置在将在后面讨论。
参见《ANSYS Commands Reference》。
如果用户选择自己的设置而不是ANSYS的缺省设置,或希望用以前版本的ANSYS的输入列表,则可用/ SOLU 模块的SOLCONTROL ,OFF命令,或在/ BATCH 命令后用/ CONFIG ,NLCONTROL,OFF命令。
参见SOLCONTROL 命令的详细描述。
ANSYS对下面的分析激活自动求解控制单场的非线性或瞬态结构以及固体力学分析,在求解自由度为UX、UY、UZ、ROTX、ROTY、ROTZ 的结合时;单场的非线性或瞬态热分析,在求解自由度为TEMP时;注意--本章后面讨论的求解控制对话框,不能对热分析做设置。
用户必须应用标准的ANSYS 求解命令或GUI来设置。
2.2 非线性静态分析步骤尽管非线性分析比线性分析变得更加复杂,但处理基本相同。
只是在非线形分析的过程中,添加了需要的非线形特性。
非线性静态分析是静态分析的一种特殊形式。
如同任何静态分析,处理流程主要由以下主要步骤组成:建模;设置求解控制;设置附加求解控制;加载;求解;考察结果。
2.2.1 建模这一步对线性和非线性分析基本上是一样的,尽管非线性分析在这一步中可能包括特殊的单元或非线性材料性质,参考§4《材料非线性分析》,和§6.1《单元非线性》。
如果模型中包含大应变效应,应力─应变数据必须依据真实应力和真实(或对数)应变表示。
ANSYS 几何非线性Ansys中有多种变形定义:1。
LARGE STRAIN:从ANSYS理论说明书的推导来看,引入了对数应变,在原贴中认为这种对数应变是真应变,笔者觉得这一点值得探讨,我们知道应变定义一般分为三种,即工程应变,或称为柯西应变,即小变形情况下通常的应变定义。
大变形情况下,以初始构形为基础,可定义格林应变,以现实构形为基础,可定义阿耳曼西应变。
笔者认为,对数应变实际上是为简化大变形分析,定义的一种应变形态。
如果在ansys中打开大变形效应,那么对数应变的引入是不必要的,因为无论格林应变还是阿耳曼西应变都考虑了现实构形和初始构形的区别。
在小变形分析中,也无必要引入对数应变,因为此时初始构形和现实构形差别甚微,不必进行区分。
但如果实际变形较大(超过10%),而在分析中又未打开大变形效应,则此时使用对数应变是有必要的。
(因为柯西应变为非可加应变,对数应变是可加的)。
所以说对数应变是在未打开大变形效应开关的前提下为考虑较大应变而人为设计的一种应变形态。
2。
LARGE ROTATION 大转动问题。
著名的例子就是钓鱼竿问题。
钓鱼竿是大变形,小应变问题。
钓鱼竿在变形时,能引起应变的变形很小,横向大变形是因为刚体转动而引起的。
但是,在连续体力学中,对刚体转动的处理是非常关键的。
为什么说很关键呢,因为刚体运动是不产生应变的,但对于某些应变定义(Strain Measures),比如上述的工程应变或者对数应变,直接代入计算后应变不为零(通过位移求得应变),也就会计算出应力来。
因此,需要消除元刚体运动的影响。
以上一段是sjtu79编辑所写,笔者要补充的是在这种情况下最好的办法就是引入格林应变或阿耳曼西应变,它们均是可以消除刚体运动影响的(通过在本构方程中使用格林或阿耳曼西本构速率),与上述两种应变对应的应力是克希霍夫应力和欧拉应力。
之所以在有限元分析中采用格林应变,是因为我们所知的初始条件均是针对初始构形的。
ANSYS求解非线性问题牛顿一拉森方法ANSYS程序的方程求解器计算一系列的联立线性方程来预测工程系统的响应。
然而,非线性结构的行为不能直接用这样一系列的线性方程表示。
需要一系列的带校正的线性近似来求解非线性问题。
逐步递增载荷和平衡迭代一种近似的非线性救求解是将载荷分成一系列的载荷增量。
可以在几个载荷步内或者在一个载步的几个子步内施加载荷增量。
在每一个增量的求解完成后,继续进行下一个载荷增量之前程序调整刚度矩阵以反映结构刚度的非线性变化。
遗憾的是,纯粹的增量近似不可避免地随着每一个载荷增量积累误差,导种结果最终失去平衡,如图1所示所示。
(a)纯粹增量式解(b)全牛顿-拉普森迭代求解图1 纯粹增量近似与牛顿-拉普森近似的关系ANSYS程序通过使用牛顿-拉普森平衡迭代克服了这种困难,它迫使在每一个载荷增量的末端解达到平衡收敛(在某个容限范围内)。
图1(b)描述了在单自由度非线性分析中牛顿-拉普森平衡迭代的使用。
在每次求解前,NR方法估算出残差矢量,这个矢量是回复力(对应于单元应力的载荷)和所加载荷的差值。
程序然后使用非平衡载荷进行线性求解,且核查收敛性。
如果不满足收敛准则,重新估算非平衡载荷,修改刚度矩阵,获得新解。
持续这种迭代过程直到问题收敛。
ANSYS程序提供了一系列命令来增强问题的收敛性,如自适应下降,线性搜索,自动载荷步,及二分等,可被激活来加强问题的收敛性,如果不能得到收敛,那么程序或者继续计算下一个载荷前或者终止(依据你的指示)。
对某些物理意义上不稳定系统的非线性静态分析,如果你仅仅使用NR方法,正切刚度矩阵可能变为降秩短阵,导致严重的收敛问题。
这样的情况包括独立实体从固定表面分离的静态接触分析,结构或者完全崩溃或者“突然变成”另一个稳定形状的非线性弯曲问题。
对这样的情况,你可以激活另外一种迭代方法,弧长方法,来帮助稳定求解。
弧长方法导致NR平衡迭代沿一段弧收敛,从而即使当正切刚度矩阵的倾斜为零或负值时,也往往阻止发散。
ANSYS几何非线性概述
一、什么是非线性
什么是非线性(non-linear)?按照百度百科的解释,非线性是指变量之间的数学关系不是直线而是曲线、曲面或不确定的属性。
而对于工程结构而言,非线性或者说非线性行为,是指外部荷载引起工程结构刚度显著改变的一种行为。
如果绘制一个非线性结构的荷载-位移曲线,则力与位移的曲线为非线性函数。
ANSYS非线性主要分为以下三大类:
1、几何非线性
大应变、大位移、大旋转
2、材料非线性
塑性、超弹性、粘弹性、蠕变
3、状态改变非线性
接触、单元生死
其中几何非线性和材料非线性是土木工程结构计算中最为常见的两种类型。
二、结构几何非线性概念理解
如果一个结构在受荷的过程经历了大变形,则变化后的几何形状能引起非线性行为。
例如,上述例子,
杆梢在轻微横向作用下是柔软的,
当外部横向荷载加大时,杆的几何形状发生改变
,力矩臂减小,引起杆的刚化响应。
几何非线性主要分为如下三种现象:
1.单元的形状改变(面积、厚度),其单独的单元刚度也将改变
2.单元的取向发生转动,其局部刚度在转化为全局分量时将会发生变化。
3.单元应变产生较大的平面内应力状态引起平面法向刚度的改变。
随着垂直挠度UY 的增加,较大的膜应力SX 将会导致刚化效应。
上述三种情况的关系如下:
应力刚化
三、ANSYS
几何非线性注意事项
1、建模注意事项 (a )单元选择注意事项
在定义单元类型时,应明白如果分析的过程中有几何非线性,应确保所选单元类型支持相应的几何非线性效应。
例如shell63单元支持应力刚化和大挠度,但不支持大应变;而shell181则支持所有的三类几何非线性,可在单元描述的特殊特征列表中找到类似信息。
特别是在选择接触单元的时候应慎重,有的接触单元是没有任何非线性能力,例如CONTAC52.
同时应注意剪切锁定以及体积锁定等不可压缩性所带来的收敛困难。
(b )预见网格扭曲
ANSYS 在第一迭代之前,会检查网格的质量;在大应变分析中,迭代计算过后的网格或许会变得严重扭曲,为防止出现不良形状,可以预见网格扭曲从而修改原始网格。
(c )足够的网格密度
为防止网格离散化错误,必须有足够的网格密度,否则就很容易造成等值线图不连续,同时如果要捕捉弯曲响应,壳和梁单元的网格密度应足够多,计算中不应有角度超过30度的单元。
一分为二,作为 三角形,形状保 持较好。
(d)耦合和约束方程
ANSYS在几何非线性计算过程中,节点坐标系不会因为考虑大转动而修正,耦合和约束方程总是作用在原始方向。
因此在大位移分析中一般要避免耦合和约束方程,连接转动和位移自由度的约束方程是基于线性和小挠度理论,因此在使用耦合和约束方程之前应仔细考虑。
2.求解注意事项
(a)何时应选择大变形?ANSYS如何打开大变形?
大变形效应可以改善求解精度,但相应地,其求解代价也加倍增加。
如果在分析问题的过程中,100%确定大变形效应不重要,这时可只选择小变形分析以使得求解效率最快,但如何对求解的问题有任何疑问,则可以始终采用大变形。
ANSYS中采用Nlgeom命令打开与关闭大变形效应。
(b)加载与边界条件
当考虑结构大变形效应时,载荷在很多情况下方向将保持不变,但是针对于某些特别载荷,当单元经历大转动时,载荷方向跟随单元而改变。
ANSYS可以根据所施加的载荷类型来模拟不同情况。
注意,在大应变分析中,压力施加于更新的面,因此由压力产生的总载将随免得伸长或缩短而变化。
在大变形分析中,一般应避免单点约和单点力,不然很容易造成应力奇异的现象。
(c)求解
求解过程中时间步长应足够小,使得在任何一个子步内,没有超过10度转动的单元。
如果在反复二分之后,模型还不能在全载荷处收敛,原因可能是实际物理上的不稳定(屈曲或者全塑性截面),这时候可以画出载荷—挠度曲线,看切向刚度是否为0。