A novel approach to the dynamic response analysis of Euler-Bernoulli beams resting on a Winkler soil model and subjected to impact loads

This work presents a novel approach to the dynamic response analysis of a Euler-Bernoulli beam resting on a Winkler soil model and subjected to an impact loading.The approach considers that damping has much less importance in controlling the maximum response to impulsive loadings because the maximum response is reached in a very short time,before the damping forces can dissipate a significant portion of the energy input into the system.The development of two sine series solutions,relating to different types of impulsive loadings,one involving a single concentrated force and the other a distributed line load,are presented.This study revealed that when a simply supported Euler-Bernoulli beam,resting on a Winkler soil model,is subject to an impact load,the resulting vertical displacements,bending moments and shear forces produced along the span of the beam are considerably affected.In particular,the quantification of this effect is best observed,relative to the corresponding static solution,via an amplification factor.The computed impact amplification factors,for the sub-grade moduli used in this study,were in magnitude greater than 2,thus confirming the multiple-degree-of-freedom nature of the problem.

Mechanical behavior and failure mechanisms of rock bolts subjected to static-dynamic loads

This study explores the effects of dynamic and static loading on rock bolt performance a key factor in maintaining the structural safety of coal mine roadways susceptible to coal bursts.Employing a housemade load frame to simulate various failure scenarios,pretension-impact-pull tests on rock bolts were conducted to scrutinize their dynamic responses under varied static load conditions and their failure traits under combined loads.The experimental results denote that with increased impact energy,maximum and average impact loads on rock bolts escalate significantly under pretension,initiating plastic deformation beyond a certain threshold.Despite minor reductions in the yield load due to impactinduced damage,pretension aids in constraining post-impact deformation rate and fluctuation degree of rock bolts.Moreover,impact-induced plastic deformation causes internal microstructure dislocation,fortifying the stiffness of the rock bolt support system.The magnitude of this fortification is directly related to the plastic deformation induced by the impact.These findings provide crucial guidance for designing rock bolt support in coal mine roadway excavation,emphasizing the necessity to consider both static and dynamic loads for improved safety and efficiency.

Numerical analysis of geosynthetic-reinforced embankment performance under moving loads

The performance of geosynthetic-reinforced embankments under traffic moving loads is always a hotspot in the geotechnical engineering field.A three-dimensional(3D)model of a geosynthetic-reinforced embankment without drainage consolidation was established using the finite element software ABAQUS.In this model,the traffic loads were simulated by two moving loads of rectangular pattern,and their amplitude,range,and moving speed were realized by a Fortran subroutine.The embankment fill was simulated by an equivalent linear viscoelastic model,which can reflect its viscoelasticity.The geogrid was simulated by the truss element,and the geocell was simulated by the membrane element.Infinite elements were utilized to weaken the boundary effect caused by the model geometry at the boundaries.Validation of the established numerical model was conducted by comparing the predicted deformations in the cross-section of the geosynthetic-reinforced embankment with those from the existing literature.On this basis,the dynamic stress and strain distribution in the pavement structure layer of the geosynthetic-reinforced embankment under a moving load was also analyzed.Finally,a parametric study was conducted to examine the influences of the different types of reinforcement,overload,and the moving load velocity on the geosynthetic-reinforced embankment.

Energy mechanism of bolt supporting effect to fissured rock under static and dynamic loads in deep coal mines

The stability control of fissured rock is difficult,especially under static and dynamic loads in deep coal mines.In this paper,the dynamic mechanical properties,strain rate evolution and energy dissipation of fissured and anchored rocks were respectively obtained by SHPB tests.It was found that bolt can provide supporting efficiency-improving effect for fissured rock against dynamic disturbance,and this effect increased quadratically with decrease in anchoring angles.Then,the energy dissipation mechanism of anchored rock was obtained by slipping model.Furthermore,bolt energy-absorbing mechanism by instantaneous tensile-shear deformation was expressed based on material mechanics,which was the larger the anchoring angle,the smaller the energy absorption,and the less the contribution to supporting efficiency improvement.On this basis,the functional relationship between energy dissipation of anchored rock and energy absorption of bolt was established.Taking the coal-gangue separation system of Longgu coal mine as an example,the optimal anchoring angle can be determined as 57.5°–67.5°.Field monitoring showed fissured rock with the optimal anchoring angle,can not only effectively control the deformation,but also fully exert the energy-absorbing and efficiency-improving effect of bolt itself.This study provides guidance to the stability control and supporting design for deep engineering under the same or similar conditions.

