Search results

This paper aims to investigate the performance of two novel numerical methods, the face-based smoothed finite element method (FS-FEM) and the edge-based smoothed finite element method (ES-FEM), which employ linear tetrahedral elements, for the purpose of strength assessment of a high-speed train hollow axle.

The calculation of stress for the wheelset, comprising an axle and two wheels, is facilitated through the application of the European axle strength design standard. This standard assists in the implementation of loading and boundary conditions and is exemplified by the typical CRH2 high-speed train wheelset. To evaluate the performance of these two methods, a hollow cylinder cantilever beam is first used as a benchmark to compare the present methods with other existing methods. Then, the strength analysis of a real wheelset model with a hollow axle is performed using different numerical methods.

The results of deflection and stress show that FS-FEM and ES-FEM offer higher accuracy and better convergence than FEM using linear tetrahedral elements. ES-FEM exhibits a superior performance to that of FS-FEM using linear tetrahedral elements, showing accuracy and convergence close to FEM using hexahedral elements.

This study channels the novel methods (FS-FEM and ES-FEM) in the static stress analysis of a railway wheelset. Based on the careful testing of FS-FEM and ES-FEM, both methods hold promise as more efficient tools for the strength analysis of complex railway structures.

In this paper, the smoothed point interpolation method (SPIM) is used to model the slope deformation. However, the computational efficiency of SPIM is not satisfying when modeling the large-scale nonlinear deformation problems of geological bodies.

In this paper, the SPIM is used to model the slope deformation. However, the computational efficiency of SPIM is not satisfying when modeling the large-scale nonlinear deformation problems of geological bodies.

A simple slope model with different mesh sizes is used to verify the performance of the efficient face-based SPIM. The first accelerating strategy greatly enhances the computational efficiency of solving the large-scale slope deformation. The second accelerating strategy effectively improves the convergence of nonlinear behavior that occurred in the slope deformation.

The designed efficient face-based SPIM can enhance the computational efficiency when analyzing large-scale nonlinear slope deformation problems, which can help to predict and prevent potential geological hazards.

This article aims to present a direct solution to handle linear constraints in finite element (FE) analysis without penalties or the Lagrange multipliers introduced.

First, the system of linear equations corresponding to the linear constraints is solved for the leading variables in terms of the free variables and the constants. Then, the reduced system of equilibrium equations with respect to the free variables is derived from the finite-dimensional virtual work equation. Finally, the algorithm is designed.

The proposed procedure is promising in three typical cases: (1) to enforce displacement constraints in any direction; (2) to implement local refinements by allowing hanging nodes from element subdivision and (3) to treat non-matching grids of distinct parts of the problem domain. The procedure is general and suitable for 3D non-linear analyses.

The algorithm is fitted only to the Galerkin-based numerical methods.

The proposed procedure does not need Lagrange multipliers or penalties. The tangential stiffness matrix of the reduced system of equilibrium equations reserves positive definiteness and symmetry. Besides, many contemporary Galerkin-based numerical methods need to tackle the enforcement of the essential conditions, whose weak forms reduce to linear constraints. As a result, the proposed procedure is quite promising.

Fused Filament Fabrication (FFF) is an extrusion-based manufacturing process using fused thermoplastics. Despite its low cost, the FFF is not extensively used in high-value industrial sectors mainly due to parts' anisotropy (related to the deposition strategy) and residual stresses (caused by successive heating cycles). Thus, this study aims to investigate the process improvement and the optimization of the printed parts.

In this work, a meshless technique – the Radial Point Interpolation Method (RPIM) – is used to numerically simulate the viscoplastic extrusion process – the initial phase of the FFF. Unlike the FEM, in meshless methods, there is no pre-established relationship between the nodes so the nodal mesh will not face mesh distortions and the discretization can easily be modified by adding or removing nodes from the initial nodal mesh. The accuracy of the obtained results highlights the importance of using meshless techniques in this field.

Meshless methods show particular relevance in this topic since the nodes can be distributed to match the layer-by-layer growing condition of the printing process.

Using the flow formulation combined with the heat transfer formulation presented here for the first time within an in-house RPIM code, an algorithm is proposed, implemented and validated for benchmark examples.

This study aims to propose a numerical modeling procedure for response analysis of elastic body floating in water and submitted to regular waves. An equivalent simplified mechanical single-degree-of-freedom system allowing to reproduce the heave movements is first developed, then the obtained lumped characteristics are used for elastic analysis of the floating body in heave motion.

First, a two-dimensional numerical model of a rigid floating body in a wave tank is implemented under DualSPHysics, an open source computational fluid dynamics (CFD) code based on smoothed particle hydrodynamics method. Then, the obtained results are exploited to derive an equivalent mechanical mass-spring-damper model. Finally, estimated equivalent characteristics are used in a structural finite element modeling of the considered body assuming elastic behavior.

Obtained results concerning the floating body displacements are represented and validated using existing experimental data in the literature. Wave forces acting on the body are also evaluated. It was found that for regular waves, it is possible to replace the complex CFD refined model by an equivalent simplified mechanical system which makes easy the use of structural finite element analysis.

The originality of this work lies in the proposed procedure to evaluate the mechanical properties of the equivalent elastic system. This allows to couple two different software tools and to take advantages of their features.

In structural health monitoring, localization of multiple slight damage without baseline data is significant and difficult. The purpose of this paper is to discuss these issues.

Damage in the structure causes singularities of displacement modes, which in turn reveals damage. Methods based on the displacement modes may fail to accurately locate the slight damage because the slight damage in engineering structure results in a relatively small variation of the displacement modes. In comparison with the displacement modes, the strain modes are more sensitive to the slight damage because the strain is the derivative of the displacement. As a result, the slight variation in displacement data will be magnified by the derivative, leading to a significant variation of the strain modes. A novel method based on strain modes is proposed for the purpose of accurately locating the multiple slight damage.

In the two bay beam and steel fixed-fixed beams, the numerical simulations and the experimental cases, respectively, illustrate that the proposed method can achieve more accurate localization in comparison with the one based on the displacement modes.

The paper offers a practical approach for more accurate localization of multiple slight damage without baseline data. And the robustness to measurement noise of the proposed method is evaluated for increasing levels of artificially added white Gaussian noise until its limit is reached, defining its range of practical applicability.

This study implements the fourth-order phase field method (PFM) for modeling fracture in brittle materials. The weak form of the fourth-order PFM requires C¹ basis functions for the crack evolution scalar field in a finite element framework. To address this, non-Sibsonian type shape functions that are nonpolynomial types based on distance measures, are used in the context of natural neighbor shape functions. The capability and efficiency of this method are studied for modeling cracks.

The weak form of the fourth-order PFM is derived from two governing equations for finite element modeling. C⁰ non-Sibsonian shape functions are derived using distance measures on a generalized quad element. Then these shape functions are degree elevated with Bernstein-Bezier (BB) patch to get higher-order continuity (C¹) in the shape function. The quad element is divided into several background triangular elements to apply the Gauss-quadrature rule for numerical integration. Both fourth-order and second-order PFMs are implemented in a finite element framework. The efficiency of the interpolation function is studied in terms of convergence and accuracy for capturing crack topology in the fourth-order PFM.

It is observed that fourth-order PFM has higher accuracy and convergence than second-order PFM using non-Sibsonian type interpolants. The former predicts higher failure loads and failure displacements compared to the second-order model due to the addition of higher-order terms in the energy equation. The fracture pattern is realistic when only the tensile part of the strain energy is taken for fracture evolution. The fracture pattern is also observed in the compressive region when both tensile and compressive energy for crack evolution are taken into account, which is unrealistic. Length scale has a certain specific effect on the failure load of the specimen.

Fourth-order PFM is implemented using C¹ non-Sibsonian type of shape functions. The derivation and implementation are carried out for both the second-order and fourth-order PFM. The length scale effect on both models is shown. The better accuracy and convergence rate of the fourth-order PFM over second-order PFM are studied using the current approach. The critical difference between the isotropic phase field and the hybrid phase field approach is also presented to showcase the importance of strain energy decomposition in PFM.

The present paper aims to explore a computational homogenisation procedure to investigate the full geometric representation of yield surfaces for isotropic porous ductile media. The effects of cell morphology and imposed boundary conditions are assessed. The sensitivity of the yield surfaces to the Lode angle is also investigated in detail.

The microscale of the material is modelled by the concept of Representative Volume Element (RVE) or unit cell, which is numerically simulated through three-dimensional finite element analyses. Numerous loading conditions are considered to create complete yield surfaces encompassing high, intermediate and low triaxialities. The influence of cell morphology on the yield surfaces is assessed considering a spherical cell with spherical void and a cubic RVE with spherical void, both under uniform strain boundary condition. The use of spherical cell is interesting as preferential directions in the effective behaviour are avoided. The periodic boundary condition, which favours strain localization, is imposed on the cubic RVE to compare the results. Small strains are assumed and the cell matrix is considered as a perfect elasto-plastic material following the von Mises yield criterion.

