• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.04a
• /
• pp.309-314
• /
• 2001

In this study, the nonlinear soil structure interaction analysis method based on finite element and boundary element method is developed. In the seismic region, the nonlinearity of near field soil has to be considered for more exact reflection of soil-structure interaction effect. Thus, nonlinear finite element program coupled with boundary elements is developed for nonlinear soil-structure interaction analysis. Using the developed numerical algorithm, the nonlinear soil-structure interaction analysis is performed and responses due to dynamic forces and seismic excitation are investigated. The developed method is verified by comparing with previous studies.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.2
• /
• pp.161-169
• /
• 2002

Based on the analytic expression for the elasto-dynamic behavior of rotating shaft, the transfer matrix is formulated for the shaft element with uniform cross-section. Timoshenko beam theory is Introduced for modeling the behavior of shaft. Complex variables representing the displacement, slope, moment and shear force are used for deriving the transfer matrix between both ends of the shaft element. Simulation result obtained by applying the transfer matrix to a general rotor model is compared with the reference result and proved to be exact. Natural frequencies and the corresponding modes are analyzed with varying the bearing: stiffness. The generally used bearings are considered for discussions. and the bearing stiffness is shown to affect the vibration characteristics of rotor.

• Structural Engineering and Mechanics
• /
• v.91 no.6
• /
• pp.591-606
• /
• 2024

This article proposes a new methodology for identifying beam damage based on changes in modal parameters using the Double Stage Extended Improved Particle Swarm Optimization (DSEIPSO) technique. A finite element code is first developed in MATLAB to model an ideal beam structure based on classical beam theory. An experimental study is then performed on a laboratory-scale beam, and the modal parameters are extracted. An improved version of the PSO algorithm is employed to update the finite element model based on the experimental measurements, representing the real structure and forming the baseline model for all further damage detection. Subsequently, structural damages are introduced in the experimental beam. The DSEIPSO algorithm is then utilized to optimize the objective function, formulated using the obtained mode shapes and the natural frequencies from the damaged and undamaged beams to identify the exact location and extent of the damage. Experimentally obtained resultsfrom a simple cantilever beam are used to validate the effectiveness of the proposed method. The illustrated results show the effectiveness of the proposed method for structural damage detection in the SHM field.

• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.11
• /
• pp.55-61
• /
• 1999

The torsional and axial vibrations of shaft system have been calculated independently because of both the limitation of computing time and the complexity of crankshaft model. In actual system, however, the torsional and axial vibrations are coupled. Therefore, in recent, many works in the coupled vibration analysis have been done to find out the more exact dynamic behavior of shaft system. The crankshaft model is very important in the vibration analysis of shaft system because most of excitation forces act on the crankshaft. It is, however, difficult to establish an exact model of crankshaft since its shape is very complex. In this work, an efficient method is proposed to construct the stiffness matrix of crankshaft using a finite element model of half crankthrow. The proposed and existing methods are compared by applying to both a simple thick beam with circular cross section and an actual crankshaft.

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

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.13 no.2
• /
• pp.10-15
• /
• 1984

The aims of this study are to obtain a suitable method and a proper mesh for investigation of the temperature distribution and heat transfer. The relative errors of temperature distribution and heat transfer for each mesh are acquired in accordance with linear finite element (FEM 3), square finite element (FEM 6), cubic finite element (FEM 10), and finite difference method (FDM). It has been found that FEM 10 is the most accurate measure to obtain the temperature distribution and heat transfer. However, no significant results have been obtained successfully, because when higher order finite element methods are used the more computational efforts are necessary due to the distribution of elements. The results of this study are as follows ; 1 . In case of a=b=L, meshes for less than $1\%$ relative errors (temperature distribution) acquired in various methods to exact solution are $2\times2,\;4\times4,\;8\times8\;and\;8\tiems8$ for each FEM 10, FEM 6, FEM 3 and FDM and a=L, b=1/2L are $10\times5$ for each FEM 3 and FDM. And the tendency of results acquired of heat transfer is similar to those above. 2 . In computational efforts (a=b=L), FEM 6 has taken 21 times. and FEM 10 154times FEM 3 and FDM and FEM 3 is the sane as FDM.

• Journal of the Korean Society for Precision Engineering
• /
• v.18 no.4
• /
• pp.64-71
• /
• 2001

In general, precise linear actuators using piezoelectric element are driven by friction force. Exact understanding of friction plays an important role in analysis and control of a motor. In this research, we designed a precise linear actuator using piezoelectric elements and observed its dynamic characteristics. By varying phase angle difference and amplitudes of the sinusoidal waves that are driving inputs, we can know that it is possible to control moving direction and distance of the slider. As preload is increased, its moving distance is decreased. And also, we have modeled a precise linear actuator using stick slip friction models such as classical, Karnopp. and reset integrator. Finally, by comparing the results of simulation and experiment, it was verified that the model is well designed.

• International Journal of Precision Engineering and Manufacturing
• /
• v.3 no.4
• /
• pp.64-71
• /
• 2002

Recently, mechanical systems such as a high-speed vehicles and railway trains moving on flexible beam structures have become a very important issue to consider. Using the general approach proposed in the first part of this paper, it is possible to predict motion of the constrained mechanical system and the elastic structure, with various kinds of foundation supporting conditions. Combined differential-algebraic equation of motion derived from both multibody dynamics theory and finite element method can be analyzed numerically using a generalized coordinate partitioning algorithm. To verify the validity of this approach, results from the simply supported elastic beam subjected to a moving load are compared with the exact solution from a reference. Finally, parametric study is conducted for a moving vehicle model on a simply supported 3-span bridge.

• Coupled systems mechanics
• /
• v.7 no.3
• /
• pp.353-371
• /
• 2018

The Wavenumber or more accurately Wavenumber-FD approach was initially introduced for two-dimensional dynamic analysis of concrete gravity dam-reservoir systems. The technique was formulated in the context of pure finite element programming in frequency domain. Later on, a variation of the method was proposed which was referred to as Wavenumber-TD approach suitable for time domain type of analysis. Recently, it is also shown that Wavenumber-FD approach may be applied for three-dimensional dynamic analysis of concrete arch dam-reservoir systems. In the present study, application of its variation (i.e., Wavenumber-TD approach) is investigated for three-dimensional problems. The method is initially described. Subsequently, the response of idealized Morrow Point arch dam-reservoir system is obtained by this method and its special cases (i.e., two other well-known absorbing conditions) for time harmonic excitation in stream direction. All results for various considered cases are compared against the exact response for models with different values of normalized reservoir length and reservoir base/sidewalls absorptive conditions.