Next Article in Journal
Ultrasonic Modification of Ag Nanowires and Their Applications in Flexible Transparent Film Heaters and SERS Detectors
Previous Article in Journal
Crack Healing and Mechanical Properties Recovery in SA 508–3 Steel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Magnetorheological Elastomers Using the Elastic–Plastic Model with Kinematic Hardening

by
Daniel Lewandowski
Department of Mechanics, Material Science and Engineering, Wroclaw University of Science and Technology, 50-372 Wroclaw, Poland
Materials 2019, 12(6), 892; https://doi.org/10.3390/ma12060892
Submission received: 8 February 2019 / Revised: 11 March 2019 / Accepted: 14 March 2019 / Published: 18 March 2019
(This article belongs to the Section Smart Materials)

Abstract

:
The paper describes the modeling process of the magnetomechanical properties of magnetorheological elastomers. Unlike the primary sources in which a viscoelastic model is commonly used, in this work, the elastic–plastic model with linear hardening was adopted. Its parameters were determined from data obtained in an experiment. Measurement data encompassed a range of various strain rates and amplitudes, as well as a range of magnetic field values. Model parameters as a function of a magnetic field were obtained as a result of identification. The correctness of data correlation was shown by comparing hysteresis loops τ ( γ ) . Satisfactory consistency between experimental and model research was achieved on the assumption that the model was applied only to higher strain rates, above the boundary value γ ˙ A = 0.025 ( 1 / s ) .

1. Introduction

The subject matter of this work is the analysis of the properties of a selected magnetorheological elastomer using a mathematical model. The material chosen for the study belongs to a large group of materials in the literature referred to as Smart Materials. Their essential characteristic is the ability to change their properties and parameters resulting from the effect of an external magnetic field [1,2]. In the case of magnetorheological composites and elastomers belonging to the same group of materials, mainly the change of their mechanical properties occurs. Due to their intended use, e.g., in mechanical constructions, mainly stiffness changes [3] or an ability to dissipate energy and damp vibration [4,5] arouse interest. It is the dissipation energy mechanisms of these materials, controlled and changed by a magnetic field, that allow using them in this way. They can also be used as actuators, which requires taking into account their magnetostrictive properties [6]. Stimulation with a magnetic field allows for a continuous change of material parameters. This characteristic means that it is possible to use magnetorheological elastomers in adaptive machines and equipment which can conform to changing operating conditions. However, this way of using materials requires the knowledge of their characteristics, i.e. a model. The description of the constitutive equation, the dependence of stress and strain as a function of a magnetic field, allows correctly constructing devices with these materials. The model is of crucial significance here also from the perspective of control and changes in a given construction. An example to illustrate this is vehicle suspensions that require conformity with current driving conditions, driving style and loads. The implementations of a magnetorheological fluid in semi-active control are presented in paper [7]. An example of using the Bouc–Wen model to describe the operation of dampers with a magnetorheological fluid in a car suspension can be found in [8].
The analysis of literature data on manufacturing, testing and modeling the properties of magnetorheological elastomers shows that a model that is frequently used in their description is the linear, viscoelastic Kelvin–Voight model [9,10]. Its popularity may be explained by the fact that it is commonly applied in the rheology of plastics where the warp of a magnetorheological elastomer is made of a polymer [11]. Unfortunately, this model is substantially simplified and does not allow correctly projecting some types of behavior, e.g., the influence of strain rate. In some scientific publications presenting research results related to these materials, their authors avoid showing stress and strain as a function of time. They offer only the so-called averaged analysis with a storage and loss modulus [12,13,14]. Unfortunately, these parameters are varied depending on frequency excitation, which is a mistake.
The proper description of magnetorheological elastomers with the use of the Kelvin–Voight model is related to their internal structure. Their crucial component is a ferromagnetic filling in the form of a magnetically soft powder, e.g., carbonyl iron [2] ASC 300 [15,16]. Mutual interaction between metallic particles with sizes of a few micrometers [17] is of crucial influence on its global properties, called the magnetorheological effect. It is generally defined as the dependence of stiffness and energy damping as a function of a magnetic field. Using the viscous element occurring in the Kelvin–Voight model as the factor responsible for an increase in stress linearly dependent on deformation speed is not sufficient. A correct description also requires the use of other compounds describing energy dissipation inside a material. A possible solution is the use of ideal-plastic elements [18,19].
This paper is divided into sections. In the initial Section 2, the MRE material and testing samples are discussed. Next, the description of the measurement stand and the test procedure are presented in Section 3. The most important part of the article is the model and process for the identification of its parameters, taking into account the influence of the magnetic field—Section 5 and Section 6. The work was completed by comparing the hysteresis loops obtained from the model and experimental data.

