Next Article in Journal
Structural Low-Level Dynamic Response Analysis Using Deviations of Idealized Edge Profiles and Video Acceleration Magnification
Next Article in Special Issue
Suppression of the Hydrodynamic Noise Induced by the Horseshoe Vortex through Mechanical Vortex Generators
Previous Article in Journal
Ultrafiltration/Granulated Active Carbon-Biofilter: Efficient Removal of a Broad Range of Micropollutants
Previous Article in Special Issue
Development of Multi-Staged Adaptive Filtering Algorithm for Periodic Structure-Based Active Vibration Control System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantification of Active Structural Path for Vibration Reduction Control of Plate Structure under Sinusoidal Excitation

School of Mechanical Engineering, Yeungnam University, 280 Daehak-ro, Gyeongsan-si, Gyeongsangbuk-do 712-749, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(4), 711; https://doi.org/10.3390/app9040711
Submission received: 21 January 2019 / Revised: 12 February 2019 / Accepted: 13 February 2019 / Published: 18 February 2019
(This article belongs to the Special Issue Active and Passive Noise Control)

Abstract

:
The engines of electric and hybrid vehicles cause vibration and noise with complex frequency spectra. This tendency is observed especially for mid-frequency components. Therefore, there are limitations in achieving vibration attenuation using only structural changes, which is the conventional isolation method. A smart structure-based active engine mounting system is a core technology that can continuously improve the noise, vibration, and harshness performance under various operating conditions by continuously controlling the dynamic characteristics of the mount. It can replace the existing mount technology which supports the engine by realizing both static and dynamic stiffness. This study focuses on vibration reduction for a typical mid-frequency excitation. Based on a source–path–receiver structure, a mathematical model for a 3D plate structure with three active paths is proposed. The amplitude and phase of the actuator were calculated for the reduction of given vibrations on the basis of the model. When controlling with the proposed method, it was necessary to perform a large amount of computation and to newly define the modeling according to the structural change. To solve this inconvenience, the NLMS (normalized least mean squares) algorithm was applied, and the results were compared. It is shown that the application of the NLMS algorithm to perform the overall vibration reduction is more effective than the previous method.

1. Introduction

Recently, the downsizing and manufacturing of high-efficiency vehicles are being pursued to produce a lightweight vehicle body and improve engine performance. This causes problems of noise, vibration, and harshness (NVH) performance deterioration, furthered by the many restrictions on the existing method of vibration reduction. To solve these problems, there is a need for active engine mounting systems. These systems are attracting attention as a technology that can continuously control mount dynamics and improve the NVH performance under various operating conditions.
Many studies have been conducted on noise and vibration control using smart materials. Particularly, because the vibrations generated in a powertrain cause interior noise due to their structural transmission. Active mounts capable of attenuating transmitted vibrations have been developed, and this will be a focus in the future. Minak et al. [1] reviewed the overall benefits of designing a solar car. In particular, it was proven that optimizing the weight of the vehicle has a positive influence on the performance. Work by Hosseini et al. [2] presented a nonlinear modeling of solenoid actuators to be used as vehicle engine mounts. Experiments on dynamic stiffness changes have shown that nonlinear equations can be linearized in specific fields. Tang et al. [3] examined the torsional vibration characteristics of a compound planetary hybrid propulsion system using the torsional vibration dynamic model. It showed good tracking performance with low-frequency vibrations. Kraus et al. [4,5] manufactured and tested an engine mount by connecting an actuator and a viscous damper in series. Using the FX-LMS algorithm, the actuator input signal was generated from the mass, and the experiment showed a reduction in noise and vibration. Hausberg et al. [6] theoretically and experimentally analyzed the nonlinearity, component temperature, and preload change of the secondary path change of active engine mounts (AEMs). The dynamic characteristics of AEMs were predicted using the FX-LMS algorithm. Yang et al. [7] applied two adaptive multi-notch filtered algorithms through four hydraulic actuators to isolate engine vibrations. However, the harmonic component increased due to nonlinearities. Sun et al. [8] suggested the in-wheel switched reluctance motor model using the Fourier fitting method and analyzed the unbalanced residual lateral force. To solve the vibration problem, a control method was proposed and showed that vibration could be effectively reduced.
Chae et al. [9] used a magnetorheological damper to reduce the vibration caused by vertical, rolling, and pitching operations. Jeon et al. [10] developed a new controllable engine mounting system used for vibration control. For new controllable engine mounting systems, right-hand mounts were produced using piezoelectric stack actuators, and roll mounts were produced using magnetorheological fluid actuators. Fakhari et al. [11] applied a robust model reference adaptive control (MRAC) to control vibrations using an electromagnetic engine mount. MRAC was applied to determine vibration isolation performance through experiments and simulations. Elahinia et al. [12] proposed a new semi-active actuator for shock and vibration isolation using magnetorheological and electrorheological fluids. Wu et al. [13] developed a vibration insulator using a magnetic spring with negative stiffness (MS-NS) with high-static–low-dynamic stiffness characteristics and showed that it can reduce the natural frequency of vibration. Truong [14] proposed a new type of hydraulic engine mount capable of inertial track control. Mathematical modeling was performed, and the simulation showed that the resonance peak and resonance frequency were changed through the optimal tuning of the inertia track region.
Kamada et al. [15] controlled active vibration by combining a column of a structure and a stack-type piezoelectric actuator. They showed that the vibration of the structure can be effectively isolated by using a piezoelectric actuator. Loukil et al. [16] studied the vibration reduction for an input applied using an accumulator and when converting the mechanical vibration energy into electrical energy and feeding it to the input of a piezoelectric actuator. They showed effective vibration reduction ability. Sui et al. [17] made an automotive engine mount actuator using PZT (Lead zirconate titanate) with low power and quick reaction characteristics, which demonstrated vibration reduction. Bartel et al. [18,19] proposed a new type of engine mount that isolates engine vibration. This was designed to withstand the dynamic forces of the structure, thus reducing power requirements by using small actuators. Liette et al. [20] conducted studies to quantify active path interactions owing to passive system dynamic characteristics in a source–path–receiver system and vibration reduction for structures with path-coupled piezoelectric stack actuators and rubber mounts. Based on the source–path–receiver structure, a model with six degrees of freedom (6-DOF) motion equations was simulated and experimented to demonstrate the reduction of vibration.
Three mounts support an actual vehicle, and the 6-DOF, introduced in [20], expresses all the three paths but is limited to a two-dimensional motion. However, as the engine vibrates, the movement occurs not only on the front but also on the sides. Hence, the model is expanded to three paths to model conditions similar to those encountered by an actual automotive engine.
This study discusses a case of typical single-frequency excitation in the mid-frequency bands occurring in engines of electric and hybrid vehicles. The overall structure is based on the source–path–receiver structure and is modeled on the basis of the schematic diagram shown in Figure 1. The source is the vehicle engine, the path is the active path, and the receiver is the sub-frame. It is assumed here that the active path consists of a combination of a piezo stack actuator and a rubber mount with no mass. This study discusses the following three points: (1) mathematical 3D modeling of structures with three active paths combining piezo stack actuators and rubber mounts; (2) numerical calculation of the amplitude and phase of each actuator input of the active path based on the mathematical modeling and the investigation of vibration reduction; (3) application of the NLMS algorithm to improve the inconvenience of the numerical method and discussion of the effect of vibration reduction on the source and receiver.

