Search results

The purpose of this paper is to introduce a new method to model a stranded wire efficiently in 3D finite element simulations.

In this method, the stranded wires are numerically approximated with the Cauer ladder network (CLN) model order reduction method in 2D. This approximates the eddy current effect such as the skin and proximity effect for the whole wire. This is then projected to a mesh which does not include each strand. The 3D fields are efficiently calculated with the CLN method and are projected in the 3D geometry to be used in simulations of electrical components with a current vector potential and a homogenized conductivity at each time step.

In applications where the stranded wire geometry is known and does not change, this homogenization approach is an efficient and accurate method, which can be used with any stranded wire configuration, homogenized stranded wire mesh and any input signal dependent on time steps or frequencies.

In comparison to other methods, this method has no direct frequency dependency, which makes the method usable in the time domain for an arbitrary input signal. The CLN can also be used to interconnected stranded cables arbitrarily in electrical components.

Pores and shrink holes are unavoidable defects in the die-casting mass production process which may significantly influence the strength, fatigue and fracture behaviour as well as the life span of structures, especially if they are subjected to high static and dynamic loads. Such defects should be considered during the design process or after production, where the defects could be detected with the help of computed tomography (CT) measurements. However, this is usually not done in today's mass production environments. This paper deals with the stress analysis of die-cast structural parts with pores found from CT measurements or that are artificially placed within a structure.

In this paper the authors illustrate two general methodologies to take into account the porosity of die-cast components in the stress analysis. The detailed geometry of a die-cast part including all discontinuities such as pores and shrink holes can be included via STL data provided by CT measurements. The first approach is a combination of the finite element method (FEM) and the finite cell method (FCM), which extends the FEM if the real geometry cuts finite elements. The FCM is only applied in regions with pores. This procedure has the advantage that all simulations with different pore distributions, real or artificial, can be calculated without changing the base finite element mesh. The second approach includes the pore information as STL data into the original CAD model and creates a new adapted finite element mesh for the simulation. Both methods are compared and evaluated for an industrial problem.

The STL data of defects which the authors received from CT measurements could not be directly applied without repairing them. Therefore, for FEM applications an appropriate repair procedure is proposed. The first approach, which combines the FEM with the FCM, the authors have realized within the commercial software tool Abaqus. This combination performs well, which is demonstrated for test examples, and is also applied for a complex industrial project. The developed in-house code still has some limitations which restrict broader application in industry. The second pure FEM-based approach works well without limitations but requires increasing computational effort if many different pore distributions are to be investigated.

A new simulation approach which combines the FEM with the FCM has been developed and implemented into the commercial Abaqus FEM software. This approach the authors have applied to simulate a real engineering die-cast structure with pores. This approach could become a preferred way to consider pores in practical applications, where the porosity can be derived either from CT measurements or are artificially adopted for design purposes. The authors have also shown how pores can be considered in the standard FEM analysis as well.

The goal of this work is to create a computational finite element model to perform thermoelastic stress analysis (TSA) with the usage of a non-ideal load frequency, containing the effects of the material thermal properties.

Throughout this document, the methodology of the model is presented first, followed by the procedure and results. The last part is reserved to results, discussion and conclusions.

This work had the main goal to create a model to perform TSA with the usage of non-ideal loading frequencies, considering the materials’ thermal properties. Loading frequencies out of the ideal range were applied and the model showed capable of good results. The created model reproduced acceptably the TSA, with the desired conditions.

This work creates a model to perform TSA with the usage of non-ideal loading frequencies, considering the materials’ thermal properties.

This study aims to numerically explore the heat transfer characteristics in turbulent two-degree-of-freedom vortex-induced vibrations (VIVs) of three elastically mounted circular cylinders.

The cylinders are at the vertices of an isosceles triangle with a base and height that are the same. The finite volume technique is used to calculate the Reynolds-averaged governing equations, whereas the structural dynamics equations are solved using the explicit integration method. Simulations are performed for three different configurations, constant mass ratio and natural frequency, as well as distinct reduced velocity values.

