The Symmetric Galerkin Boundary Element Method is advantageous for the linear elastic fracture and crackgrowth analysis of solid structures,because only boundary and crack-surface elements are needed.However,for engin...The Symmetric Galerkin Boundary Element Method is advantageous for the linear elastic fracture and crackgrowth analysis of solid structures,because only boundary and crack-surface elements are needed.However,for engineering structures subjected to body forces such as rotational inertia and gravitational loads,additional domain integral terms in the Galerkin boundary integral equation will necessitate meshing of the interior of the domain.In this study,weakly-singular SGBEM for fracture analysis of three-dimensional structures considering rotational inertia and gravitational forces are developed.By using divergence theorem or alternatively the radial integration method,the domain integral terms caused by body forces are transformed into boundary integrals.And due to the weak singularity of the formulated boundary integral equations,a simple Gauss-Legendre quadrature with a few integral points is sufficient for numerically evaluating the SGBEM equations.Some numerical examples are presented to verify this approach and results are compared with benchmark solutions.展开更多
A family of stable mixed finite elements for the linear elasticity on tetrahedral grids are constructed,where the stress is approximated by symmetric H(div)-Pk polynomial tensors and the displacement is approximated b...A family of stable mixed finite elements for the linear elasticity on tetrahedral grids are constructed,where the stress is approximated by symmetric H(div)-Pk polynomial tensors and the displacement is approximated by C-1-Pk-1polynomial vectors,for all k 4.The main ingredients for the analysis are a new basis of the space of symmetric matrices,an intrinsic H(div)bubble function space on each element,and a new technique for establishing the discrete inf-sup condition.In particular,they enable us to prove that the divergence space of the H(div)bubble function space is identical to the orthogonal complement space of the rigid motion space with respect to the vector-valued Pk-1polynomial space on each tetrahedron.The optimal error estimate is proved,verified by numerical examples.展开更多
A posteriori error estimators for the symmetric mixed finite element methods for linear elasticity problems with Dirichlet and mixed boundary conditions are proposed. Reliability and efficiency of the estimators are p...A posteriori error estimators for the symmetric mixed finite element methods for linear elasticity problems with Dirichlet and mixed boundary conditions are proposed. Reliability and efficiency of the estimators are proved. Numerical examples are presented to verify the theoretical results.展开更多
The design of mixed finite element methods in linear elasticity with symmetric stress approximations has been a longstanding open problem until Arnold and Winther designed the first family of mixed finite elements whe...The design of mixed finite element methods in linear elasticity with symmetric stress approximations has been a longstanding open problem until Arnold and Winther designed the first family of mixed finite elements where the discrete stress space is the space of H（div,Ω;S）-Pk＋1 tensors whose divergence is a Pk-1 polynomial on each triangle for k ≥ 2. Such a two dimensional family was extended, by Arnold, Awanou and Winther, to a three dimensional family of mixed elements where the discrete stress space is the space of H（div, Ω; S）-Pk＋2 tensors, whose divergence is a Pk-1 polynomial on each tetrahedron for k ≥ 2. In this paper, we are able to construct, in a unified fashion, mixed finite element methods with symmetric stress approximations on an arbitrary simplex in R^n for any space dimension. On the contrary, the discrete stress space here is the space of H（div,Ω; S）-Pk tensors, and the discrete displacement space here is the space of L^2（Ω; R^n）-Pk＋1 vectors for k ≥ n＋ 1. These finite element spaces are defined with respect to an arbitrary simplicial triangulation of the domain, and can be regarded as extensions to any dimension of those in two and three dimensions by Hu and Zhang.展开更多
基金support of the National Natural Science Foundation of China(12072011).
文摘The Symmetric Galerkin Boundary Element Method is advantageous for the linear elastic fracture and crackgrowth analysis of solid structures,because only boundary and crack-surface elements are needed.However,for engineering structures subjected to body forces such as rotational inertia and gravitational loads,additional domain integral terms in the Galerkin boundary integral equation will necessitate meshing of the interior of the domain.In this study,weakly-singular SGBEM for fracture analysis of three-dimensional structures considering rotational inertia and gravitational forces are developed.By using divergence theorem or alternatively the radial integration method,the domain integral terms caused by body forces are transformed into boundary integrals.And due to the weak singularity of the formulated boundary integral equations,a simple Gauss-Legendre quadrature with a few integral points is sufficient for numerically evaluating the SGBEM equations.Some numerical examples are presented to verify this approach and results are compared with benchmark solutions.
基金supported by National Natural Science Foundation of China(Grant Nos.11271035,91430213 and 11421101)
文摘A family of stable mixed finite elements for the linear elasticity on tetrahedral grids are constructed,where the stress is approximated by symmetric H(div)-Pk polynomial tensors and the displacement is approximated by C-1-Pk-1polynomial vectors,for all k 4.The main ingredients for the analysis are a new basis of the space of symmetric matrices,an intrinsic H(div)bubble function space on each element,and a new technique for establishing the discrete inf-sup condition.In particular,they enable us to prove that the divergence space of the H(div)bubble function space is identical to the orthogonal complement space of the rigid motion space with respect to the vector-valued Pk-1polynomial space on each tetrahedron.The optimal error estimate is proved,verified by numerical examples.
基金supported by National Science Foundation of USA(Grant No.DMS-1418934)the Sea Poly Project of Beijing Overseas Talents,National Natural Science Foundation of China(Grant Nos.11625101,91430213,11421101,11771338,11671304 and 11401026)+1 种基金Zhejiang Provincial Natural Science Foundation of China Projects(Grant Nos.LY17A010010,LY15A010015 and LY15A010016)Wenzhou Science and Technology Plan Project(Grant No.G20160019)
文摘A posteriori error estimators for the symmetric mixed finite element methods for linear elasticity problems with Dirichlet and mixed boundary conditions are proposed. Reliability and efficiency of the estimators are proved. Numerical examples are presented to verify the theoretical results.
文摘The design of mixed finite element methods in linear elasticity with symmetric stress approximations has been a longstanding open problem until Arnold and Winther designed the first family of mixed finite elements where the discrete stress space is the space of H（div,Ω;S）-Pk＋1 tensors whose divergence is a Pk-1 polynomial on each triangle for k ≥ 2. Such a two dimensional family was extended, by Arnold, Awanou and Winther, to a three dimensional family of mixed elements where the discrete stress space is the space of H（div, Ω; S）-Pk＋2 tensors, whose divergence is a Pk-1 polynomial on each tetrahedron for k ≥ 2. In this paper, we are able to construct, in a unified fashion, mixed finite element methods with symmetric stress approximations on an arbitrary simplex in R^n for any space dimension. On the contrary, the discrete stress space here is the space of H（div,Ω; S）-Pk tensors, and the discrete displacement space here is the space of L^2（Ω; R^n）-Pk＋1 vectors for k ≥ n＋ 1. These finite element spaces are defined with respect to an arbitrary simplicial triangulation of the domain, and can be regarded as extensions to any dimension of those in two and three dimensions by Hu and Zhang.
