Adaptive Damping of a Double-Beam Structure Based on Magnetorheological Elastomer

A method of vibration reduction based on activation of an MRE block that couples twin cantilever beams at their free ends is investigated. Four types of magnetorheological elastomers have been manufactured, and their rheological properties in a range of magnetic field intensities are established. Free vibrations of several double-beam structures with controllable damping members made of these MREs are investigated, and a method of semiactive control of such structures is proposed. -e effects of compression of the elastomers and alignment of the magnets used to activate them are reported. -e mathematical modeling of the system is verified experimentally.


Introduction
Systems of two parallel elastic continua are used in various devices and mechanical, civil, or aircraft structures.Examples include aircraft wing spars, bridge spans, and linear guideways used in plotters or in some technological processes.In these systems, a soft viscoelastic layer or patch may be placed between two rigid facesheets to ensure efficient vibration suppression.
Basically, there are two methods to dissipate mechanical energy in such systems.When the initial conditions imposed on both beams are not the same or when the beams have different masses (stiffnesses), one may exploit the relative transverse motion of the beams, which leads to an axial deformation of the viscoelastic element placed between them.However, the beams eventually synchronize their motion, and this way of damping becomes ineffective-the same as in the case of identical beams and initial conditions [1].en, the shear deformation of the soft middle layer becomes the predominant mechanism of energy dissipation in a sandwich beam [2][3][4][5][6][7].
In order to control these vibrations, smart materials can be used in the soft middle layer.Magnetorheological elastomers (MREs) are composite materials that are considered to be smart materials [8].Solid-state MREs overcome drawbacks of MR fluids, such as the leakage of the liquid and sedimentation of the particles.ey are classified as intelligent materials well suited for providing effective structural control.ey possess tunable rheological properties and can be used in adaptive structures for aerospace [9], automotive [10,11], railway [12,13], civil [14], and electrical engineering applications for linear and rotational motion and vibration control [15,16].
Magnetorheological elastomers exhibit interesting mechanical properties.ese materials change their rheological properties under an applied external magnetic field.is behavior is similar to magnetorheological fluids, but MREs are their solid equivalents.In magnetorheological elastomers, magnetoactive particles (e.g., iron and cobalt) are embedded within an elastic polymer matrix [17].e matrix enables relatively high strains and maintains the overall shape of the element.e magnetorheological effect is based on dipole interactions between ferromagnetic particles [18].Under a magnetic field, the particles tend to align themselves in parallel with the lines of the magnetic field [19].Due to their rapid response [20,21], the elastic moduli of these flexible smart composites can be effectively controlled in almost real-time according to varying external excitations.us, the MREs provide greater potential for vibration suppression in an active manner in a wider range of frequencies [22][23][24].
is tendency has been used to obtain several magnetorheological elastomers with aligned chain-like microstructures.e polymer matrix preserves the alignment of the particles after a curing process.From our previous research [25], it is known that in MR composites, an aligned microstructure provides a significantly greater magnetorheological effect than an isotropic arrangement of the particles [26,27].To manufacture the magnetorheological elastomers considered in our investigations described in the paper, we followed the methods of Boczkowska et al. [28].Applying this method, it has been found that the rheological properties strongly depend on the MRE microstructure, in particular on the size/shape of the particles and their arrangement.By optimizing the particles' size, shape, and alignment, the stiffness of the MREs under an applied magnetic field has been increased.
e application of magnetorheological materials to the mitigation of the vibrations of sandwich beams is a relatively new idea.ese materials can be placed along the entire length of the system in the form of an adaptive viscoelastic core [29,30], as a single patch [31], or as several patches [32].
e dynamics of these systems when subjected to various excitations has been investigated quite thoroughly, both experimentally and theoretically [33,34].
e aforementioned studies are related to passive methods of vibration reduction of double-beam structures, which are not always sufficiently effective.
In addition to beam structures, MREs can be used in plate systems.Yeh [35] presented a vibration analysis of a rectangular sandwich plate composed of an MRE core layer and a constraining layer.
e MR material exhibited a significant effect on the vibrational characteristics of the plate and a rapid time response.Vemuluri et al. [24] studied the dynamic performance of a tapered sandwich plate partially treated with an MRE layer.e laminated composite plate with an MR fluid was investigated by Ramamoorthy et al. [36].de Souza Eloy et al. [22] investigated numerically and experimentally sandwich panels with carbon/epoxy composite skins and a magnetorheological elastomer (MRE) honeycomb core.
e obtained experimental results showed good performance in the attenuation of the vibration level, especially the fundamental vibration mode of the structure.e magnetic field applied on the free end of the test panel was able to reduce significantly its first natural frequency value.
Few of these works tackle the problem of the semiactive control of such systems using magnetorheological materials.Ndemanou et al. independently proposed a method based on the Lyapunov stability theory [37,38].Namely, the controller aims at the maximization of the instantaneous energy dissipation of the system.In the first of these works, a double-beam structure is stabilized through the relative transverse motion of the beams, which leads to an axial deformation of a controllable MR damper.In the second, the shear deformation of several damping blocks made of an MRE is exploited.In another approach based on the latter mechanism, the linear quadratic regulator was applied [39].
ese studies are confined to numerical simulations of semiactively controlled systems.For example, Rajamohan et al. claim that "the implementation of the controller, however, would require further developments in compact electromagnets to apply the controlled magnetic field on the surfaces of the treated structures" [39].Indeed, in most of the studies, the electromagnetic devices were placed near the vibrating structure to control the rheological properties of the material [3,6].However, such an approach requires large and heavy machinery built around the system, which may be cumbersome.Moreover, attracting the MRE by external magnets causes the effect of negative stiffness, which is in general undesirable.As far as the knowledge of the present authors extends, the research [40] is the only one in which the electromagnets were attached directly to the vibrating structure.is study showed a rather mediocre performance of the control in the system with classical electromagnets and a damping block made of an isotropic MRE. e rest of this paper is structured as follows.Four types of magnetorheological elastomers of have been manufactured.eir rheological properties in a range of magnetic field intensities have been established, and the free vibrations of several double-beam structures with damping blocks made of these materials are investigated.Permanent magnets are then used to identify a type of MRE, providing decent damping capabilities.Apart from the increased shear stiffness and damping, additional phenomena that affect the performance of the method are reported.A dynamic model of the system is presented and used to validate a stabilizing semiactive control strategy.Finally, the theoretical findings are verified experimentally.

