Patient-specific modeling of left ventricle mechanics

We present a computational framework for the study of cardiac motion.The bio-mechanical model captures the passive and active properties of the cardiac tissue as well as the fiber architecture.We focus on the analysis of deformations of a beating left ventricle(LV),comparing numerical simulations with real data acquired by echocardiography.The goal is to determine the clinical relevance of the LV strains pattern and to investigate the relationships between that pattern and the arrangement of myocardial fibers.The proposed framework could in principle be used for a wide range of clinical applications.

A promising approach for quantifying focal stroke modeling and assessing stroke progression:optical resolution photoacoustic microscopy photothrombosis

To investigate the mechanisms underlying the onset and progression of ischemic stroke,some methods have been proposed that can simultaneously monitor and create embolisms in the animal cerebral cortex.However,these methods often require complex systems and the effect of age on cerebral embolism has not been adequately studied,although ischemic stroke is strongly age-related.In this study,we propose an optical-resolution photoacoustic microscopy-based visualized photothrombosis methodology to create and monitor ischemic stroke in mice simultaneously using a 532 nm pulsed laser.We observed the molding process in mice of different ages and presented age-dependent vascular embolism differentiation.Moreover,we integrated optical coherence tomography angiography to investigate age-associated trends in cerebrovascular variability following a stroke.Our imaging data and quantitative analyses underscore the differential cerebrovascular responses to stroke in mice of different ages,thereby highlighting the technique's potential for evaluating cerebrovascular health and unraveling age-related mechanisms involved in ischemic strokes.

Modeling and Comprehensive Review of Signaling Storms in 3GPP-Based Mobile Broadband Networks:Causes,Solutions,and Countermeasures

Control signaling is mandatory for the operation and management of all types of communication networks,including the Third Generation Partnership Project(3GPP)mobile broadband networks.However,they consume important and scarce network resources such as bandwidth and processing power.There have been several reports of these control signaling turning into signaling storms halting network operations and causing the respective Telecom companies big financial losses.This paper draws its motivation from such real network disaster incidents attributed to signaling storms.In this paper,we present a thorough survey of the causes,of the signaling storm problems in 3GPP-based mobile broadband networks and discuss in detail their possible solutions and countermeasures.We provide relevant analytical models to help quantify the effect of the potential causes and benefits of their corresponding solutions.Another important contribution of this paper is the comparison of the possible causes and solutions/countermeasures,concerning their effect on several important network aspects such as architecture,additional signaling,fidelity,etc.,in the form of a table.This paper presents an update and an extension of our earlier conference publication.To our knowledge,no similar survey study exists on the subject.

Computerised tomography vs magnetic resonance imaging for modeling of patient-specific instrumentation in total knee arthroplasty

AIM: To summarise and compare currently available evidence regarding accuracy of pre-operative imaging, which is one of the key choices for surgeons contemplating patient-specific instrumentation(PSI) surgery.METHODS: The MEDLINE and EMBASE medical literature databases were searched, from January 1990 to December 2013, to identify relevant studies. The data from several clinical studies was assimilated to allow appreciation and comparison of the accuracy of each modality. The overall accuracy of each modality was calculated as proportion of outliers > 3% in the coronal plane of both computerised tomography(CT) or magnetic resonance imaging(MRI). RESULTS: Seven clinical studies matched our inclusion criteria for comparison and were included in our study for statistical analysis. Three of these reported series using MRI and four with CT. Overall percentage of outliers > 3% in patients with CT-based PSI systems was 12.5% vs 16.9% for MRI-based systems. These results were not statistically significant. CONCLUSION: Although many studies have been undertaken to determine the ideal pre-operative imaging modality, conclusions remain speculative in the absence of long term data. Ultimately, information regarding accuracy of CT and MRI will be the main determining factor. Increased accuracy of pre-operative imaging could result in longer-term savings, and reduced accumulated dose of radiation by eliminating the need for post-operative imaging and revision surgery.

