Three-dimensional Finite Element Research Articles

Influence of connection design and material properties on stress distribution and fatigue lifetime of zygomatic implants: A finite element analysis

Zygomatic implants (ZIs) were developed as a graftless alternative to rehabilitate severely reabsorbed maxillae. This study aims to employ three-dimensional finite element analysis (FEA) to simulate the impact of external hexagonal implant connection (EHC) and internal hexagonal implant connection (IHC) on the stress distribution and fatigue lifetime within the ZI systems using parameters defined in ISO 14801:2016. Two ZI assemblies (Nobel Biocare and Noris Medical) were scanned in a micro-CT scanner and reconstructed using Nrecon software. Three-dimensional models were generated by Simpleware ScanIP Medical software. All models were exported to FEA software (ABAQUS) and subsequently to a fatigue analysis software (Fe-safe). A compressive 150 N load was applied at a 40° angle on the cap surface. A 15 Hz frequency was applied in the in silico cyclic test. The implant components had material properties of commercially pure grade 4 titanium (CPTi) and Titanium-6Aluminum-4Vanadium alloy (Ti64). Von Mises stress data, contour plots, and fatigue limits were collected and analyzed. EHC models exhibited higher peak stresses in implant components for both materials compared to IHC models. However, simulated bone support results showed the opposite trend, with higher stresses on IHCthan EHC models. The fatigue analysis revealed that assemblies with both designs exceeded ISO 14801:2016 number of cycles limits using Ti64, while CPTi groups exhibited comparatively lower worst life-repeats. In conclusion, ZIs with IHC were found to have a more homogeneous and advantageous stress distribution within both materials tested. Ti64 demonstrates a prolonged service life for both design connections.

Numerical study of steady-state dynamic characteristic of non-pneumatic tire with local structural damage

Non-pneumatic tires (NPTs) fundamentally avoid the risk of tire blowout of traditional pneumatic tires, but the overall and key component performance inevitably degrades due to factors such as high temperature, variable load, variable working conditions and impact in its service process. Once this leads to failure, it will significantly impact on vehicle safety. To this end, this paper studied the influence of local structural damage on the dynamic response of NPTs, which lays the foundation for realizing health monitoring of intelligent non-pneumatic tires (INPTs). Firstly, a three-dimensional nonlinear finite element model of NPTs was established, and the failure locations of NPTs were determined by fracture mechanics and maximum strain energy density; Secondly, the influence of structural damage on the static and dynamic performance of NPTs was analyzed; Finally, the sensitivity of the acceleration signal and sensor position of the tire inner liner to the local structural damage was studied. The research results show that structural damage will cause the stress of the spokes and shear layer to increase, and as the number of broken spokes increases, the shear layer will bear a larger load. Compared with the circumferential and lateral acceleration, the radial acceleration has the highest sensitivity to the damage of NPTs. The sensor closest to the damage location is the most sensitive to the damage. The research results provide a reference for the structural optimization design and health monitoring of INPTs.

Finite Element Modeling of the Dynamic Response of Plywood.

Modeling the dynamic properties of wood and wood-based composites is a challenging task due to naturally growing structure and moisture-dependent material properties. This paper presents the finite element modeling of plywood panels' dynamic properties. Two panels differing in thickness were analyzed: (i) 18 mm and (ii) 27 mm. The developed models consisted of individual layers of wood, which were discretized using three-dimensional finite elements formulated using an orthotropic material model. The models were subjected to an updating procedure based on experimentally determined frequency response functions. As a result of a model updating relative errors for natural frequencies obtained numerically and experimentally were not exceeding 2.0%, on average 0.7% for 18 mm thick panel and not exceeding 2.6%, on average 1.5% for 27 mm thick panel. To prove the utility of the method and at the same time to validate it, a model of a cabinet was built, which was then subjected to experimental verification. In this case, average relative differences for natural frequencies of 6.6% were obtained.

Development and experimental verification of an advanced 3D elastoplastic progressive damage model for composite materials

This study develops a sophisticated three-dimensional elastoplastic progressive damage finite element analysis (FEA) model for stiffened composite materials under axial compression. The model integrates the Hashin and Ye criteria to capture the inception of damage within the composite and introduces a cohesive zone model (CZM) to simulate debonding between stiffeners and panels. Furthermore, damage evolution based on energy release rate and plastic flow rules grounded in Hill’s criteria are employed to thoroughly explore the progressive failure behavior of stiffened composite materials. Implemented within Abaqus/Explicit through user-defined ‘VUMAT’ subroutines, the model ensures computational accuracy and reliability. Comparative evaluations against standard linear elastic progressive damage models demonstrate a statistically significant improvement of at least 5% in predicting ultimate load capacities. Additionally, this model accurately captures different failure modes, providing a nuanced understanding of stiffened composite material failure under compressive loads.

