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.

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.

In vivo imaging and computational modeling of nonlinear shear waves in living skeletal muscles

How the state of living muscles modulates the features of nonlinear elastic waves generated by external dynamic loads remains unclear because of the challenge of directly observing and modeling nonlinear elastic waves in skeletal muscles in vivo,considering their active deformation behavior.Here,this important issue is addressed by combining experiments performed with an ultrafast ultrasound imaging system to track nonlinear shear waves(shear shock waves)in muscles in vivo and finite element analysis relying on a physically motivated constitutive model to study the effect of muscle activation level.Skeletal muscle was loaded with a deep muscle stimulator to generate shear shock waves(SSWs).The particle velocities,second and third harmonics,and group velocities of the SSWs in living muscles under both passive and active states were measured in vivo.Our experimental results reveal,for the first time,that muscle states have a pronounced effect on wave features;a low level of activation may facilitate the occurrence of both the second and third harmonics,whereas a high level of activation may inhibit the third harmonic.Finite element analysis was further carried out to quantitatively explore the effect of active muscle deformation behavior on the generation and propagation of SSWs.The simulation results at different muscle activation levels confirmed the experimental findings.The ability to reveal the effects of muscle state on the features of SSWs may be helpful in elucidating the unique dynamic deformation mechanism of living skeletal muscles,quantitatively characterizing diverse shock wave-based therapy instruments,and guiding the design of muscle-mimicking soft materials.

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.

Evaluating the mechanical properties of anisotropic shale containing bedding and natural fractures with discrete element modeling

Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.

Application of 1D/3D finite elements coupling for structural nonlinear analysis

An element coupling model (ECM) method was proposed to simulate the global behavior and local damage of a structure.In order to reflect the local damage and improve the computational efficiency,three-dimensional (3D) solid elements and one-dimensional (1D) beam element were coupled by the multi-point constraint equations.A reduced scale 1?8 model test was simulated by the ECM and a full three dimensional model (3DM) contrastively.The results show that the global behavior and local damages of ECM agree well with the test and 3DM.It is indicated that the proposed method can be used in the structural nonlinear analysis accurately and efficiently.

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.

A review of discrete modeling techniques for fracturing processes in discontinuous rock masses

The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.

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.

Determination of the Normal Contact Stiffness and Integration Time Step for the Finite Element Modeling of Bristle-Surface Interaction

In finite element modeling of impact,it is necessary to define appropriate values of the normal contact stiffness,Kn,and the Integration Time Step(ITS).Because impacts are usually of very short duration,very small ITSs are required.Moreover,the selection of a suitable value of Kn is a critical issue,as the impact behavior depends dramatically on this parameter.In this work,a number of experimental tests and finite element analyses have been performed in order to obtain an appropriate value of Kn for the interaction between a bristle of a gutter brush for road sweeping and a concrete surface.Furthermore,a suitable ITS is determined.The experiments consist of releasing a steel bristle that is placed vertically at a certain distance from a concrete surface and tracking the impact.Similarly,in the finite element analyses,a beam is modeled in free fall and impacting a surface;contact and target elements are attached to the beam and the surface,respectively.The results of the experiments and the modeling are integrated through the principle of conservation of energy,the principle of linear impulse and momentum,and Newton’s second law.The results demonstrate that,for the case studied,Kn and the impact time tend to be independent of the velocity just before impact and that Kn has a very large variation,as concrete is a composite material with a rough surface.Also,the ratio between the largest height of the bristle after impact and the initial height tends to be constant.

Analysis of Do D inkjet printhead performance for printable electronics fabrication using dynamic lumped element modeling and swarm intelligence based optimal prediction

The major challenge in printable electronics fabrication is to effectively and accurately control a drop-on-demand(Do D) inkjet printhead for high printing quality. In this work, an optimal prediction model, constructed with the lumped element modeling(LEM) and the artificial bee colony(ABC) algorithm, was proposed to efficiently predict the combination of waveform parameters for obtaining the desired droplet properties. For acquiring higher simulation accuracy, a modified dynamic lumped element model(DLEM) was proposed with time-varying equivalent circuits, which can characterize the nonlinear behaviors of piezoelectric printhead. The proposed method was then applied to investigate the influences of various waveform parameters on droplet volume and velocity of nano-silver ink, and to predict the printing quality using nano-silver ink. Experimental results show that, compared with two-dimension manual search, the proposed optimal prediction model perform efficiently and accurately in searching the appropriate combination of waveform parameters for printable electronics fabrication.

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.

Vibration Response Analysis of a New Scientific Research Ship Based on Finite Element Modeling

To control the vibration level of ships under construction,MSC Software’s Patran&Nastran modeling solutions can be used to establish a detailed finite element model of a new manned submersible support mother ship based on a line drawing,including the deck layout,bulkhead section,and stiffener distribution.After a comprehensive analysis of the ship simulation conditions,boundaries,and excitation forces of the main operating equipment,modal analysis and calculation of the ship vibration can be conducted.In this study,we calculated and analyzed the vibration response of key points in the stern area of the ship’s main deck and the submersible warehouse area under design loading working conditions.We then analyzed the vibration response of typical decks(including the compass deck,steering deck,captain’s deck,forecastle deck,and main deck)under the main excitation forces and moments(such as the full swing pod and generator sets).The analysis results showed that under DESIDEP working conditions,the vibration of each deck and key areas of the support mother ship could meet the vibration code requirements of the ship’s preliminary design(using the pod excitation and generator sets).Similarly,the vibration response of a scientific research ship under other loading conditions also met the requirements of the code and provided data support for a comprehensive understanding of the ship’s vibration and noise levels.Using actual vibration measurements,the accuracy of the vibration level simulations using finite element modeling was verified,the vibration of each area of the ship comfortably meeting the requirements of the China Classification Society.

