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.

Design of Industrial Floors TR34 and Finite Element Analysis (Part 2)

The design of industrial floors will be presented in this paper. In the first part of this article the calculation methods of the TR34 British guideline will be discussed. In the second part the state of the art design methods using advanced finite element methods will be presented. The design itself may seem as slow considering the actual computer efficiency, however comparing the results to theoretical analysis and to designing methods, precision and economical nature of the method can be justified. A large number of foreign industrial floor designs were made by this method;some of them will be shown as reference at the end of the article.

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.

A MULTI-COUPLED FINITE ELEMENT ANALYSIS OF RESISTANCE SPOT WELDING PROCESS

A two-dimensional axisymmetric finite element model is developed to analyze the transient thermal and mechanical behaviors of the Resistance Spot Welding （RSW） process using commercial software ANSYS. Firstly a direct-coupled electrical-thermal Finite Element Analysis （FEA） is performed to analyze the transient thermal characteristics of the RSW process. Then based on the thermal results a sequential coupled thermo-elastic-plastic analysis is conducted to determine the mechanical features of the RSW process. The thermal history of the whole process and the temperature distribution of the weldment are obtained through the analysis. The mechanical features, including the distributions of the contact pressure at both the faying surface and the electrode-workpiece interface, the stress and strain distributions in the weldment and their changes during the RSW process, the deformation of the weldment and the electrode displacement are also calculated.

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.

Design and Structural Analysis of Flexible Wearable Chair Using Finite Element Method

The flexible wearable chair is like a light weight mobile exoskeleton that allows people to sit any-where in any working position. The traditional chair is difficult to move to different working locations due to its large size, heavy weight (~5 - 7 kg) and rigid structure and thus, they are inappropriate for workplaces where enough space is not available. Flexible wearable chair has a gross weight of 3 kg as it utilizes light-weight aluminium alloy members. Unlike the traditional chair, it consists of kinematic pairs which enable taking halts between continuous movements at any working position and thus, it is capable of reducing the risk of the physical musculoskeletal disorder substantially among workers. The objective of this paper is to focus on the mechanical design and finite element analysis (FEA) of the mechanism using ANSYS® software. In the present work, all the parts of the mechanism are designed under static load condition. The results of the analysis indicate that flexible wearable chair satisfies equilibrium and stability criterion and is capable of reducing fatigue during working in an assembly line/factory.

Finite Element Structural Analysis UsingObject-Oriented Techniques

This paper describes the object-oriented implementational method of finite element structural analysis, gives the basic concepts of the object-oriented method and objectoriented programming, develops a complete class hierarchy structure of object-oriented finite element structural analysis, and gives a part C+ + code description.

Finite Element Analysis of Dynamic Damper for CV Joint

Constant Velocity [CV] Joints are one of the most important components of front wheel drive axles. It is subjected to various stresses such as bending stress, shear stress and bearing stress. Apart from these stresses, it is also subjected to vibrations, due to out of balance tire or wheel and an out of round tire or wheel, or a bent rim. The main objective of this work is to reduce the stiffness of the damper, so that the damper can withstand within the required constraints [i.e. the forced frequency range of 80 Hz to 150 Hz]. The free vibrational and forced vibrational effects are investigated to predict the resonance phenomenon of the damper. Finite Element Analysis in ANSYS-11 software was performed to predict the dynamic behavior of the system under the required vibrational frequencies ranging from 80 Hz to 150 Hz at given loading conditions.

Stress Analysis of Buried Pipeline Installed by Horizontal Directional Drilling Using ANSYS Finite Element Software

In designing a horizontal directional drilling (HDD) pipeline project, designers face the challenge of determining the regions of maximum and minimum stresses on pipelines, ensuring the stability of the bore-hole from collapse and minimizing the stresses induced on the pipeline due to the bore-profile. This study analyses the stress induced on an HDD pipeline system using the ANSYS Version 18, mechanical APDL finite element (FE) software. The pipeline used as the case study was a gas transmission pipeline installed in south-west Nigeria. A macro-file for ANSYS Version 18, mechanical APDL used to model the pipeline was developed. The results showed that the maximum and minimum stresses induced on the HDD pipeline were at the top and bottom of the pipe, respectively;while the stresses on the sides were uniform (≈888 kg/cm2) all through the pipeline, irrespective of element number. The maximum stress occurred at the curvature point with the highest entry angle (10°), resulting in a maximum deflection at this point. The model stress validation performed by comparing results with theoretical solutions, both with respect to radius of curvature and internal pressure, showed percentage difference (errors) less than 10%. The cross sectional area validation showed a percentage difference of 0.059%.

Conventional versus minimally-invasive cervical discectomy for treatment of severe degenerative disease at C5-C6: a biomechanical comparison using a model of the full cervical spine and finite element analysis

