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.

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.

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.

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.

Finite element modelling of the geodynamic processes of the Central Andes subduction zone:A Reference Model

This paper presents preliminary results of three-dimensional thermomechanical finite-element models of a parameter study to compute the current temperature and stress distribution in the subduction zone of the central Andes （16°S-26°S） up to a depth of 400 km, the bottom of the asthenosphere. For this purpose a simulation running over c. 50,000 years will be realized based on the geometry of a generic subduction zone and an elasto-viscoplastic Drucker-Prager rheology. The kinematic and thermal boundary conditions as well as the rheological parameters represent the current state of the study area. In future works the model will be refined using a systematic study of physical parameters in order to estimate the influence of the main parameters （e.g. viscosity, fault friction, velocity, shear heating） on the results of the reference model presented here. The reference model is kept as simple as possible to be able to estimate the influence of the parameters in future studies in the best possible way, whilst minimizing comnutational time.

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.

Simulation-based Estimation of Thermal Behavior of Direct Feed Drive Mechanism with Updated Finite Element Model

Linear motors generate high heat and cause significant deformation in high speed direct feed drive mechanisms.It is relevant to estimate their deformation behavior to improve their application in precision machine tools.This paper describes a method to estimate its thermal deformation based on updated finite element（FE）model methods.Firstly,a FE model is established for a linear motor drive test rig that includes the correlation between temperature rise and its resulting deformation.The relationship between the input and output variables of the FE model is identified with a modified multivariate input/output least square support vector regression machine.Additionally,the temperature rise and displacements at some critical points on the mechanism are obtained experimentally by a system of thermocouples and an interferometer.The FE model is updated through intelligent comparison between the experimentally measured values and the results from the regression machine.The experiments for testing thermal behavior along with the updated FE model simulations is conducted on the test rig in reciprocating cycle drive conditions.The results show that the intelligently updated FE model can be implemented to analyze the temperature variation distribution of the mechanism and to estimate its thermal behavior.The accuracy of the thermal behavior estimation with the optimally updated method can be more than double that of the initial theoretical FE model.This paper provides a simulation method that is effective to estimate the thermal behavior of the direct feed drive mechanism with high accuracy.

Autologous nerve graft repair of different degrees of sciatic nerve defect:stress and displacement at the anastomosis in a three-dimensional finite element simulation model

In the repair of peripheral nerve injury using autologous or synthetic nerve grafting, the mag- nitude of tensile forces at the anastomosis affects its response to physiological stress and the ultimate success of the treatment. One-dimensional stretching is commonly used to measure changes in tensile stress and strain; however, the accuracy of this simple method is limited. There- fore, in the present study, we established three-dimensional finite element models of sciatic nerve defects repaired by autologous nerve grafts. Using PRO E 5.0 finite element simulation software, we calculated the maximum stress and displacement of an anastomosis under a 5 N load in 10-, 20-, 30-, 40-mm long autologous nerve grafts. We found that maximum displacement increased with graft length, consistent with specimen force. These findings indicate that three-dimensional finite element simulation is a feasible method for analyzing stress and displacement at the anas- tomosis after autologous nerve grafting.

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.

Finite Element Analysis in Combination with Perfectly Matched Layer to the Numerical Modeling of Acoustic Devices in Piezoelectric Materials

The characterization of finite length Surface Acoustic Wave (SAW) and Bulk acoustic Wave (BAW) resonators is addressed here. The Finite Element Analysis (FEA) induces artificial wave reflections at the edges of the mesh. In fact, these ones do not contribute in practice to the corresponding experimental response. The Perfectly Matched Layer (PML) method, allows to suppress the boundary reflections. In this work, we first demonstrate the basis of PML adapted to FEA formalism. Next, the results of such a method are depicted allowing a discussion on the behavior of finite acoustic resonators.

A Finite Element Model of Locked Plating in Femoral Shaft Fractures

Introduction: The Locking Compression Plate (LCP) system is a versatile technology that can be used either through conventional compression plating techniques or as an internal fixator with locking head screws. There have been only a few biomechanical studies examining the role of locked screw configuration on construct stability with most recommendations based upon empirical evidence or data from compression plating. This study will attempt to determine how different locked screw configurations, fracture gaps (distance between bone fragments), and interface gaps (distance between plate and bone) will affect the peak stress(von Mises stress) experienced by the plate-screw construct and, thereby, look at ways to minimize the risk of hardware failure. Materials Methods: A finite element model (FEM) was developed of a transverse mid shaft femoral fracture bridged by an eight-hole titanium LCP. Seven different screw configurations were investigated. Three different fracture gaps and three different interface gaps were studied as well. Results: The 1368 configuration was found to experience the least peak stress of 2.10 GPa while the 2367, 2457, and all filled configurations were found to have the highest peak stress (25.29 GPa, 22.78 GPa, and 23.54 GPa, respectively). Peak stress increased when the interface gap increased. Peak stress also increased as the fracture gap increased, with the largest jump between the 1 mm and 2 mm gaps. Conclusions: Every fracture is unique, and has a vast amount of parameters that must be considered when the surgeon is developing a treatment plan. For transverse femoral shaft fractures, the results of this study suggest that a working length of 2 screw holes on either side of the fracture may also lead to lower peak stress. In addition, our results demonstrate that minimizing the fracture gap and interface gap will lead to decreased stress in the plate-screw construct.

