• Composites Research
• /
• v.23 no.4
• /
• pp.7-13
• /
• 2010

A volume integral equation method (VIEM) is introduced for the solution of elastostatic problems in an unbounded isotropic elastic solids containing interacting multiple anisotropic inclusions subject to remote uniaxial tension. The method is applied to two-dimensional problems involving long parallel cylindrical inclusions. A detailed analysis of stress field at the interface between the matrix and the central inclusion is carried out for square and hexagonal packing of the inclusions. Effects of the number of anisotropic inclusions and various fiber volume fractions on the stress field at the interface between the matrix and the central inclusion are also investigated in detail. The accuracy of the method is validated by solving the single inclusion problem for which solutions are available in the literature.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.65-73
• /
• 1989

Continuum formulations for the expressions of dynamic energy release rates and computational methods for dynamic stress intensity factors are developed for the analysis of dynamic fracture problems subjected to stress wave loading. Explicit volume integral expressions for instantaneous dynamic energy release rates are derived by modeling virtual crack extensions with the dynamic Eulerian-Lagrangian kinematic description. In the finite element applications a finite region around a crack-tip is modeled by using quarter-point singular isoparametric elements, and the volume integrals are evaluated for each crack-tip element during virtual crack extensions while the singularity is maintained. It is shown that the use of the present method is more reliable and accurate for the dynamic fracture analysis than that of other path-independent integral methods when the effects of stress waves are significant.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.4
• /
• pp.437-442
• /
• 2012

To analyze a multibody system, this paper proposes an implicit numerical integration method for joint coordinates subsystem synthesis method. To verify the proposed method, a multibody model for an unmanned robot vehicle, which consists of six identical independent suspension systems, is developed. The symbolic method is applied to compute the system Jacobian matrix for the implicit integration method. The proposed method is also verified by performing rough terrain run-over simulation in comparison with the conventional implicit integration method. In addition, to evaluate the efficiency of the proposed method, the CPU time obtained by using this method is compared with that obtained by using the conventional implicit method.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.13 no.9
• /
• pp.937-946
• /
• 2002

In this paper, we present a stable solution of the transient electromagnetic scattering from the conducting objects. This method does not utilize the conventional marching-on in time (MOT) solution. Instead we solve the time domain integral equation by expressing the transient behavior of the induced current in terms of weighted Laguerre polynomials. By using this basis functions for the temporal variation, the time derivative in the integral equation can be handled analytically. Since these temporal basis functions converge to zero as time progresses, the transient response of the induced current does not have a late time oscillation. To show the validity of the proposed method, we solve a time domain electric feld integral equation and compare the results of MOT, Mie solution, and the inverse discrete Fourier transform (IDFT) of the solution obtained in the frequency domain.

• Proceedings of the Acoustical Society of Korea Conference
• /
• 1998.06c
• /
• pp.309.1-312
• /
• 1998

일반적으로 음향 문제에 상용되는 경계요소법은 Kirchhoff-Helmholtz 적분 방정식에 약특이성과 강특이성의 커널을 갖고 있어, 경계면에 매우 근접한 음장을 해석할 때 수치 적분 과정에서 큰 오차를 유발한다. 본 연구에서는 평면파 성분을 이용하여 약특이성 방정식 및 특이성이 제거된 음장 음압의 과도한 오차는 약특이성 경계 적분 방정식의 적용으로 제거될 수 있었다. 부드러운 경계면을 가진 경우는 모든 특이성의 제거가 가능하여 특이성 처리를 위한 특별한 처리가 불필요하게 되었다. 제안된 방법을 검증하기 위하여 몇 가지 단순한 모델에 대하여 경계 요소 계산을 수행하였고, 경계면 부근의 근접 음장에서 음압 예측의 정확도가 향상되는 결과를 얻었다.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.06a
• /
• pp.340-345
• /
• 2000