Evaluation of the dynamic sealing performance of cap rocks of underground gas storage under multi-cycle alternating loads

The static sealing of underground gas storage(UGS),including the integrity of cap rocks and the stability of faults,is analyzed from a macro perspective using a comprehensive geological evaluation method.Changes in pore structure,permeability,and mechanical strength of cap rocks under cyclic loads may impact the rock sealing integrity during the injection and recovery phases of UGS.In this work,the mechanical deformation and failure tests of rocks,as well as rock damage tests under alternating loads,are conducted to analyze the changes in the strength and permeability of rocks under multiple-cycle intense injection and recovery of UGS.Additionally,this study proposes an evaluation method for the dynamic sealing performance of UGS cap rocks under multi-cycle alternating loads.The findings suggest that the failure strength(70%)can be used as the critical value for rock failure,thus providing theoretical support for determining the upper limit of operating pressure and the number of injection-recovery cycles for the safe operation of a UGS system.

Modelling dynamic pantograph loads with combined numerical analysis

Appropriate interaction between pantograph and catenary is imperative for smooth operation of electric trains.Changing heights of overhead lines to accommodate level crossings,overbridges,and tunnels pose significant challenges in maintaining consistent current collection performance as the pantograph aerodynamic profile,and thus aerodynamic load changes significantly with operational height.This research aims to analyse the global flow characteristics and aerodynamic forces acting on individual components of an HSX pantograph operating in different configurations and orientations,such that the results can be combined with multibody simulations to obtain accurate dynamic insight into contact forces.Specifically,computational fluid dynamics simulations are used to investigate the pantograph component loads in a representative setting,such as that of the recessed cavity on a Class 800 train.From an aerodynamic perspective,this study indicates that the total drag force acting on non-fixed components of the pantograph is larger for the knuckle-leading orientation rather than the knuckle-trailing,although the difference between the two is found to reduce with increasing pantograph extension.Combining the aerodynamic loads acting on individual components with multibody tools allows for realistic dynamic insight into the pantograph behaviour.The results obtained show how considering aerodynamic forces enhance the realism of the models,leading to behaviour of the pantograph-catenary contact forces closely matching that seen in experimental tests.

Influence of railway wheel tread damage on wheel-rail impact loads and the durability of wheelsets

Dynamic wheel-rail contact forces induced by a severe form of wheel tread damage have been measured by a wheel impact load detector during full-scale field tests at different vehicle speeds.Based on laser scanning,the measured three-dimensional damage geometry is employed in simulations of dynamic vehicle-track interaction to calibrate and verify a simulation model.The relation between the magnitude of the impact load and various operational parameters,such as vehicle speed,lateral position of wheel-rail contact,track stiffness and position of impact within a sleeper bay,is investigated.The calibrated model is later employed in simulations featuring other forms of tread damage;their effects on impact load and subsequent fatigue impact on bearings,wheel webs and subsurface initiated rolling contact fatigue of the wheel tread are assessed.The results quantify the effects of wheel tread defects and are valuable in a shift towards condition-based maintenance of running gear,and for general assessment of the severity of different types of railway wheel tread damage.

Research on Demand Response Potential of Adjustable Loads in Demand Response Scenarios

To address the issues of limited demand response data,low generalization of demand response potential evaluation,and poor demand response effect,the article proposes a demand response potential feature extraction and prediction model based on data mining and a demand response potential assessment model for adjustable loads in demand response scenarios based on subjective and objective weight analysis.Firstly,based on the demand response process and demand response behavior,obtain demand response characteristics that characterize the process and behavior.Secondly,establish a feature extraction and prediction model based on data mining,including similar day clustering,time series decomposition,redundancy processing,and data prediction.The predicted values of each demand response feature on the response day are obtained.Thirdly,the predicted data of various characteristics on the response day are used as demand response potential evaluation indicators to represent different demand response scenarios and adjustable loads,and a demand response potential evaluation model based on subjective and objective weight allocation is established to calculate the demand response potential of different adjustable loads in different demand response scenarios.Finally,the effectiveness of the method proposed in the article is verified through examples,providing a reference for load aggregators to formulate demand response schemes.

Nonlinear response analysis and experimental verification for thin-walled plates to thermal-acoustic loads