Fabrication and Rheological Properties of the MRE Samples
Two types of elastomer matrices were used to produce the MRE samples, referred to as PU 70/30 and U1440.e PU 70/30 is a soft polyether-based polyurethane that is obtained from two polyether polyols: 14922 and VORALUX ® HF 505 (in weight ratio 70 : 30) in a mixture with the isocyanate compound HB 6013, supplied by the Dow Chemical Company.e second matrix, referred to as U1440, was synthesized using the two-component polyurethane elastomer Biresin ® U1404, supplied by SIKA.e elastomer was produced by a reaction of isocyanate prepolymer resin and amine hardener in 80 : 100 weight ratio.
In both cases, the mixing of the substrates and the curing process were conducted at room temperature.e PU 70/30 polyurethane is characterized by a low density (1.03 g/cm 3 ), low viscosity before curing (ca.1600 mPa•s), and low hardness (below 10 °ShA).e U1440 polyurethane is also characterized by a low density (1.05 g/cm 3 ), relatively low viscosity before curing (ca.3000 mPa•s), and a hardness (40 °ShA) four times higher than that of PU 70/30.ese properties allow us to obtain a relatively high magnetorheological effect in the composites based on both types of 2 Shock and Vibration matrices.e relatively low viscosity during the processing of the MRE makes it easy to obtain an alignment of the particles under an external magnetic field.Carbonyl iron powder, produced by Fluka, was used for the magnetoactive particles for the fabrication of the samples.e particles are spherical, and their average diameter is 6-9 μm.Two volume shares of the carbonyl iron particles were used.e first was 11.5 vol.%.In our previous studies [41], we obtained the highest magnetorheological effect for this amount of particles.e second amount was equal to 33 vol.%. e samples were subjected to a magnetic field during curing to obtain aligned iron particle chains within the elastomer.A magnetic field strength of 240 kA/m was used.
e characteristics of all the MRE samples are given in Table 1.
e rheological measurements were carried out using an ARES rheometer from TA Instruments under a preselected magnetic field (Figure 1).e plate-plate system with diameters of 20 mm and with a 2 mm gap between plates was used.e magnetic field was varied over the range 0-320 kA/ m. e experiments were conducted at an ambient temperature of 25 °C.e excitation was applied in a sinusoidal form with angular frequency ω � 19.9 rad/s and at constant strain of 0.1%.
MR elastomers are normally operated with small deformations in the preyield regime of the linear viscoelastic region.ey are intended to be used as structural materials in applications where the load is often dynamic.Under a cyclic dynamic loading, the material deforms and returns to its original form during one cycle.In viscoelastic materials, some part of the deformation energy input is stored and recovered during each cycle and some is dissipated as heat.
e storage modulus, G ′ , represents the ability of the vis- coelastic material to store the energy of deformation, which contributes to the material's stiffness.e loss modulus, G ″ , represents the ability of the material to dissipate the energy of deformation.
In this study, both the elastic (storage) modulus and the loss modulus were measured as functions of the magnetic field at room temperature.As shown in Figure 2, the storage modulus increases significantly with the strength of the applied magnetic field for all samples.However, this growth strongly depends on the matrix material and the amount of particles within the matrix.
e change of the modulus under a magnetic field is called the magnetorheological effect.e magnetorheological effect can be expressed by the relative change of the elastic modulus ΔG ′ /G 0 ′ under the magnetic field.is effect can be used for variable stiffness devices because it is fully reversible.It is related to the tendency of the composites to attain the minimal energy.
e soft matrix enables some elastic displacement of the particles.However, such displacements generate a strain in the matrix and a stress.As a result, the material's properties are changed.e softer matrix PU 70/30 exhibited a larger MR effect than U1404 due to larger strengthening effect of the particles under the magnetic field.
e same conclusions regarding the behavior of the elastomers under the magnetic field can be drawn for the loss modulus (see Figure 3).
During our experiments with a passively activated double-beam structure, the values of the magnetic field measured in the centers of damping elements amounted to 103, 159, and 263 kA/m, respectively, with all three types of magnets used in the study, and without the brass plate.e brass plate was intended to prevent compression of the elastomers (see Section 5.4).When the brass plate was inserted, the intensity dropped by approximately 8 kA/m.

