1. School of Materials Science and Engineering, Harbin Institute of Technology, Harbin, China.

*             Correspondence:yfzhou@hit.edu.cn

Keywords: model, effective, magnetic properties, magnetorheological elastomers, dipoles

1. INTRODUCTION

The field of smart materials has witnessed explosive growth over the past few decades, driven by the demand for materials that can sense and respond to their environment in a controlled and predictable manner. Among these, magnetorheological (MR) materials stand out for their ability to undergo rapid and significant changes in their mechanical properties upon the application of a magnetic field [1]. While magnetorheological fluids (MRFs) have been known and utilized for some time, their solid-state counterparts, magnetorheological elastomers (MREs), have emerged as a highly promising class of smart composites. MREs are essentially elastomers, such as silicone rubber or polyurethane, filled with micron-sized magnetizable particles, most commonly carbonyl iron particles [2]. In the absence of a magnetic field, these materials behave like conventional elastomers. However, when a magnetic field is applied, the embedded particles become magnetized and interact with each other, leading to a significant increase in the material’s stiffness and damping properties. This change, known as the magnetorheological effect, is rapid, reversible, and can be precisely controlled by tuning the strength of the applied magnetic field [3]. The unique properties of MREs make them ideal candidates for a wide range of applications, including:

• Vibration and shock absorption: MRE-based dampers and isolators can adapt their stiffness in real-time to effectively suppress vibrations over a broad frequency range.
• Sensors and actuators: The magneto-mechanical coupling in MREs can be exploited to create novel sensors that detect magnetic fields or mechanical deformations, as well as actuators that provide controlled motion.
• Haptic devices: MREs can provide tunable tactile feedback in virtual reality and telerobotics applications.
• Soft robotics: The ability of MREs to change shape and stiffness on demand makes them a key enabling technology for the development of soft, compliant robots [4].

The performance of any MRE-based device is fundamentally linked to the material’s magnetic properties. A thorough understanding of how the magnetic particles interact with each other and with the surrounding matrix under the influence of an external field is paramount for designing MREs with tailored functionalities [5]. This has led to the development of various theoretical models aimed at predicting the effective magnetic properties of these complex composite materials.

2. THE CLASSICAL DIPOLE MODEL: A NECESSARY FIRST STEP BUT AN INCOMPLETE PICTURE

For many years, the classical dipole model has been the go-to theoretical framework for understanding the magnetic behavior of MREs. This model simplifies the complex reality of MREs by making a series of key assumptions:

• Point Dipoles: The magnetic particles are treated as ideal point dipoles located at the center of each particle. This approximation is valid only when the distance between particles is significantly larger than their diameter.
• Uniform Magnetization: The magnetization within each particle is assumed to be uniform.
• Negligible Multipolar Interactions: The model considers only the dipole-dipole interactions between particles and neglects higher-order multipolar interactions, which can become significant at high particle concentrations [6].
• Regular Particle Arrangement: The particles are often assumed to be arranged in a regular lattice, such as a simple cubic or body-centered cubic lattice. This does not accurately represent the random or field-induced anisotropic microstructures found in real MREs.

The dipole model provides a straightforward way to calculate the magnetic interaction energy between particles and, consequently, the magnetic field-induced changes in the MRE’s properties [7]. While it has been instrumental in providing a qualitative understanding of the MR effect and has been successful in predicting the behavior of MREs with low particle concentrations, its limitations become glaringly obvious as we move towards more technologically relevant materials with higher particle loadings. The primary shortcomings of the dipole model include:

• Underestimation of Magnetic Forces: At high particle concentrations, where particles are in close proximity, the dipole model significantly underestimates the magnetic interaction forces. This is because the non-uniform magnetization and higher-order multipolar interactions, which the model ignores, play a much larger role.
• Inability to Capture Microstructure Effects: The model’s assumption of a regular particle arrangement fails to account for the complex microstructures that arise during the curing process of MREs, especially when cured in the presence of a magnetic field (anisotropic MREs). The formation of particle chains and clusters dramatically alters the magnetic response of the material [8].
• Neglect of Magneto-Mechanical Coupling: The dipole model is purely a magnetic model and does not inherently account for the mechanical properties of the elastomeric matrix. The deformation of the matrix under the influence of magnetic forces can, in turn, affect the particle arrangement and the overall magnetic response [9].

