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.

Numerical three-dimensional modeling of earthen dam piping failure

A physically-based numerical three-dimensional earthen dam piping failure model is developed for homogeneous and zoned soil dams.This model is an erosion model,coupled with force/moment equilibrium analyses.Orifice flow and two-dimensional(2D)shallow water equations(SWE)are solved to simulate dam break flows at different breaching stages.Erosion rates of different soils with different construction compaction efforts are calculated using corresponding erosion formulae.The dam's real shape,soil properties,and surrounding area are programmed.Large outer 2D-SWE grids are used to control upstream and downstream hydraulic conditions and control the boundary conditions of orifice flow,and inner 2D-SWE flow is used to scour soil and perform force/moment equilibrium analyses.This model is validated using the European Commission IMPACT(Investigation of Extreme Flood Processes and Uncertainty)Test#5 in Norway,Teton Dam failure in Idaho,USA,and Quail Creek Dike failure in Utah,USA.All calculated peak outflows are within 10%errors of observed values.Simulation results show that,for a V-shaped dam like Teton Dam,a piping breach location at the abutment tends to result in a smaller peak breach outflow than the piping breach location at the dam's center;and if Teton Dam had broken from its center for internal erosion,a peak outflow of 117851 m'/s,which is 81%larger than the peak outflow of 65120 m3/s released from its right abutment,would have been released from Teton Dam.A lower piping inlet elevation tends to cause a faster/earlier piping breach than a higher piping inlet elevation.

Numerical Models and Methods of Atmospheric Parameters Originating in the Formation of the Earth’s Climatic Cycle

Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model of key atmospheric parameters. The distribution of these parameters across the entire planet Earth is the origin of the formation of the climatic cycle, which is a normal climatic variation. To do this, the Earth is divided into eight (8) parts according to the number of key parameters to be defined in a physical representation of the model. Following this distribution, numerical models calculate the constants for the formation of water, vapor, ice, dryness, thermal energy (fire), heat, air, and humidity. These models vary in complexity depending on the indirect trigonometric direction and simplicity in the sum of neighboring models. Note that the constants obtained from the equations yield 275.156˚K (2.006˚C) for water, 273.1596˚K (0.00963˚C) for vapor, 273.1633˚K (0.0133˚C) for ice, 0.00365 in/s for atmospheric dryness, 1.996 in2/s for humidity, 2.993 in2/s for air, 1 J for thermal energy of fire, and 0.9963 J for heat. In summary, this study aims to define the main parameters and natural phenomena contributing to the modification of planetary climate. .

Cellular interplay to 3D in vitro microphysiological disease model:cell patterning microbiota-gut-brain axis

The microbiota-gut-brain axis(MGBA)has emerged as a key prospect in the bidirectional communication between two major organ systems:the brain and the gut.Homeostasis between the two organ systems allows the body to function without disease,whereas dysbiosis has long-standing evidence of etiopathological conditions.The most common communication paths are the microbial release of metabolites,soluble neurotransmitters,and immune cells.However,each pathway is intertwined with a complex one.With the emergence of in vitro models and the popularity of three-dimensional(3D)cultures and Transwells,engineering has become easier for the scientific understanding of neurodegenerative diseases.This paper briefly retraces the possible communication pathways between the gut microbiome and the brain.It further elaborates on three major diseases:autism spectrum disorder,Parkinson’s disease,and Alzheimer’s disease,which are prevalent in children and the elderly.These diseases also decrease patients’quality of life.Hence,understanding them more deeply with respect to current advances in in vitro modeling is crucial for understanding the diseases.Remodeling of MGBA in the laboratory uses many molecular technologies and biomaterial advances.Spheroids and organoids provide a more realistic picture of the cell and tissue structure than monolayers.Combining them with the Transwell system offers the advantage of compartmentalizing the two systems(apical and basal)while allowing physical and chemical cues between them.Cutting-edge technologies,such as bioprinting and microfluidic chips,might be the future of in vitro modeling,as they provide dynamicity.

Exploring mechanism of hidden,steep obliquely inclined bedding landslides using a 3DEC model:A case study of the Shanyang landslide in Shaanxi Province,China