The purpose of this study was to determine the dif-ferences in biomechanical responses of tissues in the cervical spine when pain and other problems secon-dary to severe disc degeneration disease are surgi-cally treated by conventional discectomy (CONDIS) compared to minimally-invasive discectomy (MIVDIS). A validated three-dimensional model of an intact, healthy, adult full cervical spine (C1-C7) (INT) was constructed. This model was then modified to create two models, one simulating each of the above-men-tioned two techniques for discectomy of the severely degenerated C5-C6 disc. For each of these three models, we used the finite element analysis method to obtain three biomechanical parameters at various tissues in the model, under seven different physio-logically relevant loadings. For each of the biome-chanical parameters, the results were expressed as relative change in its value when a specified combi-nation of simulated discectomy model and applied loading was used, with respect to the corresponding value in the intact model. We then computed the value of a composite biomechanical performance in-dex (CBPI) for CONDIS and MIVDIS models, with this value incorporating all of the aforementioned relative changes. We found that CBPI was marginally lower for MIVDIS model. This trend is the same as that reported for the relative complications rate and outcome measures following conventional and mini-mally-invasive discectomies in the lumbar spine. From a healthcare perspective, one implication of our finding is that minimally-invasive cervical discectomy should be considered an attractive option provided that detailed patient selection criteria are clearly de-fined and strictly followed.

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.

Relationships between femoral strength evaluated by nonlinear finite element analysis and BMD,material distribution and geometric morphology

Precisely quantifying the strength of the proximal femur and accurately assessing hip fracture risk would enable those at high risk to be identified so that preventive interventions could be taken.Development of better measures of femoral strength using the clinically

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.

Finite Element Simulation on the Spin-forming of the 3D Non-axisymmetric Thin-Walled Tubes

The roller movement trace for the 3D non-axisymmetric thin-walled tubes is a complex space curve. Besides the roller rotation caused by contact with the blank, the roller rotates around the workpiece together with the main spindle, and also moves simultaneously along the direction of the revolution radius. The method to correctly establish the finite element （FE） models of the metal spinning is based on the MSC. MARC software was introduced. The calculation formulas considering both the revolution and rotation of the roller were obtained by the mathematical deduction. The saving calculation points m should be a multiple of 4 for one revolution of the roller around the workpiece to obtain the maximum forming force for the spinning of the 3D non-axisymmetric thin-walled tubes. The simulation results conform well to the experimental ones for several spinning methods; the maximum error is less than ±15%.

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.

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.

Stochastic analysis of excavation-induced wall deflection and box culvert settlement considering spatial variability of soil stiffness

In this study, a three-dimensional (3D) finite element modelling (FEM) analysis is carried out to investigate the effects of soil spatial variability on the response of retaining walls and an adjacent box culvert due to a braced excavation. The spatial variability of soil stiffness is modelled using a variogram and calibrated by high-quality experimental data. Multiple random field samples (RFSs) of soil stiffness are generated using geostatistical analysis and mapped onto a finite element mesh for stochastic analysis of excavation-induced structural responses by Monte Carlo simulation. It is found that the spatial variability of soil stiffness can be described by an exponential variogram, and the associated vertical correlation length is varied from 1.3 m to 1.6 m. It also reveals that the spatial variability of soil stiffness has a significant effect on the variations of retaining wall deflections and box culvert settlements. The ignorance of spatial variability in 3D FEM can result in an underestimation of lateral wall deflections and culvert settlements. Thus, the stochastic structural responses obtained from the 3D analysis could serve as an effective aid for probabilistic design and analysis of excavations.

Virtual ballistic impact testing of Kevlar soft armor:Predictive and validated finite element modeling of the V0-V100 probabilistic penetration response

This works presents the first fully validated and predictive capability to model the V_0-V_(100) probabilistic penetration response of a woven fabric using a yarn-level fabric finite element model. The V_0-V_(100) curve describes the probability of complete fabric penetration as a function of projectile impact velocity. The exemplar case considered in this paper comprises of a single-layer, fully-clamped, plain-weave Kevlar fabric impacted at the center by a 17-gr, 0.22 cal FSP or fragment-simulating projectile. Each warp and fill yarn in the fabric is individually modeled using 3 D finite elements and the virtual fabric microstructure is validated in detail against the experimental fabric microstructure. Material and testing sources of statistical variability including yarn strength and modulus, inter-yarn friction, precise projectile impact location, and projectile rotation are mapped into the finite element model. A series of impact simulations at varying projectile impact velocities is executed using LS-DYNA on the fabric models, with each model comprising unique mappings. The impact velocities together with the outcomes(penetration, nonpenetration) are used to generate the numerical V_0-V_(100) curve which is then validated against the experimental V_0-V_(100) curve. The numerical Vi-Vrdata(impact, residual velocities) is also validated against the experimental Vi-Vrdata. For completeness, this paper also reports the experimental characterization data and its statistical analysis used for model input, viz. the Kevlar yarn tensile strengths, moduli, and inter-yarn friction, and the experimental ballistic test data used for model validation.