Stress Direction Research Articles

A DFT study of the CaMg2As2 material for photovoltaic applications

The research focused on investigating the structural, elastic, electronic, optical, and thermoelectric characteristics of the CaMg2As2 compound. To inspect these properties, the density functional theory (DFT) method has been applied via the Wien2K code. The electronic analysis revealed that the CaMg2As2 compound demonstrates semiconductor behavior with an indirect band gap of 1.781 eV. Additionally, its elastic properties were examined, revealing a 1.980 % AVR, indicating a variation in elasticity depending on the direction of load or stress applied. We have also studied various optical properties of CaMg2As2, including the refractive index, extinction coefficient, electron energy loss, dielectric tensor, and optical conductivity. the absorption coefficient value is zero for energy levels below 2.60 eV, meaning there is no energy absorption by the material at lower energy levels. However, as the energy surpasses 2.60 eV, the absorption coefficient increases and shows multiple peaks. This signifies that as the incident radiation's energy increases, the material starts absorbing more energy at specific energy levels. In addition, one significant finding is that the thermal conductivity of the lattice (κL) in CaMg2As2 decreases exponentially as the temperature rises. This means that as the material gets hotter, its ability to conduct heat becomes less effective. In simple terms, the material becomes less efficient at transferring heat at higher temperatures. These findings suggest that CaMg2As2 has potential as a thermoelectric material with interesting properties that deserve further investigation.

Fatigue Behaviour and Life Prediction of YSZ Thermal Barrier Coatings at Elevated Temperature under Cyclic Loads

The concentration of interfacial normal stress at the free edges of thermal barrier coatings (TBCs) can result in coating spallation. Fatigue cracking is one of the main reasons for creating free edges under complex loads. It is crucial to investigate the fatigue cracking of coatings under cyclic loads to assess potential coating failure. To address this issue, a novel model was proposed to predict the fatigue life of the YSZ topcoat under stress parallel to the interface. Firstly, this study conducted uniaxial and tensile-torsional fatigue tests at elevated temperatures on specimens with atmospheric plasma-sprayed TBCs. The test results revealed that fatigue cracks appeared in the topcoat under cyclic loads, but these cracks did not propagate into the bondcoat or substrate immediately. The number of cycles before the topcoat cracked was found to be associated with the magnitude of the cyclic load. Secondly, this study analyzed the test conditions using the finite element method. Simulations indicated that the crack direction in the topcoat under complex loading conditions was aligned with the first principal stress direction. Finally, the fatigue life prediction model of the topcoat was established based on experiments and simulations. The predicted results fell within a fourfold scatter band.

Analysis of mechanical properties of rice particle beds in confined compression tests and stress transfer predictive model

High-stress breakage during storage and processing is one of the primary factors in the decline of rice quality. To reveal the mechanical response law under this condition, confined compression and numerical simulations were used to analyze the mechanical evolution laws of rice particle beds under high-pressure in this work. The results show that the compaction process progresses from top to bottom. The total number of contact forces in the particle beds increases, the strong contact forces exhibit a diffusion in direction. The maximum principal stress is consistent with the direction of the strong force chain, particle breakage direction corresponds to the maximum principal stress direction. Finally, based on the analysis of the particle beds, a uniaxial confined compression axial stress transfer model for rice particle beds is proposed. The results are helpful for understanding the stress transfer characteristics of rice particle beds, provide basis for the improvement of related equipment.

Re-orientation of principal stresses during tunnel excavation: Model development and support strategy

The safety and stability during the excavation and production of the tunnel is a key issue that must be addressed in tunneling engineering, which is pronounced under principal stress re-orientation. A mechanical model is therefore developed to determine the boundary curve of the failure zone of a supported circular tunnel under rotating the principal stress direction. Meanwhile, the effects of various environmental factors on the pattern of the failure zone are quantitatively analyzed. Subsequently, two critical morphological indexes are determined to elucidate the uniformity of the failure zone and the functional relationship between these two morphological indices and the environmental factors, respectively. Furthermore, technical approaches and control measures for maintaining the surrounding rock are put forward. In particular, the “long and short” coordinated hierarchical support technology has been recommended and is validated through a comparison of the designed schemes and several engineering cases. The results indicate that the proposed support technology resulted in a gradual reduction of the “butterfly wings” pattern and a more consistent failure zone in the surrounding rock. These findings are useful to ensure the overall stability and anti-interference performance of the surrounding rock of a tunnel under the complicated stress field, especially with regards to the deflection of principal stress.

