3D Finite Element Model Research Articles

Experimental and numerical investigations on the impact of openings on the stiffness of CLT shear walls

Cross-laminated timber (CLT) exhibits high in-plane strength and stiffness, positioning it as a viable material for shear walls in seismic regions. However, the influence of openings on the in-plane stiffness of CLT shear walls remains inadequately studied. In this study, 43 CLT shear walls with different opening sizes were tested under in-plane loading. The stiffness reduction was less than 29% for openings 800×800mm (53% of panel width), but exceeded 60% for openings 1200×1200mm (80% of panel width). When assessed with identical opening sizes, 139mm and 175mm thick CLT shear walls displayed consistent stiffness reductions. Panels with an aspect ratio of 3:1 exhibited a more pronounced drop than 2:1 aspect ratio panels. A validated 3D finite element model was used to investigate the influence of panel layup and opening shape and location on the in-plane stiffness reduction. For a distance between the opening and panel edge exceeding 300mm, the shape and location of openings exerted negligible effects. Finally, the Dujič’s model was shown to be a suitable predictor to evaluate the in-plane stiffness reduction of CLT shear walls with openings.

A Physics Informed Neural Network for Solving the Inverse Heat Transfer Problem in Gas Turbine Rotating Cavities

Abstract Recently, Physics-Informed Neural Networks have displayed great potential in delivering swift and accurate solutions for inverse problems. Here, a Physics-Informed Neural Network (PINN) was used to solve the inverse heat conduction problem in the rotating cavities of a aeroengine high-pressure compressor internal air system. The neural network was designed to receive experimentally captured radially distributed temperature profiles as inputs and predict the associated surface heat fluxes. The correctness of these predicted heat fluxes are assessed by numerically solving the direct heat conduction equation using a 2D Finite-Element model, thereby recovering the original temperature profiles. A comparative analysis is conducted between the predicted temperature profiles and the initial inputs. The physics informed neural network is trained using noise free synthetic data, created from a range of radial temperature curve fit coefficients, and subsequently tested on noisy experimental data at engine representative conditions. The predicted temperature values exhibit some good agreement with their respective actual counterparts. Furthermore, the sensitivity of model hyperparameters are explored to showcase the capability of the proposed approach. The results show that Physics- Informed Neural Networks exhibit reduced susceptibility to experimental uncertainties when addressing inverse problems, in contrast to inverse solution methods, and offer a possible new approach for analysis of experimental data if trained on a sufficiently large dataset.

Electromagnetic responses on microstructures of duplex stainless steels based on 3D cellular and electromagnetic sensor finite element models

Electromagnetic responses on microstructures of duplex stainless steels based on 3D cellular and electromagnetic sensor finite element models

Noise generation and contribution of a single-gap bridge expansion joint: An experimental and numerical study

Noise generated from bridge expansion joints during vehicle passing has caused localized environmental problems around the world. The purpose of this paper is to investigate the vibration and noise radiation from the single-gap expansion joint induced by passing vehicles. First, a vibration and noise field test was performed for single-gap expansion joints in the highway ramps. Based on the tested results, two numerical models, involving a 3D tire-bridge expansion joint hybrid finite element model for the vibration analysis, as well as an acoustic boundary element model for the prediction of the noise radiation, were established and validated. The results of the numerical simulations and the field tests are well matched, which offers opportunities to predict the vibration and noise induced by vehicles passing through using the proposed methodology. Finally, the acoustic contribution and sound pressure distribution of different sound sources were discussed. The research showed that, for the single-gap expansion joint, the contribution of tire noise to the total noise level is dominant, and the contribution of structure-borne noise of the expansion joint and girder to the total noise level is insignificant. The findings emphasize the importance of tire noise and provide new ideas for noise control of single-gap expansion joints.

Numerical analysis of axially loaded concrete-filled steel tubular columns strengthened with CFRP grid-reinforced ECC

