The numerical analysis of static deformations and free vibrations of shear flexible plates often uses the differential quadrature method [1]. In addition, the mechanical properties of isotropic and laminated plates are determined by the collocation method of an asymmetric radial basis function [2]. Not only a laminated composite plate, but also piezo-laminated samples and laminate samples with distributed actuators can serve as working models for this type of analysis. The advanced computational models presented in this paper have the main objective of exploiting the potential associated with the use of wavelet bases and multiscale methods in the analysis of composite plate problems [3, 4]. The use of such bases leads, on the one hand, to the development of methods of arbitrarily high order, and on the other hand, they exhibit local behavior. Therefore, they can take into account the presence of singularities by locally refining the approximation space. Moreover, they provide a singularity detection tool that has been successfully used in developing efficient adaptive schemes for the numerical solution of several classes of differential equations.
The method used in this paper for the numerical solution is a collocation method based on the Deslauriers-Dubuc interpolation basis in a hierarchical form, which is the first necessary step to apply the adaptive wavelet collocation method to this class of problems. This collocation algorithm can be considered as a very efficient meshless method. The numerical model used a collocation method based on interpolating Delaurier-Dubuc wavelets. Uniform discretization was considered for the model. However, it should be noted that although the collocation method does not require a uniform grid a priori, it can be easily adapted to the case of non-uniform grids of dyadic points. The theory of first-order shear deformation allows to introduce both transverse displacement and angular displacements for rotation around perpendicular axes for the description of vibration characteristics of laminated composites. Transverse and angular displacements are independently interpolated due to the break of the connection between in-plane displacements and bending displacements for symmetrically laminated plates. The equations can be easily obtained from the following. The calculation model is based on the equations of motion for free vibrations of laminated plates, which are modified in the next stage into equations for isotropic plates of laminated composites.
The presence of oscillatory processes in the volume of the laminated composite was investigated under the following boundary conditions: an arbitrary edge with simply supported, clamped or free. The governing equations were written for normal and tangential directions of loading at the edge. The following were subsequently determined: normal bending moment, torque, shear force at the edge of the plate; rotation about axes parallel to the tangential and normal directions at the edge of the plate. At the layer boundaries, continuity of transverse shear stress is required for laminates with different materials along the thickness direction. According to the first-order shear strain assumptions, the transverse shear strain is constant across the thickness, which is a rough approximation to the actual variation even for a uniform cross-section. For uniform cross-sections, the shear strain is usually assumed to be a parabolic function of the normal coordinate. Therefore, a shear correction factor was used to approximate the transverse strain energy on an average basis.
In the first computational study, a square isotropic composite plate with simple supports was considered, which were subjected to a uniform force load. The central node of the plate was characterized by a set of finite number of normalized displacements and normal mechanical stresses. For the second numerical study, a square plate with simple supports was chosen. The plate consisted of a composite material and included four identical layers oriented at an angle of [00/900/900/00]. The plate was under a sinusoidal vertical pressure. In-plane displacements, transverse displacements, normal stresses, and in-plane transverse shear stresses were presented in normalized form.
The improved numerical method produces close to exact displacements for thinner laminated composite plates. The results for thick plates differ from the exact solution due to the absence of normal transverse strain. It should be noted that the obtained solutions are in good agreement with the finite element solutions. The method also produces highly accurate normal stresses. It should be noted that the presented study of static deformations and free vibrations of shear flexible isotropic and layered composite plates corresponds to the results of the first-order shear deformation theory. The numerical method developed in this work allowed us to analyze vibrations in local volumes of isotropic laminated composite plates, as well as in sandwich plates. The numerical experiments conducted allow us to conclude that this method can be used to describe static deformations and free vibrations in composite and sandwich plates.
References:
1. Amand Robinson M.T.., and Adali S. Dynamic stability of viscoelastic plates under axial flow by differential quadrature method. Engineering Computations. 2017. Vol. 34(4). Pp. 1240-21256. https://doi.org/ 10.1108/EC-03-2016-0113
2. Rodrigues J.D., Roque C.M.C., and Ferreira A.J.M Analysis of isotropic and laminated plates by an affine space decomposition for asymmetric radial basis functions collocation. Engineering analysis with boundary elements. 2012. Vol. 36(5). Pp. 709-715.
https://doi.org/10.1016/j.enganabound/2011.11.021
3. Zuo H. et al. Analysis of laminated composite plates using wavelet finite element method for higher-order plate theory. Composite Structures. 2015. Vol. 131. Pp. 248-258. https://doi.org/10.1016/j.compstruct.201504.064
4. Zhang X. et al. Analysis of laminated plates and shells using B-spline wavelet on interval finite element. International Journal of Structural Stability and Dynamics. 2017. Vol. 17(06). P. 1750062. https://doi.org/10.1142/S0219455417500626
|