第三章几何非线性分解.ppt

第三章第三章 几何非线性几何非线性大应变，大位移与大旋转特性大应变，大位移与大旋转特性September 30, 19983-1Structural Nonlinearities - Revision 5.5 (001156) 什么是几何非线性？什么是几何非线性？什么是几何非线性？什么是几何非线性？•变形体几何形态的改变将明显影响物体的载荷－位移（如刚度）变形体几何形态的改变将明显影响物体的载荷－位移（如刚度）特性•几何非线性并不只是指大位移，而且还包括几何状态改变所引起几何非线性并不只是指大位移，而且还包括几何状态改变所引起的任何结构响应的变化它包括大应变、大位移和大旋转的任何结构响应的变化它包括大应变、大位移和大旋转October 17, 20002几何非线性 – 5.7版本 几何非线性特性几何非线性特性几何非线性特性几何非线性特性•如果一个单元的形状发生改变（面积、厚度等），它本身的单如果一个单元的形状发生改变（面积、厚度等），它本身的单元矩阵会发生改变元矩阵会发生改变•如果一个单元的取向改变，它的单元刚度向整体刚度的转换矩如果一个单元的取向改变，它的单元刚度向整体刚度的转换矩阵将发生变化。

阵将发生变化XYXYOctober 17, 20003几何非线性 – 5.7版本 几何非线性特性（续几何非线性特性（续几何非线性特性（续几何非线性特性（续) )•如果单元应变产生了明显的面内应力如果单元应变产生了明显的面内应力 (膜应力膜应力)，垂直于面的，垂直于面的刚度会明显受影响刚度会明显受影响XYP•随着垂直位移的增加随着垂直位移的增加 (Y)，，大的膜应力大的膜应力(SX)导致刚化响应导致刚化响应UYPOctober 17, 20004几何非线性 – 5.7版本 几何非线性例子几何非线性例子几何非线性例子几何非线性例子•这是一个大应变分析的例子，一个轴对称的橡胶密封件受压缩分析包这是一个大应变分析的例子，一个轴对称的橡胶密封件受压缩分析包括接触，当密封件折叠时会发生自身接触括接触，当密封件折叠时会发生自身接触October 17, 20005几何非线性 – 5.7版本 几何非线性例子几何非线性例子几何非线性例子几何非线性例子( (续续续续) )•此例显示了绕轴线捆扎一根钢条将金属弯曲成不同的形状是生产中常此例显示了绕轴线捆扎一根钢条将金属弯曲成不同的形状是生产中常见的操作。

在此例中，应变达到见的操作在此例中，应变达到25%，，顶端旋转接近顶端旋转接近270 度！度！ 初始形状初始形状 变形形状变形形状 October 17, 20006几何非线性 – 5.7版本 几何非线性问题的复杂性几何非线性问题的复杂性几何非线性问题的复杂性几何非线性问题的复杂性•追随力追随力•非线性应力应变度量非线性应力应变度量•一致的非线性刚度矩阵一致的非线性刚度矩阵•材料非线性材料非线性•不可压缩性不可压缩性•网格扭曲网格扭曲October 17, 20007几何非线性 – 5.7版本 注意当结构经历大位移与旋转时，载荷发生了什么变化注意当结构经历大位移与旋转时，载荷发生了什么变化•在许多情况下，载荷将在变形过程中保持一致的方向在许多情况下，载荷将在变形过程中保持一致的方向•在另一些情况下，承受大旋转时，力将在另一些情况下，承受大旋转时，力将“ 追随追随”单元改变单元改变方向这两种情况这两种情况ANSYS都可模拟，取决于施加载荷的类型都可模拟，取决于施加载荷的类型追随力与载荷方向追随力与载荷方向追随力与载荷方向追随力与载荷方向October 17, 20008几何非线性 – 5.7版本 载荷方向载荷方向载荷方向载荷方向 加速度与集中力：加速度与集中力：•保持它们的初始方向，忽略单元取向。

保持它们的初始方向，忽略单元取向表面压力载荷：表面压力载荷：•随单元旋转，因此总是垂直于变形单元的表面（这些是真正的追随单元旋转，因此总是垂直于变形单元的表面（这些是真正的追随力）•更新压力表面以计算大应变效应因此，对于施加的常压力，总更新压力表面以计算大应变效应因此，对于施加的常压力，总压力载荷将随表面积的改变而改变压力载荷将随表面积的改变而改变October 17, 20009几何非线性 – 5.7版本 载荷方向载荷方向载荷方向载荷方向( (续续续续) )载荷载荷位移前的方向位移前的方向位移后的方向位移后的方向加速度加速度节点力节点力单元压力单元压力October 17, 200010几何非线性 – 5.7版本 什么是大应变？什么是大应变？什么是大应变？什么是大应变？•大应变分析设定应变不再是无限小的，而是有限的或相当大的大应变分析设定应变不再是无限小的，而是有限的或相当大的•当应变超过一定百分比及不能忽视几何形状的改变时，可认为是当应变超过一定百分比及不能忽视几何形状的改变时，可认为是大应变•大应变理论考虑了形状的改变（例如厚度，面积等等）及任何大大应变理论考虑了形状的改变（例如厚度，面积等等）及任何大旋转。