Carbon fiber-reinforced polymer (CFRP) and engineered cementitious composites (ECC) are commonly used to reinforce structural elements, with proven effectiveness in enhancing the performance of concrete-filled steel tube (CFST) columns as shown in prior experimental studies. However, the complex interaction among the components of the reinforced columns, including CFRP, ECC, the steel tube, and the concrete core, as well as the underlying mechanical mechanisms, remain insufficiently understood. This research aims to fill these gaps by developing a 3D finite element (FE) model of circular CFST short columns reinforced with CFRP grid-reinforced ECC for comprehensive numerical analysis. Sensitivity analyses are conducted to determine the governing parameters in the numerical simulation, including mesh size and critical plasticity parameters for the concrete damage model, such as the ratio of the second stress invariant on the tensile meridian to that on the compressive meridian (Kc), and the dilation angle (ψ) for both the concrete core and ECC. The model also accounts for the effects of concrete confinement provided by the CFRP grid-reinforced ECC layer and the circular steel tube. The accuracy of the FE model is validated against existing experimental data, and the model is subsequently used to investigate the behavior of reinforced CFST columns under varying geometric and material properties. The study examines failure modes, axial load-displacement responses, stress distributions, and the contribution of each material component to the load-bearing capacity of the strengthened columns. The results indicate an increase in ultimate axial strength ranging from 13% to 20% for the reinforced CFST columns. While the compressive strength of the ECC has a moderate effect on axial resistance, increasing ECC thickness, concrete strength, and steel tube yield strength significantly enhances the columns’ load-bearing capacity. Finally, a novel design model is proposed and validated, which accounts for the confinement effects provided by the CFRP grid and ECC on the concrete core.

Designing Multi-axis Compliant Mechanisms With Lockable Decoupled Inputs: A Tip–Tilt Case Study

Abstract This research is about designing multi-degree-of-freedom (multi-DOF) compliant mechanisms with decoupled inputs that can be independently locked/unlocked using bistable switches to achieve different combinations of DOFs. A case study mechanism achieving two decoupled rotational DOFs (tip and tilt) is designed, fabricated, and characterized. It can be triggered using two pairs of bistable switches, achieving drastically different states of torsional stiffness for each DOF in four sets of DOF combinations—no DOFs, a tip DOF, a tilt DOF, and both tip and tilt DOFs. Bistability and stiffness cancelation principles are exploited to achieve the desired changes in stiffness. Two flexure elements can be identified—the switch providing a negative stiffness and the cross-axis-flexural-pivot (CAFP) producing a positive stiffness. The mechanism is tuned to achieve static balancing, reaching a near-zero stiffness over much of its range. The pseudo-rigid body model and two-dimensional (2D) finite element model (FEM) are combined defining a fast method to dimension the system. The 3D FEM is simulated to validate the obtained results. For each DOF, the system is tested in two configurations (stiff and compliant) for three cycles over a ±10 deg rotation, achieving a stiffness reduction of around 99%. Comparable stiffness values were measured after triggering the switches more than once, repetitively reaching the required two states of stiffness, confirming the system's usability in practical applications. The positive stiffness provided by the CAFP is measured and compared to the device's overall stiffness, highlighting the stiffness cancelation concept.

An analytical model for predicting the magnetization loss in HTS sector-shaped conductors for fusion

Abstract Within the framework of magnetic confinement fusion, several projects worldwide are demonstrating the possibility of integrating high-temperature superconductors (HTS) in the coil systems. HTS-based technologies are highly attractive for practical applications because they can extend the operating margins of fusion coils in terms of higher temperatures, transport currents and magnetic fields. Based on the results achieved with the twisted-stacked tape cable, we have designed a novel low-loss HTS sector cable-in-conduit conductor, with a target of 60 kA at 4.5 k, 18 T, which is presently of interest for the DEMO Central Solenoid coil. In HTS cables, the AC losses can represent a significant limiting factor, therefore they must be taken into consideration both in the design phase and in the assessment of the overall magnet thermal budget. In this work, to assess the loss behavior and to optimize the cable design, we have explored different aspect ratios and arrangements of the stacked tapes within the cable layout. The magnetization losses are calculated with a 2D finite-element model based on the T-A formulation and analytical approximations based on the Brandt-Halse critical state model. Specifically, we have developed an analytical formulation that allows for the calculation of the instantaneous power losses in HTS stacked cables with a limited number of tapes per stack, achieving sufficient accuracy at high fields. The analytical model enables a sufficiently accurate assessment of the heat deposited on the conductor during those particular instants of a plasma scenario where the variation of the field is very high, such as during the critical initial discharge period of the plasma initiation.

