Clinical significance of three dimensional finite element analysis on humerus fracture

Objective: To treat humerus fracture with three dimensional pattern and finite element analysis, providing mechanical basis for treating humerus fracture. Methods: Humerus pattern was established based on the CT images, and calculation was done by ANSYS5. 6 software. Three dimensional ten-node tetrahedron unit was selected and were divided into 2 729 nodes, 49 041 units. Distribution and amount of axial compression of humerus were analyzed when clip angle was 30°, 45°, 90° between fracture face and axial line with fixed X, Y, Z directions. Results: The distribution of stress was greatly different between fracture face and non fracture face. Stress in fracture part was fairly concentrated with incomplete symmetric distribution around the center of fracture face; Greater stress distributed in the regions 10 mm from fracture face, which was 2-3 times that of other stress regions. Conclusion: Required load must be estimated under various conditions as to select the suitable internal fixation implants during the treatment of humerus fracture, which can provide helpful stress environment for fracture healing.

Three dimensional finite element analysis of anatomic distal radius Nitinol memory connector treating distal radius fracture

Objective: To study the memory biomechanical character of anatomic distal radius Nitinol memory connector (DRMC) in treating distal radius fracture. Methods: Establishing three dimensional model and finite element analysis, we calculated the stress in and around the fracture faces when distal radius fracture was fixated with DRMC. Results: Axial holding stress produced by holding part of DRMC on distal radius was 14.66 MPa. The maximum stress of holding part was 40-70 MPa, the minimum stress was 3-7 MPa,and the stress of compression part was 20-40 MPa. Conclusion: The distribution of stress produced by DRMC around the fracture line is reasonable, and axial holding stress can help stabilize fracture during earlier period. The existence of longitudal compression and memory effect can transfer fixated disused section into developed section and enhance fracture healing.

Three dimensional finite element analysis of acetabulum loaded by static stress and its biomechanical significance

Objective:To explore the mechanical behavior of acetabulum loaded by static stress and provide the mechanical basis for clinical analysis and judgement on acetabular mechanical distribution and effect of static stress.Methods:By means of computer simulation, acetabular three dimensional model was input into three dimensional finite element analysis software ANSYS7.0. The acetabular mechanical behavior was calculated and the main stress value, stress distribution and acetabular unit displacement in the direction of main stress were analyzed when anterior wall of acetabulum and acetabular crest were loaded by 1 000 N static stress. Results:When acetabular anterior wall loaded by X direction and Z direction composition force, the stress passed along 4 directions: (1)from acetabular anterior wall to pubic symphysis along superior branch of pubis firstly, (2)from acetabular anterior wall to cacroiliac joint along pelvic ring,(3)in the acetabulum,(4)from the suffered point to ischium. When acetabular crest loaded by X direction and Y direction composition force, the stress transmitted to 4 directions: (1)from acetabular crest to ilium firstly, (2)from suffered point to cacroiliac joint along pelvic ring,(3) in the acetabulum ,(4)along the pubic branch,but no stress transmitted to the ischium branch.Conclusion:Analyzing the stress distribution of acetabulum and units displacement when static stress loaded can provide internal fixation point for acetabular fracture treatment and help understand the stress distribution of acetabulum.

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 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.

3D finite element analysis on pile-soil interaction of passive pile group

The interaction between pile and soft soil of the passive pile group subjected to soil movement was analyzed with three-dimensional finite element model by using ANSYS software. The soil was assumed to be elastic-plastic complying with the Drucker-Prager yield criterion in the analysis. The large displacement of soil was considered and contact elements were used to evaluate the interaction between pile and soil. The influences of soil depth of layer and number of piles on the lateral pressure of the pile were investigated, and the lateral pressure distributions on the (2×1) pile group and on the (2×2) pile group were compared. The results show that the adjacent surcharge may result in significant lateral movement of the soft soil and considerable pressure on the pile. The pressure acting on the row near the surcharge is higher than that on the other row, due to the "barrier" and arching effects in pile groups. The passive load and its distribution should be taken into account in the design of the passive piles.