旋转October 17, 200011几何非线性 – 5.7版本 定义应变的要求定义应变的要求定义应变的要求定义应变的要求•应变是我们描述物体应变是我们描述物体变形变形 的手段尽管在一定程度上应变的数学的手段尽管在一定程度上应变的数学定义可以任意，但它仍必须满足一些要求定义可以任意，但它仍必须满足一些要求October 17, 200012几何非线性 – 5.7版本 定义应变的要求定义应变的要求定义应变的要求定义应变的要求( (续续续续) )•无变形时应变为零，有变形时不为零无变形时应变为零，有变形时不为零•通过材料的应力应变关系联系应力和应变通过材料的应力应变关系联系应力和应变•如果材料旋转大，但发生的是刚体位移，则应变为零如果材料旋转大，但发生的是刚体位移，则应变为零October 17, 200013几何非线性 – 5.7版本 应力与应变定义相对应应力与应变定义相对应应力与应变定义相对应应力与应变定义相对应•对应于度量物体变形的应变，应力度量了物体单位面积上的力对应于度量物体变形的应变，应力度量了物体单位面积上的力•尽管使用的应变可以是任意形式的，但使用的应力应与已定义的尽管使用的应变可以是任意形式的，但使用的应力应与已定义的应变相关。

应变相关October 17, 200014几何非线性 – 5.7版本 应力应变的共轭性应力应变的共轭性应力应变的共轭性应力应变的共轭性•在非线性大应变分析中使用的应力必须与应变共在非线性大应变分析中使用的应力必须与应变共轭轭 •共轭性共轭性是指应力应变相乘时，得到的是一个标量（应变能），其是指应力应变相乘时，得到的是一个标量（应变能），其值与选择的应力应变无关值与选择的应力应变无关October 17, 200015几何非线性 – 5.7版本 一维问题中的应力应变定义一维问题中的应力应变定义一维问题中的应力应变定义一维问题中的应力应变定义 •我们将通过一个简单的一维例题检查不同的应力应变我们将通过一个简单的一维例题检查不同的应力应变定义FOctober 17, 200016几何非线性 – 5.7版本 工程应变工程应变工程应变工程应变•工程应变工程应变 是是小（无穷小）应变度量，小（无穷小）应变度量，使用初始的几何使用初始的几何构形计算构形计算•工程应变取决于初始长度，如工程应变取决于初始长度，如 ，初始长度是事先已知，初始长度是事先已知的，所以工程应变是线性的的，所以工程应变是线性的。

•工程应变的应用局限于材料小旋转，中等大小的刚体旋工程应变的应用局限于材料小旋转，中等大小的刚体旋转将导致非零应变转将导致非零应变October 17, 200017几何非线性 – 5.7版本 •工程应变工程应变 的共轭应力为工程应力的共轭应力为工程应力 ，工程应力的计，工程应力的计算是当前的力算是当前的力 除以初始面积除以初始面积 工程应力工程应力工程应力工程应力October 17, 200018几何非线性 – 5.7版本 一维问题中，对数应变一维问题中，对数应变 由下面公式计算：由下面公式计算：对数应变是对数应变是非线性非线性 应变，因为它是未知的最终长度应变，因为它是未知的最终长度 的非线性函的非线性函数它同样可称为数它同样可称为log应变应变 三维等效对数应变是三维等效对数应变是Hencky应变在大应变问题中，对数应变并不能自动适应任意大的旋转在大应变问题中，对数应变并不能自动适应任意大的旋转对数应变对数应变对数应变对数应变October 17, 200019几何非线性 – 5.7版本 真实应力或真实应力或真实应力或真实应力或 Cauchy Cauchy 应力应力应力应力与对数应变与对数应变 共轭的一维应力是真实应力共轭的一维应力是真实应力 ,真实应力的计算是当真实应力的计算是当前的力前的力 除以当前（或变形的）面积除以当前（或变形的）面积 ：：真实应力通常也称为真实应力通常也称为 Cauchy 应力。

应力October 17, 200020几何非线性 – 5.7版本 Green-Lagrange Green-Lagrange 应变应变应变应变一维的一维的Green-Lagrange应变应变 由下面的公式计算：由下面的公式计算：此应变是非线性的，因为它取决于未知的更新长度此应变是非线性的，因为它取决于未知的更新长度 的平方 此种应变优于对数应变或此种应变优于对数应变或Hencky 应变应变 之处在于，它可自动适应之处在于，它可自动适应大应变问题中的任意大旋转大应变问题中的任意大旋转October 17, 200021几何非线性 – 5.7版本 第二第二第二第二 Piola-Kirchhoff Piola-Kirchhoff 应力应力应力应力Green-Lagrange 应变应变 的共轭应力是第二的共轭应力是第二Piola-Kirchhoff应应力力 它的一维计算公式为：它的一维计算公式为：需要注意的是此种应力几乎无实际意义为了输出，需要注意的是此种应力几乎无实际意义为了输出，ANSYS总总是将其转化为是将其转化为Cauchy 应力或真应力应力或真应力 输出。

