The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an ...The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an orthogonal function system mixed a weight function is defined. Next the improved moving least-square approximation is discussed in detail. The improved method has higher computational efficiency and precision than the old method, and cannot form an ill-conditioned equation system. A boundary element-free method (BEFM) for elastodynamics problems is presented by combining the boundary integral equation method for elastodynamics and the improved moving least-square approximation. The boundary element-free method is a meshless method of boundary integral equation and is a direct numerical method compared with others, in which the basic unknowns are the real solutions of the nodal variables and the boundary conditions can be applied easily. The boundary element-free method has a higher computational efficiency and precision. In addition, the numerical procedure of the boundary element-free method for elastodynamics problems is presented in this paper. Finally, some numerical examples are given.展开更多
A method based on Holmholtz equation is presented to predict thc acoustic radiation of a vibrating box by combining the finite clement method (FEM) with boundary element method (BEM). The prediction formulation is val...A method based on Holmholtz equation is presented to predict thc acoustic radiation of a vibrating box by combining the finite clement method (FEM) with boundary element method (BEM). The prediction formulation is valid for all frequencies if an extra equation is added. The acoustic field of a vibrating box has been calculated and measured. The effectiveness of the method has been proved by test results.展开更多
This paper presents mechanical quadrature methods for solving first-kind boundary integral equations on polygonal regions, which possesses high accuracy O(h0^3)and low computing complexities. Moreover, the multivariat...This paper presents mechanical quadrature methods for solving first-kind boundary integral equations on polygonal regions, which possesses high accuracy O(h0^3)and low computing complexities. Moreover, the multivariate asymptotic expansion of the error with hi^3(i = 1,…,d) power is shown. Using the multi-parameter asymptotic expansion, we not only get a high precisioin approximation solution by means of the splitting extrapolation, but also derive a posteriori estimation.展开更多
The paper begins by discussing the interpolating moving least-squares (IMLS) method. Then the formulae of the IMLS method obtained by Lancaster are revised. On the basis of the boundary element-free method (BEFM), com...The paper begins by discussing the interpolating moving least-squares (IMLS) method. Then the formulae of the IMLS method obtained by Lancaster are revised. On the basis of the boundary element-free method (BEFM), combining the boundary integral equation method with the IMLS method improved in this paper, the interpolating boundary element-free method (IBEFM) for two-dimensional elasticity problems is presented, and the corresponding formulae of the IBEFM for two-dimensional elasticity problems are obtained. In the IMLS method in this paper, the shape function satisfies the property of Kronecker δ function, and then in the IBEFM the boundary conditions can be applied directly and easily. The IBEFM is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution to the nodal variables. Thus it gives a greater computational precision. Numerical examples are presented to demonstrate the method.展开更多
In order to improve rolled strip quality, precise plate shape control theory should be established. Roll flat- tening theory is an important part of the plate shape theory. To improve the accuracy of roll flattening c...In order to improve rolled strip quality, precise plate shape control theory should be established. Roll flat- tening theory is an important part of the plate shape theory. To improve the accuracy of roll flattening calculation based on semi infinite body model, especially near the two roll barrel edges, a new and more accurate roll flattening model is proposed. Based on boundary integral equation method, an analytical model for solving a finite length semi infinite body is established. The lateral surface displacement field of the finite length semi-infinite body is simulated by finite element method (FEM) and lateral surface displacement decay functions are established. Based on the boundary integral equation method, the numerical solution of the finite length semi-infinite body under the distribu ted force is obtained and an accurate roll flattening model is established. Different from the traditional semi-infinite body model, the matrix form of the new roll flattening model is established through the mathematical derivation. The result from the new model is more consistent with that by FEM especially near the edges.展开更多
For the discretization of higher order elements, this paper presents a modified integral domain method to remove the irregular frequencies inherited in the integral equation of wave diffraction and radiation from a su...For the discretization of higher order elements, this paper presents a modified integral domain method to remove the irregular frequencies inherited in the integral equation of wave diffraction and radiation from a surface-piercing body. The set of over-determined linear equations obtained from the method is modified into a normal set of linear equations by superposing a set of linear equations with zero solutions. Numerical experiments have also been carried out to find the optimum choice of the size of the auxiliary domain and the discretization on it.展开更多
提出基于边界面法(Boundary Face Method,BFM)的完整实体应力分析方法.在该分析中,避免对结构作几何上的简化,结构的所有局部细节都按实际形状尺寸作为三维实体处理.以边界积分方程为理论基础的BFM是完整实体应力分析的自然选择.在该方...提出基于边界面法(Boundary Face Method,BFM)的完整实体应力分析方法.在该分析中,避免对结构作几何上的简化,结构的所有局部细节都按实际形状尺寸作为三维实体处理.以边界积分方程为理论基础的BFM是完整实体应力分析的自然选择.在该方法中,边界积分和场变量插值都在实体边界曲面的参数空间里实现.高斯积分点的几何数据,如坐标、雅可比和外法向量都直接由曲面算得,而不是通过单元插值近似获得,从而避免几何误差.该方法的实现直接基于边界表征的CAD模型,可做到与CAD软件的无缝连接.线弹性问题的应用实例表明,该方法可以简单有效地模拟具有细小特征的复杂结构,并且计算结果的应力精度比边界元法(Boundary Element Method,BEM)和有限元法(Finite Element Method,FEM)高.展开更多
A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to...A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to the boundary integral equation method (BIEM) for current diffraction from arbitrary 3D bodies. The boundary integral equation is discretized by higher order elements, the FMM is applied to avoid the matrix/vector product, and the resulting algebraic equation is solved by the Generalized Conjugate Residual method (GCR). Numerical examination shows that the FMM is more efficient than the direct evaluation method in computational cost and storage of computers.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.10571118)the Shanghai Leading Academic Discipline Project(Grant No.Y0103).
文摘The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an orthogonal function system mixed a weight function is defined. Next the improved moving least-square approximation is discussed in detail. The improved method has higher computational efficiency and precision than the old method, and cannot form an ill-conditioned equation system. A boundary element-free method (BEFM) for elastodynamics problems is presented by combining the boundary integral equation method for elastodynamics and the improved moving least-square approximation. The boundary element-free method is a meshless method of boundary integral equation and is a direct numerical method compared with others, in which the basic unknowns are the real solutions of the nodal variables and the boundary conditions can be applied easily. The boundary element-free method has a higher computational efficiency and precision. In addition, the numerical procedure of the boundary element-free method for elastodynamics problems is presented in this paper. Finally, some numerical examples are given.
文摘A method based on Holmholtz equation is presented to predict thc acoustic radiation of a vibrating box by combining the finite clement method (FEM) with boundary element method (BEM). The prediction formulation is valid for all frequencies if an extra equation is added. The acoustic field of a vibrating box has been calculated and measured. The effectiveness of the method has been proved by test results.
文摘This paper presents mechanical quadrature methods for solving first-kind boundary integral equations on polygonal regions, which possesses high accuracy O(h0^3)and low computing complexities. Moreover, the multivariate asymptotic expansion of the error with hi^3(i = 1,…,d) power is shown. Using the multi-parameter asymptotic expansion, we not only get a high precisioin approximation solution by means of the splitting extrapolation, but also derive a posteriori estimation.
基金supported by the National Natural Science Foundation of China (Grant No. 10871124)the Innovation Program of Shanghai Municipal Education Commission (Grant No. 09ZZ99)the ShanghaiLeading Academic Discipline Project (Grant No. J50103)
文摘The paper begins by discussing the interpolating moving least-squares (IMLS) method. Then the formulae of the IMLS method obtained by Lancaster are revised. On the basis of the boundary element-free method (BEFM), combining the boundary integral equation method with the IMLS method improved in this paper, the interpolating boundary element-free method (IBEFM) for two-dimensional elasticity problems is presented, and the corresponding formulae of the IBEFM for two-dimensional elasticity problems are obtained. In the IMLS method in this paper, the shape function satisfies the property of Kronecker δ function, and then in the IBEFM the boundary conditions can be applied directly and easily. The IBEFM is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution to the nodal variables. Thus it gives a greater computational precision. Numerical examples are presented to demonstrate the method.
基金Item Sponsored by National Natural Science Foundation of China(51075353)
文摘In order to improve rolled strip quality, precise plate shape control theory should be established. Roll flat- tening theory is an important part of the plate shape theory. To improve the accuracy of roll flattening calculation based on semi infinite body model, especially near the two roll barrel edges, a new and more accurate roll flattening model is proposed. Based on boundary integral equation method, an analytical model for solving a finite length semi infinite body is established. The lateral surface displacement field of the finite length semi-infinite body is simulated by finite element method (FEM) and lateral surface displacement decay functions are established. Based on the boundary integral equation method, the numerical solution of the finite length semi-infinite body under the distribu ted force is obtained and an accurate roll flattening model is established. Different from the traditional semi-infinite body model, the matrix form of the new roll flattening model is established through the mathematical derivation. The result from the new model is more consistent with that by FEM especially near the edges.
基金This work is a part of the research project financially supported by the National Natural Science Foundation of China
文摘For the discretization of higher order elements, this paper presents a modified integral domain method to remove the irregular frequencies inherited in the integral equation of wave diffraction and radiation from a surface-piercing body. The set of over-determined linear equations obtained from the method is modified into a normal set of linear equations by superposing a set of linear equations with zero solutions. Numerical experiments have also been carried out to find the optimum choice of the size of the auxiliary domain and the discretization on it.
文摘提出基于边界面法(Boundary Face Method,BFM)的完整实体应力分析方法.在该分析中,避免对结构作几何上的简化,结构的所有局部细节都按实际形状尺寸作为三维实体处理.以边界积分方程为理论基础的BFM是完整实体应力分析的自然选择.在该方法中,边界积分和场变量插值都在实体边界曲面的参数空间里实现.高斯积分点的几何数据,如坐标、雅可比和外法向量都直接由曲面算得,而不是通过单元插值近似获得,从而避免几何误差.该方法的实现直接基于边界表征的CAD模型,可做到与CAD软件的无缝连接.线弹性问题的应用实例表明,该方法可以简单有效地模拟具有细小特征的复杂结构,并且计算结果的应力精度比边界元法(Boundary Element Method,BEM)和有限元法(Finite Element Method,FEM)高.
文摘A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to the boundary integral equation method (BIEM) for current diffraction from arbitrary 3D bodies. The boundary integral equation is discretized by higher order elements, the FMM is applied to avoid the matrix/vector product, and the resulting algebraic equation is solved by the Generalized Conjugate Residual method (GCR). Numerical examination shows that the FMM is more efficient than the direct evaluation method in computational cost and storage of computers.