• Journal of the Earthquake Engineering Society of Korea
• /
• v.5 no.4
• /
• pp.27-38
• /
• 2001

A mechanical lumped parameter model is proposed for the dynamic modeling of a semi-infinite reservoir. A semi-analytic transmitting boundary is derived for a semi-infinite 2-D reservoir of constant depth. The characteristics of the solution are examined in both frequency and time domains. Mass, damping and spring coefficients of the mechanical model are obtained to preserve the major features of the solution such as eigenfrequencies and the shapes of Bessel functions that appear as kernels in the convolution integrals. The lumped parameter model in its final form consists of two masses, a spring and two dampers for each eigenfrequency. Application examples demonstrated that the new lumped parameter model could be used for the time domain analysis of dam-reservoir systems.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.121-128
• /
• 1999

This paper presents a time domain method for soil-structure interaction analysis for seismic loadings. It is based on the finite element formulation incorporating analytical frequency-dependent infinite elements for the far field soil. The dynamic stiffness matrices of the far field region formulated using the present method in frequency domain can be easily transformed into the corresponding matrices in time domain. At first, the equivalent earthquake forces are evaluated along the interface between the near and the far fields from the free-field response analysis carried out in frequency domain, and the results are transformed into the time domain. An efficient procedure is developed for the convolution integrals to evaluate the interaction force along the interface, which depends on the response on the interface at the past time instances as well as the concurrent instance. Then, the dynamic responses are obtained for the equivalent earthquake force and the interaction force using Newmark direct integration technique. Since the response analysis is carried out in time domain, it can be easily extended to the nonlinear analysis. Example analysis has been carried out to verify the present method in a multi-layered half-space.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.167-171
• /
• 2004

The Spar platform with deep draft is characterized as effective structure in extreme wave condition, which has larger natural period than that of waves in sea. In this paper, the time simulation of motion responses of Spar with catenary mooring line is presented in irregular waves. The memory effect is modeled by added mass at infinite frequency and convolution integrals in terms of wave damping coefficients. The added mass, wave damping coefficients and wave exciting forces are obtained from three-dimensional panel method in the frequency domain. The motion equations are consisted of forces for inetia, memory effect, hydrostatic restoring, wave exciting and mooring line. The forces of mooring line are modeled as quasi-static catenary cable.

• Communications of the Korean Mathematical Society
• /
• v.27 no.4
• /
• pp.721-732
• /
• 2012

The principal aim of the paper is to investigate new integral expression $${\int}_0^{\infty}x^{s+1}e^{-{\sigma}x^2}L_m^{(\gamma,\delta)}\;({\zeta};{\sigma}x^2)\;L_n^{(\alpha,\beta)}\;({\xi};{\sigma}x^2)\;J_s\;(xy)\;dx$$, where $y$ is a positive real number; $\sigma$, $\zeta$ and $\xi$ are complex numbers with positive real parts; $s$, $\alpha$, $\beta$, $\gamma$ and $\delta$ are complex numbers whose real parts are greater than -1; $J_n(x)$ is Bessel function and $L_n^{(\alpha,\beta)}$ (${\gamma};x$) is generalized Laguerre polynomials. Some integral formulas have been obtained. The Maple implementation has also been examined.

• Structural Engineering and Mechanics
• /
• v.2 no.3
• /
• pp.245-256
• /
• 1994

A new axisymmetric crack model is proposed on the basis of p-version of the finite element method limited to theory of small scale yielding. To this end, axisymmetric stress element is formulated by integrals of Legendre polynomial which has hierarchical nature and orthogonality relationship. The virtual crack extension method has been adopted to calculate the stress intensity factors for 3-D axisymmetric cracked bodies where the potential energy change as a function of position along the crack front is calculated. The sensitivity with respect to the aspect ratio and Poisson locking has been tested to ascertain the robustness of p-version axisymmetric element. Also, the limit value that is an exact solution obtained by FEM when degree of freedom is infinite can be estimated using the extrapolation equation based on error prediction in energy norm. Numerical examples of thick-walled cylinder, axisymmetric crack in a round bar and internal part-thorough cracked pipes are tested with high precision.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.5
• /
• pp.1-7
• /
• 2004

