• Structural Engineering and Mechanics
• /
• v.75 no.6
• /
• pp.747-769
• /
• 2020

This paper provides a semi-analytical approach to investigate the variations of 3D displacement components, electric potential, stresses, electric displacements and transverse vibration frequencies in laminated piezoelectric composite plates based on the scaled boundary finite element method (SBFEM) and the precise integration algorithm (PIA). The proposed approach can analyze the static and dynamic responses of multilayered piezoelectric plates with any number of laminae, various geometrical shapes, boundary conditions, thickness-to-length ratios and stacking sequences. Only a longitudinal surface of the plate is discretized into 2D elements, which helps to improve the computational efficiency. Comparing with plate theories and other numerical methods, only three displacement components and the electric potential are set as the basic unknown variables and can be represented analytically through the transverse direction. The whole derivation is built upon the three dimensional key equations of elasticity for the piezoelectric materials and no assumptions on the plate kinematics have been taken. By virtue of the equilibrium equations, the constitutive relations and the introduced set of scaled boundary coordinates, three-dimensional governing partial differential equations are converted into the second order ordinary differential matrix equation. Furthermore, aided by the introduced internal nodal force, a first order ordinary differential equation is obtained with its general solution in the form of a matrix exponent. To further improve the accuracy of the matrix exponent in the SBFEM, the PIA is employed to make sure any desired accuracy of the mechanical and electric variables. By virtue of the kinetic energy technique, the global mass matrix of the composite plates constituted by piezoelectric laminae is constructed for the first time based on the SBFEM. Finally, comparisons with the exact solutions and available results are made to confirm the accuracy and effectiveness of the developed methodology. What's more, the effect of boundary conditions, thickness-to-length ratios and stacking sequences of laminae on the distributions of natural frequencies, mechanical and electric fields in laminated piezoelectric composite plates is evaluated.

• Structural Engineering and Mechanics
• /
• v.42 no.4
• /
• pp.471-488
• /
• 2012

An efficient one dimensional finite element model has been presented for the dynamic analysis of composite laminated beams, using the efficient layerwise zigzag theory. To meet the convergence requirements for the weak integral formulation, cubic Hermite interpolation is used for the transverse displacement ($w_0$), and linear interpolation is used for the axial displacement ($u_0$) and shear rotation (${\psi}_0$). Each node of an element has four degrees of freedom. The expressions of variationally consistent inertia, stiffness matrices and the load vector are derived in closed form using exact integration. The formulation is validated by comparing the results with the 2D-FE results for composite symmetric and sandwich beams with various end conditions. The employed finite element model is free of shear locking. The present zigzag finite element results for natural frequencies, mode shapes of cantilever and clamped-clamped beams are obtained with a one-dimensional finite element codes developed in MATLAB. These 1D-FE results for cantilever and clamped beams are compared with the 2D-FE results obtained using ABAQUS to show the accuracy of the developed MATLAB code, for zigzag theory for these boundary conditions. This comparison establishes the accuracy of zigzag finite element analysis for dynamic response under given boundary conditions.

• Structural Engineering and Mechanics
• /
• v.9 no.5
• /
• pp.449-467
• /
• 2000

The natural frequencies and the corresponding mode shapes of a non-uniform beam carrying multiple point masses are determined by using the analytical-and-numerical-combined method. To confirm the reliability of the last approach, all the presented results are compared with those obtained from the existing literature or the conventional finite element method and close agreement is achieved. For a "uniform" beam, the natural frequencies and mode shapes of the "clamped-hinged" beam are exactly equal to those of the "hinged-clamped" beam so that one eigenvalue equation is available for two boundary conditions, but this is not true for a "non-uniform" beam. To improve this drawback, a simple transformation function ${\varphi}({\xi})=(e+{\xi}{\alpha})^2$ is presented. Where ${\xi}=x/L$ is the ratio of the axial coordinate x to the beam length L, ${\alpha}$ is a taper constant for the non-uniform beam, e=1.0 for "positive" taper and e=1.0+$|{\alpha}|$ for "negative" taper (where $|{\alpha}|$ is the absolute value of ${\alpha}$). Based on the last function, the eigenvalue equation for a non-uniform beam with "positive" taper (with increasingly varying stiffness) is also available for that with "negative" taper (with decreasingly varying stiffness) so that half of the effort may be saved. For the purpose of comparison, the eigenvalue equations for a positively-tapered beam with five types of boundary conditions are derived. Besides, a general expression for the "normal" mode shapes of the non-uniform beam is also presented.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.4
• /
• pp.734-742
• /
• 2016

Global warming and the depletion of fossil fuels have been caused by decades of reckless development. Wind energy is one form of renewable energy and is considered a future energy source. The wind tower is designed with a fundamental frequency in the soft-stiff design between the 1P and 3P range to avoid resonance. Usually, to perform natural frequency analysis of a wind tower, the boundary condition is set to the Fixed-End, and soil-pile interaction is not considered. In this study, consideration of the effect of soil-pile interaction on the wind tower was included and the difference in the natural frequency was studied. The fixed boundary condition was not affected by the soil condition and depth of the pile and the coupled spring boundary condition was unaffected by the depth of pile but affected by the depth of the pile, and the Winkler spring boundary condition is affected by both the soil condition and the depth of the pile. Therefore, the coupled spring boundary condition should be used in shallow depth soil conditions because the soil condition does not take the shallow depth soil into consideration.