Effect of the hot forming with synchronous quenching process on high cycle fatigue properties of the 2A97 Al-Li alloy

The 2A97 Al-Li alloys have been used in aircraft to achieve lightweight, its high cycle fatigue performance is one of the important factors that determines aircraft safety. The 2A97-T8 Al-Li alloy is obtained through the 2A97-T3 Al-Li alloy undergoes the hot forming with synchronous quenching process. The 2A97-T3 and 2A97-T8 Al-Li alloys are tested for high cycle fatigue under different stress amplitudes. The fatigue crack propagation path of the 2A97-T3 Al-Li alloy is not only along the direction of maximum shear stress at 45° to the stress direction but also along the direction perpendicular to the stress axis. The fatigue crack propagation path of the 2A97-T8 Al-Li alloy propagates along the direction perpendicular to the stress axis. Under high amplitude stress, the characteristics of quasi-cleavage fracture are found in the fracture morphology of 2A97-T3 Al-Li alloy. However, the characteristics of cleavage fracture are found on the fracture morphology of 2A97-T3 Al-Li alloy under low amplitude stress. Under high amplitude stress, the quasi-cleavage characteristics and the dimples are found on the fracture morphology of 2A97-T8 Al-Li alloy. Under low amplitude stress, the small voids and the dimples are found on the fracture morphology of 2A97-T8 Al-Li alloy. The precipitated phases of the 2A97-T3 Al-Li alloy are multitudinous, fine, and dispersed δ’ phases and a few T1 phases. The precipitates of the 2A97-T8 Al-Li alloy are large, fine, and acicular T1 phases and a few δ’ phases. The hot forming with synchronous quenching process significantly improves the high cycle fatigue properties of the 2A97 Al-Li alloy. In addition, the fatigue cycles, the fatigue life rise and fall plots, the S-N curves, and the fatigue crack growth behaviors of the 2A97-T3 and 2A97-T8 Al-Li alloys are obtained.

Research on crack distribution characteristics and control technology of surrounding rock in soft rock roadway under different lateral pressure coefficients

AbstractTo solve the problem of controlling large deformation of surrounding rock in deep soft rock roadway, the distribution characteristics and deformation mechanism of surrounding rock cracks in soft rock roadway under different lateral pressure coefficients are studied using numerical simulation, theoretical analysis, and field measurement. The results show that under different lateral pressure coefficients, the range of surrounding rock cracks shows three forms: round, oval, and butterfly. No matter what lateral pressure coefficient the roadway is in, the surrounding rock cracks always appear in the plastic zone, and there is a high correlation between the surrounding rock crack range and the plastic zone. The stress characteristics of the surrounding rock in the plastic zone include two main aspects. One is that the direction of the principal stress of the surrounding rock is deflected, which is manifested as an annular distribution of the direction of the maximum principal stress around the roadway. The direction of the minimum principal stress in the upper part of the roadway points to the center of the roadway, and the direction of the minimum principal stress in the lower part of the roadway deviates from the center of the roadway. Second, the ratio of the maximum to minimum principal stress in the surrounding rock is large. Under this stress characteristic, the surrounding rock in the plastic zone has strong shear dilation. The shear dilation makes the crack of the surrounding rock open so that the surrounding rock is squeezed into the roadway space, and then the roadway produces large deformation. Due to the large range of cracks in the butterfly‐shaped plastic zone, the shear dilation deformation produced by the butterfly‐shaped plastic zone is far more than that of the round/oval plastic zone. According to the crack range of roadway surrounding rock under different lateral pressure coefficients, the corresponding support scheme is put forward. Field experiments show that the support scheme can effectively control the deformation of surrounding rock and meet the requirements of roadway use.

Stress–strain Behavior of Silt Sand Soil Under Isotropically Consolidated Drained Conditions by Principal Stress Direction and Intermediate Stress Coefficient Parameters

Stress–strain Behavior of Silt Sand Soil Under Isotropically Consolidated Drained Conditions by Principal Stress Direction and Intermediate Stress Coefficient Parameters

Impact of scan strategy on principal stresses in laser powder bed fusion

