Three-dimensional(3D)parametric measurements of individual gravels in the Gobi region using point cloud technique

Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materials constituting the Gobi result in notable differences in saltation processes across various Gobi surfaces.It is challenging to describe these processes according to a uniform morphology.Therefore,it becomes imperative to articulate surface characteristics through parameters such as the three-dimensional(3D)size and shape of gravel.Collecting morphology information for Gobi gravels is essential for studying its genesis and sand saltation.To enhance the efficiency and information yield of gravel parameter measurements,this study conducted field experiments in the Gobi region across Dunhuang City,Guazhou County,and Yumen City(administrated by Jiuquan City),Gansu Province,China in March 2023.A research framework and methodology for measuring 3D parameters of gravel using point cloud were developed,alongside improved calculation formulas for 3D parameters including gravel grain size,volume,flatness,roundness,sphericity,and equivalent grain size.Leveraging multi-view geometry technology for 3D reconstruction allowed for establishing an optimal data acquisition scheme characterized by high point cloud reconstruction efficiency and clear quality.Additionally,the proposed methodology incorporated point cloud clustering,segmentation,and filtering techniques to isolate individual gravel point clouds.Advanced point cloud algorithms,including the Oriented Bounding Box(OBB),point cloud slicing method,and point cloud triangulation,were then deployed to calculate the 3D parameters of individual gravels.These systematic processes allow precise and detailed characterization of individual gravels.For gravel grain size and volume,the correlation coefficients between point cloud and manual measurements all exceeded 0.9000,confirming the feasibility of the proposed methodology for measuring 3D parameters of individual gravels.The proposed workflow yields accurate calculations of relevant parameters for Gobi gravels,providing essential data support for subsequent studies on Gobi environments.

Three-dimensional numerical analysis of plant-soil hydraulic interactions on pore water pressure of vegetated slope under different rainfall patterns

Understanding the pore water pressure distribution in unsaturated soil is crucial in predicting shallow landslides triggered by rainfall,mainly when dealing with different temporal patterns of rainfall intensity.However,the hydrological response of vegetated slopes,especially three-dimensional(3D)slopes covered with shrubs,under different rainfall patterns remains unclear and requires further investigation.To address this issue,this study adopts a novel 3D numerical model for simulating hydraulic interactions between the root system of the shrub and the surrounding soil.Three series of numerical parametric studies are conducted to investigate the influences of slope inclination,rainfall pattern and rainfall duration.Four rainfall patterns(advanced,bimodal,delayed,and uniform)and two rainfall durations(4-h intense and 168-h mild rainfall)are considered to study the hydrological response of the slope.The computed results show that 17%higher transpiration-induced suction is found for a steeper slope,which remains even after a short,intense rainfall with a 100-year return period.The extreme rainfalls with advanced(PA),bimodal(PB)and uniform(PU)rainfall patterns need to be considered for the short rainfall duration(4 h),while the delayed(PD)and uniform(PU)rainfall patterns are highly recommended for long rainfall durations(168 h).The presence of plants can improve slope stability markedly under extreme rainfall with a short duration(4 h).For the long duration(168 h),the benefit of the plant in preserving pore-water pressure(PWP)and slope stability may not be sufficient.

Three-dimensional stability calculation method for high and large composite slopes formed by mining stope and inner dump in adjacent open pits

