Three-dimensional Finite Element Model Research Articles

Cyclic performance of RC frame joints with drilled openings

The construction industry commonly employs reinforced concrete (RC) beams with utility pipe openings to facilitate the passage of pipes through structures. These openings can be made through pre-construction or post-construction processes. Pre-construction openings in RC members are included in the design, and in most cases, designers include reinforcements around the opening to deal with the stress concentration and failure mechanism caused by the presence of the opening. However, post-construction or drilled openings are added after the structure is built, potentially disrupting established load paths, and developing stress concentrations that are not anticipated by the designer. In this study, three-dimensional nonlinear numerical finite element models (FEM) models were developed and validated according to four sets of experimental data. These models formed the basis for a parametric study exploring different design parameters. Parameters such as the width and depth of the specimen and the size of the opening were investigated. The results suggest that the presence of the opening leads to a brittle behavior and a sudden drop in the strength and stiffness of the specimen after the ultimate load. Increasing the depth of the beam results in more brittle behavior in specimens with openings and decreased the drift ratio where the specimen failed by up to 75%. However, increasing the width of the specimen increases the ductility of these specimens. Finally, opening size was found to lead the specimen to behave more brittle in all cases, and decrease the maximum drift ratio from 5% to 1.5%, and decrease the load capacity in 1.5% drift ratio by up to 61%. This study will help engineers and designers to understand the effect of beam-column joints with post-construction openings under cyclic loading, whether these openings are intentional or by human error.

Study on the Seismic Response of a Water-Conveyance Tunnel Considering Non-Uniform Longitudinal Subsurface Geometry and Obliquely Incident SV-Waves

The longitudinal seismic response characteristics of a shallow-buried water-conveyance tunnel under non-uniform longitudinal subsurface geometry and obliquely incident SV-waves was studied using the numerical method, where the effect of the non-uniform longitudinal subsurface geometry due to the existence of a local one-sided rock mountain on the seismic response of the tunnel was focused on. Correspondingly, a large-scale three-dimensional (3D) finite-element model was established, where different incidence angles and incidence directions of the SV-wave were taken into consideration. Also, the non-linearity of soil and rock and the damage plastic of the concrete lining were incorporated. In addition, the wave field of the site and the acceleration response as well as damage of the tunnel were observed. The results revealed the following: (i) a local inclined subsurface geometry may focus an obliquely incident wave due to refraction or total reflection; (ii) a tunnel in a site adjacent to a rock mountain may exhibit a higher acceleration response than a tunnel in a homogeneous plain site; and (iii) damage in the tunnel in the site adjacent to a rock mountain may appear in multiple positions, and the effect of the incidence angle on the mode of compressive deformation and damage of the lining is of significance.

Finite element analysis of the thermal and thermo-mechanical coupling problems in the dry friction clutches using functionally graded material

Abstract The friction clutch is an important machine part in the torque transmission system, where it is responsible for gently transmitting the power from the engine to the gearbox and damping the tensional vibration. At the beginning of the engagement between the clutch’s elements, a large amount of heat is released due to frictional sliding between the elements of the clutch system, generating very high temperatures. In this work, the two- and three-dimensional finite element models were developed from scratch to explore the clutch disk’s temperature and contact pressure behaviors using an alternative frictional material such as functionally graded material (FGM). The two assumptions for the thermal loads were considered: uniform wear and uniform pressure assumptions. The thermal–structural problem for the new frictional material was solved numerically and then presented as the temperature distributions and the contact pressure at any instant of sliding time. Also, a comparison was made between the clutch system behaviors when using FGM with different frictional materials. The results showed that the FGM has superior results to the other available frictional materials.

Shear failure analysis of UHPC-NC composites interface with U-shaped studs

