Shear behavior of intact granite under thermo-mechanical coupling and three-dimensional morphology of shear-formed fractures

The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.

A thermo-mechanical damage constitutive model for deep rock considering brittleness-ductility transition characteristics

This paper developed a statistical damage constitutive model for deep rock by considering the effects of external load and thermal treatment temperature based on the distortion energy.The model parameters were determined through the extremum features of stress−strain curve.Subsequently,the model predictions were compared with experimental results of marble samples.It is found that when the treatment temperature rises,the coupling damage evolution curve shows an S-shape and the slope of ascending branch gradually decreases during the coupling damage evolution process.At a constant temperature,confining pressure can suppress the expansion of micro-fractures.As the confining pressure increases the rock exhibits ductility characteristics,and the shape of coupling damage curve changes from an S-shape into a quasi-parabolic shape.This model can well characterize the influence of high temperature on the mechanical properties of deep rock and its brittleness-ductility transition characteristics under confining pressure.Also,it is suitable for sandstone and granite,especially in predicting the pre-peak stage and peak stress of stress−strain curve under the coupling action of confining pressure and high temperature.The relevant results can provide a reference for further research on the constitutive relationship of rock-like materials and their engineering applications.

Study on inhomogeneous cooling behavior of extruded profile with unequal and large thicknesses during quenching using thermo-mechanical coupling model

The interfacial heat transfer coefficient between hot profile surface and cooling water was determined by using inverse heat conduction model combined with end quenching experiment. Then, a Deform-3 D thermo-mechanical coupling model for simulating the on-line water quenching of extruded profile with unequal and large thicknesses was developed. The temperature field, residual stress field and distortion of profile during quenching were investigated systematically. The results show that heat transfer coefficient increases as water flow rate increases. The peak heat transfer coefficient with higher water flow rates appears at lower interface temperatures. The temperature distribution across the cross-section of profile during quenching is severe nonuniform and the maximum temperature difference is 300 ℃ at quenching time of 3.49 s. The temperature difference through the thickness of different parts of profile first increases sharply to a maximum value, and then gradually decreases. The temperature gradient increases obviously with the increase of thickness of parts. After quenching, there exist large residual stresses on the inner side of joints of profile and the two ends of part with thickness of 10 mm. The profile presents a twisting-type distortion across the cross-section under non-uniform cooling and the maximum twisting angle during quenching is 2.78°.

Multiscale Nonlinear Thermo-Mechanical Coupling Analysis of Composite Structures with Quasi-Periodic Properties

This paper reports a multiscale analysis method to predict the thermomechanical coupling performance of composite structures with quasi-periodic properties.In these material structures,the configurations are periodic,and the material coefficients are quasi-periodic,i.e.,they depend not only on the microscale information but also on the macro location.Also,a mutual interaction between displacement and temperature fields is considered in the problem,which is our particular interest in this study.The multiscale asymptotic expansions of the temperature and displacement fields are constructed and associated error estimation in nearly pointwise sense is presented.Then,a finite element-difference algorithm based on the multiscale analysis method is brought forward in detail.Finally,some numerical examples are given.And the numerical results show that the multiscale method presented in this paper is effective and reliable to study the nonlinear thermo-mechanical coupling problem of composite structures with quasiperiodic properties.

Formation Mechanism in Alloy Steel Rolling Process Using Thermo-mechanical Coupling Method

Based on the theory of elastic-plastic finite element method, the high-speed hot continuous rolling process of a billet is simulated and analyzed in vertical and horizontal passes. The billet is dragged into the passes by contact friction force between the billet and rollers. The rollers and billet are represented by respectively rigid and deformable bodies, and three-dimensional models are developed for the billet and rollers. The distribution of deformation field, effective strain, rolling force and temperature field are accurately calculated for the whole rolling process （including unstable and stable stages）. In addition, the rolling pressure on the width symmetry center is compared with that in the in-situ experimental measurements. It is revealed that various heat exchange phenomena among the billet, rollers and surroundings can result in unbalanced temperature distribution on the cross section. Rolling force and strain can change significantly when the billet is moved towards or away from the roller gap, and keep almost invariable in the stable stage. It is expected that the simulation results would be useful for practical manufacture and provide the theoretical foundation for improvement of process planning and optimization of process parameters.

