Earticle

In this paper, the features of fatigue crack growth under ModeⅠ loading are studied. Double cantilever beam specimens were used. The effects of adherend thickness, rubber modification and adhesive thickness on fatigue crack growth were examined. The experimental results show that some of these parameters do apparently affect fatigue crack growth. Resistance to ModeⅠ fatigue crack growth are increased by rubber modification. The effects of adhesive thickness and rubber content on fatigue crack growth were explained by von Mises's equivalent stress using BEM analysis.

The present study is to investigate the effect of wave-structure interaction such as wave oscillation. The theoretical method is based upon the linear diffraction theory obtained by the boundary element method. The water depth and incident wave period in fluid region are assumed to be constant. To investigate the wave interaction due to offshore structures, the numerical program has been developed and the simulation has been carried out by varying the conditions of distance and width of offshore structures. This study can effectively be utilized for safety assessment to various breakwater systems and layout of offshore breakwater in the ocean and coastal field. It can give information for the safety to construct offshore structure and revetment in coastal region.

2차원 경계요소-유한요소 연계법을 사용하여 구형 수조구조물의 동적거동 특성을 파악하였다. 비점성, 비압축성 이상유체의 운동은 경계요소법으로 모델링하였으며 구조물의 운동은 유한요소법으로 모델링하였다. 경계요소법의 적용시 Singularity를 소거한 경계적분식을 사용하였다. 유체와 구조물의 경계면에서 적합조건과 평형조건을 만족시킴으로써 경계요소와 유한요소를 연계하였다. 유체-구조물 상호작용의 영향은 유체의 부가질량행렬과 유체 유동에 의한 강성행렬로 연계된 유체-구조물계에 반영된다. 연계된 유체-구조물계의 고유치 문제로부터 수조구조물의 고유진동수외 고유모드를 구하였으며 수평 및 수직 자바운동에 의하여 수조구조물에 발생하는 유체동압력과 자유수면의 유동도 산출하였다.

Dynamic behavior of flexible rectangular liquid containers is analyzed by a two-dimensional coupled boundary element-finite element method. The irrotational motion of inviscid and incompressible ideal fluid is modeled by boundary elements and the motion of structure by finite elements. A singularity free integral formulation is employed for the implementation of boundary element method. Coupling is performed by using compatibility and equilibrium conditions along the interface between the fluid and structure. The fluid-structure interaction effects are reflected into the coupled equation of motion as added fluid mass matrix and sloshing stiffness matrix. By solving the eigen-problem for the coupled equation of motion, natural frequencies and mode shapes of coupled system are obtained. The free surface sloshing motion and hydrodynamic pressure developed in a flexible rectangular container due to horizontal and vertical ground motions are computed in time domain.

Errors included in solutions obtained by the boundary element method are generally larger than those by the finite element method in the case that the number of discreted elements is small. One of the reasons is supposed to be attributed to the error which will be produced in the numerical integration of the singular functions in two dimensional elastic problem. Then, treatment of analytical integration to reduce computing time and to decrease errors of boundary element method are proposed.

해안구조물에 작용하는 설계파압(設計波壓)을 수치적으로 계산하기 위해서는 먼저 쇄파한계(碎波限界)에 가까운 큰 중복파랑을 수치적으로 재현할 필요가 있다. 이를 위해서는 지배방정식(支配方程式)과 비선형항(非線形項)을 포함하는 경계조건을 효과적으로 반영해야 하며 특히, 자유표면(自由表面) 경계조건(境界條件)에서의 속도의 제곱항의 처리가 중요하다. 본 연구에서는 Newton 방법을 이용하여 제곱항을 충실히 반영하므로써 일반적인 셜계파 성향에 거의 상응하는 중복파랑을 재현하였으며 기존의 섭동법(攝動法) 또는 Fourier 전개 기법 및 수리실험 결과와 비교하여 그 정확도를 검토하였다.

To calculate the wave pressure acting on coastal structures under the design wave condition, it is often necessary to numerically reproduce the big standing wave profiles close to wave breaking condition. For this, the governing equation and all nonlinear terms occurring in boundary conditions should be effectively considered in the numerical wave profile. In particular, the velocity square term in the free surface boundary condition is very important. A boundary element method is applied here to calculate the standing wave profile with the velocity square term fully treated by Newton iterative method. In order to check the validity of the method, the numerical wave profiles are compared to ones calculated by the perturbation method, the Fourier approximation method and the hydraulic experiment.

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.

묵호항내(墨湖港內) 부진동특성(副振動特性)을 조사(調査)하기 위하여 7초부터 600초 범위(範圍)사이의 입사파(入射波)에 대하여 경계요소법(境界要素法)을 이용한 수식계산(數植計算)을 실시 하였다. 또 사용된 수식계산(數植計算) 방법(方法)은 기존(旣存)에 실시 하였던 수리모형실험(水理模型實驗) 결과(結果)를 이용하여 적용성(適用性)과 타당성(妥當性)을 검토(檢討)하였다. 부진동(副振動) 조사결과(調査結果) 11초와 3분 정도(程度)의 파가 입사(入射) 할때, 높은 증폭율(增幅率)이 나타났으며 과거의 관측보고(觀測報告)와도 유사한 결과(結果)를 나타냈다. 이로서 본(本) 수치계산(數値計算) 방법(方法)은 묵호항(墨湖港)의 정온도평가(靜穩度評價), 항만계획(港灣計劃)등에 유용(有用)하게 사용될 수 있다.

