بررسی عددی اثر نسبت ناهمسانگردی، زوایه‌ی لایه‌بندی و طول ترک بر ضریب شدت تنش در سنگ فیلیت

نوع مقاله : مقاله پژوهشی

نویسندگان

1 گروه مکانیک سنگ، دانشکده فنی مهندسی، دانشگاه تربیت مدرس، تهران، ایران

2 گروه ژئومکانیک، دانشکده مهندسی علوم زمین، دانشگاه صنعتی اراک، اراک، ایران

چکیده

بررسی رفتار سنگ تحت تنش یکی از موارد مهم در طراحی سازه‌های سنگی است. ناهمسانگردی یکی از شاخصه‌های سنگ‌های رسوبی و دگرگونی است که در پایداری سازه‌های روباز و زیرزمینی تأثیرگذار است. سنگ‌های ناهمسانگرد خواص مکانیکی متفاوتی در جهت‌های مختلف دارند که در طراحی سازه‌های سنگی بایستی مدنظر قرار گیرد. به جهت اهمیت این موضوع، بررسی‌های گسترده‌ای برای شناسایی اثر ناهمسانگردی در رفتار مکانیکی سنگ شده است. روش‌های عددی که قادر به شبیه‌سازی فرایند ترک‌خوردگی در سنگ هستند، می‌توانند شکست، تمرکز و گسترش ریزترک‌ها را در اطراف فضای زیرزمینی پیش‌بینی کنند. بررسی تغییرات ضریب شدت تنش در سنگ‌های ناهمسانگرد و زاویه ترک با محور بارگذاری که در آن کشش و برش خالص و حالت ترکیبی رخ می‌دهد، سؤال اصلی این تحقیق است. در این تحقیق از روش المان محدود و نرم‌افزار آباکوس به دلیل توانایی آن در شبیه‌سازی سنگ ناهمسانگرد برای بررسی تغییرات ضریب شدت تنش استفاده شد، بدین منظور اثر نسبت ناهمسانگردی بر ضریب شدت تنش نمونه‌های ایزوتروپ عرضی در زوایای مختلف لایه‌بندی موردتحقیق قرار گرفت. با در دست داشتن مشخصات مدل‌های دیسکی CCNBD متعلق به سنگ فیلیت مطابق با استاندارد‌های ISRM نمونه‌هایی با نسبت ناهمسانگردی E/E'=1.8) و(E/E'=4  در زوایای مختلف لایه‌بندی (ψ=0-30°-45°-60°-75°-90°) و در زوایای مختلف ترک با امتداد بارگذاری (β) و در طول ترک‌های مختلف (12-24-36 میلی‌متر) با مدل‌سازی عددی موردبررسی قرار گرفت. در همه حالت‌های بررسی‌شده در طول ترک‌های مختلف و در نسبت ناهمسانگردی‌های مختلف برای یک زاویه لایه‌بندی یکسان ضریب شدت تنش مود کشش در مقدار β=0 حداکثر بوده و در β=90° حداقل است و ضریب شدت تنش مود برشی در β=0 و β=90° حداقل و β بین 30 تا 40 درجه حداکثر است. تغییر زاویه ترک نسبت به محور بارگذاری (β) در همه حالت‌های طول ترک و نسبت ناهمسانگردی تأثیر محسوسی در ضریب شدت تنش مود کشش و برش دارد. درواقع حالت کشش خالص در زوایای مختلف لایه‌بندی لزوماً در زاویه صفر درجه با محور بارگذاری (β=0) رخ نمی‌دهد و این زاویه کاملاً متأثر از زوایای لایه‌بندی (Ѱ شیستوزیته) است.

کلیدواژه‌ها

موضوعات


عنوان مقاله [English]

Numerical investigating on the effect of anisotropy ratio, layering angle, and crack length in stress intensity factor in phyllite rock

نویسندگان [English]

  • Vahid Abbasi 1
  • Morteza Ahmadi 1
  • Ehsan Mohtarami 2
