• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.89-98
• /
• 1993

As the construction of nuclear power plants are increased, nuclear wastes disposal has been faced as a serious problem. If nuclear wastes are to be buried in the underground stratum, thermo-mechanical behavior of stratum must be analyzed, because high temperature distribution has a significant effect on tunnel and surrounding stratum. In this study, in order to analyze the structural behavior of the underground which is subject to concentrated heat sources, a coupling method of nonlinear finite elements and linear boundary elements is proposed. The nonlinear finite elements (NFE) are applied in the vicinity of nuclear depository where thermo-mechanical stress is concentrated. The boundary elements are also used in infinite domain where linear behavior is expected. Using the similar method as for the problem in mechanical field, the coupled nonlinear finite element-boundary element (NFEBE) is developed. It is found that NFEBE method is more efficient than NFE which considers nonlinearity in the whole domain for the nuclear wastes depository that is expected to exhibit local nonlinearity behavior. The effect of coefficients of the rock mass such as Poisson's ratio, elastic modulus, thermal diffusivity and thermal expansion coefficient is investigated through the developed method. As a result, it is revealed that the displacements around tunnel are largely dependent on the thermal expansion coefficients.

• The Journal of the Acoustical Society of Korea
• /
• v.20 no.6
• /
• pp.94-103
• /
• 2001

The goal of using numerical methods in this study is two-fold: to replicate a set of measured, individualized HRTFs by a computer simulation, and also to visualise the resultant sound field around the head. Two methods can be wed: the Boundary Element Method (BEM) and the Infinite-Finite Element Method (IFEM). This paper presents the results of a preliminary study carried out on a KEMAR dummy-head, the geometry of which was captured with a high accuracy 3-D laser scanner and digitiser. The scanned computer model was converted to a few valid BEM and IFEM meshes with different polygon resolutions, enabling us to optimise the simulation for different frequency ranges. The results show a good agreement between simulations and measurements of the sound pressure at the blocked ear-canal of the dummy-head. The principle of reciprocity provides an effect method to simulate HRTF database. The BEM was also used to investigate the total sound field around the head, providing a tool to visualise the sound field for different arrangements of virtual acoustic imaging systems.

• Tunnel and Underground Space
• /
• v.4 no.1
• /
• pp.17-23
• /
• 1994

This study consists of two procedures on back analysis and forward analysis which is a basic tool of the former. For a safe and economical construction of underground structures, it is required to identify the structural parameters and analyze the structural behavior as exactly as possible. In this paper, a boundary element method to analyze the behavior of multi-alyered underground structures is studied, in which body forces and initial stresses are considered. That is, each layer is discritized into subregions using infinite fundamental solutions, and terms of body forces and initial stresses are transformed into boundary integral where the applied direct integral method is used. And the system of equations containing body forces and initial stresses are considered. That is, each layer is discritized into subregions using infinite fundamental solutions, and terms of body forces and initial stresses are transformed into boundary integral where the applied direct integral method is used. And the system of equations containing body forces and initial stresses are composed, then the method to solve unknowns is used with applying compatibility and equilibrium conditions between interfaces. As well, the direct search method is applied in back analysis problems. By Powell's method as a technique to search unknown parameters, assuming displacements calculated from boundary element analysis as in-situ displacements, elastic moduli and initial stresses are presumed. As consequences of this study, the results of boundary element analysis of the behavior of multilayered structure considering body forces and initial stresses are agreed with those of finite element analysis. And results of back analysis of elastic moduli and initial stresses in each layers are agreed with exact values with a little difference. Therefore, it is known that this study can be efficiently applied for analyzing the behavior of underground structures including back analysis problems.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1991.10a
• /
• pp.216-232
• /
• 1991