2. Magnetoreheological Elastomer and Specimens

The type of magnetorheological elastomer used and analyzed in this work was developed a few years ago. It consists of two main components: polymer matrix made of thermoplastic elastomer SEBS and a ferromagnetic filling—ASC300 iron type. Its detailed manufacturing procedure was described in a previous work by the author [20], where the experimental research of some properties of this type of materials was conducted. For the analysis of the model investigated in this work, a separate group of cylindrical shape samples with a height of 25 mm and a diameter of 25 mm was made. In the manufacturing process, the material was not polarized with an external magnetic field, and its internal structure can be considered anisotropic.

3. Test Stand and Magnetic Field Stimulation

The data used to identify model parameters were recorded during an experiment specially prepared for this purpose. A sample of a magnetorheological elastomer was cyclically loaded with simultaneous magnetic field simulation. Loading was conducted using a hydraulic device with a feedback system, which allowed obtaining a set of constant strain rate values. The magnetic field was generated using the so-called Halbach matrix constructed using a system of permanent magnets [21]. The matrix scheme and the method used to place the tested sample inside it are presented in Figure 1, showing: (1) the internal magnetic ring (fixed); (2) the external magnetic ring (rotating); (3) the cylindrical material sample; (4) the movable top plate attached to the sample (to transfer external force); and (5) the base. Red and blue vectors indicate action directions inside particular magnets, while the yellow ones represent a resultant magnetic field inside the matrix. The construction of the two-ring array allowed obtaining a directed magnetic field whose intensity could be adjusted by hanging the rotation angle. The measurement of magnetic field strength was conducted inside the matrix. The obtained characteristic is presented in Figure 2b. A change of the rotation angle between the rings causes the strengthening or weakening of the magnetic field inside it. This is due to the total sum of the magnetic field of all individual magnets that make up the matrix. This way of loading the sample and applying the magnetic field was described in detail in one of the previous works of the author [22].

4. Experimental Data

The samples made of a magnetorheological elastomer were deformed in a process similar to pure shearing. As shown in Figure 1, the top surface was moved parallel to the basis, and the investigated material was attached between them. The specimens were glued to the plates with a cyanoacrylate glue dedicated to rubber bonding—Loctite 480. The displacement measurement allowed calculating shear transformation in the form of angle γ . Stress values were calculated from the measurements of a force necessary to move the top surface and the sample cross-section. The character of deformation was a linear function cyclically changing at the constant speed. A sample waveform is presented in Figure 3. Hysteresis loops in the τ ( γ ) coordination system resemble an ellipsis; however, they have a sharp end and a straight line segment. Figure 4a shows hysteresis loops obtained for various frequency values (hence, also strain rates γ ˙ ) at a constant deformation amplitude and Figure 4b presents different deformation amplitude values at a constant strain rate.

5. Selection of Material Model