Patient-specific modeling of facial soft tissue based on radial basis functions transformations of a standard three-dimensional finite element model

Background An important purpose of orthodontic treatment is to gain the harmonic soft tissue profile. This article describes a novel way to build patient-specific models of facial soft tissues by transforming a standard finite element （FE） model into one that has two stages： a first transformation and a second transformation, so as to evaluate the facial soft tissue changes after orthodontic treatment for individual patients. Methods The radial basis functions （RBFs） interpolation method was used to transform the standard FE model into a patient-specific one based on landmark points. A combined strategy for selecting landmark points was developed in this study： manually for the first transformation and automatically for the second transformation. Four typical patients were chosen to validate the effectiveness of this transformation method. Results The results showed good similarity between the transformed FE models and the computed tomography （CT） models. The absolute values of average deviations were in the range of 0.375-0.700 mm at the lip-mouth region after the first transformation, and they decreased to a range of 0.116-0.286 mm after the second transformation. Conclusions The modeling results show that the second transformation resulted in enhanced accuracy compared to the first transformation. Because of these results, a third transformation is usually not necessary.

Patient-Specific Echo-Based Fluid-Structure Interaction Modeling Study of Blood Flow in the Left Ventricle with Infarction and Hypertension

Understanding cardiac blood flow behaviors is of importance for cardiovascular research and clinical assessment of ventricle functions.Patient-specific Echo-based left ventricle(LV)fluid-structure interaction(FSI)models were introduced to perform ventricle mechanical analysis,investigate flow behaviors,and evaluate the impact of myocardial infarction(MI)and hypertension on blood flow in the LV.Echo image data were acquired from 3 patients with consent obtained:one healthy volunteer(P1),one hypertension patient(P2),and one patient who had an inferior and posterior myocardial infarction(P3).The nonlinear Mooney-Rivlin model was used for ventricle tissue with material parameter values chosen to match echo-measure LV volume data.Using the healthy case as baseline,LV with MI had lower peak flow velocity(30%lower at beginejection)and hypertension LV had higher peak flow velocity(16%higher at begin-filling).The vortex area(defined as the area with vorticity>0)for P3 was 19%smaller than that of P1.The vortex area for P2 was 12%smaller than that of P1.At peak of filling,the maximum flow shear stress(FSS)for P2 and P3 were 390%higher and 63%lower than that of P1,respectively.Meanwhile,LV stress and strain of P2 were 41%and 15%higher than those of P1,respectively.LV stress and strain of P3 were 36%and 42%lower than those of P1,respectively.In conclusion,FSI models could provide both flow and structural stress/strain information which would serve as the base for further cardiovascular investigations related to disease initiation,progression,and treatment strategy selections.Large-scale studies are needed to validate our findings.

The relationship between compartment models and their stochastic counterparts:A comparative study with examples of the COVID-19 epidemic modeling

Deterministic compartment models(CMs)and stochastic models,including stochastic CMs and agent-based models,are widely utilized in epidemic modeling.However,the relationship between CMs and their corresponding stochastic models is not well understood.The present study aimed to address this gap by conducting a comparative study using the susceptible,exposed,infectious,and recovered(SEIR)model and its extended CMs from the coronavirus disease 2019 modeling literature.We demonstrated the equivalence of the numerical solution of CMs using the Euler scheme and their stochastic counterparts through theoretical analysis and simulations.Based on this equivalence,we proposed an efficient model calibration method that could replicate the exact solution of CMs in the corresponding stochastic models through parameter adjustment.The advancement in calibration techniques enhanced the accuracy of stochastic modeling in capturing the dynamics of epidemics.However,it should be noted that discrete-time stochastic models cannot perfectly reproduce the exact solution of continuous-time CMs.Additionally,we proposed a new stochastic compartment and agent mixed model as an alternative to agent-based models for large-scale population simulations with a limited number of agents.This model offered a balance between computational efficiency and accuracy.The results of this research contributed to the comparison and unification of deterministic CMs and stochastic models in epidemic modeling.Furthermore,the results had implications for the development of hybrid models that integrated the strengths of both frameworks.Overall,the present study has provided valuable epidemic modeling techniques and their practical applications for understanding and controlling the spread of infectious diseases.