This Finite Element Program(TAFEM) has been developed to be able to carry out the structural analsis of tunnel section and simulate the surrounding ground behaviour due to New Austrian Tunnelling Method, of which main support is the surrounding ground, itself. The Elasto-plastic theory has been applied. The used finite elements are 8-noded isoparametric element(rock & shotcrete), 2 or 3-noded rod element(rock bolt) and infinite boundary element. The load incremental method and tangential stiffness method has been used. Associated flow rule was applied to plastic flow and yield criteria inclued not only Mohr-Coulomb but also Drucker-Prager. In this paper, Drucker-Prager yield criterion has been used. The relationship between plastic strain and stress is based on the incremental strain concept and stress-strain equation on the basis of the stress path of each gauss point has been adopted. It may be rational that rock is considered to be no-tension material, so that no-tension analysis has been adopted in accordance with the brittle fracture constitutive equation.

• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.7
• /
• pp.158-167
• /
• 1999

Heat and moisture transfer associated with porous materials are investigated. The heat and moisture transfer in porous materials caused by the interaction of moisture gradient, temperature gradient, conduction, and evaporation are considered. The variations of temperature and moisture not only change the volume but also induce the hygro-thermal stress. The finite element formulation for solving the temperature and moisture transfer as well as the associated hygro-thermal stresses is developed. In order to verify the finite element formulation, the heat and moisture moving boundary problem in a half space and the hygro-thermo-mechanical problem in an infinite plate with a circular hole are analyzed. Temperature profile, moisture profile, and hygro-thermal stresses are compared with those of analytic solution and other investigator. Good agreements are examined

• Journal of KSNVE
• /
• v.6 no.5
• /
• pp.617-624
• /
• 1996

Finite Element Method(FEM) is one of the most popularly used method in analyzing the dynamic behaviors of structures. But unless the number of finite elements is large enough, the results from FEM are somewhat different form exact analytical solutions, especially at high frequency range. On the other hand, as the Spectral Element Method(SEM) deals directly with the governing equations of structures, the results from this method cannot but be exact regardless of any frequency range. However, despite two dimensional structures are more general, the SEM has been applied only to the analysis of one dimensional structures so far. In this paper, therefore, new methodologies are introduced to analyze the two dimensional plate structure using SEM. The results from this new method are compared with the exact analytical solutions by letting the two dimensional plate structure be one dimensional and showed the dynamic responses of two dimensional plate by including various waves propagated into x-direction.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.234-239
• /
• 2001

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isometric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• Journal of KSNVE
• /
• v.9 no.1
• /
• pp.103-112
• /
• 1999

A new modeling and analysis technique for one-dimensional distributed parameter systems is presented. First. discretized equations of motion in Laplace domain are derived by applying discretization methods for partial differential equations of a one-dimensional structure with respect to spatial coordinate. Secondly. the z and inverse z transformations are applied to the discretized equations of motion for obtaining a dynamic matrix for a uniform element. Four different discretization methods are tested with an example. Finally, taking infinite on the number of step for a uniform element leads to an exact dynamic matrix for the uniform element. A generalized modal analysis procedure for eigenvalue analysis and modal expansion is also presented. The resulting element dynamic matrix is tested with a numerical example. Another application example is provided to demonstrate the applicability of the proposed method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.1
• /
• pp.45-53
• /
• 1991

The underground structure, which has infinite or semi-infinite boundary conditions, is subjected by body forces and in-situ stresses. It also has stress concentration, which causes material nonlinear behavior, in the vicinity of the excavated surface. In this paper, some methods which can be used to transform domain integrals into boundary integrals are reviewed in order to analyze the effect of the body forces and the in-situ stresses. First, the domain integral of the body force is transformed into boundary integral by using the Galerkin tensor and divergence theorem. Second, it is transformed by writing the domain integral in cylindrical coordinates and using direct integration. The domain integral of the in-situ stress is transformed into boundary integral applying the direct integral method in cylindrical coordinates. The methodology is verified by comparing the results from the boundary element analysis with those of the finite element analysis. Coupling the above boundary elements with finite elements, the nonlinear behavior that occurs locally in the vicinity of the excavation is analyzed and the results are verified. Thus, it is concluded that the domain integrals of body forces and in-situ stresses could be performed effectively by transforming them into the boundary integrals, and the nonlinear behavior can be reasonably analyzed by coupled nonlinear finite element and boundary element method. The result of this research is expected to he used for the analysis of the underground structures in the effective manner.

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