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 analysis on effects of tunnel construction on nearby pile foundation

A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.

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.

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.

Finite element analysis of anchor plates using non-coaxial models

The non-coaxial model simulating the non-coincidence between the principal stresses and the principal plastic strain rates is employed within the framework of finite element method（FEM） to predict the behaviors of anchors embedded in granular material.The non-coaxial model is developed based on the non-coaxial yield vertex theory,and the elastic and conventional coaxial plastic deformations are simulated by using elasto-perfectly plastic Drucker-Prager yield function according to the original yield vertex theory.Both the horizontal and vertical anchors with various embedment depths are considered.Different anchor shapes and soil friction and dilation angles are also taken into account.The predictions indicate that the use of non-coaxial models leads to softer responses,compared with those using conventional coaxial models.Besides,the predicted ultimate pulling capacities are the same for both coaxial and non-coaxial models.The non-coaxial influences increase with the increasing embedment depths,and circular anchors lead to larger non-coaxial influences than strip anchors.In view of the fact that the design of anchors is mainly determined by their displacements,ignoring the non-coaxiality in finite element numerical analysis can lead to unsafe results.

Experimental and finite element analysis of tibial stress fractures using a rabbit model

AIM: To determine if rabbit models can be used to quantify the mechanical behaviour involved in tibial stress fracture(TSF) development.METHODS: Fresh rabbit tibiae were loaded under compression using a specifically-designed test apparatus. Weights were incrementally added up to a load of 30 kg and the mechanical behaviour of the tibia was analysed using tests for buckling, bone strain and hysteresis. Structural mechanics equations were subsequently employed to verify that the results were within the range of values predicted by theory. A finite element(FE) model was developed using cross-sectional computer tomography(CT) images scanned from one of the rabbit bones, and a static load of 6 kg(1.5 times the rabbit's body weight) was applied to represent running. The model was validated using the experimental strain gauge data, then geometric and elemental convergence tests were performed in order to find the minimum number of cross-sectional scans and elements respectively required for convergence. The analysis was then performed using both the model and the experimental results to investigate the mechanical behaviour of the rabbit tibia under compressive load and to examine crack initiation.RESULTS: The experimental tests showed that un der a compressive load of up to 12 kg, the rabbit tibia demonstrates linear behaviour with little hysteresis Up to 30 kg, the bone does not fail by elastic buckling however, there are low levels of tensile stress which predominately occur at and adjacent to the anterio border of the tibial midshaft: this suggests that fatigue failure occurs in these regions, since bone under cycli loading initially fails in tension. The FE model predic tions were consistent with both mechanics theory and the strain gauge results. The model was highly sensi tive to small changes in the position of the applied load due to the high slenderness ratio of the rabbit s tibia. The modelling technique used in the curren study could have applications in the development o human FE models of bone, where, unlike rabbit tibia the model would be relatively insensitive to very sma changes in load position. However, the rabbit mode itself is less beneficial as a tool to understand the me chanical behaviour of TSFs in humans due to the sma size of the rabbit bone and the limitations of human scale CT scanning equipment.CONCLUSION: The current modelling technique could be used to develop human FE models. However, the rabbit model itself has significant limitations in under standing human TSF mechanics.

Three-dimensional forward modeling for magnetotelluric sounding by finite element method

A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.

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.

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.

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.

The finite element analysis of articular cartilage fiber-reinforced composite model under rolling load

Articular cartilage is a layer of low-friction,load-bearing soft hydrated tissue covering bone-ends in diarthrosis,which plays an important role in spreading the load,reducing the joint contact stress,joint friction and wear during exercise.The vital mechanical function

Seismic performance of double-skin steel-concrete composite box piers: Part Ⅱ—Nonlinear finite element analysis

An accurate finite element （ FE） model was constructed to examine the hysteretic behavior of double-skin steel-concrete composite box （ DSCB） piers for further understanding the seismic performance of DSCB piers; where the local buckling behavior of steel tubes, the confinement of the in-filled concrete and the interface action between steel tube and in-filled concrete were considered. The accuracy of the proposed FE model was verified by the bidirectional cyclic loading test results. Based on the validated FE model, the effects of some key parameters, such as section width to steel thickness ratio, slenderness ratio, aspect ratio and axial load ratio on the hysteretic behavior of DSCB piers were investigated. Finally, the skeleton curve model of DSCB piers was proposed. The numerical simulation results reveal that the peak strength and elastic stiffness decrease with the increase of the section width to steel thickness ratio. Moreover, the increase of the slenderness ratio may result in a significant reduction in the peak strength and elastic stiffness while the ultimate displacement increases. The proposed skeleton curve model can be taken as a reference for seismic performance analyses of the DSCB piers.

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 analysis of ionic liquid gel soft actuator

A new type of soft actuator material-ionic liquid gel （ILG）, which consists of HEMA, BMIMBF4, and TiO2, can be transformed into gel state under the irradiation of ultraviolet （UV） light. In this paper, Mooney-Rivlin hyperelastic model of finite element method is proposed for the first time to study the properties of the ILG. It has been proved that the content of TiO2 has a great influence on the properties of the gel, and Young＇s modulus of the gel increases with the increase of its content, despite of reduced tensile deformation. The results in this work show that when the TiO2 content is 1.0 wt%, a large tensile deformation and a strong Young＇s modulus can be obtained to be 325% and 7.8 kPa, respectively. The material parameters of ILG with TiO2 content values of 0.2 wt%, 0.5 wt%, 1.0 wt%, and 1.5 wt% are obtained, respectively, through uniaxial tensile tests, including C10, C01, C20, C11, C02, C30, C21, C12, and C03 elements. In this paper, the large-scaled general finite element software ANSYS is used to simulate and analyze the ILG, which is based on SOLID186 element and nonlinear hyperelastic Mooney-Rivlin model. The finite element simulation analysis based stress-strain curves are almost consistent with the experimental stress-strain curves, and hence the finite element analysis of ILG is feasible and credible. This work presents a new direction for studying the performance of soft actuator for the ILG, and also contributes to the design of soft robot actuator.

Three-Dimensional Finite Element Numerical Simulation and Physical Experiment for Magnetism-Stress Detecting in Oil Casing

The casing damage has been a big problem in oilfield production. The current detection methods mostly are used after casing damage, which is not very effective. With the rapid development of China's offshore oil industry, the number of offshore oil wells is becoming larger and larger. Because the cost of offshore oil well is very high, the casing damage will cause huge economic losses. What's more, it can also bring serious pollution to marine environment. So the effective methods of detecting casing damage are required badly. The accumulation of stress is the main reason for the casing damage. Magnetic anisotropy technique based on counter magnetostriction effect can detect the stress of casing in real time and help us to find out the hidden dangers in time. It is essential for us to prevent the casing damage from occurring. However, such technique is still in the development stage. Previous studies mostly got the relationship between stress and magnetic signals by physical experiment, and the study of physical mechanism in relative magnetic permeability connecting the stress and magnetic signals is rarely reported. The present paper uses the ANSYS to do the three-dimensional finite element numerical simulation to study how the relative magnetic permeability works for the oil casing model. We find that the quantitative relationship between the stress' s variation and magnetic induction intensity's variation is: Δδ =K* ΔB, K = 8.04×109, which is proved correct by physical experiment.

Experimental and finite element analysis for fracture of coating layer of galvannealed steel sheet

Mechanical properties of galvannealed (GA) steel sheet used for automotive exposed panel and predicted failure phenomenon of its coating layer were evaluated using finite element method. V-bending test was performed to understand better the fracture of coating layer of GA steel sheet during plastic deformation. Yield strength of the coating layer was calculated by using a relative difference between hardness of coating layer measured from the nano-indentation test and that of substrate. To measure shearing strength at the interface between substrate and coating layer, shearing test with two specimens attached by an adhesive was carried out. Using the mechanical properties measured, a series of finite element analyses coupled with a failure model was performed. Results reveal that the fracture of coating layer occurs in an irregular manner at the region where compressive deformation is dominant. Meanwhile, a series of vertical cracks perpendicular to material surface are observed at the tensile stressed-region. It is found that 0.26-0.28 of local equivalent plastic strain exists at the coating and substrate at the beginning of failure. The fracture of coating layer depends on ductility of the coating layer considerably as well.

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.

Rock-soil slope stability analysis by two-phase random media and finite elements

To investigate the strong random nature of the geometric interfaces between soil and rock, a rock-soil slope is considered as a two-phase random medium. A nonlinear translation of a Gaussian field is utilized to simulate the two-phase random media, such that the soil(or rock) volume fraction and the inclination of the soil layer can be examined. The finite element method with random media incorporated as the material properties is used to determine the factor of safety of the rock-soil slope. Monte-Carlo simulations are used to estimate the statistical characteristics of the factor of safety. The failure mode of the rock-soil slope is examined by observing the maximum principal plastic strain at incipient slope failure. It is found that the critical surface of a rock-soil slope is fairly irregular, and it significantly differs from that of a pure soil slope. The factor of safety is sensitive to the soil volume faction, but it is predictable. The average factor of safety could be well predicted by the weighted harmonic average between the strength of soil and rock; the prediction model is practical and simple. Parametric studies on the inclination of the soil layer demonstrate that the most instable scenario occurs when the slope angle is consistent with the inclination of the soil layer.

Dynamic Characteristics of Inflated Torus Using Finite Element Analysis

In recent years, inflatable structures have been a subject of interest for space applications such as communication antenna, solar thermal propulsion and entry/landing systems. The inflatable structures characterized by high strength-to-mass ratios, minimal stowage volume, which makes them suitable for cost-effective large space structures. A typical example for the inflatable structure is the inflated torus which often used in order to provide structure support. In this study, our main focus is to understand the dynamic characteristics of an inflated torus in order to formulate an accurate mathematical model suitable for active vibration control applications. A commercial finite element package, ANSYS, is used to model the inflated torus. To verify the model the obtained frequencies and mode shapes are compared with the published results, which are derived using analytical approach, the verification shows a good agreement between the FEM and the analytical results. Based on the verified model, parametric study was investigated; the material thickness increase causes the natural frequencies decrease, while the increase of the inflation pressure simply results in stiffening the ring, which means that the natural frequency increased. The FEM analysis gives an easy and fast way for the vibration analysis of the structures compared with the complicated analytical solutions.