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.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
A lightweight aggregate concrete-filled steel tube(LACFST) spatial truss beam was tested under bending load. The performance was studied by the analysis of the beam deflection and strains in its chords and webs. Accor...A lightweight aggregate concrete-filled steel tube(LACFST) spatial truss beam was tested under bending load. The performance was studied by the analysis of the beam deflection and strains in its chords and webs. According to the test results, several assumptions were made to deduce the bearing capacity calculation method based on the force balance of the whole section. An optimal dimension relationship for the truss beam chords was proposed and verified by finite element analysis. Results show that the LACFST spatial truss beam failed after excessive deflection. The strain distribution agreed with Bernoulli-Euler theoretical prediction. The truss beam flexural bearing capacity calculation results matched test evidence with only a 3% difference between the two. Finite element analyses with different chord dimensions show that the ultimate bearing capacity increases as the chord dimensions increase when the chords have a diameter smaller than optimal one; otherwise, it remains almost unchanged as the chord dimensions increase.展开更多
High quality mesh plays an important role for finite element methods in science computation and numerical simulation.Whether the mesh quality is good or not,to some extent,it determines the calculation results of the ...High quality mesh plays an important role for finite element methods in science computation and numerical simulation.Whether the mesh quality is good or not,to some extent,it determines the calculation results of the accuracy and efficiency.Different from classic Lloyd iteration algorithm which is convergent slowly,a novel accelerated scheme was presented,which consists of two core parts:mesh points replacement and local edges Delaunay swapping.By using it,almost all the equilateral triangular meshes can be generated based on centroidal Voronoi tessellation(CVT).Numerical tests show that it is significantly effective with time consuming decreasing by 40%.Compared with other two types of regular mesh generation methods,CVT mesh demonstrates that higher geometric average quality increases over 0.99.展开更多
Shape optimization of mechanical components is one of the issues that have been considered in recent years. Different methods were presented such as adaptive biological for reducing costs and increasing accuracy. The ...Shape optimization of mechanical components is one of the issues that have been considered in recent years. Different methods were presented such as adaptive biological for reducing costs and increasing accuracy. The effects of step factor, the number of control points and the definition way of control points coordinates in convergence rate were studied. A code was written using ANSYS Parametric Design Language (APDL) which receives the studied parameters as input and obtains the optimum shape for the components. The results show that for achieving successful optimization, step factor should be in a specific range. It is found that the use of any coordinate system in defining control points coordinates and selection of any direction for stimulus vector of algorithm will also result in optimum shape. Furthermore, by increasing the number of control points, some non-uniformities are created in the studied boundary. Achieving acceptable accuracy seems impossible due to the creation of saw form at the studied boundary which is called "saw position".展开更多
基金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.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Project(51208176)supported by the National Natural Science Foundation of ChinaProjects(2012M511187,2013T60493)supported by the China Postdoctoral Science FoundationProject(2015B17414)supported by the Fundamental Research Funds for the Central Universities,China
文摘A lightweight aggregate concrete-filled steel tube(LACFST) spatial truss beam was tested under bending load. The performance was studied by the analysis of the beam deflection and strains in its chords and webs. According to the test results, several assumptions were made to deduce the bearing capacity calculation method based on the force balance of the whole section. An optimal dimension relationship for the truss beam chords was proposed and verified by finite element analysis. Results show that the LACFST spatial truss beam failed after excessive deflection. The strain distribution agreed with Bernoulli-Euler theoretical prediction. The truss beam flexural bearing capacity calculation results matched test evidence with only a 3% difference between the two. Finite element analyses with different chord dimensions show that the ultimate bearing capacity increases as the chord dimensions increase when the chords have a diameter smaller than optimal one; otherwise, it remains almost unchanged as the chord dimensions increase.
基金Project(11002121) supported by the National Natural Science Foundation of ChinaProject(09QDZ09) supported by Doctor Foundation of Xiangtan University, China+2 种基金Project(2009LCSSE11) supported by Hunan Key Laboratory for CSSE, ChinaProject(2011FJ3231) supported by Planned Science and Technology Project of Hunan Province,ChinaProject(12JJ3054) supported by the Provincial Natural Science Foundation of Hunan,China
文摘High quality mesh plays an important role for finite element methods in science computation and numerical simulation.Whether the mesh quality is good or not,to some extent,it determines the calculation results of the accuracy and efficiency.Different from classic Lloyd iteration algorithm which is convergent slowly,a novel accelerated scheme was presented,which consists of two core parts:mesh points replacement and local edges Delaunay swapping.By using it,almost all the equilateral triangular meshes can be generated based on centroidal Voronoi tessellation(CVT).Numerical tests show that it is significantly effective with time consuming decreasing by 40%.Compared with other two types of regular mesh generation methods,CVT mesh demonstrates that higher geometric average quality increases over 0.99.
文摘Shape optimization of mechanical components is one of the issues that have been considered in recent years. Different methods were presented such as adaptive biological for reducing costs and increasing accuracy. The effects of step factor, the number of control points and the definition way of control points coordinates in convergence rate were studied. A code was written using ANSYS Parametric Design Language (APDL) which receives the studied parameters as input and obtains the optimum shape for the components. The results show that for achieving successful optimization, step factor should be in a specific range. It is found that the use of any coordinate system in defining control points coordinates and selection of any direction for stimulus vector of algorithm will also result in optimum shape. Furthermore, by increasing the number of control points, some non-uniformities are created in the studied boundary. Achieving acceptable accuracy seems impossible due to the creation of saw form at the studied boundary which is called "saw position".
