• The Journal of the Korea Contents Association
• /
• v.11 no.1
• /
• pp.404-410
• /
• 2011

For the simulation of one-dimensional unsteady flow, the Taylor-Galerkin finite element method was adopted to the discretization of the Saint Venant equation. The model was applied to the backwater problem in a single channel and the flood routing in dendritic channel networks. The numerical solutions were compared with previously published results of finite difference and finite element methods and good agreement was observed. The model solves the continuity and the momentum equations in a sequential manner and this leads to easy implementation. Since the final system of matrix is tri-diagonal with a few additional entry due to channel junctions, the tri-diagonal matrix solution algorithm can be used with minor modification. So it is fast and economical in terms of memory for storing matrices.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.11
• /
• pp.1538-1551
• /
• 1997

Kinetic energy conservation for fixed and moving grids is examined in time-accurate finite element computation of fully unsteady inviscid flows. As numerical algorithms, fractional step method (FSM) and modified SIMPLE are used. To simulate the flow in moving grid system, arbitrary Lagrangian-Eulerian (ALE) method is adopted. In the present study, the energy conserving time integration rule for finite element algorithm is proposed and discussed schematically. It is shown that the discretization by Crank-Nicolson in time and Galerkin (central difference) in space must be used to ensure energy conservation. The developed code has been tested for a standing vortex in fixed or moving grid system, sloshing in a tank and propagation of a solitary wave, and has been shown to be a completely energy conserving algorithm.

• Steel and Composite Structures
• /
• v.27 no.3
• /
• pp.255-271
• /
• 2018

This paper deals with the transient dynamic analysis and elastic wave propagation in a functionally graded graphene platelets (FGGPLs)-reinforced composite thick hollow cylinder, which is subjected to shock loading. A micromechanical model based on the Halpin-Tsai model and rule of mixture is modified for nonlinear functionally graded distributions of graphene platelets (GPLs) in polymer matrix of composites. The governing equations are derived for an axisymmetric FGGPLs-reinforced composite cylinder with a finite length and then solved using a hybrid meshless method based on the generalized finite difference (GFD) and Newmark finite difference methods. A numerical time discretization is performed for the dynamic problem using the Newmark method. The dynamic behaviors of the displacements and stresses are obtained and discussed in detail using the modified micromechanical model and meshless GFD method. The effects of the reinforcement of the composite cylinder by GPLs on the elastic wave propagations in both displacement and stress fields are obtained for various parameters. It is concluded that the proposed micromechanical model and also the meshless GFD method have a high capability to simulate the composite structures under shock loadings, which are reinforced by FGGPLs. It is shown that the modified micromechanical model and solution technique based on the meshless GFD method are accurate. Also, the time histories of the field variables are shown for various parameters.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1992.10a
• /
• pp.117-123
• /
• 1992

In recent years, the finite element method has become one of the most popular numerical technique for obtaining solutions of engineering science problems. However, there exist various uncertainties in modeling the problems, such as the dimensions(geometry shape), the material properties, boundary conditions, etc. The consideration for the uncertainties inherent in the problems can be made by understanding the influences of uncertain parameters[1]. Determining the influences of uncertainties as statistical quantities using the standard finite element method requires enormous computing time, while the probabilistic finite element method is realized as an efficient scheme[2,3] yielding statistical solution with just a few direct computations. In this paper, a formulation of the probabilistic fluid-structure interaction problem accounting for the first order perturbation of geometric shape is derived, and especially probabilistical acoustic pressure scattering from the structure with surrounding fluid is focused on. In Section 2, governing equations for the fluid-structure problems are given. In Section 3, a finite element formulation, based on the functional, is presented. First order perturbation of geometric shape with randomness is incorporated into the finite element formulation in conjunction with discretization of the random fields in Section 4 and 5. Finally, the proposed formulation is applied to a acoustic pressure scattering problem from an infinitely long cylindrical shell structure with randomness of radial perturbation.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.06a
• /
• pp.309-313
• /
• 2000

A structural-acoustic coupling problem involving fluid in a cavity divided with flexible walls and porous materials is investigated in this paper. In many practical problems, the use of finite elements to discretize the fluid region leads to large stiffness and mass matrices. But, since the acoustic boundary element discretization requires to put elements only on the surface of structure, the size of matrices is reduced considerably. Here, we developed a numerical analysis program for the structural-acoustic coupling problems of the multi-region cavity, using boundary elements for the fluid regions and finite elements for the structure. By considering sound transmission through layered systems placed in a cavity, the accuracy of the coupled acoustical-structural finite element model has been verified by comparing its transmission loss predictions with analytical sloutions. Example problems are included to investigate the characteristics of the multi-region structural-acoustic coupling system with porous material.

