随着当代科技的迅猛发展，大批性能优越的智能复合材料不断涌现。目前用于主动控制阻尼器和谐振器的智能复合材料主要为电磁流变类材料，其材料性能满足主动控制阻尼器和谐振器中对的刚度和阻尼改变、较快响应速率等的特殊要求。

作为磁流变材料中的一员，磁流变弹性体(Magnetorheological Elastomers, MREs)与更早出现的磁流变液材料相比，不仅继承了其可控性、可逆性、响应迅速等高技术特征,还具有稳定性好、结构设计简单、制备成本低等独特的优点,因而拥有巨大的应用潜力。由于基体聚合物材料的高延展性，磁流变弹性体通常在较大荷载（外力与外磁场）作用下会发生大变形行为。其中，无磁场下磁流变弹性体的大变形力学响应作为其基础性力学行为，需要优先开展深入研究。磁流变弹性体类属的聚合物基复合材料大变形行为分析是当今学术界研究的热点方向。

针对聚合物基复合材料非线性大变形行为分析，现有的常用理论分析方法有细观力学方法和分子链网络模型两大类。鉴于复合材料的有效行为通常依赖于其局部微观结构和性能的平均响应，并不需要建立复杂的分子链网络模型，因此，细观方法在分析这类问题时显得格外有效。目前国际上已经有许多学者针对复合材料在大变形行为开展了相关的理论研究， Bergstrom和Boyce ^[1]用实验以及有限元法研究了硬质颗粒对橡胶基复合材料宏观力学性能的影响。关于复合材料的有限变形理论研究的实质性进展，直到Ponte Castaneda^[2][3]开发出二阶均化方法后才得以有效解决，并采用优化计算方法来解决能量最小问题，从而较为准确的确定不同阶段间的变形相容条件；.Bouchart^[4][5]等人提出了研究可压缩超弹性复合材料的另一种等效分析方法。 Omnès^[6]等人开发了一个广义的自洽方法，构建三相模型预测分析了填充橡胶的宏观力学性能。Yin ^[7]基于Eshelby^[8]张量建立了颗粒增强复合材料的本构模型。 Nemat-Nasser ^[9][10]首先推广了细观力学，将修正的Eshelby张量应用于有限形变分析。Abdi^[11]基于均匀化理论，预测了非均质材料在大变形时的剪切响应。 Wang等人^[12]考虑到聚合物基体的聚集程度和超弹性力学性质的随机分布，使用2D微观体积单元作模拟分析复合材料中的变形和应变场。由此看来，若将用于分析材料线性变形的细观力学方法用于非线性材料大变形分析，这就要求对细观力学模型进行修正，提出更先进的理论方法，进而较为准确的预测磁流变弹性体的大变形力学响应。L. W. Zhang^[13]对载荷均匀分布的纳米复合材料进行了几何非线性分析,采用弧长迭代算法求解纳米复合板的非线性响应。Tuan Le-Manh; Tuan Luu-Anh^[14]采用等阶分析与一阶剪切变形理论相结合的方法导出有限元模型。 得到了冯-卡门理论意义上的应力—应变关系。李剑锋和龚兴龙^[15]对磁流变弹性体的动态阻尼进行了研究和分析，分析了不同测试条件下的应变幅值.介绍了磁流变材料的数值强度分析。G S#322;awiński ， T Niezgoda ^[16]研究了羰基铁颗粒的体积分数对材料的力学性能的影响，特别是对杨氏模量的影响。

本文将基于细观力学中的自洽法方法，运用双夹杂（DI）模型，通过基体的大变形行为分段线性化处理的方式，分析各向同性磁流变弹性体的大变形非线性力学行为。本研究的开展将会丰富磁流变弹性体基础行为的分析理论，为工程应用中磁流变弹性大变形行为予以理论指导。

2. 研究的基本内容与方案

研究的基本内容：鉴于各向同性磁流变弹性体是一种颗粒随机分布颗粒增强复合材料，本文运用细观力学理论，结合分析非线性行为的数值计算方法，分析各向同性磁流变弹性体的大变形非线性力学行为。

目标；修正细观力学方法的自洽模型，定量分析各向同性磁流变弹性体的大变形非线性力学行为。

拟采用的技术方案：

鉴于细观力学方法仅适用于线性变形准静态加载方式下的力学行为，本文运用变形阶段细化的方式，将大变形行为过程分解为若干线性小变形阶段进行分析，针对每个阶段使用细观力学模型分析，使用自洽方法，分析磁流变弹性体的力学响应，并考虑该响应对对下一个阶段的材料基体性能的影响，最后将分为若干个小变形阶段进行整合，获得各向同性磁流变弹性体大变形过程中的力学响应曲线。