Search results

Structural topology optimization is computationally expensive due to the involvement of high-resolution mesh and repetitive use of finite element analysis (FEA) for computing the structural response. Since FEA consumes most of the computational time in each optimization iteration, a novel GPU-based parallel strategy for FEA is presented and applied to the large-scale structural topology optimization of 3D continuum structures.

A matrix-free solver based on preconditioned conjugate gradient (PCG) method is proposed to minimize the computational time associated with solution of linear system of equations in FEA. The proposed solver uses an innovative strategy to utilize only symmetric half of elemental stiffness matrices for implementation of the element-by-element matrix-free solver on GPU.

Using solid isotropic material with penalization (SIMP) method, the proposed matrix-free solver is tested over three 3D structural optimization problems that are discretized using all hexahedral structured and unstructured meshes. Results show that the proposed strategy demonstrates 3.1× –3.3× speedup for the FEA solver stage and overall speedup of 2.9× –3.3× over the standard element-by-element strategy on the GPU. Moreover, the proposed strategy requires almost 1.8× less GPU memory than the standard element-by-element strategy.

The proposed GPU-based matrix-free element-by-element solver takes a more general approach to the symmetry concept than previous works. It stores only symmetric half of the elemental matrices in memory and performs matrix-free sparse matrix-vector multiplication (SpMV) without any inter-thread communication. A customized data storage format is also proposed to store and access only symmetric half of elemental stiffness matrices for coalesced read and write operations on GPU over the unstructured mesh.

The purpose of this study is to develop a reliable finite element algorithm based on the transmission line method (TLM) to solve the nonlinear problem in electromagnetic field calculation.

In this paper, the TLM has been researched and applied to solve nonlinear iteration in FEM. LU decomposition method and the Jacobi preconditioned conjugate gradient method have been investigated to solve the equations in transmission line finite element method (FEM-TLM). The algorithms have been developed in C++ language. The algorithm is applied to analyze the magnetic field of a long straight current-carrying wire and a single-phase transformer.

FEM-TLM is more effective than traditional FEM in nonlinear iteration. The results of FEM-TLM have been compared and analyzed under different calculation scales. It is found that the LU decomposition method is more suitable for FEM-TLM because there is no need to repeatedly assemble the global coefficient matrix in the iterative solution process and it is not affected by the disturbance of the right-hand vector.

An effective algorithm is provided for solving nonlinear problems in the electromagnetic field, which can save a lot of computing costs. The efficiency of LU decomposition and CG method in FEM-TLM nonlinear iteration is investigated, which also makes a preliminary exploration for the research of FEM-TLM parallel algorithms.