• Structural Engineering and Mechanics
• /
• v.4 no.3
• /
• pp.209-226
• /
• 1996

One of the major divisions in the mathematical modelling of a tubular structure is to include the effect of the transverse shear stress and rotary inertia in vibration of members. During the past three decades, problems of vibration of tubular structures have been considered by some authors, and special attention has been devoted to the Timoshenko theory. There have been considerable efforts, also, to apply the method of spectral analysis to vibration of a structure with rectangular section beams. The purpose of this paper is to compare the results of the spectrally formulated finite element analyses for the Timoshenko theory with those derived from the conventional finite element method for a tubular structure. The spectrally formulated finite element starts at the same starting point as the conventional finite element formulation. However, it works in the frequency domain. Using a computer program, the proposed formulation has been extended to derive the dynamic response of a tubular structure under an impact load.

• Proceedings of the KSME Conference
• /
• 2001.06c
• /
• pp.402-407
• /
• 2001

Sensitivity analysis scheme is developed in the elasto-plastic finite element method with explicit time integration using direct differentiation method. The direct differentiation is concerned with the time integration, constitutive relation, shell element with reduced integration and the contact scheme. Sensitivity analysis results are mainly examined with the highly nonlinear and quasi-static problem with the complicated contact condition. The result shows stable sensitivity especially in the sheet metal forming analysis.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.8 no.3
• /
• pp.151-162
• /
• 2000

This study is concerned with computerized multi-level substructuring methods and stress analysis of coil springs. The purpose of substructuring methods is to reduce computing time and capacity of computer memory by multiple level reduction of the degrees of freedom in large size problems that are modeled by three dimensional continuum finite elements. In this paper, the spring super element developed is investigated with tension, torsion, and bending of a cylindrical bar in order to verify its accuracy and efficiency for the multi-level substructuring method. And then the algorithm is applied to finite element analysis of coil springs. The result demonstrates the validity of the multi-level substructuring method and the efficiency in computing time and memory by providing good computational results in coil spring analysis.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.859-865
• /
• 2007

This paper presents a method to analyze the unbalance response of a high speed polygon mirror scanner motor supported by sintered bearing and flexible supporting structures by using the finite element method and the mode superposition method. The appropriate finite element equations for polygon mirror are described by rotating annular sector element using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. The rotating components except for the polygon mirror are modeled by Timoshenko beam element including the gyroscopic effect. The flexible supporting structures are modeled by using a 4-node tetrahedron element and 4-node shell element with rotational degrees of freedom. Finite element equations of each component of the polygon mirror scanner motor and the flexible supporting structures are consistently derived by satisfying the geometric compatibility in the internal boundary between each component. The rigid link constraints are also imposed at the interface area between sleeve and sintered bearing to describe the physical motion at this interface. A global matrix equation obtained by assembling the finite element equations of each substructure is transformed to a state-space matrix-vector equation, and both damped natural frequencies and modal damping ratios are calculated by solving the associated eigenvalue problem by using the restarted Arnoldi iteration method. Unbalance responses in time and frequency domain are performed by superposing the eigenvalues and eigenvectors from the free vibration analysis. The validity of the proposed method is verified by comparing the simulated unbalance response with the experimental results. This research also shows that the flexibility of supporting structures plays an important role in determining the unbalance response of the polygon mirror scanner motor.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.20-25
• /
• 2008

The finite element analysis for porous media is severe job because constituents have different physical peoperties, and element's continuity and stability should be considered. Thus, we propose the new mixed finite element method in order to overcome the problems. In this method, multi time step, remeshing step, and sub iteration step are introduced. The multi time step and remeshing step make it possible to satisfy a stability and an accuracy during sub iteration in which global time is determined. Finally, the proposed method is compared with the ABAQUS(2007) software and exact solution(Schiffman 1967) through two dimensional consolidation model.

• Earthquakes and Structures
• /
• v.9 no.1
• /
• pp.173-194
• /
• 2015