The 2D limit equilibrium method is widely used for slope stability analysis.However,with the advancement of dump engineering,composite slopes often exhibit significant 3D mechanical effects.Consequently,it is of significant importance to develop an effective 3D stability calculation method for composite slopes to enhance the design and stability control of open-pit slope engineering.Using the composite slope formed by the mining stope and inner dump in Baiyinhua No.1 and No.2 open-pit coal mine as a case study,this research investigates the failure mode of composite slopes and establishes spatial shape equations for the sliding mass.By integrating the shear resistance and sliding force of each row of microstrip columns onto the bottom surface of the strip corresponding to the main sliding surface,a novel 2D equivalent physical and mechanical parameters analysis method for the strips on the main sliding surface of 3D sliding masses is proposed.Subsequently,a comprehensive 3D stability calculation method for composite slopes is developed,and the quantitative relationship between the coordinated development distance and its 3D stability coefficients is examined.The analysis reveals that the failure mode of the composite slope is characterized by cutting-bedding sliding,with the arc serving as the side interface and the weak layer as the bottom interface,while the destabilization mechanism primarily involves shear failure.The spatial form equation of the sliding mass comprises an ellipsoid and weak plane equation.The analysis revealed that when the coordinated development distance is 1500 m,the error rate between the 3D stability calculation result and the 2D stability calculation result of the composite slope is less than 8%,thereby verifying the proposed analytical method of equivalent physical and mechanical parameters and the 3D stability calculation method for composite slopes.Furthermore,the3D stability coefficient of the composite slope exhibits an exponential correlation with the coordinated development distance,with the coefficient gradually decreasing as the coordinated development distance increases.These findings provide a theoretical guideline for designing similar slope shape parameters and conducting stability analysis.

Phase II Clinical Study of Three-Dimensional Printed Coplanar Template Combined with CT-Guided Percutaneous Core Needle Biopsy of Pulmonary Nodules in Elderly Patients

Background: As the population age structure gradually ages, more and more elderly people were found to have pulmonary nodules during physical examinations. Most elderly people had underlying diseases such as heart, lung, brain and blood vessels and cannot tolerate surgery. Computed tomography (CT)-guided percutaneous core needle biopsy (CNB) was the first choice for pathological diagnosis and subsequent targeted drugs, immune drugs or ablation treatment. CT-guided percutaneous CNB requires clinicians with rich CNB experience to ensure high CNB accuracy, but it was easy to cause complications such as pneumothorax and hemorrhage. Three-dimensional (3D) printing coplanar template (PCT) combined with CT-guided percutaneous pulmonary CNB biopsy has been used in clinical practice, but there was no prospective, randomized controlled study. Methods: Elderly patients with lung nodules admitted to the Department of Oncology of our hospital from January 2019 to January 2023 were selected. A total of 225 elderly patients were screened, and 30 patients were included after screening. They were randomly divided into experimental group (Group A: 30 cases) and control group (Group B: 30 cases). Group A was given 3D-PCT combined with CT-guided percutaneous pulmonary CNB biopsy, Group B underwent CT-guided percutaneous pulmonary CNB. The primary outcome measure of this study was the accuracy of diagnostic CNB, and the secondary outcome measures were CNB time, number of CNB needles, number of pathological tissues and complications. Results: The diagnostic accuracy of group A and group B was 96.67% and 76.67%, respectively (P = 0.026). There were statistical differences between group A and group B in average CNB time (P = 0.001), number of CNB (1 vs more than 1, P = 0.029), and pathological tissue obtained by CNB (3 vs 1, P = 0.040). There was no statistical difference in the incidence of pneumothorax and hemorrhage between the two groups (P > 0.05). Conclusions: 3D-PCT combined with CT-guided percutaneous CNB can improve the puncture accuracy of elderly patients, shorten the puncture time, reduce the number of punctures, and increase the amount of puncture pathological tissue, without increasing pneumothorax and hemorrhage complications. We look forward to verifying this in a phase III randomized controlled clinical study. .

Three-dimensional gravity inversion based on 3D U-Net++

The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity inversion method based on 3D U-Net++.Compared with two-dimensional gravity inversion,three-dimensional(3D)gravity inversion can more precisely describe the density distribution of underground space.However,conventional 3D gravity inversion method input is two-dimensional,the input and output of the network proposed in our method are three-dimensional.In the training stage,we design a large number of diversifi ed simulation model-data pairs by using the random walk method to improve the generalization ability of the network.In the test phase,we verify the network performance by using the model-data pairs generated by the simulation.To further illustrate the eff ectiveness of the algorithm,we apply the method to the inversion of the San Nicolas mining area,and the inversion results are basically consistent with the borehole measurement results.Moreover,the results of the 3D U-Net++inversion and the 3D U-Net inversion are compared.The density models of the 3D U-Net++inversion have higher resolution,more concentrated inversion results,and a clearer boundary of the density model.

Three-Dimensional Mesoscopic Investigation on the Impact of Specimen Geometry and Bearing Strip Size on the Splitting-Tensile Properties of Coral Aggregate Concrete

The use of coral aggregate concrete(CAC)as a novel construction material has attracted significant attention for the construction of reef engineering structures.To investigate the static splitting-tensile behaviors of CAC under the influence of two factors,namely specimen geometry and bearing strip size,a three-dimensional(3D)mesoscale modeling approach with consideration for aggregate randomness in shape and distribution was adopted in this study.We established 12 different specimen models with two specimen shapes(i.e.,a cube with an edge length of 150 mm and a cylinder with dimensions ofφ150 mm×300 mm)and six strip widths(i.e.,6,9,12,15,18,and 20 mm)for calculation.The effects of specimen geometry and strip width on the splitting-tensile properties of CAC,such as failure processes,final failure patterns,and splitting-tensile strength(fst),are analyzed and discussed systematically.The results indicate the high reliability of the developed mesoscale modeling approach and reveal the optimal computational parameters for simulating and predicting the splitting-tensile properties of CAC.The fstvalues of CAC are associated with both the specimen geometry and width of the bearing strip.The fstvalues of the cube model are slightly higher than those of the cylinder model for the same bearing strip size,representing geometry effects that can be explained by differences in fracture area.Additionally,the fstvalue of CAC gradually increases with the relative width of the bearing strip ranging from 0.04 to 0.13.Based on the elastic solution theory,the variation area of CAC fstvalues with the relative width of the bearing strip was determined preliminarily,which has great significance for studying the tensile performance of CAC.

A Three-Dimensional (3D) Environment to Maintain the Integrity of Mouse Testicular Can Cause the Occurrence of Meiosis

Adhesions between different cells and extracellular matrix have been studied extensively in vitro, but little is known about their functions in testicular tissue counterparts. Spermatogonia and their companion somatic cells maintain a close association throughout spermatogenesis and this association is necessary for normal spermatogenesis. In order to keep the relative integrity of the testicular tissues, and to detect the development in vitro, culture testicular tissues in a three- dimensional （3D） agarose matrix was examined. Testicular tissues isolated from 6.5 d postpartum （dpp） mouse were cultured on the top of the matrix for 26 d with a medium height up to 4/5 of the 3D agarose matrix. The results showed that in this 3D culture environment, each type of testicular cells kept the same structure, localization and function as in vivo and might be more biologically relevant to living organisms. After culture, germ cell marker VASA and meiosis markers DAZL and SCP3 showed typical positive analysed by immunofluorescence staining and RT-PCR. It demonstrated that this 3D culture system was able to maintain the number of germ cells and promote the meiosis initiation of male germ cells.

A Three-Dimensional Numerical Model for Predicting the Weld Bead Geometry Characteristics in Laser Overlap Welding of Low Carbon Galvanized Steel

Laser welding (LW) becomes one of the most economical high quality joining processes. LW offers the advantage of very controlled heat input resulting in low distortion and the ability to weld heat sensitive components. To exploit efficiently the benefits presented by LW, it is necessary to develop an integrated approach to identify and control the welding process variables in order to produce the desired weld characteristics without being forced to use the traditional and fastidious trial and error procedures. The paper presents a study of weld bead geometry characteristics prediction for laser overlap welding of low carbon galvanized steel using 3D numerical modelling and experimental validation. The temperature dependent material properties, metallurgical transformations and enthalpy method constitute the foundation of the proposed modelling approach. An adaptive 3D heat source is adopted to simulate both keyhole and conduction mode of the LW process. The simulations are performed using 3D finite element model on commercial software. The model is used to estimate the weld bead geometry characteristics for various LW parameters, such as laser power, welding speed and laser beam diameter. The calibration and validation of the 3D numerical model are based on experimental data achieved using a 3 kW Nd:Yag laser system, a structured experimental design and confirmed statistical analysis tools. The results reveal that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects on the weld quality. The results show great concordance between predicted and measured values for weld bead geometry characteristics, such as depth of penetration, bead width at the top surface and bead width at the interface between sheets, with an average accuracy greater than 95%.

Fabrication of three-dimensional islet models by geometry-controlled hanging drop method

Hanging-drop method has been widely used to fabricate three-dimensional (3D) in vitro tissue models due to its advantages such as being easy to perform, inexpensive, and permitting precise control of cell spheroid formation. The geometry of hanging drop may play a critical role on the formation of cell spheroids, which, however, has not been explored. In this study, we developed a modified hanging-drop platform that enables the production of cell spheroids in a high-throughput manner by controlling hanging drop geometry with defined spreading ring. The surface tension force is proportional to the spreading ring and gravitational force is determined by the droplet volume, and the geometry can be determined by the balance between surface tension and gravity.β-TC-6 cell spheroids with optimized diameters were fabricated as 3D in vitro islet models. The models show morphology similar to primary islets and have functionality that more closely resembles primary islets than two-dimensional cell culture. The developed platform holds great potential for engineering well-controlled in vitro tissue models for various applications such as physiological and pathological studies, drug screening, as well as transplantation for treatment purpose.

A COMPUTER CODE FASTOR-3D FOR TRANSIENT THREE-DIMENSIONAL SIMULATION ON FLOW AND HEAT TRANS-FER WITH POROSITY APPROXIMATION

Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer. These equations were solved numerically with the finite differential method and the primitive variable approach. This method uses staggered grid and pressure correction schemes. A computer code FASTOR3D integrated the aforementioned algorithm. The preliminary results have been compared with conventional benchmark solutions. With auxiliary software DV, the numerical results were visualized in colorful images to demonstrate the variation of flow patterns and temperature profiles during the transient process. The results of the simulation code for the fluid flows and heat transfer in the sodium pool of a fast breeder reactor are acceptable.

Resolution assessment of the three-dimensional acquisition geometry of the Anhui experiment in"Geoscience Yangtze Plan"

According to the actual observation conditions of the Yangtze River valley from Anqing city to Maanshan city,we designed the 3D acquisition geometry,and applied the multi-scale checkerboard semblance analysis to assess the preliminary resolution of the designed observation.The checkerboard semblance tests use the refraction and reflection travel-time simultaneous inversion algorithm to quantitatively provide both resolution assessment of velocity structure and Moho topography.The multi-scale checker-board semblance recovery results show that while the chec-kerboard semblance threshold value is 0.5,the preliminary resolution of the designed acquisition geometry is better than 10 km in the upper crust(the depth is less than 10 km),around 15 km in the mid-crust(the depth is 10?25 km),and better than 20 km in the lower crust(the depth is 25?33 km).The preliminary tomographic resolution for the Moho topography is about 20 km in the ray-path coverage area beneath the acquisition geometry.While the checkerboard semblance threshold value is 0.75,the preliminary resolution is 20 km in the upper crust,around 20?25 km in the mid-crust and 25 km in the lower crust.And the preliminary tomographic resolution for the Moho topography is better than 30 km in the ray-path coverage area beneath the acquisition geometry.These non-linear checkerboard tests reveal that the designed acquisition geometry is suitable to image the crustal velocity structure of the Yangtze River valley in the Anhui province.

3D DEM simulation of hard rock fracture in deep tunnel excavation induced by changes in principal stress magnitude and orientation

To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.

Modeling time-dependent mechanical behavior of hard rock considering excavation-induced damage and complex 3D stress states

To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main components:a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state,a quantitative relationship between critical irreversible deformation and complex stress state,and evolution characteristics of strength parameters.The proposed model is implemented in a self-developed numerical code,i.e.CASRock.The reliability of the model is validated through experiments.It is indicated that the time-dependent fracturing potential index(xTFPI)at a given time during the attenuation creep stage shows a negative correlation with the extent of excavationinduced damage.The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress,thereby highlighting the 3D stress-dependent characteristic of the model.Finally,the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed.

Three-Dimensional Ocean Sensor Networks:A Survey

The past decade has seen a growing interest in ocean sensor networks because of their wide applications in marine research,oceanography,ocean monitoring,offshore exploration,and defense or homeland security.Ocean sensor networks are generally formed with various ocean sensors,autonomous underwater vehicles,surface stations,and research vessels.To make ocean sensor network applications viable,efficient communication among all devices and components is crucial.Due to the unique characteristics of underwater acoustic channels and the complex deployment environment in three dimensional(3D) ocean spaces,new efficient and reliable communication and networking protocols are needed in design of ocean sensor networks.In this paper,we aim to provide an overview of the most recent advances in network design principles for 3D ocean sensor networks,with focuses on deployment,localization,topology design,and position-based routing in 3D ocean spaces.

A state-of-the-art review of automated extraction of rock mass discontinuity characteristics using three-dimensional surface models

In the last two decades,significant research has been conducted in the field of automated extraction of rock mass discontinuity characteristics from three-dimensional(3D)models.This provides several methodologies for acquiring discontinuity measurements from 3D models,such as point clouds generated using laser scanning or photogrammetry.However,even with numerous automated and semiautomated methods presented in the literature,there is not one single method that can automatically characterize discontinuities accurately in a minimum of time.In this paper,we critically review all the existing methods proposed in the literature for the extraction of discontinuity characteristics such as joint sets and orientations,persistence,joint spacing,roughness and block size using point clouds,digital elevation maps,or meshes.As a result of this review,we identify the strengths and drawbacks of each method used for extracting those characteristics.We found that the approaches based on voxels and region growing are superior in extracting joint planes from 3D point clouds.Normal tensor voting with trace growth algorithm is a robust method for measuring joint trace length from 3D meshes.Spacing is estimated by calculating the perpendicular distance between joint planes.Several independent roughness indices are presented to quantify roughness from 3D surface models,but there is a need to incorporate these indices into automated methodologies.There is a lack of efficient algorithms for direct computation of block size from 3D rock mass surface models.

Two kinds of contact problems in three-dimensional icosahedral quasicrystals

Two kinds of contact problems, i.e., the frictional contact problem and the adhesive contact problem, in three-dimensional （3D） icosahedral quasicrystals are dis- cussed by a complex variable function method. For the frictional contact problem, the contact stress exhibits power singularities at the edge of the contact zone. For the adhe- sive contact problem, the contact stress exhibits oscillatory singularities at the edge of the contact zone. The numerical examples show that for the two kinds of contact problems, the contact stress exhibits singularities, and reaches the maximum value at the edge of the contact zone. The phonon-phason coupling constant has a significant effect on the contact stress intensity, while has little impact on the contact stress distribution regu- lation. The results are consistent with those of the classical elastic materials when the phonon-phason coupling constant is 0. For the adhesive contact problem, the indentation force has positive correlation with the contact displacement, but the phonon-phason cou- pling constant impact is barely perceptible. The validity of the conclusions is verified.

Effects of dry-wet cycles on three-dimensional pore structure and permeability characteristics of granite residual soil using X-ray micro computed tomography