A multiscale 3D finite element analysis of fluid/solute transport in mechanically loaded bone

The transport of fluid, nutrients, and signaling molecules in the bone lacunar-canalicular system （LCS） is critical for osteocyte survival and function. We have applied the fluorescence recovery after photobleaching （FRAP） approach to quantify load-induced fluid and solute transport in the LCS in situ, but the measurements were limited to cortical regions 30-50 μm underneath the periosteum due to the constrains of laser penetration. With this work, we aimed to expand our understanding of load-induced fluid and solute transport in both trabecular and cortical bone using a multiscaled image-based finite element analysis （FEA） approach. An intact murine tibia was first re-constructed from microCT images into a three-dimensional （3D） linear elastic FEA model, and the matrix deformations at various locations were calculated under axial loading. A segment of the above 3D model was then imported to the biphasic poroelasticity analysis platform （FEBio） to predict load-induced fluid pressure fields, and interstitial solute/fluid flows through LCS in both cortical and trabecular regions. Further, secondary flow effects such as the shear stress and/or drag force acting on osteocytes, the presumed mechano-sensors in bone, were derived using the previously developed ultrastructural model of Brinkman flow in the canaliculi. The material properties assumed in the FEA models were validated against previously obtained strain and FRAP transport data measured on the cortical cortex. Our results demonstrated the feasibility of this computational approach in estimating the fluid flux in the LCS and the cellular stimulation forces （shear and drag forces） for osteocytes in any cortical and trabecular bone locations, allowing further studies of how the activation of osteocytes correlates with in vivo functional bone formation. The study provides a promising platform to reveal potential cellular mechanisms underlying the anabolic power of exercises and physical activities in treating patients with skeletal deficiencies.

Characteristics of an Axial-flux Permanent Magnet Synchronous Machine with Contra-rotating Rotors under Unbalanced Load Condition from 3-D Finite Element Analysis

During recent years,the axial-flus PMSM with contra-rotating rotors has become a hot topic in academic research due to its high efficiency and simple structure.However,its back-EMF may be distorted under the condition of different angular positions.This paper investigates characteristics of the novel motor used for contra-propeller driving.Considering the torque ripple and current oscillation under unbalanced load condition,this paper analyzes the distorted back-EMF of the machine when its two rotors get different angular positions during rotating.The analysis results are validated by transient-magnetic 3-D FEA method,which the 3-D FEA software is used to model this motor and transient simulations are carried out to obtain its magnetic characteristic and main performances.A main focus is put on the back-EMF characteristic with different angular positions between the two rotors.Furthermore,the characteristic of torque production under unbalanced load is investigated.Finally,a prototype motor is fabricated to validate the analyses of this paper.

Optimal Design Strategies of Femur Tumor Hyperthermia Based on Finite Element Analysis of Temperature Field

A 3D femoral model was built to obtain the three-dimensional temperature distribution of femur and its surrounding tissues and provide references for clinical applications. According to the relationship between gray-value and material properties,the model was assigned with various materials to make sure that it is more similar to the real femur in geometry and physical properties. 3D temperature distribution is obtained by using finite element analysis software ANSYS 11. 0 on the basis of heat conduction theory,Laplace equation,Pennes bio-heat transfer equation,thermo physical parameters of bone tissues,the boundary condition,and initial conditions. Taken the asymmetry of the 3D distribution of temperature into account,it is necessary to adopt the heating method with multiple heat sources. This method can ensure that the temperature fields match well with the tumor tissues and kill the tumor cells efficiently under the condition of protecting the normal tissues from damage. The analysis results supply important guidance for determining the needle position and the needle number and controlling the intensity of heating.

3D Finite Element Analysis of a Man Hip Joint Femur under Impact Loads