输出October 17, 200022几何非线性 – 5.7版本 AA0P将非线性应变定义扩展至一般的三维情况将非线性应变定义扩展至一般的三维情况将非线性应变定义扩展至一般的三维情况将非线性应变定义扩展至一般的三维情况•在二维或三维问题中，当物体承受大应变变形时，不只长度发生在二维或三维问题中，当物体承受大应变变形时，不只长度发生改变，而且厚度、面积与体积都发生改变改变，而且厚度、面积与体积都发生改变October 17, 200023几何非线性 – 5.7版本 运动与变形运动与变形运动与变形运动与变形•当物体承受一些外载时，它将移动和变形当物体承受一些外载时，它将移动和变形•如果我们观察物体上一个点的运动，它的初始位置是如果我们观察物体上一个点的运动，它的初始位置是 ，最终，最终位置是位置是 ，它运动的量，它运动的量 为为XYOctober 17, 200024几何非线性 – 5.7版本 变形梯度变形梯度变形梯度变形梯度变形梯度是物体变形多少的一个度量，它的定义是：变形梯度是物体变形多少的一个度量，它的定义是：变形梯度变形梯度 包含的信息有：包含的信息有：•体积的改变体积的改变•旋转旋转•由于应变造成的形状改变由于应变造成的形状改变October 17, 200025几何非线性 – 5.7版本 变形梯度变形梯度•注意定义时，注意定义时，变形梯度变形梯度 消除了平动；这是应变定义的必要条消除了平动；这是应变定义的必要条件。

件• 在定义应变时，同样希望在定义应变时，同样希望排除旋转部分排除旋转部分（因为它对应变无贡献），（因为它对应变无贡献），并孤立出形状改变部分可使用并孤立出形状改变部分可使用极分解理论极分解理论完成这项工作完成这项工作October 17, 200026几何非线性 – 5.7版本 极分解极分解极分解极分解•变形梯度变形梯度 可使用极分解理论分解成一个旋转部分和一个应变可使用极分解理论分解成一个旋转部分和一个应变部分：部分：=旋转矩阵包含的信息为材料点作为刚体旋转的大小旋转矩阵包含的信息为材料点作为刚体旋转的大小与方向与方向= 拉拉伸矩阵包含的信息为物体在材料点处的应变伸矩阵包含的信息为物体在材料点处的应变October 17, 200027几何非线性 – 5.7版本 •在已知在已知拉拉伸矩阵伸矩阵 的情况下，可导出一维对数应变的情况下，可导出一维对数应变 与一维与一维Green-Lagrange 应变应变 的三维一般形式的三维一般形式在定义应变时使用在定义应变时使用在定义应变时使用在定义应变时使用 October 17, 200028几何非线性 – 5.7版本 Hencky Hencky 应变应变应变应变对数（对数（ Hencky ）应变可由下式计算：）应变可由下式计算：这里这里 是以矩阵形式表示的应变张量。

是以矩阵形式表示的应变张量在这种情况下，在这种情况下， 是一维对数应变或真实应变是一维对数应变或真实应变 的三的三维等效量维等效量October 17, 200029几何非线性 – 5.7版本 Green-Lagrange Green-Lagrange 应变应变应变应变在三维问题中，在三维问题中，Green-Lagrange 应变可直接由应变可直接由拉拉伸矩阵伸矩阵 计算计算出，如下式所示：出，如下式所示：这种这种应变在计算时直接忽略了旋转矩阵应变在计算时直接忽略了旋转矩阵 可以变形梯度可以变形梯度的形式写出，如下式所示：的形式写出，如下式所示：前两项是线性小应变项，最后一项是应变的非线性项前两项是线性小应变项，最后一项是应变的非线性项October 17, 200030几何非线性 – 5.7版本 ANSYS ANSYS 使用何种非线性应变？使用何种非线性应变？使用何种非线性应变？使用何种非线性应变？•在在ANSYS程序中，使用什么样的应变（如程序中，使用什么样的应变（如Hencky应变或应变或 Green-Lagrange 应变）应变） 主要取决于材料定律：主要取决于材料定律： –对于大应变塑性分析，对于大应变塑性分析，ANSYS使用对数（使用对数（ Hencky ）应变；）应变；应力应变数据以真应力－对数应变形式给出。

应力应变数据以真应力－对数应变形式给出–对于大应变超弹性与粘－弹性分析，对于大应变超弹性与粘－弹性分析，ANSYS使用使用Hencky 应应变或变或Green-Lagrange 应变，这取决于使用的单元类型应变，这取决于使用的单元类型October 17, 200031几何非线性 – 5.7版本 将非线性应力的定义扩展至三维问题将非线性应力的定义扩展至三维问题将非线性应力的定义扩展至三维问题将非线性应力的定义扩展至三维问题•正如在一维问题中一样，对应于二维、三维问题中定正如在一维问题中一样，对应于二维、三维问题中定义的非线性应变，都可定义与其共轭的应力义的非线性应变，都可定义与其共轭的应力October 17, 200032几何非线性 – 5.7版本 •Cauchy 或真应力张量或真应力张量 （以矩阵形式写出）（以矩阵形式写出） 给出了在变形后给出了在变形后构形每单位变形面积上的当前力如果我们让构形每单位变形面积上的当前力如果我们让•在三维问题中，在三维问题中，Cauchy 应力张量应力张量 通过下式联系通过下式联系 与与 •Cauchy 应力是一个有实际意义的量。

