Search results

Gives introductory remarks about chapter 1 of this group of 31 papers, from ISEF 1999 Proceedings, in the methodologies for field analysis, in the electromagnetic community. Observes that computer package implementation theory contributes to clarification. Discusses the areas covered by some of the papers ‐ such as artificial intelligence using fuzzy logic. Includes applications such as permanent magnets and looks at eddy current problems. States the finite element method is currently the most popular method used for field computation. Closes by pointing out the amalgam of topics.

The purpose of this paper is to clarify the usefulness of characteristics analysis of an interior permanent magnet (IPM) motor using the 3‐D finite element method (FEM) with prismatic elements from the view point of the accuracy and the calculation time.

The authors analyzed characteristics of an IPM motor applied by the pulse width modulation (PWM) voltage source using the 3‐D FEM with prismatic elements and with tetrahedral elements; then compared the calculated results and the calculation time using prismatic elements with those using tetrahedral elements.

It was found that the calculated current and torque, and eddy current loss using prismatic elements are almost the same as those using tetrahedral elements. The calculation time using prismatic elements is shorter than that using tetrahedral elements.

This paper verifies the usefulness of the 3‐D FEM with prismatic elements through characteristics analysis of an IPM motor driven by the PWM voltage source.

The purpose of this paper is to numerically study blood flow through a subject‐specific carotid artery with a moderately severe stenosis, also to thoroughly analyse the wall shear stress (WSS), oscillatory shear index (OSI) and WSS angular deviation (WSSAD). One of the important aspects of this study is the investigation on the influence of the extensions attached to the domain outlets.

The segmentation of the carotid artery is carried out using a deformable model based on a level set method. A geometric potential force (GPF) is employed to deform the level set to obtain the carotid artery geometry. The initial surface meshing is generated using an advanced marching cubes (MC) method, before improving the quality of the surface mesh via a number of mesh cosmetic steps. The volume mesh generation has two parts. In the first part, a quasi‐structured, boundary layer mesh is generated in the vicinity of the geometry walls. The second part of the meshing involves unstructured tetrahedral meshing of the inner part of the geometry. After the meshing stage, the flow boundary conditions are generated by numerically solving the Helmholtz equation in both space and time. Finally, the explicit characteristic‐based split (CBS) method is employed in a parallel environment to produce a detailed analysis of wall quantities.

In general, WSS is very high in the vicinity of the carotid artery apex and in the proximity of the stenosis. From the results obtained, it is clear that the influence of outlet domain extension is marginal. While the peak instantaneous WSS differs by a maximum of 5.7 per cent, the time‐averaged WSS difference due to extended domain is only 1.3 per cent. Two other derived parameters are also examined in the paper, the oscillating shear index and the WSSAD. Both these quantities also display minor or negligible differences due to domain extension.

It has been perceived that domain extension is essential to avoid wrong application of boundary conditions. The results obtained, however, conclusively show that the outlet domain extension has only a moderate influence on WSS. Thus, outlet extension to the domains may not be essential for arterial blood flows. It is also observed that the dramatic values of peak WSS obtained near the stenosis is the result of high resolution mesh along with boundary layers used in this study. Both the outcomes represent the originality of this paper.

This paper aims to propose a hybrid method, called finite element method‐Dirichlet boundary condition iteration (FEM‐DBCI), for the computation of time‐harmonic eddy current problems inside a conductor heated by coils in 3D open‐boundary geometry.

The method assumes the electrical field as unknown on a mesh of tetrahedral edge elements. The heating power density inside the conductor is then computed and a steady‐state thermal analysis is performed on the same mesh of nodal tetrahedra to calculate the temperature distribution inside the heated piece, taking radiation and convection into account. A numerical example is also provided.

The method couples a differential equation for the interior problem in terms of the electric fields with an integral equation for the exterior one. The global algebraic system is efficiently solved in an iterative way.

The paper illustrates the computation of time‐harmonic eddy current problems inside a conductor heated by coils.

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.

This work aims to describe the physical and numerical modeling of a CFD solver for hypersonic flows in thermo-chemical non-equilibrium. This paper is the second of a two-part series that concerns the application of the solver introduced in Part I to adaptive unstructured meshes.

The governing equations are discretized with an edge-based stabilized finite element method (FEM). Chemical non-equilibrium is simulated using a laminar finite-rate kinetics, while a two-temperature model is used to account for thermodynamic non-equilibrium. The equations for total quantities, species and vibrational-electronic energy conservation are loosely coupled to provide flexibility and ease of implementation. To accurately perform simulations on unstructured meshes, the non-equilibrium flow solver is coupled with an edge-based anisotropic mesh optimizer driven by the solution Hessian to carry out mesh refinement, coarsening, edge swapping and node movement.