진동하는 구조물의 음향 방사 예측에는 키르히호프-헬름홀쯔 적분 방정식에 근본을 둔 경계 요소 해석이 널리 사용된다. 이 경계 요소 해석은 익히 알고 있듯이 구조물의 동적 거동이 정량적으로 표현될 수 있는 경우는 매우 높은 정확도의 예측 결과를 제공한다. 그러나 실제 현상에서 접할 수 있는 복잡한 구조물의 음향 방사 예측에는 많은 변수들로 인해 예측의 정확도가 감소됨은 확실하다. 다른 방법으로는 실험을 통한 임의의 음장 예측 방법인 근음장 음향 홀로그래피(nearfield acoustical holography) 방법을 들 수 있다. 이 방법은 실제로 발생되는 음향 방사로부터 마이크로폰을 이용하여 홀로그램면의 음압 또는 속도를 측정하고 키르히호프-헬름홀쯔 적분 방정식에 적용하여 임의의 홀로그램면에 투사(mapping)시켜 음장을 예측하는 방법이다. 근음장 음향 홀로그래피는 탁월한 정확성을 갖고 있으나, 측정의 복잡성과 홀로그램면을 형성하기 위한 많은 이산점(절점)의 필요성 등의 단점을 갖고 있다. 본 논문에서는 또 다른 음장 예측 방법인 실험의 장점과 유한 요소 해석의 장정을 복합시킨 모드 확장 방법(modal expansion method)을 사용하여 단순 구조물인 평판의 진동에 의한 음장을 예측해 보았다. 모드 확장 방법은 구조물의 동적 거동은 모드의 선형 조합으로 표현될 수 있다는 것에 그 원리를 둔다. 본 논문은 단순 평판을 대상으로 유한 요소 해석으로 구한 모드 정보와 실험에 의해 얻은 입의 가진 주파수에 대한 진동 표면의 속도 분포를 조합하여 속도 경계 조건을 구성, 경계 요소 해석으로 음장 예측을 수행하였으며 모드 확장 방법을 사용함에 있어 고려해야할 몇 가지 사항에 대해 다루었다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.1
• /
• pp.72-80
• /
• 1996

Two analytical solutions for oscillations in a rectangular harbour are presented. In this paper, the correct solution is obtained by use of matched asymptotic expansion method, which was first derived by Mei(1989). The other solution derived from eigenfunction expansion method is also presented, in which more accurate numerical integration is employed. In order to check the solutions, amplification factors inside the harbor are calculated and plotted by both analytical methods and numerical boundary integral equation method.

• Journal of Korean Society of Disaster and Security
• /
• v.6 no.3
• /
• pp.1-8
• /
• 2013

Based on random vibration theory, a procedure for calculating the dynamic response of the tall building to time-dependent random excitation is developed. In this paper, the fluctuating along- wind load is assumed as time-dependent random process described by the time-independent random process with deterministic function during a short duration of time. By deterministic function A(t)=1-exp($-{\beta}t$), the absolute value square of oscillatory function is represented from author's studies. The time-dependent random response spectral density is represented by using the absolute value square of oscillatory function and equivalent wind load spectrum of Solari. Especially, dynamic mean square response of the tall building subjected to fluctuating wind loads was derived as analysis function by the Cauchy's Integral Formula and Residue Theorem. As analysis examples, there were compared the numerical integral analytic results with the analysis fun. results by dynamic properties of the tall uilding.

• Composites Research
• /
• v.26 no.4
• /
• pp.235-244
• /
• 2013

A volume integral equation method (VIEM) is introduced for the solution of elastostatic problems in an unbounded isotropic elastic solid containing multiple isotropic elliptical fibers of arbitrary orientation subject to uniform stress at infinity. The fibers are assumed to be long parallel elliptical cylinders composed of isotropic elastic material perfectly bonded to the isotropic matrix. The solid is assumed to be under plane strain on the plane normal to the cylinders. A detailed analysis of the stress field at the matrix-fiber interface for square and hexagonal packing of the fibers is carried out for different values of the number, orientation angles and concentration of the elliptical fibers. The accuracy and efficiency of the method are examined through comparison with results obtained from analytical and finite element methods.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.5
• /
• pp.483-491
• /
• 2010

In this paper, ${\Delta}J$-integral is used to analyze fatigue crack growth of viscoelastic material. Using analytical integral value, the J-integral is calculated. So, reduction of calculation time and increase of accuracy are made possible. The stress intensity factor is calculated using the finite element method code. In difference with existed fatigue crack analysis using ${\Delta}K$, we were successfully able to analyze various load amplitude and the fatigue crack of load cycle only with two fatigue crack growth parameters and creep compliance. The analysis gives N-${\alpha}$ curve for simulation of crack growth, and the curve almost corresponds with test results.