应力是一个有实际意义的量 变形体变形体中定义中定义单元面积单元面积的矢量的矢量= 对应的作用于面积对应的作用于面积 上的单元力上的单元力 Cauchy Cauchy 应力应力应力应力 October 17, 200033几何非线性 – 5.7版本 第二第二第二第二Piola-Kirchhoff Piola-Kirchhoff 应力应力应力应力•让让 代表从位移中导出的力代表从位移中导出的力•让让•第二第二 Piola-Kirchhoff 应力张量应力张量 使用下式联系使用下式联系 与与• 是对称的应力张量，它经常用于有限应变弹性公式，它是是对称的应力张量，它经常用于有限应变弹性公式，它是Green-Lagrange应变应变 的共轭应力张量的共轭应力张量 是一个无实际是一个无实际意义的应力张量（意义的应力张量（假想应力假想应力张量） 在在未变形未变形物体中定义物体中定义单元面积单元面积的矢量，这里的矢量，这里October 17, 200034几何非线性 – 5.7版本 与与与与 之间的关系之间的关系之间的关系之间的关系•有实际意义的有实际意义的Cauchy应力应力 可通过下式与无实际意义的第二可通过下式与无实际意义的第二 Piola-Kirchhoff 假想应力假想应力 建立直接的联系：建立直接的联系：October 17, 200035几何非线性 – 5.7版本 在非线性分析中在非线性分析中在非线性分析中在非线性分析中ANSYSANSYS使用何种应力？使用何种应力？使用何种应力？使用何种应力？•为了确保应力与应变共轭，在为了确保应力与应变共轭，在ANSYS非线性分析中使用的应力非线性分析中使用的应力度量（如度量（如Cauchy 或第二或第二Piola-Kirchhoff应力应力 ）取决于应变度）取决于应变度量。

量–当使用当使用Hencky （对数）应变时使用（对数）应变时使用 Cauchy 应力–当使用当使用Green-Lagrange 应变时使用第二应变时使用第二Piola-Kirchhoff 应力October 17, 200036几何非线性 – 5.7版本 输入时输入时输入时输入时ANSYSANSYS需要何种应力应变？需要何种应力应变？需要何种应力应变？需要何种应力应变？•对于大应变塑性分析（对于大应变塑性分析（ NLGEOM,ON ），），ANSYS需要真需要真实实应力应力应变曲线，而对于小应变分析（应变曲线，而对于小应变分析（ NLGEOM,OFF ），），ANSYS需需要工程应力应变数据要工程应力应变数据•但是，对于小应变响应，工程应变与对数（真）应变几乎一样但是，对于小应变响应，工程应变与对数（真）应变几乎一样真真实实应力与对数应变数据可用于一般塑性分析应力与对数应变数据可用于一般塑性分析•对于超弹性分析对于超弹性分析ANSYS需要工程应力应变数据来计算需要工程应力应变数据来计算Mooney-Rivlin 常数October 17, 200037几何非线性 – 5.7版本 应力应变间的转化应力应变间的转化应力应变间的转化应力应变间的转化•对于对于单轴单轴应力应变数据，工程应力应变可通过下式转化为真应力应变数据，工程应力应变可通过下式转化为真实实应力与对数应变应力与对数应变•注意上面的应力转化假设承受大应变的材料是不可压缩的或近注意上面的应力转化假设承受大应变的材料是不可压缩的或近似不可压缩的。

这种假设对于大塑性应变或超弹性材料适用似不可压缩的这种假设对于大塑性应变或超弹性材料适用October 17, 200038几何非线性 – 5.7版本 ANSYS ANSYS 输出采用何种应力形式？输出采用何种应力形式？输出采用何种应力形式？输出采用何种应力形式？•不管数值计算中使用的哪种应力（不管数值计算中使用的哪种应力（ Cauchy应力或第二应力或第二Piola-Kirchhoff 应力应力 ），），ANSYS总以实际上可解释的总以实际上可解释的Cauchy 应力形应力形式输出所有的应力结果式输出所有的应力结果October 17, 200039几何非线性 – 5.7版本 哪些单元支持大应变？哪些单元支持大应变？哪些单元支持大应变？哪些单元支持大应变？BEAM188BEAM189HYPER56HYPER58HYPER74HYPER84HYPER86PLANE2PLANE13PLANE42PLANE82PLANE182SHELL43SHELL91SHELL93SHELL181SOLID45SOLID92SOLID95SOLID185VISCO106VISCO107VISCO108October 17, 200040几何非线性 – 5.7版本 非线性全一致切向刚度矩阵非线性全一致切向刚度矩阵非线性全一致切向刚度矩阵非线性全一致切向刚度矩阵•众所周知在大多数非线性结构问题中，使用一致的非线性切向刚众所周知在大多数非线性结构问题中，使用一致的非线性切向刚度矩阵可迅速增加度矩阵可迅速增加Newton-Raphson法的收敛速度。