1 Dept. of Rock Mechanics , Faculty of Engineering, Tarbiat Modares University, Tehran, Iran
2 Dept of Geomechanics, Faculty of Earth Sciences Engineering, Arak University of Technology, Arak, Iran
چکیده [English]

Summary
Geomechanical properties of some sedimentary and metamorphic rocks are important because they are anisotropic.  These properties affect on stability of rock construction. In this research, the changes in the stress intensity factor in anisotropic rock have been investigated based on the changes in the crack angle with the loading axis and anisotropy ratio. In this research, the finite element method and ABAQUS software are used because of their ability to simulate anisotropic rock. The CCNBD disc sample was used according to the standards of ISRM. It concluded that the stress intensity factor is dependent on the anisotropy ratio, layering angle, and crack length and pure tension (mode one) does not necessarily occur at a zero-crack angle with the loading axis (β=0), and this angle is completely affected by the layering angles (Ѱ).
Introduction
Investigating the behavior of rock under load and determining its mechanical properties are important to designing rock structures. The importance increases when the rock is anisotropic. In the design two mechanical properties; the first rock strength and the second deformation should be considered. The anisotropy rocks require the simultaneous investigation of the resistance behavior and the failure mechanism [1]. A rock is called anisotropic when the value of mechanical properties is different in two different directions. Fracture toughness is the critical stress intensity factor of a sharp crack where propagation of the crack suddenly becomes rapid and unlimited [2].
For materials that have linear elastic behavior, this property is usually expressed in terms of a quantity called the critical value of the stress intensity factor. This parameter is indicated by the symbol K [3] In this research, the main question is to investigate the changes in the stress intensity factor in anisotropic rocks by changing the angle of the crack with the loading axis, anisotropy ratio, crack length, and layering angle.
Methodology and Approaches
For investigation of research question numerical method was considered. The software used Abaqus finite element(FEM) [4]. Abaqus is a set of powerful engineering simulation programs based on the FEM and can solve a wide range of problems from a relatively simple linear analysis to very complex nonlinear analysis. The rock sample was an anisotropic Phyllite belonging to the abutment of the Azad Dam in Kurdistan and the mechanical properties of the sample were taken from Abadat (2013) [5]. CCNBD disk numerical models (D=54 mm and t=30 mm) of Phyllite according to the standard of the International Society of Rock Mechanics (ISRM) with anisotropy ratio (E/E'=1.8 and E/E'=4) in different layering angles (ψ=0-30°-45°-60°-75°-90°) and in different crack angles with loading extension (β) and in different crack lengths (a=12-24-36 mm) were built and executed by Abaqus. After execution, the Abacus determined the stress intensity factor of each model, and data were collected and analyzed.
Results and Conclusions
Sensitivity analysis by variation of four parameters which were anisotropy ratio, different layering angles, different crack angles with loading extension, and different crack lengths was carried out. 180 models were built and executed by Abacus software. Three parameters were kept constant the last was changed and the stress intensity factor was evaluated. Figures 2, 3, and 4 show the effect of the anisotropy ratio variation on the angle of the crack with the loading axis in the pure tensile (mode one), pure shear (mode two), and mixed-mode in the crack length of 12 mm respectively.
1- Variations of layer angle (ψ), crack angle with loading axis (β), crack length(a), and change of anisotropy ratio have a noticeable effect on the values of the stress intensity factor.
2- In different crack lengths and in different anisotropy ratios, for a constant layering angle, the stress intensity factor of tensile mode is maximum at β=0 and minimum at β=90° and the stress intensity factor of shear mode is minimum at β=0 and β=90̊ and maximum at β between 30̊-40̊.
3- Mode one(pure tensile) in different angles of layering does not necessarily occur at a zero-degree angle with the loading axis (β=0) and this angle is completely affected by the angles of layering (Ѱ).
4- According to the results in different layering angles and in different anisotropy ratios, the β angle is different, so that the pure tensile in the layering angle of 30° and the anisotropy ratio equal to 4 occur in β =6°. However, this value is equal to 2° for the anisotropy ratio equal to 1.8