Site response analysis is an important topic in earthquake engineering. A time-domain numerical method called as one-dimensional (1D) finite element artificial boundary method is proposed to simulate the homogeneous plane elastic wave propagation in a layered half space subjected to the obliquely incident plane body wave. In this method, an exact artificial boundary condition combining the absorbing boundary condition with the inputting boundary condition is developed to model the wave absorption and input effects of the truncated half space under layer system. The spatially two-dimensional (2D) problem consisting of the layer system with the artificial boundary condition is transformed equivalently into a 1D one along the vertical direction according to Snell's law. The resulting 1D problem is solved by the finite element method with a new explicit time integration algorithm. The 1D finite element artificial boundary method is verified by analyzing two engineering sites in time domain and by comparing with the frequency-domain transfer matrix method with fast Fourier transform.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.366-371
• /
• 2007

In many CAD systems, NURBS has been employed to construct exact geometries. Recently, NURBS finite element analysis methods were proposed by some authors for convenient connection between CAD and finite element analysis. Additional advantages of NURBS FEA, such as exact geometry and no mesh generation, are obtained. However, NURBS is inefficient in local refinement and merging patches. For refinement of local region in interest, additional control points should be inserted into the entire row or column which contains the local region. There is another inefficiency of NURBS during merging patches into a large structure due to propagation of control points. In order to overcome these inefficiencies of NURBS, T-spline was proposed by Sederberg. In this work, T-spline based finite element method is proposed for efficient local refinement and merging patches. At first, accuracy and efficiency of NURBS FEA is verified and efficiency of T-spline FEA is verified by comparing with NURBS FEA.

• Biomaterials and Biomechanics in Bioengineering
• /
• v.3 no.2
• /
• pp.71-84
• /
• 2016

The finite element method is wide used in simulation in the biomechanical structures, but a lack of studies concerning finite element mesh quality in biomechanics is a reality. The present study intends to analyze the importance of the mesh quality in the finite element model results from humeral structure. A sensitivity analysis of finite element models (FEM) is presented for the humeral bone and cartilage structures. The geometry of bone and cartilage was acquired from CT scan and geometry reconstructed. The study includes 54 models from same bone geometry, with different mesh densities, constructed with tetrahedral linear elements. A finite element simulation representing the glenohumeral-joint reaction force applied on the humerus during $90^{\circ}$ abduction, with external load as the critical condition. Results from the finite element models suggest a mesh with 1.5 mm, 0.8 mm and 0.6 mm as suitable mesh sizes for cortical bone, trabecular bone and humeral cartilage, respectively. Relatively to the higher minimum principal strains are located at the proximal humerus diaphysis, and its highest value is found at the trabecular bone neck. The present study indicates the minimum mesh size in the finite element analyses in humeral structure. The cortical and trabecular bone, as well as cartilage, may not be correctly represented by meshes of the same size. The strain results presented the critical regions during the $90^{\circ}$ abduction.

• Tribology and Lubricants
• /
• v.8 no.1
• /
• pp.38-43
• /
• 1992

In solving the Reynolds equation in dynamically loaded bearing problem, it is almost impossible to find the squeeze velocity and the cavitation region by analytical method. Finite Element Method was applied to the piston-ring lubrication analysis to solve the complementary problem. The method was very efficient and any convergence problem was not encountered.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.624-634
• /
• 1994

A numerical method which combines equal-order velocity-pressure formulation originated from SIMPLE algorithm and streamline upwinding method has been developed. To verify the proposed numerical method, we considered the lid-driven cavity flow and backward facing step flow. The trend of convergence history is stable up to the error criterion beyond which the maximum value of error is oscillatory due4 to the round-off error. In the present study, all results were obtained with the single precision calculation up to the given error criterion and it was found to be sufficient for our purpose. The present results were then compared with existing experimental results using laser doppler velocimetry and numerical results using finite difference method and mixed interpolation finite element method. It has been shown that the present method gives accurate results with less memories and execution time than the coventional finite element method.