法的收敛速度•通常在迭代求解过程中使用一致的或全切向刚度矩阵将得到平方通常在迭代求解过程中使用一致的或全切向刚度矩阵将得到平方的收敛速度的收敛速度October 17, 200041几何非线性 – 5.7版本 •一致的非线性刚度矩阵一致的非线性刚度矩阵 是对离散的有限元方程求偏导得出的，是对离散的有限元方程求偏导得出的，它是单元内部力矢量它是单元内部力矢量 与单元外部力矢量与单元外部力矢量 的函数什么是一致的非线性刚度矩阵？什么是一致的非线性刚度矩阵？什么是一致的非线性刚度矩阵？什么是一致的非线性刚度矩阵？October 17, 200042几何非线性 – 5.7版本 离散的非线性静态有限元方程离散的非线性静态有限元方程离散的非线性静态有限元方程离散的非线性静态有限元方程离散的非线性静态有限元方程单元级的公式如下所示：离散的非线性静态有限元方程单元级的公式如下所示：这里这里 = 总体单元数量总体单元数量= 单元坐标系中单元内部力矢量单元坐标系中单元内部力矢量= 将将 转换到整体坐标系的转化矩阵转换到整体坐标系的转化矩阵= 总体坐标系中施加在单元上的载荷矢量总体坐标系中施加在单元上的载荷矢量October 17, 200043几何非线性 – 5.7版本 单元内部力矢量单元内部力矢量单元内部力矢量单元内部力矢量单元内部力矢量单元内部力矢量 可由下式给出：可由下式给出：= 单元单元 应变－节点应变－节点 位移矩阵位移矩阵= 单元应力矢量单元应力矢量= 单元体积单元体积这里这里使用上面定义的内部力形式，离散的非线性有限元（力平衡）方程可使用上面定义的内部力形式，离散的非线性有限元（力平衡）方程可重写成下式：重写成下式：October 17, 200044几何非线性 – 5.7版本 推导增量形式的非线性刚度矩阵推导增量形式的非线性刚度矩阵推导增量形式的非线性刚度矩阵推导增量形式的非线性刚度矩阵对离散的有限元方程求偏导可得到一致的非线性刚度对离散的有限元方程求偏导可得到一致的非线性刚度矩阵矩阵 如下所示如下所示这里这里October 17, 200045几何非线性 – 5.7版本 一致的非线性刚度矩阵一致的非线性刚度矩阵一致的非线性刚度矩阵一致的非线性刚度矩阵= 主切向矩阵主切向矩阵=考虑了考虑了 应力刚化应力刚化 作用的初始应力矩阵作用的初始应力矩阵= 在刚度关系中考虑了几何改变效应的初始位移－在刚度关系中考虑了几何改变效应的初始位移－ 旋转矩阵旋转矩阵 =在刚度关系中考虑了载荷取向改变（跟随力）效应在刚度关系中考虑了载荷取向改变（跟随力）效应的初始载荷矩阵的初始载荷矩阵October 17, 200046几何非线性 – 5.7版本 •下列梁单元与壳单元具有包含初始载荷矩阵（压力载荷刚度）的下列梁单元与壳单元具有包含初始载荷矩阵（压力载荷刚度）的全一致的切向刚度矩阵全一致的切向刚度矩阵 ：：–梁单元梁单元: 有限应变梁；有限应变梁；Beam188, Beam189–壳单元壳单元: 有限应变壳；有限应变壳； Shell181一致的非线性刚度矩阵一致的非线性刚度矩阵一致的非线性刚度矩阵一致的非线性刚度矩阵October 17, 200047几何非线性 – 5.7版本 压力载荷刚度压力载荷刚度压力载荷刚度压力载荷刚度•对于其他三维实体单元或壳单元，要包含压力载荷刚度计算需使用对于其他三维实体单元或壳单元，要包含压力载荷刚度计算需使用SURF154*。

•*对于二维分析，可在模型边部使用对于二维分析，可在模型边部使用Beam188 单元或单元或Beam189 单单元计算压力载荷刚度元计算压力载荷刚度 Surface effect element; SURF154 对于压力载荷刚度可对于压力载荷刚度可使用一个不对称刚度使用一个不对称刚度选项只有在收敛速选项只有在收敛速度慢时才使用度慢时才使用October 17, 200048几何非线性 – 5.7版本 材料非线性材料非线性材料非线性材料非线性•在分析大应变范围的结构行为时，需要处理非线性材料特性，如在分析大应变范围的结构行为时，需要处理非线性材料特性，如非线性弹性（超弹性）与塑性塑性与超弹性都将在其各自的章非线性弹性（超弹性）与塑性塑性与超弹性都将在其各自的章节中讨论节中讨论October 17, 200049几何非线性 – 5.7版本 不可压缩性不可压缩性不可压缩性不可压缩性•如前所述，在分析大应变范围的结构行为时需要处理非线性材如前所述，在分析大应变范围的结构行为时需要处理非线性材料特性，如非线性弹性（超弹性）与塑性料特性，如非线性弹性（超弹性）与塑性•类橡胶的超弹性材料与象塑性实体一样流动的材料经常显示出类橡胶的超弹性材料与象塑性实体一样流动的材料经常显示出一种称之为不可压缩性的实际现象（既产生应变但不发生体积一种称之为不可压缩性的实际现象（既产生应变但不发生体积改变）。

