• Bulletin of the Society of Naval Architects of Korea
• /
• v.24 no.3
• /
• pp.1-8
• /
• 1987

In this paper, the localized finite element method(LFEM) is applied to 3-dimensional ship motion problems in water of infinite depth. The LFEM used here is based on the functional constructed by Bai & Yeung(1974). To test the present numerical scheme, a few vertical axisymmetric bodies are treated by general 3-dimensional formulation. The computed results of hydrodynamic coefficients for a few vertical spheroids and vertical circular cylinders show good agreement with results obtained by others. The advantages of the present numerical method compared with the method of integral equation are as follows; (i) The cumbersome existence of irregular frequencies in the method of conventional integral equation is removed. (ii) The final matrix is banded and symmetric and the computation of the matrix elements is comparatively easier, whereas the size of the matrix in the present scheme is much larger. (iii) In the future research, it is possible to accommodate with the nonlinear exact free surface boundary condition in the localized finite element subdomain, whereas the linear solution is assumed in the truncated(far field) subdomain.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.6 s.249
• /
• pp.705-712
• /
• 2006

A finite element method based on the penalized subdomain-interface framework is proposed for fully-coupled, nonlinear thermomechanical analyses with thermal contact anuor radiation boundaries. In the variational formulation, a well-known penalty functional scheme is adopted for connecting subdomains and interfaces that satisfy various continuity requirements. As a logical consequence, the whole domain can be arbitrarily divided into independently-modeled subdomains without considering the conformity of meshes along their interfaces. Since the nonlinearities due to the contact and radiation boundaries can be localized within a few subdomains, the computational efficiency of the present method is greatly increased with appropriate solution algorithms. By solving some numerical problems, these advantageous features are confirmed carefully.

• International Journal of Railway
• /
• v.6 no.2
• /
• pp.59-63
• /
• 2013

The purpose of this paper was to evaluate the effectiveness of geogrid for conventional ballasted track and asphalt concrete underlayment track using PLAXIS finite element program. Geogrid element was modeled at various locations that include subballast/subgrade, subballast/ballast interfaces, middle of the ballast, and one-third depth of the ballast. The results revealed that the effectiveness of geogrid reinforcement appeared to be larger for ballasted track structure compared to asphalt concrete underlayment track. Particularly, in case of installing geogrid at one-third depth of ballast layer in a conventional ballasted track, the most effectiveness of geogrid reinforcement was achieved. The influence of geogrid axial stiffness on track substructure response was not clear to conclude. Further validations using a discrete element method along with experimental investigation are considered as a future study. The effect of asphalt concrete layer modulus was evaluated. The results exhibited that higher layer modulus seems to be effective in controlling displacement and strain of track substructure. However it also yields slightly higher stresses within track substructure. It infers that further validations are required to come up with optimum asphalt concrete mixture design to meet economical and functional criteria.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.449-456
• /
• 1998

A consistent finite element formation and analytic solutions are presented for spatial stability of thin-walled circular arch. The total potential energy is derived by applying the principle of linearized virtual work and including second order terms of finite semitangential rotations. As a result the energy functional corresponding to the semitangential rotation is obtained, in which the elastic strain energy terms are considered restrained warping effects. We have obtained analytic solution for the lateral buckling of monosymmetric thin-walled curved beam subjected to pure bending or uniform compression and it's boundary conditions are simply supported. For finite element analysis, the two node cubic Hermitian polynomials are utilized as shape Auctions. In order to illustrate the accuracy of this study, parameter studies for lateral buckling problems of circular arch are presented and compared with available solutions and numerical results analyzed by the FEM using straight beam element.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.465-472
• /
• 1998