For large deflection strongly nonlinear response problem of thin-walled structure to thermal-acoustic load, thermal-acoustic excitation test and corresponding simulation analysis for clamped metallic thin-walled plate have been implemented. Comparing calculated values with experimental values shows the consistency and verifies the effectiveness of calculation method and model for thin-walled plate subjected to thermal-acoustic load. Then this paper further completes dynamic response calculation for the cross reinforcement plate under different thermalacoustic load combinations. Based on the obtained time-domain displacement response, analyses about structure vibration forms are mainly focused on three typical motions of post-buckled plate,indicating that the relative strength between thermal load and acoustic load determines jump forms of plate. The Probability spectrum Density Functions（PDF） of displacement response were drawn and analyzed by employing statistical analysis method, and it clearly shows that the PDF of postbuckled plate exhibits bimodal phenomena. Then the Power Spectral Density（PSD） functions were used to analyze variations of response frequencies and corresponding peaks with the increase of temperatures, as well as how softening and hardening areas of the plate are determined. In the last section, this paper discusses the change laws of tensile stress and compressive stress in pre/post buckling areas, and gives the reasons for N glyph trend of the stress Root Mean Square（RMS）.

Fatigue life research and experimental verification of superalloy thin-walled structures subjected to thermal-acoustic loads

Rapid alternating stress is formed in structure subjected to harsh thermal-acoustic loads,which will affect fatigue performance and reduce fatigue life seriously.First,fatigue experiment of superalloy thin-walled structure was carried out to obtain fatigue damage location and failure time of the experiment specimen,and S-iN curves of superalloy thin-walled structure at 723 K were fitted.Then,dynamic response simulation of superalloy thin-walled structure under the same load as experiment was implemented,and fatigue life was estimated based on the fatigue life prediction model which mainly included:improved rain-flow counting method,Morrow average stress model and Miner linear cumulative damage theory.Further,comparisons between simulation solutions and experimental results achieved a consistency,which verified the validity of the Fatigue Life Prediction Model(FLPM).Moreover,taking a rectangle plate as the analysis object,the distributions of Fain-low circulation blocks and damage levels of the structure were discussed respectively.Finally,current research indicates that in pre-buckling the structure is in softened area and fatigue life decreases with the increase of temperature;in post-buckling the structure is in hardened area and fatigue life increases with the increase of temperature within a certain range.

Efficacy of incremental loads of cow's milk as a treatment for lactose malabsorption in Japan

BACKGROUND Lactose intolerance(LI)is commonly seen in East Asian countries.Several studies showed that lactose or milk loading has been used as a treatment for lactose malabsorption(LM)in Western countries,but there have been no reports regarding this type of treatment in Japan.As lactose or milk loading requires ingestion of large amounts of lactose within a short period,this is considered to be too harsh for Japanese people because of their less habitual milk consumption(175 mL per day in average)than Western people.In this study,we demonstrated lactose tolerance acquisition in a suitable way for Japanese.AIM To examine the efficacy of lactose(cow’s milk)loading treatment in patients with LM.METHODS Individuals with abdominal symptoms induced by milk or dairy products(LI symptoms)were identified with a questionnaire.A 20 g lactose hydrogen breath test(LHBT)was carried out to confirm LM diagnosis and to evaluate co-existence of small intestinal bacterial overgrowth(SIBO).Respondents diagnosed with LM were selected as study subjects and were treated with incremental loads of cow’s milk,starting from 30 mL and increasing up to 200 mL at 4-7 d intervals.After the treatment,changes in symptoms and LM diagnostic value of 20 g LHBT were investigated.Stool samples pre-and post-treatment were examined for changes in intestinal microbiota using 16S rRNA sequencing.Informed consent was obtained prior to each stage of the study.RESULTS In 46 subjects with LI symptoms(10-68 years old,mean age 34 years old)identified with the questionnaire,35(76.1%)were diagnosed with LM by 20 g LHBT,and 6 had co-existing SIBO.The treatment with incremental cow’s milk was carried out in 32 subjects diagnosed with LM(14-68 years old,median age 38.5 years old).The mean period of the treatment was 41±8.6 d.Improvement of symptoms was observed in 29(90.6%;95%confidence interval:75.0%-98.0%)subjects.Although 20 g LHBT indicated that 10(34.5%)subjects had improved diagnostic value of LM,no change was observed in 16(55.2%)subjects.Analysis of the fecal intestinal microbiota showed a significant increase in Blautia in 7 subjects who became symptom-free after the treatment(P=0.0313).CONCLUSION LM was diagnosed in approximately 75%of the subjects who had LI.Incremental loads of cow’s milk is regarded as a useful treatment for LM without affecting everyday life.