Catastrophic geological disasters frequently occur on slopes with obliquely inclined bedding structures(also referred to as obliquely inclined bedding slopes),where the apparent dip sliding is not readily visible.This phenomenon has become a focal point in landslide research.Yet,there is a lack of studies on the failure modes and mechanisms of hidden,steep obliquely inclined bedding slopes.This study investigated the Shanyang landslide in Shaanxi Province,China.Using field investigations,laboratory tests of geotechnical parameters,and the 3DEC software,this study developed a numerical model of the landslide to analyze the failure process of such slopes.The findings indicate that the Shanyang landslide primarily crept along a weak interlayer under the action of gravity.The landslide,initially following a dip angle with the support of a stable inclined rock mass,shifted direction under the influence of argillization in the weak interlayer,moving towards the apparent dip angle.The slide resistance effect of the karstic dissolution zone was increasingly significant during this process,with lateral friction being the primary resistance force.A reduction in the lateral friction due to karstic dissolution made the apparent dip sliding characteristics of the Shanyang landslide more pronounced.Notably,deformations such as bending and uplift at the slope’s foot suggest that the main slide resistance shifts from lateral friction within the karstic dissolution zone to the slope foot’s resistance force,leading to the eventual buckling failure of the landslide.This study unveils a novel failure mode of apparent dip creep-buckling in the Shanyang landslide,highlighting the critical role of lateral friction from the karstic dissolution zone in its failure mechanism.These insights offer a valuable reference for mitigating risks and preventing disasters related to obliquely inclined bedding landslides.

3D bioprinting of in vitro porous hepatoma models:establishment,evaluation,and anticancer drug testing

Traditional tumor models do not tend to accurately simulate tumor growth in vitro or enable personalized treatment and are particularly unable to discover more beneficial targeted drugs.To address this,this study describes the use of threedimensional(3D)bioprinting technology to construct a 3D model with human hepatocarcinoma SMMC-7721 cells(3DP-7721)by combining gelatin methacrylate(GelMA)and poly(ethylene oxide)(PEO)as two immiscible aqueous phases to form a bioink and innovatively applying fluorescent carbon quantum dots for long-term tracking of cells.The GelMA(10%,mass fraction)and PEO(1.6%,mass fraction)hydrogel with 3:1 volume ratio offered distinct pore-forming characteristics,satisfactorymechanical properties,and biocompatibility for the creation of the 3DP-7721 model.Immunofluorescence analysis and quantitative real-time fluorescence polymerase chain reaction(PCR)were used to evaluate the biological properties of the model.Compared with the two-dimensional culture cell model(2D-7721)and the 3D mixed culture cell model(3DM-7721),3DP-7721 significantly improved the proliferation of cells and expression of tumor-related proteins and genes.Moreover,we evaluated the differences between the three culture models and the effectiveness of antitumor drugs in the three models and discovered that the efficacy of antitumor drugs varied because of significant differences in resistance proteins and genes between the three models.In addition,the comparison of tumor formation in the three models found that the cells cultured by the 3DP-7721 model had strong tumorigenicity in nude mice.Immunohistochemical evaluation of the levels of biochemical indicators related to the formation of solid tumors showed that the 3DP-7721 model group exhibited pathological characteristics of malignant tumors,the generated solid tumors were similar to actual tumors,and the deterioration was higher.This research therefore acts as a foundation for the application of 3DP-7721 models in drug development research.

3D Road Network Modeling and Road Structure Recognition in Internet of Vehicles

Internet of Vehicles (IoV) is a new system that enables individual vehicles to connect with nearby vehicles,people, transportation infrastructure, and networks, thereby realizing amore intelligent and efficient transportationsystem. The movement of vehicles and the three-dimensional (3D) nature of the road network cause the topologicalstructure of IoV to have the high space and time complexity.Network modeling and structure recognition for 3Droads can benefit the description of topological changes for IoV. This paper proposes a 3Dgeneral roadmodel basedon discrete points of roads obtained from GIS. First, the constraints imposed by 3D roads on moving vehicles areanalyzed. Then the effects of road curvature radius (Ra), longitudinal slope (Slo), and length (Len) on speed andacceleration are studied. Finally, a general 3D road network model based on road section features is established.This paper also presents intersection and road section recognition methods based on the structural features ofthe 3D road network model and the road features. Real GIS data from a specific region of Beijing is adopted tocreate the simulation scenario, and the simulation results validate the general 3D road network model and therecognitionmethod. Therefore, thiswork makes contributions to the field of intelligent transportation by providinga comprehensive approach tomodeling the 3Droad network and its topological changes in achieving efficient trafficflowand improved road safety.