Objective： The biomechanical characters of the bone fracture of the man femoral hip joint under impact loads are explored. Methods ：A biosystem model of the man femoral hip joint by using the GE （ General Electric） lightspeed multi-lay spiral CT is conducted. A 3D finite element model is established by employing the finite element software ANSYS. The FE analysis mainly concentrates on the effects of the impact directions arising from intense movements and the parenchyma on the femoral hip joint on the stress distributions of the proximal femur. Results：The parenchyma on the hip joint has relatively large relaxation effect on the impact loads. Conclusion：Effects of the angle δ of the impact load to the anterior direction and the angle γ of the impact load to the femur shaft on the bone fracture are given;δ has larger effect on the stress and strain distributions than the angle γ,which mainly represents the fracture of the upper femur including the femoral neck fracture when the posterolateral femur is impacted, consistent with the clinical resuits.

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.

3D Finite Element Simulation of Tunnel Boring Machine Construction Processes in Deep Water Conveyance Tunnel

Applying stiffness migration method,a 3D finite element mechanical model is established to simulate the excavation and advance processes.By using 3D nonlinear finite element method,the tunnel boring machine(TBM) excavation process is dynamically simulated to analyze the stress and strain field status of surrounding rock and segment.The maximum tensile stress of segment ring caused by tunnel construction mainly lies in arch bottom and presents zonal distribution.The stress increases slightly and limitedly in the course of excavation.The maximum and minimum displacements of segment,manifesting as zonal distribution,distribute in arch bottom and vault respectively.The displacements slightly increase with the advance of TBM and gradually tend to stability.

Numerical and Dimensional Analysis for Prediction of Line Heating Residual Deformations

Line heating process is a very complex phenomenon as a variety of factors affects the amount of residual deformations. Numerical thermal and mechanical analysis of line heating for prediction of residual deformation is time consuming. In the present work dimensional analysis has been presented to obtain a new relationship between input parameters and resulting residual deformations during line heating process. The temperature distribution and residual deformations for 6 mm, 8 mm, 10 mm and 12 mm thick steel plates were numerically estimated and compared with experimental and published results. Extensive data generated through a validated FE model were used to find co-relationship between the input parameters and the resulting residual deformation by multiple regression analysis. The results obtained from the deformation equations developed in this work compared well with those of the FE analysis with a drop in the computation time in the order of 100 (computational time required for FE analysis is around 7 200 second to 9 000 seconds and where the time required for getting the residual deformation by developed equations is only 60 to 90 seconds).

A New Finite Element Model with Manufactured Error for Additive Manufacturing

Additive manufacturing(AM),adding materials layer by layer,can be used to produce objects of almost any shape or geometry.However,AM techniques cannot accurately build parts with large overhangs,especially for the large features close to horizontal,hanging over the void.The overhangs will make the manufactured model deviate from the design model,which will result in the performance of the manufactured model that cannot satisfy the design requirements.In this paper,we will propose a new finite element(FE)analysis model that includes the manufacturing errors by mimicking the AM layer by layer construction process.In such FE model,an overhang coefficient is introduced to each FE,which is defined by the support elements in the lower layer.By mimicking the AM process from the bottom layer to the top layer,all the FE properties are updated based on their overhang coefficients,which makes the computational model be able to predict the manufactured model with manufacturing errors.The proposed model can be used to predict the performance of the AM objects in the design stage,which will help the designers to improve their design by the simulation results.

A 3-Node Co-Rotational Triangular Finite Element for Non-Smooth,Folded and Multi-Shell Laminated Composite Structures

