Search

To investigate the effect of grinding speed on the surface/subsurface damage of monocrystalline silicon grinding, experiments of grinding monocrystalline silicon wafers with different grinding speeds under the condition of equal abrasive cutting depth were designed based on the workpiece rotation method. The surface morphology, roughness, surface cracks, deformations, and residual stresses after machining were analyzed through comparative experiments, and the reasons for the effect of different machining speeds on the number of cracks and depth of subsurface damage were analyzed in combination with molecular dynamics simulation models. The results indicated that with the increase in grinding speed, the surface quality gradually becomes better. Furthermore, when the grinding speed is 2399 rpm, the surface roughness Ra value formed by Nos. 600 and 5000 grinding wheels is reduced by 22.5% and 31.2%, respectively, compared with 239 rpm. Besides, when grinding with the No. 600 grinding wheel, the higher the grinding speed, the fewer the number of subsurface cracks and the smaller the crack depth. Besides, when the grinding speed is 2399 rpm, the number of cracks is reduced by about 50%, and the depth of cracks is reduced by about 25% compared with 239 rpm. Moreover, when grinding with the No. 5000 grinding wheel, the deformation and residual stress of the silicon wafer gradually decrease with the increase in grinding speed. The deformation formed at a grinding speed of 2399 rpm is about 44% lower than that at 239 rpm. With the increase in grinding speed, the number of cracks, deformation, and residual stress of silicon wafer are reduced. Moreover, the increase in grinding speed will help reduce the damage degree of the crystal structure and the thickness of the damage layer of the silicon wafer, thus improving the quality of the grinding surface.

The curvature variation regularity of the blade surface is complex and requires a high level of machining accuracy. In order to research the influence of interpolation tolerance on the machining accuracy of blade surface at different positions, the non-uniform rational b-spline was used to build the blade surface and analyze the curvature variation regularity at different positions, and the distribution regularity of interpolation tolerance at different curvature positions on interpolation moving points was explored based on the NC (Numerical Control) program. In addition, the NC interpolation principle is used to analyze the effect of different interpolation tolerance parameters on the machining accuracy of blade surfaces, and the appropriate interpolation moving point coordinates were obtained through the five-axis coordinate system transformation theory, and the five-axis cutting simulation and experiment were carried out. The results show that there are many interpolation moving points at the position with large curvature. With interpolation tolerances of 0.1 mm, 0.03 mm and 0.003 mm, the interpolation errors obtained in the regions with large curvature are 0.005 mm, 0.0034 mm and 0.0025 mm, respectively. The smaller the interpolation tolerance, the smaller the interpolation error. What is more, the smaller the interpolation tolerance is, the more obvious the improvement of interpolation accuracy will be when the interpolation tolerance is reduced. Besides, when sampling length of 250 μm, the surface roughness (Ra) detected at the small curvature position and the large curvature position are 0.434 and 1.070μm respectively. The surface roughness is better in position with smaller curvature, however, it is worse in position with larger curvature. The analytical regularity is in good agreement with both simulation and experiment values. It provides a reference for improving NC machining accuracy for similar blade surfaces.

The fragility curve expresses the probability that an asset exceeds some serviceability state for a given level of environmental perturbation or other loadings. It is an important component in the quantitative risk analysis and resilience evaluation of infrastructure exposed to natural hazards. Incidences of over-settlement of embankments are increasingly reported due to more intense and longer-duration rainfall events. This paper develops fragility curves for the rainfall-induced embankment settlement. For this purpose, an embankment incorporating enhanced seepage and displacement analysis within unsaturated soil conditions is modelled based on a reported case history. A Monte Carlo simulation is used for rainfall infiltration and embankment deformation analysis under various rainfall scenarios. Probability values are obtained to achieve three levels of damage states in terms of road embankment settlement. The parametric analysis produces the exceedance probability curves for various rainfall intensities, saturated permeabilities and embankment slope angles. This work offers an efficient tool for assessing fragility to rainfall-induced excessive settlement of embankments.