2. Modeling of Mounting System

2.1. 9-DOF Modeling

The overall schematic of the 9-DOF mounting system is based on Figure 1 and modeled as shown in Figure 2.
The 3D model was created on the basis of the lumped parameter model and expressed using orthographic projection. The floor plan shows the location of the active structural paths; m s , m r , and m s t i represent the source, receiver, and actuator mass, respectively, that were actually measured and used; I s , y and I r , y represent the inertia of the source and receiver, respectively, in the y-direction. I s f , x and I r f , x represent the inertia of the source and receiver, respectively, in the x-direction and were derived from calculations; l s i and l r i represent the distance between each path on the front of the source and receiver, respectively; l f i and l r f i indicate the distance between each path on the sides of the source and receiver, respectively; d and d 2 are the distances of the shaker relative to each direction; θ s , y , θ r , y , θ s f , x , and θ r f , x represent the rotational displacement of the source and receiver relative to the front and side, respectively; k m i and k b r i represent the complex stiffness for each path; ε s , g 1 , z and ε r , g 2 , z indicate the translational motion of the source and receiver, respectively; ε s t i , z and f s t i represent the motion and force, respectively, of the ith actuator; w z represents the force of the shaker.
In this study, the translational motion assumed only the z-direction. The rotational displacement in the y-direction was assumed to be on the front. Only the x-direction was assumed to be on the sides. Moreover, the rotational displacement was assumed to be very small. The values of each parameter are summarized in Table 1. Equations (1)–(9) express the equations of motion of the front and sides:
m s ε ¨ s , g 1 , z + ( k m 1 + k m 3 + k m 2 ) ε s , g 1 , z k m 1 ε s t 1 , z k m 2 ε s t 2 , z k m 3 ε s t 3 , z + ( k m 2 l s 2 k m 1 ε s 1 ) θ s , y + ( k m 1 l f 1 + k m 2 l f 1 k m 3 l f 2 ) θ s f , x = w z ( t )
m r ε ¨ r , g 1 , z + ( k m 4 + k m 5 + k m 6 + k b r 1 + k b r 2 + k b r 3 + k b r 4 ) ε r , g 1 , z k m 4 ε s t 1 , z k m 5 ε s t 2 , z k m 6 ε s t 3 , z + ( k m 5 l s 1 k m 4 l s 1 k b r 1 l r 1 + k b r 2 l r 2 k b r 3 l r 3 + k b r 4 l r 4 ) θ r , y + ( k m 4 l f 1 + k m 5 l f 1 k m 6 l f 2 k b r 1 l r f 1 k b r 2 l r f 2 + k b r 3 l r f 3 + k b r 4 l r f 4 ) θ r f , x = 0
m s t 1 ε ¨ s t 1 , g 1 , z + ( k m 1 + k m 4 ) ε s t 1 , z k m 1 ε s , g 1 , z k m 4 ε r , g 1 , z + k m 1 l s 1 θ s , y + k m 4 l s 1 θ r , y k m 1 l f 1 θ s f , x k m 4 l f 1 θ r f , x = f s t 1 , z ( t )
m s t 2 ε ¨ s t 2 , z + ( k m 2 + k m 5 ) ε s t 2 , z k m 2 ε s , g 1 , z k m 5 ε r , g 1 , z k m 2 l s 2 θ s , y k m 5 l s 2 θ r , y k m 2 l f 1 θ s f , x k m 5 l f 1 θ r f , x = f s t 2 , z ( t )
m s t 3 ε ¨ s t 3 , z + ( k m 3 + k m 6 ) ε s t 3 , z k m 3 ε s , g 1 , z k m 6 ε r , g 1 , z k m 3 l f 2 θ s f , x k m 6 l f 2 θ r f , x = f s t 3 , z ( t ) ,
I s θ ¨ s , y + ( k m 2 l s 2 k m 1 l s 1 ) ε s , g 1 , z + k m 1 l s 1 ε s t 1 , z k m 2 l s 2 ε s t 2 , z + ( k m 2 l s 2 2 + k m 1 l s 1 2 ) θ s , y = w z ( t ) d
I r θ ¨ r , y + ( k m 5 l s 2 k m 4 l s 2 + k b r 2 l r 2 k b r 1 l r 1 ) ε r , g 1 , z + k m 4 l s 1 ε s t 1 , z k m 5 l s 2 ε s t 2 , z + ( k m 4 l s 1 2 + k m 5 l s 2 2 + k b r 1 l r 1 2 + k b r 2 l r 2 2 + k b r 3 l r 3 2 + k b r 4 l r 4 2 ) θ r , y = 0
I s θ ¨ s f , x + ( k m 1 l f 1 + k m 2 l f 1 k m 3 l f 2 ) ε s , g 1 , z + k m 3 l f 2 ε s t 3 , z k m 1 l f 1 ε s t 1 , z k m 2 l f 1 ε s t 2 , z + ( k m 1 l f 1 2 + k m 2 l f 1 2 + k m 3 l f 2 2 ) θ s f , x = w z ( t ) d 2
I r θ ¨ r f , x + ( k m 4 l f 1 + k m 5 l f 1 k m 6 l f 2 k b r 1 l r f 1 k b r 2 l r f 3 + k b r 3 l r f 2 + k b r 4 l r f 4 ) ε r , g 1 , z + k m 6 l f 2 ε s t 3 , z k m 4 l f 1 ε s t 1 , z k m 5 l f 1 ε s t 2 , z + ( k m 4 l f 1 2 + k m 5 l f 1 2 + k m 6 l f 2 2 + k b r 1 l r f 1 2 + k b r 2 l r f 3 2 + k b r 3 l r f 2 2 + k b r 4 l r f 4 2 ) θ r f , x = 0
To express the equation of motion in a simple matrix form, the mass, stiffness, forces, and displacement vector are defined in Equations (10)–(14).
M = d i a g ( { m s m r m s t 1 m s t 2 m s t 3 I s , y I r , y I s , x I r , x } )
K = [ k m 2 + k m 3 + k m 4 0 k m 1 k m 2 k m 3 0 k m 4 + k m 5 + k m 6 + k b r 1 + k b r 2 k m 4 k m 5 k m 6 k m 1 k m 4 k m 1 + k m 4 0 0 k m 2 k m 5 0 k m 2 + k m 5 0 k m 3 k m 6 0 0 k m 3 + k m 6 k m 3 l s 2 k m 1 l s 1 0 k m 1 l s 1 k m 2 l s 2 0 0 k m 5 l s 2 k m 4 l s 1 k b r 1 l r 2 + k b r 2 l r 2 k m 4 l s 1 k m 5 l s 2 0 k m 1 l f 1 k m 2 l f 1 k m 3 l f 2 0 k m 1 l f 1 k m 2 l f 1 k m 3 l f 2 0 k m 4 l f 1 k m 5 l f 1 k m 6 l f 2 k m 4 l f 1 k m 5 l f 1 k m 6 l f 2 k m 3 l s 2 k m 1 l s 1 0 k m 1 l f 1 k m 2 l f 1 k m 3 l f 2 0 0 k m 5 l s 2 k m 4 l s 1 k b r 1 l r 2 + k b r 2 l r 2 0 k m 4 l f 1 k m 5 l f 1 k m 6 l f 2 k m 1 l s 1 k m 4 l s 1 k m 1 l f 1 k m 4 l f 1 k m 2 l s 2 k m 5 l s 2 k m 2 l f 1 k m 5 l f 1 0 0 k m 3 l f 2 k m 6 l f 2 k m 4 l s 1 2 + k m 2 l s 2 2 0 0 0 0 k b r 1 l r 1 2 + k b r 2 l r 2 2 + k m 4 l s 1 2 + k m 5 l s 2 2 0 0 0 0 k m 1 l f 1 2 + k m 2 l f 1 2 + k m 3 l f 2 2 0 0 0 0 k m 4 l f 1 2 + k m 5 l f 1 2 + k m 6 l f 2 2 ]
q ( t ) = { ε s , g 1 , z ( t ) ε r , g 1 , z ( t ) ε s t 1 , z ( t ) ε s t 2 , z ( t ) ε s t 3 , z ( t ) θ s , y ( t ) θ r , y ( t ) θ s f , x ( t ) θ r f , x ( t ) }
w ( t ) = { w z ( t ) 0 0 0 0 w z ( t ) d 0 w z ( t ) d 2 0 }
F ( t ) = { 0 0 f s t 1 , z ( t ) f s t 2 , z ( t ) f s t 3 , z ( t ) 0 0 0 0 }
The focus of this study was on the vibration reduction effect by the active path. To investigate the vibration reduction performance of the part where each active path is located, a coordinate transformation was performed using the transfer matrix defined in Equation (15).
= [ l s 1 l s 1 + l s 2 l s 2 l s 1 + l s 2 l f 2 l f 1 + l f 2 0 0 0 0 0 0 0 0 0 0 0 0 l s 1 l s 1 + l s 2 l s 2 l s 1 + l s 2 l f 2 l f 1 + l f 2 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 1 l s 1 + l s 2 1 l s 1 + l s 2 0 0 0 0 0 0 0 0 0 0 0 0 0 1 l s 1 + l s 2 1 l s 1 + l s 2 0 1 l s 1 + l s 2 1 l s 1 + l s 2 1 l f 1 + l f 2 0 0 0 0 0 0 0 0 0 0 0 0 1 l s 1 + l s 2 1 l s 1 + l s 2 1 l f 1 + l f 2 ]
The transformed displacement is defined as q = q , and the displacement matrix q for each active path is redefined as shown in Equation (16). If the equation of motion is rewritten, it is expressed as Equation (17). The schematic representation of the transformed coordinates is shown in Figure 3.
q = { ξ s t 1 , g 1 , z ( t ) ξ s t 2 g 1 , z ( t ) ξ s t 3 , g 1 , z ( t ) ε s t 1 , z ( t ) ε s t 2 , z ( t ) ε s t 3 , z ( t ) ξ s t 1 , g 2 , z ( t ) ξ s t 2 , g 2 , z ( t ) ξ s t 3 , g 2 , z ( t ) }
M q ¨ + K q = F ( x ) + W ( t )

