Modularized and Parametric Modeling Technology for Finite Element Simulations of Underground Engineering under Complicated Geological Conditions

The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.

Intelligent architecture modeling for multilevel fluvial reservoirs based on architecture interface and elements

At present,the architecture modeling method of fluvial reservoirs are still developing.Traditional methods usually use grids to characterize architecture interbeds within the reservoir.Due to the thin thickness of this type of the interlayers,the number of the model grids must be greatly expanded.The number of grids in the tens of millions often makes an expensive computation;however,upscaling the model will generate a misleading model.The above confusion is the major reason that restricts the largescale industrialization of fluvial reservoir architecture models in oilfield development and production.Therefore,this paper explores an intelligent architecture modeling method for multilevel fluvial reservoirs based on architecture interface and element.Based on the superpositional relationship of different architectural elements within the fluvial reservoir,this method uses a combination of multilevel interface constraints and non-uniform grid techniques to build a high-resolution 3D geological model for reservoir architecture.Through the grid upscaling technology of heterogeneous architecture elements,different upscaling densities are given to the lateral-accretion bedding and lateral-accretion bodies to simplify the model gridding.This new method greatly reduces the number of model grids while ensuring the accuracy of lateral-accretion bedding models,laying a foundation for large-scale numerical simulation of the subsequent industrialization of the architecture model.This method has been validated in A layer of X oilfield with meandering fluvial channel sands as reservoirs and B layer of Y oilfield with braided river sands as reservoirs.The simulation results show that it has a higher accuracy of production history matching and remaining oil distribution forecast of the targeted sand body.The numerical simulation results show that in the actual development process of oilfield,the injected water will not displace oil in a uniform diffusive manner as traditionally assumed,but in a more complex pattern with oil in upper part of sand body being left behind as residual oil due to the influences of different levels of architecture interfaces.This investigation is important to guiding reservoir evaluation,remaining oil analysis,profile control and potential tapping and well pattern adjustment.

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.

In silico optimization of actuation performance in dielectric elastomercomposites via integrated finite element modeling and deep learning

Dielectric elastomers(DEs)require balanced electric actuation performance and mechanical integrity under applied voltages.Incorporating high dielectric particles as fillers provides extensive design space to optimize concentration,morphology,and distribution for improved actuation performance and material modulus.This study presents an integrated framework combining finite element modeling(FEM)and deep learning to optimize the microstructure of DE composites.FEM first calculates actuation performance and the effective modulus across varied filler combinations,with these data used to train a convolutional neural network(CNN).Integrating the CNN into a multi-objective genetic algorithm generates designs with enhanced actuation performance and material modulus compared to the conventional optimization approach based on FEM approach within the same time.This framework harnesses artificial intelligence to navigate vast design possibilities,enabling optimized microstructures for high-performance DE composites.

Multi-scale Modeling and Finite Element Analyses of Thermal Conductivity of 3D C/SiC Composites Fabricating by Flexible-Oriented Woven Process

Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.

Discrete Element Modelling of Damage Evolution of Concrete Considering Meso-Structure of ITZ

The mechanical properties of interfacial transition zones(ITZs)have traditionally been simplified by reducing the stiffness of cement in previous simulation methods.A novel approach based on the discrete element method(DEM)has been developed for modeling concrete.This new approach efficiently simulates the meso-structure of ITZs,accurately capturing their heterogeneous properties.Validation against established uniaxial compression experiments confirms the precision of thismodel.The proposedmodel canmodel the process of damage evolution containing cracks initiation,propagation and penetration.Under increasing loads,cracks within ITZs progressively accumulate,culminating in macroscopic fractures that traverse themortarmatrix,forming the complex,serpentine path of cracks.This study reveals four distinct displacement patterns:tensile compliant,tensile opposite,mixed tensile-shear,and shear opposite patterns,each indicative of different stages in concrete’s damage evolution.The widening angle of these patterns delineates the progression of cracks,with the tensile compliant pattern signaling the initial crack appearance and the shear opposite pattern indicating the concrete model’s ultimate failure.

Modeling asymmetric fracture mechanics of Mg alloy wire in drawing process

In this study,a numerical analysis was conducted on the ductile fracture of a 2-mm diameter Mg-1Zn-0.5Mn-0.5Sr-0.1Ca alloy wire during drawing.The hexagonally close-packed crystal structure of Mg alloys causes asymmetric fracture behavior,especially in the compression region.The aim of this study is to develop a comprehensive damage model for Mg alloy wire that accurately predicts ductile fracture,with a focus on the compression region.A novel experimental method was introduced to measure the ductile fracture of Mg alloy wires under different stress states.The wire drawing process was simulated using the Generalized Incremental Stress-State dependent damage(GISSMO)Model and the Semi-Analytical Model for Polymers(SAMP)model.The damage model's prediction and the experimental results were found to be in excellent agreement,especially in determining crack initiation.Computational analysis established a safe zone diagram for die angle and reduction ratio,and experimental validation confirmed the feasibility of this approach.The proposed damage model can provide a practical and reliable analysis for optimizing the drawing process of Mg alloy wire.

