This paper studies high order compact finite volume methods on non-uniform meshes for one-dimensional elliptic and parabolic differential equations with the Robin boundary conditions.An explicit scheme and an implicit...This paper studies high order compact finite volume methods on non-uniform meshes for one-dimensional elliptic and parabolic differential equations with the Robin boundary conditions.An explicit scheme and an implicit scheme are obtained by discretizing the equivalent integral form of the equation.For the explicit scheme with nodal values,the algebraic system can be solved by the Thomas method.For the implicit scheme with both nodal values and their derivatives,the system can be implemented by a prediction-correction procedure,where in the correction stage,an implicit formula for recovering the nodal derivatives is introduced.Taking two point boundary value problem as an example,we prove that both the explicit and implicit schemes are convergent with fourth order accuracy with respect to some standard discrete norms using the energy method.Two numerical examples demonstrate the correctness and effectiveness of the schemes,as well as the indispensability of using non-uniform meshes.展开更多
A new model,called object model,for the simulation of cold roll-forming of tubes is presented.The model inherits the advantages of old models and is the embodiment of forming process that the strip is rolled step by s...A new model,called object model,for the simulation of cold roll-forming of tubes is presented.The model inherits the advantages of old models and is the embodiment of forming process that the strip is rolled step by step from feed rollers to last rolling pass.The elastic-plastic large deformation spline finite strip method based on updated Lagrangian method has been developed by improving the stiffness and transition matrix.Combined theory formulas and new analytical model,the forming process of a tube has been simulated successfully as an example.The analytical results are submitted and indicate that the proposed simulation method and new model are applicable.展开更多
In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference me...In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.展开更多
A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D F...A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D FEM of 7 nodes and 10 nodes are constructed based on the basic formula. Using these proposed elements, the multiscale numerical model for foundation subjected to harmonic periodic load, the foundation model excited by external and internal dynamic load are studied. The results show the pro- posed finite elements have higher precision than the tradi- tional elements with 4 nodes. The proposed finite elements can describe the propagation of stress waves well whenever the foundation model excited by extemal or intemal dynamic load. The proposed finite elements can be also used to con- nect the multi-scale elements. And the proposed finite elements also have high precision to make multi-scale analysis for structure.展开更多
In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant...In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant. Polygonal elements can do well in simulation of the materials behavior and provide greater flexibility for the meshing of complex geometries. Hence, the study on the polygonal element is a very useful and necessary part in the finite element method. In this paper, an n-sided polygonal element based on quadratic spline interpolant, denoted by PS2 element, is presented using the triangular area coordinates and the B-net method. The PS2 element is conforming and can exactly model the quadratic field. It is valid for both convex and non-convex polygonal element, and insensitive to mesh distortions. In addition, no mapping or coordinate transformation is required and thus no Jacobian matrix and its inverse are evaluated. Some appropriate examples are employed to evaluate the performance of the proposed element.展开更多
Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions...Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.展开更多
The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh disto...The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previous work, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the B- net method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian co- ordinates. In this paper, a thin plate spline element is devel- oped based on the spline element L8 and the refined tech- nique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.展开更多
The near crack line analysis method was used to investigate a centric crack loaded by two pairs of point shear forces in a finite plate, and the analytical solution was obtained. The solution includes the unit normal ...The near crack line analysis method was used to investigate a centric crack loaded by two pairs of point shear forces in a finite plate, and the analytical solution was obtained. The solution includes the unit normal vector of the elastic-plastic boundary near the crack line, the elastic-plastic stress fields near the crack line, the variations of the length of the plastic zone along the crack line with an external load, and the bearing capacity of a finite plate with a centric crack loaded by two pairs of point shear forces. The results are sufficiently precise near the crack line because the assumptions of small scale yielding theory have not been made and no other assumptions are taken.展开更多
In this work, a design procedure extending the B-spline based finite cell method into shape optimization is developed for axisymmetric solids involving the centrifugal force effect. We first replace the traditional co...In this work, a design procedure extending the B-spline based finite cell method into shape optimization is developed for axisymmetric solids involving the centrifugal force effect. We first replace the traditional conforming mesh in the finite element method with structured cells that are fixed during the whole design process with a view to avoid the sophisticated re-meshing and eventual mesh distortion.Then, B-spline shape functions are further implemented to yield a high-order continuity field along the cell boundary in stress analysis. By means of the implicit description of the shape boundary, stress sensitivity is analytically derived with respect to shape design variables. Finally, we illustrate the efficiency and accuracy of the proposed protocol by several numerical test cases as well as a whole design procedure carried out on an aeronautic turbine disk.展开更多
In this paper,the spline finite element method is developed to investigate free vibration problems of beams.The cubic B-spline functions are used to construct the displacement field.The assembly of elements and the in...In this paper,the spline finite element method is developed to investigate free vibration problems of beams.The cubic B-spline functions are used to construct the displacement field.The assembly of elements and the introduction of boundary conditions follow the standard finite element procedure.The results under various boundary conditions are compared with those obtained by the exact method and the finite difference method.It shows that the results are in excellent agreement with the analytical results and much more accurate than the results obtained by the finite difference method,especially for higher order modes.展开更多
In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, w...In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.展开更多
The choice of self-concordant functions is the key to efficient algorithms for linear and quadratic convex optimizations, which provide a method with polynomial-time iterations to solve linear and quadratic convex opt...The choice of self-concordant functions is the key to efficient algorithms for linear and quadratic convex optimizations, which provide a method with polynomial-time iterations to solve linear and quadratic convex optimization problems. The parameters of a self-concordant barrier function can be used to compute the complexity bound of the proposed algorithm. In this paper, it is proved that the finite barrier function is a local self-concordant barrier function. By deriving the local values of parameters of this barrier function, the desired complexity bound of an interior-point algorithm based on this local self-concordant function for linear optimization problem is obtained. The bound matches the best known bound for small-update methods.展开更多
For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan＇s element energy projection （EEP） method has ...For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan＇s element energy projection （EEP） method has the accuracy O（h^min{2k,k＋4}） The theoretical analysis coincides the reported numerical results.展开更多
In this paper,a 13-node pyramid spline element is derived by using the tetrahedron volume coordinates and the B-net method,which achieves the second order completeness in Cartesian coordinates.Some appropriate example...In this paper,a 13-node pyramid spline element is derived by using the tetrahedron volume coordinates and the B-net method,which achieves the second order completeness in Cartesian coordinates.Some appropriate examples were employed to evaluate the performance of the proposed element.The numerical results show that the spline element has much better performance compared with the isoparametric serendipity element Q20 and its degenerate pyramid element P13 especially when mesh is distorted,and it is comparable to the Lagrange element Q27.It has been demonstrated that the spline finite element method is an efficient tool for developing high accuracy elements.展开更多
Tlie rigid-flexible coupling dynamic modeling theory and the discretization methods of a rotating flexible rectangular thin plate are investigated in this paper.Based on the continuum mechanics,the rigid-flexible coup...Tlie rigid-flexible coupling dynamic modeling theory and the discretization methods of a rotating flexible rectangular thin plate are investigated in this paper.Based on the continuum mechanics,the rigid-flexible coupling dynamic model is established for the flexible rectangular thin plate undergoing large overall rotation,and the coupling term of the deformation which is caused by transverse deformation is considered.Assumed mode method(AMM),spline finite point method(SFPM)and Beizer finite point method(BFPM)are used to describe the deformation of the flexible rectangular plate,and then the dynamic equations of a rotating flexible rectangular thin plate undergoing overall motion are derived by Lagrange^equation of the second kind.The dynamics of a cantilever plate undergoing large overall rotation is simulated via using different dynamic models,and the simulation results of the first order approximation model are compared with those of the traditional zero-order approximation model.It is shown that the first order approximation model with the dynamic stiffening terms can correctly describe the dynamic behavior of the system undergoing large overall rotation,while the zero-order approximation model cannot get the correct results.And AMM.SFPM.BFPM can well describe the deformation of a rotating flexible rectangular plate.展开更多
文摘This paper studies high order compact finite volume methods on non-uniform meshes for one-dimensional elliptic and parabolic differential equations with the Robin boundary conditions.An explicit scheme and an implicit scheme are obtained by discretizing the equivalent integral form of the equation.For the explicit scheme with nodal values,the algebraic system can be solved by the Thomas method.For the implicit scheme with both nodal values and their derivatives,the system can be implemented by a prediction-correction procedure,where in the correction stage,an implicit formula for recovering the nodal derivatives is introduced.Taking two point boundary value problem as an example,we prove that both the explicit and implicit schemes are convergent with fourth order accuracy with respect to some standard discrete norms using the energy method.Two numerical examples demonstrate the correctness and effectiveness of the schemes,as well as the indispensability of using non-uniform meshes.
基金the National Natural Science Foundation of China (No. 50375135)the Talent Foundation of Beijing Jiaotong University (No. 2003RC059)
文摘A new model,called object model,for the simulation of cold roll-forming of tubes is presented.The model inherits the advantages of old models and is the embodiment of forming process that the strip is rolled step by step from feed rollers to last rolling pass.The elastic-plastic large deformation spline finite strip method based on updated Lagrangian method has been developed by improving the stiffness and transition matrix.Combined theory formulas and new analytical model,the forming process of a tube has been simulated successfully as an example.The analytical results are submitted and indicate that the proposed simulation method and new model are applicable.
基金heprojectissupportedbyNNSFofChina (No .1 9972 0 39) .
文摘In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.
基金supported by the National Natural Science Foundation of China (51109029,51178081,51138001,and 51009020)the State Key Development Program for Basic Research of China (2013CB035905)
文摘A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D FEM of 7 nodes and 10 nodes are constructed based on the basic formula. Using these proposed elements, the multiscale numerical model for foundation subjected to harmonic periodic load, the foundation model excited by external and internal dynamic load are studied. The results show the pro- posed finite elements have higher precision than the tradi- tional elements with 4 nodes. The proposed finite elements can describe the propagation of stress waves well whenever the foundation model excited by extemal or intemal dynamic load. The proposed finite elements can be also used to con- nect the multi-scale elements. And the proposed finite elements also have high precision to make multi-scale analysis for structure.
基金supported by the National Natural Science Foundation of China (60533060, 10672032, 10726067)Science Foundation of Dalian University of Technology (SFDUT07001)
文摘In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant. Polygonal elements can do well in simulation of the materials behavior and provide greater flexibility for the meshing of complex geometries. Hence, the study on the polygonal element is a very useful and necessary part in the finite element method. In this paper, an n-sided polygonal element based on quadratic spline interpolant, denoted by PS2 element, is presented using the triangular area coordinates and the B-net method. The PS2 element is conforming and can exactly model the quadratic field. It is valid for both convex and non-convex polygonal element, and insensitive to mesh distortions. In addition, no mapping or coordinate transformation is required and thus no Jacobian matrix and its inverse are evaluated. Some appropriate examples are employed to evaluate the performance of the proposed element.
基金supported by the National Natural Science Foundation of China(11001037,11102037 and 11290143)the Fundamental Research Funds for the Central Universities
文摘Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.
基金supported by the National Natural Science Foundation of China(11001037,11102037,11290143)the Fundamental Research Funds for the Central Universities(DUT13LK07)
文摘The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previous work, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the B- net method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian co- ordinates. In this paper, a thin plate spline element is devel- oped based on the spline element L8 and the refined tech- nique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.
基金Key Project(2004BA901A02) supported by the Ministry of Science and Technology of China
文摘The near crack line analysis method was used to investigate a centric crack loaded by two pairs of point shear forces in a finite plate, and the analytical solution was obtained. The solution includes the unit normal vector of the elastic-plastic boundary near the crack line, the elastic-plastic stress fields near the crack line, the variations of the length of the plastic zone along the crack line with an external load, and the bearing capacity of a finite plate with a centric crack loaded by two pairs of point shear forces. The results are sufficiently precise near the crack line because the assumptions of small scale yielding theory have not been made and no other assumptions are taken.
基金supported by the National Natura Science Foundation of China (Grant 51275424)973 Program (Gran2011CB610304)+1 种基金Research Fund for the Doctoral Program of Higher Education of China (Grant 20126102130003)the opening project (Grant KFJJ13-6M) of the State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology)
文摘In this work, a design procedure extending the B-spline based finite cell method into shape optimization is developed for axisymmetric solids involving the centrifugal force effect. We first replace the traditional conforming mesh in the finite element method with structured cells that are fixed during the whole design process with a view to avoid the sophisticated re-meshing and eventual mesh distortion.Then, B-spline shape functions are further implemented to yield a high-order continuity field along the cell boundary in stress analysis. By means of the implicit description of the shape boundary, stress sensitivity is analytically derived with respect to shape design variables. Finally, we illustrate the efficiency and accuracy of the proposed protocol by several numerical test cases as well as a whole design procedure carried out on an aeronautic turbine disk.
文摘In this paper,the spline finite element method is developed to investigate free vibration problems of beams.The cubic B-spline functions are used to construct the displacement field.The assembly of elements and the introduction of boundary conditions follow the standard finite element procedure.The results under various boundary conditions are compared with those obtained by the exact method and the finite difference method.It shows that the results are in excellent agreement with the analytical results and much more accurate than the results obtained by the finite difference method,especially for higher order modes.
基金The Project was supported by National Natural Science Foundation of China
文摘In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.
基金supported by the National Natural Science Foundation of China (Grant No.10771133)the Shanghai Leading Academic Discipline Project (Grant No.S30101)the Research Foundation for the Doctoral Program of Higher Education (Grant No.200802800010)
文摘The choice of self-concordant functions is the key to efficient algorithms for linear and quadratic convex optimizations, which provide a method with polynomial-time iterations to solve linear and quadratic convex optimization problems. The parameters of a self-concordant barrier function can be used to compute the complexity bound of the proposed algorithm. In this paper, it is proved that the finite barrier function is a local self-concordant barrier function. By deriving the local values of parameters of this barrier function, the desired complexity bound of an interior-point algorithm based on this local self-concordant function for linear optimization problem is obtained. The bound matches the best known bound for small-update methods.
基金Project supported by the National Natural Science Foundation of China (Nos. 10571046, 10371038)
文摘For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan＇s element energy projection （EEP） method has the accuracy O（h^min{2k,k＋4}） The theoretical analysis coincides the reported numerical results.
基金The project was supported by the National Natural Science Foundation of China(11001037,11102037,11072156)the Fundamental Research Funds for the Central Universities of China(DUT10ZD112,DUT10JS02)
文摘In this paper,a 13-node pyramid spline element is derived by using the tetrahedron volume coordinates and the B-net method,which achieves the second order completeness in Cartesian coordinates.Some appropriate examples were employed to evaluate the performance of the proposed element.The numerical results show that the spline element has much better performance compared with the isoparametric serendipity element Q20 and its degenerate pyramid element P13 especially when mesh is distorted,and it is comparable to the Lagrange element Q27.It has been demonstrated that the spline finite element method is an efficient tool for developing high accuracy elements.
基金the National Natural Science Foundation of China(Nos.11502098 and 11772158)the Natural Science Foundation of the Higher Education Institutions of Jiangsu Province(No.15KJB130003)+1 种基金the Doctoral Scientific Research Foundation of Jiangsu University of Science and Technology(No.120140003)the Fundamental Research Funds for Central Universities of China(No.30917011103)。
文摘Tlie rigid-flexible coupling dynamic modeling theory and the discretization methods of a rotating flexible rectangular thin plate are investigated in this paper.Based on the continuum mechanics,the rigid-flexible coupling dynamic model is established for the flexible rectangular thin plate undergoing large overall rotation,and the coupling term of the deformation which is caused by transverse deformation is considered.Assumed mode method(AMM),spline finite point method(SFPM)and Beizer finite point method(BFPM)are used to describe the deformation of the flexible rectangular plate,and then the dynamic equations of a rotating flexible rectangular thin plate undergoing overall motion are derived by Lagrange^equation of the second kind.The dynamics of a cantilever plate undergoing large overall rotation is simulated via using different dynamic models,and the simulation results of the first order approximation model are compared with those of the traditional zero-order approximation model.It is shown that the first order approximation model with the dynamic stiffening terms can correctly describe the dynamic behavior of the system undergoing large overall rotation,while the zero-order approximation model cannot get the correct results.And AMM.SFPM.BFPM can well describe the deformation of a rotating flexible rectangular plate.
