司马煤矿3号煤层底板采动破坏深度数值模拟研究

根据井下底板、含水层条件和生产安全需求,文章对司马煤矿3号煤层开采过程中底板采动破坏深度进行了数值模拟研究,研究结果表明:开采引起的底板裂隙最大破坏深度为17.5 m,但3号煤层距离奥陶系中统石灰岩岩溶裂隙极强富水含水层较远,所以绝大区域3号煤层的开采是安全的。

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.

基于PSoC的矿用传感器模拟电路设计

本文设计了一款基于PSoC的矿用传感器模拟电路,由PSoC最小系统板、数模转换模块、电流环模块三部分组成。该电路可以输出高达四路0-3.3V或0-5V电压、八路4-20m A电流环信号,以及62路开关量信号。该矿用传感器模拟电路可以模拟常用的矿用传感器的输出模拟信号,数字信号。操作方便,可通过简单的配置生成所需要的多种传感器输出的预处理信号,避免了购买高昂的传感器。

Operation optimization of prefabricated light modular radiant heating system:Thermal resistance analysis and numerical study

The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.

Case study on the mechanics of NPR anchor cable compensation for large deformation tunnel in soft rock in the Transverse Mountain area,China

A study was conducted to analyze the deformation mechanism of strongly weathered quartz schist in the Daliangshan Tunnel,located in the western Transverse Mountain area.A large deformation problem was experienced during the tunnel construction.To mitigate this problem,a support system was designed incorporating negative Poisson ratio(NPR)anchor cables with negative Poisson ratio effect.Physical model experiments,field experiments,and numerical simulation experiments were conducted to investigate the compensation mechanical behavior of NPR anchor cables.The large deformations of soft rocks in the Daliangshan Tunnel are caused by a high ground stress,a high degree of joint fracture development,and a high degree of surrounding rock fragmentation.A compensation mechanics support system combining long and short NPR anchor cables was suggested to provide sufficient counter-support force(approximately 350 kN)for the surrounding rock inside the tunnel.Comparing the NPR anchor cable support system with the original support system used in the Daliangshan tunnel showed that an NPR anchor cable support system,combining cables of 6.3 m and 10.3 m in length,effectively prevented convergence of surrounding rock deformation,and the integrated settlement convergence value remained below 300 mm.This study provides an effective scientific basis for resolving large deformation problems in deeply buried soft rocks in western transverse mountain areas.

Resistivity tomography based on multichannel electrodes

With the constantly changing engineering construction sector,the detection accuracy of conventional electrical resistivity tomography(ERT)is no longer suffi cient.A multichannel electrode design(MERT)-based ERT is introduced in this paper to address the growing need for resolution.The imaging accuracy of the ERT method is improved through the collection of apparent resistivity data in various directions by measuring the potential diff erence between diff erent channels.Numerical simulation results of the inclined high-resistivity anomaly model reveal that MERT is a precise representation of the shape,inclined direction,and buried depth of the anomaly,with thoroughfare M2N2 producing the most precise forward and inverse results.Based on the analysis results of the model resolution matrix,when the buried depth of power supply points and the gap between potential acquisition points are 30%-90%and 30%-60%of the electrode distance,respectively,the MERT approach yields superior detection outcomes.The detection eff ect of the MERT method on anomalous bodies with diff erent burial depths under the optimal parameters also indicates that the MERT method can obtain richer potential change information with higher resolution in deep areas compared to the ERT method.With the implementation of the MERT approach,the scope of applications for ERT is expanded,the accuracy of ERT detection is increased,and the progress of near-surface fi ne detection is positively infl uenced.

A Hybrid Dung Beetle Optimization Algorithm with Simulated Annealing for the Numerical Modeling of Asymmetric Wave Equations

In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.

Optimization of Channel Structure of Alkaline Water Electrolyzer by Using an Expanded Mesh as a Bipolar Plate