2.2. Calculation of Actuator Amplitude and Phase

To define the relationship between harmonic excitation and actuator response, each input is defined by Equations (18) and (19):
w z ( t ) = W z e i ω t
f s t i , z ( t ) = F s t i , z e i ( ω t + ϕ s t i , z )
where W Z is the amplitude of the shaker, and F s t i and ϕ s t i stand for the amplitude and phase of the actuator in the ith path, respectively. The displacement of the system according to the excitation force is defined for each active path and expressed by Equation (20):
ξ i ( t ) = ξ s t i , g 1 ( t )
ξ s t i , g 1 ( t ) is the displacement of the ith actuator with inertia g 1 and is affected by the excitation force of the shaker and each actuator. It can be defined by the Equations (21)–(23) as follows:
ξ 1 ( t ) = ( Ξ s 1 , 1 + Ξ s 1 , s t 1 e i ϕ s t 1 + Ξ s 1 , s t 2 e i ϕ s t 2 + Ξ s 1 , s t 3 e i ϕ s t 3 ) e i ω t
ξ 2 ( t ) = ( Ξ s 2 , 1 + Ξ s 2 , s t 1 e i ϕ s t 1 + Ξ s 2 , s t 2 e i ϕ s t 2 + Ξ s 2 , s t 3 e i ϕ s t 3 ) e i ω t
ξ 3 ( t ) = ( Ξ s 3 , 1 + Ξ s 3 , s t 1 e i ϕ s t 1 + Ξ s 3 , s t 2 e i ϕ s t 2 + Ξ s 3 , s t 3 e i ϕ s t 3 ) e i ω t
Ξ s i , 1 is the amplitude of the shaker, and Ξ s i , s t 1 , Ξ s i , s t 2 , and Ξ s i , s t 3 represent the amplitudes of the actuator for the ith active path. To calculate the amplitude and phase, a compliance matrix, represented by Equation (24), is used to define the amount of deformation caused by a unit load.
H = [ H 11 H 12 H 13 H 14 H 15 H 16 H 17 H 18 H 19 H 21 H 22 H 23 H 24 H 25 H 26 H 27 H 28 H 29 H 31 H 32 H 33 H 34 H 35 H 36 H 37 H 38 H 39 H 41 H 42 H 43 H 44 H 45 H 46 H 47 H 48 H 49 H 51 H 52 H 53 H 54 H 55 H 56 H 57 H 58 H 59 H 61 H 62 H 63 H 64 H 65 H 66 H 67 H 68 H 69 H 71 H 72 H 73 H 74 H 75 H 76 H 77 H 78 H 79 H 81 H 82 H 83 H 84 H 85 H 86 H 87 H 88 H 89 H 91 H 92 H 93 H 94 H 95 H 96 H 97 H 98 H 99 ]
The amplitude displacement vector for each active path is defined in Equation (25):
Q = { Ξ s t 1 , g 1 Ξ s t 2 , g 1 Ξ s t 3 , g 1 s t 1 s t 2 s t 3 Ξ s t 1 , g 2 Ξ s t 2 , g 2 Ξ s t 3 , g 2 }
The amplitude is calculated by Equation (26):
Q e i ω t = H { F + W } e i ω t
The amplitudes generated from the active paths are defined by Equations (27)–(29).
Ξ s 1 , 1 = ( H 11 + H 16 d + H 18 d ) W z Ξ s 1 , s t 1 = H 13 F s t 1 , Ξ s 1 , s t 2 = H 14 F s t 2 , Ξ s 1 , s t 3 = H 15 F s t 3
Ξ s 2 , 1 = ( H 21 + H 26 d + H 28 d ) W z Ξ s 2 , s t 1 = H 23 F s t 1 , Ξ s 2 , s t 2 = H 24 F s t 2 , Ξ s 2 , s t 3 = H 25 F s t 3
Ξ s 3 , 1 = ( H 31 + H 36 d + H 38 d ) W z Ξ s 3 , s t 1 = H 33 F s t 1 , Ξ s 3 , s t 2 = H 34 F s t 2 , Ξ s 3 , s t 3 = H 35 F s t 3
The phases, defined by Equations (30)–(32), are defined on the basis of amplitude:
β s 1 , 1 = ( H 11 + H 16 d + H 18 d ) β s 1 , s t 1 = H 13 , β s 1 , s t 2 = H 14 , β s 1 , s t 3 = H 15
β s 2 , 1 = ( H 21 + H 26 d + H 28 d ) β s 2 , s t 1 = H 23 , β s 2 , s t 2 = H 24 , β s 2 , s t 3 = H 25
β s 3 , 1 = ( H 31 + H 36 d + H 38 d ) β s 3 , s t 1 = H 33 , β s 3 , s t 2 = H 34 , β s 3 , s t 3 = H 35
where stands for the phase operator, and β s n , 1 represents the phase of the amplitude generated by the shaker in the nth active path; β s n , s t m represents the phase of the amplitude generated by the mth actuator in the nth active path. Using Equations (27)–(32), the magnitude and phase terms are defined by Equations (33)–(35):
ξ 1 ( t ) = ( | Ξ s 1 , 1 | e i β s 1 , 1 + | Ξ s 1 , s t 1 | e i ( β s 1 , s t 1 + ϕ s t 1 ) + | Ξ s 1 , s t 2 | e i ( β s 1 , s t 2 + ϕ s t 2 ) + | Ξ s 1 , s t 3 | e i ( β s 1 , s t 3 + ϕ s t 3 ) ) e i ω t
ξ 2 ( t ) = ( | Ξ s 2 , 1 | e i β s 2 , 1 + | Ξ s 2 , s t 1 | e i ( β s 2 , s t 1 + ϕ s t 1 ) + | Ξ s 2 , s t 2 | e i ( β s 2 , s t 2 + ϕ s t 2 ) + | Ξ s 2 , s t 3 | e i ( β s 2 , s t 3 + ϕ s t 3 ) ) e i ω t
ξ 3 ( t ) = ( | Ξ s 3 , 1 | e i β s 3 , 1 + | Ξ s 3 , s t 1 | e i ( β s 3 , s t 1 + ϕ s t 1 ) + | Ξ s 3 , s t 2 | e i ( β s 3 , s t 2 + ϕ s t 2 ) + | Ξ s 3 , s t 3 | e i ( β s 3 , s t 3 + ϕ s t 3 ) ) e i ω t
These equations have seven-phase terms in one equation. Each phase has four phases that are generated by the shaker and actuator and three phases that occur only in the actuator. Because of the presence of various phases, it is difficult to control the motion. Therefore, we can assume that the phases generated by the shaker match each other. By this definition, it is possible to create a phase difference motion in each active path. The expression of phase for active path 1 can be written as β s 1 , 1 = β s 1 , s t 1 + ϕ s t 1 . Similarly, the phases for active paths 2 and 3 are β s 1 , 1 = β s 1 , s t 2 + ϕ s t 2 and β s 1 , 1 = β s 1 , s t 3 + ϕ s t 3 , respectively. The phases generated by the shaker are represented by Equations (36)–(38):
ϕ s t 1 = β s 1 , 1 β s 1 , s t 1
ϕ s t 2 = β s 1 , 1 β s 1 , s t 2
ϕ s t 2 = β s 1 , 1 β s 1 , s t 2
From Equations (36)–(38), the phase of the actuator in each active path caused an out-of-phase motion owing to the difference between the phase generated by the shaker and that of each path generated by the actuator. In the case of active paths 2 and 3, the phase of each actuator is defined as ϕ s t 1 = β s 2 , 1 β s 2 , s t 1 , ϕ s t 2 = β s 2 , 1 β s 2 , s t 2 , and ϕ s t 3 = β s 2 , 1 β s 2 , s t 3 ; and ϕ s t 1 = β s 3 , 1 β s 3 , s t 1 , ϕ s t 2 = β s 3 , 1 β s 3 , s t 2 , and ϕ s t 3 = β s 3 , 1 β s 3 , s t 3 , respectively. Equations (39)–(41) are obtained by substituting the equation of the actuator phase into Equations (33)–(35) and summarizing them to determine the phase of the shaker:
ξ 1 ( t ) = ( | Ξ s 1 , 1 | + | Ξ s 1 , s t 1 | + | Ξ s 1 , s t 2 | + | Ξ s 1 , s t 3 | ) e i ( ω t + β s 1 , 1 )
ξ 2 ( t ) = ( | Ξ s 2 , 1 | + | Ξ s 2 , s t 1 | + | Ξ s 2 , s t 2 | + | Ξ s 2 , s t 3 | ) e i ( ω t + β s 2 , 1 )
ξ 3 ( t ) = ( | Ξ s 3 , 1 | + | Ξ s 3 , s t 1 | + | Ξ s 3 , s t 2 | + | Ξ s 3 , s t 3 | ) e i ( ω t + β s 3 , 1 )
To reduce the motion to zero, the amplitudes F s t 1 , z , F s t 2 , z , and F s t 3 , z of each actuator are calculated assuming that their magnitude is zero. Hence, we define Equations (42)–(44):
| Ξ s 1 , 1 | + | Ξ s 1 , s t 1 | + | Ξ s 1 , s t 2 | + | Ξ s 1 , s t 3 | = 0
| Ξ s 2 , 1 | + | Ξ s 2 , s t 1 | + | Ξ s 2 , s t 2 | + | Ξ s 2 , s t 3 | = 0
| Ξ s 3 , 1 | + | Ξ s 3 , s t 1 | + | Ξ s 3 , s t 2 | + | Ξ s 3 , s t 3 | = 0
To determine F s t 1 , F s t 2 , and F s t 3 , Equation (45) is given:
[ | H 13 | | H 14 | | H 15 | | H 23 | | H 24 | | H 25 | | H 33 | | H 34 | | H 35 | ] { F s t 1 F s t 2 F s t 3 } = W z [ | H 11 + d H 16 + d H 18 | | H 21 + d H 26 + d H 28 | | H 31 + d H 36 + d H 38 | ]
Equation (46) is used to calculate the amplitude of each actuator:
{ F s t 1 F s t 2 F s t 3 } = W z | H 13 | | H 24 | | H 35 | | H 13 | | H 25 | | H 34 | | H 14 | | H 23 | | H 35 | + | H 14 | | H 25 | | H 33 | + | H 15 | | H 23 | | H 34 | | H 15 | | H 24 | | H 33 | × [ ( | H 14 | | H 34 | | H 15 | | H 34 | ) ( | H 21 + d H 26 + d H 28 | ) ( | H 24 | | H 35 | | H 25 | | H 34 | ) ( | H 11 + d H 16 + d H 18 | ) ( | H 14 | | H 25 | | H 15 | | H 24 | ) ( | H 31 + d H 36 + d H 38 | ) ( | H 23 | | H 35 | | H 25 | | H 33 | ) ( | H 11 + d H 16 + d H 18 | ) ( | H 14 | | H 25 | | H 15 | | H 23 | ) ( | H 31 + d H 36 + d H 38 | ) ( | H 13 | | H 35 | | H 15 | | H 33 | ) ( | H 21 + d H 26 + d H 28 | ) ( | H 13 | | H 34 | | H 14 | | H 33 | ) ( | H 21 + d H 26 + d H 28 | ) ( | H 23 | | H 34 | | H 24 | | H 33 | ) ( | H 11 + d H 16 + d H 18 | ) ( | H 13 | | H 24 | | H 14 | | H 23 | ) ( | H 31 + d H 36 + d H 38 | ) ]

