Active vibration control analysis in smart composite structures using ANSYS

Latest revision as of 14:18, 24 February 2021

Abstract

Piezoelectric Macro-Fiber Composite (MFC) utilization is increasing in engineering fields due to its strong actuation forces and high flexibility. In this paper, piezoelectric ( ${\textstyle {d}_{33}MFC}$ ) type (M8528-P1) patches are applied for active control of a cantilever composite beam. A linear coupled ﬁnite element model for piezoelectric MFC actuation of the composite beam was developed based on the APDL-ANSYS codes by using linear piezoelectric constitutive equations to study smart composite beam behavior in open and closed-loop cases. Active vibration control response of laminated composite beams with various stacking sequence configurations was examined, and the results were compared. In this, the proportional type ( $kp$ ) of control algorithms is utilized, and in the summary of work refined two sets of finite element models. The first set their laminates to have an orientation (0 with 90), and others have an orientation (0 with 45). When the control signal applied with the gain ( $kp$ ) to a system increases, the rise time generally decreases, the actuator voltages for different control gain for all cases are observed with the increase in the proportional constant. The findings of this study indicate that the composite beams composed from (0/90)stacking sequence sets had reached stability state faster than the composite beams composed from (0/45) stacking sequence sets.

Keywords: MFC piezoelectric, smart composite beam, FEM (APDL Codes), active vibration

1. Introduction

Vibration is present in countless daily life applications, and more often, vibration is highly undesirable. Undesirable vibration may decrease the performance of the product and cause safety or economic issues [1]. Vibration damping includes active, passive, and hybrid damping techniques. It is important to understand that passive damping can be very efficient in damping out high-frequency excitation, while active damping can be used to control low-frequency vibrations [2]. Active vibration-reducing techniques turn out to be a reality for a wide variety of applications [3]. An active method integrates sensors and actuators with a flexible structure, operated by a control scheme. An essential issue in active control systems is sensor and actuator selection [4]. The ever-increasing use of smart composite materials in advanced areas such as aerospace, automotive, and sports equipment has motivated researchers to explore behavior and performance characteristics [5].

There are some previous studies reported in this work. The simple cantilever beams active vibration control is investigated in [6-10]. Piezoelectric patches are mounted on the beams as actuators. With the purpose of sensing the vibration level, a strain-gage or another piezoelectric patch can be used. The identification system and placement control technique was utilized in [6]. The piezo patches and beam structure’s finite element model is constructed, and the closed-loop control is applied in [7,8]. Also, utilized the model of the beam and piezo-patches finite element; however, modal control techniques were employed Singh et al. [9]. Xu and Koko used the commercial finite element package ANSYS to state the findings [10]. The control design was conducted in the space state that was set in the analysis of the finite element model. For the control design, MATLAB as a system toolbox was utilized in their experiment. The sensors/actuators influence site was investigated. It was noticed that the site near the fixed end for the vibration control was better. Lim investigated several modes for vibration control of a fixed square plate by using devices of location discrete sensor/actuator at maximum strain points [11]. Quek et al. introduced a typical placement method of pairs of piezoelectric sensor /actuator for the vibration control sealed compound plates [12]. Xiamin et al. used four-bar linkage to investigate active vibration control [13]. The model of the finite element is applied and the reduced mode, robust ${\textstyle {H}_{\infty }}$ , and standard H∞ control methods are analyzed. Halim and Moheimani investigated a set of proposals to a criterion for ideal placement of a collected piezoelectric sensor-actuator pair on their supported plate by using modal and locative controllability [14]. Bendine et al. studied the modeling of the active vibration control of a smart functionally graded (FGM) beam based on a higher-order shear deformation theory [15]. The locations of piezoelectric actuators and sensors have a great influence on the reliability and control system efficiency. Malgaca studied optimal parameters with respect to positive position feedback used in the real-time control algorithm determined with the objective functions [16]. Flexible single-link composite manipulator with [45/-45] and [0/90] lay-ups are performed to reduce endpoint vibrations. Takacs and Rohal-Ilkiv [17] proposed a direct closed-loop prototype method in the active vibration control system modeled in the finite element method using APDL language for the clamped-free aluminum beam piezoelectric actuation using a linear quadratic controller (LQ) with a Kalman-ﬁlter. Vashist and Chhabra [18] found the optimal piezoelectric actuator’s locations by using an integer coded genetic algorithm GA on their plate where the singular value decomposition’s control matrix was used as a fitness function. Similar results of the frequency responses between experimental and Finite Element Method were found. Yavuz et al. analyzed a system of active vibration control’s closed-loop with the purpose of suppressing the vibrations of the end-point [19]. They established the system’s mathematical model by using the Lagrange equations. The digital solution is achieved through mixing the Newmark technique with the control action. They applied PID control action with the purpose of finding the value of the actuator signal at the appropriate time. Liu X. et al. studied the piezoelectric actuators’ placement optimization and the membrane structure’s active vibration control [20]. They designed the classical linear quadratic regulator controllers with the aim of suppressing the undesired vibration. Simulation findings indicated that the piezoelectric actuators’ optimal locations are influenced by the actuator’s additional mass and stiffness. Van et al. developed a finite element model according to the theory of first-order shear deformation for active vibration control and optimal placement of laminated composite plates with linked pairs of distributed piezoelectric sensor/actuator [21]. Baghaee et al. developed the aero elastic panel flutter of the laminate composite plate to include two piezoelectric (MFC-type) [22]. They determined MFC actuator’s input voltage through the velocity or proportional feedback control algorithms according to the sensor output. Usman et al. used piezoelectric MFC energy harvesting on the cantilever beam connected to a downstream cylinder and placed stationary upstream cylinder before applying wind flow on both cylinders [23]. Optimum spacing between two cylinders was identified as being three times the circular cylinder diameter. Heganna and Joglekar [24] studied the forced vibration control analysis of the beam structure experimentally. The smart structure considered ﬁnite and ﬂat cantilever beam-like ﬁber structure with piezoelectric PZT patches attached on the surface.They also analyzed the vibration control effect of adjusting the location of PZT patches. Xie et al. proposed the fractional-order PD $\mu$ algorithm on the active vibration control of a lattice grid beam [25]. The vibration responses and dynamic model was based on the third-order shear deformation theory. A new active vibration control method is implemented using the PD $\mu$ algorithm in fractional order. Rimašauskienė et al. developed and experimentally analyzed the passive and active vibration control techniques for a thin-walled composite beam [26]. After the experiments with active vibration control, it was noticed that the vibration damping of the analyzed system was far higher at its resonant frequency in all cases. However, it could be inferred from the experimental results that passive and active vibration control systems were effective for thin-walled composite beam vibration amplitude control.

This work deals with the active vibration control analysis in smart composite beam structures by using ANSYS APDL codes. All nodes in the composite structure beam are subjected to the boundary condition (exactly 1260 boundary condition), and different proportional control gain is used to control vibration suppression to two sets of the composite beams, (0/90) sets and (0/45) sets, to reach the steady-state. From the obtained results, it can be concluded that the technique is able to dump the vibration of the composite beam structure.

2. Piezoelectric constitutive equations

Under the influence of stresses and voltage, the magnitude of the strain result and amount of electric charge accumulated by piezoelectric material specifies the actuation and sensing characteristics of these materials. Both constituent relationships under a piezoelectric material may be used to quantify these properties. Macroscopically, the piezoelectric materials display a field–strain relation [27]. The relationship is almost linear for the low electrical field, which can provide many advantages when using piezoelectric materials for system modeling and control. Nevertheless, the polarization saturates in a high electrical field, and domains extend and switch. This induces significant non-linear behavior that may be detrimental to the use of piezoelectric materials In the control implementation associated with a high electrical field [28]. In several applications of piezoelectric materials to continuous structures, the linear behavior is adopted despite the nonlinear behavior at a high electric field. The direct and converse effect of the piezoelectric phenomena involves an interaction between the mechanical behaviors of a material; Linear constitutive equations contain two mechanical variables and two electrical variables [27]. The equations governing the direct and converse piezoelectric effect in the matrix form are expressed respectively as

(1)

(2)

where ${\textstyle \left\{D\right\}}$ is the electric displacement, ${\textstyle {[e]}^{T}}$ is the transpose of the dielectric permittivity matrix ${\textstyle \left[e\right]}$ , ${\textstyle \left\{S\right\}}$ is the mechanical strain vector, ${\textstyle \left[{\epsilon }^{S}\right]}$ is the dielectric permittivity matrix at constant mechanical strain, ${\textstyle \left\{E\right\}\,}$ is the electric intensity vector, ${\textstyle \left\{T\right\}}$ is the mechanical stress vector and ${\textstyle {[c}^{E}]}$ is the matrix of elastic stiffness at constant electric field strength. The same relationship can be written in three other forms, depending on which variables are chosen to be independent. The direct relationship given by Eq.1 is normally used when modeling the sensing capability of the piezoelectric material, while the actuator capability is modeled using the converse effect relationship given by Eq. (2) [29].

3. Simulation method

In this section, ANSYS software simulates active vibration control in the smart composite. Figure 1 shows the block diagram of the obtained analysis. ANSYS/Multiphysics and ANSYS/Mechanical programs can be used for modeling structural fields with piezoelectric smart patches, whereas ${\textstyle {F}_{e}\,}$ is the force of vibration generated. The instant value of vibration producing force ${\textstyle {F}_{e}}$ can be described at each time step. In the analysis below, it is taken as ${\textstyle {F}_{o}}$ at ${\textstyle t\,=\Delta t}$ and zero for other time steps. At a piezoelectric sensor location, the strain ( $\varepsilon$ ) is calculated. Reference input for vibration cancellation is zero. ${\textstyle \,{K}_{p}}$ , ${\textstyle {K}_{v}}$ , and ${\textstyle {K}_{s}}$ are respectively, the variables the control, power, and sensor amplification factors. ${\textstyle {K}_{v}\,}$ is taken as 200, and ${\textstyle {K}_{p}\,}$ is changed in the analysis below. Only the proportional control is applied. The proportional constant for the actuator voltage ${\textstyle {V}_{a}}$ was obtained from the multiplication of ${\textstyle {K}_{p}{\,K}_{v}\,{K}_{s}}$ . Therefore, the results are not affected by changing the values of ${\textstyle {K}_{s}}$ , ${\textstyle {K}_{p}}$ and ${\textstyle {K}_{v}}$ and keeping their multiplication the same. To evaluate the performance of the vibration control, the calculated deflection at a location, ${\textstyle dt}$ , is observed.

Figure 1. Block diagram of the analysis

First, the structures of the beam type are considered. The configuration of the beam composite structure is shown in Figure 2.

Figure 2. Configuration for the structure of the beam type

3.1 Material Properties

Finite Element Analysis is used as the numerical modeling tool to simulate the dynamic characteristics of the composite structures by ANSYS macro codes, and the macro file has been written by using APDL. The first step in the coding process defined model geometry by using direct mesh generation of smart structure. In first our system, we used a composite smart beam of size (400 mm $\times$ 28 mm ) and various thickness (1.59 mm, 2.12mm and 3.18mm) according to number of laminates, with ( ${\textstyle {d}_{33}MFC-}$ $P1$ ) piezoelectric patch actuator type (M8528-P1). These typically consist of an active area (85mm $\times$ 28 mm) with (0.3 mm) thickness, and it is embedded on the beam from (10 mm) on the fixed point, as shown in Table 1 for all cases.

3.2 Simulation procedures

In the first step, the beam structures type is taken and considered. Figure 2 displays the configuration of the structure with sensor locations. The displacement value is taken from the displacement sensor, and similarly, the strain value is taken from the strain sensor as the feedback. Here, the product of ANSYS /Multiphysics (version 18) was applied to the structures of the smart composites model. The combination of the ANSYS modeling and solution into the control action codes were achieved by ANSYS [30].

The composite smart structure beam is studied, as indicated in Table 1.

Table 1. Dimensions and distances of Composite beam structure

Case No.	Composite layers	Dimensions of the structure (mm)'
1	(0/90/0)	$400\times 28\times 1.59$
2	(0/90/0/90)	$400\times 28\times 2.12$
3	(0/90/90/0)	$400\times 28\times 2.12$
4	(0/90/0/0/90/0)	$400\times 28\times 3.18$
5	(0/45/0)	$400\times 28\times 1.59$
6	(0/45/0/45)	$400\times 28\times 2.12$
7	(0/45/45/0)	$400\times 28\times 2.12$
8	(0/45/0/0/45/0)	$400\times 28\times 3.18$

The macro file creates nodes and finite elements to the structure. In the first part, composite materials have used element type SOLID185 layered structural elements for composite beam or plate part, and element type (SOLID5) elements with (MASS21) were used in the second structure part for the piezoelectric part in the structure [31]. Overall, these macro codes as shown in Table 2 indicate the defined element types of various materials.

Table 2. Element definition codes

/PREP7	! start the preprocessing
! Define Element type #1
ET,1,185	! chooses SOLID185 element for analysis composite material
KEYOPT,1,3,1	! Key option #3 = 1 LAYERED SOLID ELEMENTS
KEYOPT,1,2,2	! Key option #2 = 2 ENHANCED STRAIN FORMULATION
KEYOPT,1,8,1	! Key option #8 = 1 STORE DATA FOR ALL LAYERS

! Define Element type #2
ET,2,solid5,3	! chooses SOLID5 element for analysis piezoelectric material
ET,3,mass21

Therefore, an essential next step to confirming in three-dimensional composite material properties is assigned. As shown in Figure 2, the structure is divided into two parts. The first part deals with composite material, and the other part is actuator and all properties, as shown in Table 3, and all that has been gained from [32,33,34].

Table 3. Material properties of MFC and CFRP

Material Properties	CFRP
Young’s modulus ${\textstyle {\boldsymbol {E}}_{11}}$	140 GPa	30.34 GPa
Young’s modulus ${\textstyle {\boldsymbol {E}}_{22}}$	l8.75 GPa	15.86 GPa
Young’s modulus ${\textstyle {\boldsymbol {E}}_{33}}$	8.75 GPa	9.03 GPa
Shear modulus ${\textstyle {\boldsymbol {G}}_{12}}$	7.2 GPa	5.52 GPa
Poisson’s ratio ${\textstyle {\boldsymbol {v}}_{12}}$	0.3	0.31
Poisson’s ratio ${\textstyle {\boldsymbol {v}}_{13}}$	0.3	0.30
Density $\rho$	1550 kg/m $^{3}$	5020 kg/m $^{3}$
Piezoelectric Stress constant ( ${\textstyle {\boldsymbol {e}}_{32}}$ )		-3.47 C/m $^{2}$
Piezoelectric Stress constant ( ${\textstyle {\boldsymbol {e}}_{33}}$ )		12.36 C/m $^{2}$

The ANSYS macro of the material properties is given in Table 4. The first material is represented by composite material (Carbon fiber reinforced polymer), and the second material is represented by the piezoelectric (MFC) actuator material.

Table 4. Composite material properties codes

/PREP7		! Define Material 1# : ${\textstyle E1=140\times {10}^{9}\,Pa}$
MP,EX,1,140E9
MP,EY,1,8.75E9		! Define Material 1# : ${\textstyle E2=8.75\times {10}^{9}\,Pa}$
MP,EZ,1,8.75E9		! Define Material 1# : ${\textstyle EZ=EY\,;\,E3=8.75\times {10}^{9}\,Pa}$
MP,GXY,1,7.2E9		! Define Material 1# : shear modulus xy-plane ${\textstyle {G}_{xy}=}$ $7.2\times {10}^{9}\,Pa$
MP,GYZ,1,6.21E9		! Define Material 1# : shear modulus xz-plane ${\textstyle {G}_{xy}=6.2\times {10}^{9}\,Pa}$
MP,GXZ,1,7.2E9		! Define Material 1# : shear modulus yz-plane ${\textstyle {G}_{xy}=7.2\times {10}^{9}\,Pa}$
MP,PRXY,1,0.3		! Define Material 1# : Poisson coefficient ${\textstyle {\upsilon }_{1}=0.3}$
MP,PRYZ,1,0.3		! Define Material 1# : Poisson coefficient ${\textstyle {\upsilon }_{2}=0.3}$
MP,PRXZ,1,0.3		! Define Material 1# : Poisson coefficient ${\textstyle {\upsilon }_{3}=0.3}$
MP,DENS,1,1550		! Define Material 1# : Density for composite material

The number of composite layers, the thickness of layers, and orientation angles of fibers are assigned as macro APDL ANSYS codes for the first material, as given below in Table 5.

Table 5. Orientation angles’ Composite layers codes

SECTYPE,1,SHELL	set #1, [0/90/0]
SECDATA,0.53e-3,1,0	! 1st LAYER 1: Thk=0.53e-3m, THETA = 0°
SECDATA,0.53e-3,1,90	! 2nd LAYER 2: Thk=0.53e-3m, THETA = 90°
SECDATA,0.53e-3,1,0	! 3rd LAYER 4: Thk=0.53e-3m, THETA = 0°

The second material piezoelectric actuator properties are given in Table 6 as macro ANSYS.

Table 6. MFC properties codes

