The finite-element modeling and simulations of the intra-body communication （IBC） were investigated to provide a theoretical basis for biomedical monitoring. A finite-element model for the whole human body was devel...The finite-element modeling and simulations of the intra-body communication （IBC） were investigated to provide a theoretical basis for biomedical monitoring. A finite-element model for the whole human body was developed to simulate the IBC. The simulation of galvanic coupling IBC and electrostatic coupling IBC were implemented along with different signal transmission paths, and their attenuations were calculated. Our study showed that the position near the signal electrode had higher potential than other positions in the two types of IBC, while the potential generally decreased along the axis of the body parts. Both signal attenuations of the two types IBC increased with increasing signal transmission distance, and the electrostatic coupling IBC had comparatively higher receiving potential than the galvanic coupling IBC. The results indicated that the proposed modeling method could be used for the research of biomedical monitoring based on IBC technology.展开更多
Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and ...Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and with the nucleus removed were developed and validatedusing published experimental and clinical data. Then the models with a stainless steel nucleusprosthesis implanted and with polymer nucleus prostheses of various properties implanted were usedfor the 3D finite-element biomechanical analysis. All the above simulation and analysis were carriedout for the L4/L5 disc under a human worst--daily compression load of 2000 N. The results show thatthe polymer materials with Young's modulus of elasticity E = 0.1-100 MPa and Poisson's ratio v=0.35-0.5 are suitable to produce artificial nucleus prosthesis in view of biomechanicalconsideration.展开更多
The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor...Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor CSAMT forward modeling using the vector finite-element method. To verify the feasibility of the algorithm, we calculate the electric field, magnetic field, and tensor impedance of the 3D CSAMT far-zone field in layered media and compare them with theoretical solutions. In addition, a three-dimensional anomaly in half-space is also simulated, and the response characteristics of the impedance tensor and the apparent resistivity and impedance phase are analyzed. The results suggest that the vector finite-element method produces high-precision electromagnetic field and impedance tensor data, satisfies the electric field discontinuity, and does not require divergence correction using the vector finite-element method.展开更多
The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such...The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.展开更多
Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficien...Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.展开更多
In recent years,there has been an increasing need for climate information across diverse sectors of society.This demand has arisen from the necessity to adapt to and mitigate the impacts of climate variability and cha...In recent years,there has been an increasing need for climate information across diverse sectors of society.This demand has arisen from the necessity to adapt to and mitigate the impacts of climate variability and change.Likewise,this period has seen a significant increase in our understanding of the physical processes and mechanisms that drive precipitation and its variability across different regions of Africa.By leveraging a large volume of climate model outputs,numerous studies have investigated the model representation of African precipitation as well as underlying physical processes.These studies have assessed whether the physical processes are well depicted and whether the models are fit for informing mitigation and adaptation strategies.This paper provides a review of the progress in precipitation simulation overAfrica in state-of-the-science climate models and discusses the major issues and challenges that remain.展开更多
(Quasi-)closed-form results for the statistical properties of unmanned aerial vehicle(UAV)airto-ground channels are derived for the first time using a novel spatial-vector-based method from a threedimensional(3-D)arbi...(Quasi-)closed-form results for the statistical properties of unmanned aerial vehicle(UAV)airto-ground channels are derived for the first time using a novel spatial-vector-based method from a threedimensional(3-D)arbitrary-elevation one-cylinder model.The derived results include a closed-form expression for the space-time correlation function and some quasi-closed-form ones for the space-Doppler power spectrum density,the level crossing rate,and the average fading duration,which are shown to be the generalizations of those previously obtained from the two-dimensional(2-D)one-ring model and the 3-D low-elevation one-cylinder model for terrestrial mobile-to-mobile channels.The close agreements between the theoretical results and the simulations as well as the measurements validate the utility of the derived channel statistics.Based on the derived expressions,the impacts of some parameters on the channel characteristics are investigated in an effective,efficient,and explicable way,which leads to a general guideline on the manual parameter estimation from the measurement description.展开更多
In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical propert...In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical properties of rocks,the cracking processes of pre-cracked rocks have been extensively studied using numerical modeling methods.The peridynamics(PD)exhibits advantages over other numerical methods due to the absence of the requirements for remeshing and external crack growth criterion.However,for modeling pre-cracked rock cracking processes under impact,current PD implementations lack generally applicable rock constitutive models and impact contact models,which leads to difficulties in determining rock material parameters and efficiently calculating impact loads.This paper proposes a non-ordinary state-based peridynamics(NOSBPD)modeling method integrating the Drucker-Prager(DP)plasticity model and an efficient contact model to address the above problems.In the proposed method,the Drucker-Prager plasticity model is integrated into the NOSBPD,thereby equipping NOSBPD with the capability to accurately characterize the nonlinear stress-strain relationship inherent in rocks.An efficient contact model between particles and meshes is designed to calculate the impact loads,which is essentially a coupling method of PD with the finite element method(FEM).The effectiveness of the proposed NOSBPD modeling method is verified by comparison with other numerical methods and experiments.Experimental results indicate that the proposed method can effectively and accurately predict the 3D cracking processes of pre-cracked cracks under impact loading,and the maximum principal stress is the key driver behind wing crack formation in pre-cracked rocks.展开更多
Glassy polymers are widely used in biomedical applications in a solvent environment,yet their long-term performance is governed by the competing effects of physical aging and solvent-induced plasticization.Here,we dev...Glassy polymers are widely used in biomedical applications in a solvent environment,yet their long-term performance is governed by the competing effects of physical aging and solvent-induced plasticization.Here,we develop a constitutive model that explicitly couples the solvent concentration,structural relaxation,and mechanical response.This framework is built on a multiplicative decomposition of deformation and an Eyring-type flow rule,with structural evolution described by an effective temperature.A generalized shift factor is introduced to quantify how the solvent concentration and effective temperature jointly affect the relaxation time,thereby integrating physical aging and plasticization.The model is subsequently applied to methacrylate(MA)-based copolymer networks immersed in phosphate-buffered saline for up to nine months.Simulations accurately capture key experimental features,including the strong softening of highly swellable networks,the partial recovery due to aging,and the mitigating role of hydrophobic crosslinking in reducing solvent uptake.While the current single-mode description cannot reproduce the full relaxation spectrum,it establishes an efficient framework for predicting the long-term mechanical performance under coupled environmental and mechanical loading.This study provides a constitutive description of solvent-swollen glassy polymers,offering mechanistic insight into the interplay between plasticization and aging.Beyond biomedical MA networks,this framework establishes a foundation for predicting the long-term performance of polymer glasses under coupled aqueous environmental and mechanical loading.展开更多
Activation pruning reduces neural network complexity by eliminating low-importance neuron activations,yet identifying the critical pruning threshold—beyond which accuracy rapidly deteriorates—remains computationally...Activation pruning reduces neural network complexity by eliminating low-importance neuron activations,yet identifying the critical pruning threshold—beyond which accuracy rapidly deteriorates—remains computationally expensive and typically requires exhaustive search.We introduce a thermodynamics-inspired framework that treats activation distributions as energy-filtered physical systems and employs the free energy of activations as a principled evaluation metric.Phase-transition-like phenomena in the free-energy profile—such as extrema,inflection points,and curvature changes—yield reliable estimates of the critical pruning threshold,providing a theoretically grounded means of predicting sharp accuracy degradation.To further enhance efficiency,we propose a renormalized free energy technique that approximates full-evaluation free energy using only the activation distribution of the unpruned network.This eliminates repeated forward passes,dramatically reducing computational overhead and achieving speedups of up to 550×for MLPs.Extensive experiments across diverse vision architectures(MLP,CNN,ResNet,MobileNet,Vision Transformer)and text models(LSTM,BERT,ELECTRA,T5,GPT-2)on multiple datasets validate the generality,robustness,and computational efficiency of our approach.Overall,this work establishes a theoretically grounded and practically effective framework for activation pruning,bridging the gap between analytical understanding and efficient deployment of sparse neural networks.展开更多
Kinetic impact is the most practical planetary-defense technique,with momentum-transfer efficiency central to deflection design.We present a Monte Carlo photometric framework that couples ejecta sampling,dynamical evo...Kinetic impact is the most practical planetary-defense technique,with momentum-transfer efficiency central to deflection design.We present a Monte Carlo photometric framework that couples ejecta sampling,dynamical evolution,and image synthesis to compare directly with HST,LICIACube,ground-based and Lucy observations of the DART impact.Decomposing ejecta into(1)a highvelocity(~1600 m/s)plume exhibiting Na/K resonance,(2)a low-velocity(~1 m/s)conical component shaped by binary gravity and solar radiation pressure,and(3)meter-scale boulders,we quantify each component’s mass and momentum.Fitting photometric decay curves and morphological evolution yields size-velocity distributions and,via scaling laws,estimates of Dimorphos’bulk density,cratering parameters,and cohesive strength that agree with dynamical constraints.Photometric ejecta modeling therefore provides a robust route to constrain momentum enhancement and target properties,improving predictive capability for kinetic-deflection missions.展开更多
Xylogenesis,the process through which wood cells are formed,results in the long-term storage of carbon in woody biomass,making it a key component of the global carbon cycle.Understanding how environmental drivers infl...Xylogenesis,the process through which wood cells are formed,results in the long-term storage of carbon in woody biomass,making it a key component of the global carbon cycle.Understanding how environmental drivers influence xylogenesis during the growing season is therefore of great interest.However,studying shortterm drivers of wood production using xylogenetic data is complicated by the usual sampling scheme and the influence of eccentric growth,i.e.,heterogeneous growth around the stem.In this study,we improve xylogenesis research by introducing a statistical approach that explicitly considers seasonal phenology,short-term growth rates,and growth eccentricity.To this end,we developed Bayesian models of xylogenesis and compared them with a conventional method based on the use of Gompertz functions.Our results show that eccentricity generated high temporal autocorrelation between successive samples,and that explicitly taking it into account improved both the representativeness of phenology and intra-ring variability.We observed consistent short-term patterns in the model residuals,suggesting the influence of an unaccounted-for environmental variable on cell production.The proposed models offer several advantages over traditional methods,including robust confidence intervals around predictions,consistency with phenology,and reduced sensitivity to extreme observations at the end of the growing season,often linked to eccentric growth.These models also provide a benchmark for mechanistic testing of short-term drivers of wood formation.展开更多
Industrial robot dynamics lay the foundation for high-precision and high-speed control, and accurate identification of dynamic parameters is essential for precise dynamic calculations. The choice of friction models is...Industrial robot dynamics lay the foundation for high-precision and high-speed control, and accurate identification of dynamic parameters is essential for precise dynamic calculations. The choice of friction models is a critical component in the identification of industrial robot dynamics. Traditional static friction models struggle to capture the hysteresis effects caused by robot joint elasticity and clearances, leading to large torque prediction errors when the joint velocity crosses zero. Due to the presence of hysteresis effects, the joint velocity crosses zero in the forward direction, and the reverse direction will have different friction patterns. Although the hysteresis effects can be modeled as an ordinary differential equation(ODE), it is difficult to determine the ODE structure that achieves both generalization and accuracy to describe the hysteresis effects of the friction model. To address this issue, we propose the neural hysteresis friction(NHF), which uses neural ODE to model the hysteresis effects in a data-driven manner, thereby mitigating the current inadequacies in the study of dynamic friction characteristics. The experiments on a real 6-axis industrial robot demonstrate that our proposed method can accurately model the friction dynamics during directional switching and outperform other modeling methods. Velocity tracking control experiments show that NHF can effectively reduce tracking errors when the velocity crosses zero.展开更多
Mount Kandil is situated in the eastern sector of the EAHP(Eastern Anatolian High Plateau),to the south of the Lesser Caucasus.The mountain lies at the westernmost end of the Aras Mountains,which extends approximately...Mount Kandil is situated in the eastern sector of the EAHP(Eastern Anatolian High Plateau),to the south of the Lesser Caucasus.The mountain lies at the westernmost end of the Aras Mountains,which extends approximately 80 km along a NW-SE axis.With a summit reaching~3214 m(a.s.l.),Mount Kandil is a stratovolcano that,like many other peaks within the EAHP and the Lesser Caucasus,experienced significant environmental changes during Late Pleistocene.Among these,glacial processes stand out as the most profound,having distinctly shaped the mountains geomorphic landscape.This study presents,for the first time,a comprehensive analysis of the glacial morphology of Mount Kandil based on multiple datasets.Field-based morphological observations indicate that an area of approximately 32.62 km^(2)has been sculpted by glacial activity.Within six glaciated regions on Mount Kandil,25 cirques and 6 glacial valleys have been identified.In addition,moraines in various locations exhibit characteristic morphologies.Furthermore,valley glaciers are inferred to have descended to altitudes as low as~2000 m.The paleoequilibrium line(p ELA)was estimated to use AABR method within GIS,yielding a mean pELA of~2730 m.Ice thickness modelling indicates that the thickness of glaciers in the Kandil Mountain valleys reaches up to~350 m.Due to its orographic extension,Mount Kandil is exposed to humid northwest winds and receives substantial frontal precipitation(~686 mm annually).The compiled geomorphic,cartographic and morphometric parameters suggest that the glaciation dynamics of Mount Kandil—situated between the Southeastern Taurus and the Lesser Caucasus—closely resemble those observed in the Lesser Caucasus.This indicates that glaciation was primarily governed by northern atmospheric systems with additional influences from southerly or westerly winds.The integrated data also underscores the role of multiple atmospheric systems in controlling the glaciation regime around the Lesser Caucasus.Additionally,findings on regional pELA question the common belief that pELA increases eastward in EAHP.展开更多
In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicti...In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicting long-term roadway stability,necessitating the development of a reliable constitutive creep model and numerical simulation approach.In this study,creep experiments were conducted on pre-damaged rock with varying initial damage levels to investigate the time-dependent mechanical properties.Based on the experimental results,an accelerated-creep criterion was proposed,and an elastic-viscoplastic creep damage model(EVPCD)was established that simultaneously considers the effects of time-dependent damage and instantaneous damage caused by stress disturbances on rock creep behavior.Subsequently,the effectiveness of the proposed creep model was verified using experimental data,and the secondary development of the EVPCD model was completed based on the FLAC3D platform.Following this,a long-term stability analysis method of deep surrounding rock that accounts for excavation-and mining-induced disturbances was proposed.Using the main roadway of Xutuan Coal Mine as a case study,numerical simulations were carried out to investigate the time-dependent deformation and failure characteristics of the surrounding rock following excavation and mining disturbance.Combined with on-site monitoring of the surrounding rock damage areas,the results indicate that the EVPCD outperforms the CVISC and Nishihara models in predicting the time-dependent behavior of deep surrounding rock.展开更多
Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and ...Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and finite element analyses.Three tests of the slope deformation under earthquake and post-rainfall earthquakes are first studied using image analysis techniques.Then,based on an elastoplastic constitutive model,numerical simulations are carried out using the finite element method and compared with the centrifuge test results.Finally,a parametric study is performed to clarify the effects of antecedent rainfall on earthquake-induced slope deformation.The results show that slope deformation caused by post-rainfall earthquakes differs from that caused by earthquakes without antecedent rainfall.The seepage flow and soil strength of the slope are affected by previous rainfall conditions,such as intensity and duration,which directly influence the slope deformation caused by the subsequent earthquake.Soil displacement and strain become greater and the slip surface is more noticeable during the post-rainfall earthquake of higher intensity.In addition,the time interval between the rainfall and the earthquake has a considerable impact on the detailed characteristics of the slope deformation,and the significant deformation occurs at the time of lowest soil strength when seepage flow reaches the lower part of the slope.Moreover,the repeated intermittent rainfall greatly affects the subsequent earthquake-induced slope deformation,the main characteristics of which are closely related to the changes in saturation and strength of the slope.However,with the prolonged time gap between each round of rainfall,the earthquake-induced slope deformation becomes insignificant.展开更多
Hepatitis B Virus(HBV)infection and heavy alcohol consumption are the two primary pathogenic causes of liver cirrhosis.In this paper,we proposed a deterministic mathematical model and a logistic equation to investigat...Hepatitis B Virus(HBV)infection and heavy alcohol consumption are the two primary pathogenic causes of liver cirrhosis.In this paper,we proposed a deterministic mathematical model and a logistic equation to investigate the dynamics of liver cirrhosis progression as well as to explain the implications of variations in alcohol consumption on chronic hepatitis B patients,respectively.The intricate interactions between liver cirrhosis,recovery,and treatment dynamics are captured by the model.This study aims to show that alcohol consumption by Hepatitis B-infected individuals accelerates liver cirrhosis progression while treatment of acutely infected individuals reduces it.We proved that a unique solution of the proposed model exists,which is positive and bounded.Using the next-generation matrix approach,two basic reproductive numbers R_(A_(0))and R_(A_(max))are calculated to identify future recurrence.The equilibrium points are calculated,and both equilibria are proved locally and globally asymptotically stable when R_(0)is below and above one,respectively.It is shown that bifurcation exists at R_(0)=1 and a detailed proof for forward bifurcation is given.Furthermore,we performed the sensitivity analysis of the model parameters on R_(0).For the confirmation of analytical work,we performed numerical simulations,and the results indicate that the treatment and the inhibitory effects reduce the risk of developing liver cirrhosis in individuals,while heavy alcohol consumption accelerates markedly the liver cirrhosis progression in patients with chronic hepatitis B.展开更多
Accurate fault modeling is essential for understanding earthquake rupture processes and improving seismic hazard assessment.We present a unified framework that integrates geodetic data with multidisciplinary constrain...Accurate fault modeling is essential for understanding earthquake rupture processes and improving seismic hazard assessment.We present a unified framework that integrates geodetic data with multidisciplinary constraints,including relocated aftershocks,geological observations,and geophysical information,to adaptively model fault geometry and slip in diverse scenarios such as multi-segment and multi-event ruptures.The framework combines adaptive fault construction with a scenario-driven Bayesian joint inversion approach.Fault geometries are first built from prior constraints,such as surface ruptures and aftershocks,and then refined through probabilistic inference when such data are inadequate.To enhance computational efficiency,we introduce a Sequential Monte Carlo Fukuda-Johnson(SMC-FJ)strategy.This separates nonlinear parameters-including geometry,data weights,and smoothing factors-from linear slip parameters,which are conditionally solved via constrained least squares.Geometry updates follow a hierarchical adjustment scheme based on point,line,and structural units,enabling flexibility across rupture complexities.Synthetic tests and four case studies,including the 2022 Menyuan,2023 Türkiye,2022 Luding,and 2019 Ridgecrest earthquakes,demonstrate robustness under various constraints.For the Menyuan earthquake,full Bayesian inversion confirms that the fault geometry constrained by relocated aftershocks is highly accurate and needs only minor adjustment to match the observed surface deformation.The results further show that gradual changes in fault strike and dip modulated rupture arrest and postseismic stress accumulation,highlighting the critical role of inherited geometric structure in controlling rupture termination and delayed seismic activation.展开更多
Under sustained strong stochastic impact loads,floating-supported friction plates are susceptible to the formation of fatigue cracks that propagate along the rim.The nonlinearity and randomness introduced by the crack...Under sustained strong stochastic impact loads,floating-supported friction plates are susceptible to the formation of fatigue cracks that propagate along the rim.The nonlinearity and randomness introduced by the cracked teeth participating in the impacts significantly influence the service life and reliability of the transmission system.In this paper,an improved stiffness excitation modeling method is developed for friction plate teeth with rim cracks.It overcomes the limitations of traditional approaches that fail to accurately assess the narrow-band,large-diameter friction plate teeth with rim cracks due to constraints imposed by boundary conditions.Then,an original dynamic impact model for the floating-supported friction plate and inner hub system is proposed,incorporating the effects of bending-torsional-axial-tilting coupled motions on tooth mesh excitations and dynamic responses.This model addresses the limitations of conventional models that only consider bending-torsion coupling,thereby providing a more comprehensive representation of the system's multi-dimensional dynamic behavior.The effects of the crack propagation depth and the number of cracked teeth on the stochastic impact characteristics and vibration responses of the system are investigated.Furthermore,finite element simulations and experimental tests are conducted to validate the cracked tooth stiffness excitations and dynamic impact responses,respectively.The proposed model is anticipated to provide both a theoretical foundation and practical guidance for fault diagnosis and reliability assessment of clutch friction plates.展开更多
基金Supported by the National Natural Science Foundation of China(60801050)the Excellent Talent Fund of Beijing(2011)Excellent Young Scholars Research Fund of Beijing Institute ofTechnology(2012)
文摘The finite-element modeling and simulations of the intra-body communication （IBC） were investigated to provide a theoretical basis for biomedical monitoring. A finite-element model for the whole human body was developed to simulate the IBC. The simulation of galvanic coupling IBC and electrostatic coupling IBC were implemented along with different signal transmission paths, and their attenuations were calculated. Our study showed that the position near the signal electrode had higher potential than other positions in the two types of IBC, while the potential generally decreased along the axis of the body parts. Both signal attenuations of the two types IBC increased with increasing signal transmission distance, and the electrostatic coupling IBC had comparatively higher receiving potential than the galvanic coupling IBC. The results indicated that the proposed modeling method could be used for the research of biomedical monitoring based on IBC technology.
文摘Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and with the nucleus removed were developed and validatedusing published experimental and clinical data. Then the models with a stainless steel nucleusprosthesis implanted and with polymer nucleus prostheses of various properties implanted were usedfor the 3D finite-element biomechanical analysis. All the above simulation and analysis were carriedout for the L4/L5 disc under a human worst--daily compression load of 2000 N. The results show thatthe polymer materials with Young's modulus of elasticity E = 0.1-100 MPa and Poisson's ratio v=0.35-0.5 are suitable to produce artificial nucleus prosthesis in view of biomechanicalconsideration.
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
基金supported by the National Natural Science Foundation of China(No.41104068)the Deep Exploration in China,Sino Probe-03-05
文摘Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor CSAMT forward modeling using the vector finite-element method. To verify the feasibility of the algorithm, we calculate the electric field, magnetic field, and tensor impedance of the 3D CSAMT far-zone field in layered media and compare them with theoretical solutions. In addition, a three-dimensional anomaly in half-space is also simulated, and the response characteristics of the impedance tensor and the apparent resistivity and impedance phase are analyzed. The results suggest that the vector finite-element method produces high-precision electromagnetic field and impedance tensor data, satisfies the electric field discontinuity, and does not require divergence correction using the vector finite-element method.
基金financially supported by the National Natural Science Foundation of China(No.41204055,41164003,and 41104074)Opening Project(No.SMIL-2014-06) of Hubei Subsurface Multi-scale Imaging Lab(SMIL),China University of Geosciences(Wuhan)
文摘The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.
基金National High Technology Research and Development Program(863 Program)(No.2006AA06Z105,2007AA06Z134)
文摘Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.
基金the World Climate Research Programme(WCRP),Climate Variability and Predictability(CLIVAR),and Global Energy and Water Exchanges(GEWEX)for facilitating the coordination of African monsoon researchsupport from the Center for Earth System Modeling,Analysis,and Data at the Pennsylvania State Universitythe support of the Office of Science of the U.S.Department of Energy Biological and Environmental Research as part of the Regional&Global Model Analysis(RGMA)program area。
文摘In recent years,there has been an increasing need for climate information across diverse sectors of society.This demand has arisen from the necessity to adapt to and mitigate the impacts of climate variability and change.Likewise,this period has seen a significant increase in our understanding of the physical processes and mechanisms that drive precipitation and its variability across different regions of Africa.By leveraging a large volume of climate model outputs,numerous studies have investigated the model representation of African precipitation as well as underlying physical processes.These studies have assessed whether the physical processes are well depicted and whether the models are fit for informing mitigation and adaptation strategies.This paper provides a review of the progress in precipitation simulation overAfrica in state-of-the-science climate models and discusses the major issues and challenges that remain.
基金supported in part by the National Key Research and Development Program of China(2021YFB2900501)in part by the Shaanxi Science and Technology Innovation Team(2023-CX-TD-03)+3 种基金in part by the Science and Technology Program of Shaanxi Province(2021GXLH-Z-038)in part by the Natural Science Foundation of Hunan Province(2023JJ40607 and 2023JJ50045)in part by the Scientific Research Foundation of Hunan Provincial Education Department(23B0713 and 24B0603)in part by the National Natural Science Foundation of China(62401371,62101275,and 62372070).
文摘(Quasi-)closed-form results for the statistical properties of unmanned aerial vehicle(UAV)airto-ground channels are derived for the first time using a novel spatial-vector-based method from a threedimensional(3-D)arbitrary-elevation one-cylinder model.The derived results include a closed-form expression for the space-time correlation function and some quasi-closed-form ones for the space-Doppler power spectrum density,the level crossing rate,and the average fading duration,which are shown to be the generalizations of those previously obtained from the two-dimensional(2-D)one-ring model and the 3-D low-elevation one-cylinder model for terrestrial mobile-to-mobile channels.The close agreements between the theoretical results and the simulations as well as the measurements validate the utility of the derived channel statistics.Based on the derived expressions,the impacts of some parameters on the channel characteristics are investigated in an effective,efficient,and explicable way,which leads to a general guideline on the manual parameter estimation from the measurement description.
基金support from the National Natural Science Foundation of China(Grant Nos.42277161 and 42230709).
文摘In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical properties of rocks,the cracking processes of pre-cracked rocks have been extensively studied using numerical modeling methods.The peridynamics(PD)exhibits advantages over other numerical methods due to the absence of the requirements for remeshing and external crack growth criterion.However,for modeling pre-cracked rock cracking processes under impact,current PD implementations lack generally applicable rock constitutive models and impact contact models,which leads to difficulties in determining rock material parameters and efficiently calculating impact loads.This paper proposes a non-ordinary state-based peridynamics(NOSBPD)modeling method integrating the Drucker-Prager(DP)plasticity model and an efficient contact model to address the above problems.In the proposed method,the Drucker-Prager plasticity model is integrated into the NOSBPD,thereby equipping NOSBPD with the capability to accurately characterize the nonlinear stress-strain relationship inherent in rocks.An efficient contact model between particles and meshes is designed to calculate the impact loads,which is essentially a coupling method of PD with the finite element method(FEM).The effectiveness of the proposed NOSBPD modeling method is verified by comparison with other numerical methods and experiments.Experimental results indicate that the proposed method can effectively and accurately predict the 3D cracking processes of pre-cracked cracks under impact loading,and the maximum principal stress is the key driver behind wing crack formation in pre-cracked rocks.
基金the funding support from the Smart Medicine and Engineering Interdisciplinary Innovation Project of Ningbo University(No.ZHYG003)。
文摘Glassy polymers are widely used in biomedical applications in a solvent environment,yet their long-term performance is governed by the competing effects of physical aging and solvent-induced plasticization.Here,we develop a constitutive model that explicitly couples the solvent concentration,structural relaxation,and mechanical response.This framework is built on a multiplicative decomposition of deformation and an Eyring-type flow rule,with structural evolution described by an effective temperature.A generalized shift factor is introduced to quantify how the solvent concentration and effective temperature jointly affect the relaxation time,thereby integrating physical aging and plasticization.The model is subsequently applied to methacrylate(MA)-based copolymer networks immersed in phosphate-buffered saline for up to nine months.Simulations accurately capture key experimental features,including the strong softening of highly swellable networks,the partial recovery due to aging,and the mitigating role of hydrophobic crosslinking in reducing solvent uptake.While the current single-mode description cannot reproduce the full relaxation spectrum,it establishes an efficient framework for predicting the long-term mechanical performance under coupled environmental and mechanical loading.This study provides a constitutive description of solvent-swollen glassy polymers,offering mechanistic insight into the interplay between plasticization and aging.Beyond biomedical MA networks,this framework establishes a foundation for predicting the long-term performance of polymer glasses under coupled aqueous environmental and mechanical loading.
基金output of a research project implemented as part of the Basic Research Program at HSE University。
文摘Activation pruning reduces neural network complexity by eliminating low-importance neuron activations,yet identifying the critical pruning threshold—beyond which accuracy rapidly deteriorates—remains computationally expensive and typically requires exhaustive search.We introduce a thermodynamics-inspired framework that treats activation distributions as energy-filtered physical systems and employs the free energy of activations as a principled evaluation metric.Phase-transition-like phenomena in the free-energy profile—such as extrema,inflection points,and curvature changes—yield reliable estimates of the critical pruning threshold,providing a theoretically grounded means of predicting sharp accuracy degradation.To further enhance efficiency,we propose a renormalized free energy technique that approximates full-evaluation free energy using only the activation distribution of the unpruned network.This eliminates repeated forward passes,dramatically reducing computational overhead and achieving speedups of up to 550×for MLPs.Extensive experiments across diverse vision architectures(MLP,CNN,ResNet,MobileNet,Vision Transformer)and text models(LSTM,BERT,ELECTRA,T5,GPT-2)on multiple datasets validate the generality,robustness,and computational efficiency of our approach.Overall,this work establishes a theoretically grounded and practically effective framework for activation pruning,bridging the gap between analytical understanding and efficient deployment of sparse neural networks.
基金supported by the National Natural Science Foundation of China(Grant No.12272018)the National Key Basic Research Project(2022JCJQZD20600).
文摘Kinetic impact is the most practical planetary-defense technique,with momentum-transfer efficiency central to deflection design.We present a Monte Carlo photometric framework that couples ejecta sampling,dynamical evolution,and image synthesis to compare directly with HST,LICIACube,ground-based and Lucy observations of the DART impact.Decomposing ejecta into(1)a highvelocity(~1600 m/s)plume exhibiting Na/K resonance,(2)a low-velocity(~1 m/s)conical component shaped by binary gravity and solar radiation pressure,and(3)meter-scale boulders,we quantify each component’s mass and momentum.Fitting photometric decay curves and morphological evolution yields size-velocity distributions and,via scaling laws,estimates of Dimorphos’bulk density,cratering parameters,and cohesive strength that agree with dynamical constraints.Photometric ejecta modeling therefore provides a robust route to constrain momentum enhancement and target properties,improving predictive capability for kinetic-deflection missions.
基金supported by the Discovery Grants program of the Natural Sciences and Engineering Research Council of Canada(No.RGPIN-2021-03553)the Canadian Research Chair in dendroecology and dendroclimatology(CRC-2021-00368)+3 种基金the Ministère des Ressources Naturelles et des Forèts(MRNF,Contract no.142332177-D)the Natural Sciences and Engineering Research Council of Canada(Alliance Grant No.ALLRP 557148-20,obtained in partnership with the MRNF and Resolute Forest Products)the Fonds de recherche du Qu ebec–Nature et technologies(Partnership Research Program on the Contribution of the Forestry Sector to Climate Change MitigationGrant No.2022-0FC-309064)。
文摘Xylogenesis,the process through which wood cells are formed,results in the long-term storage of carbon in woody biomass,making it a key component of the global carbon cycle.Understanding how environmental drivers influence xylogenesis during the growing season is therefore of great interest.However,studying shortterm drivers of wood production using xylogenetic data is complicated by the usual sampling scheme and the influence of eccentric growth,i.e.,heterogeneous growth around the stem.In this study,we improve xylogenesis research by introducing a statistical approach that explicitly considers seasonal phenology,short-term growth rates,and growth eccentricity.To this end,we developed Bayesian models of xylogenesis and compared them with a conventional method based on the use of Gompertz functions.Our results show that eccentricity generated high temporal autocorrelation between successive samples,and that explicitly taking it into account improved both the representativeness of phenology and intra-ring variability.We observed consistent short-term patterns in the model residuals,suggesting the influence of an unaccounted-for environmental variable on cell production.The proposed models offer several advantages over traditional methods,including robust confidence intervals around predictions,consistency with phenology,and reduced sensitivity to extreme observations at the end of the growing season,often linked to eccentric growth.These models also provide a benchmark for mechanistic testing of short-term drivers of wood formation.
基金supported by the National Natural Science Foundation of China (Grant No.52188102)。
文摘Industrial robot dynamics lay the foundation for high-precision and high-speed control, and accurate identification of dynamic parameters is essential for precise dynamic calculations. The choice of friction models is a critical component in the identification of industrial robot dynamics. Traditional static friction models struggle to capture the hysteresis effects caused by robot joint elasticity and clearances, leading to large torque prediction errors when the joint velocity crosses zero. Due to the presence of hysteresis effects, the joint velocity crosses zero in the forward direction, and the reverse direction will have different friction patterns. Although the hysteresis effects can be modeled as an ordinary differential equation(ODE), it is difficult to determine the ODE structure that achieves both generalization and accuracy to describe the hysteresis effects of the friction model. To address this issue, we propose the neural hysteresis friction(NHF), which uses neural ODE to model the hysteresis effects in a data-driven manner, thereby mitigating the current inadequacies in the study of dynamic friction characteristics. The experiments on a real 6-axis industrial robot demonstrate that our proposed method can accurately model the friction dynamics during directional switching and outperform other modeling methods. Velocity tracking control experiments show that NHF can effectively reduce tracking errors when the velocity crosses zero.
基金supported by Van Yüzüncü Yıl University,Scientific Research Projects Coordination Unit(Project No:SDK-2025-11935)Van Yüzüncü Yıl University,Scientific Research Projects Coordination Unit for supporting the study。
文摘Mount Kandil is situated in the eastern sector of the EAHP(Eastern Anatolian High Plateau),to the south of the Lesser Caucasus.The mountain lies at the westernmost end of the Aras Mountains,which extends approximately 80 km along a NW-SE axis.With a summit reaching~3214 m(a.s.l.),Mount Kandil is a stratovolcano that,like many other peaks within the EAHP and the Lesser Caucasus,experienced significant environmental changes during Late Pleistocene.Among these,glacial processes stand out as the most profound,having distinctly shaped the mountains geomorphic landscape.This study presents,for the first time,a comprehensive analysis of the glacial morphology of Mount Kandil based on multiple datasets.Field-based morphological observations indicate that an area of approximately 32.62 km^(2)has been sculpted by glacial activity.Within six glaciated regions on Mount Kandil,25 cirques and 6 glacial valleys have been identified.In addition,moraines in various locations exhibit characteristic morphologies.Furthermore,valley glaciers are inferred to have descended to altitudes as low as~2000 m.The paleoequilibrium line(p ELA)was estimated to use AABR method within GIS,yielding a mean pELA of~2730 m.Ice thickness modelling indicates that the thickness of glaciers in the Kandil Mountain valleys reaches up to~350 m.Due to its orographic extension,Mount Kandil is exposed to humid northwest winds and receives substantial frontal precipitation(~686 mm annually).The compiled geomorphic,cartographic and morphometric parameters suggest that the glaciation dynamics of Mount Kandil—situated between the Southeastern Taurus and the Lesser Caucasus—closely resemble those observed in the Lesser Caucasus.This indicates that glaciation was primarily governed by northern atmospheric systems with additional influences from southerly or westerly winds.The integrated data also underscores the role of multiple atmospheric systems in controlling the glaciation regime around the Lesser Caucasus.Additionally,findings on regional pELA question the common belief that pELA increases eastward in EAHP.
基金funded by the National Natural Science Foundation of China(Nos.52004098,U24B2041,and 52274079)the Key Research and Development Program of Henan Province(No.251111320400)+1 种基金the Key Research Project Plan for Higher Education Institutions in Henan Province(Nos.24A570006 and 25A570002)the Scientific and Technological Research Project in Henan Province(No.242102320061).
文摘In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicting long-term roadway stability,necessitating the development of a reliable constitutive creep model and numerical simulation approach.In this study,creep experiments were conducted on pre-damaged rock with varying initial damage levels to investigate the time-dependent mechanical properties.Based on the experimental results,an accelerated-creep criterion was proposed,and an elastic-viscoplastic creep damage model(EVPCD)was established that simultaneously considers the effects of time-dependent damage and instantaneous damage caused by stress disturbances on rock creep behavior.Subsequently,the effectiveness of the proposed creep model was verified using experimental data,and the secondary development of the EVPCD model was completed based on the FLAC3D platform.Following this,a long-term stability analysis method of deep surrounding rock that accounts for excavation-and mining-induced disturbances was proposed.Using the main roadway of Xutuan Coal Mine as a case study,numerical simulations were carried out to investigate the time-dependent deformation and failure characteristics of the surrounding rock following excavation and mining disturbance.Combined with on-site monitoring of the surrounding rock damage areas,the results indicate that the EVPCD outperforms the CVISC and Nishihara models in predicting the time-dependent behavior of deep surrounding rock.
基金supported by the China Postdoctoral Science Foundation(CPSF)(Grant No.2024M762769)the Natural Science Basic Research Program of Shaanxi(Grant No.2024JC-YBQN-0333)the Postdoctoral Fellowship Program of CPSF(Grant No.GZC20232230).
文摘Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and finite element analyses.Three tests of the slope deformation under earthquake and post-rainfall earthquakes are first studied using image analysis techniques.Then,based on an elastoplastic constitutive model,numerical simulations are carried out using the finite element method and compared with the centrifuge test results.Finally,a parametric study is performed to clarify the effects of antecedent rainfall on earthquake-induced slope deformation.The results show that slope deformation caused by post-rainfall earthquakes differs from that caused by earthquakes without antecedent rainfall.The seepage flow and soil strength of the slope are affected by previous rainfall conditions,such as intensity and duration,which directly influence the slope deformation caused by the subsequent earthquake.Soil displacement and strain become greater and the slip surface is more noticeable during the post-rainfall earthquake of higher intensity.In addition,the time interval between the rainfall and the earthquake has a considerable impact on the detailed characteristics of the slope deformation,and the significant deformation occurs at the time of lowest soil strength when seepage flow reaches the lower part of the slope.Moreover,the repeated intermittent rainfall greatly affects the subsequent earthquake-induced slope deformation,the main characteristics of which are closely related to the changes in saturation and strength of the slope.However,with the prolonged time gap between each round of rainfall,the earthquake-induced slope deformation becomes insignificant.
文摘Hepatitis B Virus(HBV)infection and heavy alcohol consumption are the two primary pathogenic causes of liver cirrhosis.In this paper,we proposed a deterministic mathematical model and a logistic equation to investigate the dynamics of liver cirrhosis progression as well as to explain the implications of variations in alcohol consumption on chronic hepatitis B patients,respectively.The intricate interactions between liver cirrhosis,recovery,and treatment dynamics are captured by the model.This study aims to show that alcohol consumption by Hepatitis B-infected individuals accelerates liver cirrhosis progression while treatment of acutely infected individuals reduces it.We proved that a unique solution of the proposed model exists,which is positive and bounded.Using the next-generation matrix approach,two basic reproductive numbers R_(A_(0))and R_(A_(max))are calculated to identify future recurrence.The equilibrium points are calculated,and both equilibria are proved locally and globally asymptotically stable when R_(0)is below and above one,respectively.It is shown that bifurcation exists at R_(0)=1 and a detailed proof for forward bifurcation is given.Furthermore,we performed the sensitivity analysis of the model parameters on R_(0).For the confirmation of analytical work,we performed numerical simulations,and the results indicate that the treatment and the inhibitory effects reduce the risk of developing liver cirrhosis in individuals,while heavy alcohol consumption accelerates markedly the liver cirrhosis progression in patients with chronic hepatitis B.
基金supported by the National Natural Science Foundation of China(Grant Nos.42130101,42474002,42374003&42564002)the Jiangxi Provincial Natural Science Foundation(Grant No.20252BAC240262).
文摘Accurate fault modeling is essential for understanding earthquake rupture processes and improving seismic hazard assessment.We present a unified framework that integrates geodetic data with multidisciplinary constraints,including relocated aftershocks,geological observations,and geophysical information,to adaptively model fault geometry and slip in diverse scenarios such as multi-segment and multi-event ruptures.The framework combines adaptive fault construction with a scenario-driven Bayesian joint inversion approach.Fault geometries are first built from prior constraints,such as surface ruptures and aftershocks,and then refined through probabilistic inference when such data are inadequate.To enhance computational efficiency,we introduce a Sequential Monte Carlo Fukuda-Johnson(SMC-FJ)strategy.This separates nonlinear parameters-including geometry,data weights,and smoothing factors-from linear slip parameters,which are conditionally solved via constrained least squares.Geometry updates follow a hierarchical adjustment scheme based on point,line,and structural units,enabling flexibility across rupture complexities.Synthetic tests and four case studies,including the 2022 Menyuan,2023 Türkiye,2022 Luding,and 2019 Ridgecrest earthquakes,demonstrate robustness under various constraints.For the Menyuan earthquake,full Bayesian inversion confirms that the fault geometry constrained by relocated aftershocks is highly accurate and needs only minor adjustment to match the observed surface deformation.The results further show that gradual changes in fault strike and dip modulated rupture arrest and postseismic stress accumulation,highlighting the critical role of inherited geometric structure in controlling rupture termination and delayed seismic activation.
基金supported by the National Natural Science Foundation of China(Grant Nos.52505101,52475087,52475089,52365010)the Early-Career Young Scientists and Technologists Project of Jiangxi Province(Grant No.20252BEJ730175)。
文摘Under sustained strong stochastic impact loads,floating-supported friction plates are susceptible to the formation of fatigue cracks that propagate along the rim.The nonlinearity and randomness introduced by the cracked teeth participating in the impacts significantly influence the service life and reliability of the transmission system.In this paper,an improved stiffness excitation modeling method is developed for friction plate teeth with rim cracks.It overcomes the limitations of traditional approaches that fail to accurately assess the narrow-band,large-diameter friction plate teeth with rim cracks due to constraints imposed by boundary conditions.Then,an original dynamic impact model for the floating-supported friction plate and inner hub system is proposed,incorporating the effects of bending-torsional-axial-tilting coupled motions on tooth mesh excitations and dynamic responses.This model addresses the limitations of conventional models that only consider bending-torsion coupling,thereby providing a more comprehensive representation of the system's multi-dimensional dynamic behavior.The effects of the crack propagation depth and the number of cracked teeth on the stochastic impact characteristics and vibration responses of the system are investigated.Furthermore,finite element simulations and experimental tests are conducted to validate the cracked tooth stiffness excitations and dynamic impact responses,respectively.The proposed model is anticipated to provide both a theoretical foundation and practical guidance for fault diagnosis and reliability assessment of clutch friction plates.
