Shear Stress Research Articles

A novel method for investigating ice-substrate interfacial rheology and long-term adfreeze strength through loading-unloading creep test

Ice cementation and ice-substrate adfreeze force are the primary contributors to the high bearing capacity of pile foundations in cold regions and the stability of frozen walls in areas subjected to artificial freezing. Given the significant temperature sensitivity of ice’s shear rheology, engineering structures in ice or ice-rich soils continue to deform even under constant external loads. A thorough understanding of shear creep and the long-term adfreeze force at the ice-substrate interface is essential for predicting the continuous deformation of these structures. However, research into the shear creep behavior at frozen interfaces has historically been constrained by the precision of temperature control in experimental settings and the complexity of load paths in shear testing devices. In this study, a temperature- and stress-control device for interface shear creep is assembled firstly, and multilevel loading-unloading creep tests on steel pipes embedded in layered frozen ice were conducted. Through the decoupling of deformation progression, the viscoelastic and viscoplastic shear behaviors at the steel-ice interface under various temperatures and shear stresses were characterized, the principle of sustainable interfacial shear creep along with its underlying physical mechanism were proposed. Subsequently, with the aid of a modified nonlinear Burger model, various interfacial shear creep parameters were derived. Results reveal that the interfacial generalized shear modulus continuously improves but with a gradually weakening degree until a point of accelerating creep is reached. Additionally, the long-term adfreeze force is found to be less than half of the short-term strength, which significantly decreases as the temperature approaches the water phase transition zone. Interestingly, the stress exponent associated with the interfacial steady creep rate is considerably smaller than that predicted by Glen’s law. This research provides a theoretical basis instrumental in the engineering design in cold regions and those structures employing artificial freezing techniques.

Comparative Analysis of Circular RNAs Expression and Function between Aortic and Intracranial Aneurysms.