Displacement field reconstruction in landslide physical modeling by using a terrain laser scanner e Part 2:Application and large strain/displacement and water effect analysis

Deformation analysis is fundamental in geotechnical modeling.Nevertheless,there is still a lack of an effective method to obtain the deformation field under various experimental conditions.In this study,we introduce a processebased physical modeling of a pileereinforced reservoir landslide and present an improved deformation analysis involving large strains and water effects.We collect multieperiod point clouds using a terrain laser scanner and reconstruct its deformation field through a point cloud processing workflow.The results show that this method can accurately describe the landslide surface deformation at any time and area by both scalar and vector fields.The deformation fields in different profiles of the physical model and different stages of the evolutionary process provide adequate and detailed landslide information.We analyze the large strain upstream of the pile caused by the pile installation and the consequent violent deformation during the evolutionary process.Furthermore,our method effectively overcomes the challenges of identifying targets commonly encountered in geotechnical modeling where water effects are considered and targets are polluted,which facilitates the deformation analysis at the wading area in a reservoir landslide.Eventually,combining subsurface deformation as well as numerical modeling,we comprehensively analyze the kinematics and failure mechanisms of this complicated object involving landslides and pile foundations as well as water effects.This method is of great significance for any geotechnical modeling concerning large-strain analysis and water effects.

Influence of Clip Locations on Intraaneurysmal Flow Dynamics in Patient-specific Anterior Communicating Aneurysm Models with Different Aneurysmal Angle

To improve aneurysm treatment,this study examined the influence of clip locations on hemodynamic factors in patient-specific anterior communicating artery(ACoA)aneurysms with different aneurysmal angle.We proposed a simplified classification of ACoA aneurysms using aneurysmal angle,defined by the angle of pivot of the aneurysmal dome and the virtual two-dimensional plane created by both proximal A2 segments of anterior cerebral artery(ACA).ACoA aneurysms with three different aneurysmal angles,which are 15°,80°and 120°,were analyzed in our study.In this work,we obtained hemodynamics before and after clipping surgery with three clip locations based on clinical clipping strategies in three ACoA aneurysms with different aneurysm angles.Results showed that local high pressure occurs at impingement region of the ACoA aneurysm before clipping and new impingement region close to the clipping location after clipping treatment.For clipping the aneurysm with aneurysmal angle 15°and a wide neck,wall shear stress(WSS)distribution is more uniform when the clipping angle of two clips close to 180°comparing with other two angles.In addition,for clipping the aneurysm with aneurysmal angle 80°and 120°,local high pressure appears on new impingement region and high WSS distributes around the clipping location when the clip plane is normal to the direction of inflow of aneurysm from the dominance of A1 segment of ACA.Hence,we should avoid the impingement of inflow from the A1 segment and choose a favorable clipping location for the fastness of clip.The results of our study could preoperatively give a useful information to the decision of surgical plan.

Predictive modeling for postoperative delirium in elderly patients with abdominal malignancies using synthetic minority oversampling technique

BACKGROUND Postoperative delirium,particularly prevalent in elderly patients after abdominal cancer surgery,presents significant challenges in clinical management.AIM To develop a synthetic minority oversampling technique(SMOTE)-based model for predicting postoperative delirium in elderly abdominal cancer patients.METHODS In this retrospective cohort study,we analyzed data from 611 elderly patients who underwent abdominal malignant tumor surgery at our hospital between September 2020 and October 2022.The incidence of postoperative delirium was recorded for 7 d post-surgery.Patients were divided into delirium and non-delirium groups based on the occurrence of postoperative delirium or not.A multivariate logistic regression model was used to identify risk factors and develop a predictive model for postoperative delirium.The SMOTE technique was applied to enhance the model by oversampling the delirium cases.The model’s predictive accuracy was then validated.RESULTS In our study involving 611 elderly patients with abdominal malignant tumors,multivariate logistic regression analysis identified significant risk factors for postoperative delirium.These included the Charlson comorbidity index,American Society of Anesthesiologists classification,history of cerebrovascular disease,surgical duration,perioperative blood transfusion,and postoperative pain score.The incidence rate of postoperative delirium in our study was 22.91%.The original predictive model(P1)exhibited an area under the receiver operating characteristic curve of 0.862.In comparison,the SMOTE-based logistic early warning model(P2),which utilized the SMOTE oversampling algorithm,showed a slightly lower but comparable area under the curve of 0.856,suggesting no significant difference in performance between the two predictive approaches.CONCLUSION This study confirms that the SMOTE-enhanced predictive model for postoperative delirium in elderly abdominal tumor patients shows performance equivalent to that of traditional methods,effectively addressing data imbalance.

Volume-averaged modeling of multiphase solidification with equiaxed crystal sedimentation in a steel ingot

Macrosegregation is a critical factor that limits the mechanical properties of materials.The impact of equiaxed crystal sedimentation on macrosegregation has been extensively studied,as it plays a significant role in determining the distribution of alloying elements and impurities within a material.To improve macrosegregation in steel connecting shafts,a multiphase solidification model that couples melt flow,heat transfer,microstructure evolution,and solute transport was established based on the volume-averaged Eulerian-Eulerian approach.In this model,the effects of liquid phase,equiaxed crystals,columnar dendrites,and columnar-to-equiaxed transition(CET)during solidification and evolution of microstructure can be considered simultaneously.The sedimentation of equiaxed crystals contributes to negative macrosegregation,where regions between columnar dendrites and equiaxed crystals undergo significant A-type positive macrosegregation due to the CET.Additionally,noticeable positive macrosegregation occurs in the area of final solidification in the ingot.The improvement in macrosegregation is beneficial for enhancing the mechanical properties of connecting shafts.To mitigate the thermal convection of molten steel resulting from excessive superheating,reducing the superheating during casting without employing external fields or altering the design of the ingot mold is indeed an effective approach to control macrosegregation.

Reactive transport modeling constraints on the complex genesis of a lacustrine dolomite reservoir:A case from the Eocene Qaidam Basin,China

Reactive transport modeling(RTM)is an emerging method used to address geological issues in diagenesis research.However,the extrapolation of RTM results to practical reservoir prediction is not sufficiently understood.This paper presents a case study of the Eocene Qaidam Basin that combines RTM results with petrological and mineralogical evidence.The results show that the Eocene Xiaganchaigou Formation is characterized by mixed siliciclastic-carbonate-evaporite sedimentation in a semiclosed saline lacustrine environment.Periodic evaporation and salinization processes during the syngeneticpenecontemporaneous stage gave rise to the replacive genesis of dolomites and the cyclic enrichment of dolomite in the middle-upper parts of the meter-scale depositional sequences.The successive change in mineral paragenesis from terrigenous clastics to carbonates to evaporites was reconstructed using RTM simulations.Parametric uncertainty analyses further suggest that the evaporation intensity(brine salinity)and particle size of sediments(reactive surface area)were important rate-determining factors in the dolomitization,as shown by the relatively higher reaction rates under conditions of higher brine salinity and fine-grained sediments.Combining the simulation results with measured mineralogical and reservoir physical property data indicates that the preservation of original intergranular pores and the generation of porosity via replacive dolomitization were the major formation mechanisms of the distinctive lacustrine dolomite reservoirs(widespread submicron intercrystalline micropores)in the Eocene Qaidam Basin.The results confirm that RTM can be effectively used in geological studies,can provide a better general understanding of the dolomitizing fluid-rock interactions,and can shed light on the spatiotemporal evolution of mineralogy and porosity during dolomitization and the formation of lacustrine dolomite reservoirs.