3. Validation Using the Numerical Simulation

3.1. Simulation Overview

Simulations were conducted to determine the validity of the mathematical model. The input and sampling frequencies were set to u ( t ) = 10 sin ( 800 π t ) and 15 kHz, respectively. The equations of motion were expressed using the state space equations and are represented by Equations (47) and (48):
x ˙ ( t ) = A x ( t ) + B u ( t )
y ( t ) = C x ( t ) + D u ( t )
where A is a matrix representing the state of the system, B is an input matrix, and C is an output matrix. Each matrix is expressed by Equations (49)–(52). In addition, the output y ( t ) for the system has six outputs owing to three source and three receiver parts for the active path and is expressed by Equation (53):
A = [ O n × n I n × n K M C M ]
B = [ O n × n 1 M ]
C = [ I 3 × 3 O 3 × 3 O 3 × 3 O 3 × 3 O 3 × 3 I 3 × 3 O 9 × 9 ]
D = [ O 6 × 9 ]
y ( t ) = { a p 1 , g 1 a p 2 , g 1 a p 3 , g 1 a p 1 , g 2 a p 2 , g 2 a p 3 , g 2 } T
where a p n , g m represents the acceleration in the nth path in mass with the inertia of m.
Although the state of the system that can be measured is limited, a total state value for the system is required to construct the feedback control system. Hence, a state observer wa designed as shown in Figure 4.
Using the real system matrices, A, B, and C, the G, H, and C matrices wre constructed for the estimates. Because the state space equations are defined in continuous time, they were transformed into discrete time matrices through bilinear transformation. The observer gain L was selected appropriately, and the agreement between the estimated and actual output was confirmed, as shown in Figure 5. Therefore, the validity of the model was verified, and it was applied to the feedback control system.

