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.

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.

Development and Application of a Power Law Constitutive Model for Eddy Current Dampers

Eddy current dampers (ECDs) have emerged as highly desirable solutions for vibration control due to theirexceptional damping performance and durability. However, the existing constitutive models present challenges tothe widespread implementation of ECD technology, and there is limited availability of finite element analysis (FEA)software capable of accurately modeling the behavior of ECDs. This study addresses these issues by developing anewconstitutivemodel that is both easily understandable and user-friendly for FEAsoftware. By utilizing numericalresults obtained from electromagnetic FEA, a novel power law constitutive model is proposed to capture thenonlinear behavior of ECDs. The effectiveness of the power law constitutive model is validated throughmechanicalproperty tests and numerical seismic analysis. Furthermore, a detailed description of the application process ofthe power law constitutive model in ANSYS FEA software is provided. To facilitate the preliminary design ofECDs, an analytical derivation of energy dissipation and parameter optimization for ECDs under harmonicmotionis performed. The results demonstrate that the power law constitutive model serves as a viable alternative forconducting dynamic analysis using FEA and optimizing parameters for ECDs.

Novel damage constitutive models and new quantitative identification method for stress thresholds of rocks under uniaxial compression

Four key stress thresholds exist in the compression process of rocks,i.e.,crack closure stress(σ_(cc)),crack initiation stress(σ_(ci)),crack damage stress(σ_(cd))and compressive strength(σ_(c)).The quantitative identifications of the first three stress thresholds are of great significance for characterizing the microcrack growth and damage evolution of rocks under compression.In this paper,a new method based on damage constitutive model is proposed to quantitatively measure the stress thresholds of rocks.Firstly,two different damage constitutive models were constructed based on acoustic emission(AE)counts and Weibull distribution function considering the compaction stages of the rock and the bearing capacity of the damage element.Then,the accumulative AE counts method(ACLM),AE count rate method(CRM)and constitutive model method(CMM)were introduced to determine the stress thresholds of rocks.Finally,the stress thresholds of 9 different rocks were identified by ACLM,CRM,and CMM.The results show that the theoretical stress−strain curves obtained from the two damage constitutive models are in good agreement with that of the experimental data,and the differences between the two damage constitutive models mainly come from the evolutionary differences of the damage variables.The results of the stress thresholds identified by the CMM are in good agreement with those identified by the AE methods,i.e.,ACLM and CRM.Therefore,the proposed CMM can be used to determine the stress thresholds of rocks.

Numerical analysis of high‑speed railway slab tracks using calibrated and validated 3D time‑domain modelling

Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a novel slab trackform for high-speed railways is investigated using three-dimensional finite element modelling in Abaqus.It is then compared to the performance of a ballasted track.First,slab and ballasted track models are developed to replicate the full-scale testing of track sections.Once the models are calibrated with the experimental results,the novel slab model is developed and compared against the calibrated slab track results.The slab and ballasted track models are then extended to create linear dynamic models,considering the track geodynamics,and simulating train passages at various speeds,for which the Ledsgard documented case was used to validate the models.Trains travelling at low and high speeds are analysed to investigate the track deflections and the wave propagation in the soil,considering the issues associated with critical speeds.Various train loading methods are discussed,and the most practical approach is retained and described.Moreover,correlations are made between the geotechnical parameters of modern high-speed rail and conventional standards.It is found that considering the same ground condition,the slab track deflections are considerably smaller than those of the ballasted track at high speeds,while they show similar behaviour at low speeds.

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.

Nonlinear constitutive models of rock structural plane and their applications

Structural planes play an important role in controlling the stability of rock engineering,and the influence of structural planes should be considered in the design and construction process of rock engineering.In this paper,mechanical properties,constitutive theory,and numerical application of structural plane are studied by a combination method of laboratory tests,theoretical derivation,and program development.The test results reveal the change laws of various mechanical parameters under different roughness and normal stress.At the pre-peak stage,a non-stationary model of shear stiffness is established,and threedimensional empirical prediction models for initial shear stiffness and residual stage roughness are proposed.The nonlinear constitutive models are established based on elasto-plastic mechanics,and the algorithms of the models are developed based on the return mapping algorithm.According to a large number of statistical analysis results,empirical prediction models are proposed for model parameters expressed by structural plane characteristic parameters.Finally,the discrete element method(DEM)is chosen to embed the constitutive models for practical application.The running programs of the constitutive models have been compiled into the discrete element model library.The comparison results between the proposed model and the Mohr-Coulomb slip model show that the proposed model can better describe nonlinear changes at different stages,and the predicted shear strength,peak strain and shear stiffness are closer to the test results.The research results of the paper are conducive to the accurate evaluation of structural plane in rock engineering.

Systematic Elastostatic Stiffness Model of Over-Constrained Parallel Manipulators Without Additional Constraint Equations

The establishment of an elastostatic stiffness model for over constrained parallel manipulators(PMs),particularly those with over constrained subclosed loops,poses a challenge while ensuring numerical stability.This study addresses this issue by proposing a systematic elastostatic stiffness model based on matrix structural analysis(MSA)and independent displacement coordinates(IDCs)extraction techniques.To begin,the closed-loop PM is transformed into an open-loop PM by eliminating constraints.A subassembly element is then introduced,which considers the flexibility of both rods and joints.This approach helps circumvent the numerical instability typically encountered with traditional constraint equations.The IDCs and analytical constraint equations of nodes constrained by various joints are summarized in the appendix,utilizing multipoint constraint theory and singularity analysis,all unified within a single coordinate frame.Subsequently,the open-loop mechanism is efficiently closed by referencing the constraint equations presented in the appendix,alongside its elastostatic model.The proposed method proves to be both modeling and computationally efficient due to the comprehensive summary of the constraint equations in the Appendix,eliminating the need for additional equations.An example utilizing an over constrained subclosed loops demonstrate the application of the proposed method.In conclusion,the model proposed in this study enriches the theory of elastostatic stiffness modeling of PMs and provides an effective solution for stiffness modeling challenges they present.

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.

A framework for dynamic modelling of railway track switches considering the switch blades,actuators and control systems

The main contribution of this paper is the development and demonstration of a novel methodology that can be followed to develop a simulation twin of a railway track switch system to test the functionality in a digital environment.This is important because,globally,railway track switches are used to allow trains to change routes;they are a key part of all railway networks.However,because track switches are single points of failure and safety-critical,their inability to operate correctly can cause significant delays and concomitant costs.In order to better understand the dynamic behaviour of switches during operation,this paper has developed a full simulation twin of a complete track switch system.The approach fuses finite element for the rail bending and motion,with physics-based models of the electromechanical actuator system and the control system.Hence,it provides researchers and engineers the opportunity to explore and understand the design space around the dynamic operation of new switches and switch machines before they are built.This is useful for looking at the modification or monitoring of existing switches,and it becomes even more important when new switch concepts are being considered and evaluated.The simulation is capable of running in real time or faster meaning designs can be iterated and checked interactively.The paper describes the modelling approach,demonstrates the methodology by developing the system model for a novel“REPOINT”switch system,and evaluates the system level performance against the dynamic performance requirements for the switch.In the context of that case study,it is found that the proposed new actuation system as designed can meet(and exceed)the system performance requirements,and that the fault tolerance built into the actuation ensures continued operation after a single actuator failure.

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.

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.

Development and Validation of a Scaled Electric Combat Vehicle Tire Model

Pneumatic tire modeling and validation have been the topic of several research papers, however, most of these papers only deal with pneumatic passenger and truck tires. In recent years, wheeled-scaled vehicles have gained lots of attention as a feasible testing platform, nonetheless up to the authors’ knowledge there has been no research regarding the use of scaled tires and their effect on the overall vehicle performance characteristics. This paper presents a novel scaled electric combat vehicle tire model and validation technique. The pro-line lockdown tire size 3.00 × 7.35 is modeled using the Finite Element Analysis (FEA) technique and several materials including layered membrane, beam elements, and Mooney-Rivlin for rubber. The tire-rim assembly is then described, and the rigid body analysis is presented. The tire is then validated using an in-house custom-made static tire testing machine. The tire test rig is made specifically to test the pro-line tire model and is designed and manufactured in the laboratory. The tire is validated using vertical stiffness and footprint tests in the static domain at different operating conditions including several vertical loads. Then the tire is used to perform rolling resistance and steering analysis including the rolling resistance coefficient and the cornering stiffness. The analysis is performed at different operating conditions including longitudinal speeds of 5, 10, and 15 km/h. This tire model will be further used to determine the tractive and braking performance of the tire. Furthermore, the tire test rig will also be modified to perform cornering stiffness tests.

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.

A Modified Model for Soil–Structure Interface Considering Random Damage of Mesoscopic Contact Elements

The interaction between soil and marine structures like submarine pipeline/pipe pile/suction caisson is a complicated geotechnical mechanism process.In this study,the interface is discretized into multiple mesoscopic contact elements that are damaged randomly throughout the shearing process due to the natural heterogeneity.The evolution equation of damage variable is developed based on the Weibull function,which is able to cover a rather wide range of distribution shapes by only two parameters,making it applicable for varying scenarios.Accordingly,a statistical damage model is established by incorporating Mohr–Coulomb strength criterion,in which the interfacial residual strength is considered whereby the strain softening behavior can be described.A concept of“semi-softening”characteristic point on shear stress–displacement curve is proposed for effectively modeling the evolution of strain softening.Finally,a series of ring shear tests of the interfaces between fine sea sand and smooth/rough steel surfaces are conducted.The predicted results using the proposed model are compared with experimental data of this study as well as some results from existing literature,indicating that the model has a good performance in modeling the progressive failure and strain softening behavior for various types of soil–structure interfaces.

Constitutive modelling of fabric effect on sand liquefaction

Sand liquefaction under static and dynamic loading can cause failure of embankments,slopes,bridges and other important infrastructure.Sand liquefaction in the seabed can also cause submarine landslides and tsunamis.Fabric anisotropy related to the internal soil structure such as particle orientation,force network and void space is found to have profound influence on sand liquefaction.A constitutive model accounting for the effect of anisotropy on sand liquefaction is proposed.Evolution of fabric anisotropy during loading is considered according to the anisotropic critical state theory for sand.The model has been validated by extensive test results on Toyoura sand with different initial densities and stress states.The effect of sample preparation method on sand liquefaction is qualitatively analysed.The model has been used to investigate the response of a sand ground under earthquake loading.It is shown that sand with horizontal bedding plane has the highest resistance to liquefaction when the sand deposit is anisotropic,which is consistent with the centrifuge test results.The initial degree of fabric anisotropy has a more significant influence on the liquefaction resistance.Sand with more anisotropic fabric that can be caused by previous loading history or compaction methods has lower liquefaction resistance.

Numerical Prediction of Ride Comfort of Tracked Vehicle Equipped with Novel Flexible Road Wheels

Enhancing ride comfort has always constituted a crucial focus in the design and research of modern tracked vehicles,heavily reliant on the driving system's performance.While the road wheel is a key component of the driving system,traditional road wheels predominantly adopt a solid structure,exhibiting subpar adhesion performance and damping effects,thereby falling short of meeting the demands for high-speed,stable,and long-distance driving in tracked vehicles.Addressing this issue,this paper proposes a novel type of flexible road wheel(FRW)characterized by a catenary construction.The study investigates the ride comfort of tracked vehicles equipped with flexible road wheels by integrating finite element and vehicle dynamic.First,three-dimensional(3D)finite element(FE)models of both flexible and rigid road wheels are established,considering material and contact nonlinearities.These models are validated through a wheel radial loading test.Based on the validated FE model,the paper uncovers the relationship between load and radial deformation of the road wheel,forming the basis for a nonlinear mathematical model.Subsequently,a half-car model of a tracked vehicle with seven degrees of freedom is established using Newton's second law.A random road model,considering the track effect and employing white noise,is constructed.The study concludes by examining the ride comfort of tracked vehicles equipped with flexible and rigid road wheels under various speeds and road grades.The results demonstrate that,in comparison to the rigid road wheel(RRW),the flexible road wheel enhances the ride comfort of tracked vehicles on randomly uneven roads.This research provides a theoretical foundation for the implementation of flexible road wheels in tracked vehicles.

Physics-based seismic analysis of ancient wood structure:fault-to-structure simulation

Based on the domain reduction method,this study employs an SEM-FEM hybrid workflow which integrates the advantages of the spectral element method(SEM)for flexible and highly efficient simulation of seismic wave propagation in a three-dimensional(3D)regional-scale geophysics model and the finite element method(FEM)for fine simulation of structural response including soil-structure interaction,and performs a physics-based simulation from initial fault rupture on an ancient wood structure.After verification of the hybrid workflow,a large-scale model of an ancient wood structure in the Beijing area,The Tower of Buddhist Incense,is established and its responses under the 1665 Tongxian earthquake and the 1730 Yiheyuan earthquake are simulated.The results from the simulated ground motion and seismic response of the wood structure under the two earthquakes demonstrate that this hybrid workflow can be employed to efficiently provide insight into the relationships between geophysical parameters and the structural response,and is of great significance toward accurate input for seismic simulation of structures under specific site and fault conditions.