Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA novel method in solving seepage problems implementation in Abaqus based on the polygonal scaled boundary finite element method
126
0
0.0
(
0
)
Added by
Yang Yang
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The scaled boundary finite element method (SBFEM) is a semi-analytical computational scheme, which is based on the characteristics of the finite element method (FEM) and combines the advantages of the boundary element method (BEM). This paper integrates the scaled boundary finite element method (SBFEM) and the polygonal mesh technique into a new approach to solving the steady-state and transient seepage problems. The proposed method is implemented in Abaqus using a user-defined element (UEL). The detailed implementations of the procedure, defining the UEL element, updating the RHS and AMATRX, and solving the stiffness/mass matrix by the eigenvalue decomposition are presented. Several benchmark problems from seepage are solved to validate the proposed implementation. Results show that the polygonal element of PS-SBFEM has a higher accuracy rate than the standard FEM element in the same element size. For the transient problems, the results between PS-SBFEM and the FEM are in excellent agreement. Furthermore, the PS-SBFEM with quadtree meshes shows a good effect for solving complex geometric in the seepage problem. Hence, the proposed method is robust accurate for solving the steady-state and transient seepage problems. The developed UEL source code and the associated input files can be downloaded from https://github.com/yangyLab/PS-SBFEM.
rate research
Read More
This paper presents a steady-state and transient heat conduction analysis framework using the polygonal scaled boundary finite element method (PSBFEM) with polygon/quadtree meshes. The PSBFEM is implemented with commercial finite element code Abaqus by the User Element Sub-routine (UEL) feature. The detailed implementation of the framework, defining the UEL element, and solving the stiffness/mass matrix by the eigenvalue decomposition are presented. Several benchmark problems from heat conduction are solved to validate the proposed implementation. Results show that the PSBFEM is reliable and accurate for solving heat conduction problems. Not only can the proposed implementation help engineering practitioners analyze the heat conduction problem using polygonal mesh in Abaqus, but it also provides a reference for developing the UEL to solve other problems using the scaled boundary finite element method.
The scaled boundary finite element method (SBFEM) is capable of generating polyhedral elements with an arbitrary number of surfaces. This salient feature significantly alleviates the meshing burden being a bottleneck in the analysis pipeline in the standard finite element method (FEM). In this paper, we implement polyhedral elements based on the SBFEM into the commercial finite element software ABAQUS. To this end, user elements are provided through the user subroutine UEL. Detailed explanations regarding the data structures and implementational aspects of the procedures are given. The focus of the current implementation is on interfacial problems and therefore, element-based surfaces are created on polyhedral user elements to establish interactions. This is achieved by an overlay of standard finite elements with negligible stiffness, provided in the ABAQUS element library, with polyhedral user elements. By means of several numerical examples, the advantages of polyhedral elements regarding the treatment of non-matching interfaces and automatic mesh generation are clearly demonstrated. Thus, the performance of ABAQUS for problems involving interfaces is augmented based on the availability of polyhedral meshes. Due to the implementation of polyhedral user elements, ABAQUS can directly handle complex geometries given in the form of digital images or stereolithography (STL) files. In order to facilitate the use of the proposed approach, the code of the UEL is published open-source and can be downloaded from https://github.com/ShukaiYa/SBFEM-UEL.
In this paper, we examine the effectiveness of classic multiscale finite element method (MsFEM) (Hou and Wu, 1997; Hou et al., 1999) for mixed Dirichlet-Neumann, Robin and hemivariational inequality boundary problems. Constructing so-called boundary correctors is a common technique in existing methods to prove the convergence rate of MsFEM, while we think not reflects the essence of those problems. Instead, we focus on the first-order expansion structure. Through recently developed estimations in homogenization theory, our convergence rate is provided with milder assumptions and in neat forms.
We introduce a new efficient algorithm for Helmholtz problems in perforated domains with the design of the scheme allowing for possibly large wavenumbers. Our method is based upon the Wavelet-based Edge Multiscale Finite Element Method (WEMsFEM) as proposed recently in [14]. For a regular coarse mesh with mesh size H, we establish O(H) convergence of this algorithm under the resolution assumption, and with the level parameter being sufficiently large. The performance of the algorithm is demonstrated by extensive 2-dimensional numerical tests including those motivated by photonic crystals.
In this work, we describe a simple finite element approach that is able to resolve weak discontinuities in interface problems accurately. The approach is based on a fixed patch mesh consisting of quadrilaterals, that will stay unchanged independent of the position of the interface. Inside the patches we refine once more, either in eight triangles or in four quadrilaterals, in such a way that the interface is locally resolved. The resulting finite element approach can be considered a fitted finite element approach. In our practical implementation, we do not construct this fitted mesh, however. Instead, the local degrees of freedom are included in a parametric way in the finite element space, or to be more precise in the local mappings between a reference patch and the physical patches. We describe the implementation in the open source C++ finite element library deal.II in detail and present two numerical examples to illustrate the performance of the approach. Finally, detailed studies of the behavior of iterative linear solvers complement this work.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