Dynamic Model
is section describes the mathematical model of the system which will by used for development of the control strategy.It is based on the model already presented and verified experimentally [42], which was further simplified in [1].
e double-beam structure to be investigated is presented in Figure 4.It consists of two elastic facesheets coupled at their free ends by a smart damping element made of MRE.To prevent buckling of the structure, a distance roll modeled as a linear spring of a very high stiffness k is placed between the beams at a distance ξ k .Both beams have length L, width b, thickness d, and mass per unit length ρS and are mounted in a clamped support.
We will study the vibrations of the structure in the plane of symmetry ξ − η. e displacements of the beams from their equilibria in the η-direction are denoted by w 1,2 (ξ, t).
e beams are slender, and their deflections are small, so the linear Bernoulli-Euler beam theory is valid.
e external damping of the beams is neglected.e effect of gravity, which acts perpendicularly to the plane of vibrations, is also not taken into account.
e damping block has a rectangular shape 2a × b × 2h and mass 2 M. Its center is located at a distance of ξ � L − a from the support.e rheological parameters of this element are controlled using electromagnetic devices, which are attached at the tip of the structure and have a total mass of M a .e constitutive law of the MRE in a dynamic test can be written in the following form: where τ is the shear stress, φ � φ(t) is the angle of shear deformation, and G ′ � G ′ (u, ω) and G ″ � G ″ (u, ω) are the storage and loss moduli.Both these parameters are functions of the applied magnetic field u and the frequency of angular deformation ω.However, a typical change of storage or loss modulus due to the applied magnetic field yields only a slight change in this frequency.e values of G ′ and G ″ are nearly constant within such a narrow range of ω, so we may treat the values of both rheological parameters as independent of it.Shock and Vibration Moreover, the external devices are only able to control the rheological parameters in a limited range, where linear approximations are justi ed.us, we can write where u ∈ [0, 1] is the normalized value of the magnetic eld intensity, G 0 ′ and G ″ 0 are the storage and loss modulus without any magnetic eld applied, and r 1 and r 2 are positive constants determined experimentally.
e shear deformation of the damping element yields bending moments as well as axial forces acting on both beams.Moreover, the relative motion of the beams generates transverse forces.Having assumed that the block is not too long (2a ≪ L) and treating the parts of the beams over the damping block as rigid bodies, we may reduce the forces to the pole, which is the point over the center of the block (see Figure 5).
We now present the formulas for the bending moments and forces acting on Beam 1, which have been developed using the truss model proposed in [42].Assume that the angles of de ection of the beams at the point ξ L − a are the same and equal to φ (zw/zξ)(L − a, t) and the relative displacements of the beams at this point are Δw w 1 (L − a, t) − w 2 (L − a, t).We then have the following expressions.
For the bending moment M, and axial forces N, e in uence of the axial force appears to be of minor importance and will be disregarded.is is mainly because the axial forces act in opposite directions (i.e., when Beam 1 is compressed, Beam 2 is stretched) and partly because these forces are too small to disturb the dynamics of the rigid facesheets.With this simpli cation, we end up with a linear dynamical model in the form of the following system of coupled partial di erential equations governing the dynamics of both beams: where M and F are given by formulas ( 3) and ( 4), δ L−a is Dirac's delta at ξ L − a, and (d/dξ)δ L−a is the unit doublet at ξ L − a.Note that the positive sign of the transverse force F corresponds to beam 1, and the negative sign corresponds to the other beam.
To develop a control strategy, we need to introduce an energy-related function.
e most natural choice for this function, which also turns out to be e ective, is to de ne V as a function related to the total mechanical energy of the system, similar to [1].We neglect the kinetic energy of rotation of the damping element and the kinetic energy of its shear deformation because they make only a minor contribution to the total mechanical energy of the system.After some calculations, we obtain the following formulas: (1) e combined kinetic energy of the damping element and the actuator: (2) e elastic strain energy of the damping element (note that (zw 1 /zξ)(L − a, t) (zw 2 /zξ)(L − a, t) by assumption): (3) e kinetic energy of the beams:

Shock and Vibration
(4) e elastic strain energy of the beams: erefore, the energy function, closely related to the total mechanical energy of the system, is given by the following equation: We are going to analyze the approximate solution to the partial equations governing the system in the form of a linear combination of cantilever beam eigenfunctions: where λ i are the roots of the characteristic equation cos λ cosh λ � −1, and C i are the constants which normalize the mode shapes with respect to the norm induced by the scalar product  L 0 W i (ξ)W j (ξ)dξ.Applying Galerkin's procedure, we obtain a set of 2n second-order ordinary equations for the vectors of generalized coordinates T 1 � [T 1  1 , . . ., T 1 n ] T and T 2 � [T 2 1 , . . ., T 2 n ] T .We then reduce the order of the problem by substituting is leads to the following normal form of the set of 4n ordinary differential equations governing the dynamics: For the set of admissible controls, we assume Following [42], eight basis functions appear to be sufficient to obtain the desired precision for the approximation.
e energy function (11) for the discrete system ( 14) can be written as a quadratic function of the state variables as follows: where the energy matrix is a function of the control variable, i.e., Q � Q(u).However, we have found that Q can be treated as independent of u. us, we assume

Semiactive Control Strategy
We will now study the stabilization of the semiactive system (14) with respect to perturbations of the initial state x 0 .e problem of stabilization consists in choosing a feedback control function u ∈ U such that the solution x of the resulting feedback system associated to (14) converges in some sense to the equilibrium state x 0 within a suitable length of time.For the controller, we shall employ the method proposed in [1,43].e selection of this method is motivated by the good balance it strikes between the performance of the controller and its practical simplicity.It has been validated numerically that the method achieves a fair performance in stabilizing impulse-excited vibrations.By relying on a time-invariant state-feedback switching control law, the method can be successfully implemented using a system of displacement sensors, a signal amplifier, and a PC equipped with an analog I/O card.For the convenience of the reader, we summarize the underlying principles of the method in the following.
e aim is to design a control u * (t) ∈ U that provides the best instantaneous dissipation of the energy V as defined in (16).erefore, we require the following equation: Note that, By using Q � Q T , the insertion of the governing equations ( 14)-( 18) yields From ( 15) and (19), it follows that the solution to the problem ( 17) is given by the following equation: It should be noted that the internal structural damping guarantees that A is a Hurwitz matrix (i.e., every eigenvalue of A has strictly negative real part).Assuming the following bounds on the admissible control u min � 0, u max > 0, we can then see from Q≻0 and (19) that the control law given by (20) guarantees a permanent decrease of the global system energy _ V < 0, and thus the asymptotic stability of (14).