The model was selected based on the earlier experiments of the author and literature studies. The commonly used linear viscoelastic model, also called the Kelvin–Voigt model, despite its simplicity does not ensure providing a correct description of the behavior of magnetorheological materials. The main reason is the occurrence of some interactions inside the material structure. Magnetorheological elastomers combine two main components, the first one is an active ferromagnetic filling in the form of particles, spheres, and flakes of, e.g., iron, and the other one is a matrix, which may be made of, e.g., thermoplastic elastomers, as was the case in this work. The matrix is a binding element, however it simultaneously results in small, local dislocations in the first phase. Together, both elements constitute one uniform structure. The compact filling of space and the adherence of iron particles to each other allows, among others, for some interactions between them. As a result of material deformation on a global scale, friction takes place on the border between the particles of the ferromagnetic filling on a micro-scale. Forces generated as a consequence of these interactions depend on mutual attraction forces, which in turn are induced by an external magnetic field. Based on this type of image, an attempt was made to describe magnetorheological elastomers using dissipative elements that were independent of deformation speed and simultaneously nonlinear. The author previously used the adopted model in works depicting the so-called magnetorheological composites [23,24]. They belong to the same group of Smart Materials activated magnetically. However, they differ in their internal structure and the occurrence of a viscous liquid. The scheme of the adopted model is presented in Figure 5.
It is composed of two parallel connected branches and three elements ( τ o 2 is a shear yield point, and G 1 and G 2 are Kirchhoff’s modulus). The elastic elements cooperate with the ideally plastic element, which takes into account various types of hardening, i.e. kinematic, isotropic and mixed. The constitutive compounds describing the dependence of stress on deformation are presented below:
τ = τ 1 + τ 2 and γ = γ 1 = γ 2 ,
where:
τ 1 = 2 G 1 γ ,
τ 2 = 2 G 2 ( γ γ a ) + τ o 2 for γ ˙ < 0 and γ > τ o 2 / G 2 + γ a τ o 2 for γ ˙ < 0 and γ > τ o 2 / G 2 + γ a 2 G 2 ( γ + γ a ) τ o 2 for γ ˙ > 0 and γ < τ o 2 / G 2 γ a τ o 2 for γ ˙ > 0 and γ > τ o 2 / G 2 γ a
The adopted model, contrary to the Kelvin–Voigt model, does not consider the dependence on strain rate. It is a rather untypical approach because elastomer models are based mainly on the use of the so-called viscous liquid. In the case of the described magnetorheological elastomers, such dependencies were avoided. The characteristics obtained from the conducted experiments, which allowed observing the interaction of the material at various strain rates, are presented in Figure 6. It was found that, above a certain level of strain rate, stresses in the material state had a constant value. This allows adopting the assumption of an insignificant or non-existent influence of rate or frequency. It is justified when the used magnetorheological elastomer is employed to work with higher frequencies. In the described case, the border point can be assumed to be γ ˙ A = 0.025 ( 1 / s ) .

6. Identification of Model Parameters

The identification of the first two parameters: G 1 and τ o 2 was conducted on the basis of experimental data describing only the material loading process (without unloading). The selected parts of the hysteresis loop were used to find the intersection point with γ -axis and the straight line inclination angle. Such a solution allowed to simultaneously determining these parameters for the model. The linear function approximation was made using the linear regression method based on the least squares method. The consistency of the following function τ = γ G 1 + τ o 2 was analyzed for each loop separately, for both the top and bottom parts of the loop. Calculations were made for all available measurements of all amplitudes and strain rates. Finally, the results were averaged separately for each value of the magnetic field. The consistency criterion being the amount of dissipated energy in one cycle was used to determine parameter G 2 . It was decided to compare the size of loops obtained for the model and the experimental data. The consistency criterion was defined as follows:
δ = Δ W E Δ W M ( G 2 ) ,
where Δ W E is the area calculated from experimental data, while Δ W M ( G 2 ) is the value obtained in the model with set parameter G 2 . The calculations were made in such a way that δ was minimized and then the value of parameter G 2 was selected. Similar to the first two parameters, the calculations were made separately for each amplitude and deformation speed value and then they were averaged. The procedure was repeated subsequently for the incremental values of a magnetic field. The final identification results are presented in Figure 7. The G 1 and G 2 value order was similar. However, the character of their course as a function of a magnetic field was different. The value of parameter G 1 increased with a magnetic field, while G 2 decreased. Parameter τ o 2 behaved as expected; it grew with the increase in a magnetic field. As explained above, its increase means a rise in interaction forces between ferromagnetic particles inside the material and an increase in friction forces inside the material.
Following the analysis of the correctness of model choice and the conducted identification, a set of hysteresis loops obtained from the model are presented in Figure 8. They were calculated for the same parameters for which the experimental data shown in Figure 4b were recorded. A detailed comparison of individual loops is presented in Figure 9. The differences between the model and experimental loops are not significant and depend on the deformation amplitude.
The set of loops obtained from the model for increasing values of a magnetic field is shown in Figure 10. The model correctly maps changes in loop area and curve inclination during loading, which can be seen by comparing them with the loops presented in Figure 4. The degree to which the model reflects the loops obtained in the experiment was determined by comparing particular waveforms. For this purpose, changes in the areas of hysteresis loops and maximum stress values were analyzed. The differences related to all examined strain amplitudes, magnetic field and frequency did not exceed 10%.

