Three-dimensional Finite Element Modeling of a Maxillary Premolar Tooth Based on the Micro-CT Scanning: A Detailed Description

This study describes the details of how to construct a three-dimensional （3D） finite element model of a maxillary first premolar tooth based on micro-CT data acquisition technique, MIMICS soft- ware and ANSYS software. The tooth was scanned by micro-CT, in which 1295 slices were obtained and then 648 slices were selected for modeling. The 3D surface mesh models of enamel and dentin were created by MIMICS （STL file）. The solid mesh model was constructed by ANSYS. After the material properties and boundary conditions were set, a loading analysis was performed to demonstrate the ap- plicableness of the resulting model. The first and third principal stresses were then evaluated. The re- suits showed that the number of nodes and elements of the finite element model were 56 618 and 311801, respectively. The geometric form of the model was highly consistent with that of the true tooth, and the deviation between them was ~）.28%. The loading analysis revealed the typical stress patterns in the contour map. The maximum compressive stress existed in the contact points and the maximum tensile stress existed in the deep fissure between the two cusps. It is concluded that by using the micro-CT and highly integrated software, construction of the 3D finite element model with high quality will not be difficult for clinical researchers.

Validation and application of three-dimensional discontinuous deformation analysis with tetrahedron finite element meshed block

In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.

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.

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.

Alternating Direction Finite Volume Element Methods for Three-Dimensional Parabolic Equations

This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

CHEBYSHEV PSEUDOSPECTRAL-HYBRID FINITE ELEMENT METHOD FOR THREE-DIMENSIONAL VORTICITY EQUATION

In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability and the convergence are proved strictly.The numerical results show the advantages of this method.The technique in this paper is also applicable to other three-dimensional nonlinear problems in fluid dynamics.

THE CHARACTERISTIC FINITE ELEMENT ALTERNATING-DIRECTION METHOD AND ANALYSIS FOR THREE-DIMENSIONAL NUMERICAL RESERVOIR SIMULATION

Petroleum science has made remarkable progress in organic geochcmistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the knowledge of its evolutionary history and especially the numerical computation of fluid flow and the history of its changes under heat is vital. The mathematical model can be described as a coupled system of nonlinear partial differentical equations with initial-boundary value problems. This thesis, from actual conditions such as the effect of fluid compressibility and the three-dimensional characteristic of large-scale science-engineering computation, we put forward a kind of characteristic finite element alternating-direction schemes and obtain optimal order estimates in L^2 norm for the error in the approximate assumption.

Three-dimensional anisotropic induced polarization modeling using finite element method

With the development of geophysical exploration technology,the anisotropy of underground media has got more and more attention.At present,there are few studies on the anisotropy of the induced polarization method.This article explores the effect of anisotropy on the underground media of the induced polarization method under three-dimensional complex terrain.The research work transforms the underground electric field control equation into a variational problem,and use the unstructured finite element method to construct a large linear equation system for solving electric potentials.By the sparse matrix compression technique and symmetric successive over-relaxation preconditioned conjugate gradient algorithm(SSOR-PCG)to solve the equation system.Finally,the article uses the classic central gradient array method to obtain the forward apparent polarizability value.The calculation results of the model find that different anisotropic conditions will significantly affect the forward results which show a strong directional correlation,revealing the importance of considering anisotropy in practical work.

Evaluation value of three-dimensional finite element model analysis for bone mineral density and bone metabolism activity in patients with osteoporosis