An aneurysm is an abnormal enlargement or bulging of the wall of a blood vessel. Most often, aneurysms occur in large blood vessels - the aorta (Thoracic Aortic Aneurysm (TAA) and Abdominal Aortic Aneurysm (AAA) and brain vessels (Intracranial Aneurysm (IA)). Despite the presence of significant differences in the pathogenesis of the development and progression of IA and TAA/AAA, there are also similarities. For instance, both have been shown to be strongly influenced by shear stress, inflammatory processes, and enzymatic destruction of the elastic lamellae and extracellular matrix (ECM) proteins of the vascular wall. Moreover, although IA and TAA are predominantly considered arteriopathies with different pathological mechanisms, they share risk factors with AAA, such as hypertension and smoking. However, there is a need for a more in- -depth study of the key elements that may influence the formation and progression of a particular aneurysm to find ways of therapeutic intervention or search for a diagnostic tool. Today, it is known that the disruption of gene expression is one of the main mechanisms that contribute to the development of aneurysms. At the same time, growing evidence suggests that aberrant epigenetic regulation of gene function is strongly related to the genesis of aneurysms. Although much has been studied of the known protein-coding genes, circular RNAs (circRNAs), a relatively new and rapidly evolving large family of transcripts, have recently received much scientific attention. CircRNAs regulate gene expression through the sponging of microRNAs (miRNAs) and can also be used as therapeutic targets and biomarkers. Increasing evidence has implicated circRNAs in the pathogenesis of multiple cardiovascular diseases, including the development of aneurysms. However, the mechanism of dysregulation of certain circRNAs in a particular aneurysm remains to be studied. The discovery of circRNAs has recently advanced our understanding of the latest mode of miRNAs/target genes regulation in the development and progression of IA and TAA/AAA. The aim of this study is to compare the expression profiles of circRNAs to search for similar or different effects of certain circRNAs on the formation and progression of IA and TAA/AAA.

Analysis of composite honeycomb sandwich panels under blast loads

ABSTRACT This paper presents a novel isogeometric model to analyze the static bending, fre, and forced vibration characteristics of a sandwich of two-layer plates made of auxetic honeycomb core and laminated three-phase polymer/GNP/fiber skin subjected to a double blast load. The two layers are restricted from vertical movement by elastic pins but they can still move in the horizontal plane. The unique feature of the refined first-order shear plate theory with a distribution function g(z) is used to replace the shear correction factor, thereby allowing the shear stress at the top and bottom surfaces to be accurately calculated. The results are derived by the exact Navier solutions for rectangular plates and the IGA method with arbitrary boundary conditions for elliptical and rectangular plate structures. A comprehensive investigation was conducted to evaluate the influence of some parameters on the oscillation of sandwich two-layer plates subjected to double blast load.

In silico investigation of the reliable models of drug delivery from blood vessels to tumor using tumor spheroid-on-chip models

Due to the limitations of two-dimensional models in mimicking intratumoral drug delivery, three-dimensional tumor spheroid-on-chip models have been considered recently. Even though various models have been proposed in this field, many unknowns still exist about their reliability in imitating in vivo conditions. In addition, most of these models neglect the replication of drug transport from the bloodstream. In the current research, considering the interface of tumor and blood vessels, new designs of microfluidic devices are developed to resemble the in vivo conditions more closely. Then, drug transport in the proposed models is studied numerically using porous media theory under dynamic conditions. The simulation results of drug concentration, cell viability, fluid shear stress, and pressure distributions are compared with those of a solid tumor. As the results demonstrate, our proposed models can well mimic the in vivo conditions, including drug delivery through blood vessels. Therefore, in addition to developing tumor spheroid-on-chip models resembling tumor-vessel interaction, the present study provides a strong strategy to evaluate the reliability of in vitro models in replicating intratumoral drug delivery.

A vertically non-uniform temperature approach for the friction term computation in depth-averaged viscoplastic lava flows

Recently, depth-averaged shallow flow models have been adapted to modelling liquefied lava flows, generally characterized by a marked temperature-dependent non-Newtonian rheology. Modelling these complex flows requires to include the effects of the depth-averaged temperature gradients in the governing equations. The most significant term to correctly predict the lava mobility is the flow resistance term, which is widely estimated using the linear viscoplastic Bingham model. This non-Newtonian model allows to relate the bed shear stress to the depth-averaged lava flow features by means of a cubic equation with analytical solution when assuming a uniform temperature distribution along the vertical profile. Nevertheless, the lava temperature is non-uniform along the vertical due to the heat transfer at the bottom and the free surface, and hence the classical cubic Bingham model is not valid anymore. In this work, a depth-averaged shallow flow model is adapted for realistic lava flows considering influence of the non-uniform vertical temperature profile in the non-Newtonian resistance. This requires to modify the rheological viscoplastic models for ensuring the coupling between flow dynamics and temperature evolution. Three non-uniform temperature vertical distributions are considered: linear, piece-wise and diffusion profiles. Synthetic tests are used to show the influence of the temperature vertical profile on the numerical results. Furthermore, laboratory experimental data are used to validate this novel viscoplastic resistance formulation and to show that the calibration of its parameters is possible.

Variable stiffness performance analysis of layer jamming actuator based on bionic adhesive flaps

Soft actuators made of soft materials cannot generate precisely efficient output forces compared to rigid actuators. It is a promising strategy to equip soft actuators with variable stiffness modules of layer jamming mechanism, which could increase their stiffness as needed. Inspired by the gecko's the array of setae, bionic adhesive flaps with inclined micropillars are applied in layer jamming mechanism. In this paper, after the manufacturing process of the layer jamming actuator based on the bionic adhesive flaps is described, the equivalent stiffness models of the whole actuator are established in the unjammed and jammed states. And the shear adhesive force of a single micropillar is calculated based on the Kendall viscoelastic band model. The finite element simulation results of two bionic adhesive flaps show that the interlaminar shear stress and stiffness increase with the increase of pressure. The measurement of shear adhesive force show that the critical shear adhesive force of the bionic adhesive material is 3.2 times that of polyethylene terephthalate (PET) material, and exhibit the ability of anisotropic adhesion behavior. The variable stiffness performance of the layer jamming actuator based on bionic adhesive flaps is evaluated by three test methods, and the max stiffness reaches 8.027 N mm-1, which is 1.5 times higher than the stiffness of the layer jamming actuator based on the PET flaps. All results of simulation and experiment effectively verify the validity and superiority of applying the bionic adhesive flaps to the layer jamming mechanism to enhance the stiffness.

DEM analysis of cohesive granular shear flow using dynamic adhesion force model – Model validation for contact-dominated regime

When the DEM simulations are performed, the particle stiffness of the Discrete Element Method (DEM) is often reduced to reduce the computational cost. The Dynamic Adhesion/Cohesion Force Model (DAFM/DCFM) was proposed for adhesive/cohesive particles to make the effect of adhesion/cohesion force on the collisional motion of a particle with a reduced particle stiffness equivalent to that of the original system.This study validates the applicability of DCFM to the contact-dominated regime by means of a series of DEM simulations of two-dimensional simple shear flow of cohesive spherical particles. The results demonstrate that DCFM accurately represents the cohesive nature of the original system with respect to the fluid-to-solid phase transition of granular medium due to a reduction in shear rate, as well as the shear stress of the granular medium. It is demonstrated that the bond-breaking model proposed in the present study describes well the validity of DCFM in the contact-dominated regime.

Experimental investigation of permeability evolution and progressive damage of sandstone in the complete process of direct shear

Permeability evolution and damage behavior under shearing conditions are critical for reliable predictions of rock seepage characteristics during excavation and a correct assessment of their impacts on the occurrence of natural or induced hazards. However, publications concerned with this problem have been poorly issued so far. The primary objective of this study is to innovatively investigate permeability evolution and damage behaviors of sandstone in the complete process of direct shear. Several groups of red sandstone specimens are deformed in direct shearing using a newly modified servo-controlled direct shear apparatus to measure permeability evolution and monitor acoustic emission (AE) activities with damage under various normal stress conditions. The morphology of the rupture surfaces is scanned using a 3D laser scanner after tests. According to our experimental results, the locus of the shear stress-displacement curve is found to be closely associated with microcrack progressivity. The initial permeability in advance of shearing test and the peak permeability appearing after the peak shear stress decrease with the increasing normal stress. The minimal and maximum permeability values exhibit a lag-effect with respect to the turning and the peak points on the shear stress-displacement curve. Furthermore, a palpable surge is found in the interevent time function during shearing, which can be a potential indicator for early warning of catastrophic rupture of brittle rock materials under shearing or compression. The joint roughness coefficient (JRC) values in the direction of shearing deformation for each normal stress case are more significant than those in the direction perpendicular to the shear displacement. The impacts of our obtained results are potential ramifications for permeability control in rock engineering projects.

Mathematical derivation of a unified equations for adjoint lattice Boltzmann method in airfoil and wing aerodynamic shape optimization

Unified equations of the adjoint lattice Boltzmann method (ALBM) are derived for five applicable objective functions in 2D/3D aerodynamic shape optimization problems. The derived equations include the adjoint equation, boundary condition, terminal condition and gradient of the cost function. In this research, firstly, these relations are extracted for each objective in details and then the general form of ALBM equations are presented for all defined practical aerodynamic objective function. Five applicable cost functions which are the most important objectives in optimization of aerodynamic geometries include desired pressure and viscous shear stress (VSS) inverse design, drag and moment at fixed lift and finally lift to drag ratio at fixed angle of attack. The new extracted relations are based on the circular and spherical function scheme, and are valid for viscous/inviscid, compressible/incompressible and 2D/3D flows in all continuous flow regimes. Proof of new extracted general relations have been performed by authors.

Fluid-structure interaction analysis of a healthy aortic valve and its surrounding haemodynamics.

The opening and closing dynamics of the aortic valve (AV) has a strong influence on haemodynamics in the aortic root, and both play a pivotal role in maintaining normal physiological functions of the valve. The aim of this study was to establish a subject-specific fluid-structure interaction (FSI) workflow capable of simulating the motion of a tricuspid healthy valve and the surrounding haemodynamics under physiologically realistic conditions. A subject-specific aortic root was reconstructed from magnetic resonance (MR) images acquired from a healthy volunteer, whilst the valve leaflets were built using a parametric model fitted to the subject-specific aortic root geometry. The material behaviour of the leaflets was described using the isotropic hyperelastic Ogden model, and subject-specific boundary conditions were derived from 4D-flow MR imaging (4D-MRI). Strongly coupled FSI simulations were performed using a finite volume-based boundary conforming method implemented in FlowVision. Our FSI model was able to simulate the opening and closing of the AV throughout the entire cardiac cycle. Comparisons of simulation results with 4D-MRI showed a good agreement in key haemodynamic parameters, with stroke volume differing by 7.5% and the maximum jet velocity differing by less than 1%. Detailed analysis of wall shear stress (WSS) on the leaflets revealed much higher WSS on the ventricular side than the aortic side and different spatial patterns amongst the three leaflets.

Rheological properties of gangue-filled slurry based on L-pipe experiments

ABSTRACT The examination of rheological properties of gangue slurry is a crucial precondition for the application of gangue slurry in pipeline transportation. This study focused on the investigation of the rheological properties and hydraulic gradient of gangue slurry in a mining region located in northern Shaanxi. By employing the L-pipe test method, the research examined the impact of gangue particle size, slurry concentration, and conveying pipe diameter on the aforementioned characteristics. The findings indicated a positive correlation between the hydraulic gradient of the mine gangue slurry and the concentration of the slurry, whereas a negative correlation existed between the hydraulic gradient and both the particle size and pipe diameter. Additionally, the slurry flow rate was negatively connected with both the slurry concentration and particle size. The concentration of gangue-filled slurry had a positive correlation with both the shear stress and viscosity, whereas the particle size demonstrated a negative correlation with both properties. The optimal slurry concentration for conveying is between 70% and 75%. Furthermore, the ideal slurry particle size ratio is either 10:0 or 9:1, and the recommended pipe diameter is 125mm. The research findings offer theoretical support for the judicious selection of transport parameters when conveying gangue slurry.

Advances in Shear Stress Stimulation of Stem Cells: A Review of the Last Three Decades.

Stem cells are widely used in scientific research because of their ability to self-renew and differentiate into a variety of specialized cell types needed for body functions. However, the self-renewal and differentiation of stem cells are regulated by various stimuli, with mechanical stimulation being particularly notable due to its ability to mimic the physical environment in the body. This study systematically collected 2638 research papers published between 1994 and 2024, employing tools such as VOSviewer, CiteSpace, and GraphPad Prism to uncover research hotspots, publication trends, and collaboration networks. The results indicate a yearly increase in global research on the shear stress stimulation of stem cells, with significant contributions from the United States and China in terms of research investment and output. Future research directions include a deeper understanding of the mechanisms underlying mechanical stimulation's effects on stem cell differentiation, the development of new materials and scaffold designs to better replicate the natural cellular environment, and advancements in regenerative medicine. Despite considerable progress, challenges remain in translating basic research findings into clinical applications.

Construction nanobiotechnology approach for performance enhancement of microbially induced biomineralization (MIB) using a biopolymer encapsulated spore-based system.

The integration of green construction practices within the built environment has been significantly advanced by biotechnological innovations, among which microbially induced biomineralization (MIB), predominantly facilitated by various strains of spore-forming bacilli, emerges as a pivotal mechanism for the self-healing of concrete. However, the practical deployment of this technology faces challenges, notably the compromised viability of bacterial spores due to germination triggered by severe shear stress during concrete mixing. To address this limitation, a water-insoluble polymer (extracellular polymeric substance) produced by Cellulomonas flavigena was utilized to encapsulate and protect the spores. The encapsulation process was rigorously verified through physicochemical methodologies, including X-ray diffraction (XRD) analysis, which revealed alterations in the interlayer spacings of the extracellular polymeric substance (EPS) structure during the encapsulation process, indicating successful EPS coating of the spores. Furthermore, a proof of concept for the enhanced biomineralization capacity of EPS-coated spores was demonstrated. Standard analytical techniques confirmed the precipitation of calcite and vaterite among other minerals, underscoring the effectiveness of this novel approach. This breakthrough paves the way for the development of innovative, sustainable bioconcrete applications, aligning with broader environmental objectives and advancing the field of green construction technology.IMPORTANCEDevelopment of bioconcrete with self-healing capability through MIB constitutes an important sustainable construction biotechnology approach for restoration and repair of built environment. Like every promising technology, MIB also suffers from certain shortcomings in terms of compromised viability of the microbial cells after premature germination of the spores on exposure to shear stress caused during concrete mixing. In this study, these challenges were adequately addressed by successfully providing a protective coating of indigenously extracted EPS to the bacterial spores and elucidating the interactive mechanisms between them. The results showed stable encapsulation of the spores while providing mechanistic insights of the encapsulation phenomenon. The data also showed enhanced rate of biomineralization by encapsulated microbes when subjected to stress conditions.

Precise analysis of the static bending response of functionally graded doubly curved shallow shells via an improved finite element shear model

Doubly curved shallow shells, common in aerospace and civil engineering, pose significant challenges in predicting bending responses due to their complex geometry and material properties. Therefore, this research focused on the development of an efficient eight-node quadrilateral isoparametric element model based on a new improved first-order shear deformation theory (IFSDT) to analyze the bending deflection and stresses distribution of Functionally Graded (FG) doubly curved shallow shells. The present IFSDT obviates the necessity of a shear correction factor, which has challenges in determining for FG doubly curved shallow shells, by adopting a simpler assumption about transverse shear stresses. This assumption guarantees that the model precisely predicts the parabolic shear stresses distribution throughout the thickness of the FG shell while also meeting the requirements for free traction conditions on the top and lower surfaces of the shell. In the present analysis, five distinct types of FG doubly curved shallow shells are modeled: flat plates, spherical shells, hyperbolic parabolic shells, cylindrical shells, and elliptical paraboloid shells. A power law is employed to describe the gradual changes in material properties through the thickness direction. Numerical results for displacements and stresses are obtained for various shell types, materials, power-law factors, radii of curvature, and boundary conditions under uniform and sinusoidal loads. Comprehensive comparison and convergence studies are conducted to assess the robustness and accuracy of the proposed numerical model. This study also contributes new numerical results that are valuable for future research. Ultimately, the proposed numerical model offers a robust and straightforward tool for engineers and researchers engaged in the design and analysis of complex shell structures.

Label-Free Microfluidic Apheresis of Circulating Tumor Cell Clusters.

Screening liters of blood (i.e., apheresis) represents a generalized approach to promote the reliable access to circulating tumor cell clusters (CTCCs), which are known to be highly metastasis-competent, yet ultrarare. However, no existing CTCC sorting technology has demonstrated high throughput, high yield, low shear stress, and minimal blood dilution simultaneously as required in apheresis. Here, a label-free method is introduced termed Precision Apheresis for Non-invasive Debulking of cell Aggregates (PANDA) to continuously isolate CTCCs from undiluted blood to clean buffer through size sorting, processing 1.4 billion cells per second. The cell focusing is optimized within whole blood leveraging secondary transverse flow and margination. The PANDA chip recovers >90% of spiked ≈24 rare HeLa cell clusters from 100mL undiluted blood samples (equivalent to ≈500 billion blood cells) at 1 Lh-1 throughput, with ≤20s device residence time, ≤15Pa shear stress, and >99.9% return of blood components. The technology lays the groundwork for future routine isolation to increase the recovery of these ultrarare yet clinically significant tumor cell populations from large volumes of blood to advance cancer research, early detection, and treatment.

Strength Model for Prestressed Concrete Beams Subjected to Pure Torsion

A torsional strength model for prestressed concrete beams was proposed considering the initial crack angle, principal stress angle, and longitudinal strain, which are affected by the axial stress induced by the effective prestress. The use of the torsional effective thickness was also proposed to calculate the torsional strength of prestressed concrete beams by considering the effect of prestress. The shear element in the torsional member was simplified under the assumption that the principal tensile stress and principal compressive strain were negligible in the ultimate state. The torsional strength was determined when the principal compressive stress or shear stress at the crack surface in the shear element reached the failure criterion according to the multipotential capacity model, which considers concrete crushing and aggregate interlocking as the main resistances to the applied load. The proposed strength model was verified using test specimens collected from existing experimental studies. The proposed model accurately evaluated the torsional strength of prestressed concrete beam specimens, regardless of the key variables of the prestressed concrete specimens, where the mean value of the tested results to the calculated torsional strengths was 1.123, and the corresponding coefficient of variation was 17.7% for 104 prestressed concrete beam specimens, while the ACI 318-19 torsional design method gave the mean and coefficient of variation of 0.880 and 24.3%, respectively.

Low-velocity impact behavior of foam-based sandwich composite reinforced with warp-knitted spacer fabric; numerical and experimental study

The aim of this research is to investigate experimentally and numerically the low-velocity impact behavior of foam-based composites reinforced with warp-knitted spacer fabric (WKSF). To prepare different foam-based composites, the structural parameters of WKSFs including cell size, position, and Z-fiber height were considered. A drop weight impact test with an initial energy of 5J was carried out to examine the low-velocity impact behavior of composites, followed by experimental analyses of Mises, shear, and normal stress on various composite components. Thereafter, the impact behavior of the composites was simulated using ABAQUS/CAE software. The comparison between experimental and numerical results showed a maximum error of 9.79% in predicting the acceleration of impactor. In addition, the results revealed significant stress disparities among samples. Stress analysis showed complex patterns across samples, emphasizing structural parameter influence on stress tolerance and load-bearing capabilities. Notably, Z-fibers displayed substantial stress tolerance, while the matrix predominantly undergoes shear stress. Consequently, the ideal structure for low-velocity impact applications includes small cell size, high thickness, and non-facing hexagonal cells.

Effect of processing and formulation factors on Catalase activity in tablets

The manufacturing of tablets containing biologics exposes the biologics to thermal and shear stresses, which are likely to induce structural changes (e.g., aggregation and denaturation), leading to the loss of their activity. Saccharides often act as stabilizers of proteins in formulations, yet their stabilizing ability throughout solid oral dosage processing, such as tableting, has been barely studied. This work aimed to investigate the effects of formulation and process (tableting and spray-drying) variables on catalase tablets containing dextran, mannitol, and trehalose as potential stabilizers. Non-spray-dried and spray-dried formulations were prepared and tableted (100, 200, and 400 MPa). The enzymatic activity, number of aggregates, reflecting protein aggregation and structure modifications were studied. A principal component analysis was performed to reveal underlying correlations. It was found that tableting and spray-drying had a notable negative effect on the activity and number of aggregates formed in catalase formulations. Overall, dextran and mannitol failed to preserve the catalase activity in any unit operation studied. On the other hand, trehalose was found to preserve the activity during spray-drying but not necessarily during tableting. The study demonstrated that formulation and process variables must be considered and optimized together to preserve the characteristics of catalase throughout processing.

Statistical analysis of adhesive rod-tube joints under tensile stress for structural applications

Adhesive bonding has been replacing traditional joining methods such as welding, bolting, and riveting in the design of mechanical structures in the automotive, aerospace and aeronautic industries. This joining method has several advantages over traditional methods such as ease of manufacture, lower costs, ease of joining different materials, higher fatigue resistance, and high corrosion resistance. Although tubular adhesive joints have varying applications, such as in truss structures and vehicles, machine axles, and piping, different joint configurations exist, such as rod-tube joints (RTJ), which are not conveniently addressed in the literature. This work compares the tensile performance of adhesively bonded RTJ between aluminium alloy components (AW6082-T651), considering the variation of the main geometric parameters: overlap length (LO), tube thickness (tS), rod diameter (d), adhesive fillet angle (f), and type of adhesive. The Taguchi’s method was employed in the elaboration of the applied design of experiments (DoE). To compare the RTJ behaviour, a numerical analysis was carried out through finite element analysis (FEA) and cohesive zone modelling (CZM). Peel (σy) and shear (τxy) stresses in the adhesive layer were initially obtained by applying purely elastic models. CZM modelling made possible to obtain the damage evolution in the adhesive layer, the maximum load (Pm) and dissipated energy (U) at Pm of the adhesive joints. As a result of applying the Taguchi method, the adhesive joint that showed the best overall performance used the adhesive Araldite® AV138, LO = 40 mm, d = 20, and tS = 3 mm.

A 44-year balanced fertilizer application affected rill erosion resistance by changing humus, aggregates, and polyvalent cation

Fertilizer application can affect the physicochemical properties of soil, such as the contents of large aggregates, humus, and exchangeable cations, thereby influencing soil erosion resistance. However, the rill erosion resistance and its key influencing factors of soil following long-term balanced fertilizer application remains unclear. This study aimed to analyze the change in rill erosion resistance following 44 years of balanced fertilizer application (No changes in the type of fertilizer) to non-calcareous soils and establish a Partial Least Squares Regression (PLSR) model of rill erodibility (kd) and soil critical shear stress (τc) to changes in the physical and chemical properties of the soil. Five treatments were designed: (1) CK (no fertilizer applied), (2) N (nitrogen), (3) NP (nitrogen plus phosphorus), (4) NK (nitrogen plus potassium), and (5) NPK (nitrogen, phosphorus, and potassium). Compared to CK, the N, NP, NK, and NPK application significantly increased τc by 49.6, 96.7, 73.6, and 36.2 %, respectively. Whereas, kd increased significantly only in the NPK treatment group. The optimal partial least squares regression model showed that mean weight diameter (MWD) had the greatest positive influence on soil critical shear stress, followed by fulvic acid (FA) content, whereas water-dissolved substances had a negative influence. Long-term balanced fertilizer application can increase MWD and τc by combining micro-aggregates and humus into large aggregates. Ca2+ content had the greatest positive effect on kd. Compared with that of CK, exchangeable Ca2+ content increased significantly with NP and NPK application (7.9 and 34.3 %, respectively). Ca2+ can increase kd by binding to the polar functional groups in FA to promote the shedding of hydration shells in large aggregates. Among all treatments, the NP treatment showed the best performance for reducing kd and increasing τc. This study could contribute to the understanding of the rill erosion process and modeling in non-calcareous soils and offer a reference for agricultural erosion control treatments.