Given the extensive utilization of composite structures of ultra-high performance concrete (UHPC) and normal concrete (NC) in engineering applications, there has been a growing emphasis on studying the shear failure and mechanical properties analysis of the UHPC-NC interface. U-shaped studs assume a crucial role in reinforcing the UHPC-NC interface. A three-dimensional finite element model was constructed with different parameters to investigate the shear failure and mechanical properties of UHPC-NC incorporating U-shaped studs. The failure mode and the influences of diameter, internal span, stud number and concrete compressive strength on the shear characteristics of the UHPC-NC interface were analyzed. The numerical analysis indicates the final failure mode of the UHPC-NC interface with U-shaped stubs is the fracture of the U-shaped stud. The diameter and internal span of studs influence the interface shear capacity, while interface shear capacity is not sensitive to the compressive strength of UHPC and NC. Moreover, a calculation formula for the shear capacity of the UHPC-NC interface with U-shaped studs was proposed, considering the influence of stud diameter and internal span. An exponential function is recommended for predicting the interface shear slip and the stiffness corresponding to 0.55 times the shear capacity is adopted as the elastic shear stiffness. The interface shear slip model and the expression for elastic shear stiffness were established, providing a reference for approximate calculation without experimental data.

Experimental and numerical studies on influences of wedge reinforcement on seismic performance of loose penetrated mortise-tenon joints

Wooden wedges significantly affect the rotational stiffness and load-bearing capacity of penetrated mortise-tenon joints (PMTJ). To clarify the role of wedges on the seismic performance of loose penetrated mortise-tenon joints, three types of full-scale PMTJs, including reference, loose, and wedge-reinforced joints, were tested under quasi-static loading. The deformation characteristics, hysteresis curves, rotational stiffness degradation behavior, and energy dissipation of PMTJs with wedge were analyzed. Three-dimensional nonlinear finite element models were established to simulate the wedge effects, and the stress states of joint were analyzed. The results indicated that the ultimate bending moment capacity of reference and loose joints was reached when the splitting failure of tenon occurred, while the wedge-reinforced joint reached such capacity when the transverse compression failure occurred at the wedges. The maximum bending moment of PMTJ with a 10 mm gap and a 25 mm gap were 18.3 % and 55.3 % lower than that of reference joint, while maximum bending moments of wedge-reinforced joint with a 10 mm gap and a 25 mm gap were 12.6 % and 17.03 % larger than that of reference joint, respectively. Numerical simulation showed that wooden wedge reduced stress concentration in the tenon and alleviated compression damage of the tenon. Higher load-bearing capacity and rotational stiffness were observed in joints when hardwood wedges were used, compared with that of soft wood. Wedge reinforcement could effectively improve the load-bearing capacity of loose PMTJ. These results yielded valuable data for the performance evaluation and retrofitting of penetrated mortise-tenon joints with wedges in traditional timber structures.

Simulation of fretting wear in steel wires under variable coefficient of friction and variable wear coefficient

In this paper, the fretting wear behaviour of steel wires working in coal mining technology is studied numerically. In past studies, the fretting process of steel wires was carried out numerically considering that the coefficient of friction (COF) and wear coefficient (WC) are constant parameters. However, it has been noticed experimentally that COF increases up to a certain number of fretting cycles and then becomes constant, i.e. a steady-state stage, depending on the loading conditions. This increase in COF during the fretting process is also known as the running-up stage. The fretting wear model is modified to evaluate the influence of the varying coefficient of friction (VCOF), which is associated with the variable wear coefficient (VWC), so the influence of VWC is also considered. The subroutine UMESHMOTION used to implement the wear law is also modified to study the effect of VCOF and VWC. Therefore, in this study, the numerical results of a three-dimensional finite element (FE) model are compared, with analytical results of contact area and contact stresses, and with experimental results of peak wear depth. After validating the FE model, the wear scar, the increasing wear depth, wear volume, and the decreasing contact stress with increasing fretting cycles are determined numerically considering VCOF and VWC using cycle jump approach. The energy dissipation effect of frictional force and fretting amplitude is also studied for varying interaction properties of fretting wear models. The numerical simulations are performed by considering both elastic and plastic material properties to analyse the influence of varying interaction properties on fretting wear models at the running-up stage. The results indicate that the VWC model exhibits comparable impacts on both the elastic and plastic models. The results also show that the VWC fretting wear model leads to higher wear scar, wear volume, and wear depth values at the running-up stage as well as at the steady state stage, which are close to the experimental data.

Investigation of the Role of Osteoporotic Vertebra Degeneration on the Stability of the Lumbar Spine: In Silico Modelling under Compressive Loading.

