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.

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.

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.

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.

Finite Element Implementation of the Exponential Drucker-Prager Plasticity Model for Adhesive Joints

This paper deals with the numerical implementation of the exponential Drucker-Parger plasticitymodel in the commercial finite element software,ABAQUS,via user subroutine UMAT for adhesive joint simulations.The influence of hydrostatic pressure on adhesive strength was investigated by a modified Arcan fixture designed particularly to induce a different state of hydrostatic pressure within an adhesive layer.The developed user subroutine UMAT,which utilizes an associated plastic flow during a plastic deformation,can provide a good agreement between the simulations and the experimental data.Better numerical stability at highly positive hydrostatic pressure loads for a very high order of exponential function can also be achieved compared to when a non-associated flow is used.

A rigid-flexible coupling finite element model of coupler for analyzing train instability behavior during collision

Rail vehicles generate huge longitudinal impact loads in collisions.If unreasonable matching exists between the compressive strength of the intermediate coupler and the structural strength of the car body,the risk of car body structure damage and train derailment will increase.Herein,a four-stage rigid-flexible coupling finite element model of the coupler is established considering the coupler buckling load.The influence of the coupler buckling load on the train longitudinal-vertical-hori-zontal buckling behavior was studied,and the mechanism of the train horizontal buckling instability in train collisions was revealed.Analysis results show that an intermediate coupler should be designed to ensure that the actual buckling load is less than the compressive load when the car body structure begins to deform plastically.The actual buckling load of the coupler and the asymmetry of the structural strength of the car body in the lateral direction are two important influencing factors for the lateral buckling of a train collision.If the strength of the two sides of the car body structure in the lateral direction is asymmetrical,the deformation on the weaker side will be larger,and the end of the car body will begin to deflect under the action of the coupler force,which in turn causes the train to undergo sawtooth buckling.

Behaviour of non-ballast pre-stressed and precast track structures in high speed railway based on multiscale finite element model

In order to make further study on the mechanical property of CRTSIII type slab non-ballast track structures,which was self-designed in China,based on the method of the multiscale finite element model(FEM),the traditional FEM of slab non-ballast track structures was improved.The multiscale FEM of CRTSII type slab nonballast track structures was established based on the general finite element program ABAQUs.Then the comparative calculation was made between various FEMs,showing that the high solution precision,fast modelling speed and high solution efficiency could be obtained.Therefore,the multiscale FEM was suitable for the parametric study on mechanical behaviour of CRTSII type slab non-ballast track structures,and then the key influence factor and constructions could be optimized.

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 of the viscoelastic responses of the eye during microvolumetric changes

A linear viscoelastic finite element model was built to investigate factors that influenced the intraocular pressure (IOP) elevations due to micro-volumetric changes in the eye at three different rates. The viscoelastic properties of the cornea and the sclera, including the instantaneous modulus, equilibrium modulus, and relaxation time constants, parametrically varied to examine their effects on IOP elevations at different rates of volumetric changes. The simulated responses were in good agreement with the previously reported experimental results obtained from porcine globes, showing the general trend of higher IOP elevations at faster rates. The simulations showed that all viscoelastic properties influenced the profile of the dynamic IOP due to volumetric changes, and the relative significance of a specific parameter was highly dependent on the rate of change.

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.

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.

Conversion between solid and beam element solutions of finite element method based on meta-modeling theory:development and application to a ramp tunnel structure

In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.

A 3D finite element modeling on the relationship between crustal stress and earthquake activity in Liaoning and its vicinity

Based on the latest result in research on 3D seismic wave velocity structure of crust and uppermost mantle and taking geological setting and fracture zones into consideration, a 3D geological model for the studied region is built up. The boundary constraint and force loading boundary condition for the model are determined according to the characteristics of crustal stress field deduced from earthquake focal mechanism and in-situ stress measurement data. Using linear elastic material model a 3D finite element modeling is conducted to study the characteristics of crustal stress field. A comparison analysis between the simulated stress field and earthquake locations reveals that the moderate and strong earthquakes generally occurred in the zones with high shear stress gradient. Furthermore, the paper notices a few potential earthquake-prone regions.

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.

A non-linear finite element model of human L4-L5 lumbar spinal segment with three-dimensional solid element ligaments

This paper establishes a non-linear finite element model (NFEM) of L4-L5 lumbar spinal segment with accurate three-dimensional solid ligaments and intervertebral disc. For the purpose, the intervertebral disc and surrounding ligaments are modeled with four-nodal three-dimensional tetrahedral elements with hyper-elastic material properties. Pure moment of 10 N·m without preload is applied to the upper vertebral body under the loading conditions of lateral bending, backward extension, torsion, and forward flexion, respectively. The simulate relationship curves between generalized forces and generalized displacement of the NFEM are compared with the in vitro experimental result curves to verify NFEM. The verified results show that: (1) The range of simulated motion is a good agreement with the in vitro experimental data; (2) The NFEM can more effectively reffect the actual mechanical properties than the FE model using cable and spring elements ligaments; (3) The NFEM can be used as the basis for further research on lumbar degenerative diseases.

The Effect of Ligament Modeling Technique on Knee Joint Kinematics: A Finite Element Study

Finite element (FE) analysis has become an increasingly popular technique in the study of human joint biomechanics, as it allows for detailed analysis of the joint/tissue behavior under complex, clinically relevant loading conditions. A wide variety of modeling techniques have been utilized to model knee joint ligaments. However, the effect of a selected constitutive model to simulate the ligaments on knee kinematics remains unclear. The purpose of the current study was to determine the effect of two most common techniques utilized to model knee ligaments on joint kinematics under functional loading conditions. We hypothesized that anatomic representations of the knee ligaments with anisotropic hyperelastic properties will result in more realistic kinematics. A previously developed, extensively validated anatomic FE model of the knee developed from a healthy, young female athlete was used. FE models with 3D anatomic and simplified uniaxial representations of main knee ligaments were used to simulate four functional loading conditions. Model predictions of tibiofemoral joint kinematics were compared to experimental measures. Results demonstrated the ability of the anatomic representation of the knee ligaments (3D geometry along with anisotropic hyperelastic material) in more physiologic prediction of the human knee motion with strong correlation (r ≥ 0.9 for all comparisons) and minimum deviation (0.9° ≤ RMSE ≤ 2.29°) from experimental findings. In contrast, non-physiologic uniaxial elastic representation of the ligaments resulted in lower correlations (r ≤ 0.6 for all comparisons) and substantially higher deviation (2.6°≤ RMSE ≤ 4.2°) from experimental results. Findings of the current study support our hypothesis and highlight the critical role of soft tissue modeling technique on the resultant FE predicted joint kinematics.

Finite Element Modeling Of Low Heat Conducting Building Bricks

Heat conduction through conventional and interlocking building bricks with cavities was studied in this work. Heat transfer analysis was carried out using MATLAB? partial differential equation toolbox. Regular and staggered hole arrangements were studied. Results showed that four staggered holed interlocking bricks were effective in thermal resistance into the bricks and increasing the holes beyond four did not give any thermal resistance advantage. For the conventional bricks staggered holes did not give any thermal resistance advantage but the four-holed bricks were also adjudged to be effective in thermal resistance into the brick surface. Increasing the number of holes beyond four in conventional bricks did give some thermal resistivity advantage but very minimal. Structural strengths of holed bricks were not considered in this study.