• Geomechanics and Engineering
• /
• v.22 no.1
• /
• pp.11-23
• /
• 2020

In this study, a new analytical-numerical procedure is developed to give the stresses and strains around a circular tunnel in rock masses exhibiting different stress-strain behavior. The calculation starts from the tunnel wall and continues toward the unknown elastic-plastic boundary by a finite difference method in the annular discretized plastic zone. From the known stresses in the tunnel boundary, the strains are calculated using the elastic-plastic stiffness matrix in which three dimensional Hoek-Brown failure criterion (Jiang and Zhao 2015) and Mohr-Coulomb potential function with proper dilation angle (i.e., non-associated flow rule) are employed in terms of stress invariants. The illustrative examples give ground response curve and show correctness of the proposed approach. Finally, from the results of a great number of analyses, a simple relationship is presented to find out the closure of circular tunnel in terms of rock mass strength and tunnel depth. It can be valuable for the preliminary decision of tunnel support and for prediction of tunnel problems.

• Tunnel and Underground Space
• /
• v.27 no.1
• /
• pp.39-49
• /
• 2017

The factor of safety (FOS) is commonly used as an index to quantitatively state the degree of safety of various rock structures. Therefore it is important to understand the definition and characteristics of the adopted FOS because the calculated FOS may be different according to the definition of FOS even if it is estimated under the same stress condition. In this study, four local factors of safety based on maximum shear stress, maximum shear strength, stress invariants, and maximum principal stress were defined using the Mohr-Coulomb and Hoek-Brown failure criteria. Then, the variation characteristics of each FOS along five stress paths were investigated. It is shown that the local FOS based on the shear strength, which is widely used in the stability analysis of rock structures, results in a higher FOS value than those based on the maximum principal stress and the stress invariants. This result implies that the local FOS based on the maximum shear stress or the stress invariants is more necessary than the local FOS based on the shear strength when the conservative rock mechanics design is required. In addition, it is shown that the maximum principal stresses at failure may reveal a large difference depending on the stress path.

• Tunnel and Underground Space
• /
• v.22 no.6
• /
• pp.462-470
• /
• 2012

Implementing the generalized Hoek-Brown failure criterion in the framework of the Mohr-Coulomb criterion requires the calculation of the equivalent friction angle and cohesion. In the conventional method based on the Balmer (1952)'s theory, the tangential instantaneous friction angle and cohesion are expressed in terms of the minimum principal stress ${\sigma}_3$, which does not provide the information about the dependency of the equivalent parameters on the hydrostatic pressure and the stress path. In this study, this defect of the conventional method has been overcome by representing the equivalent parameters in terms of stress invariants. Through the example implementation of the new method, the influence of the magnitude of the hydrostatic pressure and the Lode angle on the tangential instantaneous friction angle and cohesion is investigated. It turns out that the tangential instantaneous friction angle is maximum when the stress condition is triaxial extension, while the tangential cohesion is maximum when the stress condition is triaxial compression. The dependency of the equivalent Mohr-Coulomb strength parameters on the hydrostatic pressure and the Lode angle tends to be more substantial for the favorable rockmass of larger GSI value.

• Transactions of Materials Processing
• /
• v.31 no.6
• /
• pp.351-364
• /
• 2022

The yield criterion, or called yield function, plays an important role in the study of plastic working of a sheet because it governs the plastic deformation properties of the sheet during plastic forming process. In this paper, we propose a modified version of previous anisotropic yield function (Trans. Mater. Process., 31(4) 2022, pp. 214-228) based on J2 and J3 stress invariants. The proposed anisotropic yield model has the 6th-order of stress components. The modified version of the anisotropic yield function in this study is as follows. f(J₂⁰,J₃⁰) ≡ (J₂⁰)³ + α(J₃⁰)² + β(J₂⁰)^3/2 × (J₃⁰) = k⁶ The proposed anisotropic yield function well explains the anisotropic plastic behavior of various sheets such as aluminum, high strength steel, magnesium alloy sheets etc. by introducing the parameters α and β, and also exhibits both symmetrical and asymmetrical yield surfaces. The parameters included in the proposed model are determined through an optimization algorithm from uniaxial and biaxial experimental data under proportional loading path. In this study, the validity of the proposed anisotropic yield function was verified by comparing the yield surface shape, normalized uniaxial yield stress value, and Lankford's anisotropic coefficient R-value derived with the experimental results. Application for the proposed anisotropic yield function to AA6016-T4 aluminum and DP980 sheets shows symmetrical yielding behavior and to AZ31B magnesium shows asymmetric yielding behavior, it was shown that the yield locus and yielding behavior of various types of sheet materials can be predicted reasonably by using the proposed anisotropic yield function.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.22 no.1
• /
• pp.50-55
• /
• 2013

In order to understand the biomaterial like the blood vessel of artery, there is a need to quantify the biomechanical behavior of the vessel. However, theoretical framework to describe and quantify the behaviour of blood vessel was not well established so far. For studying the biomechanical behavior of artery, Rubber-liked material which is similar to passive artery is selected since conventional theoretical interpretation is very limited to understand and predict the behavior of biomaterial. Rubber-like material is assumed to be very similar to artery and has properties of isotropy, homogeneity and is undergoing large deformation. Based on this assumption, stress developed on Rubber-like material is described by strain energy function and strain invariants which are required to understand the nonlinear elastic behavior of biomaterial. The descriptor which would be used for understanding the biomechanical behavior of artery is studied in this work.

• Elastomers and Composites
• /
• v.56 no.2
• /
• pp.79-84
• /
• 2021