Coupled Thermo-Mechanical Modeling of Long Cone-Shaped Forging Based on Rigid Visco-Plastic FEM

Based on synthetically considering the coupled thermo mechanical relations between temperature and deforming, a numerical simulation of the forging process for the special long cone shaped workpiece of Al 5.44Mg 2.15Li 0.12Zr alloy at high temperature was conducted by using the rigid visco plastic finite element method. The relations between the total load and the displacement during the forging, and the distributions of stress, strain, temperature and strain rate, which can provide useful information for the process design, are obtained.

基于声学FEM-BEM的户内变电站噪声场计算

准确计算户内变电站大型、复杂的噪声场分布,进而评价可采用降噪措施的减噪效果,是解决户内变电站噪声污染的关键问题。为此,综合声学有限元法(finite element method,FEM)求解复杂声场收敛性好及精度高的优点,及声学边界元法(boundary element method,BEM)降维求解大型声场的优势,提出了一种基于声学FEM-BEM的户内变电站噪声场求解算法。首先,建立变电站内部声源声固耦合模型,采用声学FEM求解混响噪声作用下的声固耦合响应;然后,基于声学FEM-BEM耦合理论,求解内、外耦合边界处结构单元受声固耦合激励产生的位移及应力载荷;最后,根据声压及应力载荷激发的外场声波扩散模型,基于常规Gauss数值积分法,建立外部空间声域2维BEM声学积分方程,求解外部声场。该算法在湖南某110 kV户内变电站噪声场的求解分析中得到了成功应用,与实测值的相对误差为3.61%~4.87%。

基于DEM-FEM耦合方法的碎屑流对桥墩最大冲击力的影响因素研究

碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数值模拟方法模拟了碎屑流对双柱式桥墩的冲击效应,并结合斜槽试验,验证了耦合方法的准确性,进一步分析了碎屑流冲击坡度、距离和体积密度对桥墩冲击力的影响规律。结果表明,最大冲击力与碎屑流冲击坡度、距离和体积密度分别呈幂函数(指数大于1)、幂函数(指数小于1)和线性正相关。冲击坡度、距离和体积密度对最大冲击力的敏感度值分别为3.012、0.202、0.804,在桥梁碎屑流灾害防治时需重视冲击坡度和体积密度的影响。将冲击力的数值模拟值与流体动力学模型预测值对比分析表明,流体动力学模型理论公式能较好地预测桥墩所受的最大冲击力,最大预测误差低于23.6%。相关研究结果可为山区桥梁碎屑流灾害防治与设计提供一定的参考依据。

考虑强度速率衰减效应的地震滑坡SPH-FEM模拟

基于SPH和FEM耦合的数值计算方法,引入滑动面摩擦强度的速率衰减模型,提出了一种能够模拟地震诱发滑坡破坏过程的数值模拟方法。基于所提数值方法模拟了唐家山地震滑坡,模拟结果与现场勘查结果及室内试验现象较为一致。基于模拟的滑动面上摩擦系数的演化过程,将唐家山滑坡的发生分为4个阶段:启动阶段、摩擦衰减阶段、低摩擦滑移阶段和逐步稳定阶段。模拟结果表明速度增加和摩擦强度衰减的相互促进,是触发滑体的高速运动的根本原因。提出采用滑体上作用的动摩擦力fd和动下滑力Td的比值R作为判别指标用于判断大型滑坡的启动,当首次出现R小于1时认为滑动面发生整体贯通并出现失稳启动。基于滑动面不同位置摩擦系数的演化,揭示了滑坡启动中滑动面摩擦强度衰减和滑动面的渐进贯通过程,解释了地震作用与滑动面摩擦参数速率衰减效应共同作用触发大型滑坡发生破坏的内在机理。

Thermo-mechanical coupled particle model for rock