Numerical investigation into the influence of interfacial transition zone on elastic modulus, creep, shrinkage of concrete

Interfacial transition zone (ITZ) is a weak part around aggregate in concrete, which could threaten the macroscopic properties of concrete such as strength, stiffness, deformation, and durability. However, it is difficult to quantify the ITZ effect on the macroscopic behavior of concrete only by experimental methods due to the tiny thickness of ITZ. This study employed finite element (FE) method to investigate the ITZ effect on the elastic modulus, creep, and shrinkage. A mesoscale three-phase FE model consisting of homogeneous matrix phase, aggregates and ITZ around them was developed for concrete in this study. The influences of the FE model dimension and the meshing size on the calculated results were first discussed. Then the influence of the thickness, viscoelasticity and shrinkage of ITZ on the macroscopic properties of concrete were evaluated based on the FE calculation and grey relation analysis. Finally, the measured elastic modulus, creep and shrinkage of concrete with different contents of clay on coarse aggregates were used to validate the FE model. It is found that the calculated results of the two-dimensional (2D) FE model are in good agreement with those of the three-dimensional (3D) FE model under the linear viscoelastic assumption. The meshing size has a minor effect on the calculated results when it is 0.2–2 times of the ITZ thickness. The influence of ITZ properties on the macroscopic properties of concrete will become more significant as the ITZ thickness increases. Moreover, the elastic modulus of ITZ has a greater impact on the effective properties of concrete than the thickness, creep, and shrinkage of ITZ. The findings in study can contribute to a better understanding of ITZ effect on the macroscopic properties of concrete and provide references to the reasonable consideration of ITZ in numerical studies.

Seismic response assessment of tapered piles in sandy soils: A numerical investigation

For the same volume, tapered (non-prismatic) piles have premium advantages over conventional (prismatic) cylindrical piles with respect to the axial and lateral load resistance and, thereby, are considered an economical alternative to prismatic piles. While previous investigations have studied the behavior of tapered piles subjected to different dynamic loading conditions, their seismic response has not yet received a comprehensive examination. This study aims to investigate the seismic response of tapered piles in sandy soils using non-linear three-dimensional (3D) finite element analyses. A total of 3084 cases were studied using MIDAS software to consider a wide range of parameters, including variations in soil density (Dr), taper angle (θ°), slenderness ratio (L/DAvg), and different earthquake recorded data. A parametric analysis was also undertaken to investigate the behavior of the pile lateral displacement, pile bending moment, and frictional resistance. The findings of the study reveal that tapering enhances the lateral stiffness of the pile, resulting in reduced lateral deflections and decreased sensitivity to changes in L/DAvg and ground motion intensity. Furthermore, under higher ground motion intensity, improved pile stiffness through tapering contribute to an increase in bending moments experienced by the pile. Another significant observation is the substantial improvement in frictional resistance due to soil densification. This emphasizes the critical importance of considering the soil behavior and its density in the design of pile foundations, particularly when aiming to withstand seismic loads effectively. These insights offer valuable guidance for designing pile foundations that can reliably endure the dynamic forces present during seismic events.

Comparative Evaluation of Stress Distribution and Deformation in Class II Cavities Restored With Two Different Biomimetic Restorative Materials: A Three-Dimensional Finite Element Analysis.

Aim This study aims to comparatively evaluate the stress distribution and deformation in class II cavities restored with EverX posterior and ActivaTMBioactive restorative material using three-dimensional finite element analysis. Methods On a mandibular first molar tooth model, features of a class II mesio-occlusal cavity were designed. Following the application of EverX Posterior and ActivaTMBioactive, the model was subjected to a 600 N static occlusal load. Three-dimensional finite element analysis was used to analyze the distribution of occlusal stresses among the remaining components and restorative materials. von Mises and maximum principal stresses were evaluated and compared. Results von Mises and maximum principal stresses of ActivaTMBioactive components were higher compared to the control group and EverX Posterior model. ActivaTMBioactive restorations had absorbed more stresses within the material rather than distributing it to enamel and dentin than EverX Posterior material. Also, deformation was extremely high in the ActivaTMBioactive model, which may be because of less stiffness of the material. Conclusion EverX posterior has a higher capacity to bear stresses in stress-bearingareas than ActivaTMBioactive Restorative.