Numerical simulation of fatigue crack propagation in heterogeneous geomaterials under varied loads using displacement discontinuity method

Heterogeneous brittle geomaterials are highly susceptible to cyclic loads.They contain inherent flaws and cracks that grow under fatigue loads and lead to failure.This study presents a numerical model for analyzing fatigue in these materials based on the two-dimensional(2D)boundary element method and linear elastic fracture mechanics.The process is formulated by coupling the displacement discontinuity method with the incorporation technique of dissimilar regions and the governing equations of fatigue.The heterogeneous media are assumed to consist of materials with different properties,and the interfaces are assumed to be completely bonded.In addition,the domains include multiple cracks exposed to constant and variable amplitude cyclic loads.The stress intensity factor is a crucial parameter in fatigue analysis,which is determined using the displacement field around crack tips.An incremental crack growth scheme is applied to calculating the fatigue life.The growth rate values are employed to estimate the length of crack extension when there are multiple cracks.The interaction between cracks is considered,which also includes the coalescence phenomenon.Finally,various structures under different cyclic loads are examined to evaluate the accuracy of this method.The results demonstrate the efficiency of the proposed approach in modeling fatigue crack growth and life estimation.The behavior of life curves for the heterogeneous domain was as expected.These curves illustrate the breakpoints caused by utilizing discrete incremental life equations.At these points,the trend of the curves changed with the material properties and fatigue characteristics of the new material around the crack tips.

An investigation on the strain accumulation of the lightly EICP-cemented sands under cyclic traffic loads

Industrial production of chemical cement leads to extreme emissions of greenhouse gases.Biological or bioinspired sustainable materials for soil treatment projects can be employed instead of chemical cement to heal the carbon cycle in the ecosystem.The enzyme-induced calcite precipitation(EICP)method is one of the novel bio-inspired technologies that can be employed in soil treatment projects to increase desired properties of soils.While the monotonic and cyclic behavior of the enzymatically treated sands has been investigated comprehensively,the strain accumulation pattern in these improved soils under cyclic traffic loads has not been evaluated yet.In this paper,confined and unconfined cyclic compression tests are applied to the enzymatically lightly cemented sands,and the effects of the different parameters on their strain accumulation pattern are investigated for the first time in the literature.This study uses two types of specimens with unconfined compression strengths(UCS)equal to 42 kPa and 266 kPa.It is shown that the treated specimens have a rate-dependent behavior where cyclic loads with low frequencies lead to more resilient and plastic strains in the specimens.The results show that by approaching the maximum applied stresses to the UCS of the specimens(by breaking more calcite bonds between sand particles),the rate dependency behavior of specimens will reduce.Investigation of the effects of the cementation level demonstrated that by increasing the amount of the precipitated calcite from 0.38%to 0.83%,accumulated plastic strains are reduced almost 95%under the same loading condition.Effects of the initial static loads,confining pressures,the number of cycles,and amplitudes of the cyclic loads are also evaluated.

Mechanical responses of anchoring structure under triaxial cyclic loading

Dynamic load on anchoring structures(AS)within deep roadways can result in cumulative damage and failure.This study develops an experimental device designed to test AS under triaxial loads.The device enables the investigation of the mechanical response,failure mode,instability assessment criteria,and anchorage effect of AS subjected to combined cyclic dynamic-static triaxial stress paths.The results show that the peak bearing strength is positively correlated with the anchoring matrix strength,anchorage length,and edgewise compressive strength.The bearing capacity decreases significantly when the anchorage direction is severely inclined.The free face failure modes are typically transverse cracking,concave fracturing,V-shaped slipping and detachment,and spallation detachment.Besides,when the anchoring matrix strength and the anchorage length decrease while the edgewise compressive strength,loading rate,and anchorage inclination angle increase,the failure intensity rises.Instability is determined by a negative tangent modulus of the displacement-strength curve or the continued deformation increase against the general downward trend.Under cyclic loads,the driving force that breaks the rock mass along the normal vector and the rigidity of the AS are the two factors that determine roadway stability.Finally,a control measure for surrounding rock stability is proposed to reduce the internal driving force via a pressure relief method and improve the rigidity of the AS by full-length anchorage and grouting modification.

A drug-loaded flexible substrate improves the performance of conformal cortical electrodes

Cortical electrodes are a powerful tool for the stimulation and/or recording of electrical activity in the nervous system.However,the inevitable wound caused by surgical implantation of electrodes presents bacterial infection and inflammatory reaction risks associated with foreign body exposure.Moreover,inflammation of the wound area can dramatically worsen in response to bacterial infection.These consequences can not only lead to the failure of cortical electrode implantation but also threaten the lives of patients.Herein,we prepared a hydrogel made of bacterial cellulose(BC),a flexible substrate for cortical electrodes,and further loaded antibiotic tetracycline(TC)and the anti-inflammatory drug dexamethasone(DEX)onto it.The encapsulated drugs can be released from the BC hydrogel and effectively inhibit the growth of Gram-negative and Gram-positive bacteria.Next,therapeutic cortical electrodes were developed by integrating the drug-loaded BC hydrogel and nine-channel serpentine arrays;these were used to record electrocorticography(ECoG)signals in a rat model.Due to the controlled release of TC and DEX from the BC hydrogel substrate,therapeutic cortical electrodes can alleviate or prevent symptoms associated with the bacterial infection and inflammation of brain tissue.This approach facilitates the development of drug delivery electrodes for resolving complications caused by implantable electrodes.

Reinforcing effects of polypropylene on energy absorption and fracturing of cement-based tailings backfill under impact loading

Polypropylene(PP)fiber-reinforced cement-based tailings backfill(FRCTB)is a green compound material with superior crack resistance and has good prospects for application in underground mining.However,FRCTB exhibits susceptibility to dynamic events,such as impact ground pressure and blast vibrations.This paper investigates the energy and crack distribution behavior of FRCTB under dynamic impact,considering the height/diameter(H/D)effect.Split Hopkinson pressure bar,industrial computed tomography scan,and scanning electron microscopy(SEM)experiments were carried out on six types of FRCTB.Laboratory outcomes confirmed fiber aggregation at the bottom of specimens.When H/D was less than 0.8,the proportion of PP fibers distributed along theθangle direction of80°-90°increased.For the total energy,all samples presented similar energy absorption,reflectance,and transmittance.However,a rise in H/D may cause a rise in the energy absorption rate of FRCTB during the peak phase.A positive correlation existed between the average strain rate and absorbed energy per unit volume.The increase in H/D resulted in a decreased crack volume fraction of FRCTB.When the H/D was greater than or equal to 0.7,the maximum crack volume fraction of FRCTB was observed close to the incidence plane.Radial cracks were present only in the FRCTB with an H/D ratio of 0.5.Samples with H/D ratios of 0.5 and 0.6 showed similar distributions of weakly and heavily damaged areas.PP fibers can limit the emergence and expansion of cracks by influencing their path.SEM observations revealed considerable differences in the bonding strengths between fibers and the FRCTB.Fibers that adhered particularly well to the substrate were attracted together with the hydration products adhering to surfaces.These results show that FRCTB is promising as a sustainable and green backfill for determining the design properties of mining with backfill.

Conservation genomic investigation of an endangered conifer,Thuja sutchuenensis,reveals low genetic diversity but also low genetic load

Endangered species generally have small populations with low genetic diversity and a high genetic load.Thuja sutchuenensis is an endangered conifer endemic to southwestern China.It was once considered extinct in the wild,but in 1999 was rediscovered.However,little is known about its genetic load.We collected 67 individuals from five wild,isolated T.sutchuenensis populations,and used 636,151 SNPs to analyze the level of genetic diversity and genetic load in T.sutchuenensis to delineate the conservation units of T.sutchuenensis,based on whole transcriptome sequencing data,as well as target capture sequencing data.We found that populations of T.sutchuenensis could be divided into three groups.These groups had low levels genetic diversity and were moderately genetically differentiated.Our findings also indicate that T.sutchuenensis suffered two severe bottlenecks around the Last Glaciation Period and Last Glacial Maximum.Among Thuja species,T.sutchuenensis presented the lowest genetic load and hence might have purged deleterious mutations efficiently through purifying selection.However,distribution of fitness effects analysis indicated a high extinction risk for T.sutchuenensis.Multiple lines of evidence identified three management units for T.sutchuenensis.Although T.sutchuenensis possesses a low genetic load,low genetic diversity,suboptimal fitness,and anthropogenic pressures all present an extinction risk for this rare conifer.This might also hold true for many endangered plant species in the mountains all over the world.

Damage evolution of rock-encased-backfill structure under stepwise cyclic triaxial loading