In this study, analytic solution and finite element formulation for the free vibration analysis of thin-walled circular arch, based on linearized virtual work and Vlasov's assumption, including restrained warping effect and second order terms of finite semitangential rotations, is presented. The total potential energy is derived by applying the Hellinger-Reissner principle. In this formulation, all displacement parameters of deformation are defined at the centroid axis. For the finite element formulation, the two node cubic Hermitian polynomials are utilized as shape functions. In special case, potential energy functional of thin-walled curved beam with monosymmetric cross section is derived. From this methodology, analytic solution for the free vibration of monosymmetric circular arch with simply supported is derived. In order to illustrate the accuracy of this study, various parameter studies for free vibration of circular arches are presented and compared with numerical solution analyzed by the FEM using straight beam element.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.11a
• /
• pp.173-176
• /
• 2000

Finite element model based on the Naghdi's shell theory in the general tensor-based form is formulated in the present study. Partial mixed variational functional for assumed strain is formulated in order to avoid the severe locking troubles known as transverse shear and membrane locking. The proposed assumed strain element in general tensor Naghdi's shell model provides very accurate solutions for thin shells in benchmark problems. In additions, linear elastic constitutive equations are given in the general curvilinear coordinate system including anisotropic layered structures. Thus laminated composited shell structures are easily analyzed in the present formulation.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.11a
• /
• pp.161-164
• /
• 2000

In this paper, a finite element structural analysis for thin-walled open-section composite beams with elastic couplings has been performed. The analysis includes the effects of transverse shear across beam sections, torsion warping and constrained warping. Reissner's semi-complementary energy functional is used to obtain the beam st illness coefficients The bending and torsion related warpings and the shear correct ion factors are obtained as part of the analysis. The resulting theory describes the beam kinematics in terms of the axial, flap and lag bending, flap and lag shear, torsion and torsion-warping deformations. The static response has been validated against finite element predict ions, closed form solutions, and experimental data for rectangular sol id and I-beams with elastic couplings. The free vibration results are also compared with available literature.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.04a
• /
• pp.79-86
• /
• 1997

In this study, pre-postprocessing toolbox is developed to perform elasto-plastic analyze of underground structures with transient ground water flow. This toolbox is composed of three modules. The first is the data input processor for the structural analysis. The preprocessing Is using GUI (Graphic User Interface), which is consist of dialog box, pull down, and short-cut icon, etc. The second is the structural analysis module. The analysis is based on the elasto-plastic finite element method involving additional options such as ground excavation effect, transient ground water flow, and rock bolts behavior. The last is the postprocessing module. The postprocessing is able to verify the result of the structural analysis by the graphical simulation which visualizes the element mesh, the node displacements, the element stress states, the stress contour, the ground water surface, and the rock bolt stresses. Since various options are considered separately in this toolbox, it is easy to modify the module of each processing, and to update other functional modules for the given analysis conditions.

• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.8
• /
• pp.175-184
• /
• 2003

A method of the shape design sensitivity analysis based on the boundary integral equation formulation is presented for two-dimensional inhomogeneous thermal conducting solids with multiple domains. Shape variation of the external and interface boundary is considered. A sensitivity formula of a general performance functional is derived by taking the material derivative to the boundary integral identity and by introducing an adjoint system. In numerical analysis, state variables of the primal and adjoint systems are solved by the boundary element method using quadratic elements. Two numerical examples of a compound cylinder and a thermal diffuser are taken to show implementation of the shape design sensitivity analysis. Accuracy of the present method is verified by comparing analyzed sensitivities with those by the finite difference. As application to the shape optimization, an optimal shape of the thermal diffuser is found by incorporating the sensitivity analysis algorithm in an optimization program.

• Structural Engineering and Mechanics
• /
• v.3 no.1
• /
• pp.75-89
• /
• 1995

A transition membrane element designated as CLM which has variable mid-side nodes with drilling freedoms has been presented in this paper. The functional for the linear problem, in which the drilling rotations are introduced as independent variables, has been formulated. The transition elements with variable side nodes can be efficiently used in the local mesh refinement for the in-plane structures, which have stress concentrations. A modified Gaussian quadrature is needed to be adopted to evaluate the stiffness matrices of these transition elements mainly due to the slope discontinuity of displacement within the elements. Detailed numerical studies show the excellent performance of the new transition elements developed in this study.