! Define dielectric constant
MP,DENS,3,5020	! Density for MFC piezoelectric material
MP,PERX,3,15.586E-9	! Permittivity in X direction
MP,PERY,3,15.586E-9	! Permittivity in Y direction
MP,PERZ,3,11.79E-9	! Permittivity in Z direction
TB,PIEZ,3	! Define piezoelectric table
TBDATA,16,0	! e16 piezoelectric constant
TBDATA,14,0	! e25 piezoelectric constant
TBDATA,6,-3.47	! e32 piezoelectric constant
TBDATA,9,12.36	! e33 piezoelectric constant
TB,ANEL,3	! Define structural table

After all dimensions and properties of both composite and piezoelectric properties are defined, the type of element is assigned for each type of material. SOLID185 layered structural elements are assigned for the composite part, and SOLID5 elements are used for the piezoelectric part of the structure. Then, the static solution is analyzed to evaluate the initial deflections of each node and element. The model is discretized into finite elements model, as shown in Figure 3.

Figure 3. Finite element model

The first three mode shapes and natural frequencies of the composite plates were obtained using APDL ANSYS macro codes. The obtained results were shown in Figure 4. The first three natural frequencies evaluated by the modal shape analysis method using ANSYS (15.665 Hz), (94.17Hz), and (161.63Hz), respectively, for the first case considered in the present study, are shown in Figure 4.


1st mode shape	2nd mode shape	3rd mode shape
Figure 4. Mode shapes of the first case

Fixed boundary conditions are described for the nodes ${\textstyle \,(x=0)}$ and given in Table 7 as macro ANSYS.

Table 7. Boundary conditions of composite beam codes

CP,1,volt,all	! set #1 of coupled degrees of freedom.

The degrees of freedom of the top and bottom surfaces for piezo actuator and sensor are coupled simultaneously with volt by APDL ANSYS “CP” command. Modal analysis was performed to determine parameters from vibration data measures like the inherent dynamic characteristics of the system in forms of natural frequencies, mode shapes and damping factor, and so to determine time step [35].

Householder method (reduced method) can be used for the beam structures which have the solid coupled field in ANSYS. The time step is the incremental change in time and the time step is chosen by equation formula ${\textstyle \Delta t\,=1/(20\,{f}_{h})}$ , where ${\textstyle {f}_{h}\,}$ is the highest natural frequency to be considered. Table 8 shows the three natural frequencies for the undamped system. The first mode is considered to calculate the time step and ${\textstyle \Delta t}$ is 0.0032, 0.0034, 0.0025, 0.00183, 0.00318, 0.0032, 0.0025, and 0.0018, for cases 1, 2, 3, 4, 5, 6, 7, and 8, respectively.

In this study, The strain feedback (piezoelectric sensor) is calculated at the selected sensor location, multiplied by ${\textstyle {Ks}_{1}}$ to convert the strain feedback to voltage and then subtracted from zero to calculate the error signal.

The displacement feedback is also calculated at the selected sensor location, multiplied by ${\textstyle {Ks}_{2}}$ in order to convert the displacement feedback to voltage, and also subtracted from zero to evaluate the error signal between the input and output.

Table 8. Natural frequencies of the system

Case No.	Composite layers	First (Hz)	Second (Hz)	Third (Hz)
1	(0/90/0)	15.665	94.17	161.63
2	(0/90/0/90)	14.888	87.674	193.47
3	(0/90/90/0)	19.795	120.39	196.81
4	(0/90/0/0/90/0)	27.265	165.45	217.71
5	(0/45/0)	15.701	94.379	170.62
6	(0/45/0/45)	15.589	92.137	200.82
7	(0/45/45/0)	19.945	121.29	204.18
8	(0/45/0/0/45/0)	27.697	168.02	222.82

The difference between the input reference and the sensor signal is defined as the error signal. The error value is multiplied by ${\textstyle {Kp}_{1}}$ for strain feedback and ${\textstyle Kp\,}$ for displacement feedback, and then both results are multiplied by the amplifier factor ${\textstyle Kv}$ to determine ${\textstyle Va}$ at a given time step.

The macro part that enables the calculations for the analysis of the closed loop for ${\textstyle t>\Delta t}$ is given in Table 9.

Table 9. Sensor feedback codes

do,tt,2dt,ts,dt	!time loops
*get,uz-tip,node,nldm,u,z	!sensor feedback
*err=ref-ksuz-tip**	!error signal
*Vact=kpkverr*	!actuator voltage
d,nv,volt,Vact	!Appling actuator volt
time,tt
solve
*enddo

4. Results and discussion

There are two methods to deal with active vibration control. The first method is displacement feedback signals by using the displacement sensor, and the other method is strain feedback signals by using a piezoelectric strain sensor. Displacement feedback signals have been used for the displacement to suppression of vibration.

Classical controller ( $P$ ) is the proportional control and gain ( $kp$ ) is the only variable to be changed. The present study has been used for (2,4,6, and 8) $kp$ value. The open-loop and closed-loop responses of the tip displacement are shown in Figures 5(a), (b), (c) and (d), respectively, by using APDL ANSYS. From the graph, it is found that the open-loop response vibrates a lot at a time of more than 5 sec in first Figure 5(a), 4 sec in second Figure 5(b), 3.8 in third Figure 5(c), and 1.25 sec in the last Figure 5(d) [36]. In the first case, composite laminates consist of three-ply ( $0/90/0$ ) closed-loop responses by using classical controller when we took into consideration only proportional part effects in PID control, which is shown settle at a time 1.4 sec, 0.85 sec, 0.75 sec and 0.7 sec, respectively, proportional gain ( $kp$ ) is 2, 3, 6, 8, which demonstrates that three recent responses of the system will be near each other. This reduction in the time compared with the first case is due to the increase in the number of composite layers, and, at the same time, an increase in the strength properties of the beam is observed.


(a) CFRP (0/90/0)	(b) CFRP (0/90/0/90)