改变）October 17, 200050几何非线性 – 5.7版本 非线性材料特性中不可压缩性的原因非线性材料特性中不可压缩性的原因非线性材料特性中不可压缩性的原因非线性材料特性中不可压缩性的原因•在橡胶与类橡胶材料中，泊松比接近在橡胶与类橡胶材料中，泊松比接近0.5的自然材料呈现不可压的自然材料呈现不可压缩性•在发生塑性变形的实体中，流动定律通常不允许体积的改变在发生塑性变形的实体中，流动定律通常不允许体积的改变因此，如果塑性应变很大，在大应变范围内材料响应近似不可因此，如果塑性应变很大，在大应变范围内材料响应近似不可压缩October 17, 200051几何非线性 – 5.7版本 不可压缩与网格自锁不可压缩与网格自锁不可压缩与网格自锁不可压缩与网格自锁•不论产生不可压缩条件的原因是什么，此效应意味着在大应变分不论产生不可压缩条件的原因是什么，此效应意味着在大应变分析中使用的标准的全积分有限元可能会遇到与网格自锁相关的数析中使用的标准的全积分有限元可能会遇到与网格自锁相关的数值困难October 17, 200052几何非线性 – 5.7版本 消除网格自锁消除网格自锁消除网格自锁消除网格自锁•为了处理大应变分析中与不可压缩效应相关的网格自锁问题，可为了处理大应变分析中与不可压缩效应相关的网格自锁问题，可选用不同的单元公式：选用不同的单元公式：–不协调模式不协调模式–选择的消减积分选择的消减积分–一致的消减积分一致的消减积分–混和的混和的 U-P 公式公式•我们将在单元选择一章中讨论网格自锁与不同的单元技巧。

我们将在单元选择一章中讨论网格自锁与不同的单元技巧October 17, 200053几何非线性 – 5.7版本 网格扭曲网格扭曲网格扭曲网格扭曲•大应变－大旋转分析的每一步都会使用计算出的位移更新节点坐大应变－大旋转分析的每一步都会使用计算出的位移更新节点坐标 新新”的网格将作为下一步的起始点的网格将作为下一步的起始点•如果网格严重扭曲，这一步的求解精度将成为问题这类似于使如果网格严重扭曲，这一步的求解精度将成为问题这类似于使用扭曲的网格完成一次线性分析用扭曲的网格完成一次线性分析•扭曲的网格将迅速降低分析精度有时单元会发生扭曲的网格将迅速降低分析精度有时单元会发生“ 凹陷凹陷”October 17, 200054几何非线性 – 5.7版本 减小网格扭曲效应减小网格扭曲效应减小网格扭曲效应减小网格扭曲效应•考虑后继的网格扭曲，适当划分网格考虑后继的网格扭曲，适当划分网格•避免过度约束边界上的变形避免过度约束边界上的变形•避免使用带中间节点的单元，因为此种单元更易受网格扭曲的影避免使用带中间节点的单元，因为此种单元更易受网格扭曲的影响October 17, 200055几何非线性 – 5.7版本 在大应变分析中预测网格扭曲在大应变分析中预测网格扭曲在大应变分析中预测网格扭曲在大应变分析中预测网格扭曲•在拉伸试在拉伸试件件的的颈颈缩区，缩区，预测到预测到其后其后可能可能发生发生的的网格扭曲网格扭曲后后划分划分的的初始网格。

初始网格October 17, 200056几何非线性 – 5.7版本 在大应变分析中预测网格扭曲在大应变分析中预测网格扭曲在大应变分析中预测网格扭曲在大应变分析中预测网格扭曲未变形网格未变形网格变形网格变形网格产生大的内角产生大的内角角部单元包含更好的三角部单元包含更好的三角形形状角形形状October 17, 200057几何非线性 – 5.7版本 避免过分约束边界处的变形避免过分约束边界处的变形避免过分约束边界处的变形避免过分约束边界处的变形FF•约束边界处的所有自有度，由于约束边界处的所有自有度，由于Poisson 效应将产出非常大的效应将产出非常大的应变October 17, 200058几何非线性 – 5.7版本 避免使用带中间节点的单元避免使用带中间节点的单元避免使用带中间节点的单元避免使用带中间节点的单元•在建立大应变、大位移模型时避免使用带中间节点的单元在建立大应变、大位移模型时避免使用带中间节点的单元•带中间节点的（既高阶）单元在更新几何形状时，单元的中间节带中间节点的（既高阶）单元在更新几何形状时，单元的中间节点可能穿过单元移动，导致产生负的旋转点可能穿过单元移动，导致产生负的旋转。

