• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.445-454
• /
• 1991

In this thesis, Mindlin plate element with nine nodes and three degrees-of-freedom at each node is formulated and is employed in eigen-analysis of a rectangular plates in order to alleviate locking phenomenon of eigenvalues. Eigenvalues and their modes may be locked if conventional $C_{0}$-isoparametric element is used. In order to reduce stiffness locking phenomenon, two methods (1, the general reduced and selective integration, 2, the new element that use of modified shape function) are studied. Additionally in order to reduce the error due to mass matrix, two mass matrixes (1, Gauss-Legendre mass matrix, 2, Gauss-Lobatto mass matrix) are considered. The results of eigen-analysis for two models (the square plate with all edges simply-supported and all edges built-in), computed by two methods for stiffness matrix and by two mass matrixes are compared with theoretical solutions and conventional numerical solutions. These comparisons show that the performance of the two methods with Gauss-Lobatto mass matrix is better than that of the conventional plate element. But, by considering the spurious rigid body motions, the element which employs modified shape function with full integration and Gauss-Lobatto mass matrix can elevate the accuracy and convergence of numerical solutions.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2008.11a
• /
• pp.407-413
• /
• 2008

This paper investigates the dynamics of a HDD spindle system due to the change of FDBs. Flying height of the HDD spindle system is determined through the static analysis of the FDBs, and the stiffness and damping coefficients are calculated through the dynamic analysis of the FDBs. Free vibration characteristics and shock response of the HDD spindle system are analyzed by using the finite element method and the mode superposition method. Experimental modal test is also performed to verify the accuracy of the proposed method. This research shows that the stiffness coefficients of journal heating mostly affect the rocking frequencies because their magnitude are within the range of the stiffness of supporting structure. It also shows that the damping coefficients of thrust bearing mostly affect the axial frequency because the stiffness of thrust bearing is much smaller that that of supporting structure.

• Structural Engineering and Mechanics
• /
• v.37 no.5
• /
• pp.475-487
• /
• 2011

In a flexible multi-body dynamic system the typical topological optimization method for structures cannot be directly applied, as the stiffness varies with position. In this paper, the topological optimization of the flexible multi-body dynamic system is converted into structural optimization using the equivalent static load method. First, the actual boundary conditions of the control system and the approximate stiffness curve of the mechanism are obtained from a flexible multi-body dynamical simulation. Second, the finite element models are built using the absolute nodal coordination for different positions according to the stiffness curve. For efficiency, the static reanalysis method is utilized to solve these finite element equilibrium equations. Specifically, the finite element equilibrium equations of key points in the stiffness curve are fully solved as the initial solution, and the following equilibrium equations are solved using a reanalysis method with an error controlled epsilon algorithm. In order to identify the efficiency of the elements, a non-dimensional measurement is introduced. Finally, an improved evolutional structural optimization (ESO) method is used to solve the optimization problem. The presented method is applied to the optimal design of a die bonder. The numerical results show that the presented method is practical and efficient when optimizing the design of the mechanism.

• Smart Structures and Systems
• /
• v.3 no.3
• /
• pp.245-261
• /
• 2007

Nonlinear time-domain system identification (SI) algorithm is proposed to assess damage in a shear building by synchronously estimating time-varying stiffness and damping parameters using measured acceleration data. Mass properties have been assumed as the a priori known information. Viscous damping was utilized for the current research. To chase possible nonlinear dynamic behavior under severe vibration, an incremental governing equation of vibrational motion has been utilized. Stiffness and damping parameters are estimated at each time step by minimizing the response error between measured and computed acceleration increments at the measured degrees-of-freedom. To solve a nonlinear constrained optimization problem for optimal structural parameters, sensitivities of acceleration increment were formulated with respect to stiffness and damping parameters, respectively. Incremental state vectors of vibrational motion were computed numerically by Newmark-${\beta}$ method. No model is pre-defined in the proposed algorithm for recovering the nonlinear response. A time-window scheme together with Monte Carlo iterations was utilized to estimate parameters with noise polluted sparse measured acceleration. A moving average scheme was applied to estimate the time-varying trend of structural parameters in all the examples. To examine the proposed SI algorithm, simulation studies were carried out intensively with sample shear buildings under earthquake excitations. In addition, the algorithm was applied to assess damage with laboratory test data obtained from free vibration on a three-story shear building model.

• Structural Engineering and Mechanics
• /
• v.62 no.6
• /
• pp.759-769
• /
• 2017