کلیدواژه‌ها [English]

  • Anisotropic rock
  • stress intensity factor
  • anisotropy ratio
  • layering angle
  • crack length and numerical modeling
[1]                 Xing, Yan, P. H. S. W. Kulatilake, and L. A. Sandbak. "Effect of rock mass and discontinuity mechanical properties and delayed rock supporting on tunnel stability in an underground mine." Engineering Geology 238 (2018): 62-75.‏.
[2]                 Song, Honghua, et al. "Scale effects and strength anisotropy in coal." International Journal of Coal Geology 195 (2018): 37-46.‏
[3]                 Jalalifar H., Jafari mohamadabadi B., Shahriar K. and Ahangari K., "Numerical simulation to predict the anisotropy of Faryab chromite rock using PFC software", Journal of Analytical and Numerical Methods in Mining Engineering, Volume 11, Issue 28, October 2021, Pages 23-35(In Persian).
[4]                 Behnia M., Azizian F.and Ghazvinian A.,"Numerical Modeling of Failure Mechanism of Transversely Isotopic Under Triaxial Loading By Discrete Element Method", Journal of Analytical and Numerical Methods in Mining Engineering, Volume 6, Issue 12, Mar. 2014, Pages 67-76. (In Persian).
[5]                 Javid Rad F., "Fracture mechanics and its applications in engineering", Publisher Aerospace industries, 1384(In Persian).
[6]                 Hadi H., Shahriar K., Fatehi Marji M., 2015, "Simulating the effect of disc erosion in TBM disc cutters by a semi-infinite DDM", Arabian Journal of Geosciences volume 8, pages3915–3927.
[7]                 Abdollahipour A., Fatehi Marji M., Yarahmadi Bafghi A., Gholamnejad J., 2016, "Time-dependent crack propagation in a poroelastic medium using a fully coupled hydromechanical displacement discontinuity method" International Journal of Fracture No. 199, Pages 71–87.
[8]                 Haeri H., Sarfarazi V., Ebneabbasi P., Shahbazian A., Nazari Maram A., Shahbazian A. Fatehi Marji M. and Mohamadi A. R., 2020, "XFEM and experimental simulation of failure mechanism of non-persistent joints in mortar under compression", Construction and Building Materials, Volume 236, 10 March 2020, 117500.
[9]                 Y. Tien and M. Kuo, "A failure criterion for transversely isotropic rocks," Int. J. Rock Mech. Min., vol. 38, no. January, pp. 399-412, 2001.
[10]              Sabatakakis, N., and G. Tsiambaos. "Anisotropy of Central Macedonia Phyllite and its effect on the uniaxial compressive strength." Bull PWRC 1 (1983): 2-6.
[11]              Tsidzi, Kwami Edwin Noble. "The influence of foliation on point load strength anisotropy of foliated rocks." Engineering Geology 29.1 (1990): 49-58.
[12]              Alkılıçgil, Çiğdem. "Development of a new method for mode I fracture toughness test on disc type rock specimens", MSc. thesis. 2006.‏.
[13]              Clough, R. W. (1960). "The finite element method in plane stress analysis", In Proceedings of 2nd ASCE Conference on Electronic Computation, Pittsburgh Pa., Sept. 8 and 9, 1960.
[14]              Jing, L. (2003). "A review of techniques, advances and outstanding issues in numerical modelling for rock mechanics and rock engineering", International Journal of Rock Mechanics and Mining Sciences, 40(3), 283-353.
[15]              Wang, Q., et al., "The flattened Brazilian disc specimen used for testing elastic modulus, tensile strength and fracture toughness of brittle rocks: analytical and numerical results", International Journal of Rock Mechanics and Mining Sciences, 2004. 41(2): p. 245-253.
[16]              Huang, Y., "A user-material subroutine in cooperating single crystal plasticity in the ABAQUS finite element program ", 1991: Harvard Uni.
[17]              Abadat A., "Investigating on fracture toughness changes of anisotropic rock", MSc thesis Tarbiat Modares Uni., 1393(In Persian).
[18]              Hibbitt, K.S., ABAQUS/CAE USUr's Manual. 2000