Non‐linear time history analyses of a rigid block isolated with unbonded fiber‐reinforced elastomeric isolators (UFREIs): A comparison between 3D finite element and phenomenological models

AbstractNumerical modeling represents a pivotal tool in the seismic analysis and design of structural systems, enabling the detailed prediction and examination of structural responses under seismic loading. This research conducts a comparative analysis of two numerical modeling approaches aimed at simulating the seismic response of unbonded fiber‐reinforced elastomeric isolators (UFREIs). The research focuses on a finite element (FE) model developed using Abaqus and a developed phenomenological model implemented in OpenSees, outlining the development and calibration processes for each. The FE model is developed based on simple rubber material testing data, while the phenomenological model is calibrated using experimental results from cyclic shear tests conducted on the UFREI device and the FE model. The primary objective of this study is to assess the effectiveness of these modeling approaches in predicting UFREI behavior under seismic conditions. This evaluation entails comparing model predictions with experimental data obtained from unidirectional shake table tests performed on a rigid block isolated by two UFREIs. This paper highlights the distinct advantages and limitations of each model in simulating UFREI dynamic responses during seismic events. Furthermore, it provides insights into the modeling techniques and discusses the computational demands and data requirements of each model, thereby aiding in their application to various aspects of seismic analysis and design.

Veneer-wrapped steel columns: Proof-of-concept experimental and numerical investigations

This paper presents and discusses the structural behaviour of a novel type of steel–timber composite column to be used in residential building applications. It consists of a circular hollow section wrapped by multiple beech veneer sheets and where the composite action is ensured by a structural bi-component glue applied over the contact surfaces. As part of a proof-of-concept study, five short and three long columns with lengths of 0.27 m and 2.00 m, respectively, were tested under vertical concentric load with fibres parallel to the load direction and inclined by 15°. The global response of the specimens in terms of load-displacement behaviour and failure modes was studied, along with localised strain measurements through strain gauges and a digital image correlation (DIC) system. Results showed a significant increase in buckling resistance due to the stiffness added by the veneer layers. Already with the still-moderate veneer layer thicknesses used in this study, an increase in strength of over 25 % could thus be achieved. This enhancement highlights the effectiveness of veneer wrapping in stabilizing the inner core and delaying the occurrence of flexural buckling failure. A 3D non-linear finite element model (FEM) validated these findings, showing a good correlation with experimental results in resistance and stiffness. A parametric study was conducted to explore different geometries and materials, thus creating an extended FEM-based database. A relatively good agreement with modest deviation between the FEM predictions and analytical formulations available in the European Standards was found, considering the appropriate introduction of the mechanical properties of timber.

Enhancing formability in deep drawing: A 3D finite element analysis of plastic wrinkling in AA2090 Al–Li alloy

Deep drawing is a manufacturing process used to produce billions of metal containers, crucial for the aerospace industry in forming thin-walled components efficiently. However, wrinkling due to compressive instability remains a significant challenge in sheet metal forming. This study investigates plastic wrinkling in sheet metal forming to enhance the final drawn shape. A 3D finite element model using ABAQUS/Explicit was developed to examine the effects of mesh topology, blank holder force, and sheet thickness on wrinkling in AA2090 Al–Li alloy. The finding indicates that that an increase in blank holder force increases wrinkle number and amplitude. Additionally, it was found that a 2000 N blank holder force is necessary to prevent wrinkling.

Coherent to semi-coherent transition of precipitates in Van der Merwe epitaxial films

Abstract The interplay of stresses of two coherent systems, specifically a coherent precipitate within an epitaxial film (in Van der Merve growth mode), gives rise to an interesting scenario; wherein the configurational effects and the strain energy landscape are enriched. This leads to an interdependent evolutionary pathway for the coherent to semi-coherent transition. The transition of the film-substrate interface to the semi-coherent state occurs beyond a critical thickness ( t c ) which arises from energy minimization. Using the model example of the precipitation of NbH0.5 precipitate in a Nb/Sapphire epitaxial system, the critical radii ( r film ∗ ) for model geometries have been computed. 3D finite element models have been used to compute the stress state and energetics of precipitates growing in an epitaxial film. The following important observations can be made based on the simulations. (1) The epitaxial stresses can lead to the stabilization of the coherent precipitate. (2) The transition to a semi-coherent film interface, characterized by an array of interfacial misfit dislocations, can energetically favor the semi-coherent state of the precipitate over the coherent state. The magnitude of r film ∗ depends on the spatial coordinates within the film. (3) The presence of Nb–H coherent precipitates can stabilize the coherent state of the film interface.

Mitigation and Optimization of Induced Seismicity Using Physics‐Based Forecasting

AbstractFluid injection can induce seismicity by altering stresses on pre‐existing faults. Here, we investigate minimizing induced earthquake potential by optimizing injection operations in a physics‐based forecasting framework. We built a 3D finite element model of the poroelastic crust for the Raton Basin, Central US, and used it to estimate time dependent Coulomb stress changes due to 25 years of wastewater injection in the region. Our finite element model is complemented by a statistical analysis of the seismogenic index (SI), a proxy for critically stressed faults affected by variations in the pore pressure. Forecasts of seismicity rate from our hybrid physics‐based statistical model suggest that induced seismicity in the Raton Basin, from 2001 to 2022, is still driven by wastewater injection despite declining injection rates since 2011. Our model suggests that pore pressure diffusion is the dominant cause of Coulomb stress changes at seismogenic depth, with poroelastic stress changes contributing about 5% to the driving force. Linear programming optimization for the Raton Basin reveals that it is feasible to reduce earthquake potential for a given amount of injected fluid (safety objective) or maximize fluid injection for a prescribed earthquake potential (economic objective). The optimization tends to spread out high‐rate injectors and shift them to regions of lower SI. The framework has practical importance as a tool to manage injection rate per unit field area to reduce induced earthquake potential. Our optimization framework is both flexible and adaptable to mitigate induced earthquake potential in other regions and for other types of subsurface fluid injection.

Biomechanical effects of functional clear aligners on the stomatognathic system in teens with class II malocclusion: a new model through finite element analysis

ObjectivesThe Functional Clear Aligner (FCA) is a novel orthodontic appliance designed for the treatment of Class II malocclusion with mandibular retrognathia in adolescents. The aim of this study was to investigate the biomechanical characteristics of the masticatory muscles, jawbone, and temporomandibular joint (TMJ) during mandibular advancement using either FCA or Class II elastics combined with clear aligner (Class II elastics) through finite element analysis.Materials and methodsA 3D finite element model of the ‘muscle-jawbone-TMJ-appliance’ system was constructed based on CBCT and MRI images of a boy with skeletal Class II malocclusion. Masticatory muscles included masseter, temporal, medial pterygoid, and lateral pterygoid muscles. The TMJ consists of the temporal bone’s glenoid fossa, disc, and mandibular condyle. To observe the biomechanical characteristics of the muscles and TMJ during orthodontic appliance wearing and the retention phase, two different protocols were used: Model 1: The mandibular advancement using FCA; Model 2: The mandibular advancement using Class II elastics.ResultsThe FCA group produced greater and more coordinated masticatory muscle forces compared to the Class II elastics group. Temporal and masseter muscles exhibited the most pronounced variation in muscle strength during mandibular advancement. The FCA group exhibited greater TMJ region stress compared to the Class II elastics group. Interestingly, the stress on the articular discs in both models decreased over time. Tensile stresses were observed in both the condyle and the posterior region of the articular fossa.ConclusionDuring skeletal Class II malocclusion treatment, masticatory muscle forces and stress on the TMJ were higher in the FCA group compared to the Class II elastics group. In both models, stress cushioning was provided by the articular disc.

Numerical investigation of the effect of laser peen forming process parameters on 1060 aluminium sheets under preload

Abstract Laser peen forming process changes the shape of the metal sheet or plate by plastic deformation caused by laser induced mechanical shock wave. The bending behavior of 1060 aluminium sheets under laser peen forming (LPF) are studied numerically taking into account the various laser peening process parameters for two distinct cases; with and without preload. Modal analysis has been performed on 3D finite element model (FEM) to stop post-shock oscillations. A time-dependent and space-dependent pressure distribution is applied to the laser-peened area along with preload at the free end. The numerical comparison of bending angles and surface profile of 1060 aluminium sheets under laser peen forming are studied for both the cases, with and without preload. Preloading provides better surface profile for same bending angle for a no laser shock overlap.

Parametric investigation of friction stir welding of aluminum alloy and Inconel 718 using finite element analysis

Friction Stir Welding (FSW) has emerged as a prominent technique for joining metallic materials, offering advantages in both similar and dissimilar material combinations. Despite its effectiveness, further advancements are needed to address industry demands and enhance the understanding of the FSW process. Welding nickel-based superalloys like Inconel 718 poses specific challenges, with existing literature suggesting that achieving successful welds requires high axial forces and precise process parameters, limiting its practical application in industries. To tackle these issues, this research proposes the utilization of 3D finite element models integrating multiphysics aspects, encompassing material flow and heat transfer mechanisms such as conduction, convection, and radiation. These models were validated against existing experimental data and employed to analyze temperature and strain distributions within the heat affected zone and weld nugget. The findings offer insights into the impact of various process parameters, including rotational speed, welding speed, normal force, cooling rate, and the effect of induction preheating, on key performance metrics like temperature profiles, grain size distribution, microhardness, and stress evolution. The outcomes underscore a notable correlation between process variables and performance indicators, facilitating a comprehensive understanding of the FSW process dynamics.

Strain and Deformation Analysis Using 3D Geological Finite Element Modeling with Comparison to Extensometer and Tiltmeter Observations

This study verifies the practicality of using finite element analysis for strain and deformation analysis in regions with sparse GNSS stations. A digital 3D terrain model is constructed using DEM data, and regional rock mass properties are integrated to simulate geological structures, resulting in the development of a 3D geological finite element model (FEM) using the ANSYS Workbench module. Gravity load and thermal constraints are applied to derive directional strain and deformation solutions, and the model results are compared to actual strain and tilt measurements from the Jiujiang Seismic Station (JSS). The results show that temperature variations significantly affect strain and deformation, particularly due to the elevation difference between the mountain base and summit. Higher temperatures increase thermal strain, causing tensile effects, while lower temperatures reduce thermal strain, leading to compressive effects. Strain and deformation patterns are strongly influenced by geological structures, gravity, and topography, with valleys experiencing tensile strain and ridges undergoing compression. The deformation trend indicates a southwestward movement across the study area. A comparison of FEM results with ten years of strain and tiltmeter data from JSS reveals a strong correlation between the model predictions and actual measurements, with correlation coefficients of 0.6 and 0.75 for strain in the NS and EW directions, and 0.8 and 0.9 for deformation in the NS and EW directions, respectively. These findings confirm that the 3D geological FEM is applicable for regional strain and deformation analysis, providing a feasible alternative in areas with limited GNSS monitoring. This method provides valuable insights into crustal deformation in regions with sparse strain and deformation measurement data.

Flexural behavior of bonded post-tensioned concrete beams with steel or GFRP bars

Post-tensioned concrete beams with GFRP reinforcement have gained attention in recent years due to their unique flexural behavior and potential cost benefits. Because of the absence of comprehensive guidelines in international codes and limited application in constructions, this hybrid system is adopted in the current paper as a pioneering solution for sustainable and eco-friendly construction practices. Particularly noteworthy is the corrosion resistance of GFRP bars, which significantly enhances the durability of structures. The influence of using GFRP with bonded post-tension concrete beam is experimentally investigated compared to concrete beams with prestressing and non-prestressing reinforcement. Eight concrete beams with different reinforcement types and ratios were cast and tested under three-point load. Flexural resistance, crack pattern, deflections, and strains of the beams were compared with a reference beam with prestressing and non-prestressing reinforcement. Experimental load capacities of the tested beams were compared to those obtained from various international codes. The experimental results revealed that the beams with GFRP bars had smaller failure load and larger deflection compared to those with non-prestressing steel by up to 25 % and 37 %, respectively. A 3D finite element model (FEM) was developed and validated by comparing the numerical results with those obtained from the experiments. A sensitivity analysis was performed, using the validated FEM, by changing the concrete strength, steel yield stress, and GFRP ultimate strength by ± 20 %. The sensitivity analysis included seventy-two beams which had the same dimensions, loading, and boundary conditions of the tested beams. Changing the concrete strength only led to varying the beam capacities by −7 % ∼ + 4 %. Changing the steel yield stress only led to changing of the beam capacities by −13 % ∼ + 11 %. No change was observed when the GFRP ultimate stress was varied by 20 %.