Test Stand.
As the real object of investigation, we consider a two-beam structure imitating structural elements such as aircraft wing spars, cranes, bridge spans, pipelines, and rail plotters.A photograph of the entire laboratory setup is shown in Figure 6, a schematic in Figure 4, and the fundamental parameters of this test rig in Table 2. e system was made of steel beams, each having nominal length L � 0.8 m and cross-sectional area S � 0.00003 m 2 (0.015 m × 0.002 m). e steel material was chosen because of its ferromagnetic properties, which 6 Shock and Vibration allowed us to obtain the highest magnetic field intensity.e free ends of the beams were coupled by a rectangular block made of one of the MREs that were described in Section 2. e damping element was glued using an MS-polymer adhesive.
e assembly of the MRE element in the investigated system allowed for an easy change of the magnetic flux by using various permanent magnets attached at the end of the structure.e basic parameters for the NdFeB magnets used in the investigation are presented in Table 3.
e magnets were stuck by the magnetic force generated between them and the facesheets.One of the magnets faced its north pole towards the beam while the polarization of the other one was the reverse.Consequently, the magnetic field within the damping element was as uniform as possible, and the passive action of the magnets increased both the shear stiffness and damping properties of the elastomers. is allowed a passive-adaptive damping of the structure.
e shear deformation of the damping member generates axial forces, which may produce buckling of the beam that is compressed. is complicates the measurements and analyses of the effect in which we are actually interested.To keep the gap between the beams constant over their lengths, two lightweight rolls were placed at distances of x 1 � 0.247 m and x 2 � 0.546 m from the support.e beams were tightened to the distance rolls by thin and light rubber bands.e rolls were made of aluminum tubes, so the rolling drag would be as small as possible.
e structure was fixed to a solid workshop table in a clamped configuration with two steel mounting brackets.
e oscillatory motion of the coupled two-beam system was realized by a lateral deflection in the (horizontal) y-direction.
e plane of motion was perpendicular to the force of gravity, thus the effect of gravity could be neglected.To arrive at a precise measurement of the horizontal deflection, Baumer-OADM 12I6430/S35A laser distance sensors were used.e laser sensors provided a maximum measurement resolution of 0.01 mm in less than 900 μs.Because the measurement range of the laser sensors was within the 30-130 mm interval, they were mounted near the beams at three points, marked by S 1 , S 2 , and S 3 in Figure 4.Because these sensors realized a contactless and perturbation-free signal registration, their operation did not reduce the accuracy of the remaining measurement systems in the considered setup, as confirmed in [44,45].During the tests, an NI USB-6218 analog-digital converter was used for data acquisition.All of the signals from the laser sensors were recorded by means of this converter with a sampling frequency of 1000 Hz.
e experiments were conducted according to the scheme (1A, 1B, 2A, 2B) × (thick, thin) × (no magnet, M5, M8, M15) × (compression, no compression) × (trial 1, trial 2), which resulted in 128 experimental scenarios.For the thin beam (4 mm thickness of the damping member), the results were qualitatively the same as for the thick beam (2h � 6 mm), yet quantitatively in most scenarios they are less pronounced and, therefore, they will be disregarded in this article.However, when analyzing the effect of the alignment of the magnets, the results obtained for the thin beam are easier to see, so some of them are going to be presented in Section 5.5.
e effect of the compression of the elastomer by the forces generated by the magnets is inherently present in the

Shock and Vibration
system, and for the purpose of analysis, it was excluded by inserting a brass distance plate near the damping block (see Section 5.4 for details).Two trials were performed to ensure that the results would be repeatable.Apart from a few trials when random accidents happened (such as a detachment of the sticker at which the laser sensor was aimed), the experiments were actually perfectly repeatable and, therefore, only the results for single trials will be presented.
Each trial was initiated by a de ection of the system with zero initial velocity.After releasing the structure, the vanishing free vibrations of the rst mode were recorded.e excessive amount of data was reduced by taking averages of the subsequent blocks of ve observations.Examples of time series for a trial for the thick beam, with damping member made of elastomer 1B, under the action of magnets M8, and with the e ect of compression (no brass plate), are shown in Figure 7.