Biomechanical effects of mandibular deviation on the temporomandibular joint in patients with mandibular prognathism under incisal occlusion

PurposeThe aim of this study was to investigate the biomechanical effects of mandibular deviation on the TMJ in patients with mandibular prognathism before and after orthognathic surgery using three-dimensional finite element analysis. MethodsEight patients with mandibular prognathism without deviation, eight patients with mandibular prognathism with deviation and sixteen normal subjects were recruited. Three-dimensional models of the maxillofacial were reconstructed using MIMICS. Nine muscle forces were used to simulate incisal occlusion and contact was used to simulate fossa-disc-condyle interactions. ResultsBefore surgery, the stress in the TMJ was generally greater in the Pre-MD&MP group than in the Pre-MD group; it was much greater in both groups than in the control group, ranging from about 2 to 12 times as great in the Pre-MD group and from about 5 to 64 times as great in the Pre-MD&MP group. After orthognathic surgeries, the stresses in the Post-MP&MD were significantly reduced by approximately 21.7 % to 93.4 %. And in the Post-MP group, the stresses were reduced by approximately 1.4 % to 51.1 %. ConclusionMandibular deviation exacerbated the abnormal stress distribution in the TMJ of patients with mandibular prognathism. Orthognathic surgeries could improve the stress distribution in patients with mandibular prognathism (with and without deviation). TMD was closely related to the stress levels of the TMJ.

The effects of optimized microstructured surfaces on bond strength and durability of NPJ-printed zirconia

ObjectivesThis study was to investigate the effects of optimized microstructured surfaces on bond strength and bond durability of the latest nanoparticle jetting (NPJ)-printed zirconia. MethodsZirconia microstructured surfaces with different geometries and void volume were analyzed through three-dimensional finite element analysis for surface micromorphology optimization. Zirconia disks and cylinders were additively manufactured by an NPJ 3D printer (N = 128). They were randomly divided into four groups based on surface micromorphology optimization and airborne-particle abrasion (APA) treatment before they were bonded using 10-methacryloloxydecyl dihydrogen phosphate (MDP) containing resin cement (Clearfil SA luting cement). The shear bond strengths (SBSs) were tested before and after 10,000 thermocycles and were analyzed by one-way ANOVA analysis. Failure modes were determined by optical microscopy. Zirconia surfaces were analyzed with X-ray diffraction, scanning electron microscopy, and three-dimensional interference microscopy. ResultsThe optimized microstructured surface was characterized by circular microstructures with 60 % void volume, about 20 µm of depths, about 10 µm of undercuts, and consistent beam widths. The optimized microstructured surface combined with APA treatment and MDP-containing resin cement possessed the highest SBSs both before and after thermocycling aging (P＜0.05). The greater reductions of zirconia bond strengths occurred when the zirconia were not treated with APA (P＜0.05). SignificanceThe optimized microstructured zirconia surface with circular microstructures and 60 % void volume fabricated by the latest NPJ printing technology could greatly enhance the zirconia bond strength and durability in combination with APA treatment and application of MDP-containing resin cement, which might be promising for adhesively bonded indirect restorations of NPJ-printed zirconia.

Evolution of Confinement Stress in Axially Loaded Concrete-Filled Steel Tube Stub Columns: Study on Enhancing Urban Building Efficiency

In the context of green building and sustainable urban development, understanding the mechanical behavior of structural components like concrete-filled steel tube (CFST) columns is crucial due to their improved load-bearing capacity, energy efficiency, and optimized material usage, which enhance structural resilience and sustainability. This research addresses the complex development of confining stress and its impact on the concrete core (CC) behavior within these columns, which are essential for urban infrastructure. Through extensive numerical studies, this study proposes a model to estimate the confining stress in axially loaded CFST short columns. Study findings reveal that the confinement effectiveness is influenced by variables such as compressive strength (CS) of the concrete, cross-sectional shape, and depth-to-wall thickness percentage. Additionally, the confinement is also significantly affected by the yield strain of steel εy/εc to the peak strain of unconfined concrete εc. A three-dimensional finite element model (FEM) was built for the simulation of the columns’ nonlinear behavior and was rigorously validated against experimental data. This model aids in the design and implementation of more efficient and resilient urban structures, supporting the principles of sustainable construction. The study underscores the importance of structural integrity in sustainable urban development and provides valuable insights for improving the design of green building materials.

