ansys高级非线性分析二单元技术

单元技术二AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-2单元技术本章概述•本章主要讨论18x系列的单元。超越过去的ANSYS版本,18x单元已成为非线性应用中选择的单元。–18x单元包括强大的单元公式和大量的本构模型库。•对18x单元，材料和单元技术已经分开。这就提供了一个更小的单元库，可作为一个“工具箱”，用于处理不同情况和各种本构模型。–SHELL181和BEAM188/189还具有高级的前后处理工具，这些工具是梁和壳单元特有的。AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-3单元技术...本章概述•这章的要点是:–完全积分的、传统的基于位移的连续单元在一定情形下低估位移•这称为网格锁定–因此，有不同的单元公式来处理这些问题,基于:•体积或弯曲占优的问题(结构行为)•弹性，塑性或超弹性(材料行为)•非线性求解的效率–除连续单元外，ANSYS还有庞大的壳和梁单元库•单元选择主要考虑的是基于‘薄’或‘中等厚’的壳/梁AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-4单元技术...本章概述•这章包括下列主题:A.传统的基于位移的连续单元B.连续单元中剪切和体积锁定C.选择缩减积分(B-bar)D.一致缩减积分(URI)E.增强应变公式F.混合U-P公式G.对连续单元的一般建议H.壳单元I.梁单元AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-5单元技术A.传统位移公式•无附加自由度的完全积分的低阶和高阶单元是传统的基于位移的单元的例子。–SOLID45(KEYOPT(1)=1)和PLANE42(KEYOPT(2)=1)是低阶完全积分的传统位移公式的例子。–SOLID95(KEYOPT(11)=0)是高阶完全积分传统位移公式的例子。•这实际上是14点积分公式而不是3x3x3积分方案,以后会讨论。14点积分公式比完全积分方案更有效。AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-6单元技术...传统位移公式•回顾积分点的一些重要细节:–对任何单元,自由度解{Du}是在结点求出–在积分点计算应力和应变。它们由自由度导出。例如可以由位移通过下式确定应变:[B]称为应变-位移矩阵–后处理结果时，积分点应力/应变值外推或拷贝到结点位置•右图所示为2x2积分的四结点四边形单元，红色为积分点。s,euuBDDeAdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-7单元技术...传统位移公式•传统的基于位移单元的积分点遵循Gauss积分法且和单元的阶数相同。这称为完全积分。•换句话说,完全积分意味着数值积分方法对未发生几何扭曲单元的应变能的所有分量是精确的。ElementTypeFullIntegrationOrder4NodeQuad2x28NodeQuad3x38NodeHex2x2x220NodeHex3x3x3AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-8单元技术...传统位移公式•完全积分、低阶传统位移单元易于发生剪切和体积锁定，因此很少使用。•完全积分、高阶传统位移单元也易于发生体积锁定。AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-9单元技术B.剪切和体积锁定•传统的基于位移的单元有两个问题:剪切锁定和体积锁定:–剪切锁定导致弯曲行为过分刚化(寄生剪切应力)。当细的构件承受弯曲时，这是一种几何特性。–体积锁定导致过度刚化响应。当泊松比接近或等于0.5时，这是一种材料特性。•本章重点讨论用不同单元公式解决这两个问题的方法。主要讨论连续(实体)单元。•由于非线性分析花费计算机时间太多，所以有些单元公式也提供了更有效地解决非线性问题的方法。AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-10单元技术...剪切锁定•在弯曲问题中完全积分的低阶单元呈现“过度刚化”。这个公式包含了实际并不存在的剪切应变，称为寄生剪切。(从纯弯曲的梁理论可知剪切应变gxy=0.)MMMMxy微体积纯弯曲变形中，平面截面保持平面，上下两边变成圆弧,gxy=0。完全积分的低阶单元变形中，上下两边保持直线，不再保持直角，gxy不为零。AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-11单元技术...剪切锁定实例•当长厚比增加时，模型更容易剪切锁定.–因为寄生的剪切应变/应力，所以产生的位移被低估。–下面的例子是弯曲中的梁。这种情况下剪切应力接近于零，但是如SXY等高线图中所示，发生了剪切锁定。AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-12单元技术...剪切锁定实例•这个模型呈现剪切锁定了吗?单元182(B-Bar),几乎不可压缩的Mooney-Rivlin超弹材料的平面应变答案:很意外,没有。该模型具有超弹材料属性,以B-Bar和增强应变运行,结果f非常相似.AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-13单元技术...体积锁定•材料行为是几乎或完全不可压缩时(泊松比接近或等于0.5)，在完全积分单元中发生体积锁定。–超弹材料或塑性流动可发生不可压缩(以后讨论)。–单元中产生的伪压应力导致单元对不会引起任何体积变化的变形“过度刚化”。–体积锁定也会引起收敛问题。•各种应力状态都会发生体积锁定，包括平面应变、轴对称及3-D应力。–对平面应力问题不会发生体积锁定，因为平面外应变用于满足体积不可压缩条件。AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-14单元技术...体积锁定•可把应力分解为静水压力(p)和偏差应力(s)分量:•静水压力(p)定义为体积模量(k)和体积应变(ev)的乘积:ijijijspszyxzyxvzyxvEEpssseeeeksssek2121331AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-15单元技术...体积锁定•前面幻灯片中的公式中，若泊松比接近或等于0.5,可看出:–体积模量k将很大或无穷大–体积应变ev将接近或等于零–这被称为几乎或完全不可压缩材料行为•几乎或完全不可压缩材料存在数值上的困难，且呈现出过度刚化行为。–这在体积变形问题中显而易见–从计算观点来看，对几乎不可压缩和完全不可压缩问题的处理不同。•体积锁定导致静水压力(p)的交变模式(棋盘状)，存在非线性材料时对单元可用NL,HPRES后处理静水压力。AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-16单元技术...体积锁定实例NL,HPRES的等值图如右图所示。只要有非线性材料就可得到这种输出量。用单元求解(PLESOL)后处理静水压力(NL,HPRES)使用户可以验证体积锁定是否是个问题。AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-17单元技术...练习请参考附加练习题:•练习1:剪切锁定连续单元第二章C-G节AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-19单元技术连续单元公式•后面将讨论一般的指南和建议。而下面的各部分将详细介绍用以克服剪切和体积锁定的单元技术。C.选择缩减积分(B-bar)D.一致缩减积分(URI)E.增强应变公式F.混合U-P公式AdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-20单元技术...连续单元公式•作为一个简单的解释，剪切和体积锁定是由于系统的过度约束。•利用不同的单元公式通过放松约束或引入附加的方程求解这些约束来解决这个问题。•不幸地是,没有现成的单元公式能最有效地解决锁定问题.因此在下面部分将从正反两方面来讨论每个公式。RelaxConstraintsAdditionalEquationsElementFormulation(ReducedIntegration)(AddExtraDOF)B-BarxURIxESxMixedU-PxAdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-21单元技术...连续单元公式目前在18x单元中有四个不同的单元技术:B-Bar,URI,增强应变和混合U-P。它们用于处理剪切和体积锁定:–高阶18x单元(PLANE183,SOLID186-187)通常用URI。–缺省时低阶18x单元(PLANE182,SOLID185)用B-Bar。–B-Bar和增强应变不能用于高阶单元。–混合U-P技术独立于其它技术,所以可以和B-Bar,增强应变或URI联合使用。ElementTechnologyLower-OrderElementsHigher-OrderElementsShearLocking(Bending)Nearly-Incompressible(Plasticity,Hyperelasticity)Fully-Incompressible(Hyperelasticity)B-BarY-NYNEnhancedStrainY-YYNURIYYYYNMixedU-PYYNYYAdvancedStructuralNonlinearities6.0TrainingManualSeptember30,2001Inventory#0014912-22单元技术...连续单元公式•单元选项允许用户选择合适的单元公式。–MainMenuPreprocessorElementTypeAdd/Edit/Delete…“Options”buttonindialogbox–若用命令,•KEYOPT(1)用于PLANE182的B-bar,URI和增强应变•KEYOPT(2)用于SOLID185的B-bar,URI和增强应变•K