Research Article  Open Access
M. E. Fares, M. Kh. Elmarghany, Doaa Atta, "Active Control Optimization for Minimizing the Dynamic Response of Functionally Graded Plate", International Journal of Aerospace Engineering, vol. 2014, Article ID 516303, 15 pages, 2014. https://doi.org/10.1155/2014/516303
Active Control Optimization for Minimizing the Dynamic Response of Functionally Graded Plate
Abstract
A control optimization problem for functionally graded (FG) plates is presented using a firstorder shear deformation plate theory including throughthethickness normal strain effect. The aim of the optimization is to minimize the vibrational response of FG plate with constraints on the control energy used in the damping process. An active control optimization is presented to determine the optimal level of a closed loop control function. Plate thickness and a homogeneity parameter of FG plates are used as design variables. Numerical results for the optimal control force and the total energy for a simply supported FG plate are given. The influence of throughthethickness normal strain effect on the accuracy of the obtained results is illustrated. The effectiveness of the present control and design procedure is examined.
1. Introduction
The rapid development in various industrial fields requires new materials that can serve useful functions under certain conditions. In recent years, functionally graded materials (FGMs) have gained considerable attention in many engineering applications, and they are considered as a potential structural material for future highspeed spacecraft and power generation industries. Several studies have been performed to analyze the behavior of functionally graded plates [1].
In the aerospace industry and many other engineering applications, the suppression of excessive vibrations occurring in large structures represents one of the most pressing and difficult problems facing structural designers. An effective means of suppressing excessive vibrations is by active structural control. Thus, there is need for new light materials possessing a high degree of flexibility and with very low natural damping. These factors motivated the development for more accurate tools of analysis and rigorous design methods [2â€“7]. A series of publications has appeared concerning the vibration control and design optimization of composite structures. Grandhi [8] studied the structural and control optimization of space structures. Sloss et al. [9] presented an integrated approach to solve multiobjective design and control optimization for many composite structures. Adali et al. [10, 11], Sloss et al. [12], Sadek et al. [13], and Fares et al. [14â€“16] have studied the optimal design and control optimization of composite laminated plates and shells subjected to thermomechanical loadings.
In another work of Fares et al. [17], nonlinear optimization schemes are presented to solve multiobjective design and control problems for composite laminated plates. Other studies on this topic may be found in [18, 19].
Many studies indicate that the inclusion of the normal strain and the transverse shear effects are very important for predicting accurate vibrational responses of composite plates and shells, particularly, for those made of FGMs; see, for example, [20]. In the work [20], it is showed that the neglecting of these effects may cause errors in predicting the deflections and frequencies in FG plates reaching 30%. However, most studies related to the topic of design and control optimization of composite plates and shells were carried out based on classical theories which neglect the shear deformation and normal strain effects.
The current work deals with the active control of the dynamic response of a FG orthotropic rectangular plate. The present formulation is based on a consistent firstorder plate theory including the shear deformation and throughthethickness normal strain effects [20]. The design and control objective are to minimize the dynamic response with minimum expenditure of control energy. The total elastic energy is taken as a measure of the dynamic response. The expenditure of the control energy is limited by attaching to the control objective a functional including the control force. The necessary and sufficient condition for optimal stabilization in the LiapunovBellman sense [21] is used to determine the optimal control force and controlled deflections. Numerical examples are given to assess the efficiency of the present design control approach in the damping process for FG plate.
2. Formulation and Basic Equations
Consider an elastic plate of uniform thickness , length , and width . The plate is composed of an orthotropic elastic material varying in the thickness direction. Let a Cartesian coordinates system be taken such that the midplane coincides with the plane and the plate occupies the following region:
Let the upper surface of the plate () be subjected to a transverse distributed load which may be taken as a control force.
The present formulation is based on a firstorder theory including a normal strain and shear deformation effects derived using a mixed variational approach. This theory preserves the transverse shear stresses which vanish on the top and the bottom of the plate surfaces; therefore, the condition of introducing shear correction factors is not required. The displacement field is taken in the form [20]: where (, , ) are the displacements along , , and directions, respectively, (, , ) are the displacements of a point on the midplane, and (, ) are the slopes in the and planes due to bending only (slope rotations). The constants and may take the values 0 or 1 to study the effect of the normal strain on the control process.
The infinitesimal strains associated with displacements (2) are given by where the suffix () denotes partial differentiation with respect to the coordinate .
Moreover, the stresses are taken in the following form:
The governing equations of this system are determined using the dynamic modified version of the mixed variational principle of Ressiner [20]:
Moreover, the constitutive equations are
The stress and moment resultants , , and in (5) and (6a) and (6b) can be expressed as and the stiffness coefficients , , and are given in Appendix A.
In the present control problem, all edges of the plate are taken simply supported (SSSS), and then, the following conditions must be satisfied at the plate edges:
3. The Optimal Design and Control Problem
The objective of the present control problem is to minimize the dynamic response of the plate in a specified time with the minimum possible expenditure of force , when the plate is subjected to the following initial disturbances:
The dynamic response of the plate is measured by a cost functional including the energy of the plate which is a function of displacements, its derivatives, and the velocity. The control force is introduced into the formulation by taking a performance index which comprises a weighted sum of the energy (cost functional) and penalty functional involving the control energy. Then, the mathematical formulation of the control problem can be reduced to determine the optimization variable that minimizes the following control objective: where , , and are positive constant weighting factors and and represent the strain and the kinetic energies of the plate, respectively. The functional is a penalty term involving the control function , where denotes the set of all bounded square integrable functions on the domain of solution . Note that, the functional is a quadratic and, hence, differentiable and positive definite on the domain of solution.
It is obvious that the solution of the present optimal control problem requires finding a solution of the initial boundaryvalue problem (5), (6a), (6b), (8), and (9) that minimizes the functional (10). Moreover, the design optimization procedure aims to minimize the actively controlled dynamic responses using geometry and material parameters as design variables.
4. Solution Procedure
The displacement functions (, , , , , , ) and the mechanical load may be taken for simply supported plate in the form: where , , , and are unknown functions of the time .
Using (6a) and (6b) we can get the governing equations (5) in terms of the displacements (, , , , , , ). For these equations, the inplane inertia terms are neglected. Substituting expressions (11) into the resulting equations, we get the following time equations: where is the kronecker delta and the coefficients , , , , , , , are given in Appendix B. Solving the system of (12), one gets an equation of the timedependent functions and only, where , , and are given in Appendix B.
Again, using (2), (3), (4), and (11) with the objective functional (10), we get where the coefficients () are given in Appendix B.
Now the optimal control problem is to find the control function that minimizes the functional (14); that is, that is,
The minimization of the dynamic response using the control force can be carried out using the LiapunovBellman method [21] which is considered an effective approach for the solution. The necessary and sufficient condition for minimization of the functional (14) in the LiapunovBellman sense is provided that the Liapunov function is a positive definite, then Substituting expressions (13) and (14) into (17) one gets This condition gives that
Substituting the Liapunov function (18) into the previous expression, we can get the minimum control force in the form: Again, using (18) with the expression (19), we obtain
Using (21) with the last expression and equating the coefficients of , , and by zeroes, the following system of algebraic equations is obtained: where the coefficients () are given in Appendix C.
Using (14) and (21), we obtain the following differential equation: The solution of (24) for which is given by where and are unknown coefficients which may be obtained from conditions (9) as
Using (14), (21), and (25), we can get the controlled displacement, the total energy, and the optimal control force. Then we complete the minimization process by determining the optimal design of the plate using geometric and material parameters as design variables.
5. Homogenization of Material Properties and the Design Objective
Consider that a FG plate is made of two randomly distributed isotropic constituents (ceramic and metal) and the macroscopic response of the composite is isotropic. The matrix phase, denoted by the subscript 2, is reinforced by spherical particles as a particulate phase, denoted by the subscript 1. The homogeneity of this composition is taken to vary only in the thickness direction. The effective elastic moduli of the graded composite may be evaluated from the volume fraction of the constituents and . The effective mass density at a point is given by where and are, respectively, the volume fractions of the particulate ceramic phase on the top and the bottom surfaces of the plate. The parameter indicates the volume fraction profile through the thickness. In this case, the elastic coefficients become where is Youngâ€™s modulus, is Poissonâ€™s ratio, is shear modulus, and is the effective bulk modulus.
The Mori and Tanaka [22] scheme is used to estimate the effective moduli for the regions of the graded microstructure which have a welldefined continuous matrix and a discontinuous particulate phase. The locally effective bulk modulus and the effective shear modulus of the composite are given by
The design procedure aims to determine the values of the design variable () and the index that minimize the energy ; that is, where is the optimal control force obtained from the previous control procedure.
6. Numerical Examples and Discussion
Numerical results are presented for the total energy and the control force for simply supported plates. Let the constituent materials of the FG plate be aluminum (Al) as a matrix phase and zirconia as particulate phase with material and geometrical properties given by
Tables 1 and 2 contain numerical results for the controlled total energy and the control force of simply supported Al/ZrO_{2} FG plate with their counterparts and which are calculated neglecting the effect of the normal strains (). To study the influence of the normal strain effect on the control process, the following relative errors for the control force and the total energy are presented in Tables 1, 2, and 3 for various values of the ratios , and the index . In these tables,