Rock-encased-backfill(RB)structures are common in underground mining,for example in the cut-andfill and stoping methods.To understand the effects of cyclic excavation and blasting activities on the damage of these RB structures,a series of triaxial stepwise-increasing-amplitude cyclic loading experiments was conducted with cylindrical RB specimens(rock on outside,backfill on inside)with different volume fractions of rock(VF=0.48,0.61,0.73,and 0.84),confining pressures(0,6,9,and 12 MPa),and cyclic loading rates(200,300,400,and 500 N/s).The damage evolution and meso-crack formation during the cyclic tests were analyzed with results from stress-strain hysteresis loops,acoustic emission events,and post-failure X-ray 3D fracture morphology.The results showed significant differences between cyclic and monotonic loadings of RB specimens,particularly with regard to the generation of shear microcracks,the development of stress memory and strain hardening,and the contact forces and associated friction that develops along the rock-backfill interface.One important finding is that as a function of the number of cycles,the elastic strain increases linearly and the dissipated energy increases exponentially.Also,compared with monotonic loading,the cyclic strain hardening characteristics are more sensitive to rising confining pressures during the initial compaction stage.Another finding is that compared with monotonic loading,more shear microcracks are generated during every reloading stage,but these microcracks tend to be dispersed and lessen the likelihood of large shear fracture formation.The transition from elastic to plastic behavior varies depending on the parameters of each test(confinement,volume fraction,and cyclic rate),and an interesting finding was that the transformation to plastic behavior is significantly lower under the conditions of 0.73 rock volume fraction,400 N/s cyclic loading rate,and 9 MPa confinement.All the findings have important practical implications on the ability of backfill to support underground excavations.

Liquefaction susceptibility and deformation characteristics of saturated coral sandy soils subjected to cyclic loadings-a critical review

Coral sandy soils widely exist in coral island reefs and seashores in tropical and subtropical regions.Due to the unique marine depositional environment of coral sandy soils,the engineering characteristics and responses of these soils subjected to monotonic and cyclic loadings have been a subject of intense interest among the geotechnical and earthquake engineering communities.This paper critically reviews the progress of experimental investigations on the undrained behavior of coral sandy soils under monotonic and cyclic loadings over the last three decades.The focus of coverage includes the contractive-dilative behavior,the pattern of excess pore-water pressure(EPWP)generation and the liquefaction mechanism and liquefaction resistance,the small-strain shear modulus and strain-dependent shear modulus and damping,the cyclic softening feature,and the anisotropic characteristics of undrained responses of saturated coral sandy soils.In particular,the advances made in the past decades are reviewed from the following aspects:(1)the characterization of factors that impact the mechanism and patterns of EPWP build-up;(2)the identification of liquefaction triggering in terms of the apparent viscosity and the average flow coefficient;(3)the establishment of the invariable form of strain-based,stress-based,or energy-based EPWP ratio formulas and the unique relationship between the new proxy of liquefaction resistance and the number of cycles required to reach liquefaction;(4)the establishment of the invariable form of the predictive formulas of small strain modulus and strain-dependent shear modulus;and(5)the investigation on the effects of stress-induced anisotropy on liquefaction susceptibility and dynamic deformation characteristics.Insights gained through the critical review of these advances in the past decades offer a perspective for future research to further resolve the fundamental issues concerning the liquefaction mechanism and responses of coral sandy sites subjected to cyclic loadings associated with seismic events in marine environments.

Strength and failure characteristics of marble spheres subjected to paired point loads

Failure of irregular rock samples may provide implications in the rapid estimation of rock strength,which is imperative in rock engineering practice.In this work,analytical,experimental and numerical investigations were carried out to study the mechanical properties and failure characteristics of rock spheres under paired point loads.Analytical solutions indicted that with the increase in sample size(contact angle)and decrease in Poisson’s ratio,the uneven tensile stress in theta direction decreased.Then laboratory experiments were carried out to investigate the load characteristics and failure mode of spherical marble samples with different sizes subjected to a pair of diametral point loads.The discrete element method(DEM)was adopted to study the failure process of rock spheres.The effect of the sphere diameter on the point load contact angle was examined in terms of peak load,crushed zone distribution and energy dissipation.Experimental and numerical results showed that the samples primarily fail in tension,with crushed zones formed at both loading points.With increase in the sample size,the contact angle,crushed area and total work increase.As the specimen diameter increases from 30 mm to 50 mm,the peak load on the specimen increases from 3.6 kN to 8.8 kN,and the percentage of crushed zone(ratio of crushing zone to sample radius,d/r)increased from 0.191 to 0.262.The results of the study have implications for understanding the failure of irregular rock specimens under point loading conditions and their size effects.