7. Summary and Conclusions

The identification process was completed with the determination of the model parameters. Satisfactory results have been obtained for the entire range of the magnetic field used. The course of variation of model parameters has a linear character of dependence on the magnetic field. The correctness of the obtained model was shown in comparison with the results obtained in real studies. Despite the good correlation of experimental data and the model, it should be noted that the adopted equations could not reflect the influence of strain rate changes. The loops obtained for this model are identical, regardless of the frequency of the excitation signals. Hence, the loops presented in Figure 4 are described with only one loop. The limitation of low frequencies did not seem to be a considerable simplification. Regardless of these imperfections, the description of the investigated magnetorheological elastomer obtained using this model can be considered correct. The obtained results are better than in the case of using the viscoelastic Kelvin–Voigt model where the linear viscous function could not correctly reproduce the behavior of the magnetorheological elastomer type for a wide range of strain rate.

Funding

This research was funded by the Polish Ministry of Science and Higher Education, grant No. 0401/0012/18.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Li, Y.; Li, J.; Li, W.; Du, H. A state-of-the-art review on magnetorheological elastomer devices. Smart Mater. Struct. 2014, 23, 123001. [Google Scholar] [CrossRef] [Green Version]
  2. Ahamed, R.; Choi, S.B.; Ferdaus, M.M. A state of art on magneto-rheological materials and their potential applications. J. Intell. Mater. Syst. Struct. 2018, 29, 2051–2095. [Google Scholar] [CrossRef]
  3. Kallio, M. The Elastic and Damping Properties of Magnetorheological Elastomers; VTT Publications: Espoo, Finland, 2005; pp. 3–146. [Google Scholar]
  4. Yan, Y.X.; Hui, L.X.; Yu, M.; Fu, J.; Dong, L.H. Dynamic response time of a metal foam magneto-rheological damper. Smart Mater. Struct. 2013, 22, 025026. [Google Scholar] [CrossRef]
  5. Chen, L.; Gong, X.-L.; Li, W.-H. Damping of Magnetorheological Elastomers. Chin. J. Chem. Phys. 2008, 21, 581–585. [Google Scholar] [CrossRef] [Green Version]
  6. Guan, X.; Dong, X.; Ou, J. Magnetostrictive effect of magnetorheological elastomer. J. Magn. Magn. Mater. 2008, 320, 158–163. [Google Scholar] [CrossRef]
  7. Carlson, J.D. Implementation of Semi-Active Control Using Magneto-Rheological Fluids. IFAC Proc. Volumes 2000, 33, 905–910. [Google Scholar] [CrossRef]
  8. Prabakar, R.; Sujatha, C.; Narayanan, S. Optimal semi-active preview control response of a half car vehicle model with magnetorheological damper. J. Sound Vib. 2009, 326, 400–420. [Google Scholar] [CrossRef]
  9. Kukla, M.; Górecki, J.; Malujda, I.; Talaśka, K.; Tarkowski, P. The Determination of Mechanical Properties of Magnetorheological Elastomers (MREs). Procedia Eng. 2017, 177, 324–330. [Google Scholar] [CrossRef]
  10. Agirre-Olabide, I.; Kuzhir, P.; Elejabarrieta, M. Linear magneto-viscoelastic model based on magnetic permeability components for anisotropic magnetorheological elastomers. J. Magn. Magn. Mater. 2018, 446, 155–161. [Google Scholar] [CrossRef] [Green Version]
  11. Puente-Córdova, J.; Reyes-Melo, M.; Palacios-Pineda, L.; Martínez-Perales, I.; Martínez-Romero, O.; Elías-Zúñiga, A. Fabrication and Characterization of Isotropic and Anisotropic Magnetorheological Elastomers, Based on Silicone Rubber and Carbonyl Iron Microparticles. Polymers 2018, 10, 1343. [Google Scholar] [CrossRef]
  12. Lu, X.; Qiao, X.; Watanabe, H.; Gong, X.; Yang, T.; Li, W.; Sun, K.; Li, M.; Yang, K.; Xie, H.; et al. Mechanical and structural investigation of isotropic and anisotropic thermoplastic magnetorheological elastomer composites based on poly(styrene-b-ethyleneco- butylene-b-styrene) (sebs). Rheol. Acta 2011, 51, 37–50. [Google Scholar] [CrossRef]
  13. Böse, H.; Röder, R. Magnetorheological Elastomers with High Variability of Their Mechanical Properties. J. Phys. Conf. Ser. 2009, 149, 012090. [Google Scholar] [CrossRef]
  14. Tong, Y.; Dong, X.; Qi, M. Improved tunable range of the field-induced storage modulus by using flower-like particles as the active phase of magnetorheological elastomers. Soft Matter 2018, 14, 3504–3509. [Google Scholar] [CrossRef]
  15. Zając, P.; Kaleta, J.; Lewandowski, D.; Gasperowicz, A. Isotropic magnetorheological elastomers with thermoplastic matrices: structure, damping properties and testing. Smart Mater. Struct. 2010, 19, 045014. [Google Scholar] [CrossRef]
  16. Stepanov, G.; Borin, D.; Odenbach, S. Magnetorheological effect of magneto-active elastomers containing large particles. J. Phys. Conf. Ser. 2009, 149, 012098. [Google Scholar] [CrossRef] [Green Version]
  17. Cvek, M.; Mrlik, M.; Sevcik, J.; Sedlacik, M. Tailoring Performance, Damping, and Surface Properties of Magnetorheological Elastomers via Particle-Grafting Technology. Polymers 2018, 10, 1411. [Google Scholar] [CrossRef]
  18. Melenev, P.V.; Raikher, Y.L.; Rusakov, V.V.; Stepanov, G.V. Field-induced plasticity of soft magnetic elastomers. J. Phys. Conf. Ser. 2009, 149, 012094. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, D.H.; Liao, W.H. Magnetorheological fluid dampers: A review of parametric modelling. Smart Mater. Struct. 2011, 20, 023001. [Google Scholar] [CrossRef]
  20. Kaleta, J.; Królewicz, M.; Lewandowski, D. Magnetomechanical properties of anisotropic and isotropic magnetorheological composites with thermoplastic elastomer matrices. Smart Mater. Struct. 2011, 20, 085006. [Google Scholar] [CrossRef]
  21. Halbach, K. Design of permanent multipole magnets with oriented rare earth cobalt material. Nuclear Instrum. Methods 1980, 169, 1–10. [Google Scholar] [CrossRef]
  22. Bocian, M.; Kaleta, J.; Lewandowski, D.; Przybylski, M. Tunable Absorption System based on magnetorheological elastomers and Halbach array: Design and testing. J. Magn. Magn. Mater. 2017, 435, 46–57. [Google Scholar] [CrossRef]
  23. Kaleta, J.; Lewandowski, D.; Zietek, G. Inelastic properties of magnetorheological composites: II. Model, identfication of parameters. Smart Mater. Struct. 2007, 16, 1954–1960. [Google Scholar] [CrossRef]
  24. Kaleta, J.; Lewandowski, D. Inelastic properties of magnetorheological composites: I. Fabrication, experimental tests, cyclic shear properties. Smart Mater. Struct. 2007, 16, 1948–1953. [Google Scholar] [CrossRef]
