Assessing cutter-rock interaction during TBM tunnelling in granite:Large-scale standing rotary cutting tests and 3D DEM simulations

The widespread utilisation of tunnel boring machines(TBMs)in underground construction engineering requires a detailed investigation of the cutter-rock interaction.In this paper,we conduct a series of largescale standing rotary cutting tests on granite in conjunction with high-fidelity numerical simulations based on a particle-type discrete element method(DEM)to explore the effects of key cutting parameters on the TBM cutter performance and the distribution of cutter-rock contact stresses.The assessment results of cutter performance obtained from the cutting tests and numerical simulations reveal similar dependencies on the key cutting parameters.More specifically,the normal and rolling forces exhibit a positive correlation with penetration but are slightly influenced by the cutting radius.In contrast,the side force decreases as the cutting radius increases.Additionally,the side force shows a positive relationship with the penetration for smaller cutting radii but tends to become negative as the cutting radius increases.The cutter's relative effectiveness in rock breaking is significantly impacted by the penetration but shows little dependency on the cutting radius.Consequently,an optimal penetration is identified,leading to a low boreability index and specific energy.A combined Hertz-Weibull function is developed to fit the cutter-rock contact stress distribution obtained in DEM simulations,whereby an improved CSM(Colorado School of Mines)model is proposed by replacing the original monotonic cutting force distribution with this combined Hertz-Weibull model.The proposed model outperforms the original CSM model as demonstrated by a comparison of the estimated cutting forces with those from the tests/simulations.The findings from this work that advance our understanding of TBM cutter performance have important implications for improving the efficiency and reliability of TBM tunnelling in granite.

Optimal Design of High-Speed Partial Flow Pumps using Orthogonal Tests and Numerical Simulations

To investigate the influence of structural parameters on the performances and internal flow characteristics of partial flow pumps at a low specific speed of 10000 rpm,special attention was paid to the first and second stage impeller guide vanes.Moreover,the impeller blade outlet width,impeller inlet diameter,blade inclination angle,and number of blades were considered for orthogonal tests.Accordingly,nine groups of design solutions were formed,and then used as a basis for the execution of numerical simulations(CFD)aimed at obtaining the efficiency values and heads for each design solution group.The influence of impeller geometric parameters on the efficiency and head was explored,and the“weight”of each factor was obtained via a range analysis.Optimal structural parameters were finally chosen on the basis of the numerical simulation results,and the performances of the optimized model were verified accordingly(yet by means of CFD).Evidence is provided that the increase in the efficiency and head of the optimized model was 12.11%and 23.5 m,respectively,compared with those of the original model.

Molecular simulation study of the microstructures and properties of pyridinium ionic liquid[HPy][BF_(4)]mixed with acetonitrile

The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.

Centrifugal model tests and numerical simulations for barrier dam break due to overtopping

The present study focuses on the breaching process and failure of barrier dams due to overtopping. In this work, a series of centrifugal model tests is presented to examine the failure mechanisms of landslide dams. Based on the experimental results, failure process and mechanism of barrier dam due to overtopping are analyzed and further verified by simulating the experimental overtopping failure process. The results indicate that the barrier dam will develop during the entire process of overtopping in the width direction, whereas the breach will cease to develop at an early stage in the depth direction because of the large particles that accumulate on the downstream slope. Moreover, headcut erosion can be clearly observed in the first two stages of overtopping, and coarsening on the downstream slope occurs in the last stage of overtopping. Thus, the bottom part of the barrier dam can survive after dam breaching and full dam failure becomes relatively rare for a barrier dam. Furthermore, the remaining breach would be smaller than that of a homogeneous cohesive dam under the same conditions.

A study on dynamic response of slopes under wave action using simulation tests

