This paper focuses on the direct and inverse problems for a third-order self-adjoint differential operator with non-local potential and anti-periodic boundary conditions.Firstly,we obtain the expressions for the chara...This paper focuses on the direct and inverse problems for a third-order self-adjoint differential operator with non-local potential and anti-periodic boundary conditions.Firstly,we obtain the expressions for the characteristic function and resolvent of this third-order differential operator.Secondly,by using the expression for the resolvent of the operator,we prove that the spectrum for this operator consists of simple eigenvalues and a finite number of eigenvalues with multiplicity 2.Finally,we solve the inverse problem for this operator,which states that the non-local potential function can be reconstructed from four spectra.Specially,we prove the Ambarzumyan theorem and indicate that odd or even potential functions can be reconstructed by three spectra.展开更多
An analytical scheme, which avoids using the standard Gaussian approximate quadrature to treat the boundary integrals in direct boundary element method (DBEM) of two-dimensional potential and elastic problems, is esta...An analytical scheme, which avoids using the standard Gaussian approximate quadrature to treat the boundary integrals in direct boundary element method (DBEM) of two-dimensional potential and elastic problems, is established. With some numerical results, it is shown that the better precision and high computational efficiency, especially in the band of the domain near boundary, can be derived by the present scheme.展开更多
An adaptive cell-based domain integration method（CDIM） is proposed for the treatment of domain integrals in 3D boundary element method（BEM）. The domain integrals are computed in background cells rather than volume...An adaptive cell-based domain integration method（CDIM） is proposed for the treatment of domain integrals in 3D boundary element method（BEM）. The domain integrals are computed in background cells rather than volume elements. The cells are created from the boundary elements based on an adaptive oct-tree structure and no other discretization is needed. Cells containing the boundary elements are subdivided into smaller sub-cells adaptively according to the sizes and levels of the boundary elements; and the sub-cells outside the domain are deleted to obtain the desired accuracy. The method is applied in the 3D potential and elasticity problems in this paper.展开更多
The optical potential ambiguity is a long-standing problem in the analysis of elastic scattering data.For a specific collid-ing system,ambiguous potential families can lead to different behaviors in the nearside and f...The optical potential ambiguity is a long-standing problem in the analysis of elastic scattering data.For a specific collid-ing system,ambiguous potential families can lead to different behaviors in the nearside and farside scattering components.By contrast,the envelope method can decompose the experimental data into two components with negative and positive deflection angles,respectively.Hence,a question arises as to whether the comparison between the calculated nearside(or farside)component and the derived positive-deflection-angle(or negative-deflection-angle)component can help analyze the potential ambiguity problem.In this study,we conducted a trial application of the envelope method to the potential ambiguity problem.The envelope method was improved by including uncertainties in the experimental data.The colliding systems of 16O+28Si at 215.2 MeV and 12C+12C at 1016 MeV were considered in the analyses.For each colliding system,the angular distribution experimental data were described nearly equally well by two potential sets,one of which is“surface transpar-ent”and the other is refractive.The calculated angular distributions were decomposed into nearside and farside scattering components.Using the improved envelope method,the experimental data were decomposed into the positive-deflection-angle and negative-deflection-angle components,which were then compared with the calculated nearside and farside components.The capability of the envelope method to analyze the potential ambiguities was also discussed.展开更多
Without applying any stable element techniques in the mixed methods, two simple generalized mixed element(GME) formulations were derived by combining the minimum potential energy principle and Hellinger–Reissner(H–R...Without applying any stable element techniques in the mixed methods, two simple generalized mixed element(GME) formulations were derived by combining the minimum potential energy principle and Hellinger–Reissner(H–R) variational principle. The main features of the GME formulations are that the common C0-continuous polynomial shape functions for displacement methods are used to express both displacement and stress variables, and the coefficient matrix of these formulations is not only automatically symmetric but also invertible. Hence, the numerical results of the generalized mixed methods based on the GME formulations are stable. Displacement as well as stress results can be obtained directly from the algebraic system for finite element analysis after introducing stress and displacement boundary conditions simultaneously. Numerical examples show that displacement and stress results retain the same accuracy. The results of the noncompatible generalized mixed method proposed herein are more accurate than those of the standard noncompatible displacement method. The noncompatible generalized mixed element is less sensitive to element geometric distortions.展开更多
In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the II...In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.展开更多
In the present paper, the spectrums of off-diagonal infinite-dimensional Hamiltonian operators are studied. At first, we prove that the spectrum, the continuous-spectrum, and the union of the point-spectrum and residu...In the present paper, the spectrums of off-diagonal infinite-dimensional Hamiltonian operators are studied. At first, we prove that the spectrum, the continuous-spectrum, and the union of the point-spectrum and residual- spectrum of the operators are symmetric with respect to real axis and imaginary axis. Then for the purpose of reducing the dimension of the studied problems, the spectrums of the operators are expressed by the spectrums of the product of two self-adjoint operators in state spac,3. At last, the above-mentioned results are applied to plane elasticity problems, which shows the practicability of the results.展开更多
The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （B...The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （BEFM）, combining the boundary integral equation （BIE） method with the IMLS method, the improved boundary element-free method （IBEFM） for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.展开更多
In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This st...In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.展开更多
In this paper, based on the element-free Galerkin （EFG） method and the improved complex variable moving least- square （ICVMLS） approximation, a new meshless method, which is the improved complex variable element-f...In this paper, based on the element-free Galerkin （EFG） method and the improved complex variable moving least- square （ICVMLS） approximation, a new meshless method, which is the improved complex variable element-free Galerkin （ICVEFG） method for two-dimensional potential problems, is presented. In the method, the integral weak form of control equations is employed, and the Lagrange multiplier is used to apply the essential boundary conditions. Then the corresponding formulas of the ICVEFG method for two-dimensional potential problems are obtained. Compared with the complex variable moving least-square （CVMLS） approximation proposed by Cheng, the functional in the ICVMLS approximation has an explicit physical meaning. Furthermore, the ICVEFG method has greater computational precision and efficiency. Three numerical examples are given to show the validity of the proposed method.展开更多
This study is concerned with a new,explicit approach by means of which forms of the large strain elastic potential for multiaxial rubberlike elasticity may be obtained based on data for a single deformation mode.As a ...This study is concerned with a new,explicit approach by means of which forms of the large strain elastic potential for multiaxial rubberlike elasticity may be obtained based on data for a single deformation mode.As a departure from usual studies,here for the first time errors may be estimated and rendered minimal for all possible deformation modes and,furthermore,failure behavior may be incorporated.Numerical examples presented are in accurate agreement with Treloar＇s well-known data.展开更多
Weak solution （or generalized solution） for the boundary-value problems of partial differential equations of elasticity of 3D （three-dimensional） quasicrystals is given, in which the matrix expression is used. In ...Weak solution （or generalized solution） for the boundary-value problems of partial differential equations of elasticity of 3D （three-dimensional） quasicrystals is given, in which the matrix expression is used. In terms of Korn inequality and theory of function space, we prove the uniqueness of the weak solution. This gives an extension of existence theorem of solution for classical elasticity to that of quasicrystals, and develops the weak solution theory of elasticity of 2D quasicrystals given by the second author of the paper and his students.展开更多
This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-d...This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.展开更多
By virtue of the rational interpolation procedure and logarithmic strain, a direct approach is proposed to obtain elastic potentials that exactly match uniaxial data and shear data for elastomers. This approach reduce...By virtue of the rational interpolation procedure and logarithmic strain, a direct approach is proposed to obtain elastic potentials that exactly match uniaxial data and shear data for elastomers. This approach reduces the determination of multi axial elastic potentials to that of two one-dimensional potentials, thus bypassing usual cumbersome procedures of identifying a number of unknown parameters. Predictions of the suggested potential are derived for a general biaxial stretch test and compared with the classical data given by Rivlin and Saunders （Rivlin, R. S. and Saunders, D. W. Large elastic deformation of isotropic materials. VII： experiments on the deformation of rubber. Phill. Trans. Royal Soc. London A, 243, 251-288 （1951））. Good agreement is achieved with these extensive data.展开更多
In this paper, a locking-free nonconforming rectangular finite element scheme is presented for the planar elasticity problem with pure displacement boundary condition. Meanwhile, we prove that this element is also con...In this paper, a locking-free nonconforming rectangular finite element scheme is presented for the planar elasticity problem with pure displacement boundary condition. Meanwhile, we prove that this element is also convergent for stationary Stokes problem.展开更多
The universal practices have been centralizing on the research of regularization to the direct boundary integal equations （DBIEs）. The character is elimination of singularities by using the simple solutions. However...The universal practices have been centralizing on the research of regularization to the direct boundary integal equations （DBIEs）. The character is elimination of singularities by using the simple solutions. However, up to now the research of regularization to the first kind integral equations for plane potential problems has never been found in previous literatures. The presentation is mainly devoted to the research on the regularization of the singular boundary integral equations with indirect unknowns. A novel view and idea is presented herein, in which the regularized boundary integral equations with indirect unknowns without including the Cauchy principal value （CPV） and Hadamard-finite-part （HFP） integrals are established for the plane potential problems. With some numerical results, it is shown that the better accuracy and higher efficiency, especially on the boundary, can be achieved by the present system.展开更多
In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body...In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.展开更多
This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial diffe...This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial differential equations of the above 2D problems is rewritten as an upper triangular differential system. For the associated operator matrix, the existence and the completeness of two normed orthogonal eigenfunction systems in some space are obtained, which belong to the two block operators arising in the operator matrix. Moreover, the general solution to the above 2D problem is given by the eigenfunction expansion method.展开更多
Based on the complex variable moving least-square（CVMLS） approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin（CVMLPG） method of solving two-dimensional potential proble...Based on the complex variable moving least-square（CVMLS） approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin（CVMLPG） method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square（MLS） approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin（MLPG） method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method.展开更多
The symplectic approach proposed and developed by Zhong et al. in 1990s for elasticity problems is a rational analytical method, in which ample experience is not needed as in the conventional semi-inverse method. In t...The symplectic approach proposed and developed by Zhong et al. in 1990s for elasticity problems is a rational analytical method, in which ample experience is not needed as in the conventional semi-inverse method. In the symplectic space, elasticity problems can be solved using the method of separation of variables along with the eigenfunction expansion technique, as in traditional Fourier analysis. The eigensolutions include those corresponding to zero and nonzero eigenvalues. The latter group can be further divided into α-and β-sets. This paper reformulates the form of β-set eigensolutions to achieve the stability of numerical calculation, which is very important to obtain accurate results within the symplectic frame. An example is finally given and numerical results are compared and discussed.展开更多
基金supported by the Tianjin Municipal Science and Technology Program of China(No.23JCZDJC00070)。
文摘This paper focuses on the direct and inverse problems for a third-order self-adjoint differential operator with non-local potential and anti-periodic boundary conditions.Firstly,we obtain the expressions for the characteristic function and resolvent of this third-order differential operator.Secondly,by using the expression for the resolvent of the operator,we prove that the spectrum for this operator consists of simple eigenvalues and a finite number of eigenvalues with multiplicity 2.Finally,we solve the inverse problem for this operator,which states that the non-local potential function can be reconstructed from four spectra.Specially,we prove the Ambarzumyan theorem and indicate that odd or even potential functions can be reconstructed by three spectra.
文摘An analytical scheme, which avoids using the standard Gaussian approximate quadrature to treat the boundary integrals in direct boundary element method (DBEM) of two-dimensional potential and elastic problems, is established. With some numerical results, it is shown that the better precision and high computational efficiency, especially in the band of the domain near boundary, can be derived by the present scheme.
基金Financial support for the project from the National Natural Science Foundation of China(No.51609181)
文摘An adaptive cell-based domain integration method（CDIM） is proposed for the treatment of domain integrals in 3D boundary element method（BEM）. The domain integrals are computed in background cells rather than volume elements. The cells are created from the boundary elements based on an adaptive oct-tree structure and no other discretization is needed. Cells containing the boundary elements are subdivided into smaller sub-cells adaptively according to the sizes and levels of the boundary elements; and the sub-cells outside the domain are deleted to obtain the desired accuracy. The method is applied in the 3D potential and elasticity problems in this paper.
基金This work was supported by the National Natural Science Foundation of China(Nos.12005047 and U1832105).
文摘The optical potential ambiguity is a long-standing problem in the analysis of elastic scattering data.For a specific collid-ing system,ambiguous potential families can lead to different behaviors in the nearside and farside scattering components.By contrast,the envelope method can decompose the experimental data into two components with negative and positive deflection angles,respectively.Hence,a question arises as to whether the comparison between the calculated nearside(or farside)component and the derived positive-deflection-angle(or negative-deflection-angle)component can help analyze the potential ambiguity problem.In this study,we conducted a trial application of the envelope method to the potential ambiguity problem.The envelope method was improved by including uncertainties in the experimental data.The colliding systems of 16O+28Si at 215.2 MeV and 12C+12C at 1016 MeV were considered in the analyses.For each colliding system,the angular distribution experimental data were described nearly equally well by two potential sets,one of which is“surface transpar-ent”and the other is refractive.The calculated angular distributions were decomposed into nearside and farside scattering components.Using the improved envelope method,the experimental data were decomposed into the positive-deflection-angle and negative-deflection-angle components,which were then compared with the calculated nearside and farside components.The capability of the envelope method to analyze the potential ambiguities was also discussed.
基金supported by the National Natural Science Foundation of China (Grant 11502286)
文摘Without applying any stable element techniques in the mixed methods, two simple generalized mixed element(GME) formulations were derived by combining the minimum potential energy principle and Hellinger–Reissner(H–R) variational principle. The main features of the GME formulations are that the common C0-continuous polynomial shape functions for displacement methods are used to express both displacement and stress variables, and the coefficient matrix of these formulations is not only automatically symmetric but also invertible. Hence, the numerical results of the generalized mixed methods based on the GME formulations are stable. Displacement as well as stress results can be obtained directly from the algebraic system for finite element analysis after introducing stress and displacement boundary conditions simultaneously. Numerical examples show that displacement and stress results retain the same accuracy. The results of the noncompatible generalized mixed method proposed herein are more accurate than those of the standard noncompatible displacement method. The noncompatible generalized mixed element is less sensitive to element geometric distortions.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)
文摘In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.
基金supported by the National Natural Science Foundation of China under Grant No.10562002the Specialized Research Fund for the Doctoral Program of Higher Education of China under Grant No.20070126002+1 种基金the Natural Science Foundation of Inner Mongolia under Grant No.200508010103the Inner Mongolia University Scientific Research Starting Foundation for Talented Scholars under Grant No.207066
文摘In the present paper, the spectrums of off-diagonal infinite-dimensional Hamiltonian operators are studied. At first, we prove that the spectrum, the continuous-spectrum, and the union of the point-spectrum and residual- spectrum of the operators are symmetric with respect to real axis and imaginary axis. Then for the purpose of reducing the dimension of the studied problems, the spectrums of the operators are expressed by the spectrums of the product of two self-adjoint operators in state spac,3. At last, the above-mentioned results are applied to plane elasticity problems, which shows the practicability of the results.
基金Project supported by the National Natural Science Foundation of China (Grant No 10871124)Innovation Program of Shanghai Municipal Education Commission (Grant No 09ZZ99)Shanghai Leading Academic Discipline Project (Grant No J50103)
文摘The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （BEFM）, combining the boundary integral equation （BIE） method with the IMLS method, the improved boundary element-free method （IBEFM） for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China(Grant No.10902076)the Natural Science Foundation of Shanxi Province of China(Grant No.2007011009)+1 种基金the Scientific Research and Development Program of the Shanxi Higher Education Institutions(Grant No.20091131)the Doctoral Startup Foundation of Taiyuan University of Science and Technology(Grant No.200708)
文摘In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)the Innovation Fund Project for Graduate Student of Shanghai University,China (Grant No. SHUCX112359)
文摘In this paper, based on the element-free Galerkin （EFG） method and the improved complex variable moving least- square （ICVMLS） approximation, a new meshless method, which is the improved complex variable element-free Galerkin （ICVEFG） method for two-dimensional potential problems, is presented. In the method, the integral weak form of control equations is employed, and the Lagrange multiplier is used to apply the essential boundary conditions. Then the corresponding formulas of the ICVEFG method for two-dimensional potential problems are obtained. Compared with the complex variable moving least-square （CVMLS） approximation proposed by Cheng, the functional in the ICVMLS approximation has an explicit physical meaning. Furthermore, the ICVEFG method has greater computational precision and efficiency. Three numerical examples are given to show the validity of the proposed method.
基金the support of the start-up fund from the Education Committee of China through Shanghai University(Grant S.15-B002-09-032)the fund for research innovation from Shanghai University(Grants S.10-0401-12-001)the fund from Natural Science Foundation of China(Grants 11372172,11472164)
文摘This study is concerned with a new,explicit approach by means of which forms of the large strain elastic potential for multiaxial rubberlike elasticity may be obtained based on data for a single deformation mode.As a departure from usual studies,here for the first time errors may be estimated and rendered minimal for all possible deformation modes and,furthermore,failure behavior may be incorporated.Numerical examples presented are in accurate agreement with Treloar＇s well-known data.
基金Project supported by the National Natural Science Foundation of China (Nos.10372016 and 10672022)
文摘Weak solution （or generalized solution） for the boundary-value problems of partial differential equations of elasticity of 3D （three-dimensional） quasicrystals is given, in which the matrix expression is used. In terms of Korn inequality and theory of function space, we prove the uniqueness of the weak solution. This gives an extension of existence theorem of solution for classical elasticity to that of quasicrystals, and develops the weak solution theory of elasticity of 2D quasicrystals given by the second author of the paper and his students.
基金supported by the National Natural Science Foundation of China (Grants 11571223, 51404160)Shanxi Province Science Foundation for Youths (Grant 2014021025-1)
文摘This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.
基金Project supported by the National Natural Science Foundation of China(No.11372172)the 211-Plan of the Education Committee of China(No.A.15-B002-09-032)the Research Innovation Fund of Shanghai University(No.A.10-0401-12-001)
文摘By virtue of the rational interpolation procedure and logarithmic strain, a direct approach is proposed to obtain elastic potentials that exactly match uniaxial data and shear data for elastomers. This approach reduces the determination of multi axial elastic potentials to that of two one-dimensional potentials, thus bypassing usual cumbersome procedures of identifying a number of unknown parameters. Predictions of the suggested potential are derived for a general biaxial stretch test and compared with the classical data given by Rivlin and Saunders （Rivlin, R. S. and Saunders, D. W. Large elastic deformation of isotropic materials. VII： experiments on the deformation of rubber. Phill. Trans. Royal Soc. London A, 243, 251-288 （1951））. Good agreement is achieved with these extensive data.
文摘In this paper, a locking-free nonconforming rectangular finite element scheme is presented for the planar elasticity problem with pure displacement boundary condition. Meanwhile, we prove that this element is also convergent for stationary Stokes problem.
基金Project supported by the National Natural Science Foundation of China (No.10571110)the Natural Science Foundation of Shandong Province of China (No.2003ZX12)
文摘The universal practices have been centralizing on the research of regularization to the direct boundary integal equations （DBIEs）. The character is elimination of singularities by using the simple solutions. However, up to now the research of regularization to the first kind integral equations for plane potential problems has never been found in previous literatures. The presentation is mainly devoted to the research on the regularization of the singular boundary integral equations with indirect unknowns. A novel view and idea is presented herein, in which the regularized boundary integral equations with indirect unknowns without including the Cauchy principal value （CPV） and Hadamard-finite-part （HFP） integrals are established for the plane potential problems. With some numerical results, it is shown that the better accuracy and higher efficiency, especially on the boundary, can be achieved by the present system.
基金supported by the US ARO grants 49308-MA and 56349-MAthe US AFSOR grant FA9550-06-1-024+1 种基金he US NSF grant DMS-0911434the State Key Laboratory of Scientific and Engineering Computing of Chinese Academy of Sciences during a visit by Z.Li between July-August,2008.
文摘In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.
基金Project supported by the National Natural Science Foundation of China (No. 10962004)the Special-ized Research Fund for the Doctoral Program of Higher Education of China (No. 20070126002)+1 种基金the Chunhui Program of Ministry of Education of China (No. Z2009-1-01010)the Natural Science Foundation of Inner Mongolia (No. 2009BS0101)
文摘This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial differential equations of the above 2D problems is rewritten as an upper triangular differential system. For the associated operator matrix, the existence and the completeness of two normed orthogonal eigenfunction systems in some space are obtained, which belong to the two block operators arising in the operator matrix. Moreover, the general solution to the above 2D problem is given by the eigenfunction expansion method.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China (Grant No. 11102125)
文摘Based on the complex variable moving least-square（CVMLS） approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin（CVMLPG） method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square（MLS） approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin（MLPG） method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method.
基金the National Natural Science Foundation of China (Nos. 10725210 and 10432030) the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20060335107)the Program for New Century Excellent Talents in University, MOE, China (No. NCET-05-05010)
文摘The symplectic approach proposed and developed by Zhong et al. in 1990s for elasticity problems is a rational analytical method, in which ample experience is not needed as in the conventional semi-inverse method. In the symplectic space, elasticity problems can be solved using the method of separation of variables along with the eigenfunction expansion technique, as in traditional Fourier analysis. The eigensolutions include those corresponding to zero and nonzero eigenvalues. The latter group can be further divided into α-and β-sets. This paper reformulates the form of β-set eigensolutions to achieve the stability of numerical calculation, which is very important to obtain accurate results within the symplectic frame. An example is finally given and numerical results are compared and discussed.