The Spar platform with deep draft is characterized as effective structure in extreme wave condition, which has larger natural period than that of waves in sea. In this paper, the time domain simulation of motion responses of Spar with catenary mooring line is presented in irregular waves. The memory effect is modeled by added mass at infinite frequency and convolution integrals in terms of wave damping coefficients. The added mass, wave damping coefficients and wave exciting forces are obtained from three-dimensional panel method in the frequency domain. The motion equations are consisted of forces for inertia, memory effect, hydrostatic restoring, wave exciting and mooring line. The forces of mooring line are modeled as quasi-static catenary cable.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.1 no.4
• /
• pp.59-68
• /
• 1997

In analysis of underground structures, the effects of artificial boundary conditions are considered as one of the major reasons for differences from experimental results. These phenomena can be overcome by using the boundary elements which satisfy the multi-layered half space conditions. The fundamental solutions of multi-layered half-space for boundary element method is formulated satisfying the transmission and reflection of waves at each layer interface and radiation conditions at bottom layer. The governing equations can be obtained from the displacements at each layer which are expressed in terms of harmonic functions. All types of waves can be included using the complete response from semi-infinite integrals with respect to horizontal wavenumbers using expansion of Fourier series and Hankel transformation. Two dimensional Green's functions are derived from cylindrical Navier equations and potentials performing infinite integration in y-direction. In this case, it is effective to transform into two dimensional problem using semi-analytical integration and sinusoidal Bessel function. Some verifications are given to show the accuracy and efficiency of the developed method, and numerical examples to demonstrate the dynamic behavior of underground with various properties.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.15 no.5
• /
• pp.499-508
• /
• 2004

An efficient method of moments(MoM), which can lead to the analytical evaluation of the matrix elements, is proposed to analyze microstrip structures. The present method is formulated in conjunction with use of a new closed-form spatial-domain Green's functions which are derived by use of the integral formula for semi-infinite integrals of Bessel functions. It is observed that the computational efficiency such as the amount of calculation and computation speed has been improved due to the present MoM scheme by a factor of about 4 in comparison with the previous method. To validate the proposed method, several numerical examples are presented.

• Structural Engineering and Mechanics
• /
• v.21 no.4
• /
• pp.437-450
• /
• 2005

In this study, a complete analysis of soil-structure interaction problems is presented which includes a modelling of the near surrounding of the building (near-field) and a special description of the wave propagation process in larger distances (far-field). In order to reduce the computational effort which can be very high for time domain analysis of wave propagation problems, a special approach based on similarity transformation of the infinite domain on the near-field/far-field interface is applied for the wave radiation of the far-field. The near-field is discretised with standard Finite Elements, which also allows to introduce non-linear material behaviour. In this paper, a new approach to calculate the involved convolution integrals is presented. This approximation in time leads to a dramatically reduced computational effort for long simulation times, while the accuracy of the method is not affected. Finally, some benchmark examples are presented, which are compared to a coupled Finite Element/Boundary Element approach. The results are in excellent agreement with those of the coupled Finite Element/Boundary Element procedure, while the accuracy is not reduced. Furthermore, the presented approach is easy to incorporate in any Finite Element code, so the practical relevance is high.

• JSTS:Journal of Semiconductor Technology and Science
• /
• v.10 no.3
• /
• pp.203-212
• /
• 2010

In this paper, we present an analytical model for the first eigen energy level ($E_0$) of the carriers in the inversion layer in present generation MOSFETs, having ultrathin gate oxides and high substrate doping concentrations. Commonly used approaches to evaluate $E_0$ make either or both of the following two assumptions: one is that the barrier height at the oxide-semiconductor interface is infinite (with the consequence that the wave function at this interface is forced to zero), while the other is the triangular potential well approximation within the semiconductor (resulting in a constant electric field throughout the semiconductor, equal to the surface electric field). Obviously, both these assumptions are wrong, however, in order to correctly account for these two effects, one needs to solve Schrodinger and Poisson equations simultaneously, with the approach turning numerical and computationally intensive. In this work, we have derived a closed-form analytical expression for $E_0$, with due considerations for both the assumptions mentioned above. In order to account for the finite barrier height at the oxide-semiconductor interface, we have used the asymptotic approximations of the Airy function integrals to find the wave functions at the oxide and the semiconductor. Then, by applying the boundary condition at the oxide-semiconductor interface, we developed the model for $E_0$. With regard to the second assumption, we proposed the inclusion of a fitting parameter in the wellknown effective electric field model. The results matched very well with those obtained from Li's model. Another unique contribution of this work is to explicitly account for the finite oxide-semiconductor barrier height, which none of the reported works considered.