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.

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.

The purpose of this paper is to assess the feasibility of analytical models, specifically the radial basis function method, Akbari–Ganji method and Gaussian method, in conjunction with the finite element method. The aim is to examine the impact of processing parameters on temperature history.

Through analytical investigation and finite element simulation, this research examines the influence of processing parameters on temperature history. Simufact software with a thermomechanical approach was used for finite element simulation, while radial basis function, Akbari–Ganji and Gaussian methods were used for analytical modeling to solve the heat transfer differential equation.

The accuracy of both finite element and analytical methods was validated with about 90%. The findings revealed direct relationships between thermal conductivity (from 100 to 200), laser power (from 400 to 800 W), heat source depth (from 0.35 to 0.75) and power absorption coefficient (from 0.4 to 0.8). Increasing the values of these parameters led to higher temperature history. On the other hand, density (from 7,600 to 8,200), emission coefficient (from 0.5 to 0.7) and convective heat transfer (from 35 to 90) exhibited an inverse relationship with temperature history.

The application of analytical modeling, particularly the utilization of the Akbari–Ganji, radial basis functions and Gaussian methods, showcases an innovative approach to studying directed energy deposition. This analytical investigation offers an alternative to relying solely on experimental procedures, potentially saving time and resources in the optimization of DED processes.

The high-pressure accumulator has been widely used in the hydraulic system. Failure pressure prediction is crucial for the safe design and integrity assessment of the accumulators. The purpose of this study is to accurately predict the burst pressure and location for the accumulator shells due to internal pressure.

This study concentrates the non-linear finite element simulation procedure, which allows determination of the burst pressure and crack location using extensive plastic straining criterion. Meanwhile, the full-scale hydraulic burst test and the analytical solution are conducted for comparative analysis.

A good agreement between predicted and measured the burst pressure that was obtained, and the predicted failure point coincided very well with the fracture location of the actual shell very well. Meanwhile, the burst pressure of the shells increases with wall thickness, independent of the length. It can be said that the non-linear finite element method can be employed to predict the failure behavior of a cylindrical shell with sufficient accuracy.

This paper can provide a designer with additional insight into how the pressurized hollow cylinder might fail, and the failure pressure has been predicted accurately with a minimum error below 1%, comparing the numerical results with experimental data.

Due to the miniaturization of electronic devices, the increased current density through solder joints leads to the occurrence of electromigration failure, thereby reducing the reliability of electronic devices. The purpose of this study is to propose a finite element-artificial neural network method for the prediction of temperature and current density of solder joints, and thus provide reference information for the reliability evaluation of solder joints.

The temperature distribution and current density distribution of the interconnect structure of electronic devices were investigated through finite element simulations. During the experimental process, the actual temperature of the solder joints was measured and was used to optimize the finite element model. A large amount of simulation data was obtained to analyze the neural network by varying the height of solder joints, the diameter of solder pads and the magnitude of current loads. The constructed neural network was trained, tested and optimized using this data.

Based on the finite element simulation results, the current is more concentrated in the corners of the solder joints, generating a significant amount of Joule heating, which leads to localized temperature rise. The constructed neural network is trained, tested and optimized using the simulation results. The ANN 1, used for predicting solder joint temperature, achieves a prediction accuracy of 96.9%, while the ANN 2, used for predicting solder joint current density, achieves a prediction accuracy of 93.4%.

The proposed method can effectively improve the estimation efficiency of temperature and current density in the packaging structure. This method prevails in the field of packaging, and other factors that affect the thermal, mechanical and electrical properties of the packaging structure can be introduced into the model.

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.

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.

Beginners and even experienced ones have difficulties in completing the structural fire analysis due to numerical difficulties such as convergence errors and singularity and have to spend a lot of time making many repetitive changes on the model. The aim of this article is to highlight the advantages of explicit solver which can eliminate the mentioned difficulties in finite element analysis containing highly nonlinear contacts, clearance between modeled parts at the beginning and large deflections because of high temperature. This article provides important information, especially for researchers and engineers who are new to structural fire analysis.

The finite element method is utilized to achieve mentioned purposes. First, a comparative study is conducted between implicit and explicit solvers by using Abaqus. Then, a validation process is carried out to illustrate the explicit process by using sequentially coupled heat transfer and structural analysis.

Explicit analysis offers an easier solution than implicit analysis for modeling multi-bolted connections under high temperatures. An optimum mesh density for bolted connections is presented to reflect the realistic structural behavior. Presented explicit process with the offered mesh density is used in the validation of an experimental study on multi-bolted splice connection under ISO 834 standard fire curve. A good agreement is achieved.

What makes the study valuable is that the points to be considered in the structural fire analysis are examined and it is a guide that future researchers can benefit from. This is especially true for modeling and analysis of multi-bolted connections in finite element software under high temperatures. The article can help to shorten and even eliminate the iterative debugging phases, which is a problematic and very time-consuming process for many researchers.

In the analysis of structures in a fire situation by simplified and analytical methods, one assumption is that the fire resistance time is greater than or equal to the required fire resistance time. Among the methodologies involving the fire resistance time, the most used is the tabular method, which associates fire resistance time values to structural elements based on minimum dimensions of the cross section. The tabular method is widely accepted by the technical-scientific community due to the fact that it is safe and practical. However, its main criticism is that it results in lower fire resistance times than advanced thermal and thermostructural analysis methods. The objective of this study was to evaluate the fire resistance time of reinforced concrete beams and compare it with the required fire resistance time recommended by the tabular method of NBR 15200 (ABNT, 2012).

The fire resistance time and required fire resistance time of reinforced concrete beams were evaluated using, respectively, numerical models developed based on the finite element method and the tabular method of NBR 15200 (ABNT, 2012). The influence of the following parameters was investigated: longitudinal reinforcement cover, characteristic compressive strength of concrete, beam height, longitudinal reinforcement area and arrangement of steel bars.

Among the evaluated parameters, the covering of the longitudinal reinforcement proved to be more relevant for the fire resistance time, justifying that the tabular method of NBR 15200 (ABNT, 2012) being strongly and directly influenced by this parameter. In turn, more resistant concretes, higher beams and higher steel grades have lower fire resistance time values. This is because beams in these conditions have greater resistance capacity at room temperature and, consequently, are subject to external stresses of greater magnitude. In some cases, the fire resistance time was even lower than the required fire resistance time prescribed by NBR 15200 (ABNT, 2012). Both the fire resistance time and the required fire resistance time were not influenced by the arrangement of the longitudinal reinforcements.

The present paper innovates by demonstrating the influence of other important design variables on the required fire resistance time of the NBR 15200 (ABNT, 2012). Among several conclusions, it was found that the load level to which the structural elements are subjected considerably affects their fire resistance time. For this reason, it was recommended that the methods for calculating the required fire resistance time consider the load level. In addition, the article quantifies the security degree of the tabular method and exposes some situations for which the tabular method proved to be unsafe. Moreover, in all the models analyzed, the relationship between the span and the vertical deflection associated with the failure of the beams in a fire situation was determined. With this, a span over average deflection relationship was presented in which beams in fire situations fail.

The main objective of this paper is to investigate the response of human skin to an intense temperature drop at the surface. In addition, this paper aims to evaluate the efficiency of finite difference and finite volume methods in solving the highly nonlinear form of Pennes’ bioheat equation.

One-dimensional linear and nonlinear forms of Pennes’ bioheat equation with uniform grids were used to study the behavior of human skin. The specific heat capacity, thermal conductivity and blood perfusion rate were assumed to be linear functions of temperature. The nonlinear form of the bioheat equation was solved using the Newton linearization method for the finite difference method and the Picard linearization method for the finite volume method. The algorithms were validated by comparing the results from both methods.

The study demonstrated the capacity of both finite difference and finite volume methods to solve the one-dimensional and highly nonlinear form of the bioheat equation. The investigation of human skin’s thermal behavior indicated that thermal conductivity and blood perfusion rate are the most effective properties in mitigating a surface temperature drop, while specific heat capacity has a lesser impact and can be considered constant.

This paper modeled the transient heat distribution within human skin in a one-dimensional manner, using temperate-dependent physical properties. The nonlinear equation was solved with two numerical methods to ensure the validity of the results, despite the complexity of the formulation. The findings of this study can help in understanding the behavior of human skin under extreme temperature conditions, which can be beneficial in various fields, including medical and engineering.