Alkaline water electrolysis(AWE)is the most mature technology for hydrogen production by water electrolysis.Alkaline water electrolyzer consists of multiple electrolysis cells,and a single cell consists of a diaphragm,electrodes,bipolar plates and end plates,etc.The existing industrial bipolar plate channel is concave-convex structure,which is manufactured by complicated and high-cost mold punching.This structure still results in uneven electrolyte flow and low current density in the electrolytic cell,further increasing in energy consumption and cost of AWE.Thereby,in this article,the electrochemical and flow model is firstly constructed,based on the existing industrial concave and convex flow channel structure of bipolar plate,to study the current density,electrolyte flow and bubble distribution in the electrolysis cell.The reliability of the model was verified by comparison with experimental data in literature.Among which,the electrochemical current density affects the bubble yield,on the other hand,the generated bubbles cover the electrode surface,affecting the active specific surface area and ohmic resistance,which in turn affects the electrochemical reaction.The result indicates that the flow velocity near the bottom of the concave ball approaches zero,while the flow velocity on the convex ball surface is significantly higher.Additionally,vortices are observed within the flow channel structure,leading to an uneven distribution of electrolyte.Next,modelling is used to optimize the bipolar plate structure of AWE by simulating the electrochemistry and fluid flow performances of four kinds of structures,namely,concave and convex,rhombus,wedge and expanded mesh,in the bipolar plate of alkaline water electrolyzer.The results show that the expanded mesh channel structure has the largest current density of 3330 A/m^(2)and electrolyte flow velocity of 0.507 m/s in the electrolytic cell.Under the same current density,the electrolytic cell with the expanded mesh runner structure has the smallest potential and energy consumption.This work provides a useful guide for the comprehensive understanding and optimization of channel structures,and a theoretical basis for the design of large-scale electrolyzer.

3D modeling for effect of tool eccentricity on coupled thermal and material flow characteristics during friction stir welding

A novel three-dimensional numerical model is proposed to investigate the effect of tool eccentricity on the coupled thermal and material flow characteristics in friction stir welding(FSW) process.An asymmetrical boundary condition at the tool-workpiece interface,and the dynamic mesh technique are both employed for the consideration of the tool eccentricity during tool rotating.It is found that tool eccentricity induces the periodical variation of the heat densities both at the tool-workpiece interface and inside the shear layer,but the fluctuation amplitudes of the heat density variations are limited.However,it is demonstrated that tool eccentricity results in significant variation of the material flow behavior in one tool rotating period.Moreover,the material velocity variation at the retreating side is particularly important for the formation of the periodic characteristics in FSW.The modeling result is found to be in good agreement with the experimental one.

Numerical simulation on magnetohydrodynamic power generation channel of scramjet

The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distance flight airframe power.In this paper,a computational model of the scramjet magnetohydrody⁃namic channel is developed and verified by using the commercial software Fluent.It is found that when the mag⁃netic induction intensity is 1,2,3,4 T,the power generation efficiency is 22.5%,22.3%,22.0%,21.5%,and decreases with the increase of the magnetic induction intensity,and the enthalpy extraction rate is 0.026%,0.1%,0.21%,0.34%,and increases with the increase of the magnetic induction intensity.The deceleration ef⁃fect of electromagnetic action on the airflow in the power channel increases with the increase of magnetic induc⁃tion intensity.The stronger the magnetic field intensity,the more obvious the decreasing effect of fluid Mach num⁃ber in the channel.The power generation efficiency decreases as the magnetic induction intensity increases and the enthalpy extraction rate is reversed.As the local currents gathering at inlet and outlet of the power generation area,total temperature and enthalpy along the flow direction do not vary linearly,and there are maximum and minimum values at inlet and outlet.Increasing the number of electrodes can effectively regulate the percentage of Joule heat dissipation,which can improve the power generation efficiency.

Numerical Simulation of Flow Field and Flow State Division in Thin-Film Evaporators

