• Journal of Korean Society of Steel Construction
• /
• v.14 no.4
• /
• pp.509-518
• /
• 2002

In order to perform the spatial buckling and static analysis of the nonsymmetric thin-walled beam-column element, improved exact static stiffness matrices were evaluated using equilibrium equation and force-deformation relationships. This numerical technique was obtained using a generalized linear eigenvalue problem, by introducing 14 displacement parameters and system of linear algebraic equations with complex matrices. Unlike the evaluation of dynamic stiffness matrices, some zero eigenvalues were included. Thus, displacement parameters related to these zero eigenvalues were assumed as polynomials, with their exact distributions determined using the identity condition. The exact displacement functions corresponding to three loadingcases for initial stress-resultants were then derived, by consistently combining zero and nonzero eigenvalues and corresponding eigenvectors. Finally, exact static stiffness matrices were determined by applying member force-displacement relationships to these displacement functions. The buckling loads and displacement of thin-walled beam were evaluated and compared with analytic solutions and results using ABAQUS' shell element or straight beam element.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.3007-3019
• /
• 1993

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• Advances in aircraft and spacecraft science
• /
• v.1 no.3
• /
• pp.291-310
• /
• 2014

An advanced model for the linear flutter analysis is introduced in this paper. Higher-order beam structural models are developed by using the Carrera Unified Formulation, which allows for the straightforward implementation of arbitrarily rich displacement fields without the need of a-priori kinematic assumptions. The strong form of the principle of virtual displacements is used to obtain the equations of motion and the natural boundary conditions for beams in free vibration. An exact dynamic stiffness matrix is then developed by relating the amplitudes of harmonically varying loads to those of the responses. The resulting dynamic stiffness matrix is used with particular reference to the Wittrick-Williams algorithm to carry out free vibration analyses. According to the doublet lattice method, the natural mode shapes are subsequently used as generalized motions for the generation of the unsteady aerodynamic generalized forces. Finally, the g-method is used to conduct flutter analyses of both isotropic and laminated composite lifting surfaces. The obtained results perfectly match those from 1D and 2D finite elements and those from experimental analyses. It can be stated that refined beam models are compulsory to deal with the flutter analysis of wing models whereas classical and lower-order models (up to the second-order) are not able to detect those flutter conditions that are characterized by bending-torsion couplings.

• Structural Engineering and Mechanics
• /
• v.4 no.2
• /
• pp.139-148
• /
• 1996

Vector algorithms and the relative importance of the four basic modules (computation of element stiffness matrices, assembly of the global stiffness matrix, solution of the system of linear simultaneous equations, and calculation of stresses and strains) of a finite element computer program for inelastic analysis of reinforced concrete shells are presented. Performance of the vector program is compared with a scalar program. For a cooling tower problem, the speedup factor from the scalar to the vector program is 34 for the element stiffness matrices calculation, 25.3 for the assembly of global stiffness matrix, 27.5 for the equation solver, and 37.8 for stresses, strains and nodal forces computations on a Gray Y-MP. The overall speedup factor is 30.9. When the equation solver alone is vectorized, which is computationally the most intensive part of a finite element program, a speedup factor of only 1.9 is achieved. When the rest of the program is also vectorized, a large additional speedup factor of 15.9 is attained. Therefore, it is very important that all the modules in a nonlinear program are vectorized to gain the full potential of the supercomputers. The vector finite element computer program for inelastic analysis of RC shells with layered elements developed in the present study enabled us to perform mesh convergence studies. The vector program can be used for studying the ultimate behavior of RC shells and used as a design tool.

• Structural Engineering and Mechanics
• /
• v.55 no.6
• /
• pp.1241-1260
• /
• 2015

Any rational approach to define the configuration and size of viscous fluid dampers in a structure should be based on the dynamic properties of the system with the dampers. In this paper we propose an alternative representation of the complex eigenvalues of multi degree of freedom systems with dampers to calculate new equivalent natural frequencies. Analytical expressions for the dynamic properties of a two-story building model with a linear viscous damper in the first floor (i.e. with a non-proportional damping matrix) are derived. The formulas permit to obtain the equivalent damping ratios and equivalent natural frequencies for all the modes as a function of the mass, stiffness and damping coefficient for underdamped and overdamped systems. It is shown that the commonly used formula to define the equivalent natural frequency is not applicable for this type of system and for others where the damping matrix is not proportional to the mass matrix, stiffness matrix or both. Moreover, the new expressions for the equivalent natural frequencies expose a novel phenomenon; the use of viscous fluid dampers can modify the vibration frequencies of the structure. The significance of the new equivalent natural frequencies is expounded by means of a simulated free vibration test. The proposed approach may offer a new perspective to study the effect of viscous dampers on the dynamic properties of a structure.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.2
• /
• pp.288-296
• /
• 1997

Curved beam elements with two nodes based on shallow beam geometry and strain interpolations are employed in eigenvalue analysis. In these elements, the displacement interpolation functions and mass matrices are consistent with strain fields. To assess the quality of the element mass matrix in free vibration problems, several numerical experiments are performed. In these analysis, both the inconsistent mass matrices using linear displacement interpolation function and the consistent mass matrices are used to show the difference. The numerical results demonstrate that the accuracy is closely related to the property of the mass matrix as well as that of the stiffness matrix and that the mass matrix consistent with strain fields is very beneficial to eigenvalue analysis. Also, it is proved that the strain based elements are very efficient in a wide range of element aspect ratios and curvature properties.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.7
• /
• pp.99-106
• /
• 1995

Frequency response functions are of great use in dynamic analysis of structural systems. The present paper proposes an efficient method for computation of the frequency rewponse functions of linear structural dynamic models with a sparse, non-proportional damping matrix. An exact condensation procedure is proposed which enables the present method to condense the matrices without resulting in any errors. Also, an iterative scheme is proposed to be able to avoid matrix inversion in computing frequency response matrix. The proposed method is illustrated through a numerical example.

• Structural Engineering and Mechanics
• /
• v.48 no.1
• /
• pp.77-91
• /
• 2013

The paper presents the formulation of 3-nodal line semi-analytical Mindlin-Reissner finite strip in the buckling analysis of thin-walled members, which are subjected to arbitrary loads. The finite strip is simply supported in two opposite edges. The general loading and in-plane rotation techniques are used to develop this finite strip. The linear stiffness matrix and the geometric stiffness matrix of the finite strip are given in explicit forms. To validate the proposed model and study its performance, numerical examples of some thin-walled sections have been performed and the results obtained have been compared with finite element models and the published ones.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.05a
• /
• pp.326-329
• /
• 2005

In the development of sheet-handling .machinery, it is important to predict the static and dynamic behavior of the sheets with a high degree of reliability because the sheets are fed and stacked at suck a high speed flexible media behaves geometric nonlinearity of large displacement and small strain. In this paper, static analysis of flexible media are performed by FEM considering geometric nonlinearity. Linear stiffness matrix and geometric nonlinear stiffness matrix based m the updated Lagrangian approach are derived using $C^1$ beam element and numerical simulations are performed by Updated Newton-Raphson(UNR) method.

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

In the development of sheet-handling machinery, it is important to predict the static and dynamic behavior of the sheets with a high degree of reliability because the sheets are fed and stacked at such a high speed. Flexible media behaves geometric nonlinearity of large displacement and small strain. In this paper, static and dynamic analyses of flexible media are performed by FEM considering geometric nonlinearity. Linear stiffness matrix and geometric nonlinear stiffness matrix based on the Co-rotational(CR) approach are derived and numerical simulations are performed by Updated Newton-Raphson(UNR) method and Newmark integration scheme.