• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.291-298
• /
• 1994

This paper presents a new stiffness matrix for the analysis of arbitrary thin-walled open beams in warp-restrained torsion. The element accounts for both flexural and warping torsional effects. To eliminate the ad hoc introduction of St. Venant stiffness in this $C^0$ element, the virtual work equation based on an orthogonal Cartesian coordinate system is used. The effectiveness of the derived block stiffness is addressed. The transformation matrix between two different reference systems is also shown. Numerical examples using the proposed matrix are compared with the classical solutions or other previous results in the literature.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.15 no.1
• /
• pp.77-84
• /
• 2011

This study is to develop nonlinear analysis algorithm for transfer matrix method, which can be applied to continuous beam analysis. Gauss-Lobatto integral rule is adopted and the transfer matrix is derived from stiffness matrix. In the transfer matrix method, the system equation has a constant number of unknowns regardless of number of D.O.F. Therefore, the transfer matrix method has computational efficiencies not only in linear elastic analysis but also in nonlinear analysis. To verify the developed method, the analysis results of several examples are compared with commercial code in moment-curvature, moment-displacement and load-displacement relation.

• Earthquakes and Structures
• /
• v.3 no.3_4
• /
• pp.383-411
• /
• 2012

The objective of the present paper is to review and implement the most recent developments in the Spectral Element Method (SEM), as well as improve aspects of its implementation in the study of wave propagation by numerical simulation in elastic unbounded domains. The classical formulation of the method is reviewed, and the construction of the mass matrix, stiffness matrix and the external force vector is expressed in terms of matrix operations that are familiar to earthquake engineers. To account for the radiation condition at the external boundaries of the domain, a new absorbing boundary condition, based on the Perfectly Matched Layer (PML) is proposed and implemented. The new formulation, referred to as the Multi-Axial Perfectly Matched Layer (M-PML), results from generalizing the classical Perfectly Matched Layer to a medium in which damping profiles are specified in more than one direction.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.05a
• /
• pp.660-665
• /
• 2006

This paper presents a finite element method to analyze the free and forced vibration of a hard disk drive (HDD) considering the flexibility of a spinning disk-spindle with fluid dynamic bearings (FDBs), an actuator with pivot bearings, an air bearing between head-disk interface and the base with complicated geometry. Finite element equation of each component is consistently derived with the satisfaction of the geometric compatibility of the internal boundary between each component. The spinning disk, hub and FDBs are modeled by annular sector elements, beam elements and stiffness and damping elements, respectively. The actuator am, E-block, suspension and base plate are modeled by tetrahedral elements. The pivot bearing in the actuator and the air bearing between head-disk interfaces are modeled by the stiffness element with five degrees of freedom and the axial stiffness, respectively. 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 with the restarted Arnoldi iteration method. Modal and shock testing are performed to show that the proposed method well predicts the vibration characteristics of a HDD.

• Journal of Advanced Marine Engineering and Technology
• /
• v.23 no.4
• /
• pp.504-513
• /
• 1999

Usually bellows are designed for the purpose of absorbing axial movement. To find out axial stiffness of bellows the axisymmetric shell theory using the finite element method is adopted in this paper. Bellows can be idealised by series of conical frustum-shaped elements because it is axisymmetric shell structure. The force required to deflect bellows axilly is a function of the dimensions of the bellows and the materials from which they are made. The displancements of nodal points due to small increment of force are calculated by the finite element method and the calculated nodal displacements are added to r-z cylinderical coordinates of nodal points. The new stiffness matrix of the system using the new coordinates of nodal points is adopted to calculate the another increments of nodal dis-placements that is the step by method is used in this paper. spring constant is analyzed according to the changing geometric factors of u-shaped bellows. The FEM results were agreed with experiment. Using developed FORTRAN PROGRAM spring constant can be predicted by input of a few factors.

• Proceeding of EDISON Challenge
• /
• 2016.03a
• /
• pp.195-202
• /
• 2016

In this paper, crack detection method using eigen value analysis is presented. Three methods are used: theoretical analysis, finite element method with the cracked beam elements and finite element method with three dimensional continuum elements. Finite element formulation of the cracked beam element is introduced. Additional term about stress intensity factor based on fracture mechanics theory is added to flexibility matrix of original beam to model the crack. As using calculated stiffness matrix of cracked beam element and mass matrix, natural frequencies are calculated by eigen value analysis. In the case of using continuum elements, the natural frequencies could be calculated by using EDISON CASAD solver. Several cases of crack are simulated to obtain natural frequencies corresponding the crack. The surface of natural frequency is plotted as changing with crack location and depth. Inverse analysis method is used to find crack location and depth from the natural frequencies of experimental data, which are referred by another papers. Predicted results are similar with the true crack location and depth.

• Proceedings of the Korea Concrete Institute Conference
• /
• 1998.10a
• /
• pp.450-455
• /
• 1998

In this paper, finite element with embedded discontinuous line is introduced in order to avoid the difficulties of adding new nodal points along with crack growth in discrete crack model. With the discontinuous element using discontinuous shape function, stiffness matrix of finite element is derived and dual mapping technique for numerical integration is employed. Using the finite element program made with employed algorithms, algorithm is verified and fracture analysis of simple concrete beam is performed.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.11b
• /
• pp.1066-1071
• /
• 2002

The use of frequency-dependent spectral element matrix (or exact dynamic stiffness matrix) in structural dynamics is known to provide very accurate solutions, while reducing the number of degrees-of-freedom to resolve the computational and cost problems. Thus, in the present paper, the spectral element model is formulated for the axially moving Timoshenko beam under a uniform axial tension. The high accuracy of the present spectral element is then verified by comparing its solutions with the conventional finite element solutions and exact analytical solutions. The effects of the moving speed and axial tension on the vibration characteristics, the dispersion relation, and the stability of a moving Timoshenko beam are investigated, analytically and numerically.

• Journal of KSNVE
• /
• v.10 no.1
• /
• pp.97-106
• /
• 2000

The dynamic behavior of crankshaft-bearing system in scroll compressor has been investigated using the combined methodologies of finite elements and transfer matrices. The finite element formulation is proposed including the field element for a shaft section and the point element at balancer weight locations, bearing locations, etc., whereas the conventional method is used with the elements. The Houbolt method is used to consider the time march for the integration of the system equations. The linear stiffness and damping coefficients are calculated for a finite cylindrical fluid-film bearing by solving the Reynolds equation, using finite difference method. The orbital response of crankshaft supported on the linear bearing model is obtained, considering balancer weights of motor rotor. And, the steady state displacement of crankshaft are compared with a variation in balancer weight. The loci of crankshaft at bearing locations are composed of the synchronous whirl component and the non-synchronous whirl component.

• Journal of the Korean Society for Railway
• /
• v.6 no.2
• /
• pp.129-134
• /
• 2003

This paper investigates the dynamic characteristics of a beam axially moving over multiple elastic supports. The spectral element matrix is derived first for the axially moving beam element and then it is used to formulate the spectral element matrix for the moving beam element with an interim elastic support. The moving speed dependance of the eigenvalues is numerically investigated by varying the applied axial tension and the stiffness of the elastic supports. Numerical results show that the fundamental eigenvalue vanishes first at the critical moving speed to generate the static instability.