Figure 1. Scheme of magnetorheological elastomer sample placed inside the Halbach matrix: (a) view of an array with the sample; and (b) cross-section. Notation is presented in the text.
Figure 1. Scheme of magnetorheological elastomer sample placed inside the Halbach matrix: (a) view of an array with the sample; and (b) cross-section. Notation is presented in the text.
Materials 12 00892 g001
Figure 2. (a) Photo of test setup system: 1, upper fix; 2, adjusting lever; 3, external Halbach ring; 4, external Halbach ring; 5, specimen; 6, lower fix; 7, force transducer. (b) Characterization of magnetic field dependence on rotation angle α of external ring Halbach matrix.
Figure 2. (a) Photo of test setup system: 1, upper fix; 2, adjusting lever; 3, external Halbach ring; 4, external Halbach ring; 5, specimen; 6, lower fix; 7, force transducer. (b) Characterization of magnetic field dependence on rotation angle α of external ring Halbach matrix.
Materials 12 00892 g002
Figure 3. Experimental research results: (a) stress and deformation as a function of time for selected value of magnetic field and deformation amplitude; and (b) hysteresis loops τ ( γ ) for a few selected values of magnetic field at constant values of strain amplitude and strain rate.
Figure 3. Experimental research results: (a) stress and deformation as a function of time for selected value of magnetic field and deformation amplitude; and (b) hysteresis loops τ ( γ ) for a few selected values of magnetic field at constant values of strain amplitude and strain rate.
Materials 12 00892 g003
Figure 4. Experimental research results: (a) hysteresis loops for various excitation frequencies, at the selected value of magnetic field and deformation amplitude; and (b) hysteresis loops for different strain amplitudes, at the chosen value of magnetic field and deformation amplitude.
Figure 4. Experimental research results: (a) hysteresis loops for various excitation frequencies, at the selected value of magnetic field and deformation amplitude; and (b) hysteresis loops for different strain amplitudes, at the chosen value of magnetic field and deformation amplitude.
Materials 12 00892 g004
Figure 5. The scheme of elastic–plastic model with kinematic hardening.
Figure 5. The scheme of elastic–plastic model with kinematic hardening.
Materials 12 00892 g005
Figure 6. Dependence of stress inside magnetorheological elastomer on strain rate for various strain amplitudes.
Figure 6. Dependence of stress inside magnetorheological elastomer on strain rate for various strain amplitudes.
Materials 12 00892 g006
Figure 7. Identification results. Model parameters as a function of magnetic field value: (a) stiffness G 1 and G 2 ; and (b) yield point τ o 2 .
Figure 7. Identification results. Model parameters as a function of magnetic field value: (a) stiffness G 1 and G 2 ; and (b) yield point τ o 2 .
Materials 12 00892 g007
Figure 8. Identification results. Hysteresis loops obtained from the model for three values of deformation amplitude and the set value of magnetic field.
Figure 8. Identification results. Hysteresis loops obtained from the model for three values of deformation amplitude and the set value of magnetic field.
Materials 12 00892 g008
Figure 9. Direct comparison of experimental and model results. Selected hysteresis loops for two different values of strain amplitude: (a) γ A = 0.035 ; (b) γ A = 0.067
Figure 9. Direct comparison of experimental and model results. Selected hysteresis loops for two different values of strain amplitude: (a) γ A = 0.035 ; (b) γ A = 0.067
Materials 12 00892 g009
Figure 10. Set of hysteresis loops obtained from the model for subsequently incrementing values of magnetic field.
Figure 10. Set of hysteresis loops obtained from the model for subsequently incrementing values of magnetic field.
Materials 12 00892 g010

Share and Cite

MDPI and ACS Style

Lewandowski, D. Modeling of Magnetorheological Elastomers Using the Elastic–Plastic Model with Kinematic Hardening. Materials 2019, 12, 892. https://doi.org/10.3390/ma12060892

AMA Style

Lewandowski D. Modeling of Magnetorheological Elastomers Using the Elastic–Plastic Model with Kinematic Hardening. Materials. 2019; 12(6):892. https://doi.org/10.3390/ma12060892

Chicago/Turabian Style

Lewandowski, Daniel. 2019. "Modeling of Magnetorheological Elastomers Using the Elastic–Plastic Model with Kinematic Hardening" Materials 12, no. 6: 892. https://doi.org/10.3390/ma12060892

APA Style

Lewandowski, D. (2019). Modeling of Magnetorheological Elastomers Using the Elastic–Plastic Model with Kinematic Hardening. Materials, 12(6), 892. https://doi.org/10.3390/ma12060892

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop