An empirical expression of cohesion （C） and friction angle （Ф） for layered rock was suggested. This expression was compared with a test result made by the former researchers. The constitutive relationship of a tr...An empirical expression of cohesion （C） and friction angle （Ф） for layered rock was suggested. This expression was compared with a test result made by the former researchers. The constitutive relationship of a transversely isotropic medium and Mohr-Coulomb criterion in which C and Ф vary with directions were employed, and a relative 3D elasto-plastic FEM code was developed, in which the important thing was to adopt a search-trial method to find the orientation angle （p） of shear failure plane （or weakest shear plane） with respect to the major principal stress as well as the corresponding C and Ф Taking an underground opening as the calculation object, the numerical analyses were carried out by using the FEM code for two cases of transversely isotropic rock and isotropic rock, respectively, and the computation results were compared. The results show that when the rock is a transversely isotropic one, the distributions of displacements, plastic zones and stress contours in the surrounding rock will be non-axisymmetric along the tunnel＇s vertical axis, which is very different from that of isotropic rock. The stability of the tunnel in transversely isotropic rock is relatively low.展开更多
This paper is devoted to studying the superconvergence of streamline diffusion finite element methods for convection-diffusion problems. In [8], under the condition that ε ≤ h^2 the optimal finite element error esti...This paper is devoted to studying the superconvergence of streamline diffusion finite element methods for convection-diffusion problems. In [8], under the condition that ε ≤ h^2 the optimal finite element error estimate was obtained in L^2-norm. In the present paper, however, the same error estimate result is gained under the weaker condition that ε≤h.展开更多
An element decomposition method with variance strain stabilization(EDM-VSS) is proposed. In the present EDM-VSS, the quadrilateral element is first divided into four sub-triangular cells, and the local strains in sub-...An element decomposition method with variance strain stabilization(EDM-VSS) is proposed. In the present EDM-VSS, the quadrilateral element is first divided into four sub-triangular cells, and the local strains in sub-triangular cells are obtained using linear interpolation function. For each quadrilateral element, the strain of the whole quadrilateral is the weighted average value of the local strains, which means only one integration point is adopted to construct the stiffness matrix. The stabilization item of the stiffness matrix is constructed by variance of the local strains, which can eliminate the instability of the one-point integration formulation and largely increase the accuracy of the element. Compared with conventional full integration quadrilateral element, the EDM-VSS achieves more accurate results and expends much lower computational cost. More importantly, as no mapping or coordinate transformation is involved in the present EDM-VSS, the restriction on the conventional quadrilateral elements can be removed and problem domain can be discretized in more flexible ways. To verify the accuracy and stability of the present formulation, a number of numerical examples are studied to demonstrate the efficiency of the present EDM-VSS.展开更多
基金Project(2010CB732101) supported by the National Basic Research Program of China Project(51079145) supported by the National Natural Science Foundation of China
文摘An empirical expression of cohesion （C） and friction angle （Ф） for layered rock was suggested. This expression was compared with a test result made by the former researchers. The constitutive relationship of a transversely isotropic medium and Mohr-Coulomb criterion in which C and Ф vary with directions were employed, and a relative 3D elasto-plastic FEM code was developed, in which the important thing was to adopt a search-trial method to find the orientation angle （p） of shear failure plane （or weakest shear plane） with respect to the major principal stress as well as the corresponding C and Ф Taking an underground opening as the calculation object, the numerical analyses were carried out by using the FEM code for two cases of transversely isotropic rock and isotropic rock, respectively, and the computation results were compared. The results show that when the rock is a transversely isotropic one, the distributions of displacements, plastic zones and stress contours in the surrounding rock will be non-axisymmetric along the tunnel＇s vertical axis, which is very different from that of isotropic rock. The stability of the tunnel in transversely isotropic rock is relatively low.
基金Supported by the National Natural Science Foundation of China(10471103)
文摘This paper is devoted to studying the superconvergence of streamline diffusion finite element methods for convection-diffusion problems. In [8], under the condition that ε ≤ h^2 the optimal finite element error estimate was obtained in L^2-norm. In the present paper, however, the same error estimate result is gained under the weaker condition that ε≤h.
基金supported by the National Natural Science Foundation of China(Grant Nos.11472101 and 61232014)Postdoctoral Science Foundation of China(Grant No.2013M531780)the National Laboratory for Electric Vehicles Foundations
文摘An element decomposition method with variance strain stabilization(EDM-VSS) is proposed. In the present EDM-VSS, the quadrilateral element is first divided into four sub-triangular cells, and the local strains in sub-triangular cells are obtained using linear interpolation function. For each quadrilateral element, the strain of the whole quadrilateral is the weighted average value of the local strains, which means only one integration point is adopted to construct the stiffness matrix. The stabilization item of the stiffness matrix is constructed by variance of the local strains, which can eliminate the instability of the one-point integration formulation and largely increase the accuracy of the element. Compared with conventional full integration quadrilateral element, the EDM-VSS achieves more accurate results and expends much lower computational cost. More importantly, as no mapping or coordinate transformation is involved in the present EDM-VSS, the restriction on the conventional quadrilateral elements can be removed and problem domain can be discretized in more flexible ways. To verify the accuracy and stability of the present formulation, a number of numerical examples are studied to demonstrate the efficiency of the present EDM-VSS.
