A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid m...A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.展开更多
This paper presents an adapted stabilisation method for the equal-order mixed scheme of finite elements on convex polygonal meshes to analyse the high velocity and pressure gradient of incompressible fluid flows that ...This paper presents an adapted stabilisation method for the equal-order mixed scheme of finite elements on convex polygonal meshes to analyse the high velocity and pressure gradient of incompressible fluid flows that are governed by Stokes equations system.This technique is constructed by a local pressure projection which is extremely simple,yet effective,to eliminate the poor or even non-convergence as well as the instability of equal-order mixed polygonal technique.In this research,some numerical examples of incompressible Stokes fluid flow that is coded and programmed by MATLAB will be presented to examine the effectiveness of the proposed stabilised method.展开更多
This paper investigates a polygonal finite element(PFE)to solve a two-dimensional(2D)incompressible steady fluid problem in a cavity square.It is a well-known standard benchmark(i.e.,lid-driven cavity flow)-to evaluat...This paper investigates a polygonal finite element(PFE)to solve a two-dimensional(2D)incompressible steady fluid problem in a cavity square.It is a well-known standard benchmark(i.e.,lid-driven cavity flow)-to evaluate the numerical methods in solving fluid problems controlled by the Navier-Stokes(N-S)equation system.The approximation solutions provided in this research are based on our developed equal-order mixed PFE,called Pe1Pe1.It is an exciting development based on constructing the mixed scheme method of two equal-order discretisation spaces for both fluid pressure and velocity fields of flows and our proposed stabilisation technique.In this research,to handle the nonlinear problem of N-S,the Picard iteration scheme is applied.Our proposed method’s performance and convergence are validated by several simulations coded by commercial software,i.e.,MATLAB.For this research,the benchmark is executed with variousReynolds numbers up to the maximum Re=1000.All results then numerously compared to available sources in the literature.展开更多
This paper constructs a new two-dimensional arbitrary polygonal stress hybrid dynamic(APSHD)element for structural dynamic response analysis.Firstly,the energy function is established based on Hamilton's principle...This paper constructs a new two-dimensional arbitrary polygonal stress hybrid dynamic(APSHD)element for structural dynamic response analysis.Firstly,the energy function is established based on Hamilton's principle.Then,the finite element time-space discrete format is constructed using the generalized variational principle and the direct integration method.Finally,an explicit polynomial form of the combined stress solution is give,and its derivation process is shown in detail.After completing the theoretical construction,the numerical calculation program of the APSHD element is written in Fortran,and samples are verified.Models show that the APSHD element performs well in accuracy and convergence.Furthermore,it is insensitive to mesh distortion and has low dependence on selecting time steps.展开更多
The aim of this work is to employ a modified cell-based smoothed finite element method(S-FEM)for topology optimization with the domain discretized with arbitrary polygons.In the present work,the linear polynomial basi...The aim of this work is to employ a modified cell-based smoothed finite element method(S-FEM)for topology optimization with the domain discretized with arbitrary polygons.In the present work,the linear polynomial basis function is used as the weight function instead of the constant weight function used in the standard S-FEM.This improves the accuracy and yields an optimal convergence rate.The gradients are smoothed over each smoothing domain,then used to compute the stiffness matrix.Within the proposed scheme,an optimum topology procedure is conducted over the smoothing domains.Structural materials are distributed over each smoothing domain and the filtering scheme relies on the smoothing domain.Numerical tests are carried out to pursue the performance of the proposed optimization by comparing convergence,efficiency and accuracy.展开更多
In this paper, the stress-strain curve of material is fitted by polygonal line composed of three lines. According to the theory of proportional loading in elastoplasticity, we simplify the complete stress-strain relat...In this paper, the stress-strain curve of material is fitted by polygonal line composed of three lines. According to the theory of proportional loading in elastoplasticity, we simplify the complete stress-strain relations, which are given by the increment theory of elastoplasticity. Thus, the finite element equation with the solution of displacement is derived. The assemblage elastoplastic stiffness matrix can be obtained by adding something to the elastic matrix, hence it will shorten the computing time. The determination of every loading increment follows the von Mises yield criteria. The iterative method is used in computation. It omits the redecomposition of the assemblage stiffness matrix and it will step further to shorten the computing time. Illustrations are given to the high-order element application departure from proportional loading, the computation of unloading fitting to the curve and the problem of load estimation.展开更多
In this work,we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials(FGMs).A recovery type error indicator combined with quadtree decomposition is emp...In this work,we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials(FGMs).A recovery type error indicator combined with quadtree decomposition is employed for adaptive mesh refinement.The proposed approach is capable of capturing the fracture process with a localized mesh refinement that provides notable gains in computational efficiency.The implementation is validated against experimental data and other numerical experiments on orthotropic materials with different material orientations.The results reveal an increase in the stiffness and the maximum force with increasing material orientation angle.The study is then extended to the analysis of orthotropic FGMs.It is observed that,if the gradation in fracture properties is neglected,the material gradient plays a secondary role,with the fracture behaviour being dominated by the orthotropy of the material.However,when the toughness increases along the crack propagation path,a substantial gain in fracture resistance is observed.展开更多
We propose a polygonal finite volume element method based on the mean value coordinates for anisotropic diffusion problems on star-shaped polygonalmeshes.Because the convex cells with hanging nodes are always star-sha...We propose a polygonal finite volume element method based on the mean value coordinates for anisotropic diffusion problems on star-shaped polygonalmeshes.Because the convex cells with hanging nodes are always star-shaped,the computation on themis no longer a problem.Naturally,we apply this advantage of the new polygonal finite volume element method to construct an adaptive polygonal finite volume element algorithm.Moreover,we introduce two refinement strategies,called quadtreebased refinement strategy and polytree-based refinement strategy respectively,and they all have great performance in our numerical tests.The new adaptive algorithm allows the use of hanging nodes,and the number of hanging nodes on each edge is unrestricted in general.Finally,several numerical examples are provided to show the convergence and efficiency of the proposedmethod on various polygonal meshes.The numerical results also show that the newadaptive algorithmnot only reduces the computational cost and the implementation complexity in mesh refinement,but also ensures the accuracy and convergence.展开更多
Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadrati...Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadratic serendipity element shape function is introduced from the linear generalized barycentric coordinates,and the quadratic serendipity element function space based on Wachspress coordinate is selected as the trial function space.Moreover,we construct a family of unified dual partitions for arbitrary convex polygonal meshes,which is crucial to finite volume element scheme,and propose a quadratic serendipity polygonal finite volume element method with fewer degrees of freedom.Finally,under certain geometric assumption conditions,the optimal H1 error estimate for the quadratic serendipity polygonal finite volume element scheme is obtained,and verified by numerical experiments.展开更多
Mesh-based image warping techniques typically represent image deformation using linear functions on triangular meshes or bilinear functions on rectangular meshes.This enables simple and efficient implementation,but in...Mesh-based image warping techniques typically represent image deformation using linear functions on triangular meshes or bilinear functions on rectangular meshes.This enables simple and efficient implementation,but in turn,restricts the representation capability of the deformation,often leading to unsatisfactory warping results.We present a novel,flexible polygonal finite element(poly-FEM)method for content-aware image warping.Image deformation is represented by high-order poly-FEMs on a content-aware polygonal mesh with a cell distribution adapted to saliency information in the source image.This allows highly adaptive meshes and smoother warping with fewer degrees of freedom,thus significantly extending the flexibility and capability of the warping representation.Benefiting from the continuous formulation of image deformation,our polyFEM warping method is able to compute the optimal image deformation by minimizing existing or even newly designed warping energies consisting of penalty terms for specific transformations.We demonstrate the versatility of the proposed poly-FEM warping method in representing different deformations and its superiority by comparing it to other existing state-ofthe-art methods.展开更多
A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal par...A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal particles is used in the discrete element method. Instead of a collision model of circular particles, the collision model used in our method can deal with particles of more complex shape and efficiently simulate the effects of shape on particle–particle and particle–wall interactions. For two particles falling under gravity, because of the edges and corners, different collision patterns for circular and polygonal particles are found in our simulations. The complex vortexes generated near the corners of polygonal particles affect the flow field and lead to a difference in particle motions between circular and polygonal particles. For multiple particles falling under gravity, the polygonal particles easily become stuck owing to their corners and edges, while circular particles slip along contact areas. The present method provides an efficient approach for understanding the effects of particle shape on the dynamics of non-circular particles in fluids.展开更多
文摘A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.
基金The authors would like to present our gratitude to the Flemish Government financially supporting for the VLIR-OUS TEAM Project,VN2017TEA454A103‘An innovative solution to protect Vietnamese coastal riverbanks from floods and erosion’.
文摘This paper presents an adapted stabilisation method for the equal-order mixed scheme of finite elements on convex polygonal meshes to analyse the high velocity and pressure gradient of incompressible fluid flows that are governed by Stokes equations system.This technique is constructed by a local pressure projection which is extremely simple,yet effective,to eliminate the poor or even non-convergence as well as the instability of equal-order mixed polygonal technique.In this research,some numerical examples of incompressible Stokes fluid flow that is coded and programmed by MATLAB will be presented to examine the effectiveness of the proposed stabilised method.
基金This work was supported by the VLIR-UOS TEAM Project,VN2017TEA454A 103,‘An innovative solution to protect Vietnamese coastal riverbanks from floods and erosion’funded by the Flemish Government.
文摘This paper investigates a polygonal finite element(PFE)to solve a two-dimensional(2D)incompressible steady fluid problem in a cavity square.It is a well-known standard benchmark(i.e.,lid-driven cavity flow)-to evaluate the numerical methods in solving fluid problems controlled by the Navier-Stokes(N-S)equation system.The approximation solutions provided in this research are based on our developed equal-order mixed PFE,called Pe1Pe1.It is an exciting development based on constructing the mixed scheme method of two equal-order discretisation spaces for both fluid pressure and velocity fields of flows and our proposed stabilisation technique.In this research,to handle the nonlinear problem of N-S,the Picard iteration scheme is applied.Our proposed method’s performance and convergence are validated by several simulations coded by commercial software,i.e.,MATLAB.For this research,the benchmark is executed with variousReynolds numbers up to the maximum Re=1000.All results then numerously compared to available sources in the literature.
基金funded by the National Natural Science Foundation of China(Grant No.12072135).
文摘This paper constructs a new two-dimensional arbitrary polygonal stress hybrid dynamic(APSHD)element for structural dynamic response analysis.Firstly,the energy function is established based on Hamilton's principle.Then,the finite element time-space discrete format is constructed using the generalized variational principle and the direct integration method.Finally,an explicit polynomial form of the combined stress solution is give,and its derivation process is shown in detail.After completing the theoretical construction,the numerical calculation program of the APSHD element is written in Fortran,and samples are verified.Models show that the APSHD element performs well in accuracy and convergence.Furthermore,it is insensitive to mesh distortion and has low dependence on selecting time steps.
基金support by Basic Science Research Program through the National Research Foundation(NRF)funded by Korea Ministry of Education(No.2016R1A6A1A0312812).
文摘The aim of this work is to employ a modified cell-based smoothed finite element method(S-FEM)for topology optimization with the domain discretized with arbitrary polygons.In the present work,the linear polynomial basis function is used as the weight function instead of the constant weight function used in the standard S-FEM.This improves the accuracy and yields an optimal convergence rate.The gradients are smoothed over each smoothing domain,then used to compute the stiffness matrix.Within the proposed scheme,an optimum topology procedure is conducted over the smoothing domains.Structural materials are distributed over each smoothing domain and the filtering scheme relies on the smoothing domain.Numerical tests are carried out to pursue the performance of the proposed optimization by comparing convergence,efficiency and accuracy.
文摘In this paper, the stress-strain curve of material is fitted by polygonal line composed of three lines. According to the theory of proportional loading in elastoplasticity, we simplify the complete stress-strain relations, which are given by the increment theory of elastoplasticity. Thus, the finite element equation with the solution of displacement is derived. The assemblage elastoplastic stiffness matrix can be obtained by adding something to the elastic matrix, hence it will shorten the computing time. The determination of every loading increment follows the von Mises yield criteria. The iterative method is used in computation. It omits the redecomposition of the assemblage stiffness matrix and it will step further to shorten the computing time. Illustrations are given to the high-order element application departure from proportional loading, the computation of unloading fitting to the curve and the problem of load estimation.
基金E.Martínez-Paneda acknowledges financial support from the Royal Commission for the 1851 Exhibition through their Research Fellowship programme(RF496/2018).
文摘In this work,we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials(FGMs).A recovery type error indicator combined with quadtree decomposition is employed for adaptive mesh refinement.The proposed approach is capable of capturing the fracture process with a localized mesh refinement that provides notable gains in computational efficiency.The implementation is validated against experimental data and other numerical experiments on orthotropic materials with different material orientations.The results reveal an increase in the stiffness and the maximum force with increasing material orientation angle.The study is then extended to the analysis of orthotropic FGMs.It is observed that,if the gradation in fracture properties is neglected,the material gradient plays a secondary role,with the fracture behaviour being dominated by the orthotropy of the material.However,when the toughness increases along the crack propagation path,a substantial gain in fracture resistance is observed.
基金supported by the National Natural Science Foundation of China(Nos.11871009,12271055,12371397)the Foundation of National Key Laboratory of Computational Physics for Young Scholar(No.6142A05QN23008)the Foundation of CAEP(No.CX20210044).
文摘We propose a polygonal finite volume element method based on the mean value coordinates for anisotropic diffusion problems on star-shaped polygonalmeshes.Because the convex cells with hanging nodes are always star-shaped,the computation on themis no longer a problem.Naturally,we apply this advantage of the new polygonal finite volume element method to construct an adaptive polygonal finite volume element algorithm.Moreover,we introduce two refinement strategies,called quadtreebased refinement strategy and polytree-based refinement strategy respectively,and they all have great performance in our numerical tests.The new adaptive algorithm allows the use of hanging nodes,and the number of hanging nodes on each edge is unrestricted in general.Finally,several numerical examples are provided to show the convergence and efficiency of the proposedmethod on various polygonal meshes.The numerical results also show that the newadaptive algorithmnot only reduces the computational cost and the implementation complexity in mesh refinement,but also ensures the accuracy and convergence.
基金supported by the National Natural Science Foundation of China(Nos.11871009,12271055)the Foundation of LCP and the Foundation of CAEP(CX20210044).
文摘Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadratic serendipity element shape function is introduced from the linear generalized barycentric coordinates,and the quadratic serendipity element function space based on Wachspress coordinate is selected as the trial function space.Moreover,we construct a family of unified dual partitions for arbitrary convex polygonal meshes,which is crucial to finite volume element scheme,and propose a quadratic serendipity polygonal finite volume element method with fewer degrees of freedom.Finally,under certain geometric assumption conditions,the optimal H1 error estimate for the quadratic serendipity polygonal finite volume element scheme is obtained,and verified by numerical experiments.
基金The research of Juan Cao was supported by the National Natural Science Foundation of China(Nos.61872308,61972327,and 62272402)the Xiamen Youth Innovation Funds(No.3502Z20206029)Yongjie Jessica Zhang was supported in part by NSF CMMI-1953323 and a Honda grant.
文摘Mesh-based image warping techniques typically represent image deformation using linear functions on triangular meshes or bilinear functions on rectangular meshes.This enables simple and efficient implementation,but in turn,restricts the representation capability of the deformation,often leading to unsatisfactory warping results.We present a novel,flexible polygonal finite element(poly-FEM)method for content-aware image warping.Image deformation is represented by high-order poly-FEMs on a content-aware polygonal mesh with a cell distribution adapted to saliency information in the source image.This allows highly adaptive meshes and smoother warping with fewer degrees of freedom,thus significantly extending the flexibility and capability of the warping representation.Benefiting from the continuous formulation of image deformation,our polyFEM warping method is able to compute the optimal image deformation by minimizing existing or even newly designed warping energies consisting of penalty terms for specific transformations.We demonstrate the versatility of the proposed poly-FEM warping method in representing different deformations and its superiority by comparing it to other existing state-ofthe-art methods.
基金This study was funded by the National Science Foundation of China （Grant No. 11272176）.
文摘A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal particles is used in the discrete element method. Instead of a collision model of circular particles, the collision model used in our method can deal with particles of more complex shape and efficiently simulate the effects of shape on particle–particle and particle–wall interactions. For two particles falling under gravity, because of the edges and corners, different collision patterns for circular and polygonal particles are found in our simulations. The complex vortexes generated near the corners of polygonal particles affect the flow field and lead to a difference in particle motions between circular and polygonal particles. For multiple particles falling under gravity, the polygonal particles easily become stuck owing to their corners and edges, while circular particles slip along contact areas. The present method provides an efficient approach for understanding the effects of particle shape on the dynamics of non-circular particles in fluids.