A strain energy density function is used to characterize the hyperelasticity of rubber-like materials. Conventional models, such as the Neo-Hookean, Mooney-Rivlin, and Ogden models, are widely used in automotive industries, in which the strain potential is derived from strain invariants or principal stretch ratios. A fitting procedure for experimental data is required to determine material constants for each model. However, due to the complexities of the mathematical expression, these models can only produce an accurate curve fitting in a specified strain range of the material. In this study, a hyperelastic model for Neodymium Butadiene rubber is developed by using the Artificial Neural Network. Comparing the analytical results to those obtained by conventional models revealed that the proposed model shows better agreement for both uniaxial and equibiaxial test data of the rubber.

• Transactions of Materials Processing
• /
• v.31 no.4
• /
• pp.214-228
• /
• 2022

The yield criterion, or called yield function, plays an important role in the study of plastic working of a sheet because it governs the plastic deformation properties of the sheet during plastic forming process. In this paper, we propose a novel anisotropic yield function useful for describing the plastic behavior of various anisotropic sheets. The proposed yield function includes the anisotropic version of the second stress invariant J₂ and the third stress invariant J₃. The anisotropic yield function newly proposed in this study is as follows. F(J₂)+ αG(J₃)+ βH (J₂ × J₃) = k^m The proposed yield function well explains the anisotropic plastic behavior of various sheets by introducing the parameters α and β, and also exhibits both symmetrical and asymmetrical yield surfaces. The parameters included in the proposed model are determined through an optimization algorithm from uniaxial and biaxial experimental data under proportional loading path. In this study, the validity of the proposed anisotropic yield function was verified by comparing the yield surface shape, normalized uniaxial yield stress value, and Lankford's anisotropic coefficient R-value derived with the experimental results. Application for the proposed anisotropic yield function to aluminum sheet shows symmetrical yielding behavior and to pure titanium sheet shows asymmetric yielding behavior, it was shown that the yield curve and yield behavior of various types of sheet materials can be predicted reasonably by using the proposed new yield anisotropic function.

• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.11
• /
• pp.150-157
• /
• 2003

In this paper, the large deformation of hollow silicon rubber gasket is treated. The frictional contact occurs between groove and the outer part of hollow gasket, and the frictional self-contact exists in the inner parts of hollow gasket. The silicon rubber has the nonlinear elastic behavior and its material property is approximately incompressible. Hence, the stress analysis requires an existence of a strain energy function, which is usually defined in terms of invariants or stretch ratio such as generalized Mooney-Rivlin and Ogden model. Considering large compressive deformation and friction, Mooney-Rivlin 3rd model and Coulomb's friction model are assumed. The numerical analysis is obtained by the commercial finite element program MARC. But, due to large deformation, the elements degenerate in the inner parts of hollow gasket. This means that the analysis of subsequent increments is carried out with a very poor mesh. In order to continue the analysis with a sufficient accuracy, it is necessary to use new finite element modeling by remesh. Experiments are also performed to show the validity of present method. As a conclusion, numerical results by this research have good agreements with experiments.

• Geomechanics and Engineering
• /
• v.9 no.6
• /
• pp.793-813
• /
• 2015

Characterization of discontinuous media is an endeavor that poses great challenge to engineers in practice. Since the inherent defects in cracked domains can substantially influence material resistance and govern its behavior, a lot of work is dedicated to efficiently model such effects. In order to overcome difficulties of material instability problems, one needs to comprehensively represent the geometry of cracks along with their impact on the mechanical properties of the intact material. In the present study, stress-strain results from laboratory experiments on Inada granite was used to derive crack tensor as a tool for the evaluation of fractured domain stability. It was found that the formulations proposed earlier could satisfactorily be employed to attain crack tensor via the invariants of which judgment on cracks population and induced anisotropy is possible. The earlier criteria based on crack tensor analyses were reviewed and compared to the results of the current study. It is concluded that the geometrical parameters calculated using mechanical properties could confidently be used to judge the anisotropy as well as strength of the cracked domain.

• Structural Engineering and Mechanics
• /
• v.82 no.1
• /
• pp.121-131
• /
• 2022

The main goal of this study is to prepare a program for analyzing High Strength Steel Fibrous Reinforced Concrete (HSSFRC) slabs and predict the response and strength of the slab instead of preparing a prototype and testing it in the laboratory. For this purpose, new equations are proposed to represent the material properties of High Strength Steel Fibrous Reinforced Concrete. The proposed equations obtained from performing regression analysis on many experimental results using statistical programs. The finite element method is adopted for non-linear analysis of the slabs. The eight-node "Serendipity element" (3 DoF) is chosen to represent the concrete. The layered approach is adopted for concrete elements and the steel reinforcement is represented by a smeared layer. The compression properties of the concrete are modeled by a work hardening plasticity approach and the yield condition is determined depending on the first two stress invariants. A tensile strength criterion is adopted in order to estimate the cracks propagation. many experimental results for testing slabs are compared with the numerical results of the present study and a good agreement is achieved regarding load-deflection curves and crack pattern. The response of the load deflection curve is slightly stiff at the beginning because the creep effect is not considered in this study and for assuming perfect bond between the steel reinforcement and the concrete, however, a great agreement is achieved between the ultimate load from the present study and experimental results. For the models of the tension stiffening and cracked shear modulus, the value of Bg and Bt (Where Bg and Bt are the curvature factor for the cracked shear modulus and tension stiffening models respectively) equal to 0.005 give good results compared with experimental result.