• Advances in Computational Design
• /
• v.5 no.3
• /
• pp.305-321
• /
• 2020

In this paper is presented the solution method for three-dimensional problem of transversely isotropic body's elastoplastic deformation by the finite element method (FEM). The process of problem solution consists of: determining the effective parameters of a transversely isotropic medium; construction of the finite element mesh of the body configuration, including the determination of the local minimum value of the tape width of non-zero coefficients of equation systems by using of front method; constructing of the stiffness matrix coefficients and load vector node components of the equation for an individual finite element's state according to the theory of small elastoplastic deformations for a transversely isotropic medium; the formation of a resolving symmetric-tape system of equations by summing of all state equations coefficients summing of all finite elements; solution of the system of symmetric-tape equations systems by means of the square root method; calculation of the body's elastoplastic stress-strain state by performing the iterative process of the initial stress method. For each problem solution stage, effective computational algorithms have been developed that reduce computational operations number by modifying existing solution methods and taking into account the matrix coefficients structure. As an example it is given, the problem solution of fibrous composite straining in the form of a rectangle with a system of circular holes.

• Tunnel and Underground Space
• /
• v.27 no.4
• /
• pp.243-252
• /
• 2017

Since the failure behavior of transversely isotropic rocks is significantly different from that of isotropic rocks, it is necessary to develop a transversely isotropic rock failure function in order to evaluate the stability of rock structures constructed in transversely isotropic rock masses. In this study, a spatial distribution function for strength parameters of transversely isotropic rocks is proposed, which is based on the Cassini oval curve proposed by 17th century astronomer Giovanni Domenico Cassini to model the orbit of the Sun around the Earth. The proposed distribution function consists of two model parameters which could be identified through triaxial compression tests on transversely isotropic rock samples. The original Mohr-Coulomb (M-C) failure function is extended to a three-dimensional transversely isotropic M-C failure function by employing the proposed strength parameter distribution function for the spatial distributions of the friction angle and cohesion. In order to verify the suitability of the transversely isotropic M-C failure function, both the conventional triaxial compression and true triaxial compression tests of transversely isotropic rock samples are simulated. The predicted results from the numerical experiments are consistent with the failure behavior of transversely isotropic rocks observed in the actual laboratory tests. In addition, the simulated result of true triaxial compression tests hints that the dependence of rock strength on intermediate principal stress may be closely related to the distribution of the microstructures included in the rock samples.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.13 no.6
• /
• pp.67-86
• /
• 2009

In this study, a 3D transmitting boundary in water-saturated transversely isotropic soil strata has been developed based on u-w formulation for application to general 3D analysis. Behavior in the far field region is expanded in the Fourier series, and dynamic stiffness for each term is obtained based on the u-w formulation. Transformation of the dynamic stiffness is presented to combine the transmitting boundary with the 3D finite elements for the near field region formulated in a 3D Cartesian coordinate system. The developed transmitting boundary is verified through a comparison of the dynamic behavior of a rigid circular foundation with the results from the existing numerical method. In addition, the developed transmitting boundary is applied to the analysis of the dynamic behavior of rigid foundations of diverse shapes, and the effects of the level of the groundwater table on the dynamic stiffness of a rigid rectangular foundation in the water-saturated transversely isotropic layered stratum are studied.

• Tunnel and Underground Space
• /
• v.26 no.2
• /
• pp.100-109
• /
• 2016

This study suggests three spatial distribution functions of strength parameters, which can be adopted in the derivation of failure conditions for transversely isotropic rocks. All three proposed functions, which are the oblate spheroidal function, the exponential function, and the function based on the directional projection of the strength parameter tensor, consist of two model parameters. With assumption that the cohesion and friction angle can be described by the proposed distribution functions, the transversely isotropic Mohr-Coulomb criterion is formulated and used as a failure condition in the simulation of the conventional triaxial tests. The simulation results confirm that the failure criteria incorporating the proposed distribution functions could reproduce the general trend in the variations of the axial stress at failure and the directions of failure planes with varying inclination of the weankness planes and confining pressure. Among three distribution functions, the function based on the directional projection of the strength parameter tensor yields the highest axial strength, while the axial strength estimated by the oblate spheroidal distribution function is the lowest.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.209-214
• /
• 2007

In Korean geology that crystalline rock is dominant, the properties of subsurface including the anisotropy are distributed complexly and changed abruptly. Because of such geological environments, cross-hole seismic traveltime tomography is widely used to obtain the high resolution image of the subsurface for the engineering purposes in the geotechnical sites. However, because the cross-hole tomography has a wide propagation angle coverage relatively, its data tend to include the seismic velocity anisotropy comparing with the surface seismic methods. It can cause the misinterpretation that the cross-hole seismic data including the anisotropic effects are analyzed and treated with the general processing techniques assuming the isotropy. Therefore, we need to consider the seismic anisotropy in cross-hole seismic traveltime tomography. The seismic anisotropic tomography algorithm, which is developed for evaluation of the velocity anisotropy, includes several inversion schemes in order to make the inversion process stable and robust. First of all, the set of the inversion parameters is limited to one slowness, two ratios of slowness and one direction of the anisotropy symmetric axis. The ranges of the inversion parameters are localized by the pseudo-beta transform to obtain the reasonable inversion results and the inversion constraints are controlled efficiently by ACB(Active Constraint Balancing) method. Especially, the inversion using the Fresnel volume is applied to the anisotropic tomography and it can make the anisotropic tomography more stable than ray tomography as it widens the propagation angle coverage.

• Structural Engineering and Mechanics
• /
• v.68 no.4
• /
• pp.443-457
• /
• 2018

