The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which...The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which can be solved by traditional numerical methods and an infinite domain with a crack which can be solved by Muskhelishvili’s solutions.However,this alternating method cannot be directly applied to the edge crack problems since partial crack surface of Muskhelishvili’s solutions is located outside the computational domain.In this paper,an improved alternating method,the spline fictitious boundary element alternating method(SFBEAM),based on infinite domain with the combination of spline fictitious boundary element method(SFBEM)and Muskhelishvili’s solutions is proposed to solve the edge crack problems.Since the SFBEM and Muskhelishvili’s solutions are obtained in the framework of infinite domain,no special treatment is needed for solving the problem of edge cracks.Different mixed boundary conditions edge crack problems with varies of computational parameters are given to certify the high precision,efficiency and applicability of the proposed method compared with other alternating methods and extend finite element method.展开更多
Non-singular fictitious boundary integral equations for orthotropic elastic plane problems were deduced according to boundary conditions by the techniques of singular-points-outside-domain. Then the unknown fictitious...Non-singular fictitious boundary integral equations for orthotropic elastic plane problems were deduced according to boundary conditions by the techniques of singular-points-outside-domain. Then the unknown fictitious load functions along the fictitious boundary were expressed in terms of basic spline functions, and the boundary-segment-least-squares method was proposed to eliminate the boundary residues obtained. By the above steps, numerical solutions to the integral equations can be achieved. Numerical examples are given to show the accuracy and efficiency of the proposed method.展开更多
The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a sp...The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.展开更多
Through a higher-order boundary element method based on NURBS （Non-uniform Rational B-splines）, the calculation of second-order low-frequency forces and slow drift motions is conducted for floating bodies. In the fl...Through a higher-order boundary element method based on NURBS （Non-uniform Rational B-splines）, the calculation of second-order low-frequency forces and slow drift motions is conducted for floating bodies. In the floating body＇s inner domain, an auxiliary equation is obtained by applying a Green function which satisfies the solid surface condition. Then, the auxiliary equation and the velocity potential equation are combined in the fluid domain to remove the solid angle coefficient and the singularity of the double layer potentials in the integral equation. Thus, a new velocity potential integral equation is obtained. The new equation is extended to the inner domain to reheve the irregular frequency effects; on the basis of the order analysis, the comparison is made about the contribution of all integral terms with the result in the second-order tow-frequency problem; the higher-order boundary element method based on NURBS is apphed to calculate the geometric position and velocity potentials; the slow drift motions are calculated by the spectrum analysis method. Removing the solid angle coefficient can apply NURBS technology to the hydrodynamic calculation of floating bodies with complex surfaces, and the extended boundary integral method can reduce the irregular frequency effects. Order analysis shows that free surface integral can be neglected, and the numerical results can also prove the correctness of order analysis. The results of second-order low-frequency forces and slow drift motions and the comparison with the results from references show that the application of the NURBS technology to the second-order low-frequency problem is of high efficiency and credible results.展开更多
For the computation of wave forces on structures, a B-spline expansion is applied to discretize the body surface, and represent the velocity potentials on the body surface. The expansion coefficients for the body geom...For the computation of wave forces on structures, a B-spline expansion is applied to discretize the body surface, and represent the velocity potentials on the body surface. The expansion coefficients for the body geometry are determined by the Least Square Method, and the coefficients for velocity potentials by the Galerkin method. The method can give continuous description of velocity potentials and their derivatives on the whole smooth body surface. The method has been implemented, and numerical results show that the method gives very accurate results and its convergence is fast.展开更多
In this article,we discuss a least-squares/fictitious domain method for the solution of linear elliptic boundary value problems with Robin boundary conditions.LetΩandωbe two bounded domains of R d such thatω⊂Ω.For a...In this article,we discuss a least-squares/fictitious domain method for the solution of linear elliptic boundary value problems with Robin boundary conditions.LetΩandωbe two bounded domains of R d such thatω⊂Ω.For a linear elliptic problem inΩ\ωwith Robin boundary condition on the boundaryγofω,our goal here is to develop a fictitious domain method where one solves a variant of the original problem on the fullΩ,followed by a well-chosen correction overω.This method is of the virtual control type and relies on a least-squares formulation making the problem solvable by a conjugate gradient algorithm operating in a well chosen control space.Numerical results obtained when applying our method to the solution of two-dimensional elliptic and parabolic problems are given;they suggest optimal order of convergence.展开更多
In this paper, a boundary element method is first applied to real-time animation of deformable objects and to simplify data preparation. Next, the visible external surface of the object in deforming process is represe...In this paper, a boundary element method is first applied to real-time animation of deformable objects and to simplify data preparation. Next, the visible external surface of the object in deforming process is represented by B-spline surface, whose control points are embedded in dynamic equations of BEM. Finally, the above method is applied to anatomical simulation. A pituitary model in human brain, which is reconstructed from a set of anatomical sections, is selected to be the deformable object under action of virtual tool such as scalpel or probe. It produces fair graphic realism and high speed performance. The results show that BEM not only has less computational expense than FEM, but also is convenient to combine with the 3D reconstruction and surface modeling as it enables the reduction of the dimensionality of the problem by one.展开更多
基金supported by the National Natural Science Foundation of China(51078150)the National Natural Science Foundation of China(11602087)+1 种基金the State Key Laboratory of Subtropical Building Science,South China University of Technology(2017ZB32)National Undergraduate Innovative and Entrepreneurial Training Program(201810561180).
文摘The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which can be solved by traditional numerical methods and an infinite domain with a crack which can be solved by Muskhelishvili’s solutions.However,this alternating method cannot be directly applied to the edge crack problems since partial crack surface of Muskhelishvili’s solutions is located outside the computational domain.In this paper,an improved alternating method,the spline fictitious boundary element alternating method(SFBEAM),based on infinite domain with the combination of spline fictitious boundary element method(SFBEM)and Muskhelishvili’s solutions is proposed to solve the edge crack problems.Since the SFBEM and Muskhelishvili’s solutions are obtained in the framework of infinite domain,no special treatment is needed for solving the problem of edge cracks.Different mixed boundary conditions edge crack problems with varies of computational parameters are given to certify the high precision,efficiency and applicability of the proposed method compared with other alternating methods and extend finite element method.
文摘Non-singular fictitious boundary integral equations for orthotropic elastic plane problems were deduced according to boundary conditions by the techniques of singular-points-outside-domain. Then the unknown fictitious load functions along the fictitious boundary were expressed in terms of basic spline functions, and the boundary-segment-least-squares method was proposed to eliminate the boundary residues obtained. By the above steps, numerical solutions to the integral equations can be achieved. Numerical examples are given to show the accuracy and efficiency of the proposed method.
文摘The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.
文摘Through a higher-order boundary element method based on NURBS （Non-uniform Rational B-splines）, the calculation of second-order low-frequency forces and slow drift motions is conducted for floating bodies. In the floating body＇s inner domain, an auxiliary equation is obtained by applying a Green function which satisfies the solid surface condition. Then, the auxiliary equation and the velocity potential equation are combined in the fluid domain to remove the solid angle coefficient and the singularity of the double layer potentials in the integral equation. Thus, a new velocity potential integral equation is obtained. The new equation is extended to the inner domain to reheve the irregular frequency effects; on the basis of the order analysis, the comparison is made about the contribution of all integral terms with the result in the second-order tow-frequency problem; the higher-order boundary element method based on NURBS is apphed to calculate the geometric position and velocity potentials; the slow drift motions are calculated by the spectrum analysis method. Removing the solid angle coefficient can apply NURBS technology to the hydrodynamic calculation of floating bodies with complex surfaces, and the extended boundary integral method can reduce the irregular frequency effects. Order analysis shows that free surface integral can be neglected, and the numerical results can also prove the correctness of order analysis. The results of second-order low-frequency forces and slow drift motions and the comparison with the results from references show that the application of the NURBS technology to the second-order low-frequency problem is of high efficiency and credible results.
基金National Natural Science Foundation of China under the Grant No.19732004
文摘For the computation of wave forces on structures, a B-spline expansion is applied to discretize the body surface, and represent the velocity potentials on the body surface. The expansion coefficients for the body geometry are determined by the Least Square Method, and the coefficients for velocity potentials by the Galerkin method. The method can give continuous description of velocity potentials and their derivatives on the whole smooth body surface. The method has been implemented, and numerical results show that the method gives very accurate results and its convergence is fast.
基金The first author acknowledge the support of the Institute for Advanced Study(IAS)at The Hong Kong University of Science and TechnologyThe work is partially supported by grants from RGC CA05/06.SC01 and RGC-CERG 603107.
文摘In this article,we discuss a least-squares/fictitious domain method for the solution of linear elliptic boundary value problems with Robin boundary conditions.LetΩandωbe two bounded domains of R d such thatω⊂Ω.For a linear elliptic problem inΩ\ωwith Robin boundary condition on the boundaryγofω,our goal here is to develop a fictitious domain method where one solves a variant of the original problem on the fullΩ,followed by a well-chosen correction overω.This method is of the virtual control type and relies on a least-squares formulation making the problem solvable by a conjugate gradient algorithm operating in a well chosen control space.Numerical results obtained when applying our method to the solution of two-dimensional elliptic and parabolic problems are given;they suggest optimal order of convergence.
文摘In this paper, a boundary element method is first applied to real-time animation of deformable objects and to simplify data preparation. Next, the visible external surface of the object in deforming process is represented by B-spline surface, whose control points are embedded in dynamic equations of BEM. Finally, the above method is applied to anatomical simulation. A pituitary model in human brain, which is reconstructed from a set of anatomical sections, is selected to be the deformable object under action of virtual tool such as scalpel or probe. It produces fair graphic realism and high speed performance. The results show that BEM not only has less computational expense than FEM, but also is convenient to combine with the 3D reconstruction and surface modeling as it enables the reduction of the dimensionality of the problem by one.
