【非常好】应力刚化和几何刚度

应力刚化和几何刚度

一、应力刚化ANSYS程序通过生成和使用一个称作“应力刚化矩阵”的辅助刚度矩阵来考虑应力刚化效应。

在大变形分析中何时使用应力刚化：

1、对于大多数实体单元，应力刚化的效应是与问题相关的；在大变形分析中的应用可能提高也可能降低收敛性。在大多数情况下，应该首先尝试一个应力刚化效应OFF的分析，如果遇到收敛困难时，则尝试打开应力刚化。

2、应力刚化不适应于包含由于状态改变、刚度上经历突然的不连续变化的非线性单元的结构。对于这样的结构，当应力刚化效应打开时，结构刚度上的不连续性很容易导致求解的“胀破”。

3、对于梁和壳单元，在大挠度分析中通常应该使用应力刚化。实际上，在应用这些单元进行非线性屈曲和后屈曲分析时，只有打开应力刚化才能得到精确的解。但当应用杆、梁或者壳单元来模拟刚性连杆、耦合端或者结构刚度的大变化时，不应该使用应力刚化效应。

4、无论何时使用应力刚化，务必定义一系列实际的单元实常数。使用不是人为的放大和缩小的实常数将影响对单元内部应力的计算，且将相应地降低那个单元的应力刚化效应，结果将是降低解的精度。

二、几何刚度

几何刚度矩阵表示结构在变形状态下的刚度变化，与施加的荷载有直接的关系。任意构件受到压力时，刚度有减小的倾向；反之，受到拉力时，刚度有增大的倾向。考虑几何非线性的大变形结构分析，屈曲分析等都要考虑几何刚度矩阵。例如求临界荷载P（特征值）的屈曲分析平衡方程：([K0]+P*[Kg])*{U}={0}

[K0] ：结构的弹性刚度矩阵

[Kg] ：结构的几何刚度矩阵

要使{U}有非0解，{U}的系数行列式为0，即|[K]+P*[Kg]|＝0

几何刚度矩阵又称为初应力刚度矩阵，与Ansys中称之为应力硬化的现象有关。对于梁杆体系而言，应力硬化实际上就是P-Δ效应。

应力硬化具体可参见ANSYS, Inc. Theory Reference中的3.3. Stress Stiffening。这里简述如下：应力硬化(亦称为几何硬化、增量硬化、初应力硬化和微分硬化)，是由于结构的应力状态引起结构的强化或者软化。通常存在于弯曲刚度相对轴向刚度很小的薄结构，如索、膜、梁、壳等。该效应亦可能是由大应变或者大变形引起。几何刚阵是通过前一个平衡迭代的应力状态来计算的，因此至少要迭代2次。从上可知，引起应力硬化的情况都要考虑几何刚度，如小变形条件下的P-Δ效应等。大变形情况下一般要考虑几何刚度，当然也不是必须的，Ansys 中大变形打开(NLGEOM,ON)时，同时会打开应力硬化(SSTIF,ON)，但用户也可以选择关闭。pillow兄考虑的是小变形情况下的P-Δ效应，对于无侧移刚架会出现特征值问题。但一般情况下，荷载右端项不为0，此时可以迭代计算求出荷载极值。

应力硬化理论假定单元的转动和应变是微小的，在某些结构的系统中，硬化应力仅可以通过进行大绕度分析得到。有些系统中，也可以采用小绕度或线性理论得到。如果采用小绕度或线性理论则必须在第一个载荷步中使用命令SSTIF ON。ANSYS程序通过生成和使用一个称作“应力硬化矩阵”的辅助刚度矩阵来考虑应力硬化效应。尽管刚度矩阵是使用线性理论得到的，但由于应力或应力矩阵在每次迭代之间是变化的，因此它仍旧是非线性的。在ANSYS 程序的大应变和大绕度处理中，一般都考虑到初始应力效应的影响，将其作为大应变和大绕

度理论的一个子集。对于许多实体单元和壳单元来说，当大变形效应被激活时，程序将自动包括初始硬化效应。在大变形分析中应力硬化效应的加入，是通过把应力刚度矩阵加到主刚度矩阵上以在具有大应变或大绕度性能的多数单元中产生一个“近似”的协调切向刚度矩阵。==============================================================================

几何非线性分析应力刚化随着位移增长一个有限单元已移动的坐标可以以多种方式改变结构的刚度一般来说这类问题总是是非线性的需要进行迭代获得一个有效的解大应变效应

一个结构的总刚度依赖于它的组成部件单元的方向和单刚当一个单元的结点经历位移后那个单元对总体结构刚度的贡献可以以两种方式改变，首先如果这个单元的形状改变它的单元刚度将改变看图2─1(a)，其次如果这个单元的取向改变它的局部刚度转化到全局部件的变换也将改变看图2─1b)小的变形和小的应变分析假定位移小到足够使所得到的刚度改变无足轻重这种刚度不变假定意味着使用基于最初几何形状的结构刚度的一次迭代足以计算出小变形分析中的位移什么时候使用小变形和应变依赖于特定分析中要求的精度等级相反大应变分析说明由单元的形状和取向改变，导致的刚度改变，因为刚度受位移影响且反之亦然所以在大应变分析中需要迭代求解来得到正确的位移通过发出NLGEOMONGUI 路径Main Menu>Solution>Analysis Options)来激活大应变效应这效应改变单元的形状和取向且还随单元转动表面载荷集中载荷和惯性载荷保持它们最初的方向在大多数实体单元包括所有的大应变和超弹性单元以及部分的壳单元中大应变特性是可用的在ANSYS/Linear Plus 程序中大应变效应是不可用的

图1─11 大应变和大转动

大应变处理对一个单元经历的总旋度或应变没有理论限制某些ANSYS单元类型将受到总应变的实际限制──参看下面然而应限制应变增量以保持精度因此总载荷应当被分成几个较小的步这可以NSUBSTDELTIMAUTOTS通过GUI路径Main Menu>Solution>Time/Prequent)无论何时当系统是非保守系统来自动实现如在模型中有塑性或摩擦或者有多个大位移解存在如具有突然转换现象使用小的载荷增量具有双重重要性

关于大应变的特殊建模讨论

应力─应变

在大应变求解中所有应力─应变输入和结果将依据真实应力和真实或对数应变一维时真实应变将表求为对于响应的小应变区真实应变和工程应变基本上是一致的要从小工程应变转换成对数应变使用要从工程应力转换成真实应力使用这种应力转化反对不可压缩塑性应力─应变数据是有效的

为了得到可接受的结果对真实应变超过50%的塑性分析应使用大应变单元

大应变与小

应变分析的

界定

ANSYS非线形分析指南几何非线形分析

VISCO106107及108

单元的形状

应该认识到在大应变分析的任何迭代中低劣的单元形状也就是大的纵横比过度的顶角以及具有负面积的已扭曲单元将是有害的因此你必须和注

意单元的原始形状一样注意的单元已扭曲的形状除了探测出具有负面积的单元外ANSYS程序对于求解中遇到的低劣单元形状不发出任何警告必须进行人工检查如果已扭曲的网格是不能接受的可以人工改变开始网格在容限内以产生合理的最终结果参看图2─2 图2─2 在大应变分析中避免低劣单元形状的发展具有小应变的大偏移