Rainfall-triggered waste dump instability analysis based on surface 3D deformation in physical model test

Landslide is the second largest natural disaster after earthquake. It is of significance to study the evolution laws and failure mechanism of landslides based on its surface 3D deformation information. Based on the rainfall-triggered waste dump instability model test, we studied the failure mechanisms of the waste dump by integrating surface deformation and internal slope stress and proposed novel parameters for identifying landslide stability. We developed a noncontact measurement device, which can obtain millimeter-level 3D deformation data for surface scene in physical model test;Then we developed the similar materials and established a test model for a waste dump. Based on the failure characteristics of slope surface, internal stress of slope body and displacement contours during the whole process, we divided the slope instability process in model test into four stages: rainfall infiltration and surface erosion, shallow sliding, deep sliding, and overall instability. Based on the obtained surface deformation data, we calculated the volume change during slope instability process and compared it with the point displacement on slope surface. The results showed that the volume change can not only reflect the slow-ultra acceleration process of slope failure, but also fully reflect the above four stages and reduce the fluctuations caused by random factors. Finally, this paper proposed two stability identification parameters: the volume change rate above the slip surface and the relative velocity of volume change rate. According to the calculation of these two parameters in model test, they can be used for study the deformation and failure mechanism of slope stability.

Review of Artificial Intelligence for Oil and Gas Exploration: Convolutional Neural Network Approaches and the U-Net 3D Model

Deep learning, especially through convolutional neural networks (CNN) such as the U-Net 3D model, has revolutionized fault identification from seismic data, representing a significant leap over traditional methods. Our review traces the evolution of CNN, emphasizing the adaptation and capabilities of the U-Net 3D model in automating seismic fault delineation with unprecedented accuracy. We find: 1) The transition from basic neural networks to sophisticated CNN has enabled remarkable advancements in image recognition, which are directly applicable to analyzing seismic data. The U-Net 3D model, with its innovative architecture, exemplifies this progress by providing a method for detailed and accurate fault detection with reduced manual interpretation bias. 2) The U-Net 3D model has demonstrated its superiority over traditional fault identification methods in several key areas: it has enhanced interpretation accuracy, increased operational efficiency, and reduced the subjectivity of manual methods. 3) Despite these achievements, challenges such as the need for effective data preprocessing, acquisition of high-quality annotated datasets, and achieving model generalization across different geological conditions remain. Future research should therefore focus on developing more complex network architectures and innovative training strategies to refine fault identification performance further. Our findings confirm the transformative potential of deep learning, particularly CNN like the U-Net 3D model, in geosciences, advocating for its broader integration to revolutionize geological exploration and seismic analysis.

Dynamic Modelling of a Hybrid Variable Reluctance Machine Using the 3D Finite Element Method

This paper presents the work carried out to evaluate the dynamic performance of the Hybrid Variable Reluctance Motor (HVRM). The fourth-order Runge-Kutta integration algorithm was employed to solve the equations of the dynamic model, in conjunction with the three-dimensional finite element method. The 3D numerical data was calculated using Comsol Multiphysics, which accounts for the nonlinearity of the ferromagnetic material and the 3D nature of the HVRM. The outcomes of this study are precise and accurately predict the dynamic behaviour of the HVRM in terms of rotor position response, rotational speed and torque. The distinctive contribution of this work lies in the 3D numerical modelling of the HVRM and the subsequent evaluation and analysis of its dynamic operation. Analytical and numerical 2D studies are less resource-intensive and time-consuming, and are more straightforward and rapid to analyse and interpret. However, they are constrained in their capacity to examine spatial, volumetric interactions and intricate dynamics such as flux studies where three 3D effects cannot be disregarded, winding end effects and the configuration and positioning of the interposed permanent magnet.