Experimentally Observed Vibration Period and Damping.
e experiments with free vibrations will be analyzed in this section.e oscillations at the rst natural frequency were induced.e consecutive minimal and maximal beam de ections w i and the points t i of time at which they occur were determined fully automatically by local approximations of the time series with quartic polynomials.e extremes were calculated until the amplitude of the vibrations decreased to the threshold when the signal-tonoise ratio became too low, which happened when the amplitude measured by sensor S 3 was 2 mm.An example of the time response with calculated extremes is shown in Figure 8.
To investigate the character of the vibrations, doubled times of increase and decrease τ i and logarithmic damping decrements δ i corresponding to consecutive extreme deections were calculated: No di erences were observed between the times of decrease and the times of increase; so, τ i is the estimator of the period of the vibrations.e values of the period and the decrement for consecutive cycles are shown in Figures 9 and 10.
For a classical harmonic oscillator with viscous damping, the series τ i and δ i are constant.us, the amplitude of the vibrations decreases at a constant rate, making the oscillations vanish exponentially.In the analyzed system, this is satis ed quite well in most of the analyzed scenarios.In the case of the beam with a damping member made of MRE 1B and subjected to the action of the biggest magnets, a nonlinear behavior of the material can be observed because the vibrations of high amplitudes have longer periods, and their suppression is less pronounced.e attachment of magnets generally increases the period of vibrations because of the additional inertia introduced to the system.Moreover, the damping becomes more pronounced.To analyze these e ects, an analysis analogous to the one from the previous subsection was conducted.With no magnets, the average period of vibrations is 0.390, 0.373, 0.351, and 0.324, respectively, for MRE 1A, 1B, 2A, and 2B.e logarithmic damping decrement is 0.038, 0.108, 0.067, and 0.123.e damping ratio (the quotient of the decrement and the period) is 0.097, 0.289, 0.190, and 0.378 (1/s).ese values are much higher than those theoretically predicted because structural damping is present in the real system, yet it is neglected in the model.
We will now move on to the empirical analysis of the magnetorheological e ect on vibrations (Figure 11).e action of the magnets increases the sti ness of the damping members.
e relative growth of the average vibration period is both qualitatively and quantitatively in good agreement with the theoretical predictions.For the damping (both the average decrement and the damping factor), the results are a little di erent.e relative growth of damping under the action of the magnets is lower than predicted by the model.In the real system, the energy dissipates not only through the damping member but also inside the beams, through the support, and in the form of aerodynamic drag and rolling friction of the distance rolls.
is structural and material damping is di cult to quantify, and it partly overshadows the magnetorheological e ect.However, it yields results that are more realistic because many sources of damping are inherently present in all engineering structures.8 Shock and Vibration e relative growth of the damping is much higher for both elastomers based on the PU 70/30 matrix, which con rms our theoretical predictions.Note that for both elastomers based on the U1440 matrix, the relative damping factor is almost constant, which means that within a given time period, the amplitude of vibrations vanishes by a similar amount either with the magnets or without them, so the damping turns out to be ine ective.It is only the character of the oscillations that changes: the cycles are longer and their decrement is higher.

E ect of the Attached Mass.
e additional mass of the magnets a ects the dynamics of the system.e period of vibrations increases and the damping decrement is lower.On the other hand, the magnetic eld generated by the magnets exerts the opposite e ect, because it sti ens the damping block and increases its loss modulus.To analyze these e ects separately, we apply the proposed mathematical model with the parameters the same as in experiments.e damping element is assumed to be made of MRE 1A which turned out to be the most sensitive to changing the magnetic eld of all of the tested materials.Figure 12 compares the numerically obtained tip de ections of the beam without the magnet, with an arti cial mass equal to the mass of the magnet and with the 5 mm neodymium magnet.Note that the strong eld of the permanent magnet overcomes the negative e ect of its additional mass.Still, such a positive relation of the magnetic eld intensity to the mass is di cult to be obtained using electromagnets.

Compression of the Elastomers.
When magnets are attached to the structure, they pull each other and squeeze the damping element.
is introduces an additional stress, which may be expected to have an in uence on the dynamics of the system.To analyze this e ect, we inserted a plate at one side of the damping block (see Figure 13).Brass was chosen as the material of the plate because it does not modify the distribution of the magnetic eld passing through the structure.
A comparison is shown in Figure 14. e compression of the elastomers increased their shear sti ness because the vibration period of the structure was shorter for the system without the plate than with the plate.is e ect was observed in all of the analyzed scenarios.As could be expected, it was most pronounced for softest, MRE 1A.