Based on the first-order shear deformation theory,a 3-node co-rotational triangular finite element formulation is developed for large deformation modeling of non-smooth,folded and multi-shell laminated composite structures.The two smaller components of the mid-surface normal vector of shell at a node are defined as nodal rotational variables in the co-rotational local coordinate system.In the global coordinate system,two smaller components of one vector,together with the smallest or second smallest component of another vector,of an orthogonal triad at a node on a non-smooth intersection of plates and/or shells are defined as rotational variables,whereas the two smaller components of the mid-surface normal vector at a node on the smooth part of the plate or shell(away from non-smooth intersections)are defined as rotational variables.All these vectorial rotational variables can be updated in an additive manner during an incremental solution procedure,and thus improve the computational efficiency in the nonlinear solution of these composite shell structures.Due to the commutativity of all nodal variables in calculating of the second derivatives of the local nodal variables with respect to global nodal variables,and the second derivatives of the strain energy functional with respect to local nodal variables,symmetric tangent stiffness matrices in local and global coordinate systems are obtained.To overcome shear locking,the assumed transverse shear strains obtained from the line-integration approach are employed.The reliability and computational accuracy of the present 3-node triangular shell finite element are verified through modeling two patch tests,several smooth and non-smooth laminated composite shells undergoing large displacements and large rotations.

MATHEMATICAL MODEL OF MILLING TEMPERATURE AND TEMPERATURE FIELD ANALYSIS FOR THREE-DIMEN-SIONAL COMPLEX GROOVE MILLING INSERT

The mathematical model on the temperature of the waved-edge is constructedaccording to Jaeger's theory of moving solid and based on the used temperature model of the flatinsert. It is possible to forecast the milling temperature through programming. The comparableexperiments have been done between the two new three-dimension groove inserts (waved-edge insert,great edge insert) and flat fake insert. The theoretic forecast is in good agreement with theexperimental result. According to the cutting conditions, the boundary condition of finite elementanalysis on cutting temperature field is established, and the three-dimensional temperature fieldsof inserts with grooves are analyzed by FEM, so as to offer a reference basis for the design andoptimization of insert grooves.

Biomechanical analysis of an absorbable material for treating fractures of the inferior orbital wall

AIM:To investigate the biomechanical properties and practical application of absorbable materials in orbital fracture repair.METHODS:The three-dimensional(3D)model of orbital blowout fractures was reconstructed using Mimics21.0 software.The repair guide plate model for inferior orbital wall fracture was designed using 3-matic13.0 and Geomagic wrap 21.0 software.The finite element model of orbital blowout fracture and absorbable repair plate was established using 3-matic13.0 and ANSYS Workbench 21.0 software.The mechanical response of absorbable plates,with thicknesses of 0.6 and 1.2 mm,was modeled after their placement in the orbit.Two patients with inferior orbital wall fractures volunteered to receive single-layer and double-layer absorbable plates combined with 3D printing technology to facilitate surgical treatment of orbital wall fractures.RESULTS:The finite element models of orbital blowout fracture and absorbable plate were successfully established.Finite element analysis(FEA)showed that when the Young’s modulus of the absorbable plate decreases to 3.15 MPa,the repair material with a thickness of 0.6 mm was influenced by the gravitational forces of the orbital contents,resulting in a maximum total deformation of approximately 3.3 mm.Conversely,when the absorbable plate was 1.2 mm thick,the overall maximum total deformation was around 0.4 mm.The half-year follow-up results of the clinical cases confirmed that the absorbable plate with a thickness of 1.2 mm had smaller maximum total deformation and better clinical efficacy.CONCLUSION:The biomechanical analysis observations in this study are largely consistent with the clinical situation.The use of double-layer absorbable plates in conjunction with 3D printing technology is recommended to support surgical treatment of infraorbital wall blowout fractures.

Rectangular tunnel heading stability in three dimensions and its predictive machine learning models

Tunnel heading stability in two dimensions(2D)has been extensively investigated by numerous scholars in the past decade.One significant limitation of 2D analysis is the absence of actual tunnel geometry modeling with a considerable degree of idealization.Nevertheless,it is possible to study the stability of tunnels in three dimensions(3D)with a rectangular shape using finite element limit analysis(FELA)and a nonlinear programming technique.This paper employs 3D FELA to generate rigorous solutions for stability numbers,failure mechanisms,and safety factors for rectangular-shaped tunnels.To further explore the usefulness of the produced results,multivariate adaptive regression spline(MARS)is used for machine learning of big dataset and development of design equations for practical design applications.The study should be of great benefit to tunnel design practices using the developed equations provided in the paper.

