《选择非线性分析》PPT课件.ppt

第七章,单元选择非线性分析,October18,2000,7-2,单元选取5.7版本,单元算法,传统位移协调方法Solid45KEYOPT(1)=1由于剪切锁定而很少使用非协调模式(额外形函数)Solid45缺省选项，弯曲变形选择缩减积分(B-Bar)几乎不可压缩材料，体积变形一致缩减积分(URI)几乎不可压缩材料，弯曲变形混合U-P公式不可压缩材料，超弹性,October18,2000,7-3,单元选取5.7版本,单元算法(续),为何有如此多的不同单元公式?普通非线性求解非常费时，采用不同的单元技术可更加有效地解决各种类型的非线性问题。不同材料行为(弹性、塑性、超弹性)和不同的结构行为(体积变形、弯曲)需要选择不同的单元算法。ANSYS单元库为不同类型非线性问题提供了选择适当单元的“工具条”。,October18,2000,7-4,单元选取5.7版本,传统位移协调方法,无额外形函数的完全积分低阶单元(Solid45KEYOPT(1)=1)和完全积分高阶单元(Plane2)为传统位移公式的实例。对于传统位移协调方法，数值积分对于应变能的所有组份都是数值精确的，只有节点位移(UX,UY,UZ)是基本变量。注意：由于剪切锁定和体积锁定，此公式的低阶单元极少使用。,October18,2000,7-5,单元选取5.7版本,剪切锁定,在弯曲问题中，完全积分低阶单元呈现“过分刚硬”。在弯曲中这种公式包括实际上并不存在的剪切应变，称为寄生剪切。(从纯弯曲中的梁理论可知剪切应变xy=0),微体积纯弯曲变形中，平直断面保持平直，上下两边变成圆弧，xy=0。,完全积分低阶单元变形中，上下两边保持直线，不再保持直角，xy不等于0。,October18,2000,7-6,单元选取5.7版本,剪切锁定实例,P,DeflectionofaCantileverBeam(depth/length=1/100),当长宽比增加时，在弯曲时完全积分低阶单元将发生锁定。如Plane82的二阶单元不会发生剪切锁定问题(二次形函数允许边部弯曲)。网格细分一般无助于剪切锁定!,October18,2000,7-7,单元选取5.7版本,不用于弯曲的单元公式,*在几何非线性分析中不要使用特殊形状。,October18,2000,7-8,单元选取5.7版本,当材料特性是几乎不可压缩的并使用完全积分单元时发生体积锁定(泊松比约为0.5)，不可压缩性能存在于超弹性材料或塑性流动中。在单元中出现伪压应力，从而使单元对于不会引起任何体积变化的变形“过于刚硬”。体积锁定可出现于不同的应力状态中，包括平面应变、轴对称和3-D应力。对于平面应力不会发生体积锁定。,体积锁定,October18,2000,7-9,单元选取5.7版本,体积锁定实例,厚壁轴对称圆筒的径向位移,October18,2000,7-10,单元选取5.7版本,体积锁定,体积锁定可通过压应力“棋盘状”模式(相邻单元间变化显著)检测出。在超弹性模型中可用单元等值线绘图(PLESOL)绘制静水压力(HPRES)等值线来验证此行为。如怀疑存在体积锁定，可试细分高静水压力区域的网格或改变单元类型。,October18,2000,7-11,单元选取5.7版本,非协调模式,低阶完全积分单元的形函数可被表示常曲率状态的模式所增强，这些增加的模式作为内在的自由度，因其导致网格的缝隙和重叠而被称为非协调模式。,非协调模式,无非协调模式,F,2F,F,F,2F,F,F,2F,F,F,2F,F,October18,2000,7-12,单元选取5.7版本,非协调模式,一般地非协调模式可得到精确的结果，单元间的缝隙和重叠可软化结构，阻止传统位移公式中的“过分刚硬”现象。注意，当在退化形状中生成非协调模式的单元时，非协调模式的单元将不再具有其特殊形状。例如Plane42为四边形时有非协调模式,但当退化为三节点三角形时Plane42成为常应变三角形。*警告*区域2的网格包含PLANE42三角形，这些三角形单元在弯曲中过于刚硬，如有可能采用二次(6或8节点)单元。,October18,2000,7-13,单元选取5.7版本,非协调模式,悬臂梁的挠曲(厚度/长度=1/100),非协调模式以附加公式为代价增加精度：二维:4DOF/单元三维:9DOF/单元增大.esav文件大小基于非协调模式的单元公式：5,41,42,45,43,63,143,181.为弯曲设计!,October18,2000,7-14,单元选取5.7版本,选择缩减积分(B-Bar),选择缩减积分(B-Bar)在几乎不可压缩材料中用单元平均体积应变代替积分点应变从而消除体积锁定。B-Bar把B矩阵分成两部分，即体积应变(主应变)和偏差应变。,B-Bar的优势在于其不增加求解中的自由度数。推荐将B-Bar公式用于体积变形塑性应用中，因对于B-Bar仍存在剪切锁定问题(弯曲)。,October18,2000,7-15,单元选取5.7版本,B-Bar实例,厚壁轴对称圆筒的径向位移,October18,2000,7-16,单元选取5.7版本,一致缩减积分,一致缩减积分(URI)采用比精确数值积分所需阶数低的低阶积分规则进行积分。这将导致单元变形更加容易，有助于消除剪切锁定和体积锁定。URI不需要附加的自由度。文件减小，单元计算所需CPU时间缩短(尤其对于材料非线性)。然而URI会引起应变能为零的变形模式，这被称为零能量或沙漏模式。,October18,2000,7-17,单元选取5.7版本,沙漏模式,沙漏模式是由于变形而引起零应变能的变形模式。如右图所示两例，在只有一个积分点的低阶单元中，此单个积分点未获得任何单元应变能。这可导致出现不切实际的行为。,October18,2000,7-18,单元选取5.7版本,低阶单元以及粘塑性实体单元Visco106、Visco107和Visco108。,October18,2000,7-24,单元选取5.7版本,推荐使用的实体单元,线性分析对于线性分析采用有非协调模式的一阶单元(非退化形状的Plane42,Solid45)，二阶单元适用于高应力梯度和应力集中区域中。采用Solid92进行高性能求解。具有缩减积分的Solid95单元适用于要求泊松比接近0.5的大部分问题。,October18,2000,7-25,单元选取5.7版本,推荐使用的实体单元(续),几乎不可压缩材料(塑性)可忽略弯曲的体结构变形采用Plane182,Solid185选择缩减积分(B-Bar)。对于小应变应用采用非协调模式单元Plane42,Solid45。对于大应变应用采用具有URI(特别对于大模型)的Plane182和Solid185或具有URI的Solid95。也可以采用Visco106,Visco107和Visco108单元(甚至对于与速率无关的塑性)。,October18,2000,7-26,单元选取5.7版本,推荐使用的实体单元(续),不可压缩材料(超弹性)类似橡胶材料的不可压缩条件约束需要U-P混合公式，可采用单元Hyper56,Hyper58,Hyper158或Hyper74。,October18,2000,7-27,单元选取5.7版本,壳单元-概述,当结构的总体厚度相对于典型长度很小时可使用壳单元，长度比厚度大十倍以上的问题可决定使用壳单元。公开的文献中有各种壳理论，这源于对壳位移的不同近似描述。ANSYS中的壳单元根据要求解的问题类型采用不同的公式，三个基本的壳公式包括:薄膜理论,“薄”壳理论和“厚”壳理论。,October18,2000,7-28,单元选取5.7版本,壳单元-概述(续),薄膜理论Shell41采用薄膜理论。Shell41忽略弯曲和横向剪切，只包含薄膜效应。经典Love-Kirchhoff理论Shell63是“薄”壳单元。Shell63包含弯曲和薄膜效应但忽略横向剪切变形。Reissner/Mindlin理论Shell43,143,181,91,93和99是“厚”壳单元。其包含弯曲、薄膜和横向剪切效应。横向剪切被表示为整个厚度上的常剪切应变,这种一阶近似只适用于“中等厚度”壳体。,October18,2000,7-29,单元选取5.7版本,平面变形中的壳单元,平面内壳的响应可认为是平面应力状态，因此对于壳单元不会出现体积锁定问题。(当绝对不可压缩，泊松比=0.5时Shell181支持超弹性)对于薄膜现象，壳单元的平面公式与平面实体单元的公式相似(非协调模式)。Shell41,43,63和181对于平面内变形支持非协调模式。Shell181也支持具有沙漏控制的一致缩减积分(缺省选项)。,October18,2000,7-30,单元选取5.7版本,推荐使用的壳单元,线性分析如壳的厚度非常小采用Shell63，Shell63单元不包含横向剪切效应。如横向剪切变形重要，对于均匀材料采用Shell43,Shell93或Shell143，对于复合材料采用Shell91或Shell99。注意具有一致缩减积分(缺省)的单元Shell181对大模型较快，但将需要较细的网格。,October18,2000,7-31,单元选取5.7版本,推荐使用的壳单元(续),非线性分析等向强化塑性和超弹性采用Shell181。其优势包括：较小的.esav文件，较少的CPU时间，压力载荷刚度效果，可以导入初始应力，接触分析中厚度变化。运动强化塑性，蠕变采用Shell143,Shell43和Shell93。Shell43和Shell143适用于小应变塑性，Shell93是弯曲的壳(高阶)。,October18,2000,7-32,单元选取5.7版本,梁单元-概述,梁单元可用于分析主要受侧向或横向载荷的结构，典型的梁应用包括：机器主轴，房屋构架，桥梁等。ANSYS中可用的两个梁单元公式为：Euler/Bernoulli梁Beam3和Beam4包括弯曲、轴向和扭转变形。横向剪切变形不包括于单元公式中(但可作为柔性因子应用)。Timoshenko梁Beam188和Beam189在单元公式中包括弯曲、轴向、扭转和横向剪切变形。,October18,2000,7-33,单元选取5.7版本,推荐使用的梁单元,线性分析对于线性模型采用Beam3,Beam4,Beam188或Beam189。Beam3和Beam4采用Hermitian多项式作为形函数并且在弯曲中具有三次响应。Beam188采用线性多项式作为形函数，Beam189采用二次多项式作为形函数(可作为弯曲梁)。注意有限应变梁需要更细化的网格，然而，Beam188和Beam189具有许多优秀的前后处理特色，使其非常适于处理线性模型。,October18,2000,7-34,单元选取5.7版本,推荐使用的梁单元(续),非线性分析采用Beam188和Beam189模拟各向同性强化塑性、大应变、屈曲(特征值和非线性坍塌)和/或大转动问题。用作壳加强单元，Beam188与Shell181完全兼容，并且Beam189与Shell93完全兼容。注意所有梁单元假定无约束的断面扭曲，当扭曲约束是重要的时，采用壳单元来模拟梁结构。,October18,2000,7-35,单元选取5.7版本,Beam188和Beam189,由于Beam188和Beam189对剪切变形采用一次近似(厚度方向为常剪切应力)，故其只应用于比较细长的梁。可用如下方法近似估计，(GAL2)/(EI)30,其中L是所有单元的总体长度(不是单个单元的长度)。,因弯曲响应是线性或二次的，故长度方向采用多个单元,L,