Boundary Element Parallel Computation for 3D Elastostatics Using CUDA

2011 ◽  
Author(s):  
Yingjun Wang ◽  
Qifu Wang ◽  
Gang Wang ◽  
Yunbao Huang ◽  
Yixiong Wei
Keyword(s):  
Boundary Element ◽  
Parallel Computation ◽  
Computer Aided Design ◽  
Graphics Processing Unit ◽  
Computation Method ◽  
Processing Unit ◽  
Element Analysis ◽  
Integral Methods ◽  
Cad Models ◽  
Element Method

Finite Element Method (FEM) is pervasively used in most of 3D elastostatic numerical simulations, in which Computer Aided Design (CAD) models need to be converted into mesh models first and then enriched with semantic data (e.g. material parameters, boundary conditions). The interaction between CAD models and FEM models stated above is very intensive. Boundary Element Method (BEM) has been used gradually instead of FEM in recent years because of its advantage in meshing. BEM can reduce the dimensionality of the problem by one so that the complexity in mesh generation can be decreased greatly. In this paper, we present a Boundary Element parallel computation method for 3D elastostatics. The parallel computation runs on Graphics Processing Unit (GPU) using Computing Unified Device Architecture (CUDA). Three major components are included in such method: (1) BEM theory in 3D elastostatics and the boundary element coefficient integral methods, (2) the parallel BEM algorithm using CUDA, and (3) comparison the parallel BEM using CUDA with conventional BEM and FEM respectively by examples. The dimension reduction characteristics of BEM can dispose the 3D elastostatic problem by 2D meshes, therefore we develop a new faceting function to make the ACIS facet meshes suitable for Boundary Element Analysis (BEA). The examples show that the GPU parallel algorithm in this paper can accelerate BEM computation about 40 times.

scholarly journals Boundary-Element Analysis of Magnetic Polarization Tensor for Metallic Cylinder

2021 ◽  
Author(s):  
Zhongwei Jin ◽  
ganghua qin ◽  
haidong fan ◽  
ruochen huang ◽  
ziqi chen ◽  
...  
Keyword(s):  
Boundary Element Method ◽  
Boundary Element ◽  
Scattered Field ◽  
Excitation Frequency ◽  
Computation Method ◽  
Polarization Tensor ◽  
Magnetic Polarization ◽  
Element Analysis ◽  
Metallic Sample ◽  
Element Method

The magnetic polarization tensor has a promising capability of determining the geometry and material properties of metallic samples. In this paper, a novel computation method is proposed to estimate the magnetic polarization tensors for the metallic samples using the boundary element method. In this method, the metallic sample is placed in a uniformly distributed magnetic field. Based on assumptions that the excitation frequency and/or the conductivity of the sample is very high, the metallic sample is regarded as a perfect electrical conductor (PEC). Therefore, the scattered field at a certain distance can be simulated. By utilising the boundary element method, the magnetic polarization tensor can be derived from the simulated scattered field. The theoretical calculation is presented and simulations and experiments have been carried out to validate the proposed method. The results from the simulation are matched with the analytical solution for the case of sphere samples. Moreover, there is a good agreement between the simulation results and the experimental results for the copper cylindrical samples.

Parallel Strategy of FMBEM for 3D Elastostatics and its GPU Implementation Using CUDA

2014 ◽  
Author(s):  
Zhaohui Xia ◽  
Qifu Wang ◽  
Yunbao Huang ◽  
Wei Yixiong ◽  
Wang Yingjun
Keyword(s):  
Boundary Element ◽  
Computer Aided Design ◽  
Large Scale ◽  
Near Field ◽  
Direct Calculation ◽  
Multipole Moment ◽  
Fast Multipole ◽  
Local Expansion ◽  
Cad Models ◽  
Element Method

