• Structural Engineering and Mechanics
• /
• v.11 no.5
• /
• pp.547-564
• /
• 2001

The numerical formulation of a two-phase interface element appropriate for porous lining system is presented. The formulation is isoparametric and can be applied both for 2-D and 3-D analysis. Biot's theory is utilized as the basis for the development of the element constitutive theory. In order to be capable of simulating the reinforcing characteristics of some geotextiles utilized as lining system, a reinforcement component has also been implemented into the formulation. By employing this specially developed interface finite element, the influence of soil consolidation on the stress distribution along the lining system of a reservoir and a landfill has been investigated.

• Proceedings of the Korean Society of Agricultural Engineers Conference
• /
• 2000.10a
• /
• pp.425-430
• /
• 2000

In the analysis of deformation in which the stiffness is greatly different between the adjacent materials, the desired results have been obtained by using the interface element method compared with those secured by the conventional method of the concept of continua. This study deals with the deformation analysis of soft foundation by the introduction of interface element. The physical conditions of interface element are divided into three categories by Mohr-Coulomb failure criterion ie. sliding, separation, and contact. Finally the accuracy of the program proposed in this paper is proved highly accurate by performing the comparison of the theoretical values numerical results of a model element with simplified boundary conditions.

• Structural Engineering and Mechanics
• /
• v.80 no.5
• /
• pp.625-643
• /
• 2021

In this paper, we present a new method to calculate interface warping functions for the analysis of beams with geometric and material discontinuities in the longitudinal direction. The classical Saint Venant torsion theory is extended to a three-dimensional domain by considering the longitudinal direction. The interface warping is calculated by considering both adjacent cross-sections of a given interface. We also propose a finite element procedure to simultaneously calculate the interface warping function and the corresponding twisting center. The calculated interface warping functions are employed in the continuum-mechanics based beam formulation to analyze arbitrary shape cross-section beams with longitudinal discontinuities. Compared to the previous work by Yoon and Lee (2014a), both geometric and material discontinuities are considered with fewer degrees of freedom and higher accuracy in beam finite element analysis. Through various numerical examples, the effectiveness of the proposed interface warping function is demonstrated.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.1
• /
• pp.36-43
• /
• 2015

Since the crack generation and its propagation caused by welding defects is one of the main hull damage patterns, the simulation of crack propagation process has an important significance for ship safety. Based on interface element method, a finite element model to simulate crack propagation is studied in the paper. A Lennard-Jones type potential function is employed to define potential energy of the interface element. Tensile tests of steel flat plates with initial central crack are carried out. Surface energy density and spring critical stress that are suitable for the simulation of crack propagation are determined by comparing numerical calculation and tests results. Based on a large number of simulation results, the curve of simulation correction parameter plotted against the crack length is calculated.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.5 no.6
• /
• pp.141-147
• /
• 1997

For the effective finite element analysis of the structures including material interfaces or contact surfaces, interface elements are proposed. Most of early works in this problem require not only iterative computation but also complex formulation because of the kinematic nonlinearities caused from the discontinuous behavior and the stress concentration phenomena. The proposed elements, however, are consistently formulated using relative displacements and tractions between top and bottom regular finite elements. The effectiveness of these elements are shown by solving various numerical sample problems including an leaf spring and comparing with results of general finite element analysis. As a result, more stable solutions are conveniently obtaines using interface elements than regular finite elements.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.65-92
• /
• 2019

We survey a recently developed immersed finite element method (IFEM) for the interface problems. The IFEM uses structured grids such as uniform grids, even if the interface is a smooth curve. Instead of fitting the curved interface, the bases are modified so that they satisfy the jump conditions along the interface. The early versions of IFEM [1, 2] were suboptimal in convergence order [3]. Later, the consistency terms were added to the bilinear forms [4, 5], thus the scheme became optimal and the error estimates were proven. For elasticity problems with interfaces, we modify the Crouzeix-Raviart based element to satisfy the traction conditions along the interface [6], but the consistency terms are not needed. To satisfy the Korn's inequality, we add the stabilizing terms to the bilinear form. The optimal error estimate was shown for a triangular grid. Lastly, we describe the multigrid algorithms for the discretized system arising from IFEM. The prolongation operators are designed so that the prolongated function satisfy the flux continuity condition along the interface. The W-cycle convergence was proved, and the number of V-cycle is independent of the mesh size.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2230-2241
• /
• 2006