A thermo-mechanical coupled particle model for simulation of thermally-induced rock damage based on the particle simulation method was proposed.The simulation results of three verification examples,for which the analytical solutions are available,demonstrate the correctness and usefulness of the thermo-mechanical coupled particle model.This model is applied to simulating an application example with two cases：one is temperature-independent elastic modulus and strength,while the other is temperature-dependent elastic modulus and strength.The related simulation results demonstrate that microscopic crack initiation and propagation process with consideration of temperature-independent and temperature-dependent elastic modulus and strength are different and therefore,the corresponding macroscopic failure patterns of rock are also different.On the contrary,considering the temperature-dependent elastic modulus and strength has no or little effect on the heating conduction behavior.Numerical results,which are obtained by using the proposed model with temperature-dependent elastic modulus and strength,agree well with the experimental results.This also reveals that the rock subjected to heating experiences much more cracking than the rock subjected to cooling.

Impact of cooling rate on mechanical properties and failure mechanism of sandstone under thermal-mechanical coupling effect

High geo-temperature is one of the inevitable geological disasters in deep engineering such as resource extraction,space development,and energy utilization.One of the key issues is to understand the mechanical properties and failure mechanism of high-temperature rock disturbed by low-temperature airflow after excavation.Therefore,.the experimental and numerical investigation were carried out to study the impact of cooling rate on mechanical properties and failure mechanism of high temperature sandstone.First,uniaxial compression experiments of high temperature sandstone at different real-time cooling rates were carried out to study the mechanical properties and failure modes.The experimental results indicate that the cooling rate has a significant effect on the mechanical properties and failure modes of sandstone.The peak strain,peak stress,and elastic modulus decrease with an increase in cooling rate,and the fragmentation degree after failure increases gradually.Moreover,the equivalent numerical model of heterogeneous sandstone was established using particle flow code(PFC)to reveal the failure mechanism.The results indicate that the sandstone is dominated by intragrain failure in the cooling stage,the number of microcracks is exponentially related to the cooling rate,and the higher the cooling rate,the more cracks are concentrated in the exterior region.Under axial loading,the tensile stress is mostly distributed along the radial direction,and the damage in the cooling stage is mostly due to the fracture of the radial bond.In addition,axial loading,temperature gradient and thermal stress mismatch between adjacent minerals are the main reasons for the damage of sandstone in the cooling stage.Moreover,the excessive temperature gradient in the exterior region of the sandstone is the main reason for the damage concentration in this region.

基于SPH-FEM耦合算法的台阶精准爆破延期时间研究及应用

电子雷管在大型露天矿爆破中广泛使用,利用其精准起爆特性通过设置合理延期时间可有效改善爆破质量,对矿山生产降本增效具有重要意义。基于光滑粒子流体动力学—有限单元(SPH-FEM)耦合算法,利用LS-Dyna软件建立了多个巴润矿常规台阶尺寸的三维数值模型。针对电子雷管起爆条件下不同孔间和排间延时分别建立双孔、三孔SPH-FEM三维模型,从矿岩体SPH粒子应力和移动2个维度研究了被爆矿岩体的受力破碎、抛掷位移特征,并重点分析了不同排间延时下后排炮孔爆破漏斗的偏移角度。据此确定了最佳孔间、排间延期时间范围,并通过现场验证对比试验,选出了适合该爆区的最优延时。研究结果表明:在巴润矿现行孔排距条件下利用电子雷管精准起爆,孔间延时在20~30 ms内可达到较好的破碎效果;排间延时大于35 ms时,先爆炮孔为后续炮孔提供新的自由面和抛掷空间,且大于40 ms时应力趋于稳定;在现场验证试验中,试验区内边长1 m以下的破碎矿岩体占比达91%,较非试验爆区中同尺寸岩块占比提高了14%,破碎效果得到有效改善。

基于SPH-FEM耦合方法的泥石流冲击输电塔基础的动力分析

