• Journal of the Korean Geotechnical Society
• /
• v.38 no.9
• /
• pp.45-52
• /
• 2022

Because discontinuity in the rock mass and contact of soil-structure interaction exhibits coupled thermal-hydromechanical (THM) behavior, it is necessary to develop an interface element based on the full governing equations. In this study, we derive force equilibrium, fluid continuity, and energy equilibrium equations for the interface element. Additionally, we present a stiffness matrix of the elastoplastic mechanical model for the interface element. The developed interface element uses six nodes for displacement and four nodes for water pressure and temperature in a two-dimensional analysis. The fully coupled THM analysis for fluid injection into a fault can model the complicated evolution of injection pressure due to decreasing effective stress in the fault and thermal contraction of the surrounding rock mass. However, the result of hydromechanical analysis ignoring thermal phenomena overestimates hydromechanical variables.

• Computational Structural Engineering
• /
• v.6 no.3
• /
• pp.125-137
• /
• 1993

Recently, many researches are being dealt with the adaptive method for improving the accuracy of finite element solution. This paper deals with rh-methods that are the combination of r and h-method ; r-method is to relocate the nodes for the grid optimization, h-method is to divide the elements with great error into the equal shape. As a results, rh-method has the same error decrease and convergence as h-method in the same degree of freedom, but it has more exact result of finite element in the state of restraining degrees of freedom than h-method alone.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.6
• /
• pp.527-532
• /
• 2009

In this paper simulations for the impact into ceramics and/or metal materials have been discussed. To model discrete nature for fracture and damage of brittle materials, we implemented cohesive-law fracture model with a node separation algorithm for the tensile failure and Mohr-Coulomb model for the compressive loading. The drawback of this scheme is that it requires a heavy computational time. This is because new nodes are generated continuously whenever a new crack surface is created. In order to reduce the amount of calculation, parallelization with MPI library has been implemented. For the high-speed impact problems, the mesh configuration and contact calculation changes continuously as time step advances and it causes unbalance of computational load of each processor. Dynamic load balancing technique which re-allocates the loading dynamically is used to achieve good parallel performance. Some impact problems have been simulated and the parallel performance and accuracy of the solutions are discussed.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.11 no.6
• /
• pp.2284-2291
• /
• 2010

In this paper, we investigate the vibration analysis of plates, using an 8-node shell element that accounts for the transverse shear strains and rotary inertia. The forced vibration analysis of plates subjected to arbitrary loading is investigated. In order to overcome membrane and shear locking phenomena, the assumed natural strain method is used. To improve an 8-node shell element for forced vibration analysis, the new combination of sampling points for assumed natural strain method was applied. The refined first-order shear deformation theory based on Reissner-Mindlin theory which allows the shear deformation without shear correction factor and rotary inertia effect to be considered is adopted for development of 8-node assumed strain shell element. In order to validate the finite element numerical solutions, the reference solutions of plates are presented. Results of the present theory show good agreement with the reference solution. In addition the effect of damping is investigated on the forced vibration analysis of plates.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.8 no.1
• /
• pp.129-138
• /
• 2004

In this study we determine a dynamic analysis of the existing two-span prestressed concrete box girder bridge subjected to moving vehicle loads using the experimental measurements. The moving loads applied in this paper are classified as general travelling, suddenly brake, continuous travelling, reversely travelling and reversely travelling impact loads for increasing velocities. For each travelling load, we search dynamic behaviors and characteristic in various measuring point of box girder section. In addition, the three-dimensional numerical results analyzed by the developed finite element program using flat shell element with six degrees of freedom per a node are compared with the measured experimental data. Dynamic behaviors caused impact loads by suddenly braking, reversely travelling, are bigger than by general travelling in box girder. Three-dimensional numerical results are better than one-dimensional results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1509-1517
• /
• 1993