In this article, a comparative study has been made to investigate the scattering behaviour of three layered structure model on torsional surface wave. For such model intermediate layer is taken as fiber reinforced composite, resting over a dry sandy Gibson substratum and underlying by different anelastic media. We consider two distinct mediums for topmost layer. In the first case, topmost layer has been taken as fluid saturated homogeneous porous layer, while in the second case the fluid saturated porous layer has been replaced by a transversely isotropic layer. Simple form expression for the secular equation of torsional surface wave has been worked out in both the cases by executing specific boundary conditions, which comprises Whittaker's function and its derivative, for imminent result that have been elaborated asymptotically. Some special cases have been constituted which are in excellent compliance with recorded literatures. For the sake of comparative study, numerical estimation and graphical illustration have been accomplished to identify the effects of the width ratio of the layers, Biot's gravity parameter, sandy parameter, porosity parameter and other heterogeneity parameters corresponding to the layers and half spaces, horizontal compressive and tensile initial stress on the phase velocity of torsional surface wave.

• Geophysics and Geophysical Exploration
• /
• v.16 no.2
• /
• pp.71-78
• /
• 2013

The one-way wave equation migration is much more computationally efficient comparing with reverse time migration and it can provide better image than the migration algorithm based on the ray theory. We have developed the prestack depth migration module adopting (GS) propagator designed for vertical transverse isotropic media. Since GS propagator considers the higher-order term by expanding the Taylor series of the vertical slowness in the thin slab of the phase-screen propagator, the GS migration can offer more correct image for the complex subsurface with large lateral velocity variation or steep dip. To verify the validity of the developed GS migration module, we analyzed the accuracy with the order of the GS propagator for VTI media (GSVTI propagator) and confirmed that the accuracy of the wavefield propagation with the wide angles increases as the order of the GS propagator increases. Using the synthetic seismic data, we compared the migration results obtained from the isotropic GS migration module with the anisotropic GS migration module. The results show that the anisotropic GS migration provides better images and the improvement is more evident on steeply dipping structures and in a strongly anisotropic medium.

• Tunnel and Underground Space
• /
• v.33 no.3
• /
• pp.150-168
• /
• 2023

The anisotropic behavior of rocks is primarily attributed to the directional arrangement of rock-forming minerals and the distribution characteristics of microcracks. Notably, sedimentary and metamorphic rocks often exhibit distinct transverse isotropy in terms of their strength and deformation characteristics. Consequently, it is crucial to gain accurate insights into the deformation and failure characteristics of transversely isotropic rocks during rock mechanics design processes. The deformation of such rocks is described by five independent elastic constants, which are determined through laboratory testing. In this study, the characteristics of the directional variation of apparent elastic constants in transversely isotropic rocks were investigated using experimental data reported in the literature. To achieve this, the constitutive equation proposed by Mehrabadi & Cowin was introduced to calculate the apparent elastic constants more efficiently and systematically in a rotated Cartesian coordinate system. Four transversely isotropic rock types from the literature were selected, and the influence of changes in the orientation of the weak plane on the variations of the apparent elastic modulus, apparent shear modulus, and apparent Poisson's ratio was analyzed. Based on the investigation, a new constraint on the elastic constants has been proposed. If the proposed constraint is satisfied, the directional variation of the apparent elastic constants in transversely isotropic rocks aligns with intuitive predictions of their tendencies.

• Tunnel and Underground Space
• /
• v.18 no.4
• /
• pp.252-262
• /
• 2008

In the uniaxial compressive test of a single specimen of transversely isotropic rock, its 5 independent elastic constants can not be defined since maximum 4 independent strain measurements are available theoretically. In order to solve this problem, one equation proposed by Saint Venant in 19C and confirmed by Lekhnitskii through the test experiences has been used for long time. Accordign to authors' experiences, however, this equation turned out to give erroneous elastic constants in some cases. Three new equations are suggested and their compatibilities are discussed in this paper. As the results of the analyses of the models, Lekhnitskii's suggested equation is effective for the specimen with the high dip angle whereas it results in the large erred output for that with dip angle less than $25{\sim}30$. It was found that the effectivenesses of three suggested equations and their compatibilities are subject to the dip angle and not to the amounts of elastic constants. Guide map to the selection of the compatible one of those suggested equations is presented as a result of the study.

• Journal of the Korean Geotechnical Society
• /
• v.17 no.1
• /
• pp.5-13
• /
• 2001

강도 또는 변형 이방성은 층리진 최적암, 박층모양의 변성암, 균일하게 절 리가 나있는 암석에서 흔히 볼 수 있다(Amadei, 1982, 1996). 특히, 대구지역은 퇴적암의 일종인 셰일로 구성되어 층리면의 각도에 따라 역학적인 특성이 달리 산정되어져야 한다. 이에 본 연구에서는 층리면의 각도를 달리하여 시료를 성형한 후 각종 시험을 수행하였다. 그 결과 일축 압축강도는 수평면과 층리면이 이루는 각이 0$^{\circ}$, 90$^{\circ}$에서 최고 강도를, 60$^{\circ}$에서 최저 강도를 나타내었다. 또한 간접 인장시험과 점 하중시험, 탄성파시험에서는 90$^{\circ}$에서 최고값을, 0$^{\circ}$에서 최저값을 보였다. 그리고, 층리에 따른 퇴적압의 일축 압축강도를 간접 인장강도, 점하중강도와의 상관성을 분석하여 각각의 관계를 회귀분석을 통한 상관식으로 나타내었다. 추가적으로 층리에 따른 공학적특성을 정량적으로 나타내었다.