Finite Element Method (FEM1) is pervasively used in most of 3D product design analysis, in which Computer Aided Design (CAD) models need to be converted in to mesh models first and then enriched with some material features and boundary conditions data, etc. The interaction between CAD models and FEM models is intensive. Boundary Element Method (BEM) has been expected to be advantageous in large-scale problems in recent years owing to its reduction of the dimensionality and its reduced complexity in mesh generation. However, the BEM application has so far been limited to relatively small problems due to the memory and computational complexity for matrix buildup are O(N2). The fast multipole BEM (FMBEM) combined with BEM and fast multipole method (FMM) can overcome the defect of the traditional BEM, and provides an effective method to solve the large-scale problem. Combining GPU parallel computing with FMBEM can further improve its efficiency significantly. Based on the three-dimensional elastic mechanics problems, the parallelisms of the multipole moment (ME), multipole moment to multipole moment (M2M) translation, multipole moment to local expansion (M2L) translation, local expansion to local expansion (L2L) translation and near-field direct calculation were analyzed respectively according to the characteristics of the FMM, and the parallel strategies under CUDA were presented in this paper. Three main major parts are included herein: (1) FMBEM theory in 3D elastostatics, (2) the parallel FMBEM algorithm using CUDA, and (3) comparison the GPU parallel FMBEM with BEM, FEM and FMBEM respectively by engineering examples. Numerical example results show the 3D elastostatics GPU FMBEM not only can speed up the boundary element calculation process, but also save memory which can be effective to solve the large-scale engineering problems.

Using Volume Morphing to Alter Panel Designs to Compensate Shape Distortion in Assembly

2010 ◽  
Vol 10 (2) ◽  
Author(s):  
Ramon F. Sarraga ◽  
Paul A. LeBlanc ◽  
Thomas J. Oetjens
Keyword(s):  
Computer Graphics ◽  
Computer Aided Design ◽  
Computer Software ◽  
Standard Technique ◽  
Free Form ◽  
Annual Conference ◽  
Shape Distortion ◽  
Element Analysis ◽  
Displacement Vectors ◽  
Cad Models

As automotive panels are assembled in a vehicle, they are subjected to shape distorting forces, e.g., the pressure of door seals. A standard technique for preventing shape distortions is to alter the panels’ computer aided design (CAD) in such a way that the panels assume the desired design shape under the action of the distorting forces. Volume morphing, a technique pioneered by Bézier (1978, “General Distortion of an Ensemble of Biparametric Patches,” Comput.-Aided Des., 10(2), pp. 116–120) and by Sederberg and Parry (1986, “Free-Form Deformation of Solid Geometric Models,” International Conference on Computer Graphics and Interactive Techniques, Proceedings of the 13th Annual Conference on Computer Graphics, pp. 151–160), has been extended and implemented in a computer software package called FESHAPE, created at the General Motors Research and Development Center. FESHAPE automatically modifies CAD models according to a finite set of displacement vectors obtained, e.g., from finite-element analysis or scanning tryout parts. This article discusses how FESHAPE has been successfully applied to compensate door panel distortion caused by door seals.

scholarly journals High-performance practical stiffness analysis of high-rise buildings using superfloor elements

2020 ◽  
Vol 7 (2) ◽  
pp. 211-227
Author(s):  
Ahmed A Torky ◽  
Youssef F Rashed
Keyword(s):  
Boundary Element Method ◽  
Boundary Element ◽  
Stiffness Matrix ◽  
High Performance ◽  
Parallel Processors ◽  
Element Formulation ◽  
Processing Unit ◽  
Stiffness Analysis ◽  
Industrial Level ◽  
Element Method

Abstract This study develops a high-performance computing method using OpenACC (Open Accelerator) for the stiffness matrix and load vector generation of shear-deformable plates in bending using the boundary element method on parallel processors. The boundary element formulation for plates in bending is used to derive fully populated displacement-based stiffness matrices and load vectors at degrees of freedom of interest. The computed stiffness matrix of the plate is defined as a single superfloor element and can be solved using stiffness analysis, $Ku = F$, instead of the conventional boundary element method, $Hu = Gt$. Fortran OpenACC code implementations are proposed for the computation of the superfloor element’s stiffness, which includes one serial computing code for the CPU (central processing unit) and two parallel computing codes for the GPU (graphics processing unit) and multicore CPU. As industrial level practical floors are full of supports and geometrical information, the computation time of superfloor elements is reduced dramatically when computing on parallel processors. It is demonstrated that the OpenACC implementation does not affect numerical accuracy. The feasibility and accuracy are confirmed by numerical examples that include real buildings with industrial level structural floors. Engineering computations for massive floors with immense geometrical detail and a multitude of load cases can be modeled as is without the need for simplification.

Boundary Element Analysis of Coupled Continuum and Skeletal Structures