An evaluation of the impact of osteoporosis on loss of spinal stability, with or without intervertebral disc degeneration, using computational analysis is presented. The research also investigates the correlation between osteoporosis and intervertebral disc degeneration. Three-dimensional finite element models of human lumbar spine segments were used to assess the influence of osteoporosis on spinal stability. Five different models of age-related degeneration were created using various material properties for trabecular bone and intervertebral discs. Calculation results indicate that in a spine with osteoporosis, the deformation of the intervertebral discs can increase by more than 30% when compared to a healthy spine. Thus, intervertebral disc deformation depends not only on the degree of degeneration of the discs themselves, but their deformation is also influenced by the degree of osteoporosis of the vertebrae. Additionally, the load-bearing capacity of the spine can decrease by up to 30% with osteoporosis, regardless of the degree of intervertebral disc deformation. In conclusion, osteoporosis can contribute to intervertebral disc degeneration.

Deformation of Existing Shield Tunnel Adjacent to Deep Excavations: Simulation and Monitoring Analysis

Deep excavations near subway tunnels can induce deformation, necessitating a comprehensive investigation into causal factors and mitigation strategies. Field measurements were conducted to assess both vertical and horizontal displacements of existing tunnels near a deep excavation in Shenzhen. Utilizing a validated three-dimensional finite element model that considers structure−strata interactions, this study analyzes tunnel displacements, ground movements, diaphragm wall impacts and the sensitivity of enclosure structure parameters. The results indicate that tunnel deformation correlates with enclosure structure deformation, particularly near the center of the pit. Moreover, shallow soil excavation significantly affects the vertical displacement of shallow-buried tunnels. However, the design parameters of the existing enclosure structures inadequately limit tunnel displacement. Therefore, it is crucial to intensify vertical displacement monitoring in shallow tunnels during early excavation stages and to enhance horizontal displacement monitoring during later phases. Implementing measures such as optimizing central support design or retaining soil at the pit bottom helps control maximum horizontal displacement. While support stiffness plays a greater role than retaining wall thickness, its impact on deep excavation projects is limited.

A novel semi-analytical model for laterally loaded piles

A new model for laterally loaded piles is proposed, consisting of a linear part and a non-linear iterative part. The linear part is based on the Green's function method, which strictly satisfies the displacement and stress continuity conditions at the pile-soil interface and leads to a rigorous solution. The non-linear iterative part assumes that the shear modulus of the soil around the pile decreases with the average strain of the soil layer. Only two soil parameters need to be calibrated, including the initial shear modulus of the soil and the reference shear strain at half the initial modulus. Compared to the finite element model, the advantage of this model is that it converges quickly while maintaining accuracy, which greatly improves computational efficiency. The advantage of this model over the p-y curve model is that it does not have the difficulty of selecting p-y spring parameters, especially for monopiles. The accuracy of this model is verified by comparison with three-dimensional finite element models and nine existing studies, including numerical models, field tests and centrifuge tests. These cases cover various piles with diameters of D∼0.273–10 m and aspect ratios of L/D∼3–40, installed in sand or clay by driving, jacking or vibration.

Effect of Process Parameters on Welding Residual Stress of 316L Stainless Steel Pipe.

316L stainless steel pipes are widely used in the storage and transportation of low-temperature media due to their excellent low-temperature mechanical properties and corrosion resistance. However, due to their low thermal conductivity and large coefficient of linear expansion, they often lead to significant welding residual tensile stress and thermal cracks in the weld seam. This also poses many challenges for their secure and reliable applications. In order to effectively control the crack defects caused by stress concentration near the heat-affected zone of the weld, this paper establishes a thermal elastoplastic three-dimensional finite element (FE) model, constructs a welding heat source, and simulates and studies the influence of process parameters on the residual stress around the pipeline circumference and axial direction in the heat-affected zone. Comparison and verification were conducted using simulation and experimental methods, respectively, proving the rationality of the finite element model establishment. The axial and circumferential residual stress distribution obtained by the simulation method did not have an average deviation of more than 30 MPa from the numerical values obtained by the experimental method. This study also considers the effects of welding energy, welding speed, and welding start position on the pipe's circumferential and axial residual stress laws. The results indicate that changes in welding energy and welding speed have almost no effect on the longitudinal residual stress but have a more significant effect on the transverse residual stress. The maximum transverse residual stress is reached at a welding energy of 1007.4~859.3 J/mm and a welding speed of 6.6 mm/s. Various interlayer arc-striking deflection angles can impact the cyclic phase angle of the transverse residual stress distribution in the seam center, but they do not alter its cyclic pattern. They do influence the amplitude and distribution of the longitudinal residual stress along the circumference. The residual stress distribution on the surface of the pipe fitting is homogenized and improved at 120°.