These limitations have spurred the development of more sophisticated models that aim to provide a more accurate and comprehensive description of the effective magnetic properties of MREs.

3. TOWARDS A MORE REALISTIC DESCRIPTION: BEYOND THE DIPOLE APPROXIMATION

To overcome the shortcomings of the classical dipole model, researchers have explored several avenues, incorporating more realistic physical phenomena into their theoretical frameworks. Continuum mechanics models treat the MRE as a continuous medium with effective properties that are derived by averaging the properties of the constituent phases (the particles and the matrix). These models can account for the magneto-mechanical coupling and can be used to predict the macroscopic behavior of MREs under combined mechanical and magnetic loading. However, they often rely on phenomenological parameters that need to be determined experimentally and may not provide a deep insight into the microscopic origins of the MR effect [10]. Finite element analysis has emerged as a powerful tool for simulating the behavior of MREs. By meshing the MRE into a large number of small elements and solving the governing equations of magnetism and mechanics for each element, FEA can provide a detailed picture of the local magnetic fields, stresses, and strains within the material [11]. FEA can readily handle complex microstructures and non-linear material behavior, making it a valuable tool for understanding the influence of particle arrangement and concentration on the effective magnetic properties. While FEA is highly accurate, it is also computationally expensive, especially for large-scale simulations with a high number of particles. This can make it impractical for routine design and optimization tasks. Micromechanical models bridge the gap between the overly simplified dipole model and the computationally intensive FEA. These models start from the microscopic level, considering the interactions between individual particles, but then use statistical methods to upscale these interactions to predict the macroscopic properties of the material [12]. They can account for factors such as particle size distribution, particle shape, and the formation of particle clusters.

4. A NEW PROPOSED MODEL: INTEGRATING MULTIPHYSICS FOR A COMPREHENSIVE VIEW

Building upon the strengths of these previous approaches, we propose a new, more comprehensive model for describing the effective magnetic properties of MREs. This new model is based on a Multiphysics framework that integrates key aspects of micromagnetic, continuum mechanics, and statistical physics. The core tenets of this new model are:

• Accurate Representation of Particle Magnetization: Instead of treating particles as ideal point dipoles, the new model considers the non-uniform magnetization within each particle. This can be achieved by using a multipole expansion of the magnetic field, which includes not only the dipole term but also higher-order terms such as the quadrupole and octupole terms. This is particularly important for accurately calculating the magnetic forces between closely spaced particles [13].
• Realistic Microstructure Generation: The model incorporates realistic microstructures that are representative of those found in actual MREs. This can be achieved through two primary methods:
• Image-based modeling: Using 3D imaging techniques like X-ray microtomography to obtain the actual particle distribution within a real MRE sample.
• Stochastic generation: Using statistical algorithms to generate virtual microstructures that have the same statistical properties (e.g., particle volume fraction, clustering characteristics) as real MREs [14].
• Advanced Interaction Models: The model goes beyond simple pairwise interactions and includes many-body effects. The magnetization of a given particle is not only influenced by the external field but also by the fields produced by all other particles in its vicinity. This can be captured using iterative numerical schemes or by employing effective medium theories that account for the local field corrections.
• Full Magneto-Mechanical Coupling: The new model explicitly couples the magnetic and mechanical behavior of the MRE. The magnetic forces between particles cause the elastomeric matrix to deform, which in turn changes the positions of the particles. This feedback loop is crucial for accurately predicting the magneto-induced stiffness changes and the overall effective magnetic properties [15]. This can be implemented within a finite element framework where the magnetic forces are calculated at each step of the mechanical deformation.
• Incorporation of Interfacial Effects: The model acknowledges the importance of the interface between the magnetic particles and the elastomeric matrix [16,17]. The properties of this interphase region, such as its stiffness and adhesion characteristics, can significantly influence the transmission of forces from the particles to the matrix and can affect the overall magnetic response.

