Search results

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 present paper aims to analyse the synergistic effect of pocket orientation and piezo-viscous-polar (PVP) lubrication on the performance of multi-recessed hybrid journal bearing (MHJB) system.

To simulate the behaviour of PVP lubricant in clearance space of the MHJB system, the modified form of Reynolds equation is numerically solved by using finite element method. Galerkin’s method is used to obtain the weak form of the governing equation. The system equation is solved by Gauss–Seidal iterative method to compute the unknown values of nodal oil film pressure. Subsequently, performance characteristics of bearing system are computed.

The simulated results reveal that the location of pressurised lubricant inlets significantly affects the oil film pressure distribution and may cause a significant effect on the characteristics of bearing system. Further, the use of PVP lubricant may significantly enhances the performance of the bearing system, namely.

The present work examines the influence of pocket orientation with respect to loading direction on the characteristics of PVP fluid lubricated MHJB system and provides vital information regarding the design of journal bearing system.

The peer review history for this article is available at: https://publons.com/publon/10.1108/ILT-07-2023-0241/

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.

1) The OE-MLPG penalty meshfree method is developed to solve cracked structure.(2) Smartening the numerical meshfree method by combining the particle swarm optimization (PSO) optimization algorithms and Voronoi computational geometric algorithm. (3). Selection of base functions, finding optimal penalty factor and distribution of appropriate nodal points to the accuracy of calculation in the meshless local Petrov–Galekrin (MLPG) meshless method.

Using appropriate shape functions and distribution of nodal points in local domains and sub-domains and choosing an approximation or interpolation method has an effective role in the application of meshless methods for the analysis of computational fracture mechanics problems, especially problems with geometric discontinuity and cracks. In this research, computational geometry technique, based on the Voronoi diagram (VD) and Delaunay triangulation and PSO algorithm, are used to distribute nodal points in the sub-domain of analysis (crack line and around it on the crack plane).

By doing this process, the problems caused by too closeness of nodal points in computationally sensitive areas that exist in general methods of nodal point distribution are also solved. Comparing the effect of the number of sentences of basic functions and their order in the definition of shape functions, performing the mono-objective PSO algorithm to find the penalty factor, the coefficient, convergence, arrangement of nodal points during the three stages of VD implementation and the accuracy of the answers found indicates, the efficiency of V-E-MLPG method with Ns = 7 and ß = 0.0037–0.0075 to estimation of 3D-stress intensity factors (3D-SIFs) in computational fracture mechanics.

The present manuscript is a continuation of the studies (Ref. [33]) carried out by the authors, about; feasibility assessment, improvement and solution of challenges, introduction of more capacities and capabilities of the numerical MLPG method have been used. In order to validate the modeling and accuracy of calculations, the results have been compared with the findings of reference article [34] and [35].

This study establishes a method for predicting ground vibrations caused by railway tunnels in unsaturated soils with spatial variability.

First, an improved 2.5D finite-element-method-perfect-matching-layer (FEM-PML) model is proposed. The Galerkin method is used to derive the finite element expression in the ub-pl-pg format for unsaturated soil. Unlike the ub-v-w format, which has nine degrees of freedom per node, the ub-pl-pg format has only five degrees of freedom per node; this significantly enhances the calculation efficiency. The stretching function of the PML is adopted to handle the unlimited boundary domain. Additionally, the 2.5D FEM-PML model couples the tunnel, vehicle and track structures. Next, the spatial variability of the soil parameters is simulated by random fields using the Monte Carlo method. By incorporating random fields of soil parameters into the 2.5D FEM-PML model, the effect of soil spatial variability on ground vibrations is demonstrated using a case study.

The spatial variability of the soil parameters primarily affected the vibration acceleration amplitude but had a minor effect on its spatial distribution and attenuation over time. In addition, ground vibration acceleration was more affected by the spatial variability of the soil bulk modulus of compressibility than by that of saturation.

Using the 2.5D FEM-PML model in the ub-pl-pg format of unsaturated soil enhances the computational efficiency. On this basis, with the random fields established by Monte Carlo simulation, the model can calculate the reliability of soil dynamics, which was rarely considered by previous models.

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.

In the present work, we focus on developing an in-house parallel meshless local Petrov-Galerkin (MLPG) code for the analysis of heat conduction in two-dimensional and three-dimensional regular as well as complex geometries.

The parallel MLPG code has been implemented using open multi-processing (OpenMP) application programming interface (API) on the shared memory multicore CPU architecture. Numerical simulations have been performed to find the critical regions of the serial code, and an OpenMP-based parallel MLPG code is developed, considering the critical regions of the sequential code.

Based on performance parameters such as speed-up and parallel efficiency, the credibility of the parallelization procedure has been established. Maximum speed-up and parallel efficiency are 10.94 and 0.92 for regular three-dimensional geometry (343,000 nodes). Results demonstrate the suitability of parallelization for larger nodes as parallel efficiency and speed-up are more for the larger nodes.

Few attempts have been made in parallel implementation of the MLPG method for solving large-scale industrial problems. Although the literature suggests that message-passing interface (MPI) based parallel MLPG codes have been developed, the OpenMP model has rarely been touched. This work is an attempt at the development of OpenMP-based parallel MLPG code for the very first time.

This study aims to propose and numerically assess different ways of discretising a very weak formulation of the Poisson problem.

We use integration by parts twice to shift smoothness requirements to the test functions, thereby allowing low-regularity data and solutions.

Various conforming discretisations are presented and tested, with numerical results indicating good accuracy and stability in different types of problems.

This is one of the first articles to propose and test concrete discretisations for very weak variational formulations in primal form. The numerical results, which include a problem based on real MRI data, indicate the potential of very weak finite element methods for tackling problems with low regularity.

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.

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.