3-D finite element modeling for evolution of stress field and interaction among strong earthquakes in Sichuan-Yunnan region

Based on the latest achievement about activities of geological structure, a 3-D finite-element model containing four layers of upper crust, lower crust (two layers) and upper mantle is established in the paper. By repeated tests and revisions, the boundary conditions of the model are determined. And then the background stress field, the stress field caused by fault creep and the stress field triggered by strong earthquake in Sichuan-Yunnan region, as well as their dynamic variations are calculated. The results indicate that the latter earthquake often occurs in the area with positive Coulomb rupture stress change associated with the former one, the former earthquake has a triggering effect on the latter one to a certain extent, and strong earthquake often occur in groups under the background of high stress, which is of great significance for distinguishing seismic anomalies, as well as for improving the level of earthquake prediction.

Multiple linear system techniques for 3D finite element method modeling of direct current resistivity

The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.

Discrete element modeling of acoustic emission in rock fracture

The acoustic emission （AE） features in rock fracture are simulated numerically with discrete element model （DEM）. The specimen is constructed by using spherical particles bonded via the parallel bond model. As a result of the heterogeneity in rock specimen, the failure criterion of bonded particle is coupled by the shear and tensile strengths, which follow a normal probability distribution. The Kaiser effect is simulated in the fracture process, for a cubic rock specimen under uniaxial compression with a constant rate. The AE number is estimated with breakages of bonded particles using a pair of parameters, in the temporal and spatial scale, respectively. It is found that the AE numbers and the elastic energy release curves coincide. The range for the Kaiser effect from the AE number and the elastic energy release are the same. Furthermore, the frequency-magnitude relation of the AE number shows that the value of B determined with DEM is consistent with the experimental data.

Three-Dimensional Simulations of RESET Operation in Phase-Change Random Access Memory with Blade-Type Like Phase Change Layer by Finite Element Modeling

An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell and the blade heater contactor structure by three-dimensional finite element modeling are compared with each other during RESET operation. The simulation results show that the programming region of the phase change layer in the BTL cell is much smaller, and thermal electrical distributions of the BTL cell are more concentrated on the TiN/GST interface. The results indicate that the BTL cell has the superiorities of increasing the heating efficiency, decreasing the power consumption and reducing the RESET current from 0.67mA to 0.32mA. Therefore, the BTL cell will be appropriate for high performance PCRAM device with lower power consumption and lower RESET current.

Isotropic compression behavior of granular assembly with non-spherical particles by X-ray micro-computed tomography and discrete element modeling

The particle morphological properties,such as sphericity,concavity and convexity,of a granular assembly significantly affect its macroscopic and microscopic compressive behaviors under isotropic loading condition.However,limited studies on investigating the microscopic behavior of the granular assembly with real particle shapes under isotropic compression were reported.In this study,X-ray computed tomography(mCT)and discrete element modeling(DEM)were utilized to investigate isotropic compression behavior of the granular assembly with regard to the particle morphological properties,such as particle sphericity,concavity and interparticle frictions.The mCT was first used to extract the particle morphological parameters and then the DEM was utilized to numerically investigate the influences of the particle morphological properties on the isotropic compression behavior.The image reconstruction from mCT images indicated that the presented particle quantification algorithm was robust,and the presented microscopic analysis via the DEM simulation demonstrated that the particle surface concavity significantly affected the isotropic compression behavior.The observations of the particle connectivity and local void ratio distribution also provided insights into the granular assembly under isotropic compression.Results found that the particle concavity and interparticle friction influenced the most of the isotropic compression behavior of the granular assemblies.

Hip Fracture Risk Assessment Based on Different Failure Criteria Using QCT-Based Finite Element Modeling

Precise evaluation of hip fracture risk leads to reduce hip fracture occurrence in individuals and assist to check the effect of a treatment.A subject-specific QCT-based finite element model is introduced to evaluate hip fracture risk using the strain energy,von-Mises stress,and von-Mises strain criteria during the single-leg stance and the sideways fall configurations.Choosing a proper failure criterion in hip fracture risk assessment is very important.The aim of this study is to define hip fracture risk index using the strain energy,von Mises stress,and von Mises strain criteria and compare the calculated fracture risk indices using these criteria at the critical regions of the femur.It is found that based on these criteria,the hip fracture risk at the femoral neck and the intertrochanteric region is higher than other parts of the femur,probably due to the larger amount of cancellous bone in these regions.The study results also show that the strain energy criterion gives more reasonable assessment of hip fracture risk based on the bone failure mechanism and the von-Mises strain criterion is more conservative than two other criteria and leads to higher estimate of hip fracture risk indices.

Modeling of the mechano-electrochemical effect at corrosion defect with varied inclinations on oil/gas pipelines

A 3-dimensional finite element model was built to determine the effect of inclination angle of a corro sion defect on local mechano-electrochemical(M-E)effect in a simulated soil solution.Because of the high effect of the defect inclination angle on the M-E effect,when the inclination angle is 0°(i.e.,the primary axis of the defect parallel to the longitudinal direction of the pipe),the greate st stress concentration level at the defect can be observed,which is associated with the lowest corrosion potential,the greatest anodic current density and the most serious accelerated localized corrosion.When the inclination angle is 90°,the stress concentration level reduces and the corrosion potential becomes less negative,accompanying with the decreased anodic/cathodic current densities.Besides,when the ratio(r_(ca))of the primary axial length of the defect to its secondary axial length is 1,the defect inclination does not affect the stress and the electrochemical corrosion rate at the defect.With the increase of r_(ca),the effect of the defect inclination angle is more apparent.