• Proceedings of the KSR Conference
• /
• 2005.05a
• /
• pp.344-349
• /
• 2005

The axially moving thin beam-plates exposed to sudden thermal loadings may experience severe vibrations through the thermal shock process. For accurate prediction of the thermal shock-induced vibrations, this paper develops a spectral element model for axially moving thermoelastic beam-plates. The spectral element model which is represented by spectral element matrix is formulated from the frequency-dependent dynamic shape functions which satisfy the governing equations in the frequency-domain. Thus, when compared with the classical finite element model in which simple polynomial functions are used as the shape functions, the spectral element model can provide exact solution by treating a whole uniform structure member as a single finite element, regardless of its length.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.16 no.4 s.109
• /
• pp.387-393
• /
• 2006

In this paper, a spectral element model is developed for the uniform straight pipelines conveying internal unsteady fluid flow. The spectral element matrix is formulated by using the exact frequency-domain solutions of the pipe-dynamics equations. The spectral element dynamic analysis is then conducted to evaluate the accuracy of the present spectral element model and to investigate the vibration characteristics and internal fluid characteristics of an example pipeline system.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.9
• /
• pp.137-147
• /
• 1995

The unbalance response analysis is one of the essential area in the forced vibration analysis of rotor-bearing systems. Local bearing parameters in rotor-bearing systems are the major sources which give rise to a difficulty in unbalance response computation due to the complicated dynamic properties such as rotational speed dependency and anisotropy. In the present paper, an efficient method for unbalance responses is proposed so as to easily take into account bearing parameters in computation. An exact matrix condensation procedure is proposed which enables the present method to compute unbalance responses by dealing with condensed, small matrices. The proposed method causes no errors even though the computation procedure is based on the small matrices condensed from the full matrices. The present method is illustrated through a numerical example and compared with the conventional method.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1998.04a
• /
• pp.200-205
• /
• 1998

Although discretization methods such as the transfer matrix method (TMM) and the finite element method (FEM) have played an important role in the design or analysis of rotor-bearing systems, continuous system modeling and analysis are often desirable especially for sensitivity analysis or design. The present paper proposes a comprehensive modeling procedure to obtain exact solution of general rotor-bearing systems. The proposed method considers a Timoshenko beam model and makes use of complex coordinate in the formulation. The proposed method provides exact eigensolutions and frequency response functions (FRFS) of general multi-stepped rotor-bearing systems. The first numerical example compares the proposed method with FEM. The numerical study proves that the proposed method is very efficient and useful for the analysis of rotor-bearing systems.

• Structural Engineering and Mechanics
• /
• v.22 no.6
• /
• pp.701-717
• /
• 2006

In the existing reports regarding free transverse vibrations of the Euler-Bernoulli beams, most of them studied a uniform beam carrying various concentrated elements (such as point masses, rotary inertias, linear springs, rotational springs, spring-mass systems, ${\ldots}$, etc.) or a stepped beam with one to three step changes in cross-sections but without any attachments. The purpose of this paper is to utilize the numerical assembly method (NAM) to determine the exact natural frequencies and mode shapes of the multiple-step Euler-Bernoulli beams carrying a number of lumped masses and rotary inertias. First, the coefficient matrices for an intermediate lumped mass (and rotary inertia), left-end support and right-end support of a multiple-step beam are derived. Next, the overall coefficient matrix for the whole vibrating system is obtained using the numerical assembly technique of the conventional finite element method (FEM). Finally, the exact natural frequencies and the associated mode shapes of the vibrating system are determined by equating the determinant of the last overall coefficient matrix to zero and substituting the corresponding values of integration constants into the associated eigenfunctions, respectively. The effects of distribution of lumped masses and rotary inertias on the dynamic characteristics of the multiple-step beam are also studied.

• Structural Engineering and Mechanics
• /
• v.39 no.1
• /
• pp.21-46
• /
• 2011

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.

• 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.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.18 no.2
• /
• pp.58-64
• /
• 2019

This paper describes the analysis of dynamic characteristics and prediction of the stiffness for the joint between structural members. In the process of deriving the governing equations, the stiffness values responsible for the moment and shear force were modelled by using linear and torsional springs in the middle of a clamped-clamped beam. The sensitivities of the natural frequency and modal assurance criterion were investigated as a function of the dimensionless linear and torsional spring stiffness. The reliability of the predictions for the linear and torsional stiffness values was verified by the inverse computations of the stiffness matrix. The predictive and exact theoretical stiffness values were compared for the stiffness element in the finite element formulation, and their results show an excellent correlation. It is strongly anticipated that although the proposed methodology is currently limited to the analytical utilization, it will provide a useful tool to estimate unknown joint stiffness values based on the experimental natural frequency and mode shape.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.9 no.1
• /
• pp.271-281
• /
• 2005

The effects of radiation damping is used to compensate the truncated boundary which is relatively close to the structure-fluid interface in the fluid element surrounding the submerged structures. An efficient ring element is presented to model the shell and fluid element which fully utilizes the characteristics of the axisymmetry. The computational model uses the technique which separate the meridional shape and circumferential wave mode and gets similar result with the exact solution in the eigenvalues and the earthquake analysis. The fluid-structure interaction techniques is developed in the finite element analysis of two dimensional problems using the relations between pressure, nodal unknown acceleration and added mass assuming the fluid to be invicid, incompressible and irrotational. The effectiveness and efficiency of the technique is demonstrated by analyzing the free vibration and seismic analysis using the added mass matrix considering the structural deformation effect.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.2
• /
• pp.87-96
• /
• 1999

To analyze the dynamic behavior of structure, direct integration and mode superposition may be utilized in time domain analysis. As finite number of frequencies can give relatively exact solutions, mode superposition is preferable in analyzing structural behavior. In non-linear analysis, however, mode superposition is seldom used since time-varying element stiffness changes stiffness matrix, and the change of stiffness matrix leads to the change of essential constants - natural frequencies and mode shapes. In spite of these difficulties, there are some attempts to adopt mode superposition because of low cost compared to direct integration, but the result is not satisfactory. In this paper, a method using mode superposition in non-linear analysis is presented by separating local element stiffness from global stiffness matrix with the difference between linear and non-linear restoring forces to the external force vectors included. Moreover, the hysteresis model changing with the relative deformation in each floor makes it possible to analyze non-linear behavior of structure. The proposed algorithm is applied to shear beam model and the maximum displacement is compared with the result using direct integration method.