Additive manufacturing techniques, such as laser powder bed fusion (PBF-LB), are well known for their exceptional freedom in part design. However, these techniques are also characterized by the development of large thermal gradients during production and thus residual stress (RS) formation in produced parts. In this context, neutron diffraction enables the non-destructive characterization of the bulk RS distribution. By control of the thermal gradients in the powder-bed plane by scan strategy variation we study the impact of in-process scan strategy variations on the microstructure and the three-dimensional distribution of RS. Microstructural analysis by means of electron backscatter diffraction reveals sharp microstructure transitions at the interfaces ranging from 100-200 µm. The components of the RS tensor are determined by means of neutron diffraction and the principal stress directions and magnitudes are determined by eigenvalue decomposition. We find that the distribution of RS in the powder-bed plane corresponds to the underlying scan strategy. When the alternating scan vectors align with the x- and y sample coordinate axes, the principal stress directions co-align. In the present geometry, nearly transverse isotropic stress states develop when the scan vectors are either aligned 45° between x and y or continuously rotated by 67° between each layer.

A modified SWT model for very high cycle fatigue life prediction of L-PBF Ti-6Al-4V alloy based on Single Defect: Effect of building orientation

The ultrasonic fatigue tests on laser-powder bed fused (L-PBF) Ti-6Al-4V specimens manufactured at 0°, 45°, and 90° building orientations were conducted. It is found that the very high cycle fatigue (VHCF) lives decrease with the building orientation from 0°, 45°, to 90°. All the fatigue cracks initiate from internal pores or Lack of Fusion defects. Mode Ι cracks grow perpendicular to the direction of the maximum opening stress, regardless of the building orientation. The maximum prospective defect evaluated by the extreme value statistics (EVS) method is used as the typical defect for VHCF life prediction. A modified Smith-Waston-Topper (SWT) model was proposed and combined with the Critical Distance Method to predict ultra-high cycle fatigue life. The life prediction accuracy was validated by using both experimental data obtained from the experiments conducted in this paper and experimental data taken from the literature. The results show that 84% of the predicted fatigue lives are within a scatter band of 2 standard deviation, and 16% of the predicted lives are outside a scatter band of 2 standard deviation and within a scatter band of 3 standard deviation.

A coupled thermo-mechanical model for investigating cracking and failure of composite interbedded rock

Composite interbedded rock is a special geological phenomenon that occurs when layers of rock (beds) of a given lithology lie between or alternate with other layers of a different lithology. The thermal and mechanical properties of the composite interbedded rocks exhibit strong anisotropy. Studying the thermomechanical coupling properties of natural composite interbedded rocks using existing experimental techniques is challenging. A coupled thermo-mechanical composite interbedded rock model was developed based on Fourier's heat-conduction law and Newton's second law. The thermodynamic responses and biaxial mechanical properties of composite interbedded rocks with different interlayer angles under the coupling of initial stress and thermal loading were studied for the first time. The results indicate that the thermal stresses within the granite interlayers are higher than those within the sandstone interlayers. Throughout the thermal loading process, thermal cracking occurred primarily within the (stiffer) granite interlayers. The number of thermal cracks generated during the cooling process far exceeds that during heating; at 400 °C, the number of thermal cracks during cooling is more than 35 times that during heating. When the interlayer dip angle approaches 90° and is aligned with the major principal stress direction, the phenomenon of the “double peaks” and layered failure becomes more pronounced. Within the range of variables considered, the relative impact of the three variables on peak stress in composite interbedded rock is as follows, from greatest to least: interlayer angle, thermal shock temperature, and model size. The results and developed modeling methods will contribute to the management and optimization of future engineering geology projects.

Numerical simulation of blasting behavior of rock mass with cavity under high in-situ stress