Scenario analysis of the Indonesia carbon tax impact on carbon emissions using system dynamics modeling and STIRPAT model

This study aims to develop a system dynamic(SD)forecasting model based on the STIRPAT model to forecast the effect of an IDR 30 per kg CO_(2)e carbon tax on carbon emissions,estimate future carbon emissions under ten scenarios,without and with the carbon tax,and estimate the environmental Kuznets curve(EKC)to predict Indonesia’s carbon emission peak.Carbon emission drivers in this study are decomposed into several factors,namely energy structure,energy intensity,industrial structure,GDP per capita,population,and fixed-asset investment.This study included nuclear power utilization starting in 2038.The research gaps addressed by this study compared to previous research are(1)use of the ex-ante approach,(2)inclusion of nuclear power plants,(3)testing the EKC hypothesis,and(4)contribution to government policy.The simulation results show that under the carbon tax,carbon emissions can be reduced by improving renewable energy structures,adjusting industrial structures to green businesses,and emphasizing fixed asset investment more environmentally friendly.Moreover,the result approved the EKC hypothesis.It shows an inverse U-shaped curve between GDP per capita and CO_(2)emissions in Indonesia.Indonesia’s fastest carbon emission peak is under scenario seven and is expected in 2040.Although an IDR 30 per kg CO_(2)e carbon tax and nuclear power will take decades to reduce carbon emissions,the carbon tax can still be a reference and has advantages to implement.This result can be a good beginning step for Indonesia,which has yet to gain experience with a carbon tax that can be implemented immediately and is helpful to decision-makers in putting into practice sensible measures to attain Indonesia’s carbon emission peaking.This research provides actionable insights internationally on carbon tax policies,nuclear energy adoption,EKC dynamics,global policy implications,and fostering international cooperation for carbon emission reductions.

Predicting the probability distribution of Martian rocks mechanical property based on microscale rock mechanical experiments and accurate grain-based modeling

The exploration of Mars would heavily rely on Martian rocks mechanics and engineering technology.As the mechanical property of Martian rocks is uncertain,it is of utmost importance to predict the probability distribution of Martian rocks mechanical property for the success of Mars exploration.In this paper,a fast and accurate probability distribution method for predicting the macroscale elastic modulus of Martian rocks was proposed by integrating the microscale rock mechanical experiments(micro-RME),accurate grain-based modeling(AGBM)and upscaling methods based on reliability principles.Firstly,the microstructure of NWA12564 Martian sample and elastic modulus of each mineral were obtained by micro-RME with TESCAN integrated mineral analyzer(TIMA)and nanoindentation.The best probability distribution function of the minerals was determined by Kolmogorov-Smirnov(K-S)test.Secondly,based on best distribution function of each mineral,the Monte Carlo simulations(MCS)and upscaling methods were implemented to obtain the probability distribution of upscaled elastic modulus.Thirdly,the correlation between the upscaled elastic modulus and macroscale elastic modulus obtained by AGBM was established.The accurate probability distribution of the macroscale elastic modulus was obtained by this correlation relationship.The proposed method can predict the probability distribution of Martian rocks mechanical property with any size and shape samples.

Background removal from global auroral images:Data-driven dayglow modeling

Global images of auroras obtained by cameras on spacecraft are a key tool for studying the near-Earth environment.However,the cameras are sensitive not only to auroral emissions produced by precipitating particles,but also to dayglow emissions produced by photoelectrons induced by sunlight.Nightglow emissions and scattered sunlight can contribute to the background signal.To fully utilize such images in space science,background contamination must be removed to isolate the auroral signal.Here we outline a data-driven approach to modeling the background intensity in multiple images by formulating linear inverse problems based on B-splines and spherical harmonics.The approach is robust,flexible,and iteratively deselects outliers,such as auroral emissions.The final model is smooth across the terminator and accounts for slow temporal variations and large-scale asymmetries in the dayglow.We demonstrate the model by using the three far ultraviolet cameras on the Imager for Magnetopause-to-Aurora Global Exploration(IMAGE)mission.The method can be applied to historical missions and is relevant for upcoming missions,such as the Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission.