October 17, 200059几何非线性 – 5.7版本 大弯曲大弯曲大弯曲大弯曲•作为大应变公式的一个子集，作为大应变公式的一个子集，ANSYS支持一些没有大应变特性的支持一些没有大应变特性的单元的大挠度，大旋转，小应变特性单元的大挠度，大旋转，小应变特性October 17, 200060几何非线性 – 5.7版本 大挠度理论大挠度理论大挠度理论大挠度理论•在在ANSYS中，大挠度理论（相对于大应变理论）假设变形与旋转中，大挠度理论（相对于大应变理论）假设变形与旋转大，但机械应变（产生应力）小每个单元除刚体运动外不改变大，但机械应变（产生应力）小每个单元除刚体运动外不改变形状•大挠度理论是大应变理论的子集，它适用于大挠度理论是大应变理论的子集，它适用于ANSYS的许多旧的梁的许多旧的梁单元、壳单元和其他非线性单元单元、壳单元和其他非线性单元October 17, 200061几何非线性 – 5.7版本 大挠度过程大挠度过程大挠度过程大挠度过程•大挠度理论的过程基于协转方法；既一个开始时平行于单元坐标大挠度理论的过程基于协转方法；既一个开始时平行于单元坐标系的坐标系系的坐标系“ 贴贴”在单元上。

此坐标系随单元一起旋转在单元上此坐标系随单元一起旋转•用这种方法，精确计算了材料（或单元）的旋转在旋转的坐标用这种方法，精确计算了材料（或单元）的旋转在旋转的坐标系中计算应变使用普通的小应变（工程应变）假设系中计算应变使用普通的小应变（工程应变）假设October 17, 200062几何非线性 – 5.7版本 大挠度位移，应力与应变的方向大挠度位移，应力与应变的方向大挠度位移，应力与应变的方向大挠度位移，应力与应变的方向•计算出的位移按初始方向输出，因为在大挠度分析中没有更新节计算出的位移按初始方向输出，因为在大挠度分析中没有更新节点坐标系取向点坐标系取向•大挠度（小应变）分析中的应力应变是旋转坐标系中的工程值；大挠度（小应变）分析中的应力应变是旋转坐标系中的工程值；既单元坐标系随单元旋转既单元坐标系随单元旋转October 17, 200063几何非线性 – 5.7版本 何时激活大弯曲效应？何时激活大弯曲效应？何时激活大弯曲效应？何时激活大弯曲效应？•当单元旋转当单元旋转“ 大大”到明显影响求解精度时，需激活大到明显影响求解精度时，需激活大挠度挠度效应•遗憾的是，没有明确的规定遗憾的是，没有明确的规定“ 大大”到底有多大。

到底有多大 小小”与与“ 大大”之间根据问题的不同相差非常大之间根据问题的不同相差非常大October 17, 200064几何非线性 – 5.7版本 求解选项求解选项求解选项求解选项•对支持大应变和对支持大应变和 / 或大或大变形变形的单元激活此功能的单元激活此功能Solution > Analysis Option ...•NLGEOM,ON 激活了支持大应变功能单元的此选项如果使用激活了支持大应变功能单元的此选项如果使用的单元只支持大的单元只支持大挠度挠度，， NLGEOM将激活大将激活大挠度挠度求解参照求解参照 ANSYS 单元手册单元手册October 17, 200065几何非线性 – 5.7版本 求解选项求解选项求解选项求解选项( (续续续续) )•激活求解控制时，打开非线性几何（激活求解控制时，打开非线性几何（ NLGEOM ）将在非线性刚）将在非线性刚度矩阵中缺省包含应力刚化项（在非线性刚度矩阵度矩阵中缺省包含应力刚化项（在非线性刚度矩阵[Knl]中包含中包含[K ] ）作为一个选项，你可对于一些旧单元选择在形成非线性）作为一个选项，你可对于一些旧单元选择在形成非线性刚度矩阵时不包含刚度矩阵时不包含[K ] 。

Solution > Analysis Options ... 很少关闭包含导致全一致很少关闭包含导致全一致[Knl]的应力刚化项的应力刚化项此命令对单元此命令对单元 106, 107, 108, 181, 182, 185, 188, and 189无作用无作用!October 17, 200066几何非线性 – 5.7版本 求解选项求解选项求解选项求解选项•推荐使用求解控制（缺省）推荐使用求解控制（缺省）•推荐的推荐的Newton-Raphson 选项是不带自适应下降的全选项是不带自适应下降的全 Newton-Raphson 选项（求解控制的缺省设置）选项（求解控制的缺省设置）•推荐使用自动时间步（求解控制的缺省设置）确定对自动时间推荐使用自动时间步（求解控制的缺省设置）确定对自动时间步设置足够小的最小时间步步设置足够小的最小时间步•线性搜索选项（线性搜索选项（ LNSRCH ）对收敛振荡问题有所帮助对收敛振荡问题有所帮助October 17, 200067几何非线性 – 5.7版本 注意事项注意事项注意事项注意事项•在大应变分析中预测网格扭曲，在大应变分析中预测网格扭曲，划分划分适当适当的的网格。