Detection of delamination using laser-generated Lamb waves with improved wavenumber analysis

Abstract Delamination is a typical failure mode in multilayered materials, which may potentially threaten the safety and reliability of the materials if not detected promptly. This paper proposes a novel quantitative detection method for delamination in multilayered materials based on laser-generated Lamb waves with improved wavenumber analysis. First, a three-dimensional(3D) finite element model is established to simulate the interaction between delamination and laser-generated Lamb waves. The propagation characteristics of Lamb waves in multilayered materials without defects and with delamination of three different shapes, i.e., rectangular, circular and triangular are analyzed. And the mechanism that delamination can lead to the appearance of scattered wave and new wavenumber components is investigated, facilitating the quantitative detection of delamination. Then, an improved wavenumber analysis method for quantitative detection of delamination is proposed. The new wavenumber components caused by delamination are extracted by broadband adaptive wavenumber filtering to reduce imaging interference and artifacts caused by irrelevant components, and the spatial wavenumber spectrum is obtained by fast broadband local wavenumber estimation algorithm to enhance the delamination imaging accuracy and spatial resolution. The experimental results indicate that, across various experimental scenarios, the proposed method achieves superior imaging precision compared to traditional wavenumber filtering or local wavenumber estimation methods, and can quantitatively identify delamination defects of various shapes and sizes.

Identification of a double decay length ([formula omitted]) heat flux deposition shape with embedded thermal measurement and neural network

The estimation of the heat flux density distribution profiles in tokamak devices is a very important research topic for edge plasma physics purposes and also to ensure the safety of the machine. In the radial direction, the heat flux exhibits an exponential decay that could be captured by thermal sensors distributed in the plasma facing components. Radially distributed thermal sensors based on Fiber Bragg grating technology have been embedded in the WEST lower divertor to study the heat flux deposition profiles during plasma operation. The comparison between embedded measurements and a 3D finite element model shows a small decay length (5 – 10 mm) on top of a wider heat flux with a decay length around 30 to 50 mm. A tool using neural network has been developed in order to predict the values of the different parameters describing the deposited heat flux from embedded temperature measurements in steady state. A large span of deposited heat fluxes with maximum heat flux ranging from 1 to 9 MW/m2 and decay length from 5 to 50 mm were characterized using this tool over a database of more than 250 experimental L-mode pulses performed in WEST in attached divertor configuration. The comparison of the predicted heat flux parameters values with macroscopic plasma parameters have revealed the appearance of the narrow component with the increase of the divertor power load (Pdiv) with a threshold dependant of the plasma current (IP).

Three-dimensional finite element analysis of biological signal feedback in the mechanical properties of sound production

In response to the problems of low reliability and long time consumption in traditional research on sound production, this article used electromyography (EMG) signals as input signals to build a high-precision sound production mechanics model. A Proportional-Integral-Derivative (PID) controller was used to dynamically adjust the model and develop a real-time feedback system. This article established a detailed three-dimensional (3D) finite element model including the vocal cords and throat, defined the nonlinear elastic properties and linear elastic properties of different tissues, and used tetrahedral grid partitioning technology to improve the computational accuracy of the model. Through a EMG sensor, an individual EMG was collected and filtered to remove noise. The processed EMGs were used as input parameters for the finite element model to drive the muscle units in the model. By using a PID controller to receive real-time EMG input, the error was calculated and the model was adjusted, enabling accurate simulation of the mechanical properties of vocal cord vibration under different vocal states and achieving real-time feedback. Considering the complexity of vocal cord vibration driven by biological signals, this article simulated and analyzed the modal characteristics of vocal cord vibration, and analyzed the differences in vocal cord vibration characteristics under different vocal states. The sound pressure distribution and resonance frequency were simulated and analyzed to understand the propagation characteristics of sound. Finally, by comparing and analyzing the simulated data with the actual collected data, it was found that the maximum relative error rate of the model under different sound states was 6.14%, and the overall error rate of the model was relatively low, which verified the reliability of the model. The findings demonstrated that the feedback delays of the model in different sound states were all within 100 milliseconds, indicating that the system had high real-time performance and accuracy, which was promising in practical applications.