Different morphologies for the cell imply in different yield conditions for the same load situations. The yield surfaces in correspondence to periodic boundary condition show significant differences compared to those obtained by imposing uniform strain boundary condition. The stress Lode angle has a strong influence on the geometry of the yield surfaces considering low and intermediate triaxialities.

The exhaustive computational study of the effects of cell morphologies and imposed boundary conditions fills a gap in the full representation of the flow surfaces. The homogenisation-based strategy allows us to further investigate the influence of the Lode angle on the yield surfaces.

This paper aims to investigate the responses of single piles and pile groups due to tunneling-induced ground movements in a two-layered soil system. The analyses mainly focus on the additional single pile responses in terms of bending moment, lateral deflection, axial force, shaft resistance and pile settlement. Subsequently, a series of parametric studies were carried out to better understand the responses of single piles induced by tunneling. To give further understanding regarding the pile groups, a 2 × 2 pile group with two different pile head conditions, namely, free and capped, was considered.

Using the PLAXIS three-dimensional (3D) software, a full 3D numerical modeling is performed to investigate the effects of ground movements caused by tunneling on adjacent pile foundations. The numerical model was validated using centrifuge test data found in the literature. The relevance of the 3D model is also judged by comparison with the 2D plane strain model using the PLAXIS 2D code.

The numerical test results reveal that tunneling induces significant displacements and internal forces in nearby piles. The magnitude and distribution of internal forces depend mainly on the position of the pile toe relative to the tunnel depth and the distance between the pile and the vertical axis of the tunnel. As the volume loss increases from 1% to 3%, the apparent loss of pile capacity increases from 11% to 20%. By increasing the pile length from 0.5 to 1.5 times, the tunnel depth, the maximum pile settlement and lateral deflection decrease by about 63% and 18%, respectively. On the other hand, the maximum bending moment and axial load increase by about 7 and 13 times, respectively. When the pile is located at a distance of 2.5 times the tunnel diameter (Dt), the additional pile responses become insignificant. It was found that an increase in tunnel depth from 1.5Dt to 2.5Dt (with a pile length of 3Dt) increases the maximum lateral deflection by about 420%. Regarding the interaction between tunneling and group of piles, a positive group effect was observed with a significant reduction of the internal forces in rear piles. The maximum bending moment of the front piles was found to be higher than that of the rear piles by about 47%.

Soil is a complex material that shows differently in primary loading, unloading and reloading with stress-dependent stiffness. This general behavior was not possibly being accounted for in simple elastic perfectly plastic Mohr–Coulomb model which is often used to predict the behavior of soils. Thus, in the present study, the more advanced hardening soil model with small-strain stiffness (HSsmall) is used to model the non-linear stress–strain soil behavior. Moreover, unlike previous studies THAT are usually based on the assumption that the soil is homogeneous and using numerical methods by decoupled loadings under plane strain conditions; in this study, the pile responses have been exhaustively investigated in a two-layered soil system using a fully coupled 3D numerical analysis that takes into account the real interactions between tunneling and pile foundations. The paper presents a distinctive set of findings and insights that provide valuable guidance for the design and construction of shield tunnels passing through pile foundations.

The purpose of this paper is to compare the influence of relative density (RD) and strain rates on failure mechanism and specific energy absorption (SEA) of polyamide lattices ranging from bending to stretch-dominated structures using selective laser sintering (SLS).

Three bending and two stretch-dominated unit cells were selected based on the Maxwell stability criterion. Lattices were designed with three RD and fabricated by SLS technique using PA12 material. Quasi-static compression tests with three strain rates were carried out using Taguchi's L9 experiments. The lattice compressive behaviour was verified with the Gibson–Ashby analytical model.

It has been observed that RD and strain rates played a vital role in lattice compressive properties by controlling failure mechanisms, resulting in distinct post-yielding responses as fluctuating and stable hardening in the plateau region. Analysis of variance (ANOVA) displayed the significant impact of RD and emphasised dissimilar influences of strain rate that vary to cell topology. Bending-dominated lattices showed better compressive properties than stretch-dominated lattices. The interesting observation is that stretch-dominated lattices with over-stiff topology exhibited less compressive properties contrary to the Maxwell stability criterion, whereas strain rate has less influence on the SEA of face-centered and body-centered cubic unit cells with vertical and horizontal struts (FBCCXYZ).

This comparative study is expected to provide new prospects for designing end-user parts that undergo various impact conditions like automotive bumpers and evolving techniques like hybrid and functionally graded lattices.

To the best of the authors' knowledge, this is the first work that relates the strain rate with compressive properties and also highlights the lattice behaviour transformation from ductile to brittle while the increase of RD and strain rate analytically using the Gibson–Ashby analytical model.