Influence of Various Earth-Retaining Walls on the Dynamic Response Comparison Based on 3D Modeling

The present work aims to assess earthquake-induced earth-retaining(ER)wall displacement.This study is on the dynamics analysis of various earth-retaining wall designs in hollow precast concrete panels,reinforcement concrete facing panels,and gravity-type earth-retaining walls.The finite element(FE)simulations utilized a 3D plane strain condition to model full-scale ER walls and numerous nonlinear dynamics analyses.The seismic performance of differentmodels,which includes reinforcement concrete panels and gravity-type and hollowprecast concrete ER walls,was simulated and examined using the FE approach.It also displays comparative studies such as stress distribution,deflection of the wall,acceleration across the wall height,lateral wall displacement,lateral wall pressure,and backfill plastic strain.Three components of the created ER walls were found throughout this research procedure.One is a granular reinforcement backfill,while the other is a wall-facing panel and base foundation.The dynamic response effects of varied earth-retaining walls have also been studied.It was discovered that the facing panel of the model significantly impacts the earthquake-induced displacement of ER walls.The proposed analytical model’s validity has been evaluated and compared with the reinforcement concrete facing panels,gravity-type ER wall,scientifically available data,and American Association of State Highway and Transportation Officials(AASHTO)guidelines results based on FE simulation.The results of the observations indicate that the hollow prefabricated concrete ER wall is the most feasible option due to its lower displacement and high-stress distribution compared to the two types.The methodology and results of this study establish standards for future analogous investigations and professionals,particularly in light of the increasing computational capabilities of desktop computers.

Implementation of a particle-in-cell method for the energy solver in 3D spherical geodynamic modeling

The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.

An ultrasound-guided percutaneous electrical nerve stimulation regimen devised using finite element modeling promotes functional recovery after median nerve transection

Percutaneous electrical nerve stimulation of an injured nerve can promote and accelerate peripheral nerve regeneration and improve function.When performing acupuncture and moxibustion,locating the injured nerve using ultrasound before percutaneous nerve stimulation can help prevent further injury to an already injured nerve.However,stimulation parameters have not been standardized.In this study,we constructed a multi-layer human forearm model using finite element modeling.Taking current density and activated function as optimization indicators,the optimal percutaneous nerve stimulation parameters were established.The optimal parameters were parallel placement located 3 cm apart with the injury site at the midpoint between the needles.To validate the efficacy of this regimen,we performed a randomized controlled trial in 23 patients with median nerve transection who underwent neurorrhaphy.Patients who received conventional rehabilitation combined with percutaneous electrical nerve stimulation experienced greater improvement in sensory function,motor function,and grip strength than those who received conventional rehabilitation combined with transcutaneous electrical nerve stimulation.These findings suggest that the percutaneous electrical nerve stimulation regimen established in this study can improve global median nerve function in patients with median nerve transection.

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.

Phase-field modeling of dendritic growth under forced flow based on adaptive finite element method

A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.

AUTOMATIC FINIT EELEMENT MODELING OF A WING

This paper is concerned about the automatic finite element modeling of a wing structure. The row and column method is used to identify the structure parts(ribs, spars, skins and pillars). A customization module of PCL(PATRAN Command Language under PATRAN 6.0) code from constructing airfoil curves to creating the entire wing FEM model is designed and developed. The geome tric, mesh density, material, load and boundary parameters can be easily and correctly input with the friendly interactive interface. A VFW614 wing is analyzed from creating airfoil curves to the show of stresses calculated by using NASTRAN 68 as an example. The results show that this customization module is very effective and efficient.

Discrete element modeling of migration and evolution rules of coarse aggregates in the static compaction process

To investigate migration and evolution rules of coarse aggregates in the static compaction process, an algorithm of generating digital coarse aggregates that can reflect real morphology（ such as shape, size and fracture surface） of aggregate particles, is represented by polyhedral particles based on the discrete element method（ DEM）. A digital specimen comprised of aggregates and air voids is developed. In addition,a static compaction model consisting of a digital specimen and three plates is constructed and a series of evaluation indices such as mean contact force σMCF, wall stress in direction of zcoordinate σWSZZ, porosity and coordination numbers are presented to investigate the motion rules of coarse aggregates at different compaction displacements of 7. 5, 15 and 30 mm. The three-dimensional static compaction model is also verified with laboratory measurements. The results indicate that the compaction displacements are positively related to σMCF and σWSZZ, which increase gradually with the increase in iterative steps. When the compaction proceeds, the digital specimen porosity decreases, but the coordination number increases. The variation ranges of these four indices are different at different compaction displacements. This study provides a method to analyze the compaction mechanism of particle materials such as asphalt mixture and graded broken stone.

