Search results

The purpose of this paper is to propose an effective and efficient numerical method that can consider natural frequency in multi-material topology optimization (MMTO) and which is scalable for complex three-dimensional (3D) problems.

The optimization algorithm is developed by combining custom FORTRAN code for MMTO with the open-source software Mystran, which is used as a finite element analysis (FEA) solver. The proposed algorithm allows the designer to shift the fundamental frequency of the design beyond a defined frequency spectrum from the initial designing phase. The methodology is formulated in a smooth and differentiable manner, with the sensitivity expressions, required by gradient-based optimization solvers, presented.

Natural frequency constraint has been successfully implemented into MMTO. The use of open-source software Mystran as an FEA solver in the algorithm provides ability to solve complex problems. Mystran offers powerful built-in functions for eigenvalue extraction using methods like Givens, modified Givens, inverse power and the Lanczos method, which provide the ability to solve complex models. The algorithm is successfully able to solve both two- and three-material MMTO jobs for two-dimensional and 3D geometries.

Natural frequency constraint consideration into topology optimization is very challenging due to three common issues: localized eigenmodes, mode switching and high computational cost. The proposed algorithm addresses these inherent issues, implements natural frequency constraint to MMTO and solves for complex models, which is hardly possible using conventional methods.

The purpose of this paper is to present a computationally efficient approach for the stochastic modeling of an inhomogeneous reluctivity of magnetic materials. These materials can be part of electrical machines such as a single-phase transformer (a benchmark example that is considered in this paper). The approach is based on the Karhunen–Loève expansion (KLE). The stochastic model is further used to study the statistics of the self-inductance of the primary coil as a quantity of interest (QoI).

The computation of the KLE requires solving a generalized eigenvalue problem with dense matrices. The eigenvalues and the eigenfunction are computed by using the Lanczos method that needs only matrix vector multiplications. The complexity of performing matrix vector multiplications with dense matrices is reduced by using hierarchical matrices.

The suggested approach is used to study the impact of the spatial variability in the magnetic reluctivity on the QoI. The statistics of this parameter are influenced by the correlation lengths of the random reluctivity. Both, the mean value and the standard deviation increase as the correlation length of the random reluctivity increases.

The KLE, computed by using hierarchical matrices, is used for uncertainty quantification of low frequency electrical machines as a computationally efficient approach in terms of memory requirement, as well as computation time.

The paper presents the numerical performance of the preconditioned generalized conjugate gradient (PGCG) methods in solving non‐linear convection — diffusion equations. Three non‐linear systems which describe a non‐isothermal chemical reactor, the chemically driven convection in a porous medium and the incompressible steady flow past a sphere are the test problems. The standard second order accurate centred finite difference scheme is used to discretize the models equations. The discrete approximations are solved with a double iterative process using the Newton method as outer iteration and the PGCG algorithm as inner iteration. Three PGCG techniques, which emerge to be the best performing, are tested. Laplace‐type operators are employed for preconditioning. The results show that the convergence of the PGCG methods depends strongly on the convection—diffusion ratio. The most robust algorithm is GMRES. But even with GMRES non‐convergence occurs when the convection—diffusion ratio exceeds a limit value. This value seems to be influenced by the non‐linearity type.

This paper aims to numerically study the effects of boundary conditions, pre-stress, material constants and thickness on the dynamic performance of a wrinkled thin membrane.

Based on the stability theory of plates and shells, the dynamic equations of a wrinkled thin membrane were developed, and they were solved with the Lanczos method

The effects of wrinkle-influencing factors on the dynamic performance of a wrinkled membrane are determined by the wrinkling stage. The effects are prominent when wrinkling deformation is evolving, but they are very small and can hardly be observed when wrinkling deformation is stable. Mode shapes of a wrinkled membrane are sensitive to boundary conditions, pre-stress and Poisson’s ratio, but its natural frequencies are sensitive to all these five factors.

The research work in this paper is expected to help understand the dynamic behavior of a wrinkled membrane and present access to ensuring its dynamic stability by controlling the wrinkle-influencing factors.

Very few documents investigated the dynamic properties of wrinkled membranes. No attention has yet been paid by the present literature to the global dynamic performance of a wrinkled membrane under the influences of the factors that play a pivotal role in the wrinkling deformation. In view of this, this paper numerically studied the global modes and corresponding frequencies of a wrinkled membrane and their variation with the wrinkle-influencing factors. The results indicate that the global dynamic properties of a wrinkled membrane are sensitive to these factors at the stage of wrinkling evolution.