This paper presents a sole application of boundary element method to the conjugate heat transfer problem of thermally developing laminar flow in a thick walled pipe when the fluid velocities are fully developed. Due to the coupled mechanism of heat conduction in the solid region and heat convection in the fluid region, two separate solutions in the solid and fluid regions are sought to match the solid-fluid interface continuity condition. In this method, the dual reciprocity boundary element method (DRBEM) with the axial direction marching scheme is used to solve the heat convection problem and the conventional boundary element method (BEM) of axisymmetric model is applied to solve the heat conduction problem. An iterative and numerically stable BEM solution algorithm is presented, which uses the coupled interface conditions explicitly instead of uncoupled conditions. Both the local convective heat transfer coefficient at solid-fluid interface and the local mean fluid temperature are initially guessed and updated as the unknown interface thermal conditions in the iterative solution procedure. Two examples imposing uniform temperature and heat flux boundary conditions are tested in thermally developing region and compared with analytic solutions where available. The benchmark test results are shown to be in good agreement with the analytic solutions for both examples with different boundary conditions.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.10
• /
• pp.2398-2406
• /
• 1993

Abdi's numerical method(ref.13) for representing a stress singularity by shifting the mid-side nodes of isoparametric elements is reviewed. A simple technique to obtain the optimal position of the mid-side nodes in quadratic isoparametric finite element is presented. From this technique we can directly obtain the position of the side-nodes adjacent to the crack tip. It is also observed that the present technique provides good accuracy for the expression of the opening displacement and the determination of the mid-side nodes for more wide range of material properties than that obtained by Abdicant the finite element method is applied to determine stress intensity factors for pressurized crack perpendicular to and terminating at the interface of two bonded dissimilar materials. A proper definition for stress intensity factors of a crack perpendicular to bimaterial interface is provided. It is based upon a near-tip displacement solutions on the crack surface for interface crack between two dissimilar materials. Numerical testing is carried out with the eight-node and six-node elements. The results obtained are compared with the previous solutions.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.1
• /
• pp.101-112
• /
• 1994

The measurement of residual stress distribution of $Si_3N_4/SUS304$ joint was performed on 23 specimens with the same joint condition using PSPC type X-ray stress measurement system and the two-dimensional elastoplastic analysis using finite element method was also attempted. As results, residual stress distribution near the interface on the ceramic side of the joint was revealed quantitatively. Residual stress on the ceramic side of the joint was turned out to be tensional near the interface, maximum along the edge, varying in accordance with the condition of the joint and variance to be most conspicuous for the residual stress normal to the interface characterized by the stress singularities. In the vicinity of the interface, the high stress concentration occurs and residual stress distributes three-dimensionally. Therefore, the measured stress distribution differed remarkably from the result of the two-dimensional finite-element analysis. Especially at the center of the specimen near the interface, the residual stress, $\sigma_{x}$ obtained from the finite element analysis was compressive, whereas measurement using X-ray yielded tensile $\sigma_{x}$. Here we discuss two dimensional superposition model the discrepancy between the results from the two dimensional finite element analysis and X-ray measurement.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.27 no.4
• /
• pp.163-170
• /
• 2023

In the fluid-structure interaction analysis, the finite element formulation is performed for the wave equation for dynamic fluid pressure, and the dynamic pressure is defined as a degree of freedom at the fluid nodes. Therefore, to connect the fluid to the structure, it is necessary to connect the degree of freedom of fluid dynamic pressure and the degree of freedom of structure displacement through an interface element derived from the relationship between dynamic pressure and displacement. The previously proposed fluid-structure interface elements use conformal finite element meshes in which the fluid and structure match. However, it is challenging to construct conformal meshes when complex models, such as water purification plants and wastewater treatment facilities, are models. Therefore, to increase modeling convenience, a method is required to model the fluid and structure domains by independent finite element meshes and then connect them. In this study, two fluid-structure interface elements, one based on constraints and the other based on the integration of nonsmooth functions, are proposed in nonconformal finite element meshes for structures and fluids, and their accuracy is verified.