As a numerical challenge, the super upper branch observed in the experiment is well-captured by the current numerical simulations. According to the computation findings, the vortex-shedding around the cylinders increases flow mixing and turbulence, hence enhancing heat transfer. At most reduced velocities, the Nusselt number of downstream cylinders is greater than that of upstream cylinders due to the impact of wake-induced vibration, and the maximum heat transfer improvement of these cylinders is 21% (at U_r = 16), 23% (at U_r = 5) and 20% (at U_r = 15) in the first, second and third configurations, respectively.

The main novelty of this study is inspecting the thermal behavior and turbulent flow–induced vibration of three circular cylinders in the triangular arrangement.

Hexahedral meshing is one of the most important steps in performing an accurate simulation using the finite element analysis (FEA). However, the current hexahedral meshing method in the industry is a nonautomatic and inefficient method, i.e. manually decomposing the model into suitable blocks and obtaining the hexahedral mesh from these blocks by mapping or sweeping algorithms. The purpose of this paper is to propose an almost automatic decomposition algorithm based on the 3D frame field and model features to replace the traditional time-consuming and laborious manual decomposition method.

The proposed algorithm is based on the 3D frame field and features, where features are used to construct feature-cutting surfaces and the 3D frame field is used to construct singular-cutting surfaces. The feature-cutting surfaces constructed from concave features first reduce the complexity of the model and decompose it into some coarse blocks. Then, an improved 3D frame field algorithm is performed on these coarse blocks to extract the singular structure and construct singular-cutting surfaces to further decompose the coarse blocks. In most modeling examples, the proposed algorithm uses both types of cutting surfaces to decompose models fully automatically. In a few examples with special requirements for hexahedral meshes, the algorithm requires manual input of some user-defined cutting surfaces and constructs different singular-cutting surfaces to ensure the effectiveness of the decomposition.

Benefiting from the feature decomposition and the 3D frame field algorithm, the output blocks of the proposed algorithm have no inner singular structure and are suitable for the mapping or sweeping algorithm. The introduction of internal constraints makes 3D frame field generation more robust in this paper, and it can automatically correct some invalid 3–5 singular structures. In a few examples with special requirements, the proposed algorithm successfully generates valid blocks even though the singular structure of the model is modified by user-defined cutting surfaces.

The proposed algorithm takes the advantage of feature decomposition and the 3D frame field to generate suitable blocks for a mapping or sweeping algorithm, which saves a lot of simulation time and requires less experience. The user-defined cutting surfaces enable the creation of special hexahedral meshes, which was difficult with previous algorithms. An improved 3D frame field generation method is proposed to correct some invalid singular structures and improve the robustness of the previous methods.

This study aims to assess the robustness and accuracy of the face-centred finite volume (FCFV) method for the simulation of compressible laminar flows in different regimes, using numerical benchmarks.

The work presents a detailed comparison with reference solutions published in the literature –when available– and numerical results computed using a commercial cell-centred finite volume software.

The FCFV scheme provides first-order accurate approximations of the viscous stress tensor and the heat flux, insensitively to cell distortion or stretching. The strategy demonstrates its efficiency in inviscid and viscous flows, for a wide range of Mach numbers, also in the incompressible limit. In purely inviscid flows, non-oscillatory approximations are obtained in the presence of shock waves. In the incompressible limit, accurate solutions are computed without pressure correction algorithms. The method shows its superior performance for viscous high Mach number flows, achieving physically admissible solutions without carbuncle effect and predictions of quantities of interest with errors below 5%.

The FCFV method accurately evaluates, for a wide range of compressible laminar flows, quantities of engineering interest, such as drag, lift and heat transfer coefficients, on unstructured meshes featuring distorted and highly stretched cells, with an aspect ratio up to ten thousand. The method is suitable to simulate industrial flows on complex geometries, relaxing the requirements on mesh quality introduced by existing finite volume solvers and alleviating the need for time-consuming manual procedures for mesh generation to be performed by specialised technicians.