3.2. Control Using the Calculated Amplitude and Phase

A simulation was performed to determine whether the vibration transmitted to the source and receiver was reduced using the amplitude and phase values of the calculated actuator. The equations for the actuator input are defined by Equations (54)–(56):
f s t 1 ( t ) = F s t 1 sin ( ω t + ϕ s t 1 )
f s t 2 ( t ) = F s t 2 sin ( ω t + ϕ s t 2 )
f s t 3 ( t ) = F s t 3 sin ( ω t + ϕ s t 3 )
The simulation was performed with active path 3 as a target, and it was confirmed whether the vibration transmitted to the source and receiver were reduced when the actuators in all active paths were turned on. Figure 6 and Figure 7 show the results of the vibration of the source and receiver, respectively. Table 2 shows the values of the source and receiver outputs.
The results in Table 2 show that the vibration of a s t 3 , g 1 was sharply reduced by about 98.19% because the active path 3 was used to calculate the actuator value for vibration isolation purposes. For a s t 1 , g 1 and a s t 2 , g 1 , which were untargeted, the source vibration values were greatly increased by about 450.47% and 91.28%, respectively. However, the receiver results show that the vibration of a s t 1 , g 2 was reduced by about 57% when compared to the uncontrolled results. The vibrations of a s t 2 , g 2 and a s t 3 , g 2 were decreased by about 55 and 37%, respectively.