泥石流是我国西南山区常见的地质灾害。架空输电杆塔在泥石流的冲击下往往发生基础破坏甚至会造成杆塔倒塌。首先采用光滑粒子流体动力学(smoothed particle hydrodynamics,简称SPH)方法和有限元方法(finite element method,简称FEM)相耦合的三维数值方法模拟了泥石流对杆塔基础的冲击作用;在与相关模型试验结果验证的基础上,开展了不同泥石流密度、黏度系数及初始速度条件下对输电塔基础的冲击力作用的参数分析;研究结果表明:随着泥石流初始速度的增加,冲击力峰值会随之增大;前排基础的冲击力峰值均大于后排基础;泥石流冲击过程特性受到泥石流密度和黏度系数影响。与稀性泥石流相比:黏性泥石流冲击基础后,基础下游真空区相对要小;此外,将数值模拟结果与Kwan冲击力公式及铁二院推荐的冲击压力设计公式预测值进行对比分析可以发现:Kwan冲击力公式能较好地预测出基础所受泥石流冲击力的平均趋势,最大预测误差低于30%,铁二院公式预测的稀性和黏性泥石流的冲击压力平均偏低分别约17%和28%。相关研究结果有望为泥石流频发区域输电塔基础的设计和风险评估提供一定的参考依据。

基于自适应FEM-SPH耦合算法的飞机典型部位破片冲击战伤的数值研究

针对飞机典型部位在遭到高速破片攻击后结构整体的战伤状态及破片的剩余行为开展数值模拟。应用LS-DYNA软件,结合有限单元方法(finite element method,FEM)和光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)两者的优势,建立自适应的FEM-SPH耦合模拟方法,并构建2种飞机典型部位的计算模型,采用六面体网格局部细化方法实现了核心位置的精确模拟,并进行试验来验证数值模型;开展了一系列高速冲击战伤模拟,对比了不同工况下破片高速冲击结构后形成的碎片云和破口形貌,并对破片的剩余速度和质量进行分析,确定了破片在结构蒙皮上的临界跳飞角。结果表明:自适应FEM-SPH耦合算法的计算结果与试验结果吻合良好,能够对破片高速冲击战伤进行有效准确模拟;碎片云分布形状随破片速度增加变得狭长,冲击角度会改变碎片云和结构破口形状朝向;碎片云高度和扩散速度随破片速度或角度的变化趋势基本一致并都呈线性关系;破片的速度减少量不随初始速度变化,质量减少量则与冲击速度正相关,两者与冲击角度都负相关;破片临界跳飞角与冲击速度大小基本呈线性关系。

Numerical modeling for the coupled thermo-mechanical processes and spalling phenomena in sp Pillar Stability Experiment (APSE)

In this paper, the coupled thermo-mechanical （TM） processes in the AEspoe Pillar Stability Experiment （APSE） carried out by the Swedish Nuclear Fuel and Waste Management Company （SKB） were simulated using both continuum and discontinuum based numerical methods. Two-dimensional （2D） and three- dimensional （3D） finite element method （FEM） and 2D distinct element method （DEM） with particles were used. The main objective for the large scale in situ experiment is to investigate the yielding strength of crystalline rock and the formation of the excavation disturbed/damaged zone （EDZ） during excavation of two boreholes, pressurizing of one of the boreholes and heating. For the DEM simulations, the heat flow algorithm was newly introduced into the original code. The calculated stress, displacement and temperature distributions were compared with the ones obtained from in situ measurements and FEM simulations. A parametric study for initial microcracks was also performed to reproduce the spalling phenomena observed in the APSE.

基于DEM-FEM耦合模型的质量流量对喷丸强化效果的影响

目的提出一种将有限元(FEM)与离散元(DEM)相耦合的新方法,即通过定义关键字的方式在ABAQUS软件内实现DEM-FEM耦合功能,并应用此方法建立喷丸强化模型。方法将弹丸视为等直径刚体,利用刚体动力学模拟弹丸-弹丸和弹丸-靶材之间相互作用以及整个弹丸流。通过建立的DEM-FEM耦合喷丸强化模型研究质量流量对残余应力场和表面粗糙度的影响。结果质量流量对喷丸强化效果有显著影响,最大残余应力随质量流量的增加而减小,由2 kg/min时的–596.77 MPa下降到6 kg/min时的–581.91 MPa;表面残余应力随质量流量的增加先增大后减小,由2 kg/min时的–420.86 MPa增大到4 kg/min时的472.06 MPa,随后减小到6 kg/min时的–450.50 MPa;表面粗糙度Ra随质量流量的增加减小,由2 kg/min时的11.21μm减小到6 kg/min时的9.82μm。设计了喷丸强化实验,对质量流量2 kg/min的结果进行了验证。结论实测值与模拟值吻合较好,所建立的模型可以准确反映实际喷丸强化过程。质量流量对喷丸强化效果具有显著影响。