Contribution to the Full 3D Finite Element Modelling of a Hybrid Stepping Motor with and without Current in the Coils

The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.

Effect of the mineral spatial distribution heterogeneity on the tensile strength of granite:Insights from PFC3D-GBM numerical analysis

The mechanical characteristics of crystalline rocks are affected by the heterogeneity of the spatial distribution of minerals.In this paper,a novel three-dimensional(3D)grain-based model(GBM)based on particle flow code(PFC),i.e.PFC3D-GBM,is proposed.This model can accomplish the grouping of mineral grains at the 3D scale and then filling them.Then,the effect of the position distribution,geometric size,and volume composite of mineral grains on the cracking behaviour and macroscopic properties of granite are examined by conducting Brazilian splitting tests.The numerical results show that when an external load is applied to a sample,force chains will form around each contact,and the orientation distribution of the force chains is uniform,which is independent of the external load level.Furthermore,the number of high-strength force chains is proportional to the external load level,and the main orientation distribution is consistent with the external loading direction.The main orientation of the cracks is vertical to that of the high-strength force chains.The geometric size of the mineral grains controls the mechanical behaviours.As the average grain size increases,the number of transgranular contacts with higher bonding strength in the region connecting both loading points increases.The number of high-strength force chains increases,leading to an increase in the stress concentration value required for the macroscopic failure of the sample.Due to the highest bonding strength,the generation of transgranular cracks in quartz requires a higher concentrated stress value.With increasing volume composition of quartz,the number of transgranular cracks in quartz distributed in the region connecting both loading points increases,which requires many high-strength force chains.The load level rises,leading to an increase in the tensile strength of the numerical sample.

Mathematical Wave Functions and 3D Finite Element Modelling of the Electron and Positron

The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric and Magnetic fields. Also, every moving particle has a De Broglie wavelength determined by its mass and velocity. This paper shows that all of these properties of a particle can be derived from a single wave function equation for that particle. Wave functions for the Electron and the Positron are presented and principles are provided that can be used to calculate the wave functions of all the fundamental particles in Physics. Fundamental particles such as electrons and positrons are considered to be point particles in the Standard Model of Physics and are not considered to have a structure. This paper demonstrates that they do indeed have structure and that this structure extends into the space around the particle’s center (in fact, they have infinite extent), but with rapidly diminishing energy density with the distance from that center. The particles are formed from Electromagnetic standing waves, which are stable solutions to the Schrödinger and Classical wave equations. This stable structure therefore accounts for both the wave and particle nature of these particles. In fact, all of their properties such as mass, spin and electric charge, can be accounted for from this structure. These particle properties appear to originate from a single point at the center of the wave function structure, in the same sort of way that the Shell theorem of gravity causes the gravity of a body to appear to all originate from a central point. This paper represents the first two fully characterized fundamental particles, with a complete description of their structure and properties, built up from the underlying Electromagnetic waves that comprise these and all fundamental particles.

Image processing based three-dimensional model reconstruction for cross-platform numerical simulation

Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.

Numerical study on residual current and its impact on mass transport in the Hangzhou Bay and the Changjiang Estuary I. A3-D joint model of the Hangzhou Bay and the Changjiang Estuary

A 3-D numerical model is set up in a large domain covering the Hangzhou Bay and the Changjiang Estuary based on the ECOM model in orthogonal curvilinear coordinates.The numerical schemes for baroclinic pressure gradient (BPG)terms and convective terms are improved in the paper according to the characteristics of velocity field and mass transport in the area.The model is validated by the simulations of residual current and salinity transport in the Hangzhou Bay and the Changjiang Estuary.

3D打印实物模型在颈椎疾患住院医师规范化培训教学中的应用和探讨