To improve the stretch-bend forming precision of aluminum alloy profiles with complex open sections, the relationship stress–strain distribution and the spring-back radius variation were analyzed under different stretch-bend stages. Numerical simulation was performed on the auxiliary mold stretch-bend and spring-back process for the L-section aluminum alloy profile with variable curvature. The effects of cover plate pressure and baffle gap on residual stress, section deformation and spring-back value after forming were investigated. The results show that the stretch-bend warpage can be avoided by setting the cover plate on the straight section of the profile web plate. With the increase of the cover plate pressure, the section distortion of the vertical plate increases by 82%, and the maximum spring-back increases by 71%. A baffle plate was added to the outside of the vertical plate, and the baffle gap is proportional to the amount of section distortion and inversely proportional to the spring-back value. The critical cover force (1 kN) and the critical baffle gap (1 mm) obtained by numerical analysis are used in the stretch-bend forming experiment, and the maximum error between the two results is 2.39%. The numerical simulation has a certain reference value for the rational design of the auxiliary mold stretch-bend process parameters.

Due to the limited availability of fault samples and the expensive nature of marking fault samples in Electric Multiple Unit (EMU) gearbox monitoring data, a study was conducted to simulate the degradation process of key components in the CRH5 gearbox using rigid–flexible coupling dynamics. Vibration acceleration data from the simulation were utilized to create a six-dimensional hybrid feature domain representing the degradation process. By leveraging the capabilities of the Hidden Markov Model (HMM) for handling hidden transitive probabilities in temporal data and Gated Recurrent Unit (GRU) for addressing long-distance and high-dependence temporal data, a GRU-HMM fault prediction model was developed. This model was validated using monitoring data and the six-dimensional hybrid feature domain from the CRH5 gearbox and compared against actual maintenance records. The findings indicated that the GRU-HMM fault prediction model can effectively recognize the degradation patterns of multiple components, offering higher accuracy in fault prediction compared to traditional models. These research outcomes are expected to optimize EMU maintenance schedules based on usage conditions, enhance EMU utilization rates, and reduce operational and maintenance costs, thereby providing valuable theoretical support.

Abstract To provide backup and supplementation for the Global Navigation Satellite System (GNSS), Doppler shift from Low Earth Orbit (LEO) satellites can be used as signals of opportunity to provide positioning, navigation, and timing service. In this contribution, we first investigate the model and performance of instantaneous velocity determination and positioning with LEO satellites. Given a LEO constellation with 288 satellites, we simulate Doppler shift observations at nine multi-GNSS experiment stations. Owing to the lower orbit, the performance of LEO velocity determination is much more sensitive to the initial receiver position error than that of GNSS. Statistical results show that with the initial receiver position error increased from 0.1 to 10 m, the Root Mean Square Errors (RMSEs) increase from 0.73 to 2.65 cm/s, 0.68 to 2.96 cm/s, and 1.67 to 4.15 cm/s in the east, north, and up directions, respectively. The performances with GPS are compared with GPS + LEO, and it is found that LEO Doppler shift observations contribute to GPS velocity determination. As for LEO Doppler positioning, even if more than 30 visible LEO satellites are available, the position dilution of precision values can reach several hundreds. Assuming that the error of LEO Doppler measurements is 0.01 m/s, the instantaneous Doppler positioning accuracy can achieve about a few meters, which is comparable to that of GNSS pseudorange positioning. A constant velocity model is adopted for state transition. Static LEO Doppler positioning results show that an accuracy at centimeter to decimeter level can be achieved after solution convergence. For a static simulated kinematic positioning test, the RMSEs range from a few decimeters to several meters in different regions by giving different constraints. For a dynamic positioning test, the RMSEs are about 2–3 m in high latitude region.

