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 1739 M8.0 Pingluo earthquake occurred around the Yinchuan Graben,bounded by the Helan Mountains to the west and the Ordos Block to the east.Seismological observations have shown that surface fault displacement rea...The 1739 M8.0 Pingluo earthquake occurred around the Yinchuan Graben,bounded by the Helan Mountains to the west and the Ordos Block to the east.Seismological observations have shown that surface fault displacement reaches about 2–3 m,mainly by dip-slip motion along the Helanshan Piedmont Fault.However,the documented seismic intensity is distributed predominantly within the basin area,exhibiting a sharp asymmetry across the Helanshan Piedmont Fault.Thus,the general pattern of earthquake faulting is still under debate.We built a three-dimensional elastodynamic finiteelement model to reappraise the fault mechanism.In the model,predictions from synthetic rupture models,based on available observations and the earthquake scaling law,were used as an input with the split-node technique,and the effect of basin sediment on elastic wave propagation was considered.The numerical results show that if an earthquake occurred on the Helanshan Piedmont Fault characterized by a high-angle(70°)normal fault,earthquake shaking,as predicted from the modeled peak ground velocity and peak ground acceleration,has difficulty fitting the observed result,even when the effect of sediment amplification is considered.To better fit the observed shaking pattern,the dip angle of the Helanshan Piedmont Fault must be less than about 35°between the depths of about 8–27 km,where the coseismic slip may reach about 6 m.This result leads us to conclude that the 1739 M8.0 great earthquake likely occurred on a listric normal fault at depth,in agreement with the geometry of the Helanshan Piedmont Fault,as recently evidenced by seismic reflection explorations.This conclusion means that in an intracontinental setting,a reduction in the fault dip angle along the subsurface could increase the width of the fault in the elastic crust,making misalignment between the surface rupture and the isoseismals and resulting in an increase in the upper bound of earthquake magnitude relative to simple high-angle faulting.展开更多
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.展开更多
The dissolution of MgO-refractory into the slag had an obvious influence on the steel-slag reaction and the slag property,especially for high-aluminum steels.The dissolution behavior of MgO-refractory was investigated...The dissolution of MgO-refractory into the slag had an obvious influence on the steel-slag reaction and the slag property,especially for high-aluminum steels.The dissolution behavior of MgO-refractory was investigated under various conditions,including the temperature,the initial steel composition,and the initial slag composition.A steel-slag-refractory kinetic model for high-aluminum steel was developed,which incorporated the process of MgO-refractory dissolution.The dependence of the MgO mass transfer coefficient k_(MgO)^(r)on temperature T during MgO-refractory dissolution process was established,as described by ln k_(MgO)^(r)=63,754/T+24.38524.It was indicated that the MgO dissolution rate was significantly influenced by the temperature.A higher temperature increased the dissolution rate of MgO.The initial steel composition had a slight impact on the MgO dissolution rate.Additionally,the initial slag composition strongly impacted the MgO saturation concentration and the dissolution rate.A lower initial Al_(2)O_(3)/SiO_(2)ratio increased the MgO dissolution rate.The steel-slag-refractory kinetic model accurately predicted the dissolution of MgO-refractory and the influence of dissolved MgO on the viscosity and composition change during steel-slag-refractory reactions.It was suggested that a higher temperature can hardly reduce the viscosity due to the dissolution of the MgO-refractory.展开更多
In materials science and engineering design,high-fidelity and high-efficiency numerical simulation has become a driving force for innovation and practical implementation.To address longstanding bottlenecks in the deve...In materials science and engineering design,high-fidelity and high-efficiency numerical simulation has become a driving force for innovation and practical implementation.To address longstanding bottlenecks in the development of conventional material constitutive models—such as lengthy modeling cycles and difficulties in numerical implementation—this study proposes an intelligent modeling and code generation approach powered by large languagemodels.A structured knowledge base integrating constitutive theory,numerical algorithms,and UMAT(User Material)interface specifications is constructed,and a retrieval-augmented generation strategy is employed to establish an end-to-end workflow spanning experimental data parsing,constitutive model formulation,and automatic UMAT subroutine generation.Experimental results show that the method achieves high accuracy for both a classical Johnson–Cookmodel and a physics-informed neural network(PINN)model,with key parameter identification errors below 5%.Moreover,the automatically generated UMAT subroutines yield finite element simulation results in Abaqus that are highly consistent with theoretical predictions(coefficient of determination R2>0.98)while maintaining good numerical stability.This framework is currently focused on the automatic construction of rate-dependent elastoplastic material models,and its core method also provides a clear path for extending to other constitutive categories such as hyperelasticity and viscoelasticity.This work provides an effective technical route for the rapid development and reliable numerical implementation of material constitutive models,significantly advancing the intelligence level of computational mechanics research and improving engineering application efficiency.展开更多
(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.展开更多
The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we pro...The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we propose a new approach by combining centrifugal analog modeling with numerical simulation to simulate the tectonic uplift history of the plateau based on the lower crustal flow model,and to investigate the material migration characteristics and the influence of crustal motion velocity and ductile layer viscosity on the plateau tectonic geomorphology.The models reproduce steep-sided flat-topped geomorphic features and clockwise rotation of the material at eastern Himalayan Syntaxis,verifying the rationality of the models.The results show that the greater the crustal motion velocity and the greater the ductile layer viscosity,the steeper the terrain change;and conversely,the smaller the crustal motion velocity and the smaller the ductile layer viscosity,the gentler the terrain change.This study further indicates that the weak lower crust plays an important role in the formation of geomorphic features and material migration characteristics of Qinghai-Tibet Plateau,and provides a new insight for the study of the uplift mechanism of the Tibetan Plateau.展开更多
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.展开更多
In-space cable-driven manipulators exhibit several advantages,such as a large range of motion,high dexterity,and lightweight structure.However,kinematic and dynamic analysis play an essential role in designing a cable...In-space cable-driven manipulators exhibit several advantages,such as a large range of motion,high dexterity,and lightweight structure.However,kinematic and dynamic analysis play an essential role in designing a cable-driven manipulator.In this paper,the kinematic analysis of a type of cable-driven manipulator is performed,and a motion planning scheme is conducted to actuate this manipulator.Moreover,a flexible multi-body dynamic model of a cable-driven manipulator considering the frictional contact between the cables and pulleys is established.To describe properties such as flexibility,vibration,and variable length of the cable,this paper utilizes reducedorder beam elements of the Absolute Nodal Coordinates Formulation(ANCF)in Arbitrary Lagrangian Eulerian(ALE)framework.Additionally,a virtual element is introduced to model the contact segment in the cable-pulley system.A tension decay factor is employed to account for the friction in the contact segment.To validate the proposed method,a semi-analytical model based on D'Alembert's principle is established.Cross-verification is performed to validate the accuracy of both models.The model is further applied to simulate the rotation of the cable-driven manipulator with different structural parameters and frictional factors.The results from the analyses provide valuable guidance for the design and motion control of the in-space cable-driven manipulator.Finally,a prototype of a single module is manufactured and tested.Ground experiments are carried out to verify the kinematic and dynamic models.展开更多
Predicting rock blasting outcomes in mining has been crucial since its inception.Blasting remains the most energy-and cost-efficient method for rock breaking and is often the only practical option.However,the mechanis...Predicting rock blasting outcomes in mining has been crucial since its inception.Blasting remains the most energy-and cost-efficient method for rock breaking and is often the only practical option.However,the mechanism is complex,influenced by various rock properties,explosives,and blast design parameters,making their effects difficult to quantify.Traditional stress-based models struggle with many parameters,such as stress and Poisson's ratio,which are challenging to measure in the field.Empirical models,though simpler,often oversimplify blast conditions.Both types of models are limited to simulating a few blastholes and cannot handle full-scale blasts involving hundreds of blastholes.However,modeling full-scale blasts with all blast design parameters is most required for modern mining applications.This paper presents a novel strain-based modeling approach for blasting and geomechanical applications,utilizing measurable variables such as particle velocity,strain,and displacement.By bypassing complex constitutive relations,strain-based models capture critical blasting trends and simulate full-scale blasts with full-blast design parameters with minimal calibration.The framework encompasses field strain measurements,model construction based on measurable variables,and laboratoryderived strain-failure criteria,each offering potential for future enhancement.Additionally,a standardized field test for site characterization is recommended.The approach is demonstrated through the Multiple Blasthole Fragmentation model,which simulates rock fragmentation and fragment strain during blasting,highlighting the practicality and effectiveness of strain-based modeling for multiple blasthole blasts.Moreover,this approach extends beyond blasting,with potential applications in highwall stability monitoring and other geomechanical applications.Strain-based modeling provides a simplified yet effective solution,avoiding the complexities of rock constitutive relations and field stress measurements while enabling full-blast design simulations for large-scale field blasts.展开更多
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.展开更多
Huperzine A(HupA) is a highly selective, reversible acetylcholinesterase(AChE) inhibitor that exhibits neuroprotective effects and is clinically used to manage benign memory decline.However, the specific relationship ...Huperzine A(HupA) is a highly selective, reversible acetylcholinesterase(AChE) inhibitor that exhibits neuroprotective effects and is clinically used to manage benign memory decline.However, the specific relationship between the pharmacokinetic(PK) profile of HupA and cerebral acetylcholine(ACh) dynamics remains poorly characterized. Here, we characterize the PK-pharmacodynamic(PD) properties of HupA in rats under both physiological and pathological conditions. Following a single intramuscular injection, HupA exhibits a short halflife but rapid brain penetration, while multiple dosing significantly enhances its brain exposure. In a middle cerebral artery occlusion(MCAO) rat model, HupA demonstrates increased brain distribution. Furthermore, HupA elevates ACh concentrations across multiple brain regions, concurrently modulating several monoamine neurotransmitters. Using a minimal physiologically based pharmacokinetic-pharmacodynamic(mPBPK-PD) modeling approach,cerebral ACh dynamics were accurately predicted based on the pharmacokinetics of HupA in systemic circulation. The developed mPBPK-PD model exhibits robust predictive performance and holds potential for guiding the optimization of clinical dosing regimens and improving the therapeutic efficacy of HupA.展开更多
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.展开更多
Separation bubbles forming on airfoils significantly influence aerodynamic behavior,particularly at low Reynolds numbers,making their accurate prediction a critical challenge in transition modelling.This study investi...Separation bubbles forming on airfoils significantly influence aerodynamic behavior,particularly at low Reynolds numbers,making their accurate prediction a critical challenge in transition modelling.This study investigates numerical modeling of a separation bubble and the effects of airfoil thickness and camber variation on the formation of the bubble dynamics at low Reynolds numbers.The numerical results were compared with the experimental results obtained from surface pressure distribution measurements,oil flow visualisation,and surface shear measurements to analyse the detailed flow behavior.The combination of pressure and flow visualisation techniques provided complementary insights,enabling a detailed characterisation of bubble formation.The results reveal that both the thickness and camber of the airfoil significantly influence the location,length,and stability of the bubble.At low Reynolds number flows(Re=0.5×10^(5)),particularly for highly cambered profiles,closer to the leading edge,separation and long bubbles were observed.As the Reynolds number increased,the separation point shifted to the leading edge,and reattachment became more likely.In numerical studies,transition models can accurately model the bubble initiation point;however,they often fail to model the bubble reattachment points accurately.This is due to the inadequacy of models that use empirical expressions for turbulence modelling,particularly in low Reynolds number flows,in their viscous modelling.In this study,it was concluded that transition onset terms,which specifically affect bubble formation,should be modified for more accurate modeling.展开更多
Interaction between the converter and the grid may lead to harmonic oscillations.The impedance-based method is an effective way to deal with the stability issue.In this study,the impedance-based method is used to inve...Interaction between the converter and the grid may lead to harmonic oscillations.The impedance-based method is an effective way to deal with the stability issue.In this study,the impedance-based method is used to investigate the small-signal stability of a cascaded 12-pulse line-commutated converter-based high-voltage direct current(LCC-HVDC)transmission system.In the modeling part,the impedance models of the single rectifier and inverter are established respectively with consideration to the effect of frequency coupling,which has improved the accuracy of the models.Based on the models,the AC impedance models of the cascaded LCC-HVDC transmission system are established both on the rectifier and inverter side.In the stability analysis part,the stability of the system is analyzed under different working conditions.The simulation results reveal that the established impedance model can properly represent the stability of this system.The findings of this study can provide a theoretical reference for the stability design and oscillation suppression strategy of LCC-HVDC transmission systems and LCC interconnected systems.展开更多
基金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.
基金Natural Science Foundation of China(No.42120104004)。
文摘The 1739 M8.0 Pingluo earthquake occurred around the Yinchuan Graben,bounded by the Helan Mountains to the west and the Ordos Block to the east.Seismological observations have shown that surface fault displacement reaches about 2–3 m,mainly by dip-slip motion along the Helanshan Piedmont Fault.However,the documented seismic intensity is distributed predominantly within the basin area,exhibiting a sharp asymmetry across the Helanshan Piedmont Fault.Thus,the general pattern of earthquake faulting is still under debate.We built a three-dimensional elastodynamic finiteelement model to reappraise the fault mechanism.In the model,predictions from synthetic rupture models,based on available observations and the earthquake scaling law,were used as an input with the split-node technique,and the effect of basin sediment on elastic wave propagation was considered.The numerical results show that if an earthquake occurred on the Helanshan Piedmont Fault characterized by a high-angle(70°)normal fault,earthquake shaking,as predicted from the modeled peak ground velocity and peak ground acceleration,has difficulty fitting the observed result,even when the effect of sediment amplification is considered.To better fit the observed shaking pattern,the dip angle of the Helanshan Piedmont Fault must be less than about 35°between the depths of about 8–27 km,where the coseismic slip may reach about 6 m.This result leads us to conclude that the 1739 M8.0 great earthquake likely occurred on a listric normal fault at depth,in agreement with the geometry of the Helanshan Piedmont Fault,as recently evidenced by seismic reflection explorations.This conclusion means that in an intracontinental setting,a reduction in the fault dip angle along the subsurface could increase the width of the fault in the elastic crust,making misalignment between the surface rupture and the isoseismals and resulting in an increase in the upper bound of earthquake magnitude relative to simple high-angle faulting.
基金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.
基金support from the National Key R&D Program of China(Grant No.2023YFB3709901)the National Natural Science Foundation of China(Grant No.U22A20171)+1 种基金China Baowu Low Carbon Metallurgy Innovation Foundation(Grant No.BWLCF202315)the High Steel Center(HSC)at North China University of Technology and University of Science and Technology Beijing,China.
文摘The dissolution of MgO-refractory into the slag had an obvious influence on the steel-slag reaction and the slag property,especially for high-aluminum steels.The dissolution behavior of MgO-refractory was investigated under various conditions,including the temperature,the initial steel composition,and the initial slag composition.A steel-slag-refractory kinetic model for high-aluminum steel was developed,which incorporated the process of MgO-refractory dissolution.The dependence of the MgO mass transfer coefficient k_(MgO)^(r)on temperature T during MgO-refractory dissolution process was established,as described by ln k_(MgO)^(r)=63,754/T+24.38524.It was indicated that the MgO dissolution rate was significantly influenced by the temperature.A higher temperature increased the dissolution rate of MgO.The initial steel composition had a slight impact on the MgO dissolution rate.Additionally,the initial slag composition strongly impacted the MgO saturation concentration and the dissolution rate.A lower initial Al_(2)O_(3)/SiO_(2)ratio increased the MgO dissolution rate.The steel-slag-refractory kinetic model accurately predicted the dissolution of MgO-refractory and the influence of dissolved MgO on the viscosity and composition change during steel-slag-refractory reactions.It was suggested that a higher temperature can hardly reduce the viscosity due to the dissolution of the MgO-refractory.
基金funded by the National Natural Science Foundation of China,grant number 52405341Foundation of National Key Laboratory of Computational Physics,grant number 6142A05QN24012+1 种基金Chongqing Science and Technology Committee,grant number CSTB2023NSCQ-MSX0363The Science and Technology Research Program of Chongqing Municipal Education Commission,grant number KJQN202301117.
文摘In materials science and engineering design,high-fidelity and high-efficiency numerical simulation has become a driving force for innovation and practical implementation.To address longstanding bottlenecks in the development of conventional material constitutive models—such as lengthy modeling cycles and difficulties in numerical implementation—this study proposes an intelligent modeling and code generation approach powered by large languagemodels.A structured knowledge base integrating constitutive theory,numerical algorithms,and UMAT(User Material)interface specifications is constructed,and a retrieval-augmented generation strategy is employed to establish an end-to-end workflow spanning experimental data parsing,constitutive model formulation,and automatic UMAT subroutine generation.Experimental results show that the method achieves high accuracy for both a classical Johnson–Cookmodel and a physics-informed neural network(PINN)model,with key parameter identification errors below 5%.Moreover,the automatically generated UMAT subroutines yield finite element simulation results in Abaqus that are highly consistent with theoretical predictions(coefficient of determination R2>0.98)while maintaining good numerical stability.This framework is currently focused on the automatic construction of rate-dependent elastoplastic material models,and its core method also provides a clear path for extending to other constitutive categories such as hyperelasticity and viscoelasticity.This work provides an effective technical route for the rapid development and reliable numerical implementation of material constitutive models,significantly advancing the intelligence level of computational mechanics research and improving engineering application efficiency.
基金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.
基金supported by Excellent Research Group Project for Multiphase Evolution in Hyper-Gravity of the National Natural Science Foundation of China(No.52588202)。
文摘The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we propose a new approach by combining centrifugal analog modeling with numerical simulation to simulate the tectonic uplift history of the plateau based on the lower crustal flow model,and to investigate the material migration characteristics and the influence of crustal motion velocity and ductile layer viscosity on the plateau tectonic geomorphology.The models reproduce steep-sided flat-topped geomorphic features and clockwise rotation of the material at eastern Himalayan Syntaxis,verifying the rationality of the models.The results show that the greater the crustal motion velocity and the greater the ductile layer viscosity,the steeper the terrain change;and conversely,the smaller the crustal motion velocity and the smaller the ductile layer viscosity,the gentler the terrain change.This study further indicates that the weak lower crust plays an important role in the formation of geomorphic features and material migration characteristics of Qinghai-Tibet Plateau,and provides a new insight for the study of the uplift mechanism of the Tibetan Plateau.
基金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.
基金co-supported by the National Natural Science Foundation of China(Nos.12102034 and 12125201)the Open Fund of State Key Laboratory of Robotics and Systems(HIT),China。
文摘In-space cable-driven manipulators exhibit several advantages,such as a large range of motion,high dexterity,and lightweight structure.However,kinematic and dynamic analysis play an essential role in designing a cable-driven manipulator.In this paper,the kinematic analysis of a type of cable-driven manipulator is performed,and a motion planning scheme is conducted to actuate this manipulator.Moreover,a flexible multi-body dynamic model of a cable-driven manipulator considering the frictional contact between the cables and pulleys is established.To describe properties such as flexibility,vibration,and variable length of the cable,this paper utilizes reducedorder beam elements of the Absolute Nodal Coordinates Formulation(ANCF)in Arbitrary Lagrangian Eulerian(ALE)framework.Additionally,a virtual element is introduced to model the contact segment in the cable-pulley system.A tension decay factor is employed to account for the friction in the contact segment.To validate the proposed method,a semi-analytical model based on D'Alembert's principle is established.Cross-verification is performed to validate the accuracy of both models.The model is further applied to simulate the rotation of the cable-driven manipulator with different structural parameters and frictional factors.The results from the analyses provide valuable guidance for the design and motion control of the in-space cable-driven manipulator.Finally,a prototype of a single module is manufactured and tested.Ground experiments are carried out to verify the kinematic and dynamic models.
文摘Predicting rock blasting outcomes in mining has been crucial since its inception.Blasting remains the most energy-and cost-efficient method for rock breaking and is often the only practical option.However,the mechanism is complex,influenced by various rock properties,explosives,and blast design parameters,making their effects difficult to quantify.Traditional stress-based models struggle with many parameters,such as stress and Poisson's ratio,which are challenging to measure in the field.Empirical models,though simpler,often oversimplify blast conditions.Both types of models are limited to simulating a few blastholes and cannot handle full-scale blasts involving hundreds of blastholes.However,modeling full-scale blasts with all blast design parameters is most required for modern mining applications.This paper presents a novel strain-based modeling approach for blasting and geomechanical applications,utilizing measurable variables such as particle velocity,strain,and displacement.By bypassing complex constitutive relations,strain-based models capture critical blasting trends and simulate full-scale blasts with full-blast design parameters with minimal calibration.The framework encompasses field strain measurements,model construction based on measurable variables,and laboratoryderived strain-failure criteria,each offering potential for future enhancement.Additionally,a standardized field test for site characterization is recommended.The approach is demonstrated through the Multiple Blasthole Fragmentation model,which simulates rock fragmentation and fragment strain during blasting,highlighting the practicality and effectiveness of strain-based modeling for multiple blasthole blasts.Moreover,this approach extends beyond blasting,with potential applications in highwall stability monitoring and other geomechanical applications.Strain-based modeling provides a simplified yet effective solution,avoiding the complexities of rock constitutive relations and field stress measurements while enabling full-blast design simulations for large-scale field blasts.
基金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.
基金supported by the National Key Research and Development Program of China (No. 2024YFA1308200)the National Natural Science Foundation of China (Nos. 82274009 and81973556)。
文摘Huperzine A(HupA) is a highly selective, reversible acetylcholinesterase(AChE) inhibitor that exhibits neuroprotective effects and is clinically used to manage benign memory decline.However, the specific relationship between the pharmacokinetic(PK) profile of HupA and cerebral acetylcholine(ACh) dynamics remains poorly characterized. Here, we characterize the PK-pharmacodynamic(PD) properties of HupA in rats under both physiological and pathological conditions. Following a single intramuscular injection, HupA exhibits a short halflife but rapid brain penetration, while multiple dosing significantly enhances its brain exposure. In a middle cerebral artery occlusion(MCAO) rat model, HupA demonstrates increased brain distribution. Furthermore, HupA elevates ACh concentrations across multiple brain regions, concurrently modulating several monoamine neurotransmitters. Using a minimal physiologically based pharmacokinetic-pharmacodynamic(mPBPK-PD) modeling approach,cerebral ACh dynamics were accurately predicted based on the pharmacokinetics of HupA in systemic circulation. The developed mPBPK-PD model exhibits robust predictive performance and holds potential for guiding the optimization of clinical dosing regimens and improving the therapeutic efficacy of HupA.
基金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.
基金the Scientific and Technological Research Council of Turkey(TÜB˙ITAK)for support under project number:122M826to the Scientific Research Projects Unit of Erciyes University under contract No.:FYL-2023-13162 and FYL-2024-13701.
文摘Separation bubbles forming on airfoils significantly influence aerodynamic behavior,particularly at low Reynolds numbers,making their accurate prediction a critical challenge in transition modelling.This study investigates numerical modeling of a separation bubble and the effects of airfoil thickness and camber variation on the formation of the bubble dynamics at low Reynolds numbers.The numerical results were compared with the experimental results obtained from surface pressure distribution measurements,oil flow visualisation,and surface shear measurements to analyse the detailed flow behavior.The combination of pressure and flow visualisation techniques provided complementary insights,enabling a detailed characterisation of bubble formation.The results reveal that both the thickness and camber of the airfoil significantly influence the location,length,and stability of the bubble.At low Reynolds number flows(Re=0.5×10^(5)),particularly for highly cambered profiles,closer to the leading edge,separation and long bubbles were observed.As the Reynolds number increased,the separation point shifted to the leading edge,and reattachment became more likely.In numerical studies,transition models can accurately model the bubble initiation point;however,they often fail to model the bubble reattachment points accurately.This is due to the inadequacy of models that use empirical expressions for turbulence modelling,particularly in low Reynolds number flows,in their viscous modelling.In this study,it was concluded that transition onset terms,which specifically affect bubble formation,should be modified for more accurate modeling.
基金supported in part by the National Natural Science Foundation of China under 52125704 and 51937001.
文摘Interaction between the converter and the grid may lead to harmonic oscillations.The impedance-based method is an effective way to deal with the stability issue.In this study,the impedance-based method is used to investigate the small-signal stability of a cascaded 12-pulse line-commutated converter-based high-voltage direct current(LCC-HVDC)transmission system.In the modeling part,the impedance models of the single rectifier and inverter are established respectively with consideration to the effect of frequency coupling,which has improved the accuracy of the models.Based on the models,the AC impedance models of the cascaded LCC-HVDC transmission system are established both on the rectifier and inverter side.In the stability analysis part,the stability of the system is analyzed under different working conditions.The simulation results reveal that the established impedance model can properly represent the stability of this system.The findings of this study can provide a theoretical reference for the stability design and oscillation suppression strategy of LCC-HVDC transmission systems and LCC interconnected systems.