After the erection of the Three Gorges Dam, the water level of Yangtze River will reach 175 m, and the average wave crest will be up to 1 m. Therefore the wave action cannot be neglected for the slope stability. Through simulation tests, the wave-induced dynamic response of the slope is analyzed. The soil body is taken as linear elastic body when it has a small deformation under the small wave action. Based on tests, the excess pore pressure and slope displacement under the loading in different wave period are analyzed. The ratio of dynamic strength and static strength to the breaking process of the slope is discussed. It is demonstrated that smaller wave period gives rise to a larger strain of the slope under the same stress. At different depth of water, different weakness effect on the stability of the soil slope is observed and the slope has an adaptability to the wave action to some extent.

FATIGUE LIFE PREDICTION OF CRANKSHAFT MADE OF MATERIAL 48MnV BASED ON FATIGUE TESTS,DYNAMIC SIMULATION AND FEA

S-N curve and fatigue parameters of 48MnV are obtained using small sample tests and staircase or up and down method, which paves the way for predicting fatigue life of crankshaft made of 48MnV. The fatigue life of the crankshaft of a six-cylinder engine is calculated using different damage models such as S-N method, normal strain approach, Smoth-Watson-Topper （SWT）Bannantine approach, shear strain approach, and Fatemi-Socie method based on dynamic simulation and finite element analysis （FEA） of crankshaft. The results indicate that the traditional calculation is conservative and the residual fatigue life of crankshaft is sufficient to maintain next life cycle if the crankshaft is remanufactured after its end of life.

Numerical Simulations and Model Tests of the Mooring Characteristic of A Tension Leg Platform Under Random Waves

Analyzing the dynamic response and calculating the tendon tension of the mooring system are necessary for the structural design of a tension leg platform （TLP）. The six-degree-of-freedom dynamic coupling responses and the mooring characteristics of TLP under random waves are studied by using a self-developed program. Results are verified by the 1：40 scaling factor model test conducted in the State Key Laboratory of Ocean Engineering at Shanghai JiaoTong University. The mean, range, and standard deviation of the numerical simulation and model test are compared. The influences of different sea states and wave approach angles on the dynamic response and tendon tension of the mooring system are investigated. The acceleration in the center and corner of the deck is forecasted.

Trigger mechanism of loess-mudstone landslides inferred from ring shear tests and numerical simulation

Whereas loess-mudstone landslides are widely distributed and frequently occurred at the loess Plateau,this type of landslides is hard to detect due to its particularity,and easily generates serious losses.To clarify the shear characteristics and formation mechanism of loess-mudstone landslides,field investigations,ring shear tests and numerical simulation analyses were performed on the loess specimens collected from the Dingjiagou landslide in Yan’an city,China.The test results showed that both the peak strength and residual strength of slip zone soils have a decreasing tendency with moisture content,while the increasing of normal stress caused an increase in the shear strength.These phenomena indicate that the rise in the moisture content induced by precipitation or the decreasing of normal stress due to excavation activities would result in the weakening of slip zone soils.Numerical simulations of the evolution process of slope failure using the finite element method were conducted based on the Mohr–Coulomb criterion.It was found that the heavy precipitation played a more important role in the slope instability compared with the excavation.In addition,the field investigation showed that loess soils with well-developed cracks and underlying mudstone soils provide material base for the formation of loess-mudstone landslides.Finally,the formation mechanism of this type of landslides was divided into three stages,namely,the local deformation stage,the penetration stage,the creeping-sliding stage.This study may provide a basis for understanding the sliding process of loess-mudstone landslides,as well as guidelines for the prevention and mitigation of loess-mudstone landslides.

Constitutive Modeling for Ti-6Al-4V Alloy Machining Based on the SHPB Tests and Simulation