Alignment of the Magnets.
e magnets align themselves along the lines of the magnetic eld to achieve the con guration of the lowest possible energy (Figure 15). is Shock and Vibration introduces an additional sti ness to the system.To analyze this e ect, two brass rolls were placed between the M15 magnet and the beam, in two directions (Figure 16).When the rolls were inserted vertically, the magnet could align its position with the magnetic eld.When they were in the horizontal position, this motion was prevented.

Shock and Vibration 11
Although the motion of the magnet was hardly visible-which is understandable given that the shear strains were very low-the e ect of rolling was observable, especially for the thin double-beam structure, in which the magnetic force between the magnets was higher than for the 6 mm thick structure (Figure 17).

Experimental Assessment of the Semiactive Control Strategy
6.1.Test Stand.To determine the additional functionality of the considered MRC blocks, associated with the increase in damping e ciency by applying a semiactive control strategy, an experimental investigation was carried out.For this purpose, the test stand described previously was modi ed.
Similarly to the previous test rig (Figure 6), the ends of the two steel beams were supported with mounting brackets, and the MRE 1A element was glued at the free ends of the beams (see Figure 17).e MRE 1A block was selected because of all of the tested materials, it was the softest and the most sensitive to changing the magnetic eld. is support enables horizontal displacements of the beams during the tests.In this investigation, the horizontal displacements were measured by the same laser displacement sensors OADM 12I6430/S35A.Here, four sensors in two pairs were set (see Figure 18).Each pair was used to estimate the horizontal displacement of the measuring point located between the two sensors.For this estimate, we used the average value measured by the two sensors, which allows an e cient signal ltering.
An electromagnet (GTo-30-0.5000-24VDC)was attached near the MRE 1A element, at the end of the beam structure.
is electromagnet was mounted using the MS-polymer adhesive.To direct the magnetic ux beam generated by the electromagnet perpendicularly to the MRE block, on the opposite side of the steel beam, a neodymium magnet with dimensions of 10 × 10 × 3 mm was used.e electromagnet was connected with the control unit by a current source, generating a control current proportional to a voltage signal that regulated a magnetic ux surrounded by the MRE block.e use of the current source allowed us to increase the rate of change of the magnetic ux generated by the magnet.
e electromagnet was turned on and o by a computer equipped with the semiactive control algorithm (20).e MRE block semiactive control strategy algorithm was adapted in the MATLAB environment, which was used to communicate with the measurement and control devices.For the tests, a National Instrument PCI 6251 analog-digital converter was used for data acquisition and control.All of the signals from the laser sensors were recorded by means of this converter with a sampling frequency of 1000 Hz. ese signals provided information to the algorithm, which generated the control signal to the electromagnet with a sampling frequency of 50 Hz.
e scheme of the measuring/ control system is presented in Figure 19.e detailed parameters of the test stand (see Figure 18) used in the investigations are listed in Table 4.

Free Vibration Control
Test.We aim at demonstrating the general performance of the selected MRE 1A block in  12 Shock and Vibration mitigating the structural vibration.erefore, in the present work, rather than performing a comprehensive control analysis, we limit our attention to the most typical cases of free vibration, where the de ections of the cantilever structure are concerned with the shape of its rst natural mode W 1 , as de ned in (13).For this purpose, we carried out control experiments assuming the following initial state condition: w 1,2 (ξ, 0) 50W 1 (ξ)/W 1 (L) (mm).As in the foregoing experiments, we selected the observation time to be t f 8 (s).e state-feedback switched controller (20) relies on the measurements of two displacements w 1 (0.5(ξ s1 + ξ s2 ), t) and w 1 (0.5(ξ s3 + ξ s4 ), t) estimated by the two pairs of the laser sensors.ese two measurements have been used to reconstruct the generalized coordinates T 1 1 (t) and T 1 2 (t) by using ( 12) and (13).For simplicity, we assume w 2 (0.5(ξ s1 + ξ s2 ), t) w 1 (0.5(ξ s1 + ξ s2 ), t) and w 2 (0.5(ξ s3 + ξ s4 ), t) w 1 (0.5(ξ s3 + ξ s4 ), t), and consequently T 2 1 (t) T 1 1 (t) and T 2 2 (t) T 1 2 (t).In order to justify the e ectiveness of the semiactive controller, we will use the following functional as the stabilizing performance measure: where V is the total energy function de ned as in (11).e performance of the control u u * will be examined by comparison to the worst and the best passive cases, where we assume the extreme admissible control values u u min   (obtained for the voltage of 0 V) and u u max (obtained for the voltage of 24 V), respectively.e values of the objective functional (22) in the passive and the switched control cases are listed in Table 5.Each value is normalized to the passive case u u min .e comparison of the de ection of the end point of the lower beam w 1 (L, t) for the considered cases is depicted in Figure 20(a).Regarding the overall stabilizing performance, the reader can verify that the switched control (Figure 20(b)) provides better results than both passive damping strategies.In terms of the assumed functional (22), the control outperforms the worst and best passive cases by 13.2% and 7.6%, respectively.Comparing the last extreme de ections at t ≈ 7.3 (s), the control results in a reduction of the vibration amplitude of 12.8% and 9.2%, respectively.