Progressive collapse resistance of self-resilient composite frames under fire conditions

Self-resilient composite frames exhibit robust resistance to progressive collapse and excellent seismic performance. However, their behavior under fire conditions is inadequately explored. This study investigates the collapse resistance mechanisms and response of self-resilient composite frame substructures (SRCFS) under fire-induced progressive collapse. A three-dimensional finite element model was developed, incorporating temperature-dependent thermal and mechanical properties of concrete and steel. Sequentially coupled thermal-mechanical analysis was conducted using the ABAQUS/Explicit solver in a quasi-static manner. Existing experimental data at both ambient and high temperatures were used to validate the proposed model. The study analyzed the contributions of steel beams and strands to vertical resistance, examining their impact on the collapse-resistance mechanism through a constant load-heating transient loading scheme. Structural responses under ambient and fire conditions were compared, focusing on different damage mechanisms. The effects of initial prestressing of strands, angle gauge length, angle thickness, and various heating curves on collapse resistance were evaluated. Findings suggest that the ISO-834 standard's failure criterion is overly conservative for SRCFS, as it neglects the enhanced load-carrying capacity provided by tensile catenary action. This mechanism offers additional escape time, enhancing personnel safety. The study identifies factors influencing the performance of SRCFS against progressive collapse under fire and proposes design recommendations.

Comparative analysis of internal fixation modalities for PER type IV ankle fractures: a finite element study

BackgroundThere are many options for the surgical treatment of pronation external rotation (PER) type IV ankle fractures, including the use of fibular plates and screws, the aim of this study was to investigate the biomechanical stability and safety of different internal fixation methods for PER type IV ankle fractures via finite element analysis.MethodsA three-dimensional finite element model of the ankle joint and the whole foot of a healthy 26-years-old adult male was established and validated based on computed tomography images of his lower limb, and a computer-aided design was used to produce a PER type IV ankle fracture and plate and screw model. Four different internal fixation modes were simulated, including a (all ankle fixation—utilizing a fibular plate and screws for comprehensive stabilization of the ankle), b (inferior tibiofibular joint fixation + all ankle fixation), c (inferior tibiofibular joint fixation + unfixed anterior ankle), and d (inferior tibiofibular joint fixation + unfixed anterior and posterior ankles). The results of the four different fixation methods were compared via finite element analysis, and the von Mises stresses. The displacements of the four different fixation methods were analyzed as the output indices.ResultsThere were no significant differences between the results of using fibular plates and screws and the displacement of fracture breaks among the four internal fixation modalities. The von Mises stress in the tibiotalar joint, median ankle, posterior ankle, and anterior ankle was minimized in the working condition of d, d, b, and d respectively. The von Mises stress in the fibular plate and screws was minimized in the working condition of a. The von Mises stress in the distal fibula was minimized in the working condition of a. However, the stress was mainly concentrated at the attachment point of the inferior tibiofibular anterior ligament, and in the working condition with inferior tibiofibular joint fixation, the stress was significantly concentrated in the inferior tibiofibular joint screw in all the fibular plates and screws.ConclusionsThe results of this study demonstrate the feasibility of using finite element analysis to compare the biomechanical stability and safety of four configurations of fibular plates and screws for treating PER type IV ankle fractures. All four modalities provided comparable biomechanical stability and safety, showing no significant differences. However, the current limitations of the finite element analysis methodology preclude specific clinical inferences. Further refinement of the methodology in future studies is necessary to enable reliable clinical applications.

Dentoalveolar effects of open-bite correction with the dual action vertical intra-arch technique: A finite element analysis.

