铁电问题中的参量变分原理及其应用.pdf

上海交通失擘博士后研究报告 摘要 对于铁电材料的裂尖场分析，线性压电理论是不充分的，必须考虑非线性力电耦台行为 才可能获得可信的理论和数值结果。因此，近年来铁电的宏、细观本构模型的研究有了积极 的发展。其中，一类明确借用宏观弹塑性理论中的概念和方法的本构模型对有限元分析具有 特别的吸引力。然而，将传统弹塑性有限元分析推广到此类铁电模型并不是一个显然的过程。 因缺失了传统弹塑性中能量二次型的正定性，当前普遍应用的基于电焓驻值原理和迭代法的 铁电数值算法的稳定性是难以保证的。因此，可靠的铁电数值算法依然是～个有待探讨的问 题，无论是有限元列式或非线性算法都需要深入研究。 在借用弹塑性理论基本概念的铁电模型中，通过引入一个类似屈服准则的畴变 (switching)准则，将由一个“塑性因子”^来控制材料状态在“弹性”和“塑性”之间 的变化。 fA0”plasticity“--switchi”g l^；0”elasticity”-十nDswi缸hmg 对于复杂结构的分析，材料的“弹／塑性”状态不可能预先指定，因此一般需要迭代法来计 算。现有文献中基本都采用了这一方法。然而，在弹塑性问题的计算中．已经提出了所谓的 参变量变分原理和参数二次规划算法。参变量变分原理将弹塑性问题处理为以塑性因子为控 制变量的最优控制问题，因变形能的正定性，可以应用最优控制的二次规划算法来计算弹塑 性问题。对同一个半隐式本构积分格式，该算法无须传统方法中必不可少的迭代计算，因此 数值稳定性较传统算法具有一定的优势。因Hehnholtz自由能的正定性。也可以将铁电问题 处理为以^为控制变量的最优控制问题，并以二次规划算法计算之，从而获得良好的数值稳 定性。 本文对上述～类铁电模型提出了一个参变量变分原理和参数二次规划算法，并相廊地应 的有限元具有正定性，因此在数值稳定性上具有明确的优势，对于算法稳定性分析也是有利 的。 在第一章中，建议了一个称为参变量最小Helmholtz自由能原理的变分原理。该原理适 肖j于前述的铁电宏观本构模型，以及铁电的品体塑性细观本构模型。在第一二、三章中分别计 算了单轴和多轴宏观模型的典型电滞回线、蝶形回线和铁弹曲线。计算单轴模型是为了对本 文原理和算法有一个清楚的示例，而计算多轴模型则是对算法有一个充分的考察。另外，在 第三章中还应用本文方法计算并讨论了modeE裂尖场的非线性行为．既是应用，也是对方 法的深入考察。 本报告其余部分是关于应用应力函数构造有限元的一些研究。 借助于Cossemt连续介质模型，在第四章系统地探讨了应力函数和位移对避免有限元 Cl连续性困难的互补性作用。应用几何学的外微分工具构造应力函数，显示了应力函数与 位移在几何上的相似性。通过对应力函数对偶理论的深入分析，为将应力函数列式得笠I}的余 能单元转化为具有一般位移自由度的势能单元提供了严格的理论基础，在此基础上，给出应 川廊力函数构造有限元的一般方法。 偶应力理论是有限元列式面临C1连续性困难的问题之一。平面偶应力理论的应力函数 元的同一性，而MM板位移有限元并不存在C1连续性困难。因此，在第五章研究了应川施 力函数列式构造具有一般位移自由度的平面偶应力单元的一般方法。根据这一方法，将典型 单元．数值结果表明该单元具有良好的精度和收敛性。 上海交通大学博士后研究报告 传统多变量变分原理中并没有考虑以应力函数为变分宗量。作为对应用应力函数的有限 元研究的一个深入，在第六章中研究了含应力函数和应变二类变量的Kirchhoff板弯曲有限 交列式可以构造一个含应力函数的类H．R变分原理的Kirchhoff板弯曲单元，其收敛性利计 算精度都是令人满意的。 关键词：铁电：参变量变分原理：有限元：二次规划：C‘连续性：Cosserat介质：应力函 数：对偶：比拟：偶应力：R／M板：Kirchhoff板，类H—R变分原理 上海交通大学博士后研究报告 ABSTRACT Duetononlinearelectromechanicalnearcrack isnot response tip，linearpiezoelectri