With the shift of coal seam mining to the deep, the in-situ stress of coal and rock mass increases gradually. High ground stress can limit the generation of rock cracks caused by blasting, and blasting usually shows different crushing states than low stress conditions. In order to study the blasting expansion rule of rock mass with cavity under high ground stress and the rock mass fracture state under different side stress coefficients. In this paper, the effective range of blasting and the stress distribution under blasting load are analyzed theoretically. The RHT (Riedel-Hiermaier-Thoma) model is used to numerically simulate the blasting process of rock mass with cavity under different ground stress, and the influence of ground stress and lateral pressure coefficient on the crack growth of rock mass is studied. The results show that when there is no ground stress, the damage cracks in rock mass are more concentrated in the horizontal direction and the fracture development tends to the direction where the holes are located, which confirms the guiding effect and stress concentration effect of the holes in rock mass, which helps to promote the crack penetration between the hole and the hole. The length difference of horizontal and vertical damage cracks in rock mass increases with the increase of horizontal and vertical stress difference. Under the same lateral stress coefficient, the larger the horizontal and vertical stress difference is, the stronger the inhibition effect on crack formation is. For blasting of rock mass with high ground stress, the crack formation length between gun holes decreases with the increase of stress level, and the crack extends preferentially in the direction of higher stress. Therefore, the placement of gun holes along the direction of greater stress and the shortening of hole spacing are conducive to the penetration of cracks between gun holes and empty holes. The research can provide reference for rock breaking behavior of deep rock mass blasting.

A transient radial cortical microtubule array primes cell division in Arabidopsis

Although the formation of new walls during plant cell division tends to follow maximal tensile stress direction, analyses of individual cells over time reveal a much more variable behavior. The origin of such variability as well as the exact role of interphasic microtubule behavior before cell division have remained mysterious so far. To approach this question, we took advantage of the Arabidopsis stem, where the tensile stress pattern is both highly anisotropic and stable. Although cortical microtubules (CMTs) generally align with maximal tensile stress, we detected a specific time window, ca. 3 h before cell division, where cells form a radial pattern of CMTs. This microtubule array organization preceded preprophase band (PPB) formation, a transient CMT array predicting the position of the future division plane. It was observed under different growth conditions and was not related to cell geometry or polar auxin transport. Interestingly, this cortical radial pattern correlated with the well-documented increase of cytoplasmic microtubule accumulation before cell division. This radial organization was prolonged in cells of the trm678 mutant, where CMTs are unable to form a PPB. Whereas division plane orientation in trm678 is noisier, we found that cell division symmetry was in contrast less variable between daughter cells. We propose that this "radial step" reflects a trade-off in robustness for two essential cell division attributes: symmetry and orientation. This involves a "reset" stage in G2, where an increased cytoplasmic microtubule accumulation transiently disrupts CMT alignment with tissue stress.

Influence of tectonic stress field on oil and gas migration in the Tahe Oilfield carbonate reservoir: Identification of areas favorable for reservoir formation

An understanding of the tectonic stress field distribution in the carbonate reservoirs of the Tahe Oilfield is essential for oil and gas migration and reservoir reconstruction. This study adopts a geological genesis perspective and integrates the unique attributes of the carbonate reservoirs in the Tahe Oilfield to simulate three-dimensional stress fields using the finite element method and Petrel software. The simulation, which takes pore pressure and boundary constraints into consideration, aligns with the findings of acoustic emission tests conducted by previous researchers and validates the model by corroborating the maximum horizontal principal stress direction with imaging logging interpretations. With subsequent utilization of the generalized Coulomb-Mohr shear failure criterion and the Griffith generalized maximum tensile stress criterion, the simulated stress enables the prediction of the total fracture intensity. To complement this analysis, we integrate seismic data and fault formation mechanisms to characterize variations in fracture intensity within distinct regions of the study area. Three key findings emerge from this investigation. First, due to the Hercynian thermal event, X-shaped conjugate faults display a disproportionately high degree of development compared to other primary fractures. Analyzing the mechanisms of formation and development of these X-shaped strike-slip faults reveals discernible differences in deformation characteristics between NW-trending and NE-trending faults, with the former exhibiting a greater propensity for fracture rupture and efficient oil-gas migration. Second, fault intersections are critical and efficient conduits for fluid accumulation and thereby contribute to reservoir formation. Lastly, exclusive of X-shaped conjugate faults, NE-trending faults exhibit the greatest propensity for oil and gas accumulation and migration. Collectively, these findings facilitate the identification of areas favorable for oil and gas migration and provide a crucial mechanical analysis that can inform the exploration and development of the Tahe Oilfield.

Fatigue crack propagation life analysis of propeller