To evaluate tooth displacement and periodontal stress generated by the dual action vertical intra-arch technique (DAVIT) for open-bite correction using three-dimensional finite element analysis. A three-dimensional model of the maxilla was created by modeling the cortical bone, cancellous bone, periodontal ligament, and teeth from the second molar to the central incisor of a hemiarch. All orthodontic devices were designed using specific software to reproduce their morpho-dimensional characteristics, and their physical properties were determined using Young's modulus and Poisson's coefficient of each material. A linear static simulation was performed to analyze the tooth displacements (mm) and maximum stresses (Mpa) induced in the periodontal ligament by the posterior intrusion and anterior extrusion forces generated by the DAVIT. The first and second molars showed the greatest intrusion, whereas the canines and lateral incisors showed the greatest extrusion displacement. A neutral zone of displacement corresponding to the fulcrum of occlusal plane rotation was observed in the premolar region. Buccal tipping of the molars and lingual tipping of the anterior teeth occurred with intrusion and extrusion, respectively. Posterior intrusion generated compressive stress at the apex of the buccal roots and furcation of the molars, while anterior extrusion generated tensile stress at the apex and apical third of the palatal root surface of the incisors and canines. DAVIT mechanics produced a set of beneficial effects for open-bite correction, including molar intrusion, extrusion and palatal tipping of the anterior teeth, and occlusal plane rotation with posterior teeth uprighting.

Three-dimensional unstructured adaptive finite element forward modeling simulation of magnetotellurics based on magnetic vector potential-electric scalar potential

Three-dimensional unstructured adaptive finite element forward modeling simulation of magnetotellurics based on magnetic vector potential-electric scalar potential

Analyzing the effect of undermining on suture forces during simulated skin flap surgeries with a three-dimensional finite element method

Skin flaps are common procedures used by surgeons to cover an excised area during the reconstruction of a defect. It is often a challenging task for a surgeon to come up with the most optimal design for a patient. In this paper, we set up a simulation system based on the finite element method for one of the most common flap types — the rhomboid flap. Instead of using the standard 2D planar patch, we constructed a 3D patch with multiple layers. This allowed us to investigate the impact of different undermining areas and depths. We compared the suture forces for each case and identified vertices with the largest suture force. The shape of the final suture line is also visualized for each case, which is an important clue when deciding on the most optimal skin flap orientation according to medical textbooks. We found that under the optimal undermining setup, the maximum suture force is around 0.7 N for top of the undermined layer and 1.0 N for bottom of the undermined layer. When measuring difference in final suture line shape, the maximum normalized Hausdorff distance is 0.099, which suggests that different undermining region can have significant impact on the shape of the suture line, especially in the tail region. After analyzing the suture force plots, we provided recommendations on the most optimal undermining region for rhomboid flaps.

Analysis of stress field in the head area of the Three Gorges Reservoir based on coupled fluid-solid theory

The Three Gorges Reservoir, one of the largest water conservation system in the world, has been of keen interest to scientists globally since its impoundment. After construction of the dam, there has been a significant increase in seismic activity in the head area of the reservoir. It is generally accepted that earthquakes in this region are predominantly caused by the Jiuwanxi and Xiannvshan faults. This study focused on the stress changes occurring in the research area. A three-dimensional finite element model of the reservoir area was constructed using the geological structure and digital ground elevation data of the reservoir area. The fluid-solid coupling theory was applied to calculate the dynamic spatial changes in pore pressure and Coulomb stress in the faults and surrounding rocks during reservoir impoundment. The findings indicated that the added head pore water pressure at the bottom of the reservoir had a maximum impact range of approximately −2800 m on the surrounding rock, whereas the Xiannvshan and Jiuwanxi faults had a maximum diffusion range of approximately −4300 m. Rock permeability also played a significant role in the water storage process. During the 1 56 m water impoundment stage, owing to rapid water storage activity, stress could not be transmitted to both sides in a timely manner, resulting in the formation of an extreme stress change zone at −4000 m inside the fault. This may have been the reason for the frequent earthquakes during this stage. The 17 5 m cycle water storage stage also exhibited a significant degree of seismicity, potentially attributable to the long-term infiltration of reservoir water and accumulation of stress in the previous stage. The stress in the study area at the four stages are in a process of accumulation-release-accumulation-release.

Analysis of the behavior of piled foundations in unconventional soil

Building foundations using pile groups generally assume that the piles absorb the entire load of the structure. However, it is known that pile caps rest not only on the piles but also on the surrounding soil. These foundations involve great complexity in the phenomenon of load sharing, especially when they are situated on unconventional soils, such as tropical climate soils. To study this behavior, static load tests were conducted on real foundations. These tests were also used to calibrate a three-dimensional finite element numerical model (FEA-3D) using Rocscience's RS3 software. The results showed that the settlement corresponding to the allowable load of the pile groups was significantly lower, ranging from 22 % to 73 %, compared to pile caps that exert contact pressure on the soil. For pile groups, the efficiency ranged between 92 % and 109 %, while for pile caps with contact, it ranged from 132 % to 178 %. In this sense, it becomes evident that estimating the settlement of the pile group based on the conventional design, which does not consider the contact between the cap and the soil, can significantly underestimate the settlement, and provide an overly optimistic analysis of the soil-structure interaction.