Objective: To study the evaluation value of three-dimensional finite element model analysis for bone mineral density (BMD) and bone metabolism activity in patients with osteoporosis. Methods: A total of 218 patients who were diagnosed with osteoporosis in the hospital between February 2014 and January 2017 were collected as observation group, and 100 healthy volunteers who received physical examination in the hospital during the same period were selected as normal control group. The femoral head of the two groups was analyzed by three-dimensional finite element model, and the femoral head BMD levels and serum bone metabolism index contents were measured. Pearson test was used to evaluate the evaluation value of femoral head three-dimensional finite element model for osteoporosis. Results: The cancellous bone and cortical bone Von Mises stress value of observation group were lower than those of normal control group, and femoral neck BMD value of observation group was lower than that of normal control group;serum bone metabolism index BGP content was lower than that of normal control group while NBAP, TRACP-5b and CTX-1 contents were higher than those of normal control group. Pearson test showed that the cancellous bone and cortical bone Von Mises stress value of patients with osteoporosis were directly correlated with BMD value and bone metabolism index contents. Conclusion: The three-dimensional finite element model analysis resultsof patients with osteoporosis can objectively reflect the femoral headBMD value and bone metabolism activity, and is a reliable way to evaluate the risk of long-term fractures.

Three-dimensional finite element analysls for external midface distraction after different osteotomy in patients with cleft lip and palate

Objective To study three - dimensional finite element analysis for external midface distraction after different osteotomy in patients with cleft lip and palate ( Clp) . Methods Three - dimensional Fem models of Le Fort Ⅰ,Ⅱand Ⅲ,osteotomy in Clp patients were estabolished. External midface distraction were simulated. An anteriorly and inferiorly directed 900 g force was

Combined load bearing capacity of rigid piles embedded in a crossanisotropic clay deposit using 3D finite element lower bound

In this study,an iterative-based three-dimensional finite element lower bound in association with the second-order cone programming method is adopted to evaluate the limit load of a single pile embedded in cross-anisotropic soils under general loading condition.The lower bound solutions of the pile embedded in an anisotropic soil deposit can be found by formulating the element equilibrium,equilibrium of shear and normal stresses along discontinuities,boundary conditions,yield function,and optimizing the objective function through the second-order cone programming method in conjunction with an iterative-based update procedure.A general loading condition is considered to profile the expansion of the safe load in the vertical-horizontal-moment(V-H-M)space.The results of this study are compared and validated against three different cases including an isotropic lateral loading,anisotropic end bearing capacity,and a pile embedded in an isotropic soil deposit under general loading condition.A parametric study is conducted to evaluate the impact of different influencing factors.It was found that the effect of anisotropy on the variation of lateral limit load of a single pile is more pronounced than the corresponding vertical and bending moment limit loads,whereas the interface properties have more significant effects on the vertical and bending moment limit loads in comparison to the lateral limit load.

Parametric Study of Hip Fracture Risk Using QCT-Based Finite Element Analysis

Various parameters such as age,height,weight,and body mass index(BMI)influence the hip fracture risk in the elderly which is the most common injury during the sideways fall.This paper presents a parametric study of hip fracture risk based on the gender,age,height,weight,and BMI of subjects using the subject-specific QCT-based finite element modelling and simulation of single-leg stance and sideways fall loadings.Hip fracture risk is estimated using the strain energy failure criterion as a combination of bone stresses and strains leading to more accurate and reasonable results based on the bone failure mechanism.Understanding the effects of various parameters on hip fracture risk can help to prescribe more accurate preventive and treatment plans for a community based on the gender,age,height,weight,and BMI of the population.Results of this study show an increase in hip fracture risk with the increase of age,body height,weight,and BMI in both women and men under the single-leg stance and the sideways fall configurations.

Experimental and Finite Element Analysis for Multi-lead Rubber Bearings:A Comparative Study

The mechanical properties of multi-lead rubber bearings （MLRBs） were investigated by experiment and finite element analysis. First, the vertical stiffness, horizontal stiffness and yielded shear force were tested for four MLRB specimens and two specimens of the single-lead rubber bearings （ SLRBs）. Then, the MLRBs were modeled by the explicit finite element analysis software ANSYS/ LS-DYNA, in order to evaluate the horizontal force-displacement hysteretic curves under static vertical and dynamical horizontal loadings. The disagreement between the tested and theoretical values was less than 11.4%, and MLRBs and SLRBs were similar in vertical stiffness, pre-yield stiffness and yield stiffness.

Development of a comprehensive finite element cervical spine model for studying neck injury of pilot