A numerical simulation using the Boundary Element Method is performed for investigating the seiche in the Mukho harbor. The range of the incident waves period is from 7 seconds to 10 minutes, which these surface waves usually may be found in cases of the swell and the long period storm surges caused by storm winds. It is found that the seiche may be occurred by resonating with the incident waves of which the periods are about 11 seconds and 3 minutes. The numerical method was verified by the hydraulic model experiment data which had performed before. Numerical results for the seiche in the Mukho harbor show good agreements with the field observations, so that this research may be useful to estimate harbor calmness in a harbor and the harbor planning.

This paper deals with the applicability of linear elastic fracture mechanics under centrifugal force. Stress intensity factors K are calculated as a function of the inclination crack of length 2a, the position at different angular velocities 1200rpm, 2400rpm and at different values of the inclination crack angle .phi. ( .phi. = 0 .deg. , 15 .deg. , 30 .deg. , 45 .deg. , 60 .deg. , 75 .deg. , 90 .deg. ) and are measured in models of rotation disks using a boundary element method. Especially, stress intensity factors $K_{l}$ and $K_{ll}$ obtained separately from the crack tip of the mixed mode, were used to further investigate the influence of $K_{l}$ and $K_{ll}$ on fracture in rotating disks. With the increase in the speed of rotation, the effect of K/ sub l/became larger where as that of $K_{ll}$ became small. For the increase in the inclination crack angle .phi. , a decrease in $K_{l}$ and an increase in $K_{ll}$ were observed.

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.

A two-dimensional program for the analysis of bimaterial inclusion has been developed using the bound- ary element method. In order to study the effects of circular inclusion on the stress field of the crack tip, numerical analysis was performed for the straight crack of finite length around the symmetric circular inclusion whose modulus of elasticity was different from that of the matrix material. In the case of inclusion whose stiffness was smaller than that of the matrix material, the stress intensity factor was found to increase as the crack enamated. The stress intensity factor was uninfluenced from the radial change in inclusion and remained constant for the stiffness equivalent to the matrix materials, where as it decreased for the inclusion with larger stiffness. For the vareation in the distance of the inclusion, a small increase in the stress intensity factor was observed for the case with small or equal stiffness compared with the matrix materials. The inclusion with larger stiffness showed a gradual decrease in the strss intensity factor as the crack emanated.

Boundary element solution method is introduced for radiation heat transfer problem with objects inside a 2-D enclosure, where shadow zone exists. Surfaces are assumed as diffuse gray in a transparent medium. Boundary integral and boundary element equations were derived from radiation transfer equation, and their theory is reviewed. Also the process of solution methods to implement the boundary element method is analysed and explained with consideration of shadow effects. BEM solution results are compared with two test problems and are found to be good agreements with the both analytic and ANSYS Fluent numerical solutions. Therefore the current BEM analysis for radiation heat transfer problem can be considered as verified, and their efficacy with engineering applicability is established as a result.

The boundary element alternating method is applied to two dimensional and axisymmetric racture mechanics problems. Different crack depths are considered with particular emphasis on cracks. Several methods of calculating the stress intensity factor are presented and applied to several practical examples to demonstrate their accuracy.

In this paper, a new algorithm of coupling Element-Free Galerkin Method(EFGM) and Boundary Element Method(BEM) using the variational formulation is presented. A global variational coupling formulation of EFGM-BEM is achieved by combining the variational form on each subregion. In the formulation, Lagrange multiplier method is introduced to satisfy the compatibility conditions between EFGM subregion and BEM subregion. Some numerical examples are studied to verify accuracy and efficiency of the proposed method, in which numerical performance of the method is compared with that of conventional method such as EFGM-BEM direct coupling method, EFGM and BEM. The proposed method incorporating the merits of EFGM and BEM is expected to be applied to special engineering problems such as the crack propogation problems in very large domain, and underground structures with joints.

In this study, a framework for coupling of the convolution quadrature time-domain boundary element method (CQBEM) and image-based finite element method (IMFEM) is presented for 2-D elastic wave propagation. This coupling method has three advantages: 1) the finite element modeling for heterogeneous areas can be performed without difficulties by using digital data for the analysis model, 2) wave propagation in an infinite domain can be calculated with high accuracy by using the CQBEM, and 3) a small time-step size can be used. In general, a small time-step size cannot be used in the classical time-domain boundary element method. However, the CQBEM used in this analysis can address a small time-step size. This makes it possible to couple the CQBEM and image-based FEM which require a small-time step size. In this study, the formulation and validation of the pro-posed method are described and confirmed by solving fundamental elastic wave scattering problems. As a numerical example, elastic wave scattering in inhomogeneous media is demonstrated using the proposed coupling method.

We consider a numerical scheme for solving a nonlinear boundary integral equation (BIE) obtained by reformulation the nonlinear boundary value problem (BVP). We give a simple alternative to the standard collocation method for the nonlinear BIE. This method consists of one conventional linear system and another coupled linear system resulting from an auxiliary BIE which is obtained by differentiating both side of the nonlinear interior integral equations. We obtain an analogue BIE through the perturbation of the fundamental solution of Laplace's equation. We procure the super-convergence of approximate solutions.

This paper treats the optimization of rotor construction in the permanent magnet synchronous motor (PMSM) for electric vehicle (EV). While the field system of PMSM has generally one magnet per pole, we replace the magnet into plural sub-magnets. The dimensions of each sub-magnet are determined by the concept of pulse width modulation (PWM). By adopting the proposed rotor construction, we can not only reduce the space harmonics of the air-gap field but also provide space for rotor bars (i.e., damper windings) around the direct-axis. From the investigation by hybrid EE-BE (coupled finite element and boundary element) method coupled with both electric circuit and motion equation, we verify that the construction is effective for practical use.