• Journal of the Korean Society of Safety
• /
• v.26 no.2
• /
• pp.62-69
• /
• 2011

This paper presents a new simple two-dimensional frame finite element able to accurately estimate the load-carrying capacity of reinforced concrete beams flexurally strengthened externally bonded fiber reinforced polymer (FRP) strips and plates. The proposed analysis model considers distributed plasticity with layer-discretization of the cross-sections and the bond-slip behavior of epoxy layer. The proposed model is used to predict the load-carrying capacity and the applied load-midspan deflection response of RC beams subjected to bending loading. Numerical simulations and experimental measurements are compared based on numerous tests available in the literature and published by different authors. The numerically simulated response agree remarkably well with the corresponding experimental results. Thus, the proposed model is suitable for efficient and accurate modeling and analysis of flexural strengthening of RC beams with externally bonded FRP sheets/plates and for practical use in design-oriented parametric studies.

• Ocean Systems Engineering
• /
• v.8 no.3
• /
• pp.329-343
• /
• 2018

The inclination of seabed profile (sloped seabed) is one of the known topographic features which can be observed at different seabed level in the large offshore basin. A mooring system connected between the platform and global seabed is an integral part of the floating structure which tries to keep the floating platform settled in its own position against hostile sea environment. This paper deals with an investigation of the motion responses of an FPSO platform moored on the sloped seabed under the combined action of wave, wind and current loads. A three-dimensional panel discretization method has been used to model the floating body. To introduce the connection of multi-segmented non-linear elastic catenary mooring cables with the sloped seabed, a quasi-static composite catenary model is employed. The model and analysis have been completed by using hydrodynamic diffraction code AQWA. Validation of the numerical model has been successfully carried out with an experimental work published in the latest literature. The analysis procedure in this study has been followed time domain analysis. The study involves an objective oriented investigation on platform motions, in order to identify the effects of the slopped seabed, the action of the wave, wind and current loads and the presence of riser system. In the end, an effective analysis has been performed to identify a stable mooring model in demand of reducing structural responses of the FPSO.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.3
• /
• pp.53-61
• /
• 1984

A nonlinear formulation of a beam finite element(NB6) on the total Lagrangian mode for the geometrically nonlinear analysis of two-dimensional elastic framed structures is presented. The NB6 beam element has been degenerated from the three-dimensional continuum by introducing the deep beam assumptions and consists of three reference nodes and three relative nodes. The element characteristics are derived by discretizing the beam equations of motion using the Galerkin weighted residual method and are reduced-integrated repeatedly for each loading step by the Newton-Raphson iteration techpique. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed nonlinear NB6 beam element.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.3
• /
• pp.1-10
• /
• 1988

A higher order plate bending finite element using cubic in-plane displacement profiles is proposed for geometrically nonlinear analysis of thin and thick plates. The higher order plate bending element has been derived from the three dimensional plate-like continuum by discretization of the equations of motion by Galerkin weighted residual method, together with enforcing higher order plate assumptions. Total Lagrangian formulation has been used for geometrically nonlinear analysis of plates and consistent linearization by Newton-Raphson method has been performed to solve the nonlinear equations. The element characteristics have been computed by, selective reduced integration technique using Gauss quadrature to avoid shear locking phenomenon in case of extremely thin plates. Several numerical examples were solved with FEAP macro program to demonstrate versatility and accuracy of the present higher order plate bending element.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.1-8
• /
• 1992

This paper presents a p-version finite element approach for modeling the stress distribution around a circular hole in a finite strip subjected to membrane and flexural behaviors. Also, same problem with a crack emanating from a perforated tension strip was solved by virtual crack extension method. The p-version of the finite element method based on integrals of Legendre polynomials is shown to perform very well for modeling geometries with very steep stress gradients in the vicinity of a circular cutout. Here, the transfinite mapping technique for circular boundaries was used to avoid the discretization errors. The numerical results from the proposed scheme have a good comparison with those by Nisida, Howland, Newman etc. and the conventional finite element approach.