目的探讨3 D打印实物模型在提高颈椎疾患住院医师规范化培训(简称住培)教学效果中的作用。方法选取参加骨科住培的颈椎疾患学习的住院医师80名为研究对象,随机分为两组,观察组(40人)以幻灯讲解加3 D打印实物模型教学,对照组(40人)采用幻灯讲解教学。课后即时和课后1月分别对住院医师的知识掌握情况和诊疗技术进行现场测试,利用问卷的形式对住院医师的满意度和参与兴趣(10分法)进行调查。结果课后即时观察组的知识掌握情况和对照组相当(P>0.05),但是观察组的诊疗技术测试成绩显著高于对照组(P<0.05),观察组的满意度和参与兴趣评分均显著高于对照组,两者对比差异有统计学意义(均P<0.05);课后1月时,观察组的知识掌握情况和诊疗技术测试成绩均显著高于对照组(均P<0.05)。结论3 D打印实物模型应用于颈椎疾患住培的带教中,可显著提高教学效果,提升住院医师的满意度和参与兴趣,值得在住培教学中推广。

Numerical modeling of 3D DC resistivity method in the mixed space-wavenumber domain

Forward modeling is the basis of inversion imaging and quantitative interpretation for DC resistivity exploration.Currently,a numerical model of the DC resistivity method must be finely divided to obtain a highly accurate solution under complex conditions,resulting in a long calculation time and large storage.Therefore,we propose a 3D numerical simulation method in a mixed space-wavenumber domain to overcome this challenge.The partial differential equation about abnormal potential is transformed into many independent ordinary differential equations with different wavenumbers using a 2D Fourier transform along the x axis and y axis direction.In this way,a large-scale 3D numerical simulation problem is decomposed into several 1D numerical simulation problems,which significantly reduces the computational and storage requirements.In addition,these ordinary 1D differential equations with different wavenumbers are independent of each other and high parallelelism of the algorithm.They are solved using a finite-element algorithm combined with a chasing method,and the obtained solution is modified using a contraction operator.In this method,the vertical direction is reserved as the spatial domain,then grid size can be determined flexibly based on the underground current density distribution,which considers the solution accuracy and calculation efficiency.In addition,for the first time,we use the contraction operator in the integral equation method to iterate the algorithm.The algorithm takes advantage of the high efficiency of the standard Fourier transform and chasing method,as well as the fast convergence of the contraction operator.We verified the accuracy of the algorithm and the convergence of the contraction operator.Compared with a volume integral method and goal-oriented adaptive finite-element method,the proposed algorithm has lower memory requirements and high computational efficiency,making it suitable for calculating a model with large-scale nodes.Moreover,different examples are used to verify the high adaptability and parallelism of the proposed algorithm.The findings show that the 3D numerical simulation method of DC resistivity method in a mixed space-wavenumber domain is highly efficient,precise,and parallel.

Numerical Simulation of Injection Molding Cooling Process Based on 3D Surface Model

The design of the cooling system of injection molds directly affects both productivity and the quality of the final part. Using the cooling process CAE system to instruct the mold design, the efficiency and quality of design can be improved greatly. At the same time, it is helpful to confirm the cooling system structure and optimize the process conditions. In this paper, the 3D surface model of mold cavity is used to replace the middle-plane model in the simulation by Boundary Element Method, which break the bottleneck of the application of the injection molding simulation softwares base on the middle-plane model. With the improvements of this paper, a practical and commercial simulation software of injection molding cooling process named as HsCAE3D6.0 is developed.

On numerical modeling of low-head direct chill ingot caster for magnesium alloy AZ31

A comprehensive 3D turbulent CFD study has been carried out to simulate a Low-Head(LH)vertical Direct Chill(DC)rolling ingot caster for the common magnesium alloy AZ31.The model used in this study takes into account the coupled laminar/turbulent melt flow and solidification aspects of the process and is based on the control-volume finite-difference approach.Following the aluminum/magnesium DC casting industrial practices,the LH mold is taken as 30 mm with a hot top of 60 mm.The previously verified in-house code has been modified to model the present casting process.Important quantitative results are obtained for four casting speeds,for three inlet melt pouring temperatures(superheats)and for three metal-mold contact heat transfer coefficients for the steady state operational phase of the caster.The variable cooling water temperatures reported by the industry are considered for the primary and secondary cooling zones during the simulations.Specifically,the temperature and velocity fields,sump depth and sump profiles,mushy region thickness,solid shell thickness at the exit of the mold and axial temperature profiles at the center and at three strategic locations at the surface of the slab are presented and discussed.