Optimization of Material Coefficients in the Holzapfel-Gasser-Ogden Material Model for the Main Four Ligaments of the Knee Joint-A Finite Element Study

Accurate representation of soft tissue material properties plays a crucial role in computational biomechanics. Several material models have been used for knee ligaments in finite element (FE) studies, including the neo-Hookean model (widely used) and the Holzapfel-Gasser-Ogden (HGO) model (seldom used). While the coefficients of neo-Hookean models for the knee ligaments are available in the literature, limited data exists for the HGO model. Furthermore, no peer-reviewed comparison of these two material models for the knee ligaments while including the 3D representation of the ligaments for both material models is present in the literature. We used mechanical properties from the tensile test experiments in the literature for each ligament to obtain the HGO material coefficients while accounting for the ligaments’ viscoelastic behavior. Resultant coefficients were then used in an Abaqus/explicit knee model to simulate bipedal landing from a jump. The simulations were repeated with neo-Hookean values from the literature. Knee kinematics plus ACL and MCL strains were evaluated and compared for these two material models. The outputs from the simulations with HGO properties were predominantly within 1.5 standard deviations from the mean in-vitro data. When the material properties changed to Neo-Hookean, the outputs for kinematics and strain values were higher than the HGO case, and in most instances, they were outside the experimental range for ACL and MCL strains (by up to 11.35 SD) as well as some ITR angles (by up to 2.86 SD). Reported HGO material model with optimized coefficients produces a more realistic representation of the ligaments’ material properties, and will help improve the outcomes of FE models for more accurate predictions of knee behavior.

Rapid Modeling and Finite Element Analysis of Missile Stiffener

Stiffener is the important missile structure to ensure the structural strength of the missile. In order to improve the design efficiency and the quality of the missile stiffener, the methods of missile stiffener rapid modeling and analysis are proposed. First, the problems of traditional manual modeling of the stiffener are analyzed. According to the problems and actual requirement of modeling, volume decomposition method is used to divide the stiffener into the upper section, the lower section and the web in order for feature analysis and parameter extraction. Then based on the parameters the basic unit decomposed above is created for Boolean operation to establish the stiffener. Finally, a rapid stiffener modeling and analysis program were developed based on UG Open API, the modeling and analysis result validates the feasibility of the method.

Finite element model for arch bridge vibration dynamics considering effect of suspender length adjustment on geometry stiffness matrix

In this paper, we established a finite element (FEM) model to analyze the dynamic characteristics of arch bridges. In this model, the effects of adjustment to the length of a suspender on its geometry stiffness matrix are stressed. The FEM equations of mechanics characteristics, natural frequency and main mode are set up based on the first order matrix perturbation theory. Applicantion of the proposed model to analyze a real arch bridge proved the improvement in the simulation precision of dynamical characteristics of the arch bridge by considering the effects of suspender length variation.

Computer model of the human head-neck and finite element analysis

The difficulty in establishing the finite element model of head and cervical spine is interpreted in the study. A head-neck 3D model is constructed accurately and quickly by the technology of CT scan,the automatically modeling of Mimics software and the RE technology of Geomagic software. Then the finite element model of the head-neck which is close to the real one is set up by the preprocessor of the FEM soft ware ANSYS. After the transient finite element analysis is performed on the model,the historical response of the displacement of the head is obtained. The result is validated by the result of the existed experiment. The stress,as well as the deformation,of nodes in the head and the cervical spine at any time benefits a lot to the clinic study on the injure to the head and neck caused by the impacts. And all the analysis is done by limited computer available.

Torsional Response Analysis of Flexible Pipe Based on Theory and Finite Element Method

As key components connecting offshore floating production platforms and subsea imports, offshore flexible pipes play significant roles in oil, natural gas, and water injection. It is found that torsional failure is one of the failure modes of flexible pipes during transportation and laying. In this paper, a theoretical model(TM) of a flexible pipe section mechanics is established, in which the carcass layer and the pressure armor layer are equivalent to the orthogonal anisotropic layers due to its complex cross-section structure. The calculation results of the TM are consistent with those of a finite element model(FEM), which can accurately describe the torsional response of the flexible pipe.Subsequently, the TM and FEM are used to discuss the influence of boundary conditions on the torsional response.The structure of the flexible pipe is stable when twisted counterclockwise. However, limiting the top axial displacement can improve the axial and radial instability of the tensile armor layer when twisted clockwise. Finally, it is recommended that the flexible pipe can be kept under top fixation during service or installation to avoid torsional failure.