3.3. Control Using the NLMS Algorithm

To isolate the source vibration of the target path in the previous section, the amplitude and phase of the actuator were calculated. However, it is evident that the vibration was greatly increased for the untargeted path. In addition, the calculation method cannot be applied to a signal with a real multi-frequency. Therefore, both the source and the receiver cannot control vibrations satisfactorily. Therefore, to compare with the computational isolation method, the simulation was carried out for the source and receiver vibration isolation using the NLMS algorithm. The simulation was based on the schematic representation shown in Figure 8 and was used as an input to the actuator using the signal from the source.
The inputs of the simulations were provided under the same conditions for comparison, and the signals entering the actuator were applied by appropriate adjustments to the signals from the NLMS algorithm. The simulation results of the source and receiver are shown in Figure 9 and Figure 10, respectively. Table 3 presents the values for each output.
The result of the source part of Table 3 shows that the vibration of a s t 1 , g 1 increased by about 105.62% compared to the previous value, but those of a s t 2 , g 1 and a s t 3 , g 1 decreased by 22.66 and 65.84%, respectively. In addition, the results for the receiver show that the vibrations of a s t 1 , g 2 , a s t 2 , g 2 , and a s t 3 , g 2 decreased by 10.08, 4.86, and 90.29%, respectively. The reduction effect on the receiver was generally good for the calculation. However, when the results of the calculation in the source part wee compared with the those of the control using NLMS, the vibrations of a s t 1 , g 1 and a s t 2 , g 1 tended to greatly increase in the calculation method, but the result of using the NLMS algorithm tended to decrease the vibration in a s t 2 , g 1 . Moreover, the result of the calculation of vibration for a s t 1 , g 1 showed a tendency to increase by 450.47%, which is about 5.5 times, whereas for the result using the NLMS algorithm, the vibration tended to increase by 105.62%, i.e., it doubled. Therefore, it can be seen that the use of the NLMS algorithm in the calculation method is more efficient for the whole system than the use of the actuator input.

4. Conclusions

This study was conducted to reduce vibration by coupling a piezo actuator and rubber mount between the plates of vehicle engines. The main contributions of this study are as follows: (1) Based on the source–path–receiver structure, three-dimensional motion equations for three active-path structures were established. Considering the effects on lateral rotational displacement, this enables a more realistic simulation compared to the modeling established in previous studies. (2) The NLMS algorithm was applied to the active control and compared to the calculation method of the actuator input. Through this result, it was shown that applying the NLMS algorithm is effective for overall vibration reduction.
The plate structure with active paths was modeled on the basis of the lumped parameter model and expressed using the equation of motion. When excited by a single frequency, the vibration-reduced effects were evaluated by comparing the results obtained by using the actuator input with those obtained by applying the NLMS algorithm. The computing method using the actuator input showed that the vibration of the receiver was largely reduced. However, the vibration of the non-targeted source was greatly increased. Moreover, this method is troublesome because it requires to newly define the modeling according to the structure change and cannot be applied when a multi-frequency input is applied. To solve this inconvenience, the NLMS algorithm was applied. When the NLMS algorithm was used, the vibration for the active path 1 in the source tended to increase, but those of the other active paths tended to reduce overall in the source and receiver. In particular, it was seen that the vibration of the active path 3 was largely isolated. Hence, future works should be conducted to study the relationship between each active path to control vibration more efficiently. In addition, an experimental verification of the results of this study should be carried out. Finally, our results show that vibrations are reduced for single-frequency inputs; hence, simulation and experiments for multi-frequency signals should be performed in the future.