Introduction-The cervical spine is subjected to injury frequently,especially among pilots who are usually on the condition of high acceleration.Injuries of the cervical spine will be potential risk of damage to the spinal cord,which could be result in life threatening

Importance of Three-Dimensional Piezoelectric Coupling Modeling in Quantitative Analysis of Piezoelectric Actuators

This paper demonstrates the importance of three-dimensional(3-D)piezoelectric coupling in the electromechanical behavior of piezoelectric devices using three-dimensional finite element analyses based on weak and strong coupling models for a thin cantilevered piezoelectric bimorph actuator.It is found that there is a significant difference between the strong and weak coupling solutions given by coupling direct and inverse piezoelectric effects(i.e.,piezoelectric coupling effect).In addition,there is significant longitudinal bending caused by the constraint of the inverse piezoelectric effect in the width direction at the fixed end(i.e.,3-D effect).Hence,modeling of these effects or 3-D piezoelectric coupling modeling is an electromechanical basis for the piezoelectric devices,which contributes to the accurate prediction of their behavior.

Virtual Element Formulation for Finite Strain Elastodynamics Dedicated to Professor Karl Stark Pister for his 95th birthday

The virtual element method(VEM)can be seen as an extension of the classical finite element method(FEM)based on Galerkin projection.It allows meshes with highly irregular shaped elements,including concave shapes.So far the virtual element method has been applied to various engineering problems such as elasto-plasticity,multiphysics,damage and fracture mechanics.This work focuses on the extension of the virtual element method to efficient modeling of nonlinear elasto-dynamics undergoing large deformations.Within this framework,we employ low-order ansatz functions in two and three dimensions for elements that can have arbitrary polygonal shape.The formulations considered in this contribution are based on minimization of potential function for both the static and the dynamic behavior.Generally the construction of a virtual element is based on a projection part and a stabilization part.While the stiffness matrix needs a suitable stabilization,the mass matrix can be calculated using only the projection part.For the implicit time integration scheme,Newmark-Method is used.To show the performance of the method,various two-and three-dimensional numerical examples in are presented.

Stress Distribution in Maxillary Central Incisors without Ferrules: A Finite Element Analysis of Post and Core Systems

Purpose: The purpose of this study was to identify optimal post and core materials for central incisors without ferrules using three-dimensional finite element analysis and three-point bending tests. Methods: Stress analyses were performed with six models: cast metal post and core (MP), composite resin core alone, straight fiber-reinforced post-composite resin core (FSR), tapered fiber-reinforced post-composite resin core, straight titanium post-composite resin core (TSR), and tapered titanium post-composite resin core (TTR). A 100-N load was applied to the lingual surface at a 45° angle to the long axis of the tooth. Maximum von Mises stress distributions were calculated with finite element analysis software. Five samples each of composite resin, straight fiber-reinforced post, straight titanium post, straight fiber-reinforced post and composite resin, and straight titanium post and composite resin were subjected to three-point bending tests, followed by analysis of variance and Tukey’s multiple comparison test. Results: Stress distribution was optimal on TTR. Maximum von Mises stress on the cervical side of the post was greatest in TSR (693 MPa) and TTR (556 MPa). Maximum stress on the apical side of the post was greatest in MP (110 MPa). Maximum stress in surrounding dentin was lowest in MP (203 MPa) and TTR (250 MPa). Gap distance was smallest in MP (0.09 mm) and largest in FSR (0.26 mm). Mean maximum three-point bending force was lowest in composite resin (26.9 N/mm) and highest in titanium post and composite resin (97.1 N/mm). Titanium post bending strength was consistently greater than that of the fiber-reinforced post (p < 0.01). Conclusion: These results revealed optimal stress distribution and high bending strength with the tapered titanium post and resin combination, suggesting that this combination can most effectively prevent root or post fracture in an anterior tooth without a ferrule.

THREE-DIMENSIONAL ELLIPTIC CRACK UNDER IMPACT LOADING

The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.