Three-dimensional finite element analysis of the biomechanical properties of different material implants for replacing missing teeth.

The purpose of this study was to analyze the biomechanical properties of implants made of different materials to replace missing teeth by using three-dimensional finite element analysis and provide a theoretic basis for clinical application. CBCT data was imported into the Mimics and 3-Matic to construct the three-dimensional finite element model of a missing tooth restored by an implant. Then, the model was imported into the Marc Mentat. Based on the variations of the implant materials(titanium, titanium-zirconia, zirconia and poly (ether-ether-ketone) (PEEK))and bone densities (high and low),a total of eight models were created. An axial load of 150N was applied to the crown of the implant to simulate the actual occlusal situation. Both the maximum values of stresses in the cortical bone and implant were observed in the Zr-low model. The maximum displacements of the implants were also within the normal range except for the PEEK models. The cancellous bone strains were mainly distributed in the apical area of the implant, and the maximum value (3225μstrain) was found in PEEK-low model. Under the premise of the same implant material, the relevant data from various indices in low-density bone models were larger than that in high-density bone models. From the biomechanical point of view, zirconia, titanium and titanium-zirconia were all acceptable implant materials for replacing missing teeth and possessed excellent mechanical properties, while the application of PEEK material needs to be further optimized and modified.

Estimation of load-carrying capacity of cracked RC beams using 3D digital twin model integrated with point clouds and images

Digital twin (DT) is capable of managing, simulating, and future performance prediction of existing infrastructures by integrating three-dimensional (3D) geometric model and finite element (FE) model. This study aims to automatically generate geometric DT of existing RC members using point clouds and incorporate damages detected from images into FE models. An edge-based modeling method was proposed to extract sharp features for geometry reconstruction with incomplete point clouds, and an automated procedure was established for registering images to point clouds, mapping cracks detected in images to FE models, and assigning reduced material properties. Four RC beams were tested in a laboratory to validate the proposed method. Existing cracks were found to have a significant effect on structural capacity and failure mode. The results also showed that the developed FE model produced a good agreement between experimental and predicted failure loads while providing a more efficient and time-saving assessment protocol.

Characterization of mechanical stress in the occurrence of cortical opacification in age-related cataracts using three-dimensional finite element model of the human lens and RNA-seq

Cataract is the leading cause of blindness across the world. Age-related cataract (ARC) is the most common type of cataract, but its pathogenesis is not fully understood. Using three-dimensional finite element modeling combining experimental biotechnology, our study demonstrates that external forces during accommodation cause mechanical stress predominantly in lens cortex, basically matching the localization of opacities in cortical ARCs. We identified the cellular senescence and upregulation of PIEZO1 mRNA in HLECs under mechanical stretch. This mechano-induced senescence in HLECs might be mediated by PIEZO1-related pathways, portraying a potential biomechanical cause of cortical ARCs. Our study updates the fundamental insight towards cataractogenesis, paving the way for further exploration of ARCs pathogenesis and nonsurgical treatment.

Stress, damage, and fatigue performance analysis of CFRP/Al double-sided countersunk riveted joints with variable rivet-hole clearance

The CFRP/Al countersunk rivet joints are one of the most fatigued dangerous points of the aircraft fuselage. To deeply investigate the failure mechanisms of the riveted joint structures, this paper utilizes electromagnetic riveting technology to prepare single-lap joint specimens and investigates the quasi-static tensile and fatigue failure behaviours of the joints under different rivet-hole clearances and stress levels. In addition, to assist this investigation, a three-dimensional finite element model, combining riveting and fatigue loading, is developed. This model takes into account the strain rate effects during riveting process and progressive failure behavior of various typical damage modes in CFRP laminates. The results indicate that Aluminum alloy sheets fracture dominates the fatigue failure of the joints. Increasing the clearance can reduce the stress amplitude on the shear plane and further enhance the fatigue life of the joints. Furthermore, all additional CFRP damage during the loading process is the growth of the initial riveting damage, and the damage-induced fibre bearing capacity decrease is considered the reason for the increase in stress amplitude of the Al sheet. Last, excessive clearance enlargement does not appreciably diminish the initial riveting damage, instead, it could result in a joint strength reduction.

