• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.6
• /
• pp.507-515
• /
• 2017

Flexible multibody simulations are widely used in the industry to design mechanical systems. In flexible multibody dynamics, deformation coordinates are described either relatively in the body reference frame that is floating in the space or in the inertial reference frame. Moreover, these deformation coordinates are generated based on the discretization of the body according to the finite element approach. Therefore, the formulation of the flexible multibody system always deals with a huge number of degrees of freedom and the numerical solution methods require a substantial amount of computational time. Parallel computational methods are a solution for efficient computation. However, most of the parallel computational methods are focused on the efficient solution of large-sized linear equations. For multibody analysis, we need to develop an efficient formulation that could be suitable for parallel computation. In this paper, we developed a subsystem synthesis method for a flexible multibody system and proposed efficient parallel computational schemes based on the OpenMP API in order to achieve efficient computation. Simulations of a rotating blade system, which consists of three identical blades, were carried out with two different parallel computational schemes. Actual CPU times were measured to investigate the efficiency of the proposed parallel schemes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.37-49
• /
• 2000

This paper is for the first essential study on the development of unified finite element formulations for solving problems related to the dynamics/thermoelastics behavior of solids. In the first part of formulations, the finite element method is based on the introduction of a new quantity defined as heat displacement, which allows the heat conduction equations to be written in a form equivalent to the equation of motion, and the equations of coupled thermoelasticity to be written in a unified form. The equations obtained are used to express a variational formulation which, together with the concept of generalized coordinates, yields a set of differential equations with the time as an independent variable. Using the Laplace transform, the resulting finite element equations are described in the transform domain. In the second, the Laplace transform is applied to both the equation of heat conduction derived in the first part and the equations of motions and their corresponding boundary conditions, which is referred to the transformed equation. Selections of interpolation functions dependent on only the space variable and an application of the weighted residual method to the coupled equation result in the necessary finite element matrices in the transformed domain. Finally, to prove the validity of two approaches, a comparison with one finite element equation and the other is made term by term.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.51-61
• /
• 2000

This paper discusses the homogenization method to determine effective average elastic constants of a linear structure by considering its microstructure. A detailed description on the homogenization method is given for the linear elastic material and then the finite element approximation is performed for an investigation of elastic properties. An asymptotic expansion is carried out in the cross-section area, or in the unit cell. Two and three lay-up structures made up of individual isotropic constituents are chosen for numerical examples to check discrepancies between results generated by this theoretical development and the conventional approach. Asymptotic characteristics of the process in extracting the stiffness of structure locally formed by spatial repetitions yield underestimated values of stiffness. These discrepancies are detected by the asymptotic corrective term which is ascribed to considerations of microscopic perturbations and proved in the finite element formulation. The asymptotic analysis is the more reasonable in analysing the composite material, rather than the conventional approach to calculate the macroscopic average for elastic properties.

• Journal of Korean Society of Steel Construction
• /
• v.11 no.2 s.39
• /
• pp.201-212
• /
• 1999

A beam supported by a flexible elastic support is commonly used as structural elements, e.g., braced beam, railway track, etc. The elastic support can be located in arbitrary point in the cross-section. This paper investigates the effects of support eccentricity on the elastic buckling of beams with elastic supports. The effects of stiffness of the elastic support are also studied. A beam element with elastic supports and the analysis program are developed for elastic buckling analysis using finite element formulation. The elastic support is modeled by elastic spring element. Using the offset technique, the eccentricity of support is taken into account. A beam element having 14 degrees of freedom including the warping degree of freedom is used. Various numerical example analyses show that the present formulation and analysis program accurately and effectively compute the buckling load and mode of beams with elastic supports.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2008.04a
• /
• pp.1053-1056
• /
• 2008

Permeability coefficient of concrete is a substantial key parameter for understanding the durability performance of concrete and its micro-structural densification. Many researches to deal with the issue have been accomplished, however, it is very rare to deal with the theoretical study on permeability coefficient in connection with carbonation of concrete and the effect of volumetric fraction of cement paste or aggregate on the permeability coefficient. The majority of these researches have not dealt with this issue combined with carbonation of concrete, although carbonation can significantly impact on the permeability coefficient of concrete. The purpose of this study is to establish a fundamental approach to compute the permeability coefficient of (non)carbonated concrete. When simulating micro-structural characteristics as a starting point for deriving a model for the permeability coefficient by the numerical simulation program for cementitious materials, HYMOSTRUC, a more realistic formulation can be achieved. For several compositions of cement pastes, the permeability coefficient is calculated with the analytical formulation, followed by a microstructure-based model. Emphasis is on the micro-structural changes and its effective change of the permeability coefficient of carbonated concrete. The results of micro-structural water permeability coefficient model will be compared with results achieved from permeability experiments.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.2
• /
• pp.153-160
• /
• 2008

In this study, we propose a new three-node hybrid-mixed curved beam element with the relaxed-equiribrium stress functions for static analysis. The proposed element considering shear deformation is based on the Hellinger-Reissner variational principle. The stress functions are carefully chosen from three important considerations: (i) all the kinematic deformation modes must be suppressed, and (ii) the spurious constraints must be removed in the limiting behaviors via the field-consistency, and (iii) the relaxed equilibrium conditions could be incorporated because it might be impossible to select the stress functions and parameters to fully satisfy both the equiribrium conditions and the suppression of kinematic deformation modes in the three-node curved beam hybrid-mixed formulation. Numerical examples confirm the superior and stable behavior of the proposed element regardless of slenderness ratio and curvature. Besides, the proposed element shows the outstanding performance in predicting the stress resultant distributions.