A constitutive model is critical for the prediction accuracy of a metal cutting simulation. The highest strain rate involved in the cutting process can be in the range of 104-106 s 1. Flow stresses at high strain rates are close to that of cutting are difficult to test via experiments. Split Hopkinson compression bar （SHPB） technology is used to study the deformation behavior of Ti-6Al-4V alloy at strain rates of 10 -4-10 4s- 1. The Johnson Cook （JC） model was applied to characterize the flow stresses of the SHPB tests at various conditions. The parameters of the JC model are optimized by using a genetic algorithm technology. The JC plastic model and the energy density-based ductile failure criteria are adopted in the proposed SHPB finite element simulation model. The simulated flow stresses and the failure characteristics, such as the cracks along the adiabatic shear bands agree well with the experimental results. Afterwards, the SHPB simulation is used to simulate higher strain rate（approximately 3 × 10 4 s -1） conditions by minimizing the size of the specimen. The JC model parameters covering higher strain rate conditions which are close to the deformation condition in cutting were calculated based on the flow stresses obtained by using the SHPB tests （10 -4 - 10 4 s- 1） and simulation （up to 3 × 10 4 s - 1）. The cutting simulation using the constitutive parameters is validated by the measured forces and chip morphology. The constitutive model and parameters for high strain rate conditions that are identical to those of cutting were obtained based on the SHPB tests and simulation.

Dynamic simulation of pilot thickener operation using phenomenological model with results validation for continuous and discontinuous tests

The phenomenological theory of sedimentation-thickening processes predicts the settling behavior of a flocculated suspension in dependence of two functions, the batch flux density function and the effective solid stress. These functions were determined using batch settling tests. The governing equations for sedimentation were then solved numerically for these functions and the predictions were compared to the experimental results from pilot scale thickener tests. Firstly, the continuous tests were performed in the plexiglass pilot thickener at different feed flow rates and discharge rates and the solid volume fraction of discharge, the bed height and the time were recorded for each condition. These tests were also simulated and it was observed that there is a good agreement between the results of continuous tests and the results of dynamic simulation. Secondly, the discontinuous tests were performed in the plexiglass pilot thickener at different feed flow rates with a discharge rate of zero. The bed formation rate was determined for each condition. These tests were also simulated and it was observed that there is a good agreement between the results of discontinuous tests and the results of simulation.

DISCREPANCIES IN THE REGRESSION MODELLING OF RECRYSTALLIZATION RATE AS USING THE DATA FROM PHYSICAL SIMULATION TESTS

The analysis of numerous experimental equations published in the literature reveals awide scatter in the predictions for the static recrystallization kinetics of steels. Thepowers of the deformation variables, strain and strain rate, similarly as the powerof the grain size vary in these equations. These differences are highlighted and thetypical values are compared between torsion and compression tests. Potential errorsin physical simulation testing are discussed.

Simulation of unconventional well tests with the finite volume method

The finite volume method has been successfully applied in several engineering fields and has shown outstanding performance in fluid dynamics simulation. In this paper, the general framework for the simulation ofnear-wellbore systems using the finite volume method is described. The mathematical model and the numerical model developed by the authors are presented and discussed. A radial geometry in the vertical plane was implemented so as to thoroughly describe near-wellbore phenomena. The model was then used to simulate injection tests in an oil reservoir through a horizontal well and proved very powerful to correctly reproduce the transient pressure behavior. The reason for this is the robustness of the method, which is independent of the gridding options because the discretization is performed in the physical space. The model is able to describe the phenomena taking place in the reservoir even in complex situations, i.e. in the presence of heterogeneities and permeability barriers, demonstrating the flexibility of the finite volume method when simulating non-conventional tests. The results are presented in comparison with those obtained with the finite difference numerical approach and with analytical methods, if possible.

Numerical Simulation of PMM Tests for A Ship in Close Proximity to Sidewall and Maneuvering Stability Analysis

As the maneuverability of a ship navigating close to a bank is influenced by the sidewall, the assessment of ship maneuvering stability is important. The hydrodynamic derivatives measured by the planar motion mechanism （PMM） test provide a way to predict the change of ship maneuverability. This paper presents a numerical simulation of PMM model tests with variant distances to a vertical bank by using unsteady RANS equations. A hybrid dynamic mesh technique is developed to realize the mesh configuration and remeshing of dynamic PMM tests when the ship is close to the bank. The proposed method is validated by comparing numerical results with results of PMM tests in a circulating water channel. The first-order hydrodynamic derivatives of the ship are analyzed from the time history of lateral force and yaw moment according to the multiple-run simulating procedure and the variations of hydrodynamic derivatives with the ship-sidewall distance are given. The straight line stability and directional stability are also discussed and stable or unstable zone of proportional-derivative （PD） controller parameters for directional stability is shown, which can be a reference for course keeping operation when sailing near a bank.

Laboratory Model Tests and DEM Simulations of Unloading-Induced Tunnel Failure Mechanism

Tunnel excavation is a complicated loading-unloading-reloading process characterized by decreased radial stresses and increased axial stresses.An approach that considers only loading,is generally used in tunnel model testing.However,this approach is incapable of characterizing the unloading effects induced by excavation on surrounding rocks and hence presents radial and tangential stress paths during the failure process that are different from the actual stress state of tunnels.This paper carried out a comparative analysis using laboratory model testing and particle flow code(PFC2D)-based numerical simulation,and shed light upon the crack propagation process and,microscopic stress and force chain variations during the loading-unloading process.The failure mode observed in the unloading model test is shear failure.The force chains are strongly correlated with the concrete fracture propagation.In addition,the change patterns of the radial and tangential stresses of surrounding rocks in the broken region,as well as the influence of the initial stress on failure loads are revealed.The surrounding soil of tunnel failure evolution as well as extent and shape of the damage zone during the excavation-induced unloading were also studied.

Computer Simulation of the Indentation Creep Tests on Particle-Reinforced Composites

A systematical simulation has been carried out on the indentation creep test on particle-reinforced composites. The deformation, failure mechanisms and life are analyzed by three reasonable models. The following five factors have been considered simultaneously: creep property of the particle, creep property of the matrix, the shape of the particle, the volume fraction of the particle and the size (relative size to the particle) of the indentation indenter. For all the cases, the power law respecting to the applied stress can be used to model the steady indentation creep depth rate of the indenter, and the detail expressions have been presented. The computer simulation precision is analyzed by the two-phase model and the three-phase model. Two places of the stress concentration are found in the composites. One is ahead of the indentation indenter, where the high stress state is deduced by the edge of the indenter and will decrease rapidly near to a steady value with the creep time. The other one is at the interface, where the high stress state is deduced by the misfit of material properties between the particles and matrix. It has been found that the creep dissipation energy density other than a stress parameter can be used to be the criterion to model the debonding of the interfaces. With the criterion of the critical creep dissipation energy density, a power law to the applied stress with negative exponent can be used to model the failure life deduced by the debonding of interfaces. The influences of the shape of the particles and the matching of creep properties of particle and matrix can be discussed for the failure. With a crack model, the further growth of interface crack is analyzed, and some important experimental phenomena can be predicted. The failure mechanism which the particle will be punched into matrix has been also discussed. The critical differences between the creep properties of the particles and matrix have been calculated, after a parameter has been defined. In the view of competition of failure mechanisms, the best matching of the creep properties of the two phases and the best shape of the particles are discussed for the composite design.

Numerical simulation of direct shear tests on mechanical properties of talus deposits based on self-adaptive PCNN digital image processing