(c) CFRP (0/90/90/0)	(d) CFRP (0/90/0/90/0/90
Figure 5. Open-loop and closed-loop responses of a composite beam of tip displacement

However, the evidence for this relationship is conclusive and has been retrieved from PID rules; with increasing the proportional range, the vibration converges slowly and settling time increases. On the other side, the controlling signal has been decreasing by increasing the value of the proportional constant that takes part in PID and neglected integral and derivative parts. Whereas in the second case, composite laminates consist of four plies (0/90/0/90), the closed-loop responses or all proportional gain ( $kp$ ) shown as in Figure 5(b) is close to the results of the first case because the effect of the fourth ply is smaller (90 degree) composite ply with the beam in the fourth case, which demonstrates that control responses of the system will be near each other and this reduction in the time compared with others is due to the increase in the number of composite layers and at the same time increases the strength properties of the beam. Based on the results, it has been noted that an increase in laminates of composite plate decreases the vibration and time stability of the system.

For amplitude comparison between two sets of composite materials beam structures, first set consist of (0-90) degrees of laminates and another set consist of (0-45) degrees of laminates. This based on the responses that were simulated with high accuracy using the proportional gain formulation. As shown in Figures 5 and 6, all curves have nearly the same shape, but speed to reach the stable case in the second set is smaller than the first set. The reason is attributed to the greater flexural strength for the second compared to the former.


(a) CFRP (0/45/0)	(b) CFRP (0/45/0/45)

(c) CFRP (0/45/45/0)	(d) CFRP (0/45/0/45/0/45)
Figure 6. Open-loop and closed-loop responses of the composite beam of tip displacement

Piezoelectric MFC has a maximum voltage, so it is peak to peak voltage. While peak voltage is the voltage from the zero horizontal axis line to the top of the waveform, the peak to peak voltage is the voltage of the one cycle of a wave, all the way from the peak of the negative side to the peak of the positive side [37]. The maximum voltage for MFC to be applied is (2000 volt) for the piezoelectric actuator. The uncontrolled and controlled responses obtained by actuation voltage can be seen in Figures 7(a), (b), (c), (d), (e), (f), (g), and (h). From the graph, it has been observed that in the open-loop, the response was nil because in the open-loop, there is no link between vibrate response and actuator.

On the other hand, in all cases, when the control signal applied with the gain to a system increases, the rise time generally decreases as shown in Figure 7, where the actuator voltages for different control gain for all cases are observed to increase with a proportional constant. Meanwhile, $kp=8$ has the highest value of actuator voltage observed in all cases [38].

As can be seen from Figures 7(b), (c), (d), (f), (g), and (h), the minimum value of actuator voltage is (-500V) and smaller than this volt by APDL-ANSYS macro file has been changed to (-500 Volt) because the sinusoidal voltage of piezoelectric MFC-P1 is worked between (-500 volt) to (+1500 Volt) voltage cycle at 0.1 Hz was output by a voltage amplifier [39].


(a) CFRP (0/90/0)	(b) CFRP (0/90/0/90)

(c) CFRP (0/90/90/0)	(d) CFRP (0/90/0/0/90/0)

(e) CFRP (0/45/0)	(f) CFRP (0/45/0/45)

(g) CFRP (0/45/45/0)	(h) CFRP (0/45/0/0/45/0)
Figure 7. Open-loop and closed-loop responses of the composite beam of the actuation voltage

Another important finding was the effect of fiber orientations of composite laminate layers on active vibrations, as shown in Figure 8. The composite beams composed of laminates (0/90/0), (0/90/0/90), (0/90/90/0), and (0/90/0/0/90/0) reach stability faster than the composite beams composed of laminates (0/45/0), (0/45/0/45), (0/45/45/0), and (0/45/0/0/45/0) because the first set has higher relaxation modulus than second sets [40].


(a)	(b)

(c)	(d)
Figure 8. Comparision between two sets of similar composite beam

5. Conclusions

A FEM for piezoelectric with composite beam structure was proposed to simulate model dynamic, and the electrical model of piezoelectric which permitted the simulations of type ( ${\textstyle {d}_{33}-}$ $MFC$ ) effects. The performance of the composite beam structure is improved when piezoelectric type M8528(P1) is used. It is shown that the classic controller provides in increasing control gain ( ${\textstyle kp}$ ) and better active damping than low control gain, especially at ${\textstyle kp=8}$ .

The active control capability of such a composite beam is analyzed using a simple control strategy. Overall, these controlled responses show that, although all different composite laminate plates are subjected to the same intensity of mechanical load. The observed decrease in the vibration of the composite plate might be explained by the increase in the proportional gain in PID control on all cases of the composite laminate plate, while at the same time decreasing the voltage actuator of piezoelectric. The model can be used for more investigations with general composition, applied loads, and another control method. The final objectives of analysis and modeling are the minimization of weight (or mass) of the beam structure under behavior constraints, increasing the performance of the piezoelectric actuator, and maximization of natural frequencies of the vibration modes. However, given the immersed short time for each case, caution must be taken. The difference between two similar different in orientation laminate, two refined finite element models based 0/90 laminate set and 0/45 laminate set, for the analysis of laminated beam with embedded surface bonded piezoelectric laminae have been presented. The evidence from this study suggests the composite beams composed from (0/90) stacking sequence sets had reached stability state faster than the composite beams composed from (0/45) stacking sequencesets.

References

[1] Sun L., Li W., Wu Y., Lan Q. Active vibration control of a conical shell using piezoelectric ceramics. Journal of Low Frequency Noise, Vibration and Active Control, 36(4):366-375, 2017.

[2] Baz A.M.. Active and passive vibration damping. John Wiley & Sons, 2019.

[3] Kouider B., Polat A. Numerical investigation of active vibration control of composite plate using PID controller. International Conference on Advances and Innovations in Engineering, 2017.

[4] Symans M.D., Constantinou M.C. Semi-active control systems for seismic protection of structures: A state-of-the-art review. Engineering structures, 21(6):469-487, 1999.

[5] Kumar V. Experimental investigation on active vibration control of smart composite shell structures. PhD Thesis, Jawaharlal Nehru Technological University Hyderabad, 2011.

[6] Manning W.J., Plummer A.R., Levesley M.C. Vibration control of a flexible beam with integrated actuators and sensors. Smart Materials and Structures, 9(6):932, 2000.

[7] Gaudenzi P., Carbonaro R., Benzi E. Control of beam vibrations by means of piezoelectric devices: theory and experiments. Composite Structures, 50(4):373-379, 2000.

[8] Bruant I., Coffignal G., Lene F., Verge M. Active control of beam structures with piezoelectric actuators and sensors: modeling and simulation. Smart Materials & Structures, 10(2):404, 2001.

[9] Singh S.P., Pruthi H.S., Agarwal V.P. Efficient modal control strategies for active control of vibrations. Journal of Sound and Vibration, 262(3):563-575, 2003.

[10] Xu S.X., Koko T.S. Finite element analysis and design of actively controlled piezoelectric smart structures. Finite Elements in Analysis and Design, 40(3):241-262, 2004.

[11] Lim Y.H. Finite-element simulation of closed loop vibration control of a smart plate under transient loading. Smart Materials and Structures, 12(2):272, 2003.

[12] Quek S.T., Wang S.Y., Ang, K.K. Vibration control of composite plates via optimal placement of piezoelectric patches. Journal of Intelligent Material Systems and Structures, 14(4-5):229-245, 2003.

[13] Xianmin Z., Changjian S., Erdman A.G. Active vibration controller design and comparison study of flexible linkage mechanism systems. Mechanism and Machine Theory, 37(9):985-997, 2002.

[14] Halim D., Moheimani S.R. An optimization approach to optimal placement of collocated piezoelectric actuators and sensors on a thin plate. Mechatronics, 13(1):27-47, 2003.

[15] Bendine K., Boukhoulda F.B., Haddag B., Nouari M. Active vibration control of composite plate with optimal placement of piezoelectric patches. Mechanics of Advanced Materials and Structures, 26(4):341-349, 2019.

[16] Malgaca L. Integration of active vibration control methods with finite element models of smart laminated composite structures. Composite Structures, 92(7):1651-1663, 2010.

[17] Takacs G., Rohal-Ilkiv B. Direct closed-loop active vibration control system prototype in ANSYS. In INTER-NOISE and NOISE-CON Congress and Conference Proceedings, 244(1):1-12, Institute of Noise Control Engineering, April 2012.

[18] Vashist S.K., Chhabra D. Optimal placement of piezoelectric actuators on plate structures for active vibration control using genetic algorithm. In Active and Passive Smart Structures and Integrated Systems, International Society for Optics and Photonics, 9057:9, 2014.

[19] Yavuz Ş., Malgaca L., Karagülle H. Analysis of active vibration control of multi-degree-of-freedom flexible systems by Newmark method. Simulation Modelling Practice and Theory, 69:136-148, 2016.

[20] Liu X., Cai G., Peng F., Zhang H. Piezoelectric actuator placement optimization and active vibration control of a membrane structure. Acta Mechanica Solida Sinica, 31(1):66-79, 2018.

[21] Van Tham V., Quoc T.H., Tu T.M. Optimal placement and active vibration control of composite plates integrated piezoelectric sensor/actuator pairs. Vietnam Journal of Science and Technology, 56(1):113, 2018.

[22] Baghaee M., Farrokhabadi A., Jafari-Talookolaei R.A. Modeling, analysis, and control of MFC sandwiched laminate panel flutter with general layups and arbitrary boundary conditions. Composite Structures, 223:110940, 2019.

[23] Usman M., Hanif A., Kim I.H., Jung H.J. Experimental validation of a novel piezoelectric energy harvesting system employing wake galloping phenomenon for a broad wind spectrum. Energy, 153:882-889, 2018.

[24] Heganna S.S., Joglekar J.J. Active vibration control of smart structure using PZT patches. Procedia Computer Science, 89:710-715, 2016.

[25] Xie C., Wu Y., Liu Z. Modeling and active vibration control of lattice grid beam with piezoelectric fiber composite using fractional order PDμ algorithm. Composite Structures, 198:126-134, 2018.

[26] Rimašauskienė R., Jūrėnas V., Radzienski M., Rimašauskas M., Ostachowicz W. Experimental analysis of active–passive vibration control on thin-walled composite beam. Composite Structures, 223:110975, 2019.

[27] IEEE. An American National Standard: IEEE Standard on Piezoelectricity. IEEE Transactions on Sonics and Ultrasonics, (2):1-55, 1984.

[28] Abuzaid A., Hrairi M., Dawood M.S.I. Survey of active structural control and repair using piezoelectric patches. Actuators 4(2):77-98, Multidisciplinary Digital Publishing Institute, June 2015.

[29] Heywang W., Lubitz K., Wersing W. Eds. Piezoelectricity: evolution and future of a technology. Springer Science & Business Media, Springer Series in Materials Science, Vol. 114, 2008.

[30] ANSYS, ANSYS Mechanical APDL Basic Analysis Guide/ANSYS help, release 15.0 2013. Inc. Canonsburg, PA, 2017.

[31] Kamel M.A., Ibrahim K., Ahmed A.E.M. Vibration control of smart cantilever beam using finite element method. Alexandria Engineering Journal, 58(2):591-601, 2019.

[32] Trindade M.A., Benjeddou A. Finite element characterisation of multilayer d31 piezoelectric macro-fibre composites. Composite Structures, 151:47-57, 2016.

[33] Vijay Kumar G., Raja S., Prasanna K.B., Sudha V. Finite element analysis and vibration control of a deep composite cylindrical shell using MFC actuators. Smart Materials Research, 2012.

[34] Smart Materials Corporation (Accessed Nov 2019). Actuator and sensor ,Available online: https://www.smart-material.com/MFC-product-P1.html

[35] He J., Fu Z.F., Modal Analysis. Linacre House, Jorden Hill, Oxford, 2001.

[36] Huang B., Kim H.S., Youn B.D. Active vibration control of smart composite laminates with partial debonding of actuator. International Journal of Precision Engineering and Manufacturing, 16(4):831-840, 2015.

[37] Ian Sinclair. Passive Components for Circuit Design. Linacre House, Jordan Hill, Oxford OX2 8DP, 2004.

[38] Garcia-Segura T., Yepes V., Marti J.V., Alcala J. Optimization of concrete I-beams using a new hybrid glowworm swarm algorithm. Latin American Journal of Solids and Structures, 11(7):1190-1205, 2014.

[39] Wu M., Yuan X., Luo H., Chen H., Chen C., Zhou K., Zhang D. Enhanced actuation performance of piezoelectric fiber composites induced by incorporated BaTiO3 nanoparticles in epoxy resin. Physics Letters A, 381(19):1641-1647, 2017.

[40] Das R. Stress relaxation of ehermostting prepreg composites during cure. MSc. Thesis, Bangladesh University of Engineering and Technology, 2008.

Latest revision as of 14:18, 24 February 2021

Abstract

1. Introduction

2. Piezoelectric constitutive equations

3. Simulation method

3.1 Material Properties

3.2 Simulation procedures

4. Results and discussion

5. Conclusions

References

Document information

Document Score

Share this document

Keywords

claim authorship

@@ Line 10: / Line 10: @@
 Mechanics and Mechatronics Department, College of Engineering, Salahaddin University, Erbil, Iraq. Email: [mailto:azad.saber@su.edu.krd azad.saber@su.edu.krd]
 -->
 ==Abstract==
-<span style="text-align: center; font-size: 75%;">Piezoelectric Macro-Fiber Composite (MFC) utilization is increasing in engineering fields due to its strong actuation forces and high flexibility. In this paper, piezoelectric (</span> <math display="inline">{d}_{33}MFC</math><span style="text-align: center; font-size: 75%;">) type (M8528-P1) patches are applied for active control of a cantilever composite beam. A linear coupled ﬁnite element model for piezoelectric MFC actuation of the composite beam was developed based on the APDL-ANSYS codes by using linear piezoelectric constitutive equations to study smart composite beam behavior in open and closed-loop cases. Active vibration control response of laminated composite beams with various stacking sequence configurations was examined, and the results were compared. In this, the proportional type (kp) of control algorithms is utilized,  and in the summary of work refined two sets of finite element models. The first set their laminates to have an orientation (0 with 90), and others have an orientation (0 with 45). When the control signal applied with the gain (kp) to a system increases, the rise time generally decreases, the actuator voltages for different control gain for all cases are observed with the increase in the proportional constant. The findings of this study indicate that the composite beams composed from (0/90)stacking sequence sets had reached stability state faster than the composite beams composed from (0/45) stacking sequence sets.</span>
+Piezoelectric Macro-Fiber Composite (MFC) utilization is increasing in engineering fields due to its strong actuation forces and high flexibility. In this paper, piezoelectric (<math display="inline">{d}_{33}MFC</math>) type (M8528-P1) patches are applied for active control of a cantilever composite beam. A linear coupled ﬁnite element model for piezoelectric MFC actuation of the composite beam was developed based on the APDL-ANSYS codes by using linear piezoelectric constitutive equations to study smart composite beam behavior in open and closed-loop cases. Active vibration control response of laminated composite beams with various stacking sequence configurations was examined, and the results were compared. In this, the proportional type (<math>kp</math>) of control algorithms is utilized,  and in the summary of work refined two sets of finite element models. The first set their laminates to have an orientation (0 with 90), and others have an orientation (0 with 45). When the control signal applied with the gain (<math>kp</math>) to a system increases, the rise time generally decreases, the actuator voltages for different control gain for all cases are observed with the increase in the proportional constant. The findings of this study indicate that the composite beams composed from (0/90)stacking sequence sets had reached stability state faster than the composite beams composed from (0/45) stacking sequence sets.
-==Keywords:==
+'''Keywords''': MFC piezoelectric, smart composite beam, FEM (APDL Codes), active vibration
-MFC piezoelectric, smart composite beam, FEM (APDL Codes), Active vibration.
 ==1. Introduction==
@@ Line 23: / Line 20: @@
 Vibration is present in countless daily life applications, and more often, vibration is highly undesirable. Undesirable vibration may decrease the performance of the product and cause safety or economic issues [1]. Vibration damping includes active, passive, and hybrid damping techniques. It is important to understand that passive damping can be very efficient in damping out high-frequency excitation, while active damping can be used to control low-frequency vibrations [2]. Active vibration-reducing techniques turn out to be a reality for a wide variety of applications [3]. An active method integrates sensors and actuators with a flexible structure, operated by a control scheme. An essential issue in active control systems is sensor and actuator selection [4]. The ever-increasing use of smart composite materials in advanced areas such as aerospace, automotive, and sports equipment has motivated researchers to explore behavior and performance characteristics [5].
-There are some previous studies reported in this work. The simple cantilever beams active vibration control is investigated in [6-10]. Piezoelectric patches are mounted on the beams as actuators. With the purpose of sensing the vibration level, a strain-gage or another piezoelectric patch can be used. The identification system and placement control technique was utilized in [6]. The piezo patches and beam structure’s finite element model is constructed, and the closed-loop control is applied in [7,8]. Also, Singh et al. [9] utilized the model of the beam and piezo-patches finite element; however, modal control techniques were employed. Xu and Koko [10] used the commercial finite element package ANSYS to state the findings. The control design was conducted in the space state that was set in the analysis of the finite element model. For the control design, MATLAB as a system toolbox was utilized in their experiment. The sensors/actuators influence site was investigated. It was noticed that the site near the fixed end for the vibration control was better. Lim [11] investigated several modes for vibration control of a fixed square plate by using devices of location discrete sensor/actuator at maximum strain points. Quek et al. [12] introduced a typical placement method of pairs of piezoelectric sensor /actuator for the vibration control sealed compound plates. Xiamin et al. [13] used four-bar linkage to investigate active vibration control. The model of the finite element is applied and the reduced mode, robust  <math display="inline">{H}_{\infty }</math>, and standard H∞ control methods are analyzed. Halim and Moheimani [14] investigated a set of proposals to a criterion for ideal placement of a collected piezoelectric sensor-actuator pair on their supported plate by using modal and locative controllability. Bendine et al. [15] studied the modeling of the active vibration control of a smart functionally graded (FGM) beam based on a higher-order shear deformation theory. The locations of piezoelectric actuators and sensors have a great influence on the reliability and control system efficiency. Malgaca, L [16] studied optimal parameters with respect to positive position feedback used in the real-time control algorithm determined with the objective functions. Flexible single-link composite manipulator with [45/-45] and [0/90] lay-ups are performed to reduce endpoint vibrations. Takacs, G. and Rohal-Ilkiv [17] proposed a direct closed-loop prototype method in the active vibration control system modeled in the finite element method using APDL language for the clamped-free aluminum beam piezoelectric actuation using a linear quadratic controller (LQ) with a Kalman-ﬁlter.Vashist and Chhabra [18] found the optimal piezoelectric actuator’s locations by using an integer coded genetic algorithm GA on their plate where the singular value decomposition’s control matrix was used as a fitness function. Similar results of the frequency responses between experimental and Finite Element Method were found. Yavuz et al.[19] analyzed a system of active vibration control’s closed-loop with the purpose of suppressing the vibrations of the end-point. They established the system’s mathematical model by using the Lagrange equations.  The digital solution is achieved through mixing the Newmark technique with the control action. They applied PID control action with the purpose of finding the value of the actuator signal at the appropriate time. Liu X. et al.[20] studied the piezoelectric actuators’ placement optimization and the membrane structure’s active vibration control. They designed the classical linear quadratic regulator controllers with the aim of suppressing the undesired vibration. Simulation findings indicated that the piezoelectric actuators’ optimal locations are influenced by the actuator’s additional mass and stiffness. Van et al. [21] developed a finite element model according to the theory of first-order shear deformation for active vibration control and optimal placement of laminated composite plates with linked pairs of distributed piezoelectric sensor/actuator. Baghaee et al. [22] developed the aero elastic panel flutter of the laminate composite plate to include two piezoelectric (MFC-type). They determined MFC actuator’s input voltage through the velocity or proportional feedback control algorithms according to the sensor output. Usman et al. [23] used piezoelectric MFC energy harvesting on the cantilever beam connected to a downstream cylinder and placed stationary upstream cylinder before applying wind flow on both cylinders. Optimum spacing between two cylinders was identified as being three times the circular cylinder diameter. Heganna and Joglekar [24] studied the forced vibration control analysis of the beam structure experimentally. The smart structure considered ﬁnite and ﬂat cantilever beam-like ﬁber structure with piezoelectric PZT patches attached on the surface.They also analyzed the vibration control effect of adjusting the location of PZT patches. Xie et al. [25] proposed the fractional-order PD<sup>μ</sup> algorithm on the active vibration control of a lattice grid beam. The vibration responses and dynamic model was based on the third-order shear deformation theory. A new active vibration control method is implemented using the PDμ algorithm in fractional order.Rimašauskienė et al. [26] developed and experimentally analyzed the passive and active vibration control techniques for a thin-walled composite beam. After the experiments with active vibration control, it was noticed that the vibration damping of the analyzed system was far higher at its resonant frequency in all cases. However, it could be inferred from the experimental results that passive and active vibration control systems were effective for thin-walled composite beam vibration amplitude control.
+There are some previous studies reported in this work. The simple cantilever beams active vibration control is investigated in [6-10]. Piezoelectric patches are mounted on the beams as actuators. With the purpose of sensing the vibration level, a strain-gage or another piezoelectric patch can be used. The identification system and placement control technique was utilized in [6]. The piezo patches and beam structure’s finite element model is constructed, and the closed-loop control is applied in [7,8]. Also,  utilized the model of the beam and piezo-patches finite element; however, modal control techniques were employed Singh et al. [9]. Xu and Koko  used the commercial finite element package ANSYS to state the findings [10]. The control design was conducted in the space state that was set in the analysis of the finite element model. For the control design, MATLAB as a system toolbox was utilized in their experiment. The sensors/actuators influence site was investigated. It was noticed that the site near the fixed end for the vibration control was better. Lim investigated several modes for vibration control of a fixed square plate by using devices of location discrete sensor/actuator at maximum strain points [11]. Quek et al.  introduced a typical placement method of pairs of piezoelectric sensor /actuator for the vibration control sealed compound plates [12]. Xiamin et al.  used four-bar linkage to investigate active vibration control [13]. The model of the finite element is applied and the reduced mode, robust  <math display="inline">{H}_{\infty }</math>, and standard H∞ control methods are analyzed. Halim and Moheimani investigated a set of proposals to a criterion for ideal placement of a collected piezoelectric sensor-actuator pair on their supported plate by using modal and locative controllability [14]. Bendine et al. studied the modeling of the active vibration control of a smart functionally graded (FGM) beam based on a higher-order shear deformation theory [15]. The locations of piezoelectric actuators and sensors have a great influence on the reliability and control system efficiency. Malgaca  studied optimal parameters with respect to positive position feedback used in the real-time control algorithm determined with the objective functions [16]. Flexible single-link composite manipulator with [45/-45] and [0/90] lay-ups are performed to reduce endpoint vibrations. Takacs and Rohal-Ilkiv [17] proposed a direct closed-loop prototype method in the active vibration control system modeled in the finite element method using APDL language for the clamped-free aluminum beam piezoelectric actuation using a linear quadratic controller (LQ) with a Kalman-ﬁlter. Vashist and Chhabra [18] found the optimal piezoelectric actuator’s locations by using an integer coded genetic algorithm GA on their plate where the singular value decomposition’s control matrix was used as a fitness function. Similar results of the frequency responses between experimental and Finite Element Method were found. Yavuz et al.  analyzed a system of active vibration control’s closed-loop with the purpose of suppressing the vibrations of the end-point [19]. They established the system’s mathematical model by using the Lagrange equations.  The digital solution is achieved through mixing the Newmark technique with the control action. They applied PID control action with the purpose of finding the value of the actuator signal at the appropriate time. Liu X. et al. studied the piezoelectric actuators’ placement optimization and the membrane structure’s active vibration control [20]. They designed the classical linear quadratic regulator controllers with the aim of suppressing the undesired vibration. Simulation findings indicated that the piezoelectric actuators’ optimal locations are influenced by the actuator’s additional mass and stiffness. Van et al.  developed a finite element model according to the theory of first-order shear deformation for active vibration control and optimal placement of laminated composite plates with linked pairs of distributed piezoelectric sensor/actuator [21]. Baghaee et al.  developed the aero elastic panel flutter of the laminate composite plate to include two piezoelectric (MFC-type) [22]. They determined MFC actuator’s input voltage through the velocity or proportional feedback control algorithms according to the sensor output. Usman et al.  used piezoelectric MFC energy harvesting on the cantilever beam connected to a downstream cylinder and placed stationary upstream cylinder before applying wind flow on both cylinders [23]. Optimum spacing between two cylinders was identified as being three times the circular cylinder diameter. Heganna and Joglekar [24] studied the forced vibration control analysis of the beam structure experimentally. The smart structure considered ﬁnite and ﬂat cantilever beam-like ﬁber structure with piezoelectric PZT patches attached on the surface.They also analyzed the vibration control effect of adjusting the location of PZT patches. Xie et al. proposed the fractional-order PD<math>\mu</math> algorithm on the active vibration control of a lattice grid beam [25]. The vibration responses and dynamic model was based on the third-order shear deformation theory. A new active vibration control method is implemented using the PD<math>\mu</math> algorithm in fractional order. Rimašauskienė et al.  developed and experimentally analyzed the passive and active vibration control techniques for a thin-walled composite beam [26]. After the experiments with active vibration control, it was noticed that the vibration damping of the analyzed system was far higher at its resonant frequency in all cases. However, it could be inferred from the experimental results that passive and active vibration control systems were effective for thin-walled composite beam vibration amplitude control.
 This work deals with the active vibration control analysis in smart composite beam structures by using ANSYS APDL codes. All nodes in the composite structure beam are subjected to the boundary condition (exactly 1260 boundary condition), and different proportional control gain is used to control vibration suppression to two sets of the composite beams, (0/90) sets and (0/45) sets, to reach the steady-state. From the obtained results, it can be concluded that the technique is able to dump the vibration of the composite beam structure.
-==2. Piezoelectric Constitutive Equations==
+==2. Piezoelectric constitutive equations==
-Under the influence of stresses and voltage, the magnitude of the strain result and amount of electric charge accumulated by piezoelectric material specifies the actuation and sensing characteristics of these materials. Both constituent relationships under a piezoelectric material may be used to quantify these properties. Macroscopically, the piezoelectric materials display a field–strain relation [27]. The relationship is almost linear for the low electrical field, which can provide many advantages when using piezoelectric materials for system modeling and control. Nevertheless,  the polarization saturates in a high electrical field, and domains extend and switch. This induces significant non-linear behavior that may be detrimental to the use of piezoelectric materials In the control implementation associated with a high electrical field [28]. In several applications of piezoelectric materials to continuous structures,  the linear behavior is adopted despite the nonlinear behavior at a high electric field. The direct and converse effect of the piezoelectric phenomena involves an interaction between the mechanical behaviors of a material; Linear constitutive equations contain two mechanical variables and two electrical variables [27] . The equations governing the direct and converse piezoelectric effect in the matrix form are expressed respectively as
+Under the influence of stresses and voltage, the magnitude of the strain result and amount of electric charge accumulated by piezoelectric material specifies the actuation and sensing characteristics of these materials. Both constituent relationships under a piezoelectric material may be used to quantify these properties. Macroscopically, the piezoelectric materials display a field–strain relation [27]. The relationship is almost linear for the low electrical field, which can provide many advantages when using piezoelectric materials for system modeling and control. Nevertheless,  the polarization saturates in a high electrical field, and domains extend and switch. This induces significant non-linear behavior that may be detrimental to the use of piezoelectric materials In the control implementation associated with a high electrical field [28]. In several applications of piezoelectric materials to continuous structures,  the linear behavior is adopted despite the nonlinear behavior at a high electric field. The direct and converse effect of the piezoelectric phenomena involves an interaction between the mechanical behaviors of a material; Linear constitutive equations contain two mechanical variables and two electrical variables [27]. The equations governing the direct and converse piezoelectric effect in the matrix form are expressed respectively as
-<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+{| class="formulaSCP" style="width: 100%; text-align: left;"
-<div style="text-align: right; direction: ltr; margin-left: 1em;">
-{| class="formulaSCP" style="width: 100%; text-align: center;"
 |-
-| <math>\left\{ D\right\} ={[e]}^{T}\, \left\{ S\right\} +\left[ {\epsilon }^{S}\right] \left\{ E\right\} \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \, \, \, (1)</math>
+|
+{| style="text-align: center; margin:auto;width: 100%;"
+|-
+| style="text-align: center;" | <math>\left\{ D\right\} ={[e]}^{T}\, \left\{ S\right\} +\left[ {\epsilon }^{S}\right] \left\{ E\right\}</math>
+|}
+| style="width: 5px;text-align: right;white-space: nowrap;" |(1)
 |}
-<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+{| class="formulaSCP" style="width: 100%; text-align: left;"
-<div style="text-align: right; direction: ltr; margin-left: 1em;">
-{| class="formulaSCP" style="width: 100%; text-align: center;"
 |-
-|<math display="inline">\left\{ T\right\} ={[c}^{E}]\, \left\{ S\right\} +\left[ e\right] \left\{ E\right\} \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \, \, \, (2)</math>
+|
+{| style="text-align: center; margin:auto;width: 100%;"
+|-
+| style="text-align: center;" | <math>\left\{ T\right\} ={[c}^{E}]\, \left\{ S\right\} +\left[ e\right] \left\{ E\right\}</math>
+|}
+| style="width: 5px;text-align: right;white-space: nowrap;" |(2)
 |}
-</div>
-</div>
-</div>
-</div>
-where <math display="inline">\left\{ D\right\}</math>  is the electric displacement, <math display="inline">{[e]}^{T}</math>  is the transpose of the dielectric permittivity matrix <math display="inline">\left[ e\right]</math> , <math display="inline">\left\{ S\right\}</math>  is the mechanical strain vector, <math display="inline">\left[ {\epsilon }^{S}\right]</math>  is the dielectric permittivity matrix at constant mechanical strain, <math display="inline">\left\{ E\right\} \,</math> is the electric intensity vector, <math display="inline">\left\{ T\right\}</math>  is the mechanical stress vector and <math display="inline">{[c}^{E}]</math> is the matrix of elastic stiffness at constant electric field strength. The same relationship can be written in three other forms, depending on which variables are chosen to be independent. The direct relationship given by Eq.1 is normally used when modeling the sensing capability of the piezoelectric material, while the actuator capability is modeled using the converse effect relationship given by Eq.2 [29].
+where <math display="inline">\left\{ D\right\}</math>  is the electric displacement, <math display="inline">{[e]}^{T}</math>  is the transpose of the dielectric permittivity matrix <math display="inline">\left[ e\right]</math>, <math display="inline">\left\{ S\right\}</math>  is the mechanical strain vector, <math display="inline">\left[ {\epsilon }^{S}\right]</math>  is the dielectric permittivity matrix at constant mechanical strain, <math display="inline">\left\{ E\right\} \,</math> is the electric intensity vector, <math display="inline">\left\{ T\right\}</math>  is the mechanical stress vector and <math display="inline">{[c}^{E}]</math> is the matrix of elastic stiffness at constant electric field strength. The same relationship can be written in three other forms, depending on which variables are chosen to be independent. The direct relationship given by Eq.1 is normally used when modeling the sensing capability of the piezoelectric material, while the actuator capability is modeled using the converse effect relationship given by Eq. (2) [29].
-==3. Simulation Method==
+==3. Simulation method==
-In this section, ANSYS software simulates active vibration control in the smart composite. Fig. 1 shows the block diagram of the obtained analysis. ANSYS/Multiphysics and ANSYS/Mechanical programs can be used for modeling structural fields with piezoelectric smart patches, whereas  <math display="inline">{F}_{e}\,</math> is the force of vibration generated. The instant value of vibration producing force <math display="inline">{F}_{e}</math> can be described at each time step. In the analysis below, it is taken as <math display="inline">{F}_{o}</math> at <math display="inline">t\, =</math><math>\, \Delta t</math> and zero for other time steps. At a piezoelectric sensor location, the strain ( ε )  is calculated. Reference input for vibration cancellation is zero. <math display="inline">\, {K}_{p}</math>, <math display="inline">{K}_{v}</math>, and <math display="inline">{K}_{s}</math> are respectively, the variables the control, power, and sensor amplification factors. <math display="inline">{K}_{v}\,</math> is taken as 200, and <math display="inline">{K}_{p}\,</math> is changed in the analysis below. Only the proportional control is applied. The proportional constant for the actuator voltage <math display="inline">{V}_{a}</math> was obtained from the multiplication of <math display="inline">{K}_{p}{\, K}_{v}\, {K}_{s}</math>. Therefore, the results are not affected by changing the values of <math display="inline">{K}_{s}</math>, <math display="inline">{K}_{p}</math> and <math display="inline">{K}_{v}</math> and keeping their multiplication the same. To evaluate the performance of the vibration control, the calculated deflection at a location, <math display="inline">dt</math>, is observed.
+In this section, ANSYS software simulates active vibration control in the smart composite. [[#img-1|Figure 1]] shows the block diagram of the obtained analysis. ANSYS/Multiphysics and ANSYS/Mechanical programs can be used for modeling structural fields with piezoelectric smart patches, whereas  <math display="inline">{F}_{e}\,</math> is the force of vibration generated. The instant value of vibration producing force <math display="inline">{F}_{e}</math> can be described at each time step. In the analysis below, it is taken as <math display="inline">{F}_{o}</math> at <math display="inline">t\, = \Delta t</math> and zero for other time steps. At a piezoelectric sensor location, the strain (<math>\varepsilon</math>)  is calculated. Reference input for vibration cancellation is zero. <math display="inline">\, {K}_{p}</math>, <math display="inline">{K}_{v}</math>, and <math display="inline">{K}_{s}</math> are respectively, the variables the control, power, and sensor amplification factors. <math display="inline">{K}_{v}\,</math> is taken as 200, and <math display="inline">{K}_{p}\,</math> is changed in the analysis below. Only the proportional control is applied. The proportional constant for the actuator voltage <math display="inline">{V}_{a}</math> was obtained from the multiplication of <math display="inline">{K}_{p}{\, K}_{v}\, {K}_{s}</math>. Therefore, the results are not affected by changing the values of <math display="inline">{K}_{s}</math>, <math display="inline">{K}_{p}</math> and <math display="inline">{K}_{v}</math> and keeping their multiplication the same. To evaluate the performance of the vibration control, the calculated deflection at a location, <math display="inline">dt</math>, is observed.
-{| style="width: 100%;border-collapse: collapse;"
+<div id='img-1'></div>
-|-
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;width: 100%;"|[[Image:Review_223450771672-image1.jpeg|432px]]
 |-
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Fig. 1 Block diagram of the analysis</span>
+|style="padding:10px;"| [[Image:Review_223450771672-image1.jpeg|432px]]
+|- style="text-align: center; font-size: 75%;"
+| colspan="1" style="padding-bottom:10px;"| '''Figure 1'''. Block diagram of the analysis
 |}
-<span style="text-align: left; font-size: 100%;">First, the structures of the beam type are considered. The configuration of the beam composite structure is shown in Fig. 2.</span>
+First, the structures of the beam type are considered. The configuration of the beam composite structure is shown in [[#img-2|Figure 2]].
-{| style="width: 100%;border-collapse: collapse;"
+<div id='img-1'></div>
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;width: 100%;"|[[Image:Review_223450771672-image2.png|444px]]
+|style="padding:10px;"| [[Image:Review_223450771672-image2.png|444px]]
-|-
+|- style="text-align: center; font-size: 75%;"
-| style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Fig. 2  Configuration for the structure of the beam type</span>
+| colspan="1" style="padding-bottom:10px;"| '''Figure 2'''. Configuration for the structure of the beam type
 |}
 ====3.1 Material Properties====
-<span style="text-align: center; font-size: 75%;">Finite Element Analysis is used as the numerical modeling tool to simulate the dynamic characteristics of the composite structures by ANSYS macro codes, and the macro file has been written by using APDL. The first step in the coding process defined model geometry by using direct mesh generation of smart structure. In first our system, we used a composite smart beam of size (400 mm x 28 mm ) and various thickness (1.59 mm, 2.12mm and 3.18mm) according to number of laminates, with (</span> <math display="inline">{d}_{33}MFC-</math><math>P1</math><span style="text-align: center; font-size: 75%;">) piezoelectric patch actuator type (M8528-P1). These typically consist of an active area (85mm x 28 mm) with (0.3 mm) thickness, and it is embedded on the beam from (10 mm) on the fixed point, as shown in Table 1 for all cases.</span>
+Finite Element Analysis is used as the numerical modeling tool to simulate the dynamic characteristics of the composite structures by ANSYS macro codes, and the macro file has been written by using APDL. The first step in the coding process defined model geometry by using direct mesh generation of smart structure. In first our system, we used a composite smart beam of size (400 mm <math>\times</math> 28 mm ) and various thickness (1.59 mm, 2.12mm and 3.18mm) according to number of laminates, with (<math display="inline">{d}_{33}MFC-</math><math>P1</math>) piezoelectric patch actuator type (M8528-P1). These typically consist of an active area (85mm <math>\times</math>  28 mm) with (0.3 mm) thickness, and it is embedded on the beam from (10 mm) on the fixed point, as shown in Table 1 for all cases.
 ====3.2 Simulation procedures ====
-In the first step, the beam structures type is taken and considered. Fig. 2 displays the configuration of the structure with sensor locations. The displacement value is taken from the displacement sensor, and similarly, the strain value is taken from the strain sensor as the feedback. Here, the product of ANSYS /Multiphysics (version 18) was applied to the structures of the smart composites model. The combination of the ANSYS modeling and solution into the control action codes were achieved by ANSYS [30].
+In the first step, the beam structures type is taken and considered. [[#img-2|Figure 2]] displays the configuration of the structure with sensor locations. The displacement value is taken from the displacement sensor, and similarly, the strain value is taken from the strain sensor as the feedback. Here, the product of ANSYS /Multiphysics (version 18) was applied to the structures of the smart composites model. The combination of the ANSYS modeling and solution into the control action codes were achieved by ANSYS [30].
-The composite smart structure beam is studied, as indicated in Table 1.
+The composite smart structure beam is studied, as indicated in [[#tab-1|Table 1]].
-{| style="width: 74%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+<div class="center" style="font-size: 75%;">'''Table 1'''. Dimensions and distances  of Composite beam structure</div>
+<div id='tab-1'></div>
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
+|-style="text-align:center"
+! Case No. !! Composite layers !! Dimensions of the structure (mm)'
+|-style="text-align:center"
+|  style="vertical-align: top;"|1
+|  style="text-align: center;vertical-align: top;"|(0/90/0)
+|  style="text-align: center;vertical-align: top;"|<math>400\times 28\times 1.59</math>
 |-
-|  colspan='3'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Table 1  Dimensions and distances  of Composite beam structure</span>
+|  style="text-align: center;vertical-align: top;"|2
+|  style="text-align: center;vertical-align: top;"|(0/90/0/90)
+|  style="text-align: center;vertical-align: top;"|<math>400\times 28\times 2.12</math>
 |-
-|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Case No.'''</span>
+|  style="text-align: center;vertical-align: top;"|3
-|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Composite layers'''</span>
+|  style="text-align: center;vertical-align: top;"|(0/90/90/0)
-|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Dimensions of the structure (mm)'''</span>
+|  style="text-align: center;vertical-align: top;"|<math>400\times 28\times 2.12</math>
 |-
-|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1</span>
+|  style="text-align: center;vertical-align: top;"|4
-|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/90/0)</span>
+|  style="text-align: center;vertical-align: top;"|(0/90/0/0/90/0)
-|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">400x28x1.59</span>
+|  style="text-align: center;vertical-align: top;"|<math>400\times 28\times 3.18</math>
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2</span>
+|  style="text-align: center;vertical-align: top;"|5
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/90/0/90)</span>
+|  style="text-align: center;vertical-align: top;"|(0/45/0)
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">400x28x2.12</span>
+|  style="text-align: center;vertical-align: top;"|<math>400\times 28\times  1.59</math>
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">3</span>
+|  style="text-align: center;vertical-align: top;"|6
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/90/90/0)</span>
+|  style="text-align: center;vertical-align: top;"|(0/45/0/45)
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">400x28x2.12</span>
+|  style="text-align: center;vertical-align: top;"|<math>400\times 28\times 2.12</math>
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">4</span>
+|  style="text-align: center;vertical-align: top;"|7
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/90/0/0/90/0)</span>
+|  style="text-align: center;vertical-align: top;"|(0/45/45/0)
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">400x28x3.18</span>
+|  style="text-align: center;vertical-align: top;"|<math>400\times 28\times 2.12</math>
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5</span>
+|  style="text-align: center;vertical-align: top;"|8
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/45/0)</span>
+|  style="text-align: center;vertical-align: top;"|(0/45/0/0/45/0)
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">400x28x1.59</span>
+|  style="text-align: center;vertical-align: top;"|<math>400\times 28\times 3.18</math>
-|-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">6</span>
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/45/0/45)</span>
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">400x28x2.12</span>
-|-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">7</span>
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/45/45/0)</span>
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">400x28x2.12</span>
-|-
-|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">8</span>
-|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/45/0/0/45/0)</span>
-|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">400x28x3.18</span>
-|-
-|  colspan='3'  style="border-top: 1pt solid black;vertical-align: top;"|
 |}
-<span style="text-align: center; font-size: 100%;">The macro file creates nodes and finite elements to the structure. In the first part, composite materials have used element type SOLID185 layered structural elements for composite beam or plate part, and element type (SOLID5) elements with (MASS21) were used in the second structure part for the piezoelectric part in the structure [31]. Overall, these macro codes as shown in Table 2 indicate the defined element types of various materials.</span>
+The macro file creates nodes and finite elements to the structure. In the first part, composite materials have used element type SOLID185 layered structural elements for composite beam or plate part, and element type (SOLID5) elements with (MASS21) were used in the second structure part for the piezoelectric part in the structure [31]. Overall, these macro codes as shown in [[#tab-2|Table 2]] indicate the defined element types of various materials.
+<div class="center" style="font-size: 75%;">'''Table 2'''. Element definition codes</div>
-{| style="width: 76%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+<div id='tab-2'></div>
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Table 2 Element definition codes</span>
+! style="text-align:left;" | /PREP7 !! style="text-align:left;" | ! start the preprocessing
 |-
-|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-bottom: 2pt solid white;border-right: 2pt solid white;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''/PREP7'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''! Define Element type #1'''
-|  style="border-top: 1pt solid black;border-left: 2pt solid white;border-bottom: 2pt solid white;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! start the preprocessing</span>
+|  style="vertical-align: top;"|
 |-
-|  style="border-top: 2pt solid white;border-left: 1pt solid black;border-bottom: 2pt solid white;border-right: 2pt solid white;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''! Define Element type #1'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''ET,1,185'''
-|  style="border-top: 2pt solid white;border-left: 2pt solid white;border-bottom: 2pt solid white;border-right: 1pt solid black;vertical-align: top;"|
+|  style="vertical-align: top;"|! chooses SOLID185 element for analysis composite material
 |-
-|  style="border-top: 2pt solid white;border-left: 1pt solid black;border-bottom: 2pt solid white;border-right: 2pt solid white;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''ET,1,185'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''KEYOPT,1,3,1'''
-|  style="border-top: 2pt solid white;border-left: 2pt solid white;border-bottom: 2pt solid white;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! chooses SOLID185 element for analysis composite material                       </span>
+|  style="vertical-align: top;"|! Key option #3 = 1  LAYERED SOLID  ELEMENTS
 |-
-|  style="border-top: 2pt solid white;border-left: 1pt solid black;border-bottom: 2pt solid white;border-right: 2pt solid white;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''KEYOPT,1,3,1'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''KEYOPT,1,2,2'''
-|  style="border-top: 2pt solid white;border-left: 2pt solid white;border-bottom: 2pt solid white;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Key option #3 = 1  LAYERED SOLID  ELEMENTS  </span>
+|  style="vertical-align: top;"|! Key option #2 = 2  ENHANCED STRAIN  FORMULATION
 |-
-|  style="border-top: 2pt solid white;border-left: 1pt solid black;border-bottom: 1pt solid white;border-right: 2pt solid white;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''KEYOPT,1,2,2'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''KEYOPT,1,8,1'''
-|  style="border-top: 2pt solid white;border-left: 2pt solid white;border-bottom: 1pt solid black;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Key option #2 = 2  ENHANCED STRAIN  FORMULATION</span>
+|  style="vertical-align: top;"|! Key option #8 = 1 STORE  DATA FOR ALL LAYERS
 |-
-|  style="border-top: 1pt solid white;border-left: 1pt solid black;border-bottom: 1pt solid white;border-right: 2pt solid white;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''KEYOPT,1,8,1'''</span>
+|  style="vertical-align: top;"|
-|  style="border-top: 2pt solid white;border-left: 2pt solid white;border-bottom: 1pt solid white;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Key option #8 = 1 STORE  DATA FOR ALL LAYERS</span>
+|  style="vertical-align: top;"|
 |-
-|  style="border-top: 1pt solid white;border-left: 1pt solid black;border-bottom: 1pt solid white;border-right: 2pt solid white;vertical-align: top;"|
+|   style="vertical-align: top;padding-left:10px;"|'''! Define Element type #2'''
-|  style="border-top: 1pt solid white;border-left: 2pt solid white;border-bottom: 1pt solid white;border-right: 1pt solid black;vertical-align: top;"|
+|  style="vertical-align: top;"|
 |-
-|   style="border-top: 1pt solid white;border-left: 1pt solid black;border-bottom: 1pt solid white;border-right: 2pt solid white;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''! Define Element type #2'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''ET,2,solid5,3'''
-|  style="border-top: 1pt solid white;border-left: 2pt solid white;border-bottom: 1pt solid white;border-right: 1pt solid black;vertical-align: top;"|
+|  style="vertical-align: top;"|! chooses SOLID5 element for analysis piezoelectric material
 |-
-|  style="border-top: 1pt solid white;border-left: 1pt solid black;border-bottom: 1pt solid white;border-right: 2pt solid white;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''ET,2,solid5,3'''</span>
+|  style="vertical-align: top;padding-left:10px;padding-bottom:5px;"|'''ET,3,mass21'''
-|  style="border-top: 2pt solid white;border-left: 2pt solid white;border-bottom: 1pt solid white;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! chooses SOLID5 element for analysis piezoelectric material</span>
+|  style="vertical-align: top;padding-bottom:5px;"|
-|-
-|  style="border-top: 1pt solid white;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 2pt solid white;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''ET,3,mass21'''</span>
-|  style="border-top: 1pt solid black;border-left: 2pt solid white;border-bottom: 1pt solid black;border-right: 1pt solid black;vertical-align: top;"|
 |}
-<span style="text-align: center; font-size: 75%;">Therefore, an essential next step to confirming in three-dimensional composite material properties is assigned.  As shown in Fig. 2, the structure is divided into two parts. The first part deals with composite material, and the other part is actuator and all properties, as shown in Table 3, and all that has been gained from [32][33]&[34]. </span>
+Therefore, an essential next step to confirming in three-dimensional composite material properties is assigned.  As shown in [[#img-2|Figure 2]], the structure is divided into two parts. The first part deals with composite material, and the other part is actuator and all properties, as shown in [[#tab-3|Table 3]], and all that has been gained from [32,33,34].
-{| style="width: 74%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+<div class="center" style="font-size: 75%;">'''Table 3'''. Material properties of MFC and CFRP</div>
+<div id='tab-3'></div>
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
+|-style="text-align:center"
+! colspan='2' | Material Properties !! colspan='2' |CFRP !! MFC
 |-
-|  colspan='3'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Table 3 Material properties of MFC and CFRP</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|Young’s modulus  <math display="inline">{\boldsymbol{E}}_{{11}}</math>
+|  style="text-align: center;vertical-align: top;"|140 GPa
+|  colspan='2'  style="text-align: center;vertical-align: top;"|30.34 GPa
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Material Properties'''</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|Young’s modulus  <math display="inline">{\boldsymbol{E}}_{{22}}</math>
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''CFRP'''</span>
+|  style="text-align: center;vertical-align: top;"|l8.75 GPa
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MFC'''</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|15.86 GPa
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Young’s modulus </span> <math display="inline">{\boldsymbol{E}}_{\boldsymbol{11}}</math>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|Young’s modulus  <math display="inline">{\boldsymbol{E}}_{{33}}</math>
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\boldsymbol{140\, GPa}</math>
+|  style="text-align: center;vertical-align: top;"|8.75 GPa
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">30.34 GPa</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|9.03 GPa
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Young’s modulus </span> <math display="inline">{\boldsymbol{E}}_{\boldsymbol{22}}</math>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|Shear modulus  <math display="inline">{\boldsymbol{G}}_{{12}}</math>
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\boldsymbol{8.75\, GPa}</math>
+|  style="text-align: center;vertical-align: top;"|7.2 GPa
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">15.86 GPa</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|5.52 GPa
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Young’s modulus </span> <math display="inline">{\boldsymbol{E}}_{\boldsymbol{33}}</math>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|Poisson’s ratio  <math display="inline">{\boldsymbol{v}}_{{12}}</math>
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\boldsymbol{8.75\, GPa}</math>
+|  style="text-align: center;vertical-align: top;"|0.3
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">9.03 GPa</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|0.31
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Shear modulus </span> <math display="inline">{\boldsymbol{G}}_{\boldsymbol{12}}</math>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|Poisson’s ratio  <math display="inline">{\boldsymbol{v}}_{{13}}</math>
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\boldsymbol{7.2\, GPa}</math>
+|  style="text-align: center;vertical-align: top;"|0.3
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5.52 GPa</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|0.30
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Poisson’s ratio </span> <math display="inline">{\boldsymbol{v}}_{\boldsymbol{12}}</math>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|Density <math>\rho</math>
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\boldsymbol{0.3}</math>
+|  style="text-align: center;vertical-align: top;"|1550 kg/m<math>^3</math>
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.31</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|5020  kg/m<math>^3</math>
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Poisson’s ratio </span> <math display="inline">{\boldsymbol{v}}_{\boldsymbol{13}}</math>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|Piezoelectric Stress constant (<math display="inline">{\boldsymbol{e}}_{{32}}</math>)
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\boldsymbol{0.3}</math>
+|  style="text-align: center;vertical-align: top;"|
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.30</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|-3.47  C/m<math>^2</math>
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Density ρ</span>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|Piezoelectric Stress constant (<math display="inline">{\boldsymbol{e}}_{{33}}</math>)
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\boldsymbol{1550\, kg/}{\boldsymbol{m}}^{\boldsymbol{3}}</math>
+|  style="text-align: center;vertical-align: top;"|
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5020 </span> <math display="inline">\boldsymbol{\, kg/}{\boldsymbol{m}}^{\boldsymbol{3}}</math>
+|  colspan='2'  style="text-align: center;vertical-align: top;"|12.36   C/m<math>^2</math>
-|-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Piezoelectric Stress constant (</span> <math display="inline">{\boldsymbol{e}}_{\boldsymbol{32}}</math><span style="text-align: center; font-size: 75%;">)</span>
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-3.47 </span> <math display="inline">\boldsymbol{C/}{\boldsymbol{m}}^{\boldsymbol{2}}</math>
-|-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Piezoelectric Stress constant (</span> <math display="inline">{\boldsymbol{e}}_{\boldsymbol{33}}</math><span style="text-align: center; font-size: 75%;">)</span>
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">12.36 </span> <math display="inline">\boldsymbol{C/}{\boldsymbol{m}}^{\boldsymbol{2}}</math>
 |}
-<span style="text-align: center; font-size: 100%;">The ANSYS macro of the material properties is given in Table 4. The first material is represented by composite material (Carbon fiber reinforced polymer), and the second material is represented by the piezoelectric (MFC) actuator material.</span>
+The ANSYS macro of the material properties is given in [[#tab-4|Table 4]]. The first material is represented by composite material (Carbon fiber reinforced polymer), and the second material is represented by the piezoelectric (MFC) actuator material.
-{| style="width: 74%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+<div class="center" style="font-size: 75%;">'''Table 4'''. Composite material properties codes</div>
-|-
-|  colspan='1'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Table 4 Composite material properties codes</span>
+<div id='tab-4'></div>
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
+|-style="text-align:left;"
+!  colspan='2' style="text-align:left;" | /PREP7 !!  colspan='2' style="text-align:left;" |! Define  Material 1# : <math display="inline">E1= 140\times {10}^{9}\, Pa</math>
 |-
-|  colspan='2'  style="border-top: 1pt solid black;border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''/PREP7'''</span>
+|  colspan='2'  style="vertical-align: top;padding-left:10px;"|'''MP,EX,1,140E9'''
+|colspan='2'  style="vertical-align: top;"|
-<span style="text-align: center; font-size: 75%;">'''MP,EX,1,140E9'''</span>
-|  colspan='2'  style="border-top: 1pt solid black;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# : E1=</span> <math display="inline">140x{10}^{9}\, Pa</math>
 |-
-|  colspan='2'  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,EY,1,8.75E9'''</span>
+|  colspan='2'  style="vertical-align: top;width:80px;padding-left:10px;"|'''MP,EY,1,8.75E9'''
-|  colspan='2'  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# : E2=</span> <math display="inline">8.75x{10}^{9}\, Pa</math>
+|  colspan='2'  style="vertical-align: top;"|! Define  Material 1# : <math display="inline">E2=8.75\times {10}^{9}\, Pa</math>
 |-
-|  colspan='2'  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,EZ,1,8.75E9'''</span>
+|  colspan='2'  style="vertical-align: top;padding-left:10px;"|'''MP,EZ,1,8.75E9'''
-|  colspan='2'  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# : EZ=EY ; E3=</span> <math display="inline">8.75x{10}^{9}\, Pa</math>
+|  colspan='2'  style="vertical-align: top;"|! Define  Material 1# : <math display="inline">EZ=EY \,;\, E3= 8.75\times {10}^{9}\, Pa</math>
 |-
-|  colspan='2'  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,GXY,1,7.2E9'''</span>
+|  colspan='2'  style="vertical-align: top;padding-left:10px;"|'''MP,GXY,1,7.2E9'''
-|  colspan='2'  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# : shear modulus xy-plane  </span> <math display="inline">{G}_{xy}=</math><math>7.2x{10}^{9}\, Pa</math>
+|  colspan='2'  style="vertical-align: top;"|! Define  Material 1# : shear modulus xy-plane   <math display="inline">{G}_{xy}=</math><math>7.2\times {10}^{9}\, Pa</math>
 |-
-|  colspan='2'  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,GYZ,1,6.21E9'''</span>
+|  colspan='2'  style="vertical-align: top;padding-left:10px;"|'''MP,GYZ,1,6.21E9'''
-|  colspan='2'  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# : shear modulus xz-plane   </span> <math display="inline">{G}_{xy}=</math><math>6.2x{10}^{9}\, Pa</math>
+|  colspan='2'  style="vertical-align: top;"|! Define  Material 1# : shear modulus xz-plane   <math display="inline">{G}_{xy}=6.2\times {10}^{9}\, Pa</math>
 |-
-|  colspan='2'  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,GXZ,1,7.2E9'''</span>
+|  colspan='2'  style="vertical-align: top;padding-left:10px;"|'''MP,GXZ,1,7.2E9'''
-|  colspan='2'  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# : shear modulus yz-plane   </span> <math display="inline">{G}_{xy}=</math><math>7.2x{10}^{9}\, Pa</math>
+|  colspan='2'  style="vertical-align: top;"|! Define  Material 1# : shear modulus yz-plane   <math display="inline">{G}_{xy}=7.2\times {10}^{9}\, Pa</math>
 |-
-|  colspan='2'  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,PRXY,1,0.3'''</span>
+|  colspan='2'  style="vertical-align: top;padding-left:10px;"|'''MP,PRXY,1,0.3'''
-|  colspan='2'  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# : Poisson coefficient </span> <math display="inline">{\upsilon }_{1}=</math><math>0.3</math>
+|  colspan='2'  style="vertical-align: top;"|! Define  Material 1# : Poisson coefficient <math display="inline">{\upsilon }_{1}=0.3</math>
 |-
-|  colspan='2'  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,PRYZ,1,0.3'''</span>
+|  colspan='2'  style="vertical-align: top;padding-left:10px;"|'''MP,PRYZ,1,0.3'''
-|  colspan='2'  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# : Poisson coefficient </span> <math display="inline">{\upsilon }_{2}=</math><math>0.3</math>
+|  colspan='2'  style="vertical-align: top;"|! Define  Material 1# : Poisson coefficient  <math display="inline">{\upsilon }_{2}=0.3</math>
 |-
-|  colspan='2'  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,PRXZ,1,0.3'''</span>
+|  colspan='2'  style="vertical-align: top;padding-left:10px;"|'''MP,PRXZ,1,0.3'''
-|  colspan='2'  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# : Poisson coefficient </span> <math display="inline">{\upsilon }_{3}=</math><math>0.3</math>
+|  colspan='2'  style="vertical-align: top;"|! Define  Material 1# : Poisson coefficient  <math display="inline">{\upsilon }_{3}=0.3</math>
 |-
-|  colspan='2'  style="border-left: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,DENS,1,1550'''</span>
+|  colspan='2'  style="vertical-align: top;padding-left:10px;padding-bottom:5px;"|'''MP,DENS,1,1550'''
-|  colspan='2'  style="border-bottom: 1pt solid black;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define  Material 1# :  Density for composite material</span>
+|  colspan='2'  style="vertical-align: top;"|! Define  Material 1# :  Density for composite material
 |}
-<span style="text-align: center; font-size: 100%;">The number of composite layers, the thickness of layers, and orientation angles of fibers are assigned as macro APDL ANSYS codes for the first material, as given below in Table 5.</span>
+The number of composite layers, the thickness of layers, and orientation angles of fibers are assigned as macro APDL ANSYS codes for the first material, as given below in [[#tab-5|Table 5]].
-{| style="width: 76%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+<div class="center" style="font-size: 75%;">'''Table 5'''. Orientation angles’ Composite layers codes</div>
+<div id='tab-5'></div>
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
+|-style="text-align:left;"
+! style="text-align:left;"|SECTYPE,1,SHELL !! style="text-align:left;"| set #1, [0/90/0]
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Table 5 Orientation angles’ Composite layers codes</span>
+|  style="vertical-align: top;padding-left:10px;"|'''SECDATA,0.53e-3,1,0'''
+|  style="vertical-align: top;"|! 1st  LAYER 1: Thk=0.53e-3m, THETA = 0°
 |-
-|  style="border-top: 1pt solid black;border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''SECTYPE,1,SHELL'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''SECDATA,0.53e-3,1,90'''
-|  style="border-top: 1pt solid black;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">set #1, [0/90/0]</span>
+|  style="vertical-align: top;"|! 2nd LAYER 2: Thk=0.53e-3m, THETA = 90°
 |-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''SECDATA,0.53e-3,1,0'''</span>
+|  style="vertical-align: top;padding-left:10px;padding-bottom:5px;"|'''SECDATA,0.53e-3,1,0'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! 1st  LAYER 1: Thk=0.53e-3m, THETA = 0</span> <math display="inline">\boldsymbol{\mbox{°}}</math>
+|  style="vertical-align: top;padding-bottom:5px;"|! 3rd LAYER 4: Thk=0.53e-3m, THETA = 0°
-|-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''SECDATA,0.53e-3,1,90'''</span>
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! 2nd LAYER 2: Thk=0.53e-3m, THETA = 90</span> <math display="inline">\boldsymbol{\mbox{°}}</math>
-|-
-|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''SECDATA,0.53e-3,1,0'''</span>
-|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! 3rd LAYER 4: Thk=0.53e-3m, THETA = 0</span> <math display="inline">\boldsymbol{\mbox{°}}</math>
 |}
-<span style="text-align: center; font-size: 100%;">The second material piezoelectric actuator properties are given in Table 6 as macro ANSYS.</span>
+The second material piezoelectric actuator properties are given in [[#tab-6|Table 6]] as macro ANSYS.
+<div class="center" style="font-size: 75%;">'''Table 6'''. MFC properties codes</div>
+<div id='tab-6'></div>
-{| style="width: 76%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
+|-style="text-align:left;"
+! ! Define dielectric constant !!
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Table 6  MFC properties codes</span>
+|  style="vertical-align: top;padding-left:10px;"|'''MP,DENS,3,5020'''
+|  style="vertical-align: top;"|! Density for MFC  piezoelectric material
 |-
-|  style="border-top: 1pt solid black;border-left: 1pt solid black;  vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''! Define dielectric constant'''</span>
+|  style="vertical-align: top;padding-left:10px;"| '''MP,PERX,3,15.586E-9'''
-|  style="border-top: 1pt solid black; border-right: 1pt solid black;vertical-align: top;"|
+|  style=" vertical-align: top;"|! Permittivity in X direction
 |-
-|  style=" ;border-left: 1pt solid black;  vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,DENS,3,5020'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''MP,PERY,3,15.586E-9'''
-|  style=" border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Density for MFC  piezoelectric material</span>
+|  style="vertical-align: top;"|! Permittivity in Y direction
 |-
-|  style=" ;border-left: 1pt solid black;  vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,PERX,3,15.586E-9'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''MP,PERZ,3,11.79E-9'''
-|  style=" border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Permittivity in X direction</span>
+|  style="vertical-align: top;"|! Permittivity in Z direction
 |-
-|  style="border-left: 1pt solid black; vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,PERY,3,15.586E-9'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''TB,PIEZ,3'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Permittivity in Y direction</span>
+|  style="vertical-align: top;"|! Define piezoelectric table
 |-
-|  style="border-left: 1pt solid black; vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''MP,PERZ,3,11.79E-9'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''TBDATA,16,0'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Permittivity in Z direction</span>
+|  style="vertical-align: top;"|! e16 piezoelectric constant
 |-
-|  style="border-left: 1pt solid black; vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''TB,PIEZ,3'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''TBDATA,14,0'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define piezoelectric table</span>
+|  style="vertical-align: top;"|! e25 piezoelectric constant
 |-
-|  style="border-left: 1pt solid black; vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''TBDATA,16,0'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''TBDATA,6,-3.47'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! e16 piezoelectric constant</span>
+|  style="vertical-align: top;"|! e32 piezoelectric constant
 |-
-|  style="border-left: 1pt solid black; vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''TBDATA,14,0'''</span>
+|  style="vertical-align: top;padding-left:10px;"|'''TBDATA,9,12.36'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! e25 piezoelectric constant</span>
+|  style="vertical-align: top;"|! e33 piezoelectric constant
 |-
-|  style="border-left: 1pt solid black; vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''TBDATA,6,-3.47'''</span>
+|  style="vertical-align: top;padding-left:10px;padding-bottom:5px;"|'''TB,ANEL,3'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! e32 piezoelectric constant</span>
+|  style="vertical-align: top;"|! Define structural table
-|-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''TBDATA,9,12.36'''</span>
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! e33 piezoelectric constant</span>
-|-
-|  style="border-left: 1pt solid black;border-bottom: 1pt solid black; vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''TB,ANEL,3'''</span>
-|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! Define structural table</span>
 |}
-<span style="text-align: center; font-size: 100%;">    After all dimensions and properties of both composite and piezoelectric properties are defined, the type of element is assigned for each type of material. SOLID185 layered structural elements are assigned for the composite part, and SOLID5 elements are used for the piezoelectric part of the structure. Then, the static solution is analyzed to evaluate the initial deflections of each node and element. The model is discretized into finite elements model, as shown in Fig. 3.
+After all dimensions and properties of both composite and piezoelectric properties are defined, the type of element is assigned for each type of material. SOLID185 layered structural elements are assigned for the composite part, and SOLID5 elements are used for the piezoelectric part of the structure. Then, the static solution is analyzed to evaluate the initial deflections of each node and element. The model is discretized into finite elements model, as shown in [[#img-3|Figure 3]].
+<div id='img-3'></div>
-The first three mode shapes and natural frequencies of the composite plates were obtained using APDL ANSYS macro codes. The obtained results were shown in Fig. 4. The first three natural frequencies evaluated by the modal shape analysis method using ANSYS (15.665 Hz), (94.17Hz), and (161.63Hz), respectively, for the first case considered in the present study, are shown in Fig. 4.</span>
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
-<span style="text-align: center; font-size: 100%;">     Fixed boundary conditions are described for the nodes </span> <math display="inline">\, (x=</math><math>0)</math><span style="text-align: center; font-size: 100%;">and given in Table 7 as macro ANSYS.</span>
-{| style="width: 76%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
-|-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Table 7 Boundary conditions of composite beam codes</span>
 |-
-|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''CP,1,volt,all    '''</span>
+|style="padding:10px;"|  [[Image:Review_223450771672-image3.jpeg|474px]]
-|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">! set #1 of coupled degrees of freedom.</span>
+|- style="text-align: center; font-size: 75%;"
+| colspan="1" style="padding-bottom:10px;"| '''Figure 3'''. Finite element model
 |}
-<span style="text-align: center; font-size: 100%;"> The degrees of freedom of the top and bottom surfaces for piezo actuator and sensor are coupled simultaneously with volt by APDL ANSYS “CP” command. Modal analysis was performed to determine parameters from vibration data measures like the inherent dynamic characteristics of the system in forms of natural frequencies, mode shapes and damping factor, and so to determine time step [35].</span>
+The first three mode shapes and natural frequencies of the composite plates were obtained using APDL ANSYS macro codes. The obtained results were shown in [[#img-4|Figure 4]]. The first three natural frequencies evaluated by the modal shape analysis method using ANSYS (15.665 Hz), (94.17Hz), and (161.63Hz), respectively, for the first case considered in the present study, are shown in [[#img-4|Figure 4]].
+<div id='img-4'></div>
-{| style="width: 76%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|  style="border: 1pt solid black;text-align: center;vertical-align: top;width: 100%;"|[[Image:Review_223450771672-image3.jpeg|474px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image4.jpeg|298px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image5.jpeg|298px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image6.jpeg|298px]]
 |-
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Fig. 3  Finite element model</span>
+|  style="text-align: center;font-size: 75%;"|1st mode shape
+|  style="text-align: center;font-size: 75%;"|2nd mode shape
+|  style="text-align: center;font-size: 75%;"|3rd mode shape
+|- style="text-align: center; font-size: 75%;"
+| colspan="3" style="padding:10px;"| '''Figure 4'''. Mode shapes of the first case
 |}
-<span style="text-align: center; font-size: 100%;">  Householder method (reduced method) can be used for the beam structures which have the solid coupled field in ANSYS. The time step is the incremental change in time and the time step is chosen by equation formula  </span> <math display="inline">\Delta t\, =</math><math>\, 1/(20\, {f}_{h})</math><span style="text-align: center; font-size: 100%;">, where </span> <math display="inline">{f}_{h}\,</math> <span style="text-align: center; font-size: 100%;"> is the highest natural frequency to be considered. Table 8 shows the three natural frequencies for the undamped system. The first mode is considered to calculate the time step and </span> <math display="inline">\Delta t</math><span style="text-align: center; font-size: 100%;"> is 0.0032, 0.0034, 0.0025, 0.00183, 0.00318, 0.0032, 0.0025, and 0.0018, for cases 1, 2, 3, 4, 5, 6, 7, and 8, respectively.</span>
+Fixed boundary conditions are described for the nodes  <math display="inline">\, (x=0)</math> and given in [[#tab-7|Table 7]] as macro ANSYS.
-{| style="width: 100%;border-collapse: collapse;"
+<div class="center" style="font-size: 75%;">'''Table 7'''. Boundary conditions of composite beam codes</div>
-|-
-|  style="border: 1pt solid black;text-align: center;width: 34%;"|[[Image:Review_223450771672-image4.jpeg|198px]]
+<div id='tab-7'></div>
-|  style="border: 1pt solid black;text-align: center;width: 33%;"|[[Image:Review_223450771672-image5.jpeg|198px]]
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-|  style="border: 1pt solid black;text-align: center;width: 32%;"|[[Image:Review_223450771672-image6.jpeg|198px]]
+|-style="text-align:center"
-|-
+! CP,1,volt,all    !! ! set #1 of coupled degrees of freedom.
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1st mode shape</span>
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2nd mode shape</span>
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3rd mode shape</span>
-|-
-|  colspan='3'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Fig. 4  Mode shapes of the first case</span>
 |}
-<span style="text-align: center; font-size: 100%;">In this study, The strain feedback (piezoelectric sensor) is calculated at the selected sensor location, multiplied by </span> <math display="inline">{Ks}_{1}</math><span style="text-align: center; font-size: 100%;"> to convert the strain feedback to voltage and then subtracted from zero to calculate the error signal.</span>
+The degrees of freedom of the top and bottom surfaces for piezo actuator and sensor are coupled simultaneously with volt by APDL ANSYS “CP” command. Modal analysis was performed to determine parameters from vibration data measures like the inherent dynamic characteristics of the system in forms of natural frequencies, mode shapes and damping factor, and so to determine time step [35].
-<span style="text-align: center; font-size: 100%;">    The displacement feedback is also calculated at the selected sensor location, multiplied by </span> <math display="inline">{Ks}_{2}</math><span style="text-align: center; font-size: 100%;"> in order to convert the displacement feedback to voltage, and also subtracted from zero to evaluate the error signal between the input and output.</span>
+Householder method (reduced method) can be used for the beam structures which have the solid coupled field in ANSYS. The time step is the incremental change in time and the time step is chosen by equation formula  <math display="inline">\Delta t\, = 1/(20\, {f}_{h})</math>, where <math display="inline">{f}_{h}\,</math> is the highest natural frequency to be considered. [[#tab-8|Table 8]] shows the three natural frequencies for the undamped system. The first mode is considered to calculate the time step and <math display="inline">\Delta t</math> is 0.0032, 0.0034, 0.0025, 0.00183, 0.00318, 0.0032, 0.0025, and 0.0018, for cases 1, 2, 3, 4, 5, 6, 7, and 8, respectively.
-{| style="width: 76%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+In this study, The strain feedback (piezoelectric sensor) is calculated at the selected sensor location, multiplied by  <math display="inline">{Ks}_{1}</math> to convert the strain feedback to voltage and then subtracted from zero to calculate the error signal.
-|-
-|  colspan='5'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Table 8'''  '''Natural frequencies of the system</span>
-|-
-|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Case No.'''</span>
-|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Composite layers'''</span>
-|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''First '''</span>
-<span style="text-align: center; font-size: 75%;">'''(Hz)'''</span>
+The displacement feedback is also calculated at the selected sensor location, multiplied by <math display="inline">{Ks}_{2}</math> in order to convert the displacement feedback to voltage, and also subtracted from zero to evaluate the error signal between the input and output.
-|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''Second'''</span>
-<span style="text-align: center; font-size: 75%;">''' (Hz)'''</span>
+<div class="center" style="font-size: 75%;">'''Table 8'''. Natural frequencies of the system</div>
-|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">'''Third'''</span>
-<span style="text-align: center; font-size: 75%;">''' (Hz)'''</span>
+<div id='tab-8'></div>
-|-
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-|  style="border-top: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1</span>
+|-style="text-align:center"
-|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/90/0)</span>
+! Case No. !! Composite layers !! First (Hz) !! Second (Hz) !! Third (Hz)
-|  style="border-top: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">15.665</span>
-|  style="border-top: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">94.17</span>
-|  style="border-top: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">161.63</span>
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2</span>
+|  style="text-align: center;vertical-align: top;"|1
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/90/0/90)</span>
+|  style="text-align: center;vertical-align: top;"|(0/90/0)
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">14.888</span>
+|  style="text-align: center;vertical-align: bottom;"|15.665
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">87.674</span>
+|  style="text-align: center;vertical-align: bottom;"|94.17
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">193.47</span>
+|  style="text-align: center;vertical-align: bottom;"|161.63
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">3</span>
+|  style="text-align: center;vertical-align: top;"|2
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/90/90/0)</span>
+|  style="text-align: center;vertical-align: top;"|(0/90/0/90)
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">19.795</span>
+|  style="text-align: center;vertical-align: bottom;"|14.888
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">120.39</span>
+|  style="text-align: center;vertical-align: bottom;"|87.674
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">196.81</span>
+|  style="text-align: center;vertical-align: bottom;"|193.47
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">4</span>
+|  style="text-align: center;vertical-align: top;"|3
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/90/0/0/90/0)</span>
+|  style="text-align: center;vertical-align: top;"|(0/90/90/0)
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">27.265</span>
+|  style="text-align: center;vertical-align: bottom;"|19.795
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">165.45</span>
+|  style="text-align: center;vertical-align: bottom;"|120.39
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">217.71</span>
+|  style="text-align: center;vertical-align: bottom;"|196.81
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5</span>
+|  style="text-align: center;vertical-align: top;"|4
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/45/0)</span>
+|  style="text-align: center;vertical-align: top;"|(0/90/0/0/90/0)
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">15.701</span>
+|  style="text-align: center;vertical-align: bottom;"|27.265
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">94.379</span>
+|  style="text-align: center;vertical-align: bottom;"|165.45
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">170.62</span>
+|  style="text-align: center;vertical-align: bottom;"|217.71
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">6</span>
+|  style="text-align: center;vertical-align: top;"|5
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/45/0/45)</span>
+|  style="text-align: center;vertical-align: top;"|(0/45/0)
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">15.589</span>
+|  style="text-align: center;vertical-align: bottom;"|15.701
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">92.137</span>
+|  style="text-align: center;vertical-align: bottom;"|94.379
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">200.82</span>
+|  style="text-align: center;vertical-align: bottom;"|170.62
 |-
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">7</span>
+|  style="text-align: center;vertical-align: top;"|6
-|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/45/45/0)</span>
+|  style="text-align: center;vertical-align: top;"|(0/45/0/45)
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">19.945</span>
+|  style="text-align: center;vertical-align: bottom;"|15.589
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">121.29</span>
+|  style="text-align: center;vertical-align: bottom;"|92.137
-|  style="text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">204.18</span>
+|  style="text-align: center;vertical-align: bottom;"|200.82
 |-
-|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">8</span>
+|  style="text-align: center;vertical-align: top;"|7
-|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(0/45/0/0/45/0)</span>
+|  style="text-align: center;vertical-align: top;"|(0/45/45/0)
-|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">27.697</span>
+|  style="text-align: center;vertical-align: bottom;"|19.945
-|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">168.02</span>
+|  style="text-align: center;vertical-align: bottom;"|121.29
-|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|<span style="text-align: center; font-size: 75%;">222.82</span>
+|  style="text-align: center;vertical-align: bottom;"|204.18
 |-
-|  colspan='5'  style="border-top: 1pt solid black;vertical-align: top;"|
+|  style="text-align: center;vertical-align: top;"|8
+|  style="text-align: center;vertical-align: top;"|(0/45/0/0/45/0)
+|  style="text-align: center;vertical-align: bottom;"|27.697
+|  style="text-align: center;vertical-align: bottom;"|168.02
+|  style="text-align: center;vertical-align: bottom;"|222.82
 |}
-<span style="text-align: center; font-size: 100%;">   The difference between the input reference and the sensor signal is defined as the error signal. The error value is multiplied by </span> <math display="inline">{Kp}_{1}</math><span style="text-align: center; font-size: 100%;">for strain feedback and </span> <math display="inline">Kp\,</math> <span style="text-align: center; font-size: 100%;">for displacement feedback, and then both results are multiplied by the amplifier factor </span> <math display="inline">Kv</math><span style="text-align: center; font-size: 100%;"> to determine </span> <math display="inline">Va</math><span style="text-align: center; font-size: 100%;"> at a given time step. </span>
+The difference between the input reference and the sensor signal is defined as the error signal. The error value is multiplied by <math display="inline">{Kp}_{1}</math> for strain feedback and <math display="inline">Kp\,</math> for displacement feedback, and then both results are multiplied by the amplifier factor <math display="inline">Kv</math> to determine  <math display="inline">Va</math> at a given time step.
-The macro part that enables the calculations for the analysis of the closed loop for t > ∆t is given in table 9.
+The macro part that enables the calculations for the analysis of the closed loop for <math display="inline">t > \Delta t</math> is given in [[#tab-9|Table 9]].
-{| style="width: 76%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+<div class="center" style="font-size: 75%;">'''Table 9'''. Sensor feedback codes</div>
+<div id='tab-9'></div>
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
+|-style="text-align:left;"
+! style="text-align:left;"| *do,tt,2*dt,ts,dt !! style="text-align:left;"|!time loops
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Table 9  Sensor feedback codes</span>
+|  style="vertical-align: top;padding-left:5px;"|'''*get,uz-tip,node,nldm,u,z'''
+|  style="vertical-align: top;"|!sensor feedback
 |-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''*do,tt,2*dt,ts,dt'''</span>
+|  style="vertical-align: top;padding-left:5px;"|'''err=ref-ks*uz-tip'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">!time loops</span>
+|  style="vertical-align: top;"|!error signal
 |-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''*get,uz-tip,node,nldm,u,z'''</span>
+|  style="vertical-align: top;padding-left:5px;"|'''Vact=kp*kv*err'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">!sensor feedback</span>
+|  style="vertical-align: top;"|!actuator voltage
 |-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''err=ref-ks*uz-tip'''</span>
+|  style="vertical-align: top;padding-left:5px;"|'''d,nv,volt,Vact'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">!error signal</span>
+|  style="vertical-align: top;"|!Appling actuator volt
 |-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Vact=kp*kv*err'''</span>
+|  style="vertical-align: top;padding-left:5px;"|'''time,tt'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">!actuator voltage</span>
+|  style="vertical-align: top;"|
 |-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''d,nv,volt,Vact'''</span>
+|  style="vertical-align: top;padding-left:5px;"|'''solve'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">!Appling actuator volt</span>
+|  style="vertical-align: top;"|
 |-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''time,tt'''</span>
+|  style="vertical-align: top;padding-bottom:5px;padding-left:5px;"|'''*enddo'''
-|  style="border-right: 1pt solid black;vertical-align: top;"|
+|  style="vertical-align: top;"|
-|-
-|  style="border-left: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''solve'''</span>
-|  style="border-right: 1pt solid black;vertical-align: top;"|
-|-
-|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''*enddo'''</span>
-|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;vertical-align: top;"|
 |}
+==4. Results and discussion==
-==4. Results and Discussion==
 There are two methods to deal with active vibration control. The first method is displacement feedback signals by using the displacement sensor, and the other method is strain feedback signals by using a piezoelectric strain sensor. Displacement feedback signals have been used for the displacement to suppression of vibration.
-Classical controller (P) is the proportional control and gain (kp) is the only variable to be changed. The present study has been used for (2,4,6, and 8)   kp value. The open-loop and closed-loop responses of the tip displacement are shown in Fig.s 5(a), 5(b),5(c), and 5(d), respectively, by using APDL ANSYS. From the graph, it is found that the open-loop response vibrates a lot at a time of more than 5 sec in first Fig. 5(a), 4 sec in second Fig. 5(b), 3.8 in third Fig. 5(c), and 1.25 sec in the last Fig. 5(d) [36]. In the first case, composite laminates consist of three-ply (0/90/0) closed-loop responses by using classical controller when we took into consideration only proportional part effects in PID control, which is shown settle at a time 1.4 sec, 0.85 sec, 0.75 sec and 0.7 sec, respectively, proportional gain (kp) is 2, 3, 6, 8, which demonstrates that three recent responses of the system will be near each other. This reduction in the time compared with the first case is due to the increase in the number of composite layers, and, at the same time, an increase in the strength properties of the beam is observed.
+Classical controller (<math>P</math>) is the proportional control and gain (<math>kp</math>) is the only variable to be changed. The present study has been used for (2,4,6, and 8)   <math>kp</math> value. The open-loop and closed-loop responses of the tip displacement are shown in [[#img-5|Figures 5]](a), (b), (c) and (d), respectively, by using APDL ANSYS. From the graph, it is found that the open-loop response vibrates a lot at a time of more than 5 sec in first [[#img-5|Figure 5]](a), 4 sec in second [[#img-5|Figure 5]](b), 3.8 in third [[#img-5|Figure 5]](c), and 1.25 sec in the last [[#img-5|Figure 5]](d) [36]. In the first case, composite laminates consist of three-ply (<math>0/90/0</math>) closed-loop responses by using classical controller when we took into consideration only proportional part effects in PID control, which is shown settle at a time 1.4 sec, 0.85 sec, 0.75 sec and 0.7 sec, respectively, proportional gain (<math>kp</math>) is 2, 3, 6, 8, which demonstrates that three recent responses of the system will be near each other. This reduction in the time compared with the first case is due to the increase in the number of composite layers, and, at the same time, an increase in the strength properties of the beam is observed.
+<div id='img-5'></div>
-{| style="width: 100%;border-collapse: collapse;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;"> [[Image:Review_223450771672-image7.jpeg|294px]] </span>
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image7.jpeg|394px]]
-|  colspan='2'  style="border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;"> [[Image:Review_223450771672-image8.jpeg|294px]] </span>
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image8.jpeg|394px]]
 |-
-|  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(a) CFRP (0/90/0)</span>
+|  style="text-align: center;font-size: 75%;"|(a) CFRP (0/90/0)
-|  colspan='2'  style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(b) CFRP (0/90/0/90)</span>
+|  style="text-align: center;font-size: 75%;"|(b) CFRP (0/90/0/90)
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;"> [[Image:Review_223450771672-image9.jpeg|294px]] </span>
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image9.jpeg|394px]]
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;"> [[Image:Review_223450771672-image10.jpeg|294px]] </span>
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image10.jpeg|394px]]
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(c) CFRP (0/90/90/0)</span>
+|  style="text-align: center;font-size: 75%;"|(c) CFRP (0/90/90/0)
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(d) CFRP (0/90/0/90/0/90)</span>
+|  style="text-align: center;font-size: 75%;"|(d) CFRP (0/90/0/90/0/90
-|-
+|- style="text-align: center; font-size: 75%;"
-|  colspan='3'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Fig. 5 Open-loop and closed-loop responses of a composite beam of tip displacement.</span>
+| colspan="2" style="padding:10px;"| '''Figure 5'''. Open-loop and closed-loop responses of a composite beam of tip displacement
 |}
-<span style="text-align: center; font-size: 100%;">However, the evidence for this relationship is conclusive and has been retrieved from PID rules; with increasing the proportional range, the vibration converges slowly and settling time increases. On the other side, the controlling signal has been decreasing by increasing the value of the proportional constant that takes part in PID and neglected integral and derivative parts. Whereas in the second case, composite laminates consist of four plies (0/90/0/90), the closed-loop responses or all proportional gain (kp) shown as in Fig. 5(b) is close to the results of the first case because the effect of the fourth ply is smaller (90 degree) composite ply with the beam in the fourth case, which demonstrates that control responses of the system will be near each other and this reduction in the time compared with others is due to the increase in the number of composite layers and at the same time increases the strength properties of the beam. Based on the results, it has been noted that an increase in laminates of composite plate decreases the vibration and time stability of the system.</span>
+However, the evidence for this relationship is conclusive and has been retrieved from PID rules; with increasing the proportional range, the vibration converges slowly and settling time increases. On the other side, the controlling signal has been decreasing by increasing the value of the proportional constant that takes part in PID and neglected integral and derivative parts. Whereas in the second case, composite laminates consist of four plies (0/90/0/90), the closed-loop responses or all proportional gain (<math>kp</math>) shown as in [[#img-5|Figure 5]](b) is close to the results of the first case because the effect of the fourth ply is smaller (90 degree) composite ply with the beam in the fourth case, which demonstrates that control responses of the system will be near each other and this reduction in the time compared with others is due to the increase in the number of composite layers and at the same time increases the strength properties of the beam. Based on the results, it has been noted that an increase in laminates of composite plate decreases the vibration and time stability of the system.
-For amplitude comparison between two sets of composite materials beam structures, first set consist of (0-90) degrees of laminates and another set consist of (0-45) degrees of laminates. This based on the responses that were simulated with high accuracy using the proportional gain formulation. As shown in Figs. 5 and 6, all curves have nearly the same shape, but speed to reach the stable case in the second set is smaller than the first set. The reason is attributed to the greater flexural strength for the second compared to the former.
+For amplitude comparison between two sets of composite materials beam structures, first set consist of (0-90) degrees of laminates and another set consist of (0-45) degrees of laminates. This based on the responses that were simulated with high accuracy using the proportional gain formulation. As shown in [[#img-5|Figures 5]] and [[#img-6|6]], all curves have nearly the same shape, but speed to reach the stable case in the second set is smaller than the first set. The reason is attributed to the greater flexural strength for the second compared to the former.
+<div id='img-6'></div>
-Piezoelectric MFC has a maximum voltage, so it is peak to peak voltage. While peak voltage is the voltage from the zero horizontal axis line to the top of the waveform, the peak to peak voltage is the voltage of the one cycle of a wave, all the way from the peak of the negative side to the peak of the positive side [37]. The maximum voltage for MFC to be applied is (2000 volt) for the piezoelectric actuator. The uncontrolled and controlled responses obtained by actuation voltage can be seen in Figs. 7(a), 7(b), 7(c), 7(d), 7(e), 7(f), 7(g), and 7(h). From the graph, it has been observed that in the open-loop, the response was nil because in the open-loop, there is no link between vibrate response and actuator.
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
-{| style="width: 100%;border-collapse: collapse;"
-|-
-|  style="border: 1pt solid black;text-align: center;width: 50%;"|[[Image:Review_223450771672-image11.jpeg|282px]]
-|  style="border: 1pt solid black;text-align: center;width: 50%;"|[[Image:Review_223450771672-image12.jpeg|282px]]
 |-
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(a) CFRP (0/45/0)</span>
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image11.jpeg|382px]]
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(b) CFRP (0/45/0/45)</span>
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image12.jpeg|382px]]
 |-
-|  style="border: 1pt solid black;text-align: center;width: 50%;"|[[Image:Review_223450771672-image13.jpeg|288px]]
+|  style="text-align: center;font-size: 75%;"|(a) CFRP (0/45/0)
-|  style="border: 1pt solid black;text-align: center;width: 50%;"|[[Image:Review_223450771672-image14.jpeg|282px]]
+|  style="text-align: center;font-size: 75%;"|(b) CFRP (0/45/0/45)
 |-
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(c) CFRP (0/45/45/0)</span>
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image13.jpeg|388px]]
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(d) CFRP (0/45/0/45/0/45)</span>
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image14.jpeg|382px]]
 |-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Fig. 6 Open-loop and closed-loop responses of the composite beam of tip displacement</span>
+|  style="text-align: center;font-size: 75%;"|(c) CFRP (0/45/45/0)
+|  style="text-align: center;font-size: 75%;"|(d) CFRP (0/45/0/45/0/45)
+|- style="text-align: center; font-size: 75%;"
+| colspan="2" style="padding:10px;"| '''Figure 6'''. Open-loop and closed-loop responses of the composite beam of tip displacement
 |}
-<span style="text-align: center; font-size: 100%;">On the other hand, in all cases, when the control signal applied with the gain to a system increases, the rise time generally decreases as shown in Fig. 7, where the actuator voltages for different control gain for all cases are observed to increase with a proportional constant. Meanwhile, kp=8 has the highest value of actuator voltage observed in all cases [38].</span>
+Piezoelectric MFC has a maximum voltage, so it is peak to peak voltage. While peak voltage is the voltage from the zero horizontal axis line to the top of the waveform, the peak to peak voltage is the voltage of the one cycle of a wave, all the way from the peak of the negative side to the peak of the positive side [37]. The maximum voltage for MFC to be applied is (2000 volt) for the piezoelectric actuator. The uncontrolled and controlled responses obtained by actuation voltage can be seen in [[#img-7|Figures 7]](a), (b), (c), (d), (e), (f), (g), and (h). From the graph, it has been observed that in the open-loop, the response was nil because in the open-loop, there is no link between vibrate response and actuator.
-As can be seen from Figs. 7(b), 7(c), 7(d), 7(f), 7(g), and 7(h), the minimum value of actuator voltage is (-500V) and smaller than this volt by APDL-ANSYS macro file has been changed to (-500 Volt) because the sinusoidal voltage of piezoelectric MFC-P1 is worked between (-500 volt) to (+1500 Volt) voltage cycle at 0.1 Hz was output by a voltage amplifier [39].
+On the other hand, in all cases, when the control signal applied with the gain to a system increases, the rise time generally decreases as shown in [[#img-7|Figure 7]], where the actuator voltages for different control gain for all cases are observed to increase with a proportional constant. Meanwhile, <math>kp=8</math> has the highest value of actuator voltage observed in all cases [38].
-{| style="width: 100%;border-collapse: collapse;"
+As can be seen from [[#img-7|Figures 7]](b), (c), (d), (f), (g), and (h), the minimum value of actuator voltage is (-500V) and smaller than this volt by APDL-ANSYS macro file has been changed to (-500 Volt) because the sinusoidal voltage of piezoelectric MFC-P1 is worked between (-500 volt) to (+1500 Volt) voltage cycle at 0.1 Hz was output by a voltage amplifier [39].
+<div id='img-7'></div>
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|  style="border: 1pt solid black;text-align: center;width: 44%;"|[[Image:Review_223450771672-image15.jpeg|240px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image15.jpeg|340px]]
-|  style="border: 1pt solid black;text-align: center;width: 55%;"|[[Image:Review_223450771672-image16.jpeg|240px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image16.jpeg|340px]]
 |-
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(a) CFRP (0/90/0)</span>
+|  style="text-align: center;font-size: 75%;"|(a) CFRP (0/90/0)
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(b)  CFRP (0/90/0/90)</span>
+|  style="text-align: center;font-size: 75%;"|(b) CFRP (0/90/0/90)
 |-
-|  style="border: 1pt solid black;text-align: center;width: 44%;"|[[Image:Review_223450771672-image17.jpeg|240px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image17.jpeg|340px]]
-|  style="border: 1pt solid black;text-align: center;width: 55%;"|[[Image:Review_223450771672-image18.jpeg|240px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image18.jpeg|340px]]
 |-
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(c) CFRP (0/90/90/0)</span>
+|  style="text-align: center;font-size: 75%;"|(c) CFRP (0/90/90/0)
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(d) CFRP (0/90/0/0/90/0)</span>
+|  style="text-align: center;font-size: 75%;"|(d) CFRP (0/90/0/0/90/0)
 |-
-|  style="border: 1pt solid black;text-align: center;width: 44%;"|[[Image:Review_223450771672-image19.jpeg|240px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image19.jpeg|340px]]
-|  style="border: 1pt solid black;text-align: center;width: 55%;"|[[Image:Review_223450771672-image20.jpeg|240px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image20.jpeg|340px]]
 |-
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(e) CFRP (0/45/0)</span>
+|  style="text-align: center;font-size: 75%;"|(e) CFRP (0/45/0)
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(f) CFRP (0/45/0/45)</span>
+|  style="text-align: center;font-size: 75%;"|(f) CFRP (0/45/0/45)
 |-
-|  style="border: 1pt solid black;text-align: center;width: 44%;"|[[Image:Review_223450771672-image21.jpeg|240px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image21.jpeg|340px]]
-|  style="border: 1pt solid black;text-align: center;width: 55%;"|[[Image:Review_223450771672-image22.jpeg|234px]]
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image22.jpeg|334px]]
 |-
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(g) CFRP (0/45/45/0)</span>
+|  style="text-align: center;font-size: 75%;"|(g) CFRP (0/45/45/0)
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(h) CFRP (0/45/0/0/45/0)</span>
+|  style="text-align: center;font-size: 75%;"|(h) CFRP (0/45/0/0/45/0)
-|-
+|- style="text-align: center; font-size: 75%;"
-|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Fig. 7  Open loop and closed loop responses of composite beam of actuation voltage.</span>
+| colspan="2" style="padding:10px;"| '''Figure 7'''. Open-loop and closed-loop responses of the composite beam of the actuation voltage
 |}
-<span style="text-align: center; font-size: 100%;">Another important finding was the effect of fiber orientations of composite laminate layers on active vibrations, as shown in Fig.8. The composite beams composed of laminates (0/90/0), (0/90/0/90), (0/90/90/0), and (0/90/0/0/90/0) reach stability faster than the composite beams composed of laminates (0/45/0), (0/45/0/45), (0/45/45/0), and (0/45/0/0/45/0) because the first set has higher relaxation modulus than second sets [40].</span>
+Another important finding was the effect of fiber orientations of composite laminate layers on active vibrations, as shown in [[#img-8|Figure 8]]. The composite beams composed of laminates (0/90/0), (0/90/0/90), (0/90/90/0), and (0/90/0/0/90/0) reach stability faster than the composite beams composed of laminates (0/45/0), (0/45/0/45), (0/45/45/0), and (0/45/0/0/45/0) because the first set has higher relaxation modulus than second sets [40].
+<div id='img-8'></div>
-{| style="width: 100%;border-collapse: collapse;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|  colspan='2'  style="border-top: 1pt solid black;"|
+|  style="text-align: center;padding:10px;"| [[Image:Review_223450771672-image23.jpeg|300px]]
+|  style="text-align: center;padding:10px;"|[[Image:Review_223450771672-image24.jpeg|294px]]
 |-
-|  style="border: 1pt solid black;width: 51%;"|[[Image:Review_223450771672-image23.jpeg|300px]]
+|  style="text-align: center;font-size: 75%;"|(a)
-|  style="border: 1pt solid black;width: 48%;"|[[Image:Review_223450771672-image24.jpeg|294px]]
+|  style="text-align: center;font-size: 75%;"|(b)
 |-
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(a)</span>
+|  style="text-align: center;padding:10px;"|[[Image:Review_223450771672-image25.jpeg|294px]]
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(b)</span>
+|  style="text-align: center;padding:10px;"|[[Image:Review_223450771672-image26.jpeg|294px]]
 |-
-|  style="border: 1pt solid black;width: 51%;"|[[Image:Review_223450771672-image25.jpeg|294px]]
+|  style="text-align: center;font-size: 75%;"|(c)
-|  style="border: 1pt solid black;width: 48%;"|[[Image:Review_223450771672-image26.jpeg|294px]]
+|  style="text-align: center;font-size: 75%;"|(d)
-|-
+|- style="text-align: center; font-size: 75%;"
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(c)</span>
+| colspan="2" style="padding:10px;"| '''Figure 8'''. Comparision between two sets of similar composite beam
-|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">(d)</span>
-|-
-|  colspan='2'  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Fig. 8 Comparision between two sets of similar composite beam</span>
 |}
 ==5. Conclusions==
-<span style="text-align: center; font-size: 100%;">  A FEM for piezoelectric with composite beam structure was proposed to simulate model dynamic, and the electrical model of piezoelectric which permitted the simulations of type ( </span> <math display="inline">{d}_{33}-</math><math>MFC</math><span style="text-align: center; font-size: 100%;"> ) effects. The performance of the composite beam structure is improved when piezoelectric type M8528(P1) is used. It is shown that the classic controller provides in increasing control gain (kp) and better active damping than low control gain, especially at kp=8.</span>
+A FEM for piezoelectric with composite beam structure was proposed to simulate model dynamic, and the electrical model of piezoelectric which permitted the simulations of type (<math display="inline">{d}_{33}-</math><math>MFC</math>) effects. The performance of the composite beam structure is improved when piezoelectric type M8528(P1) is used. It is shown that the classic controller provides in increasing control gain (<math display="inline">kp</math>) and better active damping than low control gain, especially at <math display="inline">kp=8</math>.
-<span id='_GoBack'>The active control capability of such a composite beam is analyzed using a simple control strategy. Overall, these controlled responses show that, although all different composite laminate plates are subjected to the same intensity of mechanical load. The observed decrease in the vibration of the composite plate might be explained by the increase in the proportional gain in PID control on all cases of the composite laminate plate, while at the same time decreasing the voltage actuator of piezoelectric. The model can be used for more investigations with general composition, applied loads, and another control method. The final objectives of analysis and modeling are the minimization of weight (or mass) of the beam structure under behavior constraints, increasing the performance of the piezoelectric actuator, and maximization of natural frequencies of the vibration modes. However, given the immersed short time for each case, caution must be taken. The difference between two similar different in orientation laminate, two refined finite element models based 0/90 laminate set and 0/45 laminate set, for the analysis of laminated beam with embedded surface bonded piezoelectric laminae have been presented. The evidence from this study suggests the composite beams composed from (0/90) stacking sequence sets had reached stability state faster than the composite beams composed from (0/45) stacking sequencesets.</span>
+The active control capability of such a composite beam is analyzed using a simple control strategy. Overall, these controlled responses show that, although all different composite laminate plates are subjected to the same intensity of mechanical load. The observed decrease in the vibration of the composite plate might be explained by the increase in the proportional gain in PID control on all cases of the composite laminate plate, while at the same time decreasing the voltage actuator of piezoelectric. The model can be used for more investigations with general composition, applied loads, and another control method. The final objectives of analysis and modeling are the minimization of weight (or mass) of the beam structure under behavior constraints, increasing the performance of the piezoelectric actuator, and maximization of natural frequencies of the vibration modes. However, given the immersed short time for each case, caution must be taken. The difference between two similar different in orientation laminate, two refined finite element models based 0/90 laminate set and 0/45 laminate set, for the analysis of laminated beam with embedded surface bonded piezoelectric laminae have been presented. The evidence from this study suggests the composite beams composed from (0/90) stacking sequence sets had reached stability state faster than the composite beams composed from (0/45) stacking sequencesets.
 ==References==
+<div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;">
-[1] Sun, L., Li, W., Wu, Y. and Lan, Q. Active vibration control of a conical shell using piezoelectric ceramics. Journal of Low Frequency Noise, Vibration and Active Control, 36(4), pp.366-375. ,2017.
-[2] Baz, A.M.. Active and Passive Vibration Damping. John Wiley & Sons. , 2019.
+[1] Sun L., Li W., Wu Y., Lan Q. Active vibration control of a conical shell using piezoelectric ceramics. Journal of Low Frequency Noise, Vibration and Active Control, 36(4):366-375, 2017.
-[3] Kouider, B. and Polat, A. Numerical investigation of active vibration control of composite plate using PID controller. 2017.
+[2] Baz A.M.. Active and passive vibration damping. John Wiley & Sons, 2019.
-[4] Symans, M.D. and Constantinou, M.C. Semi-active control systems for seismic protection of structures: a state-of-the-art review. Engineering structures, 21(6), pp.469-487, 1999.
+[3] Kouider B., Polat A. Numerical investigation of active vibration control of composite plate using PID controller. International Conference on Advances and Innovations in Engineering, 2017.
-[5] Kumar, V. Experimental investigation on active vibration control of smart composite shell structures. PhD thesis ,Jawaharlal Nehru Technological University Hyderabad, 2011.
+[4] Symans M.D., Constantinou M.C. Semi-active control systems for seismic protection of structures: A state-of-the-art review. Engineering structures, 21(6):469-487, 1999.
-[6] Manning, W.J., Plummer, A.R. and Levesley, M.C. Vibration control of a flexible beam with integrated actuators and sensors . Smart Materials and Structures, 9(6), p.932, 2000.
+[5] Kumar V. Experimental investigation on active vibration control of smart composite shell structures. PhD Thesis, Jawaharlal Nehru Technological University Hyderabad, 2011.
-[7] Gaudenzi, P., Carbonaro, R. and Benzi, E.   Control of beam vibrations by means of piezoelectric devices: theory and experiments . Composite structures, 50(4), pp.373-379, 2000.
+[6] Manning W.J., Plummer A.R., Levesley M.C. Vibration control of a flexible beam with integrated actuators and sensors. Smart Materials and Structures, 9(6):932, 2000.
-[8] Bruant, I., Coffignal, G., Lene, F. and Verge, M.  Active control of beam structures with piezoelectric actuators and sensors: modeling and simulation . Smart materials and structures, 10(2), p.404, 2001.
+[7] Gaudenzi P., Carbonaro R., Benzi E.   Control of beam vibrations by means of piezoelectric devices: theory and experiments. Composite Structures, 50(4):373-379, 2000.
-[9] Singh, S.P., Pruthi, H.S. and Agarwal, V.P. Efficient modal control strategies for active control of vibrations .  Journal of Sound and Vibration, 262(3), pp.563-575,2003.
+[8] Bruant I., Coffignal G., Lene F., Verge M.  Active control of beam structures with piezoelectric actuators and sensors: modeling and simulation. Smart Materials & Structures, 10(2):404, 2001.
-[10] Xu, S.X. and Koko, T.S.  Finite element analysis and design of actively controlled piezoelectric smart structures . Finite elements in analysis and design, 40(3), pp.241-262 , 2004.
+[9] Singh S.P., Pruthi H.S., Agarwal V.P. Efficient modal control strategies for active control of vibrations.  Journal of Sound and Vibration, 262(3):563-575, 2003.
-[11] Lim, Y.H. Finite-element simulation of closed loop vibration control of a smart plate under transient loading . Smart Materials and Structures, 12(2), p.272, 2003.
+[10] Xu S.X., Koko T.S.  Finite element analysis and design of actively controlled piezoelectric smart structures. Finite Elements in Analysis and Design, 40(3):241-262, 2004.
-[12] Quek, S.T., Wang, S.Y. and Ang, K.K. Vibration control of composite plates via optimal placement of piezoelectric patches . Journal of intelligent material systems and structures, 14(4-5), pp.229-245, 2003.
+[11] Lim Y.H. Finite-element simulation of closed loop vibration control of a smart plate under transient loading. Smart Materials and Structures, 12(2):272, 2003.
-[13] Xianmin, Z., Changjian, S. and Erdman, A.G.  Active vibration controller design and comparison study of flexible linkage mechanism systems . Mechanism and Machine Theory, 37(9), pp.985-997, 2002.
+[12] Quek S.T., Wang S.Y., Ang, K.K. Vibration control of composite plates via optimal placement of piezoelectric patches. Journal of Intelligent Material Systems and Structures, 14(4-5):229-245, 2003.
-[14] Halim, D. and Moheimani, S.R.  An optimization approach to optimal placement of collocated piezoelectric actuators and sensors on a thin plate . Mechatronics, 13(1), pp.27-47, 2003.
+[13] Xianmin Z., Changjian S., Erdman A.G.  Active vibration controller design and comparison study of flexible linkage mechanism systems. Mechanism and Machine Theory, 37(9):985-997, 2002.
-[15] Bendine, K., Boukhoulda, F.B., Haddag, B. and Nouari, M. Active vibration control of composite plate with optimal placement of piezoelectric patches. Mechanics of Advanced Materials and Structures, 26(4), pp.341-349, 2019.
+[14] Halim D., Moheimani S.R.  An optimization approach to optimal placement of collocated piezoelectric actuators and sensors on a thin plate. Mechatronics, 13(1):27-47, 2003.
-[16] Malgaca, L. Integration of active vibration control methods with finite element models of smart laminated composite structures. Composite Structures, 92(7), pp.1651-1663, 2010.
+[15] Bendine K., Boukhoulda F.B., Haddag B., Nouari M. Active vibration control of composite plate with optimal placement of piezoelectric patches. Mechanics of Advanced Materials and Structures, 26(4):341-349, 2019.
-[17] Takacs, G. and Rohal-Ilkiv, B., April. Direct closed-loop active vibration control system prototype in ANSYS. In INTER-NOISE and NOISE-CON Congress and Conference Proceedings (Vol. 244, No. 1, pp. 1-12). Institute of Noise Control Engineering , 2012.
+[16] Malgaca L. Integration of active vibration control methods with finite element models of smart laminated composite structures. Composite Structures, 92(7):1651-1663, 2010.
-[18] Vashist, S.K. & Chhabra, D. Optimal placement of piezoelectric actuators on plate structures for active vibration control using genetic algorithm , In Active and Passive Smart Structures and Integrated Systems  , International Society for Optics and Photonics,  , 9057, p. 9, 2014.
+[17] Takacs G., Rohal-Ilkiv B. Direct closed-loop active vibration control system prototype in ANSYS. In INTER-NOISE and NOISE-CON Congress and Conference Proceedings, 244(1):1-12, Institute of Noise Control Engineering, April 2012.
-[19] Yavuz, Ş., Malgaca, L. and Karagülle, H.  Analysis of active vibration control of multi-degree-of-freedom flexible systems by Newmark method . Simulation Modelling Practice and Theory, 69, pp.136-148 , 2016.
+[18] Vashist S.K.,  Chhabra D. Optimal placement of piezoelectric actuators on plate structures for active vibration control using genetic algorithm. In Active and Passive Smart Structures and Integrated Systems, International Society for Optics and Photonics, 9057:9, 2014.
-[20] Liu, X., Cai, G., Peng, F. and Zhang, H.  Piezoelectric actuator placement optimization and active vibration control of a membrane structure. Acta Mechanica Solida Sinica, 31(1), pp.66-79.  [https://doi.org/ https://doi.org/] 10.1007/s10338-018-0005-y, 2018.
+[19] Yavuz Ş., Malgaca L., Karagülle H.  Analysis of active vibration control of multi-degree-of-freedom flexible systems by Newmark method. Simulation Modelling Practice and Theory, 69:136-148, 2016.
-[21] Van Tham, V., Quoc, T.H. and Tu, T.M.  Optimal placement and active vibration control of composite plates integrated piezoelectric sensor/actuator pairs . Vietnam Journal of Science and Technology, 56(1), p.113, 2018.
+[20] Liu X., Cai G., Peng F., Zhang H.  Piezoelectric actuator placement optimization and active vibration control of a membrane structure. Acta Mechanica Solida Sinica, 31(1):66-79, 2018.
-[22] Baghaee, M., Farrokhabadi, A. and Jafari-Talookolaei, R.A. Modeling, analysis, and control of MFC sandwiched laminate panel flutter with general layups and arbitrary boundary conditions. Composite Structures, 223, p.110940, 2019.
+[21] Van Tham V., Quoc T.H., Tu T.M.  Optimal placement and active vibration control of composite plates integrated piezoelectric sensor/actuator pairs. Vietnam Journal of Science and Technology, 56(1):113, 2018.
-[23] Usman, M., Hanif, A., Kim, I.H. and Jung, H.J.  Experimental validation of a novel piezoelectric energy harvesting system employing wake galloping phenomenon for a broad wind spectrum. Energy, 153, pp.882-889, 2018.
+[22] Baghaee M., Farrokhabadi A., Jafari-Talookolaei R.A. Modeling, analysis, and control of MFC sandwiched laminate panel flutter with general layups and arbitrary boundary conditions. Composite Structures, 223:110940, 2019.
-[24] Heganna, S.S. and Joglekar, J.J., 2016. Active vibration control of smart structure using PZT patches. Procedia Computer Science, 89, pp.710-715.
+[23] Usman M., Hanif A., Kim I.H., Jung H.J.  Experimental validation of a novel piezoelectric energy harvesting system employing wake galloping phenomenon for a broad wind spectrum. Energy, 153:882-889, 2018.
-[25] Xie, C., Wu, Y. and Liu, Z., 2018. Modeling and active vibration control of lattice grid beam with piezoelectric fiber composite using fractional order PDμ algorithm. Composite Structures, 198, pp.126-134.
+[24] Heganna S.S., Joglekar J.J. Active vibration control of smart structure using PZT patches. Procedia Computer Science, 89:710-715, 2016.
-[26] Rimašauskienė, R., Jūrėnas, V., Radzienski, M., Rimašauskas, M. and Ostachowicz, W., 2019. Experimental analysis of active–passive vibration control on thin-walled composite beam. Composite Structures, 223, p.110975.
+[25] Xie C., Wu Y., Liu Z. Modeling and active vibration control of lattice grid beam with piezoelectric fiber composite using fractional order PDμ algorithm. Composite Structures, 198:126-134, 2018.
-[27] IEEE, 1984. An American National Standard: IEEE Standard on Piezoelectricity. IEEE Transactions on Sonics and Ultrasonics, (2), pp.1-55.
+[26] Rimašauskienė R., Jūrėnas V., Radzienski M., Rimašauskas M., Ostachowicz W. Experimental analysis of active–passive vibration control on thin-walled composite beam. Composite Structures, 223:110975, 2019.
-[28] Abuzaid, A., Hrairi, M. and Dawood, M.S.I., 2015, June. Survey of active structural control and repair using piezoelectric patches. In Actuators (Vol. 4, No. 2, pp. 77-98). Multidisciplinary Digital Publishing Institute.
+[27] IEEE. An American National Standard: IEEE Standard on Piezoelectricity. IEEE Transactions on Sonics and Ultrasonics, (2):1-55, 1984.
-[29] Heywang, W., Lubitz, K. and Wersing, W. eds., 2008. Piezoelectricity: evolution and future of a technology (Vol. 114). Springer Science & Business Media.
+[28] Abuzaid A., Hrairi M., Dawood M.S.I. Survey of active structural control and repair using piezoelectric patches. Actuators 4(2):77-98, Multidisciplinary Digital Publishing Institute, June 2015.
-[30] ANSYS, ANSYS Mechanical APDL Basic Analysis Guide/ANSYS help, release 15.0 2013, Inc., Canonsburg, PA, 2017.
+[29] Heywang W., Lubitz K., Wersing W. Eds. Piezoelectricity: evolution and future of a technology. Springer Science & Business Media, Springer Series in Materials Science, Vol. 114, 2008.
-[31] Kamel, M.A., Ibrahim, K. and Ahmed, A.E.M. Vibration control of smart cantilever beam using finite element method. Alexandria Engineering Journal, 58(2), pp.591-601,2019.
+[30] ANSYS, ANSYS Mechanical APDL Basic Analysis Guide/ANSYS help, release 15.0 2013. Inc. Canonsburg, PA, 2017.
-[32] Trindade, M.A. and Benjeddou, A.  Finite element characterisation of multilayer d31 piezoelectric macro-fibre composites . Composite Structures, 151, pp.47-57,2016.
+[31] Kamel M.A., Ibrahim K., Ahmed A.E.M. Vibration control of smart cantilever beam using finite element method. Alexandria Engineering Journal, 58(2):591-601, 2019.
-[33] Vijay Kumar, G., Raja, S., Prasanna, K.B. and Sudha, V. Finite element analysis and vibration control of a deep composite cylindrical shell using MFC actuators . Smart Materials Research, 2012.
+[32] Trindade M.A., Benjeddou A.  Finite element characterisation of multilayer d31 piezoelectric macro-fibre composites. Composite Structures, 151:47-57, 2016.
-[34] Smart Materials Corporation. (Accessed Nov 2019). Actuator and sensor ,Available online: [https://www.smart-material.com/MFC-product-P1.html https://www.smart-material.com/MFC-product-P1.html]
+[33] Vijay Kumar G., Raja S., Prasanna K.B., Sudha V. Finite element analysis and vibration control of a deep composite cylindrical shell using MFC actuators. Smart Materials Research, 2012.
-[35] Jimin He and Zhi-Fang Fu , Modal Analysis, Linacre House, Jorden Hill, Oxford , 2001.
+[34] Smart Materials Corporation (Accessed Nov 2019). Actuator and sensor ,Available online: [https://www.smart-material.com/MFC-product-P1.html https://www.smart-material.com/MFC-product-P1.html]
-[36] Huang, B., Kim, H.S. and Youn, B.D.  Active vibration control of smart composite laminates with partial debonding of actuator . International Journal of Precision Engineering and Manufacturing, 16(4), pp.831-840, 2015.
+[35] He J., Fu Z.F., Modal Analysis. Linacre House, Jorden Hill, Oxford, 2001.
-[37] Ian Sinclair . Passive Components for Circuit Design, Linacre House, Jordan Hill, Oxford OX2 8DP, 2004.
+[36] Huang B., Kim H.S., Youn B.D.  Active vibration control of smart composite laminates with partial debonding of actuator. International Journal of Precision Engineering and Manufacturing, 16(4):831-840, 2015.
-[38] Garcia-Segura, T., Yepes, V., Marti, J.V. and Alcala, J. Optimization of concrete I-beams using a new hybrid glowworm swarm algorithm . Latin American Journal of Solids and Structures, 11(7), pp.1190-1205, 2014.
+[37] Ian Sinclair. Passive Components for Circuit Design. Linacre House, Jordan Hill, Oxford OX2 8DP, 2004.
-[39] Wu, M., Yuan, X., Luo, H., Chen, H., Chen, C., Zhou, K. and Zhang, D. Enhanced actuation performance of piezoelectric fiber composites induced by incorporated BaTiO3 nanoparticles in epoxy resin . Physics Letters A, 381(19), pp.1641-1647, 2017.
+[38] Garcia-Segura T., Yepes V., Marti J.V., Alcala J. Optimization of concrete I-beams using a new hybrid glowworm swarm algorithm. Latin American Journal of Solids and Structures, 11(7):1190-1205, 2014.
-[40] Rony Das. Stress relaxation of ehermostting prepreg composites during cure  ,MSc. thesis , Bangladesh University of Engineering and Technology, 2008
+[39] Wu M., Yuan X., Luo H., Chen H., Chen C., Zhou K., Zhang D. Enhanced actuation performance of piezoelectric fiber composites induced by incorporated BaTiO3 nanoparticles in epoxy resin. Physics Letters A, 381(19):1641-1647, 2017.
+[40] Das R. Stress relaxation of ehermostting prepreg composites during cure. MSc. Thesis, Bangladesh University of Engineering and Technology, 2008.
+</div>