The flow field and flow state of thin-film evaporators are complex,and it is significant to effectively divide and quantify the flow field and flow state,as well as to study the internal flow field distribution and material mixing characteristics to improve the efficiency of thin-film evaporators.By using computational fluid dynamics(CFD)numerical simulation,the distribution pattern of the high-viscosity fluid flow field in the thin-film evaporators was obtained.It was found that the staggered interrupted blades could greatly promote material mixing and transportation,and impact the film formation of high-viscosity materials on the evaporator wall.Furthermore,a flow field state recognition method based on radial volume fraction statistics was proposed,and could quantitatively describe the internal flow field of thin-film evaporators.The method divides the high-viscosity materials in the thin-film evaporators into three flow states,the liquid film state,the exchange state and the liquid mass state.The three states of materials could be quantitatively described.The results show that the materials in the exchange state can connect the liquid film and the liquid mass,complete the material mixing and exchange,renew the liquid film,and maintain continuous and efficient liquid film evaporation.

Three-dimensional CSAMT Forward Modeling of Potential Landslide Sliding Surfaces Using Finite Element Method

Landslides are a type of natural disaster that can cause substantial harm to humanity.Monitoring and predicting the initiation of potential landslides is critical to avoiding losses due to disasters and economic activities.The impact of the controlled-source audio-frequency magnetotelluric method on investigating landslide surfaces is assessed through numerical simulations with a finite element approach.A Dirichlet boundary condition is selected to match the truncated boundary,resulting in a remarkable improvement in simulation effi ciency.Rederivation of the formulas for a layered medium adept to the controlled-source audiofrequency magnetotelluric method is necessary to determine the electromagnetic fi eld at any location along the truncated boundary.After the reliability evaluation of the new codes,a landslide model with a slide surface is designed,and the characteristics of its electromagnetic fi eld and the apparent resistivity are studied.Instead of the total electromagnetic fi eld,which is strongly infl uenced by topography variation,the apparent resistivity should be used for sliding surface detection.The normalized pure anomalous electromagnetic fi eld may also be employed to quickly assess the detectability of the sliding surface.Overall,this study demonstrates that the controlled-source audio-frequency magnetotelluric method can be employed for investigating landslides,and recommends survey parameters,including confi guration,frequency range,and length of survey line in landslide exploration.

Forming defects in aluminum alloy hot stamping of side-door impact beam

The forming defects, including thinning, rupture, wrinkling and springback, usually arising in producing a side-door impact beam, were investigated by trial and numerical simulation. A temperature-related constitutive model specific to the temperature range from 350 °C to 500 °C was established and used for the numerical simulation. The trial and numerical simulation were conducted to clarify the quantitative characteristics of forming defects and to analyze the effects of process parameters on the forming defects. Results show that the rupture situation is ameliorated and the springback is eliminated in the aluminum alloy hot stamping. The wrinkling severity decreases with increasing blank holder force (BHF), but the BHF greater than 15 kN causes the rupture at the deepest drawing position of workpiece. The forming defects are avoided with lubricant in the feasible ranges of process parameters: the BHF of 3 to 5 kN and the stamping speed of 50 to 200 mm/s.

Forced convection rheoforming process for preparation of 7075 aluminum alloy semisolid slurry and its numerical simulation

A self-developed forced convection rheoforming （FCR） machine for the preparation of light alloy semisolid slurry was introduced. The microstructure characteristics of 7075 aluminium alloy semisolid slurry at different stirring speeds prepared by the FCR process were analyzed. The experimental results suggest that with the increase of the stirring speed, the mean grain size of the semisolid decreases and the shape factor as well as the number of primary grains increase. Meanwhile, the preparation process of semisolid slurry was numerically simulated. The flow characteristics of the melt in the device and the effect of the stirring speed on temperature field and solid fraction of the melt were investigated. The simulated results show that during the preparation process of semisolid slurry, there is a complex convection within the FCR device that obviously changes the temperature field distribution and solid fraction of the melt. When the convection intensity increases, the scope of the undercooling gradient of the melt is reduced and temperature distribution is improved.

Numerical simulation of intermediate phase growth in Ti/Al alternate foils

To investigate the diffusion reaction between Ti/Al solid diffusion couple, Ti/Al alternate foils formed by hot pressing were annealed at 525, 550, 575 and 600 °C for time ranging from 1 to 40 h. The experimental results show that TiAl3 was the only observed phase at Ti/Al interface. The interface thermodynamics favored the preferential formation of TiAl3 in Ti/Al couple. The growth of TiAl3 layer occurred mainly towards Al foil side and exhibited a parabolic law. Using the interdiffusion coefficients calculated based on the contribution of grain boundary diffusion, the growth of TiAl3 was simulated numerically with the finite difference method, and the simulated results were in good agreement with the experimental ones.