旋提手法治疗神经根型颈椎病的应力及椎间孔形态学特征:三维有限元分析

背景:颈椎旋提手法治疗神经根型颈椎病的临床疗效确切,安全性较好,已广泛应用于临床中,但其对各节段颈椎椎体、椎间盘的生物力学及椎间孔面积的影响尚未系统阐明。目的:基于三维有限元分析技术开展相关研究,为旋提手法治疗神经根型颈椎病的效应机制提供数字化证据。方法:选择1例颈背部无异常疾患的志愿者,将其颈椎旋转至生理极限体位状态下的CT影像资料作为有限元模型材料。分别通过Mimics19.0、Geomagic Studio 2013、Hypermash 14.0及ANSYS Workbench 2020 R2软件完成对有限元模型的初步构建,基于文献完成颈椎各组织结构的弹性模量、弹性系数赋值,基于团队前期工作基础在模型上加载旋提手法的力学工况,分析颈椎旋提手法对C3-T1节段各椎体、椎间盘的力学参数及其对颈椎旋颈对侧椎间孔面积的影响。结果与结论:(1)旋提手法实施时,骨性结构的应力明显高于椎间盘等软组织;(2)手法操作时各颈椎椎体顶部椎板处应力高,底部椎板处应力低,关节突关节处及横突处应力较低;而各椎间盘顶部髓核处应力较低,底部与椎板连接处应力较高,但椎间盘应力最高点在顶部外侧纤维环处;(3)另外,加载提扳力后较之加载前,C6/C7节段椎间孔投影面积明显增大;(4)提示旋提手法操作具有改变颈椎本身应力结构的力学特征,并能扩大患者旋颈对侧的C6/7椎间孔面积,实现治疗神经根型颈椎病的目的。

3D打印个性化正畸装置对单根牙垂直移动的三维有限元分析

背景:基于牙齿垂直移动原理,通过数字化设计结合3D打印制作个性化正畸装置,可使个性化正畸装置与单个切牙形成支抗系统。借助微种植体绝对支抗作用可对单根牙在三维方向上进行精准控制。目的:基于生物力学原理设计了牙位11、12、21、22压低移动与伸长移动的个性化正畸装置,通过三维有限元法分析个性化正畸装置的安全性与其对牙齿的移动作用。方法:应用Mimics、Geomagic Wrap、SolidWorks和Ansys Workbench软件分别建立了上颌前牙区(牙位11、12、21、22)牙槽骨-牙周膜-上颌切牙-个性化悬臂微种植体-连接板-个性化托槽的三维有限元模型,每个牙位分别设置个性化压低移动与伸长移动正畸装置,加载300 g拉力或推力,分析个性化正畸装置各部件应力水平,计算牙齿位移与牙周膜应力分布情况。结果与结论:①个性化压低移动正畸装置受到的最大等效应力为162.90 MPa,个性化伸长移动正畸装置受到的最大等效应力239.57 MPa,二者的最大等效应力均位于个性化托槽加力附件的垂直部分,个性化正畸装置各部件所受等效应力均在屈服强度内,均具有良好的安全性;②在个性化正畸装置作用下,牙齿初始位移呈现整体压低或伸长移动的趋势,垂直方向上的位移远远超过在水平方向和矢状方向上的位移;牙周膜根尖或颈部位置出现等效应力峰值,在根尖区牙周膜出现等效应力集中区域;③结果显示,个性化正畸装置可以使牙位11、12、21、22近似达到压低移动或伸长移动的目的,初步证实了个性化垂直移动正畸装置的有效性。