The macro mechanical properties of materials with characteristics of large scale and complicated structural composition can be analyzed through its reconstructed meso-structures.In this work,the meso-structures of talus deposits that widely exist in the hydro-power engineering in the southwest of China were first reconstructed by small particles according to the in-situ photographs based on the self-adaptive PCNN digital image processing,and then numerical direct shear tests were carried out for studying the mechanical properties of talus deposits.Results indicate that the reconstructed meso-structures of talus deposits are more consistent with the actual situation because the self-adaptive PCNN digital image processing has a higher discrimination in the details of soil-rock segmentation.The existence and random distribution of rock blocks make the initial shear stiffness,the peak strength and the residual strength higher than those of the "pure soil" with particle size less than 1.25 cm apparently,but reduce the displacements required for the talus deposits reaching its peak shear strength.The increase of rock proportion causes a significant improvement in the internal friction angle of talus deposit,which to a certain degree leads to the characteristics of shear stress-displacement curves having a changing trend from the plastic strain softening deformation to the nonlinear strain hardening deformation,while an unconspicuous increase in cohesion.The uncertainty and heterogeneity of rock distributions cause the differences of rock proportion within shear zone,leading to a relatively strong fluctuation in peak strengths during the shear process,while movement features of rock blocks,such as translation,rotation and crossing,expand the scope of shear zone,increase the required shear force,and also directly lead to the misjudgment that the lower shear strength is obtained from the samples with high rock proportion.That,however,just explains the reason why the shear strength gained from a small amount of indoor test data is not consistent with engineering practice.

DEM Simulation to Determine the Influence on the Experimental Results of Tests of Iron Pellets When the Dimensions of the Test Device Are Varied

The current study is based on the DEM computer simulation of three experimental test devices with different dimensions to determine the difference in the results of the formation of shear and repose angles that the particles experience when grouped under the action of the gravitational force. In this respect, the experimental test devices with different height, width, and depth were geometrically modeled with iron pellet particles using morphology and a granulometric variation from 6 mm to 9 mm of equivalent diameter in its spherical shape. Depending on the results obtained, a reliable size of the experimental test device will be available to obtain the necessary data for a correct adjustment of the calibration parameters for the DEM simulation of mining-metallurgical processes that use granulated material of iron pellet.

基于Plant Simulation的化纤自动落丝系统仿真实验分析

基于Plant Simulation软件,构建一种化纤自动落丝系统的3D数字化模型,对系统效率关键输入因子进行参数化设置,将系统各工位三维数字模型导入软件中建立层次化运动机构图形结构,并运用Simltalk语言实现三维动作仿真,通过基础物理参数设置,保证了仿真数字化模型与现实系统更具一致性。通过多级实验设计分析了系统影响因子对于系统效率的影响特性曲线,并进一步通过动态参数化实验方法,计算出双输入因子对系统效能影响的敏感度。该实验结果可为化纤自动落丝系统的建设成本控制与可行性分析提供指导,具有较好的工程应用价值。

基于Simulation X的多轴转向系统功能安全仿真分析

基于ISO 26262《道路车辆功能安全》的标准要求,对特种车的多轴电液转向系统进行分析,以提高系统的安全性和可靠性。运用Simulation X软件建立详细的多轴特种车仿真模型,通过模拟故障模式注入进行仿真试验。对仿真结果和数据进行分析,评估得到故障的严重性、暴露度和可控性,从而确定相应的汽车安全完整性等级。基于故障注入仿真的汽车功能安全分析方法,在系统早期设计阶段可以作为评估架构安全性的重要手段。

Centrifuge Model Tests and Numerical Simulations of the Impact of Underwater Explosion on an Air-Backed Steel Plate

Damage and threats to hydraulic and submarine structures by underwater explosions(UNDEXs)have raised much attention.The centrifuge model test,compared to prototype test,is a more promising way to examine the problem while reducing cost and satisfying the similitude requirements of both Mach and Froude numbers simultaneously.This study used a systematic approach employing centrifuge model tests and numerical simulations to investigate the effects of UNDEXs on an air-backed steel plate.Nineteen methodical centrifuge tests of UNDEXs were conducted.The shock wave pressure,bubble oscillation pressure,acceleration and the strain of the air-backed steel plate were recorded and compared with numerical studies using the finite element analysis(FEA)commercial software ABAQUS.By implementing empirically derived and physically measured pressures into the numerical models,the effects of the shock wave and bubble oscillation on the steel plate were investigated.Generally,the numerical results were in agreement with the experimental results.These results showed that the peak pressure of an UNDEX has a significant effect on the peak acceleration of the steel plate and that the impulse of the UNDEX pressure governs the peak strain of the steel plate.