几何非线性分析
- 格式:doc
- 大小:96.00 KB
- 文档页数:18
线性分析在结构方面就是指应力应变曲线刚开始的弹性部分,也就是没有达到应力屈服点的结构分析非线性分析包括状态非线性,几何非线性,以及材料非线性,状态非线性比如就是钓鱼竿,几何比如就是物体的大变形,材料比如就是塑性材料属性。
2.非线性行为的原因引起结构非线性的原因很多,主要可分为以下3种类型。
(1)状态变化(包括接触)许多普通结构表现出一种与状态相关的非线性行为。
例如,一根只能拉伸的电缆可能是松弛的,也可能是绷紧的;轴承套可能是接触的,也可能是不接触的;冻土可能是冻结的,也可能是融化的。
这些系统的刚度由于系统状态的改变而突然变化。
状态改变或许和载荷直接有关(如在电缆情况中),也可能是由某种外部原因引起的(如在冻土中的紊乱热力学条件)。
接触是一种很普遍的非线性行为,接触是状态变化非线性类型中一个特殊而重要的子集。
(2)几何非线性(3)材料非线性非线性的应力-应变关系是结构非线性的常见原因。
许多因素可以影响材料的应力-应变性质,包括加载历史(如在弹-塑性响应状况下)、环境状况(如温度)、加载的时间总量(如在蠕变响应状况下)等。
3.非线性结构分析中应注意的问题(1)牛顿-拉普森方法ANSYS程序的方程求解器可以通过计算一系列的联立线性方程来预测工程系统的响应。
然而,非线性结构的行为不能直接用这样一系列的线性方程来表示,需要一系列的带校正的线性近似来求解非线性问题。
几何非线性分析随着位移增长,一个有限单元已移动的坐标可以以多种方式改变结构的刚度。
一般来说这类问题总是非线性的,需要进行迭代获得一个有效的解。
大应变效应大应变处理对一个单元经历的总应变没有理论限制。
然而,应限制应变增量以保持精度。
因此,总载荷应当被分成几个较小的载荷步。
无论何时(当系统是非保守系统,在模型中有塑性或摩擦,或者有多个大位移解存在时)使用小的载荷增量都具有双重重要性。
材料非线性分析1.概述(1)塑性定义塑性是指在某种给定载荷下,材料产生永久变形的特性,对于大多的工程材料来说,当其内部应力低于比例极限时,应力与应变的关系是线性的。
几何非线性
1)几何非线性一般是由于变形后的状态与变形之前相差较大,必须
在变形后的状态下重新建立物体的平衡方程,来保证问题的准确性,主要
包括——大挠度、大转动、大曲率三种:
2)几何非线性的计算方法有两种——TL方法(完全拉格朗日方法)、UL方法(更新拉格朗日方法),其中UL方法的平衡是建立在变形后坐标
系上,具有更强的适用性,如下图:
3)和材料非线性问题相比,几何非线性问题有更为复杂多样的荷载-
位移的平衡路径,存在荷载控制下的位移跳跃和位移控制下的荷载跌落,
如下图:
4)对于这种具有负斜率的荷载-位移曲线,能够解决的求解算法是——弧长法,防止了一般基于增量法的牛顿-辛普森方法的双折返现象,其
越过顶峰的示意图如下:
5)接下来的文章,将通过例如轴对称固支圆板的几何非线性分析,
就该问题做更详细的说明:。
几何非线性分析中应注意的几个问题
1.因为非线性,查看某个单独荷载的作用结果无意义。
应将各种荷载放在同一种工况下进行
分析。
2.在做斜拉桥的成桥阶段分析时,如果要对自重+二期荷载+初拉力的荷载组合进行几何非线
性分析,建议在建立索单元时,直接在Lu/L列中选择初拉力后,在右侧输入初拉力值。
3.在做几何非线性分析时,尽量不要使用释放梁端约束功能。
建议在需要释放梁端约束的节
点位置建立两个节点,节点间用弹性连接连接,在需要释放约束的方向不输入刚度值即可。
几何非线性分析中不推荐使用释放梁端约束的原因如下:
a.几何非线性分析中的单元几何刚度是使用节点的坐标计算的(不断修正)。
b.释放梁端约束后,因为梁单元的端部节点和另一个单元的节点共享一个节点,单元的
位移和节点的位移会有不一致的问题,从而造成几何非线性分析不容易收敛。
4.几何非线性分析和P-delta分析不能同时进行,几何非线性分析属于大位移分析,P-delta
分析属于小位移分析,做几何非线性分析时不必再做P-delta分析。
如果用户这两个分析都要做,建议另存模型后分别分析。
5.几何非线性分析中采用的方法有TL法、UL法、CR法等,MIDAS中采用的是CR法。
6.MIDAS中几何非线性分析中适用的单元有桁架(包括索)、梁单元、平面应力单元、板单元,
如果与其他单元(如实体单元)混合使用时,只能考虑其他单元的刚度效应,不能考虑其他单元的几何非线性效应。
7.在Civil 671版本中,施工阶段分析可同时考虑非线性累加模型和收缩和徐变分析(同时勾
选即可),即非线性分析的累加模型可以考虑收缩和徐变。
考虑非线性分析(独立模型):考虑几何非线性进行施工阶段分析。
独立模型阶段:在各个施工阶段形成独立模型来进行分析。
独
立模型的几何非线性分析和考虑时间依存特性的分析不能同时
进行。
此时,除了非线性分析控制选项之外的其它选项不能设
定。
累加模型阶段:累加各个施工阶段的结果来进行非线性分析。
进行累加模型的几何非线性分析时,可以考虑时间依存特性的
效果和索初拉力类型,还可以考虑施工阶段新激活构件的初始
切向位移(包括未闭合配合力)。
考虑平衡构件内力:考虑平衡力进行非线性分析。
只考虑P-Delta效应:考虑P-Delta效应进行施工阶段分析。
能够考虑轴力对
考虑时间依存效果(累加模型):考虑材料的徐变和收缩、抗压强度的变化。
注:
1.累加模型阶段的几何非线性分析必须是以实际位移
为基准进行的,所以必须要勾选“赋予施工阶段中新
激活构件初始切向位移”选项。
就算用户不选择此项,
程序内部自动选择。
2.进行斜拉桥的几何非线性分析时,在激活拉索的阶段,
不能同时激活其它单元、除索张力以外其它荷载、其
它边界条件。
3.施工阶段中激活支撑(或连接、弹簧等)边界条件时,
要在相应构件激活之前一个施工阶段激活。
如果构件
和边界条件同时激活,边界条件被激活时考虑了前一
个施工阶段引起的变形,结果会有很大的误差。
所以
先激活构件的情况,在激活边界条件时要选择“变形
前”选项。
几何非线性分析ANSYS非线形分析指南几何非线形分析随着位移增长,一个有限单元已移动的坐标可以以多种方式改变结构的刚度。
一般来说这类问题总是是非线性的,需要进行迭代获得一个有效的解。
大应变效应一个结构的总刚度依赖于它的组成部件(单元)的方向和单刚。
当一个单元的结点经历位移后,那个单元对总体结构刚度的贡献可以以两种方式改变变。
首先,如果这个单元的形状改变,它的单元刚度将改变。
(看图2─1(a))。
其次,如果这个单元的取向改变,它的局部刚度转化到全局部件的变换也将改变。
(看图2─1(b))。
小的变形和小的应变分析假定位移小到足够使所得到的刚度改变无足轻重。
这种刚度不变假定意味着使用基于最初几何形状的结构刚度的一次迭代足以计算出小变形分析中的位移。
(什么时候使用“小”变形和应变依赖于特定分析中要求的精度等级。
相反,大应变分析说明由单元的形状和取向改变导致的刚度改变。
因为刚度受位移影响,且反之亦然,所以在大应变分析中需要迭代求解来得到正确的位移。
通过发出NLGEOM,ON(GUI路径Main Menu>Solution>Analysis Options),来激活大应变效应。
这效应改变单元的形状和取向,且还随单元转动表面载荷。
(集中载荷和惯性载荷保持它们最初的方向。
)在大多数实体单元(包括所有的大应变和超弹性单元),以及部分的壳单元中大应变特性是可用的。
在ANSYS/Linear Plus程序中大应变效应是不可用的。
第1页ANSYS非线形分析指南几何非线形分析图1─11 大应变和大转动大应变处理对一个单元经历的总旋度或应变没有理论限制。
(某些ANSYS单元类型将受到总应变的实际限制──参看下面。
)然而,应限制应变增量以保持精度。
因此,总载荷应当被分成几个较小的步,这可以〔NSUBST,DELTIM,AUTOTS〕,通过GUI路径Main Menu>Solution>Time/Prequent)。
钢结构的几何非线性分析在结构工程设计与研究中,几何非线性分析是一项重要的任务,特别是在钢结构的设计过程中。
钢结构的几何非线性分析考虑了结构形变和位移的影响,以更准确地评估结构的性能和稳定性。
一、概述钢结构通常由大量的钢材构件组成,这些构件经受荷载作用后会发生形变和变形。
当荷载作用超过结构的弹性极限时,结构材料开始发生非弹性变形,即产生塑性变形。
这种塑性变形会导致结构的刚度和稳定性发生变化,因此在设计过程中必须考虑几何非线性效应。
二、几何非线性分析方法1. 大位移理论大位移理论是几何非线性分析的基础理论之一。
它考虑了结构在受荷载作用下发生的大位移和大变形,能够更真实地模拟结构的实际响应。
大位移理论通过引入非线性应变和非线性应力来描述结构的变形情况,从而得到更准确的分析结果。
2. 几何非线性有限元分析几何非线性有限元分析是常用的计算方法之一。
该方法将结构离散化为有限数量的单元,并在每个单元内考虑非线性效应。
通过求解非线性方程组,可以得到结构的位移和应力分布,从而评估结构的承载能力和稳定性。
三、应用领域钢结构的几何非线性分析广泛应用于工程实践中。
以下是一些典型的应用领域:1. 结构稳定性分析钢结构在受到外部荷载作用下,可能发生稳定性失效。
几何非线性分析可以考虑结构的大位移和大变形,并通过评估结构的临界载荷以判断稳定性。
2. 构件受力分析在实际工程中,钢结构的各个构件可能存在复杂的荷载作用,如弯曲、剪切和扭转等。
几何非线性分析可以考虑这些复杂的受力情况,从而准确评估构件的受力性能。
3. 地震响应分析钢结构在地震荷载下会发生较大的位移和变形,甚至可能发生破坏。
几何非线性分析可以模拟结构在地震作用下的响应,评估结构的安全性。
四、结论钢结构的几何非线性分析是设计和评估钢结构性能的重要手段。
通过考虑结构的大位移和大变形效应,可以更准确地预测结构的响应和稳定性。
在实际工程中,几何非线性分析应用广泛,涵盖了结构稳定性、构件受力分析和地震响应分析等方面。
钢结构的非线性分析钢结构作为一种重要的结构形式,在建筑和工程领域被广泛应用。
而在设计和分析这类结构时,非线性分析是不可或缺的一部分。
本文将围绕钢结构的非线性分析展开讨论,并就该主题进行全面的阐述。
一、引言钢结构的非线性分析是指在考虑结构材料和结构构件在受荷过程中的非线性特性的条件下,对结构的变形、承载力和稳定性进行分析。
与线性分析相比,非线性分析更为精确,能够更好地反映实际结构的力学行为。
因此,在实际工程设计中,钢结构的非线性分析具有重要意义。
二、非线性分析的类型1. 几何非线性分析几何非线性分析是指在受荷过程中,结构的几何形状发生较大变形时的分析方法。
在传统线性分析中,通常假设结构的变形是较小的,而几何非线性分析则能更准确地考虑结构变形对力学特性的影响。
2. 材料非线性分析材料非线性分析是指考虑结构材料在受荷过程中的非线性特性进行的分析。
钢材的应力-应变曲线在高应力水平下表现出明显的非线性特性,材料非线性分析能更真实地模拟实际情况,确保结构的安全性。
3. 接触非线性分析钢结构中的接触问题也是需要考虑的一个重要方面。
接触非线性分析是指在考虑结构构件之间接触和摩擦时进行的分析。
通过准确分析接触问题,可以更精确地确定结构的承载能力和变形情况。
三、非线性分析的数值方法为了实现钢结构的非线性分析,需要借助于数值计算方法。
目前常用的数值方法包括有限元法、非线性弹性法和塑性铰接法等。
1. 有限元法有限元法是一种将结构划分为许多小单元,通过对这些小单元的力学特性进行分析,再综合考虑整体的力学性能的分析方法。
对于钢结构的非线性分析,有限元法能够较准确地考虑结构材料和几何的非线性特性。
2. 非线性弹性法非线性弹性法是基于弹性理论的扩展,通过引入非线性材料的应力-应变关系进行分析。
该方法适用于分析较小变形下的结构非线性行为。
3. 塑性铰接法塑性铰接法是一种将钢材的塑性行为简化为铰节点模型的分析方法。
通过确定铰节点的位置和性能,可以快速而准确地分析钢结构的非线性特性。
解非线性方程是方法主要有:增量法、迭代法、增量迭代混合法。
几何非线性有限元方法:1、完全的拉格朗日列式法(T.L.Formulation)在整个分析过程中,以t=0时的位形作为参考,且参考位形保持不变,这种列式称为完全的拉格朗日列式(T.L法)对于任意应力-应变关系与几何运动方程,杆系单元的平衡方程可由虚功原理推导得到:式(1)式中各量分别为:应变矩阵,是单元应变与节点位移的关系矩阵;单元的应力向量;杆端位移向量;V是单元体积分域,对T.L列式,是变形前的单元体积域;单元杆端力向量;直接按上式建立单元刚度方程并建立结构有限元列式,称为全量列式法。
在几何非线性分析中,按全量列式法得到的单元刚度矩阵和结构刚度矩阵往往是非对称的,对求解不利,因此多采用增量列式法。
将式(1)写成微分形式变形后得:式(2)这就是增量形式T.L列式的单元平衡方程。
式中为:单元弹性刚度矩阵、单元初位移刚度矩阵或单元大位移刚度矩阵、初应力刚度矩阵、三个刚度矩阵之和,称为单元切线刚度矩阵。
2、修正的拉格朗日列式法(U.L.Formulation)在建立t+∆t时刻物体平衡方程时,如果我们选择的参照位形不是未变形状态t=0时的位形,而是最后一个已知平衡状态,即本增量步起始的t时刻位形为参照位形,这种列式法称为修正的拉格朗日列式法(U.L列式)。
增量形式的U.L列式结构平衡方程可写成:式(3)3、T.L列式与U.L列式的比较T.L列式与U.L列式是不同学派用不同的简化方程及理论导出的不同方法,但是它们在相同的荷载增量步内其线性化的切线刚度矩阵应该相同,这一点已得到多个实际例题的证明。
T.L列式与U.L列式的不同点比较内容| T.L列式| U.L列式| 注意点计算单刚的积分域| 在初始构形的体积域内进行| 在变形后的t时刻体积域内进行| U.L列式必须保留节点坐标值精度| 保留了刚度阵中所有线性与非线性项| 忽略了高阶非线性| U.L列式的荷载增量不能过大单刚组集成总刚| 用初始时刻各单元结构总体坐标系中的方向余弦形成转换阵,计算过程不变| 用变形后t时刻单元在结构总体坐标中的方向余弦形成转换阵,计算过程中不断改变| U.L列式中组集荷载向量也必须注意方向余弦的改变本构关系的处理| 在大应变时,非线性本构关系不易引入| 比较容易引入大应变非线性本构关系| U.L方法更适用于混凝土徐变分析从理论上讲,这这两种方法都可以用于各种几何非线性分析。
midas,civil可以分析材料非线性与几何非线性篇一:midaS几何非线性理论知识当结构的变形相对杆件长度已不能忽略时,为了在结构变形后的形状上建立平衡,并考虑初始缺陷对结构屈曲承载力的影响,必须对结构进行基于大挠度理论的非线性屈曲分析。
在midas中可以这样处理:对于索结构或张悬梁结构中,定义的只受拉索单元并不能进行特征值分析,因为其只能定义在几何非线性分析中。
如要进行特征值分析,那么要将只受拉索单元转换为只受拉桁架单元。
先对该结构进行几何非线性,得出自重作用下的初始索力,然后将索单元定义为只受拉桁架单元,将计算所得的索力按初始荷载加到单元中:荷载->初始荷载->小位移->初始单元内力加入张力。
1、问:在midaS中如何计算自重作用下活荷载的稳定系数(屈曲分析安全系数)?答:稳定分析又叫屈曲分析,所谓的荷载安全系数(临界荷载系数)均是对应于某种荷载工况或荷载组合的。
例如:当有自重w和集中活荷载P作用时,屈曲分析结果临界荷载系数为10的话,表示在10*(w+P)大小的荷载作用下结构可能发生屈曲。
但这也许并不是我们想要的结果。
我们想知道的是在自重(或自重+二期恒载)存在的情况下,多大的活荷载作用下会发生失稳,即想知道w+Scale*P中的Scale值。
我们推荐下列反复计算的方法。
步骤一:先按w+P计算屈曲分析,如果得到临街荷载系数S1。
步骤二:按w+S1*P计算屈曲,得临界荷载系数S2。
步骤二:按w+S1*S2*P计算屈曲,得临界荷载系数S3。
重复上述步骤,直到临街荷载系数接近于1.0,此时的S1*S2*S3*Sn 即为活荷载的最终临界荷载系数。
(参见下图)midas官方网站的说话,供大家参考:考虑几何非线性同时进行稳定分析可以实现。
方法如下:1、将进行稳定分析所用荷载定义在一个荷载工况下;2、定义非线性分析控制,选择几何非线性,在非线性分析荷载工况中添加此荷载工况,并对其定义加载步骤;3、分析;4、查看结果中的阶段步骤时程图表,查找变形发生突变的位置点,及加载系数,即可推知发生失稳的极限荷载。
midascivil可以分析材料非线性与几何非线性篇一:midaS几何非线性理论知识当结构的变形相对杆件长度已不能忽略时,为了在结构变形后的形状上建立平衡,并考虑初始缺陷对结构屈曲承载力的影响,必须对结构进行基于大挠度理论的非线性屈曲分析。
在midas中可以这样处理:对于索结构或张悬梁结构中,定义的只受拉索单元并不能进行特征值分析,因为其只能定义在几何非线性分析中。
如要进行特征值分析,那么要将只受拉索单元转换为只受拉桁架单元。
先对该结构进行几何非线性,得出自重作用下的初始索力,然后将索单元定义为只受拉桁架单元,将计算所得的索力按初始荷载加到单元中:荷载->初始荷载->小位移->初始单元内力加入张力。
1、问:在midaS中如何计算自重作用下活荷载的稳定系数(屈曲分析安全系数)?答:稳定分析又叫屈曲分析,所谓的荷载安全系数(临界荷载系数)均是对应于某种荷载工况或荷载组合的。
例如:当有自重w和集中活荷载P作用时,屈曲分析结果临界荷载系数为10的话,表示在10*(w+P)大小的荷载作用下结构可能发生屈曲。
但这也许并不是我们想要的结果。
我们想知道的是在自重(或自重+二期恒载)存在的情况下,多大的活荷载作用下会发生失稳,即想知道w+Scale*P中的Scale值。
我们推荐下列反复计算的方法。
步骤一:先按w+P计算屈曲分析,如果得到临街荷载系数S1。
步骤二:按w+S1*P计算屈曲,得临界荷载系数S2。
步骤二:按w+S1*S2*P计算屈曲,得临界荷载系数S3。
重复上述步骤,直到临街荷载系数接近于1.0,此时的S1*S2*S3*Sn 即为活荷载的最终临界荷载系数。
(参见下图)midas官方网站的说话,供大家参考:考虑几何非线性同时进行稳定分析可以实现。
方法如下:1、将进行稳定分析所用荷载定义在一个荷载工况下;2、定义非线性分析控制,选择几何非线性,在非线性分析荷载工况中添加此荷载工况,并对其定义加载步骤;3、分析;4、查看结果中的阶段步骤时程图表,查找变形发生突变的位置点,及加载系数,即可推知发生失稳的极限荷载。
ansys 非线性分析原理ANSYS中的非线性分析是指通过考虑材料的非线性行为、几何非线性和边界条件的非线性等因素,对结构进行分析和计算。
非线性分析的原理主要包括以下几个方面。
1. 材料的非线性行为:考虑到材料在受载作用下的非线性行为,一般采用弹塑性分析方法。
弹塑性材料在受力时会出现应力-应变曲线的非线性特征,这需要使用合适的本构模型来描述。
ANSYS中常用的本构模型有弹塑性模型、弹性模型等,根据问题的实际情况选择适当的本构模型进行分析。
2. 几何的非线性效应:当结构在受载作用下出现较大的变形时,就需要考虑几何非线性效应。
一般情况下,当结构的变形较小时可以忽略几何非线性,反之则需要进行几何非线性分析。
几何非线性的分析可通过使用大变形理论来描述结构的非线性变形,并进行相应的计算。
3. 边界条件的非线性效应:非线性分析还需要考虑边界条件的非线性效应。
在实际工程中,边界条件往往是随着结构的变形而变化的,如约束条件的变化、边界载荷的变化等。
这些非线性边界条件会对结构的响应产生影响,因此需要将其考虑在内进行非线性分析。
在ANSYS中进行非线性分析时,通常需要进行以下步骤:1. 定义材料的本构模型:选择合适的弹塑性模型或弹性模型,并设置相应的参数。
2. 构建几何模型:根据实际工程要求,构建结构的几何模型,并对其进行离散化,即将结构分割成有限元网格。
3. 施加边界条件和载荷:根据实际工况,为结构施加边界条件和载荷。
4. 求解非线性方程组:通过非线性方程的迭代求解方法,求解得到结构的非线性响应。
5. 分析结果的后处理:对求解得到的结果进行分析和后处理,获取所需的工程参数和信息。
总之,非线性分析在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)。
无论何时当系统是非保守系统,来自动实现如在模型中有塑性或摩擦,或者有多个大位移解存在,如具有突然转换现象,使用小的载荷增量具有双重重要性。
关于大应变的特殊建模讨论应力─应变在大应变求解中,所有应力─应变输入和结果将依据真实应力和真实(或对数)应变。
(一维时,真实应变将表求为。
对于响应的小应变区,真实应变和工程应变基本上是一致的。
)要从小工程应变转换成对数应变,使用。
要从工程应力转换成真实应力,使用。
(这种应力)转化反对不可压缩塑性应力─应变数据是有效的。
)为了得到可接受的结果,对真实应变超过50%的塑性分析,应使用大应变单元(VISCO106,107及108)。
单元的形状应该认识到在大应变分析的任何迭代中低劣的单元形状(也就是,大的纵横比,过度的顶角以及具有负面积的已扭曲单元)将是有害的。
因此,你必须和注意单元的原始形状一样注意的单元已扭曲的形状。
(除了探测出具有负面积的单元外,ANSYS程序对于求解中遇到的低劣单元形状不发出任何警告,必须进行人工检查)如果已扭曲的网格是不能接受的,可以人工改变开始网格(在容限内)以产生合理的最终结果(参看图2─2)。
图2─2在大应变分析中避免低劣单元形状的发展具有小应变的大偏移小应变大转动某些单元支持大的转动,但不支持大的形状改变。
一种称作大挠度的大应变特性的受限形式对这类单元是适用的。
在一个大挠度分析中,单元的转动可以任意地大,但是应变假定是小的。
大挠度效应(没有大的形状改变)在ANSYS/Linear Plus程序中是可用的。
(在ANSYS/Mechanical,以及ANSYS/Structural产品中,对于支持大应变特性的单元,大挠度效应不能独立于大应变效应被激活。
)在所有梁单元和大多数壳单元中,以及许多非线性单元中这个特性是可用的。
通过打开NLGEOM,ON (GUI路径Main Menu>Solution>Anolysis Options)来激活那些支持这一特性的单元中的大位移效应。
应力刚化结构的面外刚度可能严重地受那个结构中面内应力的状态的影响。
面内应力和横向刚度之间的联系,通称为应力刚化,在薄的,高应力的结构中,如缆索或薄膜中,是最明显的。
一个鼓面,当它绷紧时会产生垂向刚度,这是应力强化结构的一个普通的例子。
尽管应力刚化理论假定单元的转动和应变是小的,在某些结构的系统中(如在图2─3(a)中),刚化应力仅可以通过进行大挠度分析得到。
在其它的系统中(如图2─3(b)中),刚化应力可采用小挠度或线性理论得到。
图2─3 应力硬化梁要在第二类系统中使用应力硬化,必须在第一个载荷步中发出SSTIF,ON (GUI路径Main Menu>Solution>Analysis Options)。
ANSYS程序通过生成和使用一个称作“应力刚化矩阵”的辅助刚度矩阵来考虑应力刚化效应。
尽管应力刚度矩阵是使用线性理论得到的,但由于应力(应力刚度矩阵)在每次迭代之间是变化的这个事实因而它是非线性的。
大应变和大挠度处理包括进初始应力效应作为它们的理论的一个子集,对于许多实体和壳单元,当大变型效应被激活时〔NLGEOM,ON〕(GUI路径Main Menu>Solution>Analysis Options)自动包括进初始硬化效应。
在大变形分析中〔NLGEOM,ON〕包含应力刚化效应〔SSTIF,ON〕将把应力刚度矩阵加到主刚度矩阵上以在具有大应变或大挠度性能的大多数单元中产生一个“近似的”协调切向刚度矩阵。
例外情况包括BEAM4和SHELL63,以及不把“应力刚化”列为特殊特点的任何单元。
对于BEAM4和SHELL63,你可以通过设置KEYOPT(2)=1和NLGEOM,ON在初始求解前激活应力刚化。
当大变形效应为ON(开)时这个KEYOPT设置激活一个协调切向刚度矩阵选项。
当协调切向刚度矩阵被激活时(也就是,当KEYOPT(2)=1且NLGEOM,ON 时)SSTIF对BEAM4和SHELL63将不起作用。
在大变型分析中何时应当使用应力刚化·对于大多数实体单元,应力刚化的效应是与问题相关的;在大变型分析中的应用可能提高也可能降低收敛性。
在大多数情况下,首先应该尝试一个应力刚化效应OFF(关闭)的分析。
如果你正在模拟一个受到弯曲或拉伸载荷的薄的结构,当用应力硬化OFF(关)时遇到收敛困难,则尝试打开应力硬化。
·应力刚化不建议用于包含“不连续单元”(由于状态改变,刚度上经历突然的不连续变化的非线性单元,如各种接触单元,SOLID65,等等)的结构。
对于这样的问题,当应力刚化为ON(开)时,结构刚度上的不连续线性很容易导致求解“胀破”。
·对于桁、梁和壳单元,在大挠度分析中通常应使用应力刚化。
实际上,在应用这些单元进行非线性屈曲和后屈曲分析时,只有当打开应力刚化时才得到精确的解。
(对于BEAM4和SHELL63,你通过设置单元KEYOPT(2)=1激活大挠度分析中〔NLGEOM,ON〕的应力刚化。
)然而,当你应用杆、梁或者壳单元来模拟刚性连杆,耦合端或者结构刚度的大变化时,你不应使用应力刚化。
注意:无论何时使用应力刚化,务必定义一系列实际的单元实常数。
使用不是“成比例”(也就是,人为的放大或缩小)的实常数将影响对单元内部应力的计算,且将相应地降低那个单元的应力刚化效应。
结果将是降低解的精度。
旋转软化旋转软化为动态质量效应调整(软化)旋转物体的刚度矩阵。
在小位移分析中这种调整近似于由于大的环形运动而导致几何形状改变的效应。
通常它和预应力[PSTRES](GUI路径Main Menu>Solution>Analysis Options)一起使用,这种预应力由旋转物体中的离心力所产生。
它不应和其它变形非线性,大挠度和大应变一起使用。
旋转软化用OMEGA命令中的KPSIN来激活(GUI路径Main Menu>Preprocessor>Loads>-Loads-Apply>-Structural-Other>Angular Velotity)。
关于非线性分析的忠告和准则着手进行非线性分析通过比较小心地采用时间和方法,可以避免许多和一般的非线性分析有关的困难,下列建议对你可能是有益的了解程序的运作方式和结构的表现行为如果你以前没有使用过某一种特别的非线性特性,在将它用于大的,复杂的模型前,构造一个非常简单的模型(也就是,仅包含少量单元),以及确保你理解了如何处理这种特性。
·通过首先分析一个简化模型,以便使你对结构的特性有一个初步了解。
对于非线性静态模型,一个初步的线性静态分析可以使你知道模型的哪一个区域将首先经历非线性响应,以及在什么载荷范围这些非线性将开始起作用。
对于非线性瞬态分析,一个对梁,质量块及弹簧的初步模拟可以使你用最小的代价对结构的动态有一个深入了解。
在你着手最终的非线性瞬时动态分析前,初步非线性静态,线性瞬时动态,和/或模态分析同样地可以有助于你理解你结构的非线性动态响应的不同的方面。
·阅读和理解程序的输出信息和警告。
至少,在你尝试后处理你的结果前,确保你的问题收敛。
对于与路程相关的问题,打印输出的平衡迭代记录在帮助你确定你的结果是有效还是无效方面是特别重的。
简化·尽可能简化最终模型。
如果可以将3─D结构表示为2─D平面应力,平面应变或轴对称模型,那么这样做,如果可以通过对称或反对称表面的使用缩减你的模型尺寸,那么这样做。
(然而,如果你的模型非对称加载,通常你不可以利用反对称来缩减非线性模型的大小。
由于大位移,反对称变成不可用的。
)如果你可以忽略某个非线性细节而不影响你模型的关键区域的结果,那么这样做。
·只要有可能就依照静态等效载荷模拟瞬时动态加载。
·考虑对模型的线性部分建立子结构以降低中间载荷或时间增量及平衡迭代所需要的计算时间。
采用足够的网格密度·考虑到经受塑性变形的区域要求一个合理的积分点密度。
每个低阶单元将提供和高阶单元所能提供的一样多积分点数,因此经常优先用于塑性分析。
在重要塑性区域网格密度变得特别地重要,因为大挠度要求对于一个精确的解,个单元的变形(弯曲)不能超过30度。
·在接触表面上提供足够的网格密度以允许接触应力以一种平滑方式分布。
·提供足够用于分析应力的网格密度。
那些应力或应变关心的面与那些需要对位移或非线性解析处的面相比要求相对好的网格。
·使用足够表征最高的重要模态形式的网格密度。
所需单元数目依赖于单元的假定位移形状函数,以及模态形状本身。
·使用足够可以用来分析通过结构的任何瞬时动态波传播的网格密度。
如果波传播是重要的,那么至少提供20个单元来分析一个波长。
逐步加载·对于非保守的,与路径相关的系统,你需要以足够小的增量施加载荷以确保你的分析紧紧地跟随结构的载荷响应曲线。