Fatigue crack fracture is a common failure mode in structures. For the propeller with initial crack, the hydrodynamic effect cannot be ignored when studying the fatigue crack propagation behavior. Taking the full-size KP505 propeller as the research object, its hydrodynamic performance is simulated and verified. By using the fluid-structure interaction (FSI), the pressure distribution and stress distribution of propeller under different working conditions are analyzed comprehensively, and the position of propeller hot spot stress is obtained. According to the principal stress and its direction of the hot spot stress, the location and orientation of possible macroscopic crack initiation are determined. By using the linear elastic fracture mechanics and finite element method (FEM), the propeller fatigue crack propagation law and life are calculated. The analysis results show that the Forman-Newman-de Koning (FNK) model predicts more complete fatigue crack propagation life than Paris model. The fatigue crack propagation of propeller has two stages, steady propagation stage, and accelerated propagation stage. The crack breaking blade thickness is the dividing point between the two stages. In addition, the effects of the advance coefficient, stress ratio, and initial crack size on the computed results are also analyzed, which provide reference for ship driving and propeller maintenance.

Impact of stress regime change on the permeability of a naturally fractured carbonate buildup (Latemar, the Dolomites, northern Italy)

Abstract. Changing stress regimes control fracture network geometry and influence porosity and permeability in carbonate reservoirs. Using outcrop data analysis and a displacement-based linear elastic finite-element method, we investigate the impact of stress regime change on fracture network permeability. The model is based on fracture networks, specifically fracture substructures. The Latemar, predominantly affected by subsidence deformation and Alpine compression, is taken as an outcrop analogue for an isolated (Mesozoic) carbonate buildup with fracture-dominated permeability. We apply a novel strategy involving two compressive boundary loading conditions constrained by the study area's NW–SE and N–S stress directions. Stress-dependent heterogeneous apertures and effective permeability were computed in the 2D domain by (i) using the local stress state within the fracture substructure and (ii) running a single-phase flow analysis considering the fracture apertures in each fracture substructure. Our results show that the impact of the modelled far-field stresses at (i) subsidence deformation from the NW–SE and (ii) Alpine deformation from N–S increased the overall fracture aperture and permeability. In each case, increasing permeability is associated with open fractures parallel to the orientation of the loading stages and with fracture densities. The anisotropy of permeability is increased by the density and connectedness of the fracture network and affected by shear dilation. The two far-field stresses simultaneously acting within the selected fracture substructure at a different magnitude and orientation do not necessarily cancel each other out in the mechanical deformation modelling. These stresses affect the overall aperture and permeability distributions and the flow patterns. These effects – potentially ignored in simpler stress-dependent permeability – can result in significant inaccuracies in permeability estimation.

Geomechanical properties of laminated shale and bedding shale after water absorption: A case study of the Chang 7 shale in Ordos basin, China

There are abundant oil and gas resources in China's continental shale formations, and these formations often contain sandy lamina, tuffaceous lamina and carbonate lamina parallel to the bedding planes, resulting in complex geomechanical properties. During the drilling process, the mechanical weak surface structure inside the shale and the hydration effect of drilling fluid may easily cause wellbore instability, which reduces the safety and benefit of drilling. To explore the geomechanical properties of laminated shale and bedding shale after water absorption, a series of tests were conducted on the Chang 7 laminated shale and bedding shale in the Ordos Basin. The result indicates that the compressive strength and elastic modulus of laminated shale and bedding shale show a trend of decreasing first and then starting to drop as β (the angle between bedding/lamina and the direction of stress) keeps increasing, and laminated shale has stronger anisotropy, lower compressive strength and elastic modulus. The moisture content of laminated shale and bedding shale increases as the soaking time increases, which leads to a decrease in the compressive strength. Compared to bedding shale, laminated shale reaches water saturation faster and has a higher moisture content. Based on the macroscopic and microscopic images before and after shale hydration, it can be concluded that natural microfractures are easily formed at the interface between the lamina and shale matrix, and brittle minerals at the edges of natural microfractures are prone to detachment, thereby increase the hydration area and intensify the hydration process. Compared to bedding shale formations, laminated shale formations have a higher hole-size elongation ratio and are more prone to collapse according to on-site data. This study reveals the geomechanical properties of laminated shale under the dual influence of anisotropy and hydration, provides support for drilling design in such formations.

Shear-induced permeability evolution of natural fractures in granite: Implications for stimulation of EGS reservoirs

