A unified numerical approach for soft to hard magneto-viscoelastically coupled polymers

TL;DR

This paper introduces a comprehensive finite element framework for simulating magneto-viscoelastic polymers, capable of handling both soft and hard magnetic behaviors, compressibility, and time-dependent effects, validated through practical examples.

Contribution

A unified numerical approach for magneto-viscoelastic polymers that integrates soft and hard magnetic behaviors, compressibility, and viscoelasticity within a finite element framework.

Findings

The framework accurately models magneto-mechanical responses of MAPs.

Viscoelastic parameters significantly influence the response characteristics.

The method successfully simulates boundary value problems like beams and robotic grippers.

Abstract

The last decade has witnessed the emergence of magneto-active polymers (MAPs) as one of the most advanced multi-functional soft composites. Depending on the magnetisation mechanisms and responsive behaviour, MAPs are mainly classified into two groups: i) hard magnetic MAPs in which a large residual magnetic flux density sustains even after the removal of the external magnetic field, and ii) soft magnetic MAPs where the magnetisation of the filler particles disappear upon the removal of the external magnetic field. Polymeric materials are widely treated as fully incompressible solids that require special numerical treatment to solve the associated boundary value problem. Furthermore, both soft and hard magnetic particles-filled soft polymers are inherently viscoelastic. Therefore, the aim of this paper is to devise a unified finite element method-based numerical framework for magneto-mechanically coupled systems that can work for compressible and fully incompressible materials and from hard to soft MAPs, including the effects of the time-dependent viscoelastic behaviour of the underlying matrix. First, variational formulations for the uncoupled problem for hard MAPs and the coupled problem for soft MAPs are derived. The weak forms are then discretised with higher-order B\'ezier elements while the evolution equation for internal variables in viscoelastic models is solved using the generalised-alpha time integration scheme, which is implicit and second-order accurate. Finally, a series of experimentally-driven boundary value problems consisting of the beam and robotic gripper models are solved in magneto-mechanically coupled settings, demonstrating the versatility of the proposed numerical framework. The effect of viscoelastic material parameters on the response characteristics of MAPs under coupled magneto-mechanical loading is also studied.