Surrogate modeling for unsaturated infiltration via the physics and equality-constrained artificial neural networks

Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

Artificial intelligence-driven radiomics study in cancer:the role of feature engineering and modeling

Modern medicine is reliant on various medical imaging technologies for non-invasively observing patients’anatomy.However,the interpretation of medical images can be highly subjective and dependent on the expertise of clinicians.Moreover,some potentially useful quantitative information in medical images,especially that which is not visible to the naked eye,is often ignored during clinical practice.In contrast,radiomics performs high-throughput feature extraction from medical images,which enables quantitative analysis of medical images and prediction of various clinical endpoints.Studies have reported that radiomics exhibits promising performance in diagnosis and predicting treatment responses and prognosis,demonstrating its potential to be a non-invasive auxiliary tool for personalized medicine.However,radiomics remains in a developmental phase as numerous technical challenges have yet to be solved,especially in feature engineering and statistical modeling.In this review,we introduce the current utility of radiomics by summarizing research on its application in the diagnosis,prognosis,and prediction of treatment responses in patients with cancer.We focus on machine learning approaches,for feature extraction and selection during feature engineering and for imbalanced datasets and multi-modality fusion during statistical modeling.Furthermore,we introduce the stability,reproducibility,and interpretability of features,and the generalizability and interpretability of models.Finally,we offer possible solutions to current challenges in radiomics research.

Three-dimensional Modeling of Ore-forming Elements and Mineralization Prognosis for the Yechangping Mo Deposit,Henan Province,China

Three-dimensional geochemical modeling of ore-forming elements is crucial for predicting deep mineralization.This approach provides key information for the quantitative prediction of deep mineral localization,three-dimensional fine interpolation,analysis of spatial distribution patterns,and extraction of quantitative mineral-seeking markers.The Yechangping molybdenum(Mo)deposit is a significant and extensive porphyry-skarn deposit in the East Qinling-Dabie Mo polymetallic metallogenic belt at the southern margin of the North China Block.Abundant borehole data on oreforming elements underpin deep geochemical predictions.The methodology includes the following steps:(1)Threedimensional geological modeling of the deposit was established.(2)Correlation,cluster,and factor analyses post delineation of mineralization stages and determination of mineral generation sequence to identify(Cu,Pb,Zn,Ag)and(Mo,W,mfe)assemblages.(3)A three-dimensional geochemical block model was constructed for Mo,W,mfe,Cu,Zn,Pb,and Ag using the ordinary kriging method,and the variational function was developed.(4)Spatial distribution and enrichment characteristics analysis of ore-forming elements are performed to extract geological information,employing the variogram and w(Cu+Pb+Zn+Ag)/w(Mo+W)as predictive indicators.(5)Identifying the western,northwestern,and southwestern areas of the mine with limited mineralization potential,contrasted by the northeastern and southeastern areas favorable for mineral exploration.

Modeling time-dependent mechanical behavior of hard rock considering excavation-induced damage and complex 3D stress states

To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main components:a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state,a quantitative relationship between critical irreversible deformation and complex stress state,and evolution characteristics of strength parameters.The proposed model is implemented in a self-developed numerical code,i.e.CASRock.The reliability of the model is validated through experiments.It is indicated that the time-dependent fracturing potential index(xTFPI)at a given time during the attenuation creep stage shows a negative correlation with the extent of excavationinduced damage.The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress,thereby highlighting the 3D stress-dependent characteristic of the model.Finally,the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed.