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.

Among different types of engineering structures, plates play a significant role. Their analysis necessitates numerical modeling with finite elements, such as triangular, quadrangular or sector plate elements, owing to the intricate geometrical shapes and applied loads. The scope of this study is the development of a new rectangular finite element for thin plate bending based on the strain approach using Airy's function. It is called a rectangular plate finite element using Airy function (RPFEUAF) and has four nodes. Each node had three degrees of freedom: one transverse displacement (w) and two normal rotations (x, y).

Equilibrium conditions are used to generate the interpolation functions for the fields of strain, displacements and stresses. The evolution of the Airy function solutions yielded the selection of these polynomial bi-harmonic functions. The variational principle and the analytical integration approach are used to evaluate the basic stiffness matrix.

The numerical findings for thin plates quickly approach the Kirchhoff solution. The results obtained are compared to the analytical solution based on Kirchhoff theory.

The efficiency of the strain based approach using Airy's function is confirmed, and the robustness of the presented element RPFEUAF is demonstrated. Because of this, the current element is more reliable, better suited for computations and especially intriguing for modeling this kind of structure.

A few earlier studies presented infeasible heatline trajectories for natural convection within annular domains involving an inner circular cylinder and outer square/circular enclosure. The purpose of this paper is to revisit and illustrate the correct heatline trajectories for various test cases.

Galerkin finite element based methodology and space adaptive grid have been used to simulate natural convective flows within the annular domains. The prediction of heatlines involves derivatives at the nodes, which are evaluated based on finite element basis functions and contributions from neighboring elements.

The heatlines in the earlier work indicate infeasible heat flow paths such as heat flow from one portion to the other of isothermal hot walls and heat flow across the adiabatic walls. Current results illustrate physically consistent heat flow paths involving perpendicularly emerging heatlines from hot to cold walls for conductive transport, long heat flow paths around the closed-loop heatline cells for convective transport and parallel layout of heatlines to the adiabatic walls. Results also demonstrate complex heatlines involving multiple flow vortices and complex flow structures.

Current work translates heatfunctions from energy flux vectors, which are determined by using basis sets. This work demonstrates the expected heatline trajectories for various scenarios involving conductive and convective heat transport within enclosures with an inner hot object as a first attempt, and the results are precursors for the understanding of energy flow estimates.

The response of the Pine Flat dam–water–foundation rock system is studied by a new described approach (i.e. FE-(FE-TE)-FE). The initial part of study is focused on the time harmonic analysis. In this part, it is possible to compare the transfer functions against corresponding responses obtained by the FE-(FE-HE)-FE approach (referred to as exact method which employs a rigorous fluid hyper-element). Subsequently, the transient analysis is carried out. In that part, it is only possible to compare the results for low and high normalized reservoir length cases. Therefore, the sensitivity of results is controlled due to normalized reservoir length values.

In the present study, dynamic analysis of a typical concrete gravity dam–water–foundation rock system is formulated by the FE-(FE-TE)-FE approach. In this technique, dam and foundation rock are discretized by plane solid finite elements while, water domain near-field region is discretized by plane fluid finite elements. Moreover, the H-W (i.e. Hagstrom–Warburton) high-order condition is imposed at the reservoir truncation boundary. This task is formulated by employing a truncation element at that boundary. It is emphasized that reservoir far-field is excluded from the discretized model.

High orders of H-W condition, such as O5-5 considered herein, generate highly accurate responses for both possible excitations under both types of full reflective and absorptive reservoir bottom conditions. It is such that transfer functions are hardly distinguishable from corresponding exact responses obtained through the FE-(FE-HE)-FE approach in time harmonic analyses. This is controlled for both low and high normalized reservoir length cases (L/H = 1 and 3). Moreover, it can be claimed that transient analysis leads practically to exact results (in numerical sense) when one is employing high order H-W truncation element. In other words, the results are not sensitive to reservoir normalized length under these circumstances.

Dynamic analysis of concrete gravity dam–water–foundation rock systems is formulated by a new method. The salient aspect of the technique is that it utilizes H-W high-order condition at the truncation boundary. The method is discussed for all types of excitation and reservoir bottom conditions.

Gas metal arc-based directed energy deposition (GMA-DED) process experiences residual stress (RS) developed due to heat accumulation during successive layer deposition as a significant challenge. To address that, monitoring of transient temperature distribution concerning time is a critical input. Finite element analysis (FEA) is considered a decisive engineering tool in quantifying temperature and RS in all manufacturing processes. However, computational time and prediction accuracy has always been a matter of concern for FEA-based prediction of responses in the GMA-DED process. Therefore, this study aims to investigate the effect of finite element mesh variations on the developed RS in the GMA-DED process.

The variation in the element shape functions, i.e. linear- and quadratic-interpolation elements, has been used to model a single-track 10-layered thin-walled component in Ansys parametric design language. Two cases have been proposed in this study: Case 1 has been meshed with the linear-interpolation elements and Case 2 has been meshed with the combination of linear- and quadratic-interpolation elements. Furthermore, the modelled responses are authenticated with the experimental results measured through the data acquisition system for temperature and RS.

A good agreement of temperature and RS profile has been observed between predicted and experimental values. Considering similar parameters, Case 1 produced an average error of 4.13%, whereas Case 2 produced an average error of 23.45% in temperature prediction. Besides, comparing the longitudinal stress in the transverse direction for Cases 1 and 2 produced an error of 8.282% and 12.796%, respectively.

To avoid the costly and time-taking experimental approach, the experts have suggested the utilization of numerical methods in the design optimization of engineering problems. The FEA approach, however, is a subtle tool, still, it faces high computational cost and low accuracy based on the choice of selected element technology. This research can serve as a basis for the choice of element technology which can predict better responses in the thermo-mechanical modelling of the GMA-DED process.

This paper aims to estimate the bearing capacity of a surface strip and circular footings lying on layered sand using numerical limit analysis.

Lower and upper bound limit analysis, as well as finite elements and second-order conic programming (SOCP), are used in this analysis. The yield criterion of Mohr-Coulomb is used to model soil behavior. Using this technique, stringent lower and upper bounds on ultimate bearing capacity can be achieved by assuming an associated flow law.

The obtained results indicate that the exact collapse load is typically being bracketed to within 6% about a mean of both the bounds. The obtained results are compared with the existing literature wherever applicable.

To the best of the authors’ knowledge, no study has used lower and upper bound limit analysis, as well as finite elements and SOCP, to estimate the bearing capacity of a surface strip and circular footings lying on layered sand.

Recently, the utilization of cold-formed steel (CFS) roof truss systems and different types of other combination structural support systems, such as concrete or hot-rolled steel support, becomes more frequently used. This paper aims to identify the load transfer characteristics of three different design details for cold-formed truss to supporting system connections and to propose simplified modelling approach for practices.

Simplification modelling of connection design could be proposed for practical purpose based on the load transfer characteristics obtained from detailed study using finite element method. A cold-formed roof truss system with connection is modelled using line elements. However, the supporting system is not modelled in this work. Three types of connection involve, which are five pieces of CFS L-angle brackets, one-piece of CFS L-angle brackets and three types of bolts connection are modelled.

The results of analysis show that the connections located on the loaded side experienced higher reactions than those far from loaded side. From the result, it is also found that the option of “Fixed But” support condition in STAAD.Pro with translational degree of freedom being restrained is the most suitable way to represent the CFS L-angle brackets design for Type 1 connection for use in truss modelled using line elements.

Such increase in usage necessitates an appropriate connection detailing depending on the behaviour of the connection.

Bearings are critical components used to support loads and facilitate motion for rotating and sliding parts of the machinery. Bearing malfunctions can cause catastrophic failures. Hence, failure analysis and endeavors to improve bearing performance are essential discussions for worldwide designers, manufacturers and end users of vital machinery. This study aims to investigate a type of roller bearing from the railway industry with premature failures. The task arises because locomotives’ maintenance and service life quality are vital to railway operations while providing transportation services for the nation. To assist in maintaining the designated locomotives, the present study scrutinizes the causes of failure of heavy-duty roller bearings from locomotive bogie axleboxes.

It is intended to inspect this bearing service life and statistically scrutinize its design parameters to reveal the failures’ shortcomings and origins. The significant measures include examinations of their failures’ primary and vital factors by comparing them with a real-life service history of 16 roller bearings of the same type. The bearings come from the axleboxes of a locomotive bogie with an axle load of 20 tons. The bearing loads are estimated using the EN13104 standard document and confirmed by the finite element method using ABAQUS engineering software. To validate the finite element modeling results, the bearings’ stress analysis is performed using the Hertzian contact theory that demonstrated perfect conformity. The said methods are also used to search for the areas susceptible to failures in these bearings. With the inclusion and exploitation of the bearing maintenance conditions and the logbook recordings of the locomotives for the past seven years, the critical cause for this type of bearing’s failures is surveyed and discussed.

With the inclusion and exploitation of the bearing maintenance conditions and the logbook recordings of the locomotives for the past seven years, the critical cause for this type of bearing’s failures is surveyed and discussed. As a crucial result, it is found that deprived maintenance and inadequate lubrication are the root causes of the loss of the selected bearings.

For the designated locomotives, the origins of the heavy-duty roller bearing failures and its design shortcomings are revealed by examining and comparing them with a real-life service history of many of the same types of bearings. The novelty of the research is in using the combination of the methods mentioned above and its decent outcome.

When the reinforced concrete beams are reinforced by bonding steel plates to the bottom, excessive use of steel plates will make the reinforced concrete beams become super-reinforced beams, and there are security risks in the actual use of super-reinforced beams. In order to avoid the occurrence of this situation, the purpose of this paper is to study the calculation method of the maximum number of bonded steel plates to reinforce reinforced concrete beams.

First of all, when establishing the limit failure state of the reinforced member, this paper comprehensively considers the role of the tensile steel bar and steel plate and takes the load effect before reinforcement as the negative contribution of the maximum number of bonded steel plates that can be used for reinforcement. Through the definition of the equivalent tensile strength, equivalent elastic modulus and equivalent yield strain of the tensile steel bar and steel plate, a method to determine the relative limit compression zone height of the reinforced member is obtained. Second, based on the maximum ratio of (reinforcement + steel plate), the relative limit compression zone height and the equivalent tensile strength of the tensile steel bar and steel plate of the reinforced member, the calculation method of the maximum number of bonded steel plates is derived. Then, the static load test of the test beam is carried out and the corresponding numerical model is established, and the reliability of the numerical model is verified by comparison. Finally, the accuracy of the calculation method of the maximum number of bonded steel plates is proved by the numerical model.

The numerical simulation results show that when the steel plate width is 800 mm and the thickness is 1–4 mm, the reinforced concrete beam has a delayed yield platform when it reaches the limit state, and the failure mode conforms to the basic stress characteristics of the balanced-reinforced beam. When the steel plate thickness is 5–8 mm, the sudden failure occurs without obvious warning when the reinforced concrete beam reaches the limit state. The failure mode conforms to the basic mechanical characteristics of the super-reinforced beam failure, and the bending moment of the beam failure depends only on the compressive strength of the concrete. The results of the calculation and analysis show that the maximum number of bonded steel plates for reinforced concrete beams in this experiment is 3,487 mm². When the width of the steel plate is 800 mm, the maximum thickness of the steel plate can be 4.36 mm. That is, when the thickness of the steel plate, the reinforced concrete beam is still the balanced-reinforced beam. When the thickness of the steel plate, the reinforced concrete beam will become a super-reinforced beam after reinforcement. The calculation results are in good agreement with the numerical simulation results, which proves the accuracy of the calculation method.

This paper presents a method for calculating the maximum number of steel plates attached to the bottom of reinforced concrete beams. First, based on the experimental research, the failure mode of reinforced concrete beams with different number of steel plates is simulated by the numerical model, and then the result of the calculation method is compared with the result of the numerical simulation to ensure the accuracy of the calculation method of the maximum number of bonded steel plates. And the study does not require a large number of experimental samples, which has a certain economy. The research result can be used to control the number of steel plates in similar reinforcement designs.

This paper aims to focus on the finite element method in the frequency domain (FD-FEM) for the transient electric field in the non-sinusoidal steady state under the non-sinusoidal periodic voltage excitation.

Firstly, the boundary value problem of the transient electric field in the frequency domain is described, and the finite element equation of the FD-FEM is derived by Galerkin’s method. Secondly, the constrained electric field equation on the boundary in the frequency domain (FD-CEFEB) is also derived, which can solve the electric field intensity on the boundary and the dielectric interface with high accuracy. Thirdly, the calculation procedures of the FD-FEM with FD-CEFEB are introduced in detail. Finally, a numerical example of the press-packed insulated gate bipolar transistor under the working condition of the repetitive turn-on and turn-off is given.

The FD-CEFEB improves numerical accuracy of electric field intensity on the boundary and interfacial charge density, which can be achieved by modifying the existing FD-FEMs’ code in appropriate steps. Moreover, the proposed FD-FEM and the FD-CEFEB will only increase calculation costs by a little compared with the traditional FD-FEMs.

The FD-CEFEB can directly solve the electric field intensity on the boundary and the dielectric interface with high accuracy. This paper provides a new FD-FEM for the transient electric field in the non-sinusoidal steady state with high accuracy, which is suitable for combined insulation structure with a long time constant.