• Journal of the Korean Geotechnical Society
• /
• v.15 no.4
• /
• pp.221-234
• /
• 1999

This paper deals with the application of joint element in the finite element modeling of discontinuities encountered during rock tunneling. A nodal displacement joint element was implemented in a two dimensional finite element program GEOFE2D. The applicability of the joint element for modeling of discontinuities and the numerical stability of the implemented algorithm were examined by comparing the results of reduced small scale model tests as well as commercially available FEM program. The GEOFE2D was then used to analyze a tunnel crossed by a major discontinuity for the purpose of understanding the effect of discontinuity on the tunnel behavior. In addition, a modeling technique for the junction of discontinuity and shotcrete lining was presented. The results of analysis indicated that the stress-strain field around the tunnel is significantly altered by the presence of discontinuity, and that the stresses in the shotcrete lining considerably increase at the junction of the shotcrete lining and the discontinuity. It is therefore concluded that the major discontinuities must be carefully modeled in the finite element analysis of a tunneling problem in order to obtain more reliable results close to actual tunnel behavior.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.14 no.2
• /
• pp.174-183
• /
• 2006

The resistance spot-welded region in most current finite element crash models is characterized as a rigid beam at the location of the welded spot. The region is modeled to fail with a failure criterion which is a function of the axial and shear load at the rigid beam. The role of this rigid beam is simply to transfer the load across the welded components. The calculation of the load acting on the rigid beam is important to evaluate the failure of the spot-weld. In this paper, numerical simulation is carried out to evaluate the calculation of the load at the rigid beam. The load calculated from the precise finite element model of the spot-welded region considering the residual stress due to the thermal history during the spot welding procedure is regarded as the reference value and the value of the load is compared with the one obtained from the spot-welded model using the rigid beam with respect to the element size, the element shape and the number of imposed constraints. Analysis results demonstrate that the load acting on the spot-welded element is correctly calculated by the change of the element shape around the welded region and the location of welded constrains. The results provide a guideline for an accurate finite element modeling of the spot-welded region in the crash analysis of vehicles.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.11
• /
• pp.3441-3453
• /
• 1996

A finite element anlaysis is performed for large deformations of a felxible beam. The total Lagrangian formulation for a general large deformation, which involves finite rotations, is chosen and the exponential map is used to treat finite rotations from the Eulerian point of view. The finite elements results are confirmed for several cases of deformations through comparison to a first order elasticity solution obtained by numerical integration, and the agreement between the two is found to be excellent. For lateral buckling, the point of vanishing determinant of the resulting unsymmetric tangent stiffness is traced to examine its relationship to bifurcation points. It is found that the points of vanishing determinant is not corresponding to bifurcation points for large deformation in general, which suggests that the present unsymmetric tangent stiffness is not an exact first derivative of internal forces with respect to displacement.

• Proceedings of the KSME Conference
• /
• 2000.11a
• /
• pp.486-491
• /
• 2000

In this study, finite elements addition and removal method by stress range is applied to optimize shapes in structures, without using classical and numerical optimization methods and search methods. The program based on this algorithm is developed and compared to theoritial results with considerable accuracy. Classical methods need mesh generation for finite element analysis for every iteration, the developed method needs updated mesh data such as coordinates of nodes, elements connectivity, and loads on nodes. And other tools of finite element analysis can be in use as a black box to interface with this program.

• Structural Engineering and Mechanics
• /
• v.80 no.5
• /
• pp.553-562
• /
• 2021

In this paper, a hybrid scaled boundary finite element and finite volume method (SBFEM-FVM) is proposed for simulating hydraulic-fracture propagation in brittle concrete materials. As a semi-analytical method, the scaled boundary finite element method is introduced for modelling concrete crack propagation under both an external force and water pressure. The finite volume method is employed to model the water within the crack and consider the relationship between the water pressure and the crack opening distance. The cohesive crack model is used to analyse the non-linear fracture process zone. The numerical results are compared with experimental data, indicating that the F-CMOD curves and water pressure changes under different loading conditions are approximately the same. Different types of water pressure distributions are also studied with the proposed coupled model, and the results show that the internal water pressure distribution has an important influence on crack propagation.