Spudcan Penetration Simulation Using the Coupled Eulerian-Lagrangian Method with Thermo-Mechanical Coupled Analysis

A novel modeling technique based on the coupled Eulerian-Lagrangian(CEL) method is provided to solve the geotechnical problems with large deformations. The technique is intended to solve the update problem of soil mechanical properties during spudcan penetration in normally consolidated clay soil. In the CEL model, the normal method of assigning an increasing shear strength profile with depth(NA) is defective due to its Eulerian framework. In this paper, a new technique is proposed to update soil material properties by introducing thermo-mechanical coupled analysis(TMCA) to the CEL models. During establishment of the CEL models, the optimal penetration velocity and minimum mesh size are determined through parametric studies. Reasonability and accuracy are then verified through comparison of the preliminary results with the soil flow configuration and penetration resistance(Fv) of a centrifuge test, and the results of the proposed method are compared with those of the remeshing and interpolation technique with small strain(RITSS) method. To achieve a CEL model with satisfactory accuracy, the NA and TMCA methods implemented in the CEL models and the RITSS method are first adopted in weightless soil. Comparison of the findings with those obtained in previous studies shows that the TMCA method can update material properties and predict Fv. The TMCA method is then applied to soils with self-weight and different shear strength profiles. Results show that the proposed method is capable of accurately modeling the large deformation problem of spudcan penetration in non-homogeneous clay.

Numerical modeling of coupled thermo-mechanical response of a rock pillar

Understanding the rock mass response to excavation and thermal loading and improving the capability of the numerical models for simulating the progressive failure process of brittle rocks are important for safety assessment and optimization design of nuclear waste repositories.The international cooperative DECOVALEX-2011 project provides a platform for development,validation and comparison of numerical models,in which the sp pillar stability experiment（APSE） was selected as the modeling target for Task B.This paper presents the modeling results of Wuhan University（WHU） team for stages 1 and 2 of Task B by using a coupled thermo-mechanical model within the framework of continuum mechanics.The rock mass response to excavation is modeled with linear elastic,elastoplastic and brittle-plastic models,while the response to heating is modeled with a coupled thermo-elastic model.The capabilities and limitations of the model for representation of the thermo-mechanical responses of the rock pillar are discussed by comparing the modeling results with experimental observations.The results may provide a helpful reference for the stability and safety assessment of the hard granite host rock in China＇s Beishan preselected area for high-level radioactive waste disposal.

FEM COUPLING FIELD ITERATION AND ITS CONVERGENCE FOR A GMM ACTURATOR

The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load types is investigated. Then the calculated deformations are compared with the experimental values. The results convince that the CI of FEM is suitable for the simulation of energy coupling and transformation mechanism of the GMM. At last, the output deformation properties are studied under different input currents, showing that there is a good compromise between good linearity and large strain under the prestress 6 MPa.

基于DEM-FEM耦合的环锤式破碎机反击板可靠性分析

由于环锤式破碎机破碎物料的过程是一个极其复杂多变的过程,且与反击板所受的碰撞力是不连续、无规则的,因此,无法使用公式推导或使用单一软件对反击板受力进行计算,为了能够更好地确定破碎机破碎过程中反击板的载荷大小、分布和寿命,文中提出了一种基于离散元(DEM)和有限元(FEM)耦合计算的环锤式破碎机反击板可靠性分析方法,通过EDEM建立环锤式破碎机破碎物料仿真过程,并提取出破碎过程中反击板的载荷分布,与ANSYS Workbench建立耦合计算通道,计算得到在正常工作参数下反击板的受力大小和位置,在此基础上采用n Code Design Life疲劳分析软件对反击板进行寿命分析,得到反击板的受力云图和损伤云图,研究转子转速、反击板倾角、转子与筛板距离等不同因素下对反击板可靠性的影响,结果表明:转子转速为1000 r∕min时反击板的最大受力值为247.9 MPa、反击板角度为60°时应力最大值为206.48 MPa、转子与筛板距离为530 mm时受力最大值为268.37 MPa。