The use of recycled coarse aggregate in concrete structures promotes environmental sustainability; however, performance of these structures might be negatively impacted when it is used as a replacement to traditional aggregate. This paper aims to simulate recycled concrete beams strengthened with carbon fiber-reinforced polymer (CFRP), to advance the modeling and use of recycled concrete structures.

To investigate the performance of beams with recycled coarse aggregate concrete (RCAC), finite element models (FEMs) were developed to simulate 12 preloaded RCAC beams, strengthened with two CFRP strengthening schemes. Details of the modeling are provided including the material models, boundary conditions, applied loads, analysis solver, mesh analysis and computational efficiency.

Using FEM, a parametric study was carried out to assess the influence of CFRP thickness on the strengthening efficiency. The FEM provided results in good agreement with those from the experiments with differences and standard deviation not exceeding 11.1% and 3.1%, respectively. It was found that increasing the CFRP laminate thickness improved the load-carrying capacity of the strengthened beams.

The developed models simulate the preloading and loading up to failure with/without CFRP strengthening for the investigated beams. Moreover, the models were validated against the experimental results of 12 beams in terms of crack pattern as well as load, deflection and strain.

The purpose of this paper is to present a methodology for calculating eddy current losses in the core of a single-phase power voltage transformer, which, unlike a standard power transformer, has an open-type core (I-type core). In those apparatus, reduction of core losses is achieved by using a multipart open-type core that is created by merging a larger number of leaner cores.

3D FEM approach for calculation of eddy current losses in open-type cores based on a weak AλA formulation is presented. Method in which redundant degrees of freedom are eliminated is shown. This enables faster convergence of the simulation. The results are benchmarked using simulations with standard AVA formulation.

Results using weak AλA formulation with elimination of redundant degrees of freedom are in agreement with both simulation using only weak AλA formulation and with simulation based on AVA formulation.

The presented methodology is valid in linear cases, whereas the nonlinear case will be part of future work.

Presented procedure can be used for the optimization when designing the open-type core of apparatus like power voltage transformers.

The presented method is specifically adapted for calculating eddy currents in the open-type core. The method is based on a weak formulation for the magnetic vector potential A and the current vector potential λ, incorporating numerical homogenization and a straightforward elimination of redundant degrees of freedom, resulting in faster convergence of the simulation.

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 ability to characterise highly conducting objects, that may also be highly magnetic, by the complex symmetric rank–2 magnetic polarizability tensor (MPT) is important for metal detection applications including discriminating between threat and non-threat objects in security screening, identifying unexploded anti-personnel landmines and ordnance and identifying metals of high commercial value in scrap sorting. Many everyday non-threat items have both a large electrical conductivity and a magnetic behaviour, which, for sufficiently weak fields and the frequencies of interest, can be modelled by a high relative magnetic permeability. This paper aims to discuss the aforementioned idea.

The numerical simulation of the MPT for everyday non-threat highly conducting magnetic objects over a broad range of frequencies is challenging due to the resulting thin skin depths. The authors address this by employing higher order edge finite element discretisations based on unstructured meshes of tetrahedral elements with the addition of thin layers of prismatic elements. Furthermore, computer aided design (CAD) geometrical models of the non-threat and threat object are often not available and, instead, the authors extract the geometrical features of an object from an imaging procedure.

The authors obtain accurate numerical MPT characterisations that are in close agreement with experimental measurements for realistic physical objects. The assessment of uncertainty shows the impact of geometrical and material parameter uncertainties on the computational results.

The authors present novel computations and measurements of MPT characterisations of realistic objects made of magnetic materials. A novel assessment of uncertainty in the numerical predictions of MPT characterisations for uncertain geometry and material parameters is included.