1. Introduction
MREs are smart materials that are composed of a polymer matrix filled with magnetic particles, usually iron or iron oxide, which can respond to an external magnetic field. The magnetic particles are typically dispersed within the elastomer matrix in a random manner, forming a network of interconnected chains. When an external magnetic field is applied to the MRE, the magnetic particles within the polymer matrix align themselves with the direction of the field, causing the elastomer to stiffen and become more rigid. Conversely, when the magnetic field is removed, the particles return to their random orientation and the elastomer returns to its original soft and flexible state [
1]. The stiffness of MRE can be controlled by adjusting the strength and orientation of the magnetic field.
MRE have the potential for a range of applications in various industries. Some key examples include vibration control [
2] and damping systems [
3,
4], where MRE can reduce vibrations and enhance the stability of structures and machinery. MRE can also be utilized in the development of soft robotics [
5] and prosthetics with tunable stiffness [
6], offering improved flexibility and control. Additionally, MRE have potential for seismic protection [
7], as they can enhance the seismic resistance of structures by providing adaptive damping to reduce the impact of seismic waves. In the aerospace and defense industries, MRE can be integrated into aircraft and spacecraft design to improve vibration control and damping [
8], leading to better performance and stability. Furthermore, MRE may be applied to medical device design [
9], such as stents and catheters, to offer greater control and flexibility. As research in this field progresses, the applications of MRE continue to expand, showcasing their potential for diverse applications in different industries.
MRE changes in their mechanical properties in response to an applied magnetic field. Due to their unique behavior, accurate modeling of MRE is crucial for their effective application in various industries. Accurate modeling of MRE allows for the prediction of their behavior under different magnetic field strengths, which is essential in designing MRE-based devices and systems. It also enables the optimization of MRE for specific applications by predicting the mechanical properties, such as stiffness and damping, and the response time of MREs to changes in the magnetic field. Furthermore, accurate modeling of MREs can lead to the development of more advanced MRE-based systems that require a more precise understanding of the behavior of MRE to ensure optimal performance.
There are two methods in modeling the behavior of MRE, namely parametric and non-parametric models [
10]. Parametric models refer to models that are based on a set of predefined parameters or assumptions about the behavior of the material. These models typically involve the use of mathematical equations and require knowledge of the material's properties and characteristics, such as its stiffness and damping. On the other hand, non-parametric models do not rely on pre-defined assumptions or parameters, but instead use data-driven approaches to develop models based on observed behavior. These models can include machine learning algorithms, such as artificial neural networks or support vector machines that learn from data to predict the behavior of the material.
Both parametric and non-parametric models have been used in MRE research. Parametric models, such as micromechanical models or continuum mechanics models, have been developed to describe the behavior of MRE based on the physical properties of the material. Non-parametric models, such as artificial neural networks [
11] or fuzzy logic systems [
12], have been used to predict the behavior of MRE based on experimental data. Overall, both approaches have their strengths and weaknesses, and the choice of which model to use depends on the specific application and available data. Parametric models can provide a more fundamental understanding of the material behavior, but may require more knowledge of the material properties. Non-parametric models can be more flexible and adaptable, but may require more experimental data to develop accurate models [
13].
In this paper, an impact behavior model for MRE is proposed which takes into account the hysteresis characteristics of MRE. To model the hysteresis characteristics, multilayers exponential based Preisach model is proposed. Exponential functions are commonly used in modeling hysteresis as they can capture the nonlinear behavior of the material. The multilayer exponential function allows for a more complex model that can better capture the behavior of MRE. To optimize the parameters of the model, an optimization tool namely PSO is used. PSO is a nature-inspired optimization algorithm that simulates the social behavior of swarms in nature, such as flocks of birds or schools of fish. The algorithm is designed to find the optimal solution to a given optimization problem by iteratively adjusting a group of particles or agents in a search space.
2. Design and Fabrication of MRE based Double Acting Actuator
In traditional MRE, the actuator works mainly in compression mode. However, in this study, a double-acting MRE actuator which can work in both compression and extension modes was developed and fabricated. The schematic diagram and the prototype of the MRE actuator is shown in
Figure 1a and
Figure 1b, respectively. MRE actuators consist of a piston and cylinder containing MRE material, a coil or electromagnet and a power source. When a mechanical vibration causes the piston to move relative to the cylinder, a magnetic field is generated by the coil or electromagnet that causes the magnetic particles in the MRE to align them in the direction of the field. This changes the stiffness and damping properties of the MRE, providing a damping force that opposes the motion of the piston. The damping force can be adjusted by varying the strength and direction of the magnetic field, which is controlled by the controller.
The fabrication process of MRE typically involves several steps. First, an elastomer matrix is selected to meet specific application requirements. Next, magnetic particles are selected and their size is determined. Then, the magnetic particles are prepared by mixing them with the elastomer. The mixture is cured by pouring it into a mold. Finally, the MRE is magnetized by exposing them to a strong magnetic field, which magnetizes the particles and induces the desired mechanical properties [
14]. The coils are made of copper wire coated with a layer of enamel insulation with the diameter of 0.7 mm and the number of coils turns is 250. The parameters and the composition for MRE fabrication is shown in
Table 1.
3. Experimental Study on the Force-Displacement Characteristics of MRE for Impact Response
Upon fabricating the proposed MRE, a series of drop impact tests were conducted to evaluate the effectiveness of the MRE under impact loads by analyzing their force-displacement characteristics. The drop impact test is performed to assess the ability of a material to resist impact loading. The test involves several steps. First, the specifications for the test, such as the height of the drop, the weight of the impactor, and other relevant factors, are determined. Next, the test specimen or product is prepared to ensure that it is in suitable condition for testing. The impactor is then secured to the drop apparatus, carefully aligned and oriented. The impactor is raised to the predetermined height and released, allowing it to free-fall and strike the test specimen. The impact is observed and to ensure accurate and reliable results, the test is typically repeated several times. Finally, the results can be analyzed.
Figure 2 depicts the drop test machine that was utilized in the experiment.
The tests were carried out by subjecting the material to a sudden external forces to observe its behavior. The Instron Drop Impact Machine and CEAST Software were utilized to set the experiment variables, including the impact energy, impact velocity, falling height, total mass, and applied current. The parameters set to the Instron Drop Impact Machine during experiment is shown in the
Table 2.
Figure 2 shows the experimental setup that was conducted in the Automotive Lab at Universiti Pertahanan Nasional Malaysia (UPNM).
The experiments were conducted to investigate the behavior of the MRE under different current inputs to the coils. In order to achieve this, varying current values were applied to the coils at 0, 0.5, 1, 1.5, and 2 Amperes. To ensure the accuracy and consistency of the data collected, each experiment was repeated multiple times at each current level. This approach helps to reduce the effects of random errors and improves the reliability of the experimental results.
Figure 3 presents the experimental results displaying the force-displacement relationship with respect to the different current inputs. The figure clearly illustrates the changes in stiffness and damping properties of the MRE with varying current inputs. Upon analyzing each force-displacement curve, it was found that the upper slope of the curve represents the response of the MRE undergoing compression, while the lower slope represents the response of the MRE undergoing retraction.
4. Modeling the Hysteresis Characteristics of MREs using Multi-Layers Exponential based Preisach Model Optimized with PSO
This section describes the algorithm for a multi-layer exponential-based Preisach model that was developed using experimental data on the force-displacement characteristics of MRE under impact loading conditions for input currents injected to the coils at 0, 0.5, 1, 1.5, and 2 Amperes. Additionally, the interpolation method used to obtain the force produced by the MRE for input currents between the specified values is explained. Finally, an optimization tool called Particle Swarm Optimization (PSO) is also discussed.
4.1. Multi-Layers Exponential based Preisach Model
Preisach model is a mathematical tool used to describe hysteresis, which is a phenomenon in which the output of a system depends not only on the current input but also on its history [
16,
17]. The model consists of a set of hysterons, which are basic units that represent the behavior of the system. Each hysteron is associated with a particular input/output history and has a unique switching threshold. The model describes the behavior of the system by representing it as a distribution of hysterons, each of which contributes to the overall response of the system [
18]. The Preisach model consists of many relay hysterons connected in parallel, given weights, and summed. This can be visualized by a block diagram as shown in
Figure 4. Each of these relays has different
α and
β thresholds and is scaled by
µ. With increasing the number of
N, the true hysteresis curve can be approximated better [
19].
An exponential function is a mathematical function in the form of
. Exponential functions are used to model phenomena that grow or decay at a constant percentage rate over time, such as population growth, radioactive decay, or compound interest. They have a characteristic curve that starts out rapidly and then increases or decreases slowly. First, second and
n-th hysterons are defined as follows:
The general relationship between input-output of the Preisach model is written as follows:
Referring to
Figure 3, there are five hysteresis loops representing the force-displacement characteristics of MRE under impact loading for different values of the current injected into the coils. The currents include 0, 0.5, 1, 1.5, and 2 Amperes. Each loop has two sets up Preisach model to represent the upper and lower slopes.
In this proposed model, some controlled parameters are optimized using PSO to accurately represent the hysteresis response of MREs under impact loadings, they are:
,
, …,
;
,
, …,
;
,
, …,
;
,
, …,
; … ;
,
, …,
;
,
, …,
.
As previously mentioned, the proposed model is developed based on experimental data obtained from force-displacement measurements for input currents of 0, 0.5, 1, 1.5, and 2 Amperes. To calculate the force generated by MREs for input currents between 0 A and 0.5 A, an interpolation approach is used based on the known or specified displacement. The interpolation algorithm is detailed in
Figure 5, which uses the experimental data for 0 Ampere and 0.5 Ampere as an example. Similar algorithm is used to calculate the force generated by MREs for input currents between 0.5 A and 1 A, between 1 A and 1.5 A, and between 1.5 A and 2 A. The algorithm for calculating force produced by MREs is as follows:
where,
: Current injected to the MREs coils (between 0A to 0.5 A)
: Force produced by MRE at current i during compression
: Force produced by MRE at current i during retraction
: Force at upper slope of 0A
: Force at lower slope of 0A
: Force at upper slope of 0.5 A
: Force at lower slope of 0.5A
4.2. Optimization of Multi-Layers Sigmoidal Functions using PSO
PSO is a metaheuristic optimization algorithm that is inspired by the social behavior of bird flocking or fish schooling [
20]. In PSO, a set of particles are initialized randomly in the search space and move towards the optimal solution by updating their positions based on their own previous positions and the positions of the best performing particles in the swarm. The basic theory of PSO can be summarized in the following steps [
21]:
Initialization: A population of particles is generated randomly in the search space. Each particle represents a potential solution to the optimization problem.
Evaluation: The fitness of each particle is evaluated based on the objective function to be optimized.
Update of the particle's best position: Each particle keeps track of the best position it has visited so far, denoted as Pbest. If the fitness of the current position is better than its Pbest, the particle updates its Pbest.
Update of the swarm's best position: The best position among all the Pbest positions of the particles is denoted as Gbest. If the fitness of the current Gbest position is better than its previous value, the swarm updates its Gbest.
-
Update of the particle's velocity and position: Each particle updates its velocity and position based on its current velocity, its distance from its Pbest, and its distance from the Gbest. The velocity and position updates are given by the following equations:
Position update:
where
and
are the velocity and position of particle
i at time
t, w is the inertia weight,
and
are the acceleration coefficients, and
and
are random numbers between 0 and 1.
Termination: The algorithm terminates when a stopping criterion is met, such as reaching a maximum number of iterations or a satisfactory fitness level.
By iterating through these steps, the particles in the swarm collectively move towards the optimal solution of the optimization problem.
In this study, the parameters of multi-layers exponential based Preisach model that will be tuned using PSO are , , …, ; , , …, ; , , …, ; , , …, ; … ; , , …, ; , , …, . In PSO, there are several optimization parameters that need to be set to ensure the algorithm performs optimally. The choice of parameter values for PSO can depend on the specific problem being solved and the characteristics of the search space. A common approach is to use a parameter tuning method, such as grid search or random search, to find the optimal parameter values. These parameters include:
Swarm size: 40
Maximum number of iterations: 100
Inertia weight: 0.9
Acceleration coefficients ( and ): 1.42
5. Results and Discussions
This section presents a comparison between the simulated responses of the proposed MRE model and the corresponding experimental data. The maximum error of the predicted force will also be analyzed and discussed. Finally, the accuracy of the model will be tested and validated using input currents ranging from 0 to 0.5 A, 0.5 to 1 A, 1 to 1.5 A, and 1.5 to 2 A. Effects of varying the swarm size and the number of iterations on the PSO are also analyzed.
5.1. Comparison between the Simulated Response of Multi-Layers Exponential based Preisach Model with the experimental Data
Figure 6 compares the simulated model response to experimental data for no current injected into the coils. The model's response closely matches the experimental data, indicating a high degree of accuracy in the simulation. However, a maximum error of 5% occurred in area A, which may be attributed to limitations in the experimental setup or inaccuracies in the simulation model. Further investigation and refinement of the model could potentially reduce this error and improve the overall accuracy of the simulation.
In
Figure 7, the simulated model response is compared to experimental data for an input current of 0.5A injected into the coils. The simulation results exhibit a high degree of accuracy and closely match the experimental data, but the maximum error occurs in area B with a percentage error of approximately 10%. This discrepancy could be due to limitations in the experimental setup, variations in the material properties, or uncertainties in the simulation model. Despite this limitation, the close agreement between the simulated and experimental data validates the model's capability to predict the system's behavior accurately under different input conditions.
Figure 8 compares the simulated model response to experimental data for an input current of 1A injected into the coils. The simulation results exhibit a high degree of accuracy and closely match the experimental data, but the maximum error is observed in area C with a percentage error of approximately 4.8%. Again, this discrepancy may be due to experimental limitations, material property variations, or uncertainties in the simulation model. However, the close agreement between the simulated and experimental data provides confidence in the model's ability to accurately predict the system's behavior.
Figure 9 compares the simulated model response to experimental data for the input current of 1.5A injected into the coils. The simulation results closely match the experimental data, demonstrating the model's high degree of accuracy in predicting the system's response under this input condition. However, it is worth noting that the maximum error occurs in area D, with a percentage error of approximately 4.7%. The close agreement between the simulated and experimental data validates the model's capability to accurately predict the system's behavior.
Lastly, in
Figure 10, the simulated model response is compared to experimental data for an input current of 2A injected into the coils. The simulation results closely match the experimental data, indicating that the model accurately predicts the system's response under this input condition. However, the maximum error is observed in area E with a percentage error of approximately 8.1%. Nevertheless, the close agreement between the simulated and experimental data validates the model's capability to predict the system's behavior accurately. The maximum error of predicted force of the proposed model is summarized in
Table 3.
5.2. Validation of the Interpollated Model
Previously, it was mentioned that the force-displacement characteristics of MRE under impact loading for input currents ranging from 0-0.5, 0.5-1, 1-1.5, and 1.5-2 Amperes were predicted using an interpolation algorithm. This algorithm is commonly used to estimate data points within a range of values based on known data points. To validate the accuracy of this interpolation algorithm, experimental works were conducted using drop impact tests. These tests were carried out with the input currents injected to MRE coils of 0.3, 0.7, 1.3, and 1.7 Amperes and the resulting force-displacement data were recorded. The experimental data obtained from the drop impact tests were then compared to the corresponding model response with the same input current as shown in
Figure 11. It can be seen from the figure that the interpollation algorithm is valid. The model responses closely follow the experimental results with acceptable error. The maximum error of predicted force of the proposed model in the interpollation regions is summarized in
Table 4.
5.3. Effects of Varying Swarm Size and the Number of Iterations
The number of iterations is an important parameter in PSO, as it determines the length of time the particles are allowed to search for the optimal solution. In general, increasing the number of iterations in PSO can improve the algorithm's ability to find the global optimum, as it allows the particles more time to explore the search space and converge on the best solution. However, at a certain point, additional iterations may not lead to any further improvement in the solution, as the particles may have already converged to a local optimum. On the other hand, decreasing the number of iterations can lead to faster execution times, but at the cost of potentially missing out on better solutions. Therefore, it is important to strike a balance between the number of iterations and the desired level of performance index and execution time. Referring to Figure 12, the number of iteration selected was 100 since it has fast convergence rate and able to achieve the performance index as good as the number of iteration 120.
Figure 12.
Effect of Varying the Number of Iteration to the Performance of PSO.
Figure 12.
Effect of Varying the Number of Iteration to the Performance of PSO.
The swarm size refers to the number of particles in the population, and it determines the diversity and convergence rate of the swarm. Increasing the swarm size can lead to better global exploration, as there are more particles searching the solution space. However, it can also increase the computational cost and reduce the convergence rate, as there are more particles to communicate and update. On the other hand, reducing the swarm size can improve the convergence rate, as there are fewer particles to communicate and update. However, it may also decrease the diversity of the swarm, which can lead to premature convergence and suboptimal solutions. Referring to Figure 13, the swarm size selected in this study was 40 as it shows fast convergence rate and achieved the performance index as good as the swarm size 100.
Figure 12.
Effect of Varying the Number of Particles to the Performance of PSO.
Figure 12.
Effect of Varying the Number of Particles to the Performance of PSO.
6. Conclusions
This paper presents a comprehensive study on the hysteresis behavior modeling of magnetorheological elastomers under impact loadings. The study proposes a multi-layered exponential-based Preisach model that is enhanced with particle swarm optimization to provide a reliable and accurate framework for capturing the complex hysteresis behavior of the material. The developed model demonstrates excellent performance in capturing the dynamic response of magnetorheological elastomers under various impact loading scenarios. The results show that the model's response closely matches the experimental data, with a maximum prediction error of 10% or less. The interpolated model's response also shows good agreement with the experimental data, with a maximum percentage error of less than 8.5%. The study also examines the effects of varying the number of iterations and the number of particles on the performance of PSO. Overall, the findings suggest that the proposed model provides a promising approach for accurately predicting the hysteresis behavior of magnetorheological elastomers under impact loadings.
Author Contributions
Conceptualization: Hudha and Kadir. Design / Fabrication: Alwi. Testing: Alwi and Kadir. Data Analysis: Alwi and Amer. Matlab Programming: Alwi and Kadir. Optimization: Alwi and Amer. Paper writing: Alwi and Hudha. .
Funding
This research is funded by Ministry of Higher Education (MOHE) of Malaysia.
Acknowledgments
This research is fully supported by FRGS grant, FRGS/1/2021/TK02/UPNM/02/1 led by A/Prof. Dr. K. Hudha. The authors fully acknowledged Ministry of Higher Education (MOHE) and Nasional Defence University of Malaysia for the approved fund which makes this important research viable and effective.
Conflicts of Interest
The authors declare that there is no conflict of interests.
References
- Y. Li, J. Li, W. Li, and H. Du, “A state-of-the-art review on magnetorheological elastomer devices,” Smart Mater. Struct. 2014, 23, 12. [CrossRef]
- A. Masa’id, B. W. Lenggana, U. Ubaidillah, D. D. Susilo, and S.-B. Choi, “A Review on Vibration Control Strategies Using Magnetorheological Materials Actuators: Application Perspective,” Actuators. 2023, 12, 3, 113. [CrossRef]
- R. P. Li, C. Bin Du, F. Guo, G. J. Yu, and X. G. Lin, “Performance of Variable Negative Stiffness MRE Vibration Isolation System,” Adv. Mater. Sci. Eng. 2015. [CrossRef]
- S. Sun et al., “A magnetorheological elastomer rail damper for wideband attenuation of rail noise and vibration,” J. Intell. Mater. Syst. Struct. 2020, 31, 2, 220–228. [CrossRef]
- N. Bira, P. Dhagat, and J. R. Davidson, “A Review of Magnetic Elastomers and Their Role in Soft Robotics,” Front. Robot. 2020, 7 , 1–9. [CrossRef]
- Ubaidillah, S. A. Mazlan, J. Sutrisno, and H. Zamzuri, “Potential applications of magnetorheological elastomers,” Appl. Mech. Mater. 2014, 663, 695–699. [CrossRef]
- R. Brancati, G. Di Massa, S. Pagano, and S. Strano, “Seismic isolation system for lightweight structures based on MRE devices,” Lect. Notes Eng. Comput. Sci. 2015 , 2218, 1313–1318.
- P. Skalski and K. Kalita, “Implementation of Magnetorheological Elastomers in Transport,” Trans. Inst. Aviat. 2016, 245, 4, 189–198. [CrossRef]
- A. Munteanu et al., “Reprocessed magnetorheological elastomers with reduced carbon footprint and their piezoresistive properties,” Sci. Rep. 2022, 12, 1, 1–12. [CrossRef]
- A. Z. Bin Pokaad, K. Hudha, M. Z. B. M. Nasir, and Ubaidillah, “Simulation and experimental studies on the behaviour of a magnetorheological damper under impact loading,” Int. J. Struct. Eng.2011, 2, 164–187. [CrossRef]
- H. W. He, Z. H. Wang, Z. J. Zou, and Y. Liu, “Nonparametric modeling of ship maneuvering motion based on self-designed fully connected neural network,” Ocean Eng. 2022 ,251,111113. [CrossRef]
- S. Mahmoud Taheri and G. Hesamian, “Non-parametric statistical tests for fuzzy observations: Fuzzy test statistic approach,” Int. J. Fuzzy Log. Intell. Syst. 2017, 17, 3, 145–153. [CrossRef]
- J. T. Thorson and I. G. Taylor, “A comparison of parametric, semi-parametric, and non-parametric approaches to selectivity in age-structured assessment models,” Fish. 2014 ,158, 74–83. [CrossRef]
- R. Brancati, G. Di Massa, M. Di Vaio, S. Pagano, and S. Santini, Experimental investigation on magneto-rheological elastomers, Springer International Publishing.2019, 68.
- K. Hudha, N. Salwa Sobri, K. Sumasundram, N. Akmal Haniffah, Z. Abd Kadir, and M. Sabirin Rahmat, “Investigation on the Effect of the Ferrous Particles Size on the Impact Absorption Capability of Magnetorheological Elastomer,” 2022 IEEE 18th Int. Colloq. Signal Process. Appl. 2022 , 299–303.
- M. M. Kuczmann, “Dynamic Preisach hysteresis model,” J. Adv. Res. Phys. 2010, 1, 1–5.
- C. Natale, F. Velardi, and C. Visone, “Identification and compensation of Preisach hysteresis models for magnetostrictive actuators,” Phys. B Condens. Matter. 2001, 306, 1–4, 161–165. [CrossRef]
- S. Willerich, “on Vectorial Relay Operators,” vol. 56, no. 3, pp. 1–4, 2020.
- P. B. Nguyen, S. B. Choi, and B. K. Song, “Development of a novel diagonal-weighted Preisach model for rate-independent hysteresis,” Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2017, 231, 5, 961–976. [CrossRef]
- F. Marini and B. Walczak, “Particle swarm optimization (PSO). A tutorial,” Chemom. Intell. Lab. Syst. 2015,149, 153–165. [CrossRef]
- C. L. Chan and C. L. Chen, “A cautious PSO with conditional random,” Expert Syst. Appl., 2015, 42, 8, 4120–4125. [CrossRef]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).