5. MATHEMATICAL FORMULATION OF THE NEW MODEL

The mathematical formulation of the proposed model involves solving a coupled set of equations governing the magnetic and mechanical behavior of the MRE. Magnetic Problem: The magnetic field H and magnetic induction B are related by the constitutive equation: B = μ0​ (H + M). where μ0​ is the permeability of free space and M is the magnetization. The total magnetization of the composite is the volume average of the magnetization of the particles and the matrix. The magnetization of the particles is a non-linear function of the local magnetic field and can be described by a Langevin function or a more sophisticated micromagnetic model [18]. Mechanical Problem: The mechanical behavior of the elastomeric matrix is described by the equations of elasticity. The stress tensor σ is related to the strain tensor ε by the constitutive law of the elastomer, which can be a neo-Hookean, Mooney-Rivlin, or another hyper elastic model. Coupling: The coupling between the magnetic and mechanical problems is achieved through the magnetic forces acting on the particles. The magnetic force on a single particle can be calculated by integrating the Maxwell stress tensor over its surface. This force is then applied as a body force in the mechanical problem, causing the matrix to deform. The deformation of the matrix leads to a new particle configuration, which then requires a re-calculation of the magnetic field. This coupled problem is solved iteratively until a self-consistent solution is obtained [19].

6. EXPERIMENTAL VALIDATION AND FUTURE OUTLOOK

The development of any new theoretical model must be accompanied by rigorous experimental validation. A variety of experimental techniques can be used to measure the effective magnetic properties of MREs, including:

• Vibrating Sample Magnetometry (VSM): To measure the magnetization curves of MRE samples.
• Inductance-based methods: To measure the magnetic permeability of MREs as a function of the applied magnetic field.
• Dynamic Mechanical Analysis (DMA) with a magnetic field: To measure the magneto-induced changes in stiffness and damping.

By comparing the predictions of the new model with experimental data obtained from well-characterized MRE samples, we can assess its accuracy and predictive power. The future of MRE modeling lies in the continued integration of multiphysics phenomena and the development of computationally efficient algorithms. The proposed new model represents a significant step forward from the classical dipole model and provides a more robust framework for understanding and designing these fascinating smart materials. As our computational capabilities continue to grow, we can expect to see even more sophisticated models that can capture the full complexity of MREs, paving the way for the development of novel devices with unprecedented performance and functionality.

7. CONCLUSION

The journey from the simple elegance of the classical dipole model to the comprehensive, Multiphysics framework proposed in this review reflects the evolution of our understanding of magnetorheological elastomers. While the dipole model provided the initial foundational insights, the demand for high-performance MREs has necessitated a more nuanced and accurate approach. The new model, by integrating a more realistic description of particle magnetization, microstructure, and magneto-mechanical coupling, offers a powerful predictive tool for the scientific and engineering communities. It moves beyond simple approximations and embraces the inherent complexity of these composite materials. This allows for a more rational design of MREs with tailored magnetic and mechanical properties, accelerating their transition from laboratory curiosities to real-world technological solutions. The continued development and refinement of such advanced models, in close conjunction with experimental validation, will be the key to unlocking the full potential of magnetorheological elastomers and ushering in a new era of smart materials and adaptive structures.

REFERENCES

