• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1994년도 추계학술대회 논문집
• /
• pp.837-842
• /
• 1994

The boundary element method was used for studying singularities of an interface crack with contact zones. Because the contact zone size was extremely small in a tention field, a large number of Gaussian points were is used for numerical integration of the Kernels. In order to estimate the contact zone size, iteration method was used. The interface crack tips with contact zones showed no oscillatory behavior and J-integral values across the interface were conserved.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제23권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.

• Structural Engineering and Mechanics
• /
• 제80권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.

• Structural Engineering and Mechanics
• /
• 제5권6호
• /
• pp.705-713
• /
• 1997

The residual thermal stresses at the interface corner between the elastic fiber and the viscoelastic matrix of a two-dimensional unidirectional laminate due to cooling from cure temperature down to room temperature were studied. The matrix material was assumed to be thermorheologically simple. The time-domain boundary element method was employed to investigate the nature of stresses on the interface. Numerical results show that very large stress gradients are present at the interface corner and this stress singularity might lead to local yielding or fiber-matrix debonding.

• 대한조선학회논문집
• /
• 제58권5호
• /
• pp.314-321
• /
• 2021

In the existing whole ship analysis, topside was modeled as mass element. However recently, the topside is modeled as beam element due to the owner's requirement to improve the maturity of the whole ship FE model. To follow the owner'srequirement, detailed information for topside drawing and modeling, which may delay analysis schedule, is needed. However, it is hard to respond effectively to this matter due to the lack of study on the topside from the hull perspective. Therefore in this study, the effect of the topside on the hull is investigated when the topside is modeled as a mass element or beam element respectively. In addition, the interface modeling method is analyzed to verify modeling method used in the existing whole ship analysis. The results indicate that the interface and topside modeling method used in existing whole ship analysis are appropriate. This conclusion will be the technical basis for responding to owner's requirement about the topside modeling method.

• 한국농공학회논문집
• /
• 제50권4호
• /
• pp.67-75
• /
• 2008

In construction facilities such as bridges, the fluid boundary layer(or water film) is formed at the structure-soil interface by the inflow into the system due to rainfall or/and rising ground-water. As a result, the structure-soil interaction(SSI) state changes into the structure-fluid-soil interaction(SFSI) state. In general, construction facilities may be endangered by the inflow of water into the soil foundation. Thus, it is important to predict the dynamic SFSI responses accurately so that the facilities may be properly designed against such dangers. It is desired to have the robust tools of attaining such a purpose. However, there has not been any report of a method for the SFSI analyses. The objective of this study is to propose an efficient method of finite element modelling using the new interface element named hybrid interface element capable of giving reasonable predictions of the dynamic SFSI response. This element enables the simulation of the limited normal tensile resistance and the tangential hydro-plane behaviour, which has not been preceded in the previous studies. The hybrid interface element was tested numerically for its validity and employed in the analysis of SFSI responses of the continuous bridge subjected to seismic load under rainfall or/and rising ground-water condition. It showed that dynamic responses of the continuous bridge resting on direct foundation may be amplified under rainfall condition and consequently lead to significant variation of stresses.

• 한국정밀공학회지
• /
• 제16권4호통권97호
• /
• pp.197-204
• /
• 1999

The boundary element method was used for studying singularities of an interface crack with contact zones. The iterative procedure is applied to estimate the contact zone size. Because the contact zone size was extremely small in a tension field, a large number of Gaussian points were used for numerical integration of the Kernels. Stress extrapolation method and J-integral were used ofr determining stress intensity factors. When the interface crack was assumed to have opened tips, oscillatory singularities appear near the tips of the interface crack. But the interface crack with contact zone which Comninou suggested had no oscillatory behavior. The contact zone size under shear loading was much larger than that under tensile. The stress intensity factors computed by stress extrapolation method were close to those of Comninou's solution. And the stress intensity factor evaluated by J-integral was similar to that by stress extrapolation method.

• 대한수학회보
• /
• 제51권3호
• /
• pp.717-737
• /
• 2014

We study the effect of numerical quadrature in space on semidiscrete and fully discrete piecewise linear finite element methods for parabolic interface problems. Optimal $L^2(L^2)$ and $L^2(H^1)$ error estimates are shown to hold for semidiscrete problem under suitable regularity of the true solution in whole domain. Further, fully discrete scheme based on backward Euler method has also analyzed and optimal $L^2(L^2)$ norm error estimate is established. The error estimates are obtained for fitted finite element discretization based on straight interface triangles.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집A
• /
• pp.176-181
• /
• 2000

Stress intensity factors for a rigid line inclusion tying along a bimaterial interface are calculated by the boundary element method with the multiregion and double-Point techniques. The formula between the stress intensity factors and the inclusion surface stresses are derived. The numerical values of the stress intensity factors for the bimaterial interface rigid line inclusion in the infinite body are proved to be in good agreement within 3% when compared with the previous exact solutions. In the finite bimaterial systems, the stress intensity factors for the center and edge rigid line inclusions at interface are computed with the variation of the rigid line inclusion length and the shear modulus ratio under the biaxial and uniaxial loading conditions.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회A
• /
• pp.442-447
• /
• 2008

This paper presents a new approach to simulate fluid-solid interaction problems involving non-matching interfaces. The coupling between fluid and solid domains with dissimilar finite element meshes consisting of 4-node quadrilateral elements is achieved by using the interface element method (IEM). Conditions of compatibility between fluid and solid meshes are satisfied exactly by introducing the interface elements defined on interfacing regions. Importantly, a consistent transfer of loads through matching interface element meshes guarantees the present method to be an efficient approach of the solution strategy to fluid-solid interaction problems. An arbitrary Lagrangian-Eulerian (ALE) description is adopted for the fluid domain, while for the solid domain an updated Lagrangian formulation is considered to accommodate finite deformations of an elastic structure. The stabilized equal order velocity-pressure elements for incompressible flows are used in the motion of fluids. Fully coupled equations are solved simultaneously in a single computational domain. Numerical results are presented for fluid-solid interaction problems involving nonmatching interfaces to demonstrate the effectiveness of the methodology.