Permeability evolution of natural fractures during shearing is one of the critical issues in enhanced geothermal system (EGS). In this study, we conducted a series of numerical simulations of shear-flow tests under different normal stresses to investigate the permeability evolution and shear behavior of natural fractures during shearing. All numerical simulations in this study are executed based on a coupled hydro-mechanical pore network model (PNM) for fluid flow within the discrete element method (DEM). The numerical model uses a 33 mm × 25 mm fracture surface profile extracted from X-ray computed tomography scanning data of joints in a rock core retrieved from a 4.2-km-deep well at the Pohang EGS site. The results demonstrate a positive correlation between fracture permeability and shear displacement when the normal stress is relatively high. The central aspect of permeability evolution lies in the variation of local aperture distribution along the pathway: an increase in the variance of fracture aperture distribution leads to an increase in fracture permeability. Permeability evolution of fractures is significantly affected by normal stress, fracture roughness, and relative flow direction. Larger fracture roughness facilitates the formation of high-permeability pathways. Higher local normal stress often corresponds to lower local permeability. Moreover, high normal stress amplifies the effect of fracture roughness on permeability evolution. These findings enhance our comprehensive understanding of hydraulic-mechanical coupling effects during EGS stimulation.

Effect of true triaxial principal stress unloading rate on strain energy density of sandstone

Deep rock are often in a true triaxial stress state. Studying the impacts of varying unloading speeds on their strain energy (SE) density is highly significant for predicting rock stability. Through true triaxial unloading principal stress experiments and true triaxial stress equilibrium unloading experiments on sandstone, this paper proposes a method to compute the SE density in a true triaxial compressive unloading principal stress test. This method aims to analyze the SE variation in rocks under the action of true triaxial unloading principal stresses. Acoustic emission is used to verify the correctness of the SE density calculation method in this paper. This study found that: (1) Unloading in one principal stress direction causes the SE density to rise in the other principal stress directions. This rise in SE, depending on its reversibility, can be categorized into elastic and dissipated SE. (2)When unloading principal stresses, the released elastic SE density in the unloading direction is influence by the stress path and rate. (3) The higher the unloading speed will leads to greater increases in the input SE density, elastic SE density, and dissipative SE density in the other principal stress directions. (4) The dissipated SE generated under true triaxial compression by unloading the principal stress is positively correlated with the damage to the rock; with an increase in unloading rate, there is a corresponding increase in the formation of cracks after unloading. (5) Utilizing the stress balance unloading test, we propose a calculation method for SE density in true triaxial unloading principal stress tests.

Room temperature creep behavior of a single crystal nickel-based alloys

The deformation characteristics of a single-crystal nickel-based alloy containing Re during creep at room temperature were studied by means of creep property tests, microstructure observations and contrast analysis of the dislocation configuration. The results show that during the deformation of the alloy at room temperature, the original cubic γ′ phase transforms into a rhombus shape along the direction of maximum shearing stress, and its deformation feature is that the dislocation slips in the matrix and shears the γ′ phase. The <110> superdislocation shear into the γ′ phase can cross slip from the {111} plane to the {100} plane, forming a K-W lock, and can also be decomposed at the {111} plane. The dislocation configurations of (1/2)<110> partial dislocation plus antiphase boundary (APB) and (1/3)<112> partial dislocation plus SISF can effectively inhibit the slip and cross-slip of the dislocation and improve the deformation resistance of the alloy. At the later stage of creep, under the action of shear stress, the initial slip system is activated first to distort the γ and γ′ phases, and then the secondary slip system is activated and shears the primary slip system, resulting in a large stress concentration at the delivery point and the initiation of cracks in this area. With the alternating activation of the primary slip system and secondary slip system during creep, the initiation and expansion of cracks continue. Damage and fracture mechanisms occur in alloys during room temperature creep.

Wellbore breakouts in heavily fractured rocks: A coupled discrete fracture network-distinct element method analysis

Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely. However, the majority of such wellbore breakout analyses were based on continuum mechanics. In addition to failure in intact rocks, wellbore breakouts can also be initiated along natural discontinuities, e.g. weak planes and fractures. Furthermore, the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation. This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin, China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method (DEM) and the discrete fracture network (DFN). The developed model was validated against caliper log measurement, and its stability study was carried out by stress and displacement analyses. A parametric study was performed to investigate the effects of the characteristics of fracture distribution (orientation and length) on borehole stability by sensitivity studies. Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability. Moreover, an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress (i.e. the direction of wellbore breakouts) towards alternative directions, ultimately leading to the whole wellbore failure. These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.