Conclusions
e use of magnetorheological elastomers for the mitigation of the free vibrations of a double-beam structure was studied.Four MREs with two types of polymeric matrices and two volume shares of iron particles were manufactured and tested on a rheometer to assess their storage and loss moduli within a magnetic eld.
en, damping elements made of these materials were fabricated and glued between two cantilever beams, so that they coupled the free ends of the beams.e free vibrations of the rst mode were induced and analyzed, and a simple linear model of the structure was developed.Both theoretical and experimental studies showed that the performance of the investigated method depends markedly on the type of the material and the permanent magnets attached: (i) e PU 70/30 polymeric matrix is characterized by lower storage and loss moduli than the U1440 matrix, so its damping capabilities are not as good as the latter when no magnetic eld is applied.(ii) e magnetorheological e ect for the PU 70/30 polymeric matrix is much more visible than for the U1440 matrix, which allows easier control of the damper.When a su ciently strong magnetic eld is applied, the storage and loss moduli of MREs based on the rst type of matrix exceed those for the U1440-based elastomers.(iii) Magnetorheological elastomers having 33 vol.%Fe provide higher sti ness, damping, and magnetorheological e ects than those containing 11.5 vol.% of iron particles.(iv) MRE 1B (PU 70/30, 33 vol.%Fe) was the most e ective for vibration suppression of the investigated structure.We think that this is due to the unique microstructure and microstresses inside this material, which were obtained due to the selection of the properties of the components of this composite.(v) e compression of the elastomers under the action of the magnetic forces, as well as the e ect of the alignment of the magnets with the magnetic eld, increases the sti ness and damping of the damping members.(vi) Application of an elastomeric damping member at the end of the structure results in a sigmoidal shape of the rst mode, so additional elements are necessarily along the structure to prevent buckling deformation.(vii) e proposed model is very simple but it accurately captures the dynamics of the structure and allows development of the control strategy.
Data Availability e data used to support the ndings of this study are included within the article.Shock and Vibration

Figure 2 :Figure 3 :
Figure 2: Storage modulus G′ and relative MR e ect for samples under di erent magnetic eld strengths.

Figure 1 :
Figure 1: ARES rheometer with magnetic coil device and measurement scheme.

Figure 4 :
Figure 4: Double-beam structure with embedded intelligent damping MRE element.

Figure 5 :
Figure 5: Forces generated by the damping block reduced to the pole.

Figure 6 :
Figure 6: Laboratory two-beam construction system with MR elastomer damper.

Figure 15 :
Figure 15: Alignment of magnets with the magnetic eld.

Figure 16 :
Figure 16: e rolls placed in two directions to analyze the e ect of alignment of the magnets.

Figure 18 :Figure 19 :
Figure 18: e test stand for the investigation of the control applied to the MRE element.

Figure 20 :
Figure 20: e trajectories of the tip de ection in the passive and the switched control cases (a) and switched control signal (b).

Table 1 :
Characterization of fabricated MRE samples.

Table 2 :
Fundamental parameters of the test rig.

Table 3 :
Neodymium Iron Boron magnets used in the investigation.

Table 4 :
Fundamental parameters of the second test stand.
Parameter Value Parameters of beam material (dural) ρ 2.7 g/cm 3 , E 69 GPa, E * 40.3 MPa•s Dimensions of a single beam L 860, b 20, d 1.6 mm

Table 5 :
e values of the objective function in the passive and the switched control cases.