• Journal of Mechanical Science and Technology
• /
• v.14 no.2
• /
• pp.197-206
• /
• 2000

A boundary integral equation method in the shape design sensitivity analysis is developed for the elasticity problems with axisymmetric non-homogeneous bodies. Functionals involving displacements and tractions at the zonal interface are considered. Sensitivity formula in terms of the interface shape variation is then derived by taking derivative of the boundary integral identity. Adjoint problem is defined such that displacement and traction discontinuity is imposed at the interface. Analytic example for a compound cylinder is taken to show the validity of the derived sensitivity formula. In the numerical implementation, solutions at the interface for the primal and adjoint system are used for the sensitivity. While the BEM is a natural tool for the solution, more generalization should be made since it should handle the jump conditions at the interface. Accuracy of the sensitivity is evaluated numerically by the same compound cylinder problem. The endosseous implant-bone interface problem is considered next as a practical application, in which the stress value is of great importance for successful osseointegration at the interface. As a preliminary step, a simple model with tapered cylinder is considered in this paper. Numerical accuracy is shown to be excellent which promises that the method can be used as an efficient and reliable tool in the optimization procedure for the implant design. Though only the axisymmetric problem is considered here, the method can be applied to general elasticity problems having interface.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.18 no.1
• /
• pp.1-8
• /
• 1981

To contribute to the preventive measures against local vibrations of ship's deck panels, some investigations into the prediction method of the natural frequency of the vibration of stiffened plates were done. Firstly, an analytical method based on the orthotropic plate analogy and the Rayleigh method using eigenfunctions of the Euler beam was shown, and numerical results of a regularly stiffened plate were compared with experimental results. And then, the method was extended to stiffened plates having one or two irregular stiffeners to obtain an approximate formula showing the relation between the change of the natural frequency and the size of the irregular stiffeners. The latter case was investigated for the purpose of providing a convenient design manual applicable to cure of local resonant vibrations of ships' deck panels by additional reinforcement of one or two stiffeners. In the analytical development the boundary was assumed to be rigidly supported and elastically restrained against rotation. In the experiment, however, only an extreme case i.e. simply supported boundary was investigated. The results of the investigation show that there is a fairly good conformity between the analytical results and the experimental ones in the first case, and that the approximate formula for the second case is confirmed also to be reliable for the design purpose. Considering that actual boundary conditions of deck panels in ship structures lie mostly somewhere between the simple support and the fixed, the authors discussed problems of the joint efficiency at the boundary of deck panels from the viewpoint of the practical application of the formulae.

• Journal of Power System Engineering
• /
• v.12 no.2
• /
• pp.35-40
• /
• 2008

This paper deals with the free vibrations of conical shells with linearly variable thickness by the transfer influence coefficient method. The classical thin shell theory based upon the Flugge theory is assumed and the governing equations of a conical shell are written as a coupled set of first order matrix differential equations using the transfer matrix. The Runge-Kutta-Gill integration method is used to solve the governing differential equation. The natural frequencies and corresponding mode shapes are calculated numerically for the conical shells with linearly variable thickness and various boundary conditions at the edges. The present method is applied to conical shells with linearly varying thickness, and the effects of the semi-vertex angle, the number of circumferential waves and thickness ratio on vibration are studied.

• Smart Structures and Systems
• /
• v.20 no.3
• /
• pp.351-368
• /
• 2017

To peruse the free vibration of curved functionally graded piezoelectric (FGP) nanosize beam in thermal environment, nonlocal elasticity theory is applied for modeling the nano scale effect. The governing equations are obtained via the energy method. Analytically Navier solution is employed to solve the governing equations for simply supported boundary conditions. Solving these equations enables us to estimate the natural frequency for curved FGP nanobeam under the effect of a uniform temperature change and external electric voltage. The results determined are verified by comparing the results by available ones in literature. The effects of various parameters such as nonlocality, uniform temperature changes, external electric voltage, gradient index, opening angle and aspect ratio of curved FGP nanobeam on the natural frequency are successfully discussed. The results revealed that the natural frequency of curved FGP nanobeam is significantly influenced by these effects.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.13 no.3
• /
• pp.155-163
• /
• 2003

This paper presents the characteristics of the radial nitration of cylindrical piezoelectric transducers. The differential equations of piezoelectric radial motion have been derived in terms of the radial displacement and electric potential, which are functions of the radial and axial coordinates. Applying mechanical and electrical boundary conditions has yielded the characteristic equation of radial vibration. Numerical results of the natural frequencies have been compared with the experimental observations reported earlier for the transducers of several sizes, and have shown a good agreement for the fundamental mode. The paper discusses the dependence of the natural frequencies on the radius and thickness of the piezoelectric cylinders and the difference between Piezoelectric and elastic resonances.

• Steel and Composite Structures
• /
• v.28 no.4
• /
• pp.461-470
• /
• 2018

In this study, the effect of polyurethane foam filler, in addition to surface layer thickness and core material thickness, on vibration characteristics of sandwich structures was investigated. The manufacturing process was carried out according to the Taguchi method. The natural frequencies and damping ratios of the produced samples were determined experimentally for fixed-free boundary conditions. In addition, solid models were developed for test samples and their finite element analyses were performed with $ANSYS^{(R)}$ to obtain their natural frequencies and mode shapes. An acceptably good agreement was found with the comparison of experimental results with the numerically obtained ones. The most effective parameters on the vibration characteristics of the sandwich structure were determined by the Taguchi method.