• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.8
• /
• pp.1429-1437
• /
• 1992

A new finite element formulation based on the relaxation type hereditary integral is presented for a time-domain analysis of isotropic, linear viscoelastic problems. The semi-discrete variational approximation and elastic-viscoelastic correspondence principle are used in the theoretical development of the proposed method. In a time-stepping procedure of final, linear algebraic system equations, only a small additional computation for past history is required since the equivalent stiffness matrix is constant. The viscoelasticity matrices are derived and the stress computation algorithm is given in matrix form. The effect of time increment and Gauss point numbers on the numerical accuracy is examined. Two dimensional numerical examples of plane strain and plane stress are solved and compared with the analytical solutions to demonstrate the versatility and accuracy of the present method.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.3
• /
• pp.166-177
• /
• 2005

This study developed the jet integral model to analyze the behavior of the wastewater discharge in the near field using the fourth order Runge-Kutta method in order to numerically solve the problems of six ordinary differential equations and six unknowns. This jet integral model used the entrainment hypothesis and the manipulation of sonle shape constant. This study also conducted the hydraulic experiments fnr single horizontal buoyant Jet using LIF through the calibration procedure. The results calculated by the previous models, CORMIX 1 and VISJET, and the proposed jet integral model were compared to the hydraulic experimental results. The centerline trajectories predicted by the proposed model were in good agreements with the experimental results in the transition region whereas the trajectories calculated by the VISJET model agreed well with the measured data in the momentum and buoyancy-dominated regions. The centerline dilution calculated by the proposed model agreed generally with the measured dilution in the intial and transition regions while the centerline dilution predicted by the CORMIX 1 was in good agreements with the experimental results in the momentum and buoyancy-dominated regions.

• Economic and Environmental Geology
• /
• v.27 no.2
• /
• pp.191-199
• /
• 1994

We have developed an algorithm based on the method of integral equations to simulate the magnetotelluric (MT) responses of three-dimensional (3-D) bodies in a layered half-space. The inhomogeneities are divided into a number of cells and are replaced by an equivalent current distribution which is approximated by pulse basis functions. A matrix equation is constructed using the electric Green's tensor function appropriate to a layered earth, and is solved for the vector current in each cell. Subsequently, scattered fields are found by integrating electric and magnetic Green's tensor functions over the scattering current About a 3-D conductive body near the earth's surface, interpretation using 2-D transverse electric modeling schemes can imply highly erratic low resistivities at depth. This is why these routines do not account for the effect of boundary charges. However, centrally located profiles across elongate 3-D prisms may be modeled accurately with a 2-D transverse magnetic algorithm, which implicitly includes boundary charges in its formulation. Multifrequency calculations show that apparent resistivity and impedance phase are really two complementary parameters. Hence, they should be treated simultaneously in broadband MT interpretation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.3
• /
• pp.267-275
• /
• 2009

This paper presents the p-convergent coupling element on the basis of the ESSE(equivalent single layer shell element) and the PLLE(partial-linear layerwise element) to analyze laminated composite plates. The ESSE is formulated by the degenerated shell theory, on the other hand, the assumption of the PLLE is piecewise linear variation of the in-plane displacement and a constant value of lateral displacement across the thickness. The proposed finite element model is based on p-convergence approach. The integrals of Legendre polynomials and Gauss-Lobatto technique are chosen to interpolate displacement fields and to implement numerical quadrature, respectively. This study has been focused on the verification of p-convergent element. For this purpose, various finite element multiple models associated with the combination of ESSE and PLLE elements are tested to show numerical stability. The simple examples such as a cantilever beam subjected vertical load and a plate with tension are adopted to evaluate the performance of proposed element.

• Journal of Korea Water Resources Association
• /
• v.50 no.9
• /
• pp.637-646
• /
• 2017

This study presents a numerical model for simulating dense interflows. The governing equations are provided and the finite difference method is used with the $k-{\varepsilon}$ turbulence model. The model is used to simulate a dense interflow established in a deep ambient water, resulting velocity and excess density profiles. It is observed that velocity decreases in the longitudinal direction due to water entrainment in the vicinity of the outlet and rarely changes for increased Richardson number. Similarity collapses of velocity and excess density are obtained, but those of turbulent kinetic energy and dissipation rate are not. A shape factor for the dense interflow is obtained from the simulated profiles. The value of this shape factor can be used in the layer-averaged modeling of dense interflows. In addition, a buoyancy-related parameter ($c_{3{\varepsilon}}$) for the $k-{\varepsilon}$ model and the volume expansion coefficient (${\beta}_0$) are obtained from the simulated results.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2013.10a
• /
• pp.206-213
• /
• 2013

Transient response of a crack in a functionally graded piezoelectric material (FGPM) interface layer between two dissimilar homogeneous piezoelectric layers under anti-plane shear is analyzed using integral transform approaches. The properties of the FGPM layer vary continuously along the thickness. Laplace and Fourier transforms are used to reduce the problem to two sets of dual integral equations, which are then expressed to the Fredholm integral equations of the second kind. Numerical values on the dynamic energy release rate (DERR) are presented for the FGPM to show the effects on electric loading, gradient of the material properties, and thickness of the layers. Computed results yield following conclusions: (a) the DERR increases with the increase of the gradient of the material properties of the FGPM layer; (b) certain direction and magnitude of the electric impact loading impedes crack extension; (c) increase of the thickness of the FGPM layer and the homogeneous piezoelectric layer which has larger material properties than those of the crack plane are beneficial to increase of the resistance of transient fracture of the FGPM layer.

• Economic and Environmental Geology
• /
• v.27 no.1
• /
• pp.41-47
• /
• 1994

A fast Hankel transform (FHT) algorithm (Anderson, 1982) is applied to numerical evaluation of many Green's tensor integrals encountered in three-dimensional electromagnetic modeling using integral equations. Efficient computation of Hankel transforms is obtained by a combination of related and lagged convolutions which are available in the FHT. We express Green's tensor integrals for a layered half-space, and rewrite those to a form of related functions so that the FHT can be applied in an efficient manner. By use of the FHT, a complete or full matrix of the related Hankel transform can be rapidly and accurately calculated for about the same computation time as would be required for a single direct convolution. Computing time for a five-layer half-space shows that the FHT is about 117 and 4 times faster than conventional direct and multiple lagged convolution methods, respectively.

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

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

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.7
• /
• pp.881-889
• /
• 2010

A volume integral equation method (VIEM) is introduced for solving the elastostatic problems related to an unbounded isotropic elastic solid; this solid is subjected to remote uniaxial tension, and it contains multiple interacting isotropic inclusions. The method is applied to two-dimensional problems involving long parallel cylindrical inclusions. A detailed analysis of the stress field at the interface between the matrix and the central inclusion is carried out; square and hexagonal packing of the inclusions are considered. The effects of the number of isotropic inclusions and different fiber volume fractions on the stress field at the interface between the matrix and the central inclusion are also investigated in detail. The accuracy and efficiency of the method are clarified by comparing the results obtained by analytical and finite element methods. The VIEM is shown to be very accurate and effective for investigating the local stresses in composites containing isotropic fibers.