Study on hot deformation behavior of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy using a combination of strain-compensated Arrhenius constitutive model and finite element simulation method

Isothermal hot compression experiments were conducted on homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy to investigate hot deformation behavior at the temperature range of 673-773 K and the strain rate range of 0.001-1 s^(-1)by using a Gleeble-1500D thermo mechanical simulator.Metallographic characterization on samples deformed to true strain of 0.70 illustrates the occurrence of flow localization and/or microcrack at deformation conditions of 673 K/0.01 s^(-1),673 K/1 s^(-1)and 698 K/1 s^(-1),indicating that these three deformation conditions should be excluded during hot working of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.Based on the measured true stress-strain data,the strain-compensated Arrhenius constitutive model was constructed and then incorporated into UHARD subroutine of ABAQUS software to study hot deformation process of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.By comparison with measured force-displacement curves,the predicted results can describe well the rheological behavior of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy,verifying the validity of finite element simulation of hot compression process with this complicated constitutive model.Numerical results demonstrate that the distribution of values of material parameters(α,n,Q and ln A)within deformed sample is inhomogeneous.This issue is directly correlated to the uneven distribution of equivalent plastic strain due to the friction effect.Moreover,at a given temperature the increase of strain rate would result in the decrease of equivalent plastic strain within the central region of deformed sample,which hinders the occurrence of dynamic recrystallization(DRX).

Finite Element Simulation of Radial Tire Building and Shaping Processes Using an Elasto-Viscoplastic Model

The comprehensive tire building and shaping processes are investigated through the finite element method(FEM)in this article.The mechanical properties of the uncured rubber from different tire components are investigated through cyclic loading-unloading experiments under different strain rates.Based on the experiments,an elastoviscoplastic constitutive model is adopted to describe themechanical behaviors of the uncured rubber.The distinct mechanical properties,including the stress level,hysteresis and residual strain,of the uncured rubber can all be well characterized.The whole tire building process(including component winding,rubber bladder inflation,component stitching and carcass band folding-back)and the shaping process are simulated using this constitutive model.The simulated green tire profile is in good agreement with the actual profile obtained through 3D scanning.The deformation and stress of the rubber components and the cord reinforcements during production can be obtained fromthe FE simulation,which is helpful for judging the rationality of the tire construction design.Finally,the influence of the parameter“drum width”is investigated,and the simulated result is found to be consistent with the experimental observations,which verifies the effectiveness of the simulation.The established simulation strategy provides some guiding significance for the improvement of tire design parameters and the elimination of tire production defects.

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.

Finite element modeling of pore-fluid flow in the Dachang ore district,Guangxi,China:Implications for hydrothermal mineralization

Convective heat transfer associated with the circulation of porefluid in porous rocks and fractures within the upper crust of the Earth is substantial when the temperature gradient is sufficiently high. In order to understand the process of Snpolymetallic mineralization in the Dachang ore district of Guangxi, a finite element method has been used in this study to simulate both pore-fluid flow and heat transfer in this district. On the basis of related geological, tectonic and geophysical constraints, a computational model was established. It enables a computational simulation and sensitivity analysis to be carried out for investigating ore-forming pore-fluid flow and other key factors that may affect hydrothermal ore genesis in the district. The related simulation results have indicated that： （1） permeable fault zones in the Dacbang ore district can serve as preferential pathways for pore-fluid flow on a regional-scale; and （2） the pore-fluid flow can affect the salinity distribution. This latter factor is part of the reason why Sn-polymetallic mineralization has taken place in this district.

Finite element modeling assumptions impact on seismic response demands of MRF-buildings

Recent seismic events have raised concerns over the safety and vulnerability of reinforced concrete moment resisting frame ＂RC-MRF＂ buildings. The seismic response of such buildings is greatly dependent on the computational tools used and the inherent assumptions in the modelling process. Thus, it is essential to investigate the sensitivity of the response demands to the corresponding modelling assumption. Many parameters and assumptions are justified to generate effective structural finite element（FE） models of buildings to simulate lateral behaviour and evaluate seismic design demands. As such, the present study focuses on the development of reliable FE models with various levels of refinement. The effects of the FE modelling assumptions on the seismic response demands on the design of buildings are investigated. the predictive ability of a FE model is tied to the accuracy of numerical analysis; a numerical analysis is performed for a series of symmetric buildings in active seismic zones. The results of the seismic response demands are presented in a comparative format to confirm drift and strength limits requirements. A proposed model is formulated based on a simplified modeling approach, where the most refined model is used to calibrate the simplified model.

2.5-D modeling of cross-hole electromagnetic measurement by nite element method

A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number （Ky） domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell＇s equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in