NUMERICAL SIMULATIONS OF TEMPERATURE FIELD IN DIRECT METAL LASER SINTERING PROCESS

A mathematical model is developed for simulating the heat transferring behavior in a direct metal laser sintering process. The model considers the thermal phenomena involved in the process, including conduction, radiation, and convection. A formula for the calculation of the heat conductivity of a sintering system containing solid phase, liquid phase, and gas phase is given. Due to the continuous movement of the laser beam, a local coordinate system centered on the laser beam is used to simplify the analytical calculation. Assuming that it is approximately a Gaussian laser beam, the heat conduction equation is resolved based on the assumption of the thermal insulating boundary conditions and the fixed thermal physical parameters. The FORTRAN language is employed to compile the program to simulate the temperature field in the direct copper powder sintering process. It shows a good agreement with the preliminary experimental results.[KH3/4D]

Visualization and simulation of plastic material flow in friction stir welding of 2024 aluminium alloy plates

Thin copper sheets as marker material were embedded into weld path of 2024 aluminium alloy plates and their final position after friction stir welding was examined by metallographic techniques. Referring to the visualized material flow patterns, a three-dimensional model was developed to conduct the numerical simulation of the temperature profile and plastic material flow in friction stir welding. The calculated velocity contour of plastic flow in close proximity of the tool is generally consistent with the visualized results. As the tool rotation speed increases at a constant tool travel speed, the material flow near the pin gets stronger. The predicted shape and size of the weld nugget zone match with the experimentally measured ones.

CFD numerical simulation of flow velocity characteristics of hydrocyclone

A CFD based numerical simulation of flow velocity of hydrocyclone was conducted with different structural and operational parameters to investigate its distribution characteristics and influencing mechanism. The results show there exist several unsymmetrical envelopes of equal vertical velocities in both upward inner flows and downward outer flows in the hydrocyclone, and the cone angle and apex diameter have remarkable influence on the vertical location of the cone bottom of the envelope of zero vertical velocity. It is also found that the tangential velocity isolines exist in the horizontal planes located in the effective separation region of hydrocyclone. The increase of feed pressure has almost no effect on the distribution characteristics of both vertical velocity and tangential velocity in hydrocyclone, but the magnitude and gradient of tangential velocity are increased obviously to make the motion velocity of high density particles to the wall increased and to make the cyclonic separation effect improved.

Numerical simulation of acoustic pressure field for ultrasonic grain refinement of AZ80 magnesium alloy

Ultrasound with different intensities was applied to treating AZ80 alloy melt to improve its solidification structure.The average grain size of the alloy could be decreased from 303 to 148 μm after the ultrasound with intensity of 30.48 W/cm2 was applied.To gain insight into the mechanism of ultrasonic treatment which affected the microstructure of the alloy,numerical simulations were carried out and the effects of different ultrasonic pressures on the behaviors of cavitation bubble in the melt were studied.The ultrasonic field propagation in the melt was also characterized.The results show that samples from different positions are subjected to different acoustic pressures and the effect of grain refinement by ultrasonic treatment for these samples is different.With the increase of ultrasonic intensity,the acoustic pressure is increased and the grain size is decreased generally.

Inhomogeneity of density and mechanical properties of A357 aluminum alloy backward extruded in semi-solid state

The inhomogeneity of density and mechanical properties of A357 aluminum alloy in the semi-solid state were investigated.Numerical simulation and backward extrusion were adopted to study the preparation of cup shells.The results show that the relative density of the wall is the lowest in samples,and that of the base is the highest.With increasing the billet height,more time is needed for relative density of the corner to reach the maximum value,and the relative densities in every region improve evidently with increasing the pressure.The tensile stress was simulated to be the largest at the corner,and the hot tearings were forecasted to mainly appear at the corner too.By employing proper billet height and pressure,the extruded samples consisted of fine and uniform microstructures,and can obtain excellent mechanical properties and Brinell hardness.