• Steel and Composite Structures
• /
• v.10 no.6
• /
• pp.541-561
• /
• 2010

This paper presents finite element analyses, experimental measurements and finite element model updating of an arch type steel laboratory bridge model using semi-rigid connections. The laboratory bridge model is a single span and fixed base structure with a length of 6.1 m and width of 1.1m. The height of the bridge column is 0.85 m and the maximum arch height is 0.95 m. Firstly, a finite element model of the bridge is created in SAP2000 program and analytical dynamic characteristics such as natural frequencies and mode shapes are determined. Then, experimental measurements using ambient vibration tests are performed and dynamic characteristics (natural frequencies, mode shapes and damping ratios) are obtained. Ambient vibration tests are performed under natural excitations such as wind and small impact effects. The Enhanced Frequency Domain Decomposition method in the frequency domain and the Stochastic Subspace Identification method in the time domain are used to extract the dynamic characteristics. Then the finite element model of the bridge is updated using linear elastic rotational springs in the supports and structural element connections to minimize the differences between analytically and experimentally estimated dynamic characteristics. At the end of the study, maximum differences in the natural frequencies are reduced on average from 47% to 2.6%. It is seen that there is a good agreement between analytical and experimental results after finite element model updating. Also, connection percentages of the all structural elements to joints are determined depending on the rotational spring stiffness.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.1
• /
• pp.88-99
• /
• 1996

The explicit scheme for finite element analysis of sheet metal forming problems has been widely used for providing practical solutions since it improves the convergency problem, memory size and computational time especially for the case of complicated geometry and large element number. The explicit schemes in general use are based on the elastic-plastic modeling of material requiring large computataion time. In the present work, a basic formulation for rigid-plastic explicit finite element analysis of plain strain sheet metal forming problems has been proposed. The effect of some basic parameters involved in the dynamic analysis has been studied in detail. Thus, the effective ranges of parameters have been proposed for numerical simultion by the rigid-plastic explicit finite element method. A direct trial-and-error method is introduced to treat contact and friction. In computation, sheet material is assumed to possess normal anisotropy and rigid-plastic workhardening characteristics. In order to show the validity and effectiveness of the proposed explicit scheme, computations are carried out for cylindrical punch stretching and the computational results are compared with those by the implicit scheme as well as with a commercial code. The proposed rigid-plastic exlicit finite element method can be used as a robust and efficient computational method for analysis of sheet metal forming.

• 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.

• Transactions of Materials Processing
• /
• v.16 no.3 s.93
• /
• pp.187-192
• /
• 2007

The finite volume method for forging simulation is examined to reveal its possibility as well as its problem in this paper. For this study, the finite volume method based MSC/SuperForge and the finite element method based AFDEX are employed. The simulated results of the homogeneous compression obtained by the two softwares are compared to indicate the problems of the finite volume method while several application examples are given to show the possibility of the finite volume method fur simulation of complex hot forging processes. It is shown that the finite volume method can not predict the exact solution of the homogeneous compression especially in terms of forming load and deformed shape but that it is helpful to simulate very complex forging processes which can hardly be simulated by the conventional finite element method.

• Structural Engineering and Mechanics
• /
• v.11 no.6
• /
• pp.591-604
• /
• 2001

In this paper, a finite element with embedded displacement discontinuity which eliminates the need for remeshing of elements in the discrete crack approach is applied for the progressive fracture analysis of concrete structures. A finite element formulation is implemented with the extension of the principle of virtual work to a continuum which contains internal displacement discontinuity. By introducing a discontinuous displacement shape function into the finite element formulation, the displacement discontinuity is obtained within an element. By applying either a nonlinear or an idealized linear softening curve representing the fracture process zone (FPZ) of concrete as a constitutive equation to the displacement discontinuity, progressive fracture analysis of concrete structures is performed. In this analysis, localized progressive fracture simultaneous with crack closure in concrete structures under mixed mode loading is simulated by adopting the unloading path in the softening curve. Several examples demonstrate the capability of the analytical technique for the progressive fracture analysis of concrete structures.