• Structural Engineering and Mechanics
• /
• v.24 no.6
• /
• pp.647-664
• /
• 2006

In this study, an improved finite element formulation with a scheme of solution for the large deflection analysis of inextensible prismatic and nonprismatic slender beams is developed. For this purpose, a three-noded Lagrangian beam-element with two dependent degrees of freedom per node (i.e., the vertical displacement, y, and the actual slope, $dy/ds=sin{\theta}$, where s is the curved coordinate along the deflected beam) is used to derive the element stiffness matrix. The element stiffness matrix in the global xy-coordinate system is achieved by means of coordinate transformation of a highly nonlinear ($6{\times}6$) element matrix in the local sy-coordinate. Because of bending with large curvature, highly nonlinear expressions are developed within the global stiffness matrix. To achieve the solution after specifying the proper loading and boundary conditions, an iterative quasi-linearization technique with successive corrections are employed considering these nonlinear expressions to remain constant during all iterations of the solution. In order to verify the validity and the accuracy of this study, the vertical and the horizontal displacements of prismatic and nonprismatic beams subjected to various cases of loading and boundary conditions are evaluated and compared with analytic solutions and numerical results by available references and the results by ADINA, and excellent agreements were achieved. The main advantage of the present technique is that the solution is directly obtained, i.e., non-incremental approach, using few iterations (3 to 6 iterations) and without the need to split the stiffness matrix into elastic and geometric matrices.

• Structural Engineering and Mechanics
• /
• v.15 no.1
• /
• pp.83-114
• /
• 2003

A new simple 3-node triangular flat-shell element with standard nodal DOF (6 DOF per node) is proposed for the linear and geometrically nonlinear analysis of very thin to thick plate and shell structures. The formulation of element GT9 (Long and Xu 1994), a generalized conforming membrane element with rigid rotational freedoms, is employed as the membrane component of the new shell element. Both one-point reduced integration scheme and a corresponding stabilization matrix are adopted for avoiding membrane locking and hourglass phenomenon. The bending component of the new element comes from a new generalized conforming Kirchhoff-Mindlin plate element TSL-T9, which is derived in this paper based on semiLoof constrains and rational shear interpolation. Thus the convergence can be guaranteed and no shear locking will happen. Furthermore, a simple hybrid procedure is suggested to improve the stress solutions, and the Updated Lagrangian formulae are also established for the geometrically nonlinear problems. Numerical results with solutions, which are solved by some other recent element models and the models in the commercial finite element software ABAQUS, are presented. They show that the proposed element, denoted as GMST18, exhibits excellent and better performance for the analysis of thin-think plates and shells in both linear and geometrically nonlinear problems.

• Journal of the Korean Chemical Society
• /
• v.20 no.2
• /
• pp.97-110
• /
• 1976

Effects of direct multiple quantum excitations in vibrational energy transfer were investigated. Vibrational transition probabilities for 0${\rightarrow}$2, 0${\rightarrow}$3, and 0${\rightarrow}$4 excitations were explicitly formulated including both direct 0→n excitations and stepwise single quantum processes. For the formulation the perturbing force was derived from the exponential potential including terms up to fourth order in the vibrational amplitude. The head-on collinear collision model between a harmonic oscillator and an incident particle was employed, and the formulation was based on the semiclassical approximation. Numerical results were obtained for five different collision systems (Ar${\cdots}$O-N, He${\cdots}$H-H, He${\cdots}$H-Cl, 5${\cdots}$1-2, 2${\cdots}$12-12). Comparison between the present results and those obtained using the linearized interaction potential showed that the overall effect of including the direct multiple quantum transition is to decrease the probabilities at low collision energies and to increase them at high energies. The present results were found to be significantly different from those obtained using the linearized potential for collision systems He${\cdots}$H-H, He${\cdots}$H-Cl, and 5${\cdots}$1-2. For systems Ar${\cdots}$O-N and 2${\cdots}$12-12 the differences were negligible.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.12
• /
• pp.17-25
• /
• 2012

FEM requires the stress-strain relationship equations for numerical analyses. However, most formulations for the stress-strain relationship published up to the present are not satisfactory enough to properly express all the levels from the small strain to the peak. Tatsuoka and Shibuya (1991) suggested a new single formulation applicable not only to a wide range of geo-materials from soft clay to soft rock, but also to a wide range of strain levels from $10^{-6}$ to $10^{-2}$. The plain strain compression test is carried out to seven samples of research standard sand specimens and two samples of glass beads, which have been used at world-renowned research institutes. In this study, strains of the maximum principal stress (${\sigma}_1$) and the minimum principal stress (${\sigma}_3$) were thoroughly measured from $10^{-6}$ to $10^{-2}$, and the result, applied to Tatsuoka and Shibuya's new formulation, coincided closely with the measured data of the stress-strain relationship from the small strain to the peak.