Due to seasonal climate alterations,the microstructure and permeability of granite residual soil are easily affected by multiple dry-wet cycles.The X-ray micro computed tomography(micro-CT)acted as a nondestructive tool for characterizing the microstructure of soil samples exposed to a range of damage levels induced by dry-wet cycles.Subsequently,the variations of pore distribution and permeability due to drywet cycling effects were revealed based on three-dimensional(3D)pore distribution analysis and seepage simulations.According to the results,granite residual soils could be separated into four different components,namely,pores,clay,quartz,and hematite,from micro-CT images.The reconstructed 3D pore models dynamically demonstrated the expanding and connecting patterns of pore structures during drywet cycles.The values of porosity and connectivity are positively correlated with the number of dry-wet cycles,which were expressed by exponential and linear functions,respectively.The pore volume probability distribution curves of granite residual soil coincide with the χ^(2)distribution curve,which verifies the effectiveness of the assumption of χ^(2)distribution probability.The pore volume distribution curves suggest that the pores in soils were divided into four types based on their volumes,i.e.micropores,mesopores,macropores,and cracks.From a quantitative and visual perspective,considerable small pores are gradually transformed into cracks with a large volume and a high connectivity.Under the action of dry-wet cycles,the number of seepage flow streamlines which contribute to water permeation in seepage simulation increases distinctly,as well as the permeability and hydraulic conductivity.The calculated hydraulic conductivity is comparable with measured ones with an acceptable error margin in general,verifying the accuracy of seepage simulations based on micro-CT results.

Three-dimensional multi-constraint route planning of unmanned aerial vehicle low-altitude penetration based on coevolutionary multi-agent genetic algorithm

To address the issue of premature convergence and slow convergence rate in three-dimensional (3D) route planning of unmanned aerial vehicle (UAV) low-altitude penetration,a novel route planning method was proposed.First and foremost,a coevolutionary multi-agent genetic algorithm (CE-MAGA) was formed by introducing coevolutionary mechanism to multi-agent genetic algorithm (MAGA),an efficient global optimization algorithm.A dynamic route representation form was also adopted to improve the flight route accuracy.Moreover,an efficient constraint handling method was used to simplify the treatment of multi-constraint and reduce the time-cost of planning computation.Simulation and corresponding analysis show that the planning results of CE-MAGA have better performance on terrain following,terrain avoidance,threat avoidance (TF/TA2) and lower route costs than other existing algorithms.In addition,feasible flight routes can be acquired within 2 s,and the convergence rate of the whole evolutionary process is very fast.

A simplified three-dimensional extension of Hoek-Brown strength criterion

The Hoek-Brown(HB)strength criterion has been applied widely in a large number of projects around the world.However,this criterion ignores the intermediate principal stress s2.Many evidences have demonstrated that the rock strength is dependent on s2.Thus it is necessary to extend the HB criterion into a three-dimensional(3D)form.In this study,the effect of s2 on the strength of rocks is identified by reviewing the true triaxial tests of various rock types reported in the literature.A simple 3D strength criterion is developed.The modified criterion is verified by the true triaxial tests of 13 rock types.The results indicate that the modified criterion can achieve a good fit to most of rock types.It can represent a series of criteria as b varies.For comparisons,several existing 3D versions of the HB criterion are selected to predict the strengths of these rock types.It is indicated that the proposed criterion works better than other criteria.A substantial relationship between parameter b and the unconfined compressive strength is established,which guarantees that the proposed criterion can still work well even in the absence of true triaxial test data.

Effect of fractures on mechanical behavior of sand powder 3D printing rock analogue under triaxial compression

In practical engineering applications,rock mass are often found to be subjected to a triaxial stress state.Concurrently,defects like joints and fractures have a notable impact on the mechanical behavior of rock mass.Such defects are identified as crucial contributors to the failure and instability of the surrounding rock,subsequently impacting the engineering stability.The study aimed to investigate the impact of fracture geometry and confining pressure on the deformation,failure characteristics,and strength of specimens using sand powder 3D printing technology and conventional triaxial compression tests.The results indicate that the number of fractures present considerably influences the peak strength,axial peak strain and elastic modulus of the specimens.Confining pressure is an important factor affecting the failure pattern of the specimen,under which the specimen is more prone to shear failure,but the initiation,expansion and penetration processes of secondary cracks in different fracture specimens are different.This study confirmed the feasibility of using sand powder 3D printing specimens as soft rock analogs for triaxial compression research.The insights from this research are deemed essential for a deeper understanding of the mechanical behavior of fractured surrounding rocks when under triaxial stress state.