Impact of small-angle shield undercrossing on existing tunnels

A three-dimensional finite element model was developed to investigate the effects of small-angle tunnels on existing tunnels in terms of surface settlement, arch bottom settlement, and ellipticalization of existing tunnel linings. The model simulated two small-angle tunnels crossing two existing tunnels. The maximum surface subsidence was observed at the recently excavated tunnel centers. The presence of existing tunnels reduced surface subsidence caused by ring-by-ring excavation. Surface settlement followed a “W” curve pattern, with two local minima gradually getting close as the tunnel was excavated. The maximum value of the arch bottom settlement of the existing tunnel gradually approached the intersection of the existing tunnel and the newly excavated tunnel. When the absolute value of the rate of change of the arch bottom settlement was high, the tunnel produced a large inter-ring misalignment. The line connecting the centers of the existing tunnel and recently excavated tunnel was parallel to the major axis direction of the ring after the deformation of the lining on the existing tunnel. Therefore, a small-angle tunnel passing through an existing tunnel may result in misalignment and torsion of the existing tunnel, requiring careful monitoring of the inter-ring bolt force to ensure safety.

Modeling of laser-assisted micro-milling Inconel718

Nickel-based superalloy Inconel718 is a typical difficult-to-machine material with the characteristic of high strength and hardness. Demands for micro parts of nickel-based superalloy Inconel718 are growing in aerospace, biomedical and other fields. Laser-assisted micro-milling (LAMM) can realize precision machining of the workpiece because of its lower cutting force. However, the machining mechanism and the change law of cutting force during LAMM of nickel-based superalloy Inconel718 are still unclear. There is a crucial need for modeling of LAMM Inconel718. This paper explores LAMM of nickel-based superalloy Inconel718 and realizes the simulation of LAMM process. A three-dimensional transient heat transfer finite element (FE) model is established to predict the temperature distribution of a Gaussian heat source in LAMM. Multiple sets of simulations of laser pre-heating with various laser power are conducted to analyze the temperature fields of workpiece and determine the appropriate laser parameters. A novel 3D model of the full machining process of LAMM Inconel 718 is constructed by ABAQUS with multi-subroutine, a fully coupled thermal-mechanical process is achieved. A modified Johnson-Cook constitutive model is used to reflect the size effect and flow stress of LAMM Inconel 718. Finally, the validity of the models is verified by experiments. Based on the verified model, the appropriate laser pre-heating parameters are obtained, and the prediction of cutting forces of LAMM Inconel718 is achieved.

A numerical investigation for the development of functionally graded Ti/HA tibial implant for total ankle replacement: Influence of material gradation law and volume fraction index.

Stress shielding is one of the major concerns for total ankle replacement implants nowadays, because it is responsible for implant-induced bone resorption. The bone resorption contributes to the aseptic loosening and failure of ankle implants in later stages. To reduce the stress shielding, improvements can be made in the implant material by decreasing the elastic mismatch between the implant and the tibia bone. This study proposes a new functionally graded material (FGM) based tibial implant for minimizing the problem of stress shielding. Three-dimensional finite element (FE) models of the intact tibia and the implanted tibiae were created to study the influence of material gradation law and volume fraction index on stress shielding and implant-bone micromotion. Different implant materials were considered that is, cobalt-chromium, titanium (Ti), and FGM with Ti at the bottom and hydroxyapatite (HA) at the top. The FE models of FGM implants were generated by using different volume fractions and the rule of mixtures. The rule of mixtures was used to calculate the FGM properties based on the local volume fraction. The volume fraction was defined by using exponential, power, and sigmoid laws. For the power and sigmoid law varying volume fraction indices (0.1, 0.2, 0.5, 1, 2, and 5) were considered. The geometry resembling STAR® ankle system tibial implant was considered for the present study. The results indicate that FGMs lower stress shielding but also marginally increase implant-bone micromotion; however, the values were within the acceptable limit for bone ingrowth. It is observed that the material gradation law and volume fraction index influence the performance of FGM tibial implants. The tibial implant composed of FGM using power law with a volume fraction index of 0.1 was the preferred option because it showed the least stress shielding.