The effect of central axial load on natural frequencies of various thin-walled beams, are investigated by some researchers using different methods such as finite element, transfer matrix and dynamic stiffness matrix methods. However, there are situations that the load will be off centre. This type of loading is called eccentric load. The effect of the eccentricity of axial load on the natural frequencies of asymmetric thin-walled beams is a subject that has not been investigated so far. In this paper, the mentioned effect is studied using exact dynamic stiffness matrix method. Flexure and torsion of the aforesaid thin-walled beam is based on the Bernoulli-Euler and Vlasov theories, respectively. Therefore, the intended thin-walled beam has flexural rigidity, saint-venant torsional rigidity and warping rigidity. In this paper, the Hamilton‟s principle is used for deriving governing partial differential equations of motion and force boundary conditions. Throughout the process, the uniform distribution of mass in the member is accounted for exactly and thus necessitates the solution of a transcendental eigenvalue problem. This is accomplished using the Wittrick-Williams algorithm. Finally, in order to verify the accuracy of the presented theory, the numerical solutions are given and compared with the results that are available in the literature and finite element solutions using ABAQUS software.

• Structural Engineering and Mechanics
• /
• v.21 no.5
• /
• pp.539-551
• /
• 2005

Finite elements based on isoparametric formulation are known to suffer spurious stiffness properties and corresponding stress oscillations, even when care is taken to ensure that completeness and continuity requirements are enforced. This occurs frequently when the physics of the problem requires multiple strain components to be defined. This kind of error, commonly known as locking, can be circumvented by using reduced integration techniques to evaluate the element stiffness matrices instead of the full integration that is mathematically prescribed. However, the reduced integration technique itself can have a further drawback - rank deficiency, which physically implies that spurious energy modes (e.g., hourglass modes) are introduced because of reduced integration. Such instability in an existing stiffness matrix is generally detected by means of an eigenvalue test. In this paper we show that a knowledge of the dimension of the solution space spanned by the column vectors of the strain-displacement matrix can be used to identify the instabilities arising in an element due to reduced/selective integration techniques a priori, without having to complete the element stiffness matrix formulation and then test for zero eigenvalues.

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

• Proceedings of the KSME Conference
• /
• 2000.11a
• /
• pp.363-369
• /
• 2000

Due to environmental problem for reduction in fuel consumption, vehicle emission and etc., many automotive makers are trying to reduce the weight of the vehicle. The most effective way to reduce the weight of vehicle is to use lighter materials, aluminum, plastics. Aluminum Space Frame has many advantages in weight reduction, body stiffness, ease of model change and so on. So, most of automotive manufacturers are attempting to develope Aluminum Space Frame body. For these reasons, we have developed Aluminum Intensive Vehicle based on steel monocoque body with Hyundai Motor Company. We achieved about 30% weight reduction, the stiffness of our model was higher than that of conventional steel monocoque body. In this paper, with optimization using FEM analysis, we could get more weight reduction and body stiffness increase. In the long run, we analyzed by means of simulation using PAM-CRASH to evaluate crush and crash characteristic of Aluminum Intensive Vehicle in comparison to steel monocoque automotive.

• Journal of Magnetics
• /
• v.18 no.4
• /
• pp.443-453
• /
• 2013

Maglev wind turbine generator (MWTG) technology has been widely studied due to its low loss, low maintenance cost, and high reliability. However, the dynamics of the magnetic bearing system differ fromthe those of the traditional mechanical bearing system. A horizontal axial MWTG supported with a permanent magnetic bearing is designed in this research and the radial forces and the natural frequencies of the rotor system are studied. The results show that the generatorhas a cyclical magnetic forceand an unreasonable bearing stiffness may mean that the rotor system needs to work in the resonance region; the bearing stiffness is the key factor to avoid this problem. This is the main rule of the bearing stiffness design in the MWTG, and this rule can also be used in other maglev permanent magnet motors.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.487-492
• /
• 2007

The enhancement of the service life of damaged or cracked structures is a major issue for researchers and engineers. The hierarchical void element with the integrals of Legend polynomials is used to characterize the fracture behavior of unpatched crack as well as repaired crack with bonded composite patches by computing the stress intensity factors and stress contours at the crack tip. The numerical approach is based on the v-version degenerate shell element including the theory of anisotropic laminated composites. Since the equivalent single layer approach is adopted in this study, the proposed element is necessary to represent a discontinuous crack part as a continuum body with zero stiffness of materials. Thus the aspect ratio of this element to represent the crack should be extremely slender. The sensitivity of numerical solution with respect to energy release rate, displacement and stress has been tested to show the robustness of hierarchical void element as the aspect ratio is increased up to 2000. The stiffness derivative method and displacement extrapolation method have been applied to calculate the stress intensity factors of Mode I problem.