Author Contributions

D.H. and B.K. initiated and developed the ideas related to this research work. D.H. and B.K. developed novel methods, derived relevant formulations, and carried out performance analyses and numerical analyses. D.H. wrote the paper draft under B.K.’s guidance, and D.H. finalized the paper.

Funding

This paper was supported by the National Research Foundation of Korea (NRF) grant funded by the government of Korea (Ministry of Science, ICT & Future Planning) (No. 2016R1C1B1010888).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Minak, G.; Fragassa, C.; De Camargo, F.V. A Brief Review on Determinant Aspects in Energy Efficient Solar Car Design and Manufacturing. In Proceedings of the International Conference on Sustainable Design and Manufacturing, Bologna, Italy, 26–28 April 2017; Volume 68, pp. 847–856. [Google Scholar]
  2. Hosseini, A.M.; Arzanpour, S.; Golnaraghi, F.; Parameswaran, A.M. Solenoid actuator design and modeling with application in engine vibration isolators. J. Vib. Control 2012, 19, 1015–1023. [Google Scholar] [CrossRef]
  3. Tang, X.; Yang, W.; Hu, X.; Zhang, D. A novel simplified model for torsional vibration analysis of a series-parallel hybrid electric vehicle. Mech. Syst. Signal Process. 2017, 85, 329–338. [Google Scholar] [CrossRef]
  4. Kraus, R.; Herold, S.; Millitzer, J.; Jungblut, T. Development of Active Engine Mounts Based on Piezo Actuators. ATZ Worldw. 2014, 116, 46–51. [Google Scholar] [CrossRef]
  5. Harun, M.H.; Azhari, M.A.; Yunos, M.R.M.; Yamin, A.K.M. Characterization of a Magnetorheological Fluid Damper Applied to Semi-Active Engine Mounting System. In Proceedings of the 6th ECCOMAS Conference on Smart Structures and Materials, Glasgow, UK, 11–15 June 2018; Volume 5, pp. 2448–2459. [Google Scholar]
  6. Hausberg, F.; Scheiblegger, C.; Pfeffer, P.; Plöchl, M.; Hecker, S.; Rupp, M. Experimental and analytical study of secondart path variations in active engine mounts. J. Sound Vib. 2015, 340, 22–38. [Google Scholar] [CrossRef]
  7. Yang, T.; Suai, Z.; Sun, Y.; Zhu, M.; Xiao, Y.; Liu, X.; Du, J.; Jin, G.; Liu, Z. Active vibration isolation system for a diesel engine. Noise Control Eng. J. 2012, 60, 267–282. [Google Scholar] [CrossRef]
  8. Sun, W.; Li, Y.; Huang, J.; Zhang, N. Vibration effect and control of In-Wheel Switched Reluctance Motor for electric vehicle. J. Sound Vib. 2015, 338, 105–120. [Google Scholar] [CrossRef]
  9. Chae, H.D.; Choi, S.-B. A new vibration isolation bed stage with magnetorheological dampers for ambulance vehicles. Smart Mater. Struct. 2014, 24, 17001. [Google Scholar] [CrossRef]
  10. Jeon, J.; Han, Y.-M.; Lee, D.-Y.; Choi, S.-B. Vibration control of the engine body of a vehicle utilizing the magnetorheological roll mount and the piezostack right-hand mount. Proc. Inst. Mech. Eng. D J. Autom. Eng. 2013, 227, 1562–1577. [Google Scholar] [CrossRef]
  11. Fakhari, V.; Choi, S.-B.; Cho, C.-H. A new robust adaptive controller for vibration control of active engine mount subjected to large uncertainties. Smart Mater. Struct. 2015, 24, 45044. [Google Scholar] [CrossRef]
  12. Elahinia, M.; Ciocanel, C.; Nguyen, T.M.; Wang, S. MR- and ER-Based Semiactive Engine Mounts: A Review. Smart Mater. Res. 2013, 2013, 1–21. [Google Scholar] [CrossRef]
  13. Wu, W.; Chen, X.; Shan, Y. Analysis and experiment of a vibration isolator using a novel magnetic spring with negative stiffness. J. Sound Vib. 2014, 333, 2958–2970. [Google Scholar] [CrossRef]
  14. Truong, T.Q.; Ahn, K.K. A new type of semi-active hydraulic engine mount using controllable area of inertia track. J. Sound Vib. 2010, 329, 247–260. [Google Scholar] [CrossRef]
  15. Kamada, T.; Fujita, T.; Hatayama, T.; Arikabe, T.; Murai, N.; Aizawa, S.; Tohyama, K. Active vibration control of frame structures with smart structures using piezoelectric actuators (Vibration control by control of bending moments of columns). Smart Mater. Struct. 1997, 6, 448–456. [Google Scholar] [CrossRef]
  16. Loukil, T.; Bareille, O.; Ichchou, M.; Haddar, M. A low power consumption control scheme: Application to a piezostack-based active mount. Front. Mech. Eng. 2013, 8, 383–389. [Google Scholar] [CrossRef]
  17. Sui, L.; Xiong, X.; Shi, G. Piezoelectric Actuator Design and Application on Active Vibration Control. Phys. Procedia 2012, 25, 1388–1396. [Google Scholar] [CrossRef]
  18. Herold, S.; Kraus, R.; Millitzer, J.; de Rue, G. Vibration control of a medium-sized vehicle by a novel active engine mount. In Proceedings of the 4th PT PIESA Symposium, Nuremberg, Germany, 26 April 2013. [Google Scholar]
  19. Bartel, T.; Herold, S.; Mayer, D.; Melz, T. Development and testing of active vibration control systems with piezoelectric actuators. In Proceedings of the 6th Eccomas Conference on Smart Structures and Materials, Torino, Italy, 24–26 June 2013. [Google Scholar]
  20. Liette, J.; Dreyer, J.T.; Singh, R. Interaction between two active structural paths for source mass motion control over mid-frequency range. J. Sound Vib. 2014, 333, 2369–2385. [Google Scholar] [CrossRef]