The authors examined the numerical natural frequency analysis of a 2D functionally graded (FG) truncated thick hollow cone using 3D elasticity theory.

The material properties of the 2D-FGM (two dimensional-functionally graded materials) cone are graded along the radial and axial axes of the cone using a power–law distribution. The eigenvalue problem was solved using finite element analysis (FEA) employing graded hexahedral elements, and the verification of the finite element approach was assessed by comparing the current solution to earlier experimental studies.

The effects of semivertex angle, material distribution and the cone configuration on the natural frequencies have been analyzed. For various semivertex angles, thickness, length and power law exponents, many results in the form of natural frequencies and mode shapes are presented for the 2D-FGM cone. As a result, the effects of the given parameters were addressed, and the results were compared, demonstrating the direct efficiency of raising the power–law exponents and cone thickness on the rise of natural frequencies.

For the first time, the numerical natural frequency analysis of a 2D-FG truncated thick hollow truncated cone based on 3D equations of elasticity has been investigated. The material properties of the truncated cone have been distributed along two directions, which has not been considered before in any research for the truncated thick cone. The reason for using these innovative volume fraction functions is the lack of accurate coverage by functions that are available in the literature (Asemi et al., 2011; Babaei et al. 2021).

Computational algorithms for finite element dynamic analysis of large‐scale structural problems that exploit both concurrent and parallel features of multiple instruction multiple data streams computers are presented. A new computer program architecture is used in which large finite element domains are automatically divided into subdomains. The number of subdomains generated is equal to the number of available processors. The spatial solution is obtained using a basis of orthogonal vectors. The temporal solution is computed exactly. Discussion is focused on the concurrent generation of global Ritz vectors. Examples run on a hypercube multiprocessor confirm the potential of the proposed scheme.

As a practical engineering method, earthquake response spectra play an important role in seismic hazard assessment and in seismic design of structures. However, the computing scale and the efficiency of commercial software restricted the solution of complex structures. There is a clear need of developing large-scale and highly efficient finite element procedures for response spectrum analysis.

In this paper, the kernel theories for earthquake response spectra are deduced and the corresponding parallel solution flow via the modal superposition method is presented. Based on the algorithm and the parallel data structure of JAUMIN framework, a parallel finite element (FE) solution module is established. Using the solution procedure on a supercomputer equipped with up to thousands of processors, the correctness and parallel scalability of the algorithm are evaluated via numerical experiments of typical engineering examples.

The results show that the solution module has the same precision as the commercial FE software ANSYS; the maximum solution scale achieves 154 million degrees of freedom (DOFs) with a favorable parallel computing efficiency, going far beyond the computing ability of the commercial FE software.

The solution scale in this paper is very challenging for the large-scale parallel computing of structural dynamics and will promote the dynamic analysis ability of complex facilities greatly.

The paper illustrates how software designed and constructed in a truly modular style can rapidly be adapted to exploit new developments in analytical techniques. Given the increasing cost of writing new software, it is argued that this concept has commercial, as well as academic, attractiveness.

A hierarchically preconditioned conjugate gradient (PCG) method for finite element analysis is presented. Its use is demonstrated for the difficult problem of the non‐linear analysis of 3D reinforced concrete structures. Examples highlight the dramatic savings in computer storage and more modest savings in solution times obtained using PCG especially for large problems.

This study aims to reduce the weight of the door, improve the operating efficiency and ensure the safety of vehicle operation.

Based on traditional aluminium alloy doors, a new type of honeycomb composite material was developed. Tests were conducted to determine the honeycomb compression resistance, honeycomb and skin shear performance, plate bending, thermal conductivity and environmental protection. Eight doors were developed based on the full-side open structure, and static strength and stiffness analyses were performed simultaneously. To solve door vibration problems, modal analysis and test were carried out.

The test results showed that the weight of the door was reduced by more than 40% whilst ensuring the strength and stiffness of the vehicle. The first–sixth-order test mode of the door was increased by more than 14% compared with existing aluminium alloy doors.

A new type of honeycomb composite material was used in this study. The test results showed that the weight of the door was reduced by more than 40% whilst ensuring the strength and stiffness of the vehicle. The 1st-to-6th order test mode of the door was increased by more than 14% compared with the existing aluminium alloy door.