The paper shows, through comparisons with experimental and other numerical results, how FEM + anisotropic mesh optimization are the natural choice to accurately simulate hypersonic non-equilibrium flows on unstructured meshes. Three-dimensional test cases demonstrate how, for high-speed flows, shocks resolution, and not necessarily boundary layers resolution, is the main driver of solution accuracy at walls. Equally distributing the error among all elements in a suitably defined Riemannian space yields highly anisotropic grids that feature well-resolved shock waves. The resulting high level of accuracy in the computation of the enthalpy jump translates into accurate wall heat flux predictions. At the opposite end, in all cases examined, high-quality but isotropic unstructured meshes gave very poor solutions with severely inadequate heat flux distributions not even featuring expected symmetries. The paper unequivocally demonstrates that unstructured anisotropically adapted meshes are the best, and may be the only, way for accurate and cost-effective hypersonic flow solutions.

Although many hypersonic flow solvers are developed for unstructured meshes, few numerical simulations on unstructured meshes are presented in the literature. This work demonstrates that the proposed approach can be used successfully for hypersonic flows on unstructured meshes.

The purpose of this study is to model the dynamic characteristics of an opened supersonic disk-gap-band parachute.

A fluid-structure interaction (FSI) method with body-fitted mesh is used to simulate the supersonic parachute. The compressible flow is modeled using large-eddy simulation (LES). A contact algorithm based on the penalty function with a virtual contact domain is proposed to solve the negative volume problem of the body-fitted mesh. Automatic unstructured mesh generation and automatic mesh moving schemes are used to handle complex deformations of the canopy.

The opened disk-gap-band parachute is simulated using Mach 2.0, and the simulation results fit well with the wind tunnel test data. It is found that the LES model can successfully predict large-scale turbulent vortex in the flow. This study also demonstrates the capability of the present FSI method as a tool to predict shock oscillation and breathing phenomenon of the canopy.

The contact algorithm based on the penalty function with a virtual contact domain is proposed for the first time. This methodology can be used to solve the negative volume problem of the dynamic mesh in the flow field.

An easy and robust grid‐based approach is proposed to construct the fully hexahedral mesh in three‐dimensional case and its application for the mesh regeneration or remeshing during the finite element simulation of a metal forming process is presented to show the validity and effectiveness of the scheme. The proposed scheme enables the construction of the provisional mesh by superimposing the regular cubical grid over the object to be meshed and removing the exterior grid points and cells. Because the constructed provisional mesh has the discrete rugged boundary that is quite different from the boundary geometry of the object to be meshed, the nodes on the boundary of the provisional mesh are projected onto the object boundary. The main disadvantage of the mesh constructed by grid‐based approaches is its severely distorted elements on the boundary owing to the projection of the rugged boundary onto the object boundary. In order to improve the quality of boundary elements, some layers of elements on the boundary surface are constructed and the nodes are repositioned by mesh smoothing. Consequently, the quality of boundary elements is effectively improved.

The rostrum of a paddlefish provides hydrodynamic stability during feeding process in addition to detect the food using receptors that are randomly distributed in the rostrum. The exterior tissue of the rostrum covers the cartilage that surrounds the bones forming interlocking star shaped bones.

The aim of this work is to assess the mechanical behavior of four finite element models varying the type of formulation as follows: linear-reduced integration, linear-full integration, quadratic-reduced integration and quadratic-full integration. The paper also presents the load transfer mechanisms of the bone structure of the rostrum. The base material used in the study was steel with elastic–plastic behavior as a homogeneous material before applying materials properties that represents the behavior of bones, cartilages and tissues.

Conclusions are based on comparison among the four models. There is no significant difference between integration orders for similar type of elements. Quadratic-reduced integration formulation resulted in lower structural stiffness compared with linear formulation as seen by higher displacements and stresses than using linearly formulated elements. It is concluded that second-order elements with reduced integration are the alternative to analyze biological structures as they can better adapt to the complex natural contours and can model accurately stress concentrations and distributions without over stiffening their general response.

The use of advanced computational mechanics techniques to analyze the complex geometry and components of the paddlefish rostrum provides a viable avenue to gain fundamental understanding of the proper finite element formulation needed to successfully obtain the system behavior and hot spot locations.

The purpose of this study to analyze the plastic anisotropy of 6061 aluminum alloy with finite deformation using crystal plasticity finite element method.

A representative volume element (RVE) model was constructed by Voronoi tessellation. In this model, grain shapes, grain orientations and distribution of grains were involved. The mechanical response of the component under composite loading was tested using specify cruciform specimen. Moreover, different stress and strain states in the specific central region were analyzed to reveal the effects of complex loading.

We calculated the influence of misorientation of adjacent grains as well as the evolution of the micro structure’s plastic deformation on the macroscopic deformation of the structure. Geometry design for the cruciform specimen helps obtain a homogenous distribution of the stress and strain at the specimen center. In this process, the initial grain orientation is also an important factor, and the larger misorientations between special grains may cause greater stress concentration.

The influence of micro-scale factors on macro-scale plastic anisotropy of AA6061 is analyzed using RVE model and cruciform specimen, and they offer a reference for related research.