Investigation into the effect of deltoid ligament injury on rotational ankle instability using a three-dimensional ankle finite element model.

Injury to the lateral collateral ligament of the ankle may cause ankle instability and, when combined with deltoid ligament (DL) injury, may lead to a more complex situation known as rotational ankle instability (RAI). It is unclear how DL rupture interferes with the mechanical function of an ankle joint with RAI. To study the influence of DL injury on the biomechanical function of the ankle joint. A comprehensive finite element model of an ankle joint, incorporating detailed ligaments, was developed from MRI scans of an adult female. A range of ligament injury scenarios were simulated in the ankle joint model, which was then subjected to a static standing load of 300N and a 1.5Nm internal and external rotation torque. The analysis focused on comparing the distribution and peak values of von Mises stress in the articular cartilages of both the tibia and talus and measuring the talus rotation angle and contact area of the talocrural joint. The dimensions and location of insertion points of ligaments in the finite element ankle model were adopted from previous anatomical research and dissection studies. The anterior drawer distance in the finite element model was within 6.5% of the anatomical range, and the talus tilt angle was within 3% of anatomical results. During static standing, a combined rupture of the anterior talofibular ligament (ATFL) and anterior tibiotalar ligament (ATTL) generates new stress concentrations on the talus cartilage, which markedly increases the joint contact area and stress on the cartilage. During static standing with external rotation, the anterior talofibular ligament and anterior tibiotalar ligament ruptured the ankle's rotational angle by 21.8% compared to an intact joint. In contrast, static standing with internal rotation led to a similar increase in stress and a nearly 2.5 times increase in the talus rotational angle. Injury to the DL altered the stress distribution in the tibiotalar joint and increased the talus rotation angle when subjected to a rotational torque, which may increase the risk of RAI. When treating RAI, it is essential to address not only multi-band DL injuries but also single-band deep DL injuries, especially those affecting the ATTL.

Deformation caused by T-shape expansion of existing metro station: complex challenges and strategies

This case study presents an investigation on existing metro station deformation characteristics during T-type expansion in Nanjing, China. A three-dimensional elastic-plastic finite element model is adopted to consider the interaction of soil, retaining structure, horizontal support system and existing station structure. The results indicate that coupling deformation of asymmetric uplift deflection and existing structure rotation occurs on existing structure. Moreover, the sensitivity of influencing factors is analyzed, including the dimension of retaining structure, the support spacing, the excavation method, the soil reinforcement depth in passive area, the width of retaining soil during excavation, the relative position between foundation pit and existing structure, the excavation depth, the existing structure rigidity. Furthermore, based on the analysis of difficulties in the T-shape expansion projects, the controlling measures for the deformation of existing metro station are proposed. The achievements of this case study can be used as a reference for similar underground expansion projects.

A nonlinear analytical method for pipeline strain under soil thaw subsidence load in permafrost region

In frozen terrain, the transportation of heated oil triggers thermal soil thaw and collapse, leading to pipeline distortion and breakdown. During substantial ground settlement in such terrains, current strain prediction techniques for pipelines fall short, neglecting significant pipeline bending, non-linear material behaviors, or detailed pipeline-permafrost interplay. To address this, a model focusing on pipeline-soil dynamics is devised, incorporating material non-linearity and pronounced pipeline distortion. This model derives analytical strain solutions for pipelines by aggregating incremental distortions of pipeline components via a strain iteration technique. Concurrently, a three-dimensional solid finite element model is crafted, assessing pipeline strains across varying subsidence zone dimensions, factoring in the complicated permafrost environment. The model's efficiency is affirmed by its less than 3% discrepancy from finite element method results, showcasing its enhanced accuracy and dependability in predicting pipeline strains under extensive soil settlement scenarios.