1. Bustamante, R. (2011). Magnetorheological Elastomers: Modeling and Identification. VDM Verlag.
2. Jolly, M. R., Bender, J. W., & Carlson, J. D. (1999). Properties and applications of commercial magnetorheological fluids. Journal of Intelligent Material Systems and Structures, 10(1), 5-13.
3. Gong, X. L., Zhang, X. Z., & Zhang, P. Q. (2005). Fabrication and characterization of isotropic magnetorheological elastomers. Polymer Testing, 24(5), 669-676.
4. Kallio, M., Lindroos, T., Aalto, S., & Kallio, P. (2007). A new model for the magnetic properties of magnetorheological elastomers. Journal of Physics D: Applied Physics, 40(19), 5873.
5. Varga, Z., Filipcsei, G., & Zrínyi, M. (2006). Magnetic field induced phase transition in ferrogels. Polymer Gels and Networks, 4(2), 125-133.
6. Wood, D. A., & Wereley, N. M. (2013). Large-strain quasistatic and dynamic characterization of magnetorheological elastomers. Journal of Intelligent Material Systems and Structures, 24(8), 929-940.
7. Jing, K., Li, H., Xiang, H., & Peng, X. (2025). A New Model to Describe the Effective Magnetic Properties of Magnetorheological Elastomers. Physics, 7(2), 21.
8. Vatandoost, H., Sedaghati, R., & Rakheja, S. (2022). A novel methodology for accurate estimation of magnetic permeability of magnetorheological elastomers. Journal of Magnetism and Magnetic Materials, 560, 169669.
9. Kumar, V., Park, S. J., Lee, D. J., & Park, S. S. (2021). Mechanical and magnetic response of magneto‐rheological elastomers with different types of fillers and their hybrids. Journal of Applied Polymer Science, 138(37), 50957.
10. Saber, A., & Sedaghati, R. (2023). The modeling of magnetorheological elastomers: A state‐of‐the‐art review. Advanced Engineering Materials, 25(16), 2300182.
11. Díez, A. G., Tubio, C. R., Etxebarria, J. G., & Lanceros-Mendez, S. (2021). Magnetorheological elastomer‐based materials and devices: state of the art and future perspectives. Advanced Engineering Materials, 23(6), 2100240.
12. Snarskii, A. A., Shamonin, M., & Yuskevich, P. (2021). Effect of magnetic-field-induced restructuring on the elastic properties of magnetoactive elastomers. Journal of Magnetism and Magnetic Materials, 517, 167392.
13. Kumar, V., & Park, S. S. (2024). Stretchable Magneto-Mechanical Configurations with High Magnetic Sensitivity Based on “Gel-Type” Soft Rubber for Intelligent Applications. Gels, 10(1), 80.
14. Stepanov, G. V., Bakhtiiarov, A. V., Lobanov, D. A., Borin, D. Y., Semerenko, D. A., & Storozhenko, P. A. (2022). Magnetoresistive and magnetocapacitive effects in magnetic elastomers. SN Applied Sciences, 4(6), 178.
15. Vasilyeva, M., Nagornov, D., & Orlov, G. (2021). Research on dynamic and mechanical properties of magnetoactive elastomers with high permeability magnetic filling agent at complex magneto-temperature exposure. Materials, 14(9), 2376.
16. Silva, J. A., Gouveia, C., Dinis, G., Pinto, A. M., & Pereira, A. M. (2022). Giant magnetostriction in low-concentration magnetorheological elastomers. Composites Part B: Engineering, 243, 110125.
17. Peng, Z., Zhai, Z., Yang, R., Xu, H., & Wu, Y. (2024). Fabrication and property research of a new 3D-Printable magnetorheological elastomer (MRE). Materials Today Physics, 45, 101467.
18. Izadifard, M. H., Norouzi, M., Ghatee, M., Nazari, M. B., & Alimardan, A. (2024). Experimental characterization and fractional modelling of anisotropic magnetorheological elastomers under the influence of temperature and magnetic fields. Soft Matter, 20(7), 1629-1650.
19. Tahir, S., Usman, M., & Umer, M. A. (2022). Effect of volume fraction on shear mode properties of Fe-Co and Fe-Ni filled magneto-rheological elastomers. Polymers, 14(14), 2968.