2007 ◽  
Vol 10 (4) ◽  
pp. 415-438
Author(s):  
Youssef F. Rashed
Keyword(s):  
Boundary Element ◽  
Main Idea ◽  
Element Analysis ◽  
Skeletal Structures ◽  
Stiffness Matrices ◽  
Assembly Technique ◽  
A New Technique ◽  
Set Up ◽  
The Continuum ◽  
Element Method

This paper presents a new technique for solving coupled continuum and skeletal structures. The technique is based on employing the well-known flexibility and stiffness methods within the boundary element method (BEM). The analyzed problem is divided into: continuum parts, which are modeled using the BEM and skeletal parts which are modeled using the flexibility or stiffness methods. The main idea of the presented technique is to set up a methodology to generate flexibility or stiffness matrices for the continuum parts using the BEM. To do so, several flexibility and stiffness models are developed. The developed technique is tested on three problems. Results are compared to those obtained from the finite element method (FEM) to show the validity of the developed technique. The present technique gains both advantages of the BEM and the FEM as it allows boundary-only discretization for the continuum parts and uses the banded assembly technique of FEM for the overall structure.

Clearance Measurement of 3D Objects Using Accessibility Cone

2019 ◽  
Author(s):  
Masatomo Inui ◽  
Kouhei Nishimiya ◽  
Nobuyuki Umezu
Keyword(s):  
Graphics Processing Unit ◽  
Three Dimensional ◽  
Computation Time ◽  
Depth Information ◽  
Computation Method ◽  
Processing Unit ◽  
Three Dimensional Objects ◽  
Mechanical Products ◽  
Definition Of ◽  
Graphics Processing

Abstract Clearance is a basic parameter in the design of mechanical products, generally specified as the distance between two shape elements, for example, the width of a slot. This definition is unsuitable for evaluating the clearance during assembly or manufacturing tasks, where the depth information is also critical. In this paper, we propose a novel definition of clearance for the surface of three-dimensional objects. Unlike the typical methods used to define clearance, the proposed method can simultaneously handle the relationship between the width and depth in the clearance, and thus, obtain an intuitive understanding regarding the assembly and manufacturing capability of a product. Our definition is based on the accessibility cone of a point on the object’s surface; further, the peak angle of the accessibility cone corresponds to the clearance at this point. A computation method of the clearance is presented and the results of its application are demonstrated. Our method uses the rendering function of a graphics processing unit to compute the clearance. A large computation time necessary for the analysis is considered as a problem regarding the practical use of this clearance definition.

Using CAD Models and Their Semantics to Prepare F.E. Simulations

2005 ◽  
Author(s):  
Okba Hamri ◽  
Jean-Claude Le´on ◽  
Franca Giannini ◽  
Bianca Falcidieno
Keyword(s):  
Finite Element ◽  
Data Structure ◽  
Computer Aided Design ◽  
Simulation Models ◽  
Product Model ◽  
Element Analysis ◽  
Time Model ◽  
Usual Practice ◽  
Shape Adaptation ◽  
Cad Models

The preparation of simulation models from Computer Aided Design (CAD) models is still a difficult task since shape changes are often required to adapt a component or a mechanical system to the hypotheses and specifications of the simulation model. Detail removal or idealization operations are among the current treatments performed during the preparation of simulation models. Most of the time, model exchanges are required between the engineering office and the simulation engineers, often producing losses of information and lacking of robustness. Thus, inefficient processes and remodelling phases form the usual practice. In this paper we show that geometric models can be extracted from CAD software as well as some of their semantics. This semantics can then be transferred, used and eventually preserved during the shape adaptation process required for a given Finite Element Analysis (FEA). The software environment enabling this transfer simultaneously requires the description of the initial B-Rep NURBS model as well as that of the adapted one. The process set up is based on STandard for the Exchange of Product model data (STEP) to provide a robust link between CAD and shape adaptation environments. In order to describe the appropriate variety of shapes required for the Finite Element (FE) preparation, a specific data structure is proposed to express the corresponding topology of the models. Hence, it is shown that the operators associated to the FE preparation process can take advantage of this data structure and the semantics of the initial CAD model that can be attached to the adapted model. Examples illustrating the various process steps and corresponding operations are provided and demonstrate the robustness of the approach.

Sign in / Sign up

Export Citation Format

Share Document