• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.6 no.6
• /
• pp.525-530
• /
• 2005

In this study, it is presented analysis results of bending problems in the anisotropic thick shell and the anisotropic thin shell bending problems. In the numerical analysis of various mechanical problems involving complex partial differential equations, finite element method is used. Both Kirchoffs assumptions and Mindlin assumptions are used as the basic governing equations of bending problems in the anisotropic shells. The analysis results are compared between the anisotropic thick shells and the anisotropic thin shells for the various width-thickness ratios. The numerical method of this study will be contributed not only to analysis the bending behavior of anisotropic shells but also to design the anisotropic shells.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.35 no.1
• /
• pp.54-59
• /
• 1999

High-speed electronic digital computers have enabled engineers to employ various numerical discretization techniques for solutions of complex problems. The Finite Element Method is one of the such technique. The Finite Element Method is one of the numerical analysis based on the concepts of fundamental mathematical approximation. Three dimensional plate structures used often in partition of ship, box girder and frame are analyzed by Finite Element Method. In design of structures, the static deflections, stress concentrations and dynamic deflections must be considered. However, these problem belong to geometrically nonlinear mechanical structure analysis. The analysis of each element is independent, but coupling occurs in assembly process of elements. So, to overcome such a difficulty the shell theory which includes transformation matrix and a fictitious rotational stiffness is taken into account. Also, the Mindlin's theory which is considered the effect of shear deformation is used. The Mindlin's theory is based on assumption that the normal to the midsurface before deformation is "not necessarily normal to the midsurface after deformation", and is more powerful than Kirchoff's theory in thick plate analysis. To ensure that a small number of element can represent a relatively complex form of the type which is liable to occur in real, rather than in academic problem, eight-node quadratic isoparametric elements are used. are used.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2008.11a
• /
• pp.71-81
• /
• 2008

This paper proposes a spectral element which can represent dynamic responses in high frequency domain such as Lamb waves on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by piezoelectric layer (PZT layer) bonded on a base plate. In the two layer beam model, a PZT layer is assumed to be rigidly bonded on a base beam. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with electro mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are formulated through equations of motions converted into frequency domain. A detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through comparison results with the conventional 2-D FEM and the previously developed spectral elements.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.11
• /
• pp.1157-1169
• /
• 2008

This paper presents spectral element formulation which approximates Lamb wave propagation by PZT transducers bonded on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by a piezoelectric (PZT) layer rigidly bonded on a base plate. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Euler-Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with the electro-mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are obtained through equations of motions converted into frequency domain. Detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through numerical examples.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.2
• /
• pp.63-70
• /
• 2014

Design sensitivity analysis and topology design optimization for a piezoelectric crystal resonator are developed. The piezoelectric crystal resonator is deformed mechanically when subjected to electric charge on the electrodes, or vice versa. The Mindlin plate theory with higher-order interpolations along thickness direction is employed for analyzing the thickness-shear vibrations of the crystal resonator. Thin electrode plates are masked on the top and bottom layers of the crystal plate in order to enforce to vibrate it or detect electric signals. Although the electrode is very thin, its weight and shape could change the performance of the resonators. Thus, the design variables are the bulk material densities corresponding to the mass of masking electrode plates. An optimization problem is formulated to find the optimal topology of electrodes, maximizing the thickness-shear contribution of strain energy at the desired motion and restricting the allowable volume and area of masking plates. The necessary design gradients for the thickness-shear frequency(eigenvalue) and the corresponding mode shape(eigenvector) are computed very efficiently and accurately using the analytical design sensitivity analysis method using the eigenvector expansion concept. Through some demonstrative numerical examples, the design sensitivity analysis method is verified to be very efficient and accurate by comparing with the finite difference method. It is also observed that the optimal electrode design yields an improved mode shape and thickness-shear energy.