Low-frequency vibration assisted bend-stretch forming has a good friction reduction effect due to the “surface effect” and enhances the forming ability of materials compared to traditional bend-stretch forming. In order to further improve the friction state of the forming contact interface, the contact interface friction between AL6005A aluminum alloy and Q345B steel under low-frequency vibration was analyzed by using the elastic-plastic contact theory. The formula for calculating the sliding friction coefficient of contact pressure and low-frequency vibration parameters was established. The influence of forming contact pressure and vibration frequency on the contact stress, friction force, and friction coefficient of asperity at the contact interface is obtained through the simulation analysis of low-frequency vibration sliding friction of double rough contact surface. The results show that the application of low-frequency vibration can weaken the influence of contact pressure on the friction coefficient and reduce the friction coefficient between AL6005A aluminum alloy and Q345B steel. Within a specific range of contact force, the sliding friction coefficient decreases first and then increases with the increase of vibration frequency. There is a critical vibration frequency to minimize the sliding friction coefficient of the contact interface.

The milling force of hollow thin-walled welded structural parts varies greatly at each position, and the analytical model was less studied, which affects the prediction accuracy of milling force. In order to construct an analytical model of the milling force, a model was constructed to calculate the cutting depth of the tool slot milling based on the contact relationship between the cutter and the workpiece. Moreover, based on the structural properties of the hollow thin-walled welded structural parts, the effective cutting depth function at different locations was derived by using infinitesimal element integration method, and the dynamic milling force prediction model was solved at each location. Furthermore, the differences between the theoretical and experimental values were compared and analyzed from the viewpoint of the peak value and the trend of the milling force, and the chip thickness was verified. The results show that the milling force is minimal in all directions at the position of the upper and lower skin zone, and maximal at the position of the inclined rib zone. What’s more, due to the large cutting depth at the position of the vertical plate zone, the milling force in X direction increases first and then decreases, and continue increases at the position of the vertical plate zone, besides, the milling force in Y and Z directions increases first and then decreases. What’s more, the prediction errors of the forces in X, Y and Z directions at the position of the upper and lower skin plate zone are 7.01%, 2.04% and 8.24%, respectively, the prediction errors of the forces in X, Y and Z directions at the position of the inclined rib plate zone are 3.22%、8.76% and 7.77%, respectively, and the prediction errors of the forces in X, Y and Z directions at the position of the vertical plate zone are 7.01%, 2.04% and 8.24%, respectively.

For the friction and wear of the braking surface of the disc brake, the mechanism of friction contact deformation of rough surface is analyzed. Considering elastoplastic contact mechanics of the friction surfaces, combined with Archard wear calculation theory, a wear calculation model of fractal rough surface is established, and a wear life prediction of the brake disc is proposed. Based on high-speed and heavy-load braking conditions, the influence of different braking parameters on the surface wear of brake disc is studied. Moreover, a pin-on-disc wear test is used to verify the validity of the wear calculation model, and the average relative error of test and simulation is 4.97%. The research results show that the fractal parameters affect the complexity of the rough surface, and surface wear is mainly caused by interface plastic contact. The friction coefficient is large, which can strengthen the surface shear effect. Large load and high speed would increase the real contact area of the surface and increase the wear volume per unit time. Furthermore, the wear life of the brake disc could be estimated, which is of great significance for raising the wear life of the brake disc, improving braking life and braking efficiency.