Identification of hanger damage based on deflection influence matrix

Tied-arch bridges, a vital component of modern infrastructure, are susceptible to various forms of damage, particularly hangers. The detection and identification of such damages are crucial for maintaining structural integrity and safety. However, traditional methods face challenges in terms of accuracy and efficiency. This study aims to develop a refined method for hanger damage identification in tied-arch bridges, to address the limitations of existing techniques. By focusing on deflection changes at the anchoring points between the hangers and tie beams, we sought to enhance the precision of damage detection. We propose an innovative approach based on the concept of influence lines, introducing the ‘generalized deflection difference influence line'and the ‘deflection difference influence matrix’. Then proposed a new identification index for identifying the damaged hanger after matrix. An actual tied-arch bridge was used to validate the proposed approach. A detailed three-dimensional finite element model of the bridge was developed and calibrated using dynamic and static response data. Thirty different hanger-damage conditions were simulated to evaluate the effectiveness of the proposed method. Our findings reveal that the deflection difference influence matrix offers more detailed and comprehensive information on bridge distribution points than traditional methods. Our method proved effective in identifying hanger damage, irrespective of its location on the bridge. In additionally, the identification efficiency of the method can be improved by adjusting the magnitude of the applied load, with larger loads amplifying the detectability of damage. This study highlights the potential of the deflection difference influence matrix to revolutionize hanger damage identification for tied-arch bridges. Its adaptability, accuracy, and efficiency are significant advancements over existing methods. This study successfully demonstrates an innovative and reliable method for hanger damage identification.

Post-earthquake fire performance of partially encased composite steel and concrete columns

Considering the potential severity of losses caused by post-earthquake fires compared to earthquakes alone, this study employed a three-dimensional Finite Element (FE) model established in ABAQUS to investigate the response of partially encased composite (PEC) steel and concrete columns under post-earthquake fire, serving as preliminary analysis for subsequent experimental studies. The validity of the FE model was verified using results from existing seismic and fire loading tests. In simulating the hysteresis behavior of PEC columns, this study evaluated the impacts of concrete discrete cracking and mechanical models on the simulation outcomes, thereby refining the modeling technique. A typical PEC column was selected as the focus of this investigation. Sequential earthquake-fire simulations were conducted using a data transfer technique to assess the fire resistance performance of different seismic damaged columns with different axial compression load ratios. The results indicate that the behavior of damaged PEC columns evolves based on initial damage, and the failure exhibits overall bending accompanied by local buckling of the steel flange. As the axial load ratio increases, the fire resistance time of seismic-damaged PEC columns significantly decreases. When the axial compression ratio is increased from 0.2 to 0.5, the fire resistance time of columns with damage factor D within 0.6 decreases by 59.07% on average. This sensitivity to axial pressure is particularly pronounced when earthquake damage is severe. Moreover, it finds that the lateral drift ratio under seismic action should be controlled within 2% to ensure that the post-earthquake fire resistance will not be significantly reduced.

Three-dimensional simulation on phase change heat transfer of moving spherical particle

Investigating the motion process of paraffin particle under a mesoscopic scale is very helpful to master the heat transfer mechanism of the mushy zone in solid–liquid phase change. Prior researches predominantly focused on either particle motion or phase change, without accounting for their coupling. In this study, the three-dimensional finite element method is used to solve the motion and heat transfer of the particle, taking into account the coupling of motion and phase change. The moving grid is employed to trace the interface of moving particles and liquid using the Arbitrary Lagrangian-Euler method, which is able to calculate flow and heat transfer in fluid region and allow grid deformation in solid region. The results indicate that uneven fluid velocity around the particles results in non-uniform temperature gradients on particle’s surface. In areas with larger temperature gradient, it melts more quickly and its deformation is obvious. The maximum is located at particle orientation angle 30°. Changing the size and direction of fluid velocity will reinforce or inhibit the movement of particle. The changing time of motion direction when fluid velocity is −0.03 m/s shorten 2.7 times than when fluid velocity is −0.01 m/s. The initial temperature of the surrounding fluid has a positive relationship with the particle melting rate. The channel width affects significantly the motion and melting of particle. The initial position of particle has an important effect on it’s falling movement, and the closer is close to the wall of channel, the greater the displacement of left and right movement.