Note that the relative errors in the control force decrease with the increase of ratio , and these errors increase with increasing the index of the volume fraction profile, so that the errors reach more than 22% for plates with and . This means that the effect of the normal strain is weak for thin FG plates with . Moreover, the inclusion of the effect of normal strain in the formulation with the cases is required for obtaining more accurate results. This may be explained by the fact that, with increasing the index , the homogeneity of the plate differs rapidly through the plate thickness. This requires the inclusion of the normal strain effect in the formulation to give more accurate description for the deformation process through the thickness. Table 3 contains similar numerical results as in Tables 1 and 2 for second mode , . These results show that the relative errors increase with the increase of the mode numbers and confirming the importance of the inclusion of the normal strain effect in the formulation of the problem to obtain more accurate results. In general, the effect of the normal strain on the accuracy of the numerical results of the total energy is less significant than that on the control force .
Figures 1 and 2 contain curves for the controlled force and the total energy () versus the ratios and for different values of index for SSSS plates. These curves show that the volume fraction profile through the plate thickness has a significant effect on the control force and controlled energy, particularly for FG plates with and square plates. Moreover, the curves show the high sensitivity of the total energy and control force of the FG plates versus the change of ratios , and the index . This indicates that these design variables may play an important role in reducing the vibrational response of the FG plates. Figure 3 shows curves and curves against the time parameter . These curves show that the present design control procedure reduces the expenditure of the control energy and the time of the controlling process. Figure 4 contains curves for the control force and the total energy against the index number for different values of the thickness ratio . These curves confirm the previous discussion on the efficiency of the present procedure in damping process, where the use of the thickness ratio as design variable reduces significantly the level of the total energy and decreases the consumed control energy. Figure 5 contains curves and curves against the volume fraction . Note that the design procedure over the plate thickness is more efficient than that over the volume fractional profile, while the design procedure over both and is the most efficient. In general the present control and design procedure strongly improves the performance of the FG plates.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
7. Conclusion
A structural and control optimization procedure is presented to minimize the vibrational response of FG plates. A firstorder plate theory including the normal strain and shear deformation effects is used. Optimal level of closedloop control force, plate thickness, and the volume fractional index are determined. The influence of the normal strain effect on the accuracy of the obtained results is illustrated. It is found that the inclusion of the normal stain effect is required, so that neglecting the normal stain effect leads to errors reaching more than 22% for square plates with homogeneity index . The presented results show that the present control and design procedure has a high efficiency in improving the performance of the FG plates.
Appendices
A. The Definition of the Parameters in (6a) and (6b)
The stiffness coefficients , , and in (6a) and (6b) are
B. The Definition of the Parameters in (12) and (14)
The coefficients , , , , , , , in (12) are The coefficients () in (14) are