Gas hydrates have been considered as a new energy that could replace conventional fossil resources in the future because of their high energy density, environmental friendliness, and enormous reserves. To further analyze the potential distribution of gas hydrate stability zones (GHSZ) and the formation of a gas hydrate system in the Shenhu area of the South China Sea (SCS), a 3D petroleum simulation model (PSM) was built from 3D seismic interpretations and all available geological data. Based on the thermal calibration of the 3D model, the evolution of the GHSZ, hydrocarbon generation and migration, and the formation and accumulation of gas hydrates were simulated for the first time in the area. Thermal simulation shows that the methane source of gas hydrate originated from shallow biogenic gas and deep thermogenic gas. Most areas are dominated by shallow biogenic gas, while, only about 3% of the deep thermogenic gas derived from Enping Formation source rock and contributed to the gas hydrate formation within a few areas in the southeast. The thermogenic gas migrated vertically into the GHSZ through connecting faults, mud diapir, and/or gas chimney to form gas hydrate. The source rocks of the Wenchang Formation, a deep thermogenic gas source, began to enter the main hydrocarbon generation window at 28.4 Ma. The Enping source rock began to generate oil from 25 Ma on and gas from 16 Ma on. Since 5.3 Ma, most areas of the source rocks have generated a gas window, and only the shallower parts in the east still in the oil window, which had lasted until now. The shallow biogenic gas source rocks from the Hanjiang, Yuehai, and Wanshan formations generated gas in different periods, respectively. The Qionghai Formation began to generate hydrocarbon from 0.3 Ma and until now. Other results show that the GHSZ developed mainly during the Quaternary and Neogene (Wanshan Formation) and the GHSZ is thicker in the southern area and thinner in the northern part with a positive correlation with water depth. Starting at 11.6 Ma, the GHSZ developed in the Hanjiang Formation in the south of the Shenhu area and gradually expanded to the north to cover most of the study area at 5.3 Ma during the Yuehai Formation. From 1.8 Ma on, the GHSZ covered the entire study area. At the same time, the GHSZ in the Hanjiang Formation disappeared because of the change in temperature and pressure. At present, the GHSZ in the Yuehai Formation has disappeared, while the Quaternary and Wanshan are the two main formations for GHSZ development. The formation and distribution of gas hydrates are fundamentally controlled by the space-time coupling between the hydrocarbon generation and expulsion time and distribution of the GHSZ. The simulation results of gas hydrate accumulation and distribution were verified by drilling results and the matching rate is 84%. This is the first time that 3D simulation was successfully conducted with PSM technology in the Shenhu area and it provides important guidance for gas hydrate study in other areas of the SCS.

Titanium alloys are widely used in aerospace, the military industry, electronics, automotive fields, etc., due to their excellent properties such as low density, high strength, high-temperature resistance, and corrosion resistance. Many components need to be finished precisely after being cut in these applications. In order to achieve high-quality and high-efficiency finishing of titanium alloys, ultrasonic-assisted magnetic abrasive finishing (UAMAF) was introduced in this research. The machining performance for UAMAF of a titanium alloy was studied by experimentally comparing UAMAF and magnetic abrasive finishing (MAF). The results show that the cutting force of UAMAF can reach 2 to 4 times that of MAF, and it decreases rapidly with the increase in the machining gap due to the energy loss of ultrasonic impact in the transmission between magnetic abrasives. The surface roughness of UAMAF can reach about Ra 0.075 μm, which is reduced by about 59% compared with MAF. The main wear type of the magnetic abrasive is that the diamond grits fell off the magnetic abrasive in both UAMAF and MAF. The uniform wear of the magnetic abrasive is realized, and the utilization ratio of the magnetic abrasive is obviously improved in UAMAF.

For material fracture and severe wear in braking conditions, the discrete element method (DEM) is used to simulate the wear process of the braking interface explicitly. Based on the central difference method, particle motion equations are established considering the influence of elemental damping on particle contact. Combined with the Particle Flow Code (PFC) software, a DEM wear model of the braking interface is established using the parallel bond modeling method. The braking wear process is simulated, and the material damage process is investigated. The simulation results demonstrate that with the increase of the initial braking load and the initial braking speed, the wear depth increased by 24.75% and 16.22%, respectively. The increase in the number of detached particles leads to an increasing trend of fracture force chains, which increases the thickness of the flowing particle layer, revealing the micro-wear mechanism of the braking interface.