In this paper,the multi-agent model about shop logistics is set up.This model has 8 agents:raw materials stock agent,process agent,testing agent,transition agent,production information agent,scheduling agent,process a...In this paper,the multi-agent model about shop logistics is set up.This model has 8 agents:raw materials stock agent,process agent,testing agent,transition agent,production information agent,scheduling agent,process agent and stock agent.The scheduling agent has three subagents:manager agent(MA),resource agent(RA)and part agent(PA).MA,PA and RA are communicating equally that guarantees agility of the whole MAS system.The part tasks pass between MA,RA and PA as an integer,which can guarantee the consistency of the data.We use a detailed example about shop logistics scheduling in a semiconductor company to explain the principle.In this example,we use two scheduling strategies:FCFS and SPT.The result data indicates that the average flow time and lingering ratio are changed using different strategy.It is proves that the multi-agent scheduling is useful.展开更多
With the increasing role of social media in information dissemination,effectively simulating and analyzing public event dynamics has become a key research focus.We present an interactive visual analysis system for sim...With the increasing role of social media in information dissemination,effectively simulating and analyzing public event dynamics has become a key research focus.We present an interactive visual analysis system for simulating social media events using multi-agent models powered by large language models(LLMs).By modeling agents with diverse characteristics,the system explores how agents perceive information,adjust their emotions and stances,provide feedback,and influence the trajectory of events.The system integrates real-time interactive simulation with multi-perspective visualization,enabling users to investigate event trajectories and key influencing factors under varied configurations.Theoretical work standardizes agent attributes and interaction mechanisms,supporting realistic simulation of social media behaviors.Evaluation through indicators and case studies demonstrates the system’s effectiveness and adaptability,offering a novel tool for public event analysis across open social platforms.展开更多
Multi-agent systems often require good interoperability in the process of completing their assigned tasks.This paper first models the static structure and dynamic behavior of multiagent systems based on layered weight...Multi-agent systems often require good interoperability in the process of completing their assigned tasks.This paper first models the static structure and dynamic behavior of multiagent systems based on layered weighted scale-free community network and susceptible-infected-recovered(SIR)model.To solve the problem of difficulty in describing the changes in the structure and collaboration mode of the system under external factors,a two-dimensional Monte Carlo method and an improved dynamic Bayesian network are used to simulate the impact of external environmental factors on multi-agent systems.A collaborative information flow path optimization algorithm for agents under environmental factors is designed based on the Dijkstra algorithm.A method for evaluating system interoperability is designed based on simulation experiments,providing reference for the construction planning and optimization of organizational application of the system.Finally,the feasibility of the method is verified through case studies.展开更多
The similarities and differences between the container terminal logistics system(CTLS)and the Harvard-architecture computer system are compared in terms of organization and architecture.The mapping relation and the mo...The similarities and differences between the container terminal logistics system(CTLS)and the Harvard-architecture computer system are compared in terms of organization and architecture.The mapping relation and the modeling framework of the CTLS are presented based on multi-agent,and the successful algorithms in the computer domain are applied to the modeling framework,such as the dynamic priority and multilevel feedback scheduling algorithm.In addition,a model and simulation on a certain quay at Shanghai harbor is built up on the AnyLogic platform to support the decision-making of terminal on service cost.It validates the feasibility and creditability of the above systematic methodology.展开更多
With the new characteristics of global cooperation in supply chains being synthetically considered,a hybrid model to the cooperative negotiation process for the order distribution in supply chain is mainly studied.Aft...With the new characteristics of global cooperation in supply chains being synthetically considered,a hybrid model to the cooperative negotiation process for the order distribution in supply chain is mainly studied.After reviewing and analyzing some main domestic and overseas processes in cooperative negotiation modeling in supply chain,some problems are subsequently pointed out.For example,the traditional simple multi-agent system(MAS)frameworks which have some limitations,are not suitable for solving modeling complex systems.To solve these problems,thinking with the aid of the multi-agent structure and complex system modeling,the manufacturing supply chain is taken as an example,and a time Petri net production model is adopted to decompose the materials.And then a cooperative negotiation model for the order distribution in supply chain is constructed based on combining multi-agent techniques with time Petri net modeling.The simulation results reveal that the above model helps solve the problems of cooperative negotiation in supply chains.展开更多
With the advent of sixth-generation mobile communications(6G),space-air-ground integrated networks have become mainstream.This paper focuses on collaborative scheduling for mobile edge computing(MEC)under a three-tier...With the advent of sixth-generation mobile communications(6G),space-air-ground integrated networks have become mainstream.This paper focuses on collaborative scheduling for mobile edge computing(MEC)under a three-tier heterogeneous architecture composed of mobile devices,unmanned aerial vehicles(UAVs),and macro base stations(BSs).This scenario typically faces fast channel fading,dynamic computational loads,and energy constraints,whereas classical queuing-theoretic or convex-optimization approaches struggle to yield robust solutions in highly dynamic settings.To address this issue,we formulate a multi-agent Markov decision process(MDP)for an air-ground-fused MEC system,unify link selection,bandwidth/power allocation,and task offloading into a continuous action space and propose a joint scheduling strategy that is based on an improved MATD3 algorithm.The improvements include Alternating Layer Normalization(ALN)in the actor to suppress gradient variance,Residual Orthogonalization(RO)in the critic to reduce the correlation between the twin Q-value estimates,and a dynamic-temperature reward to enable adaptive trade-offs during training.On a multi-user,dual-link simulation platform,we conduct ablation and baseline comparisons.The results reveal that the proposed method has better convergence and stability.Compared with MADDPG,TD3,and DSAC,our algorithm achieves more robust performance across key metrics.展开更多
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.展开更多
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.展开更多
Accurate mechanical modeling is essential for robotic belt grinding(RBG), a process characterized by compliant contact mechanisms that make force prediction particularly challenging. However, existing mechanical model...Accurate mechanical modeling is essential for robotic belt grinding(RBG), a process characterized by compliant contact mechanisms that make force prediction particularly challenging. However, existing mechanical models predominantly focus on macroscale compliance while neglecting grain-scale compliant motion. Moreover, abrasive grains are typically idealized as regular shapes, overlooking the inherent stochasticity of real grain geometries. This study proposes a shapeequivalence method for modeling stochastic abrasive grains and develops a multiscale compliant force model for RBG. Specifically, an individual grain is represented as a polygonal pyramid with stochastic edges that is mathematically equivalent to a cone;this method unifies the treatment of grain geometries and streamlines the modeling process. The mathematical equivalence relationship for random grain shapes is further derived based on a grain-compliant contact model. By integrating grain geometric characteristics and progressive grain wear, an analytical mechanical model that captures both the static contact force and dynamic grinding force is established, thereby describing the transition from grain-workpiece compliant interaction to belt-workpiece elastic contact. Grinding experiments were conducted using abrasive belts with different grain shape distributions to validate the model. The results demonstrated reliable predictions of the tangential grinding force and its component characteristics. Additional analyses were performed to reveal how the tangential grinding force varies with wear time and grinding parameters.展开更多
基金Supported by the Zhejiang Province Science Foundation of China(M703022)
文摘In this paper,the multi-agent model about shop logistics is set up.This model has 8 agents:raw materials stock agent,process agent,testing agent,transition agent,production information agent,scheduling agent,process agent and stock agent.The scheduling agent has three subagents:manager agent(MA),resource agent(RA)and part agent(PA).MA,PA and RA are communicating equally that guarantees agility of the whole MAS system.The part tasks pass between MA,RA and PA as an integer,which can guarantee the consistency of the data.We use a detailed example about shop logistics scheduling in a semiconductor company to explain the principle.In this example,we use two scheduling strategies:FCFS and SPT.The result data indicates that the average flow time and lingering ratio are changed using different strategy.It is proves that the multi-agent scheduling is useful.
基金supported by the Natural Science Foundation of China(NSFC No.62472099)Ji Hua Laboratory S&T Program(X250881UG250).
文摘With the increasing role of social media in information dissemination,effectively simulating and analyzing public event dynamics has become a key research focus.We present an interactive visual analysis system for simulating social media events using multi-agent models powered by large language models(LLMs).By modeling agents with diverse characteristics,the system explores how agents perceive information,adjust their emotions and stances,provide feedback,and influence the trajectory of events.The system integrates real-time interactive simulation with multi-perspective visualization,enabling users to investigate event trajectories and key influencing factors under varied configurations.Theoretical work standardizes agent attributes and interaction mechanisms,supporting realistic simulation of social media behaviors.Evaluation through indicators and case studies demonstrates the system’s effectiveness and adaptability,offering a novel tool for public event analysis across open social platforms.
基金supported by the Key R&D Projects in Jiangsu Province(BE2021729)the Key Primary Research Project of Primary Strengthening Program(KYZYJKKCJC23001).
文摘Multi-agent systems often require good interoperability in the process of completing their assigned tasks.This paper first models the static structure and dynamic behavior of multiagent systems based on layered weighted scale-free community network and susceptible-infected-recovered(SIR)model.To solve the problem of difficulty in describing the changes in the structure and collaboration mode of the system under external factors,a two-dimensional Monte Carlo method and an improved dynamic Bayesian network are used to simulate the impact of external environmental factors on multi-agent systems.A collaborative information flow path optimization algorithm for agents under environmental factors is designed based on the Dijkstra algorithm.A method for evaluating system interoperability is designed based on simulation experiments,providing reference for the construction planning and optimization of organizational application of the system.Finally,the feasibility of the method is verified through case studies.
基金The National Key Technology R&D Program of China during the 11th Five-Year Plan Period(No.2006BAH02A06)
文摘The similarities and differences between the container terminal logistics system(CTLS)and the Harvard-architecture computer system are compared in terms of organization and architecture.The mapping relation and the modeling framework of the CTLS are presented based on multi-agent,and the successful algorithms in the computer domain are applied to the modeling framework,such as the dynamic priority and multilevel feedback scheduling algorithm.In addition,a model and simulation on a certain quay at Shanghai harbor is built up on the AnyLogic platform to support the decision-making of terminal on service cost.It validates the feasibility and creditability of the above systematic methodology.
基金The National Natural Science Foundation of China(No.70401013)the National Key Technology R&D Program of China during the 11th Five-Year Plan Period(No.2006BAH02A06)
文摘With the new characteristics of global cooperation in supply chains being synthetically considered,a hybrid model to the cooperative negotiation process for the order distribution in supply chain is mainly studied.After reviewing and analyzing some main domestic and overseas processes in cooperative negotiation modeling in supply chain,some problems are subsequently pointed out.For example,the traditional simple multi-agent system(MAS)frameworks which have some limitations,are not suitable for solving modeling complex systems.To solve these problems,thinking with the aid of the multi-agent structure and complex system modeling,the manufacturing supply chain is taken as an example,and a time Petri net production model is adopted to decompose the materials.And then a cooperative negotiation model for the order distribution in supply chain is constructed based on combining multi-agent techniques with time Petri net modeling.The simulation results reveal that the above model helps solve the problems of cooperative negotiation in supply chains.
文摘With the advent of sixth-generation mobile communications(6G),space-air-ground integrated networks have become mainstream.This paper focuses on collaborative scheduling for mobile edge computing(MEC)under a three-tier heterogeneous architecture composed of mobile devices,unmanned aerial vehicles(UAVs),and macro base stations(BSs).This scenario typically faces fast channel fading,dynamic computational loads,and energy constraints,whereas classical queuing-theoretic or convex-optimization approaches struggle to yield robust solutions in highly dynamic settings.To address this issue,we formulate a multi-agent Markov decision process(MDP)for an air-ground-fused MEC system,unify link selection,bandwidth/power allocation,and task offloading into a continuous action space and propose a joint scheduling strategy that is based on an improved MATD3 algorithm.The improvements include Alternating Layer Normalization(ALN)in the actor to suppress gradient variance,Residual Orthogonalization(RO)in the critic to reduce the correlation between the twin Q-value estimates,and a dynamic-temperature reward to enable adaptive trade-offs during training.On a multi-user,dual-link simulation platform,we conduct ablation and baseline comparisons.The results reveal that the proposed method has better convergence and stability.Compared with MADDPG,TD3,and DSAC,our algorithm achieves more robust performance across key metrics.
基金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.
基金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.
基金supported by the National Natural Science Foundation of China (Grant Nos.52505554,52575571)the Postdoctoral Fellowship Program of CPSF (Grant No.GZB20250348)。
文摘Accurate mechanical modeling is essential for robotic belt grinding(RBG), a process characterized by compliant contact mechanisms that make force prediction particularly challenging. However, existing mechanical models predominantly focus on macroscale compliance while neglecting grain-scale compliant motion. Moreover, abrasive grains are typically idealized as regular shapes, overlooking the inherent stochasticity of real grain geometries. This study proposes a shapeequivalence method for modeling stochastic abrasive grains and develops a multiscale compliant force model for RBG. Specifically, an individual grain is represented as a polygonal pyramid with stochastic edges that is mathematically equivalent to a cone;this method unifies the treatment of grain geometries and streamlines the modeling process. The mathematical equivalence relationship for random grain shapes is further derived based on a grain-compliant contact model. By integrating grain geometric characteristics and progressive grain wear, an analytical mechanical model that captures both the static contact force and dynamic grinding force is established, thereby describing the transition from grain-workpiece compliant interaction to belt-workpiece elastic contact. Grinding experiments were conducted using abrasive belts with different grain shape distributions to validate the model. The results demonstrated reliable predictions of the tangential grinding force and its component characteristics. Additional analyses were performed to reveal how the tangential grinding force varies with wear time and grinding parameters.
