• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.2A
• /
• pp.131-136
• /
• 2010

In this study, most of famous algorithms for the corner problem are listed. By comparing these with implemented codes and theoretical dissections, new algorithms are developed. These algorithms are combined by the existing auxiliary equations. All relating algorithms are numerically tested with 3 problems. Two problems have well-known analytical solutions and the result of another example is compared with the one of the published paper. The conducted research reveals the characteristics of existing algorithms and demonstrates newly developed algorithms can produce a reasonable solution by reflecting various type of boundary conditions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.109-122
• /
• 1983

An accurate thick beam element (TB4) which includes the effects of the shear deformation and rotatory inertia has been degenerated from the three dimensional continuum by employing the Timoshenko beam assumptions. The proposed TB4 element has four nodes and two degrees of freedom at each node, totally eight degrees of freedom. The transverse deflection W and plane rotation ${\theta}$ with the cubic interpolation functions are selected as nodal variables. The element characteristics are formulated by discretizing the beam equations of motion, using the Galerkin weighted residual method, and are numerically integrated by the reduced shear integration technique, using the three-point Gauss quadrature with the various shear coefficients. Several numerical examples are analyzed to demonstrate the accuracy and the monotonic convergence behavior of the proposed TB4 beam element. The result indicates that the TB4 element shows the more excellent performance and the monotonic convergence behavior than the other existing Timoshenko beam type elements for the whole range of the beam aspect ratios, in both static and free vibration analyses.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.16 no.10
• /
• pp.6673-6678
• /
• 2015

Recently, the free form buildings are constructed frequently. Exterior and interior components of these buildings have the free cross-section and a curved shape. So, There are many usages of classical finite element having tapered section and free-style shape. Some general commercial applications like ETABS, SAP2000, MIDAS are usually used for the safety evaluation of the free form structures. However, there are some limits in the accuracy of structural analysis and the length of analysis time because a very complicated finite element mesh have to be used. Therefore, In this study, a pre and post program module was developed to take advantage of general 3-D curved beam element which has a free-style curved shape and mathematical backgrounds. Pre-post processing module has been developed in this study was developed to control the curvature of the curved members by the NURBS control points. As a result, fast geometric modeling than was possible commercial applications. In addition, realistic depiction of the shape and behavior patterns were possible because of the free-form building allows visual check of the free form.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.6A
• /
• pp.577-585
• /
• 2009

A 9-node degenerated shell finite element(FE), which has been developed for assessment of ultimate pressure capacity and nonlinear analysis for nuclear containment building is described in this paper. Reissner-Midnlin(RM) assumptions are adopted to develop the shell FE so that transverse shear deformation effects is considered. Material model for concrete prior to cracking is constructed based on the equivalent stress-equivalent strain relationship. Tension stiffening model, shear transfer mechanism and compressive strength reduction model are used to model the material behavior of concrete after cracking. Niwa and Aoyagi-Yamada failure criteria have been adapted to find initial cracking point in compression-tension and tension-tension region, respectively. Finally, the performance of the developed program is tested and demonstrated with several examples. From the numerical tests, the present results show a good agreement with experimental data or other numerical results.

• 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.26 no.1
• /
• pp.39-48
• /
• 2013

This paper develops a 2-D moving least squares(MLS) difference method for Stefan problem by extending the 1-D version of the conventional method. Unlike to 1-D interfacial modeling, the complex topology change in 2-D domain due to arbitrarily moving boundary is successfully modelled. The MLS derivative approximation that drives the kinetics of moving boundary is derived while the strong merit of MLS Difference Method that utilizes only nodal computation is effectively conserved. The governing equations are differentiated by an implicit scheme for achieving numerical stability and the moving boundary is updated by an explicit scheme for maximizing numerical efficiency. Numerical experiments prove that the MLS Difference Method shows very good accuracy and efficiency in solving complex 2-D Stefan problems.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.7
• /
• pp.44-50
• /
• 2002

Low-velocity impact on composite sandwich panel has been investigated. Contact force is computed from a proposed modified Hertzian contact law. The Hertzian contact law is constructed by adjusting numerical value of the exponent and reducing the through-the- thickness elastic constant of honeycomb core. The equivalent transverse elastic constant is calculated from the rule of mixture. Nonlinear equation to calculate the contact force is solved by the Newton-Raphson method and time integration is done by the Newmark-beta method. A finite element program for the low-velocity impact analysis is coded by implementing these techniques and an 18-node assumed strain solid element. Behaviors of composite sandwich panels subjected to low-velocity impact are analyzed for various cases with different geometry and lay-ups. It has been found that the present code with the proposed contact law can predict measured contact forces and contact times for most cases within reasonable error bounds.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.6
• /
• pp.359-366
• /
• 2020

In this study, a simulation-based damage estimation method for helidecks is proposed using an artificial neural network. The structural members that share a connecting node in the helideck are regarded as a damage group, and a total of 37,400 damage scenarios are numerically generated by applying randomly assigned damage to up to three damage groups. Modal analysis is then performed for all the damage scenarios, which are selectively used as either training or validation or verification sets based on the purpose of use. An artificial neural network with three hidden layers is constructed using a PyTorch program to recognize the patterns of the modal responses of the helideck model under both damaged and undamaged states, and the network is successively trained to minimize the loss function. Finally, the estimated damage rate from the proposed artificial neural network is compared to the actual assigned damage rate using 400 verification scenarios to show that the neural network is able to estimate the location and amount of structural damage precisely.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.21 no.1
• /
• pp.72-82
• /
• 2015

In order to effectively and economically apply the previous POM(Princeton Ocean Model) WAD(Wetting And Drying) to the coastal area, the POM WAD was modified such as the water elevation input of tidal harmonics in the open boundaries was included and a CTS(Computing Time Saving) technique was introduced to the model. The modified model was tested to the standing waves in the rectangular bay and the hydraulic experiments for the flow and heat diffusion in the 3D basin. The numerical results showed a good agreement with the analytical solutions of the standing waves and the observed values by the hydraulic experiments, respectively. And also when the modified model with the CTS technique was applied to Gwangyang Bay of Korea, the computing time was decreased by as much as 39.4%.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.5 no.3
• /
• pp.49-59
• /
• 1985

Formulation of the geometric optimization for truss structures based on the elasticity theory turn out to be the nonlinear programming problem which has to deal with the Cross sectional area of the member and the coordinates of its nodes simultaneously. A few techniques have been proposed and adopted for the analysis of this nonlinear programming problem for the time being. These techniques, however, bear some limitations on truss shapes loading conditions and design criteria for the practical application to real structures. A generalized algorithm for the geometric optimization of the truss structures which can eliminate the above mentioned limitations, is developed in this study. The algorithm developed utilizes the two-phases technique. In the first phase, the cross sectional area of the truss member is optimized by transforming the nonlinear problem into SUMT, and solving SUMT utilizing the modified Newton-Raphson method. In the second phase, the geometric shape is optimized utilizing the unidirctional search technique of the Rosenbrock method which make it possible to minimize only the objective function. The algorithm developed in this study is numerically tested for several truss structures with various shapes, loading conditions and design criteria, and compared with the results of the other algorithms to examme its applicability and stability. The numerical comparisons show that the two-phases algorithm developed in this study is safely applicable to any design criteria, and the convergency rate is very fast and stable compared with other iteration methods for the geometric optimization of truss structures.