Figure 1. Powertrain with three active structural paths.
Figure 1. Powertrain with three active structural paths.
Applsci 09 00711 g001
Figure 2. Modeling of source-active structural path–receiver.
Figure 2. Modeling of source-active structural path–receiver.
Applsci 09 00711 g002
Figure 3. Schematic representation (a) before coordinate change; (b) after coordinate change.
Figure 3. Schematic representation (a) before coordinate change; (b) after coordinate change.
Applsci 09 00711 g003
Figure 4. State observer design.
Figure 4. State observer design.
Applsci 09 00711 g004
Figure 5. Result of the state observer.
Figure 5. Result of the state observer.
Applsci 09 00711 g005
Figure 6. Result of the vibration of the source related to each active path; the blue and red lines represent the vibration of the source before and after applying the control, respectively.
Figure 6. Result of the vibration of the source related to each active path; the blue and red lines represent the vibration of the source before and after applying the control, respectively.
Applsci 09 00711 g006
Figure 7. Result of the vibration of the receiver related to each active path; the blue and red lines represent the vibration of the receiver before and after applying the control, respectively.
Figure 7. Result of the vibration of the receiver related to each active path; the blue and red lines represent the vibration of the receiver before and after applying the control, respectively.
Applsci 09 00711 g007
Figure 8. Schematic representation of the simulation using the NLMS (Normalized Least Mean Squares) algorithm.
Figure 8. Schematic representation of the simulation using the NLMS (Normalized Least Mean Squares) algorithm.
Applsci 09 00711 g008
Figure 9. Result of the vibration of the source related to each active path using the NLMS algorithm; the blue and red lines represent the vibration of the source before and after applying the control, respectively.
Figure 9. Result of the vibration of the source related to each active path using the NLMS algorithm; the blue and red lines represent the vibration of the source before and after applying the control, respectively.
Applsci 09 00711 g009
Figure 10. Result of the vibration of the receiver related to each active path using the NLMS algorithm; the blue and red lines represent the vibration of the receiver before and after applying the control, respectively.
Figure 10. Result of the vibration of the receiver related to each active path using the NLMS algorithm; the blue and red lines represent the vibration of the receiver before and after applying the control, respectively.
Applsci 09 00711 g010
Table 1. Variable values and units for 9 degrees of freedom (DOF).
Table 1. Variable values and units for 9 degrees of freedom (DOF).
VariableValuesUnitsVariableValuesUnits
m s = m r 1.081 kg k m 5 0.53 ( 1 + i 0.256 ) kN mm 1
m s t 1 = m s t 3 0.067 kg k m 6 0.61 ( 1 + i 0.300 ) kN mm 1
m s t 2 0.075 kg k b r n ( n = 1 , 2 ) 0.42 ( 1 + i 0.300 ) kN mm 1
I s , y = I r , y 16.172 g m 2 l s n ( n = 1 , 2 ) 100 mm
I s f , x = I r f , x 5.940 g m 2 l r n ( n = 1 , 2 ) 140 mm
k m 1 5.46 ( 1 + i 0.034 ) kN mm 1 l f n ( n = 1 , 2 ) 80 mm
k m 2 5.46 ( 1 + i 0.034 ) kN mm 1 l r f n ( n = 1 ~ 4 ) 100 mm
k m 3 2.48 ( 1 + i 0.036 ) kN mm 1 d 50 mm
k m 4 0.61 ( 1 + i 0.300 ) kN mm 1 d 2 10 mm
Table 2. Vibration reduction when the actuator was turned on.
Table 2. Vibration reduction when the actuator was turned on.
Source a s t 1 , g 1 a s t 2 , g 1 a s t 3 , g 1
Original 0.533 1.514 1.095
All actuators turned on 2.934 ( 450.47 % ) 2.896 ( 91.28 % ) 0.019 ( 98.19 % )
Receiver a s t 1 , g 2 a s t 2 , g 2 a s t 3 , g 2
Original 0.2382 0.37 0.6122
All actuators turned on 0.1021 ( 57.13 % ) 0.165 ( 55.18 % ) 0.384 ( 37.21 % )
Table 3. Result of vibration reduction when the actuator was turned on using the NLMS signal.
Table 3. Result of vibration reduction when the actuator was turned on using the NLMS signal.
Source a s t 1 , g 1 a s t 2 , g 1 a s t 3 , g 1
Original 0.533 1.514 1.095
Control using the NLMS 1.096 ( 105.62 % ) 1.171 ( 22.66 % ) 0.374 ( 65.84 % )
Receiver a s t 1 , g 2 a s t 2 , g 2 a s t 3 , g 2
Original 0.2382 0.37 0.6122
Control using the NLMS 0.214 ( 10.08 % ) 0.352 ( 4.86 % ) 0.059 ( 90.29 % )

Share and Cite

MDPI and ACS Style

Hong, D.; Kim, B. Quantification of Active Structural Path for Vibration Reduction Control of Plate Structure under Sinusoidal Excitation. Appl. Sci. 2019, 9, 711. https://doi.org/10.3390/app9040711

AMA Style

Hong D, Kim B. Quantification of Active Structural Path for Vibration Reduction Control of Plate Structure under Sinusoidal Excitation. Applied Sciences. 2019; 9(4):711. https://doi.org/10.3390/app9040711

Chicago/Turabian Style

Hong, Dongwoo, and Byeongil Kim. 2019. "Quantification of Active Structural Path for Vibration Reduction Control of Plate Structure under Sinusoidal Excitation" Applied Sciences 9, no. 4: 711. https://doi.org/10.3390/app9040711

APA Style

Hong, D., & Kim, B. (2019). Quantification of Active Structural Path for Vibration Reduction Control of Plate Structure under Sinusoidal Excitation. Applied Sciences, 9(4), 711. https://doi.org/10.3390/app9040711

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