参考预测网网格参考预测网格扭曲的指南格扭曲的指南•使用适当的单元类型和积分准则以解决网格自锁问题在单使用适当的单元类型和积分准则以解决网格自锁问题在单元选择一章中有相关的更详细信息元选择一章中有相关的更详细信息•如果使用的单元不直接支持压力载荷刚度，可使用表面单元（如果使用的单元不直接支持压力载荷刚度，可使用表面单元（ SURF154 ）如果收敛速度慢，激活不对称选项）如果收敛速度慢，激活不对称选项KEYOPT(5)=1October 17, 200068几何非线性 – 5.7版本 注意事项注意事项注意事项注意事项•时间步大小应控制在每时间步大小应控制在每个个子步中的最大旋转度数小于子步中的最大旋转度数小于5或或10度•对大旋转分析不要使用预测对大旋转分析不要使用预测•梁单元和壳单元使用足够的网格密度；没有一个单元可承受超过梁单元和壳单元使用足够的网格密度；没有一个单元可承受超过30度的弯曲度的弯曲度•如果自动时间步重复进行二分，这可能是由于结构不稳定绘制如果自动时间步重复进行二分，这可能是由于结构不稳定绘制载荷位移响应曲线载荷位移响应曲线October 17, 200069几何非线性 – 5.7版本 练练练练 习习习习做：做：为熟悉大应变求解，分析一个轴对称的礅粗问题。

为熟悉大应变求解，分析一个轴对称的礅粗问题•此问题具有几何非线性（大变此问题具有几何非线性（大变形和大应变），材料非线性形和大应变），材料非线性（塑性）和状态非线性（接触）（塑性）和状态非线性（接触）特性•在在变形变形的工件中累积等效塑性的工件中累积等效塑性应变9.0 mmTarget Reaction = 800October 17, 200070几何非线性 – 5.7版本 练练练练 习习习习步骤步骤:1. 恢复数据库恢复数据库 upset.db 在此数据库中包含分析使用的有限元网在此数据库中包含分析使用的有限元网格和固定的边界条件格和固定的边界条件2. 进入求解器并在进入求解器并在pilot 节点上施加节点上施加-Y 方向的方向的9.0 mm 的位移在的位移在数据库中有一个数据库中有一个n_load 参数指明了这个参数指明了这个pilot 节点3. 激活非线性几何激活非线性几何 (NLGEOM,ON)4. 将将“ time”设置为施加的位移值设置为施加的位移值5. 监视监视pilot 节点节点(n_load) 处的反作用力，将其设置为监视文件的处的反作用力，将其设置为监视文件的变量变量1。

October 17, 200071几何非线性 – 5.7版本 练练练练 习习习习步骤步骤:6. 设置输出控制，写十步结果到设置输出控制，写十步结果到 .rst文件中文件中 (OUTRES)7. 设置自动时间步的子步数尝试设置自动时间步的子步数尝试NSUBST,100,1000,208. 求解9. 在通用后处理器中，查阅结果绘制变形形状与应力结果查阅在通用后处理器中，查阅结果绘制变形形状与应力结果查阅单元扭曲情况（单元扭曲情况（ PlotCtrls > Style > Edge Options ）可以根）可以根据单元扭曲水平发现一些潜在的问题吗？据单元扭曲水平发现一些潜在的问题吗？October 17, 200072几何非线性 – 5.7版本 练练练练 习习习习步骤步骤:10. 在时间历程后处理器中绘制力位移曲线在时间历程后处理器中绘制力位移曲线11. 查阅输出文件和监视文件注意对一定的子步有一个警告信息－查阅输出文件和监视文件注意对一定的子步有一个警告信息－塑性应变超限怎样修改子步数减少塑性应变增量？塑性应变超限怎样修改子步数减少塑性应变增量？October 17, 200073几何非线性 – 5.7版本 练练练练 习习习习额外练习额外练习•如果有时间重新开始如果有时间重新开始分析，可施加分析，可施加12.0 mm的位移。

绘制力的位移绘制力位移曲线为什么随位移曲线为什么随位移增大反作用力显位移增大反作用力显著增加？著增加？October 17, 200074几何非线性 – 5.7版本 练练练练 习习习习问题的答案：问题的答案：施加了施加了9.0 mm 的位移时，在工件的位移时，在工件“ 角部角部”的单元发生扭曲其内的单元发生扭曲其内部角接近部角接近180度第二个需要关注的第二个需要关注的是接近中心线的单是接近中心线的单元的长宽比元的长宽比October 17, 200075几何非线性 – 5.7版本 练练练练 习习习习问题的答案：问题的答案：•增加最小子步数以在每子步上施加更小的载荷增量（增加最小子步数以在每子步上施加更小的载荷增量（ 例如例如NSUBST,100,1000,40））•既然问题中包含塑性，材料特性接近不可压缩由大塑性应变引既然问题中包含塑性，材料特性接近不可压缩由大塑性应变引起的不可压缩性导致压缩工件所需的力增加起的不可压缩性导致压缩工件所需的力增加October 17, 200076几何非线性 – 5.7版本 练习解答练习解答练习解答练习解答/BATCHRESUME,UPSET,DB/SOLUNLGEOM,ONSOLC,ONTIME,9NSUBST,100,1000,20OUTRES,ALL,-10MONITOR,VAR1,N_LOAD,FYD,N_LOAD,UY,-9SOLVETIME,12D,N_LOAD,UY,-12SOLVEFINISH注意：注意：October 17, 200077几何非线性 – 5.7版本 。