The purpose of this study is to develop a program for the automatic derivation of the symbolic equations necessary for the kinematic and dynamic analyses of the spatial mechanism. For this purpose, a symbolic manipulation package called MCSYMA is used. Every symbolic equation is formulated using relative joint coordinate to obtain the numerically efficient system equations. These equations are produced in FORTRAN statements and linked to a FORTRAN program for numerical analysis. Several examples are taken for comparison with the commercial package called DADS which is using Cartesian coordinate approach. Also, this symbolic formulation approach is compared with a conventional numerical approach for an example. The results show that this symbolic approach with relative joint coordinate system is most efficient in computational time among three and is recommended for the derivation of macro elements frequently used.

• Journal of the Korean Recycled Construction Resources Institute
• /
• v.6 no.2
• /
• pp.100-107
• /
• 2018

The material and geometrical nonlinear finite elment analysis of UHPFRC 50M composite box girder was carried out. Constitute law in tension and compressive region of UHPFRC and HPC were modeled based on specimen test. The accuracy of nonlinear FEM analysis was verified by the experimental result of UHPFRC 50M composite girder. The UHPFRC 50M segmental composite box girder which has 1.5% steel fiber of volume fraction, 135MPa compressive strength and 18MPa tensile strength was tested. The post-tensioned UHPFRC composite girder consisted of three segment UHPFRC U-girder and High Strength Concrete reinforced slab. The parts of UHPFRC girder were modeled by 8nodes hexahedron elements and reinforcement bars and tendons were built by 2nodes linear elements by Midas FEA software. The constitutive laws of concrete materials were selected Multi-linear model both of tension and compression function under total strain crack model, which was included in classifying of smeared crack model. The nonlinearity of reinforcement elements and tendon was simulated by Von Mises criteria. The nonlinear static analysis was applied by incremental-iteration method with convergence criteria of Newton-Raphson. The validation of numerical analysis was verified by comparison with experimental result and numerical analysis result of load-deflection response, neutral axis coordinate change, and cracking pattern of girder. The load-deflection response was fitted very well with comparison to the experimental result. The finite element analysis is seen to satisfactorily predict flexural behavioral responses of post-tensioned, reinforced UHPFRC composite box girder.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.6
• /
• pp.884-894
• /
• 1987

Stress intensity factors for the bimaterial interface cracks are determined by the boundary element method employing the multiregion technique along with the double-point concept. For this purpose, the formulas relating the stress intensity factors to the crack surface displacements, which are applicable to both the homogeneous and the bimaterial systems, are derived and the accuracy of the results is discussed using the preexisting analytic solutions. Besides, the stress intensity factors for the edge-cracked bimaterial plates are computed with various crack lengths and shear modulus ratios under the biaxial and the uniaxial loadings, respectively, to demonstrate the dependence of stress intensity factors on the loading conditions and the material properties.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.2995-3006
• /
• 1993

An automatic mesh generation scheme has been developed for finite element analysis with two-dimensional, quadrilateral elements. The basic strategies of the method are to transform the analysis domain into loops with key nodes and the loops are recursively subdivided into subloops with the use of best split lines. Finally by using the basic loop operators, the meshes are completed. In this algorithm an eight-node loop operator is proposed, which is useful in the area where the change of element size is large and the splitting criteria for subdividing the loops have also been modified to the existing algorithms. Lines, arcs, and cubic spline curves are used to define the boundaries of analysis domain. Sample meshes for several geometries are presented to demonstrate the robustness of the algorithm.

• Transactions of Materials Processing
• /
• v.6 no.5
• /
• pp.408-415
• /
• 1997

The backward tracing scheme of the finite element analysis, which is counted to be unique and useful for process design in metal forming, has been developed and applied successfully in industry to several metal forming processes. Here the backward tracing scheme is implemented for process design of three-dimensional plastic deformation in metal forming, and it is applied to a precision coining process. The contact problem between the die and workpiece has been treated carefully during backward tracing simulation in three-dimensional deformation. The results confirm that the application of the developed program implemented with backward tracing scheme of the rigid plastic finite element leads to a reasonable initial piercing hole configuration. It is concluded that three-dimensional extension of the scheme appears to be successful for industrial applications.