Cosserat材料弹塑性摩擦接触分析及多尺度有限元方法.doc

Cosserat材料弹塑性摩擦接触分析及多尺度有限元方法Elastic-Plastic Frictional Contact Analysis and Multiscale Finite Element Method for Cosserat Materials大连理工大学博士学位论文格式规范- PAGE 122 -摘 要材料的微结构对于材料的宏观力学性能存在着巨大的影响。在对一些特定材料（如：颗粒材料、多孔介质和金属基复合材料等）和物理现象（如：应力集中和应变局部化等）进行数值模拟时，经典连续介质力学理论在物理及数值方面存在着一定的缺陷，其主要原因是由于控制方程内尺度参数的缺失以及局部数学特性的改变导致了相应的数值有限元解存在着网格相关性，有时甚至不能得到正确的数值结果。为了更有效的反映材料的微结构特性，研究者提出了一些广义的连续介质本构模型，例如Cosserat连续体模型。Cosserat连续体模型中引入了独立的旋转自由度，并通过在本构关系中引入材料内尺度参数。从力学观点出发，许多微孔材料（例如：金属泡沫、钢架结构、颗粒材料和海绵骨）都可以看作是内部由微梁结构组成的材料。当这类材料受到外载作用时，材料内部的弯曲作用是主要的变形机制，也就是说位移与转动同时存在。由于Cosserat连续体模型中旋转自由度以及新的材料参数的引入，使得Cosserat连续介质模型在对微孔材料进行力学分析时成为了一个合适的选择。本文基于参数变分原理，发展了Cosserat模型三维弹塑性分析的参数二次规划算法，并对工程中常见的应变局部化问题进行了模拟；构建了Cosserat材料弹塑性摩擦接触问题的参数二次规划算法，对Cosserat材料参数对接触面力学特性进行了分析；另外，提出了非均质Cosserat材料的多尺度有限元方法，避免了一般均匀化方法确定Cosserat材料参数的困难，并获得了具有较高精度的数值结果。首先，本文基于参数变分原理及其相应的参数二次规划算法，发展了Cosserat模型三维弹塑性分析的有限元模型，并采用该模型对工程中常见的应变/软化局部化问题进行了分析。由于Cosserat连续介质模型的本构方程中存在材料的内尺度参数，该模型可有效的避免经典连续介质模型进行数值模拟是遇到的网格依赖性问题。数值结果表明所发展的三维Cosserat连续介质弹塑性有限元模型可以有效的模拟应变局部化现象并且具有很好的数值稳定性，同时所获得的数值结果具有良好的非网格依赖性。其次，给出了Cosserat材料弹塑性接触问题的参数变分原理，以此为基础发展了Cosserat材料平面及空间弹塑性摩擦接触问题的参数二次规划算法。弹塑性问题和摩擦接触问题统一采用参数变分原理进行处理，大大简化了原双非线性问题的求解过程，同时采用一定的数学方法消除了接触问题中引入的惩罚因子，保证了算法的稳定性并提高了数值结果的计算精度。数值结果说明了所发展模型的重要性以及可行性，且在某些状况下接触面接触力的分布情况强烈依赖于Cosseart材料参数。最终，本文提出了周期性非均质Cosserat材料扩展多尺度有限元分析模型。其主要特点采用数值方法构造多尺度基函数有效的扑捉粗网格的细尺度特性。考虑到Cosserat模型中转角自由度的存在，基于两节点梁单元的转角与横向位移关系，构造了扩展多尺度有限元方法（EMsFEM）的指定边界条件，并发展了相应的周期性边界条件。采用上述两种边界条件，分别构造了Cosserat模型的位移与转角场的数值基函数，建立了宏观变形与微观应力-应变间的关系。该方法不需要像一般均匀化方法那样估计非均质Cosserat材料总体参数。数值算例说明了所发展算法的有效性与高效性。关键词：参数变分原理；Cosserat模型；应变局部化；接触问题；多尺度有限元法Elastic-Plastic Frictional Contact Analysis and Multiscale Finite Element Method for Cosserat MaterialsAbstractMicrostructures of materials play an important role in determining their macroscopic mechanical behaviors. It has been found that classical continuum theory is deficient both on physical and numerical aspects for some special materials (such as granular materials, porous media and metal matrix com