In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least s...In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.展开更多
In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squar...In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares （CVMLS） approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin （EFG） method and the ICVMLS approximation, the improved complex variable element-free Galerkin （ICVEFG） method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.展开更多
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.展开更多
In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved c...In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin（ICVEFG） method is presented for two-dimensional（2D） elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.展开更多
Based on the complex variable moving least-square （CVMLS） approximation, the complex variable element-free Galerkin （CVEFG） method for two-dimensional viscoelasticity problems under the creep condition is presente...Based on the complex variable moving least-square （CVMLS） approximation, the complex variable element-free Galerkin （CVEFG） method for two-dimensional viscoelasticity problems under the creep condition is presented in this paper. The Galerkin weak form is employed to obtain the equation system, and the penalty method is used to apply the essential boundary conditions, then the corresponding formulae of the CVEFG method for two-dimensional viscoelasticity problems under the creep condition are obtained. Compared with the element-free Galerkin （EFG） method, with the same node distribution, the CVEFG method has higher precision, and to obtain the similar precision, the CVEFG method has greater computational efficiency. Some numerical examples are given to demonstrate the validity and the efficiency of the method.展开更多
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.展开更多
Based on the improved interpolating moving least-squares （ⅡMLS） method and the Galerkin weak form, an improved interpolating element-free Galerkin （ⅡEFG） method is presented for two-dimensional elasticity proble...Based on the improved interpolating moving least-squares （ⅡMLS） method and the Galerkin weak form, an improved interpolating element-free Galerkin （ⅡEFG） method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares （IMLS） method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin （IEFG） method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG 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.展开更多
On the basis of the complex variable moving least-square （CVMLS） approximation, a complex variable meshless local Petrov-Galerkin （CVMLPG） method is presented for transient heat conduction problems. The method is ...On the basis of the complex variable moving least-square （CVMLS） approximation, a complex variable meshless local Petrov-Galerkin （CVMLPG） method is presented for transient heat conduction problems. The method is developed based on the CVMLS approximation for constructing shape functions at scattered points, and the Heaviside step function is used as a test function in each sub-domain to avoid the need for a domain integral in symmetric weak form. In the construction of the well-performed shape function, the trial function of a two-dimensional （2D） problem is formed with a one-dimensional （1D） basis function, thus improving computational efficiency. The numerical results are compared with the exact solutions of the problems and the finite element method （FEM）. This comparison illustrates the accuracy as well as the capability of the CVMLPG method.展开更多
By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is propose...By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.展开更多
This paper presents an interpolating element-free Galerkin(IEFG) method for solving the two-dimensional(2D) elastic large deformation problems. By using the improved interpolating moving least-squares method to form s...This paper presents an interpolating element-free Galerkin(IEFG) method for solving the two-dimensional(2D) elastic large deformation problems. By using the improved interpolating moving least-squares method to form shape function, and using the Galerkin weak form of 2D elastic large deformation problems to obtain the discrete equations, we obtain the formulae of the IEFG method for 2D elastic large deformation problems. As the displacement boundary conditions can be applied directly, the IEFG method can acquire higher computational efficiency and accuracy than the traditional element-free Galerkin(EFG)method, which is based on the moving least-squares approximation and can not apply the displacement boundary conditions directly. To analyze the influences of node distribution, scale parameter of influence domain and the loading step on the numerical solutions of the IEFG method, three numerical examples are proposed. The IEFG method has almost the same high accuracy as the EFG method, and for some 2D elastic large deformation problems the IEFG method even has higher computational accuracy.展开更多
The Burton-Miller boundary integral formulation is solved by a complex variable boundary element-free method(CVBEFM)for the boundary-only meshless analysis of acoustic problems with arbitrary wavenumbers.To regularize...The Burton-Miller boundary integral formulation is solved by a complex variable boundary element-free method(CVBEFM)for the boundary-only meshless analysis of acoustic problems with arbitrary wavenumbers.To regularize both strongly singular and hypersingular integrals and to avoid the computation of the solid angle and its normal derivative,a weakly singular Burton-Miller formulation is derived by considering the normal derivative of the solid angle and adopting the singularity subtraction procedures.To facilitate the implementation of the CVBEFM and the approximation of gradients of the boundary variables,a stabilized complex variable moving least-square approximation is selected in the meshless discretization procedure.The results show the accuracy and efficiency of the present CVBEFM and reveal that the method can produce satisfactory results for all wavenumbers,even for extremely large wavenumbers such as k=10000.展开更多
基金supported by the Graduate Student Scientific Research Innovation Project through Research Innovation Fund for Graduate Students in Shanxi Province(Project No.2024KY648).
文摘In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.
基金supported by the National Natural Science Foundation of China (Grant No.11026223)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, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares （CVMLS） approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin （EFG） method and the ICVMLS approximation, the improved complex variable element-free Galerkin （ICVEFG） method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11171208 and U1433104)
文摘In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin（ICVEFG） method is presented for two-dimensional（2D） elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)
文摘Based on the complex variable moving least-square （CVMLS） approximation, the complex variable element-free Galerkin （CVEFG） method for two-dimensional viscoelasticity problems under the creep condition is presented in this paper. The Galerkin weak form is employed to obtain the equation system, and the penalty method is used to apply the essential boundary conditions, then the corresponding formulae of the CVEFG method for two-dimensional viscoelasticity problems under the creep condition are obtained. Compared with the element-free Galerkin （EFG） method, with the same node distribution, the CVEFG method has higher precision, and to obtain the similar precision, the CVEFG method has greater computational efficiency. Some numerical examples are given to demonstrate the validity and the efficiency of the method.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)
文摘Based on the improved interpolating moving least-squares （ⅡMLS） method and the Galerkin weak form, an improved interpolating element-free Galerkin （ⅡEFG） method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares （IMLS） method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin （IEFG） method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG 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.
基金supported by the National Natural Science Foundation of China(Grant No.51078250)the Research Project by Shanxi Scholarship Council of Shanxi Province,China(Grant No.2013-096)the Scientific&Technological Innovation Program for Postgraduates of Taiyuan University of Science and Technology,China(Grant No.20125026)
文摘On the basis of the complex variable moving least-square （CVMLS） approximation, a complex variable meshless local Petrov-Galerkin （CVMLPG） method is presented for transient heat conduction problems. The method is developed based on the CVMLS approximation for constructing shape functions at scattered points, and the Heaviside step function is used as a test function in each sub-domain to avoid the need for a domain integral in symmetric weak form. In the construction of the well-performed shape function, the trial function of a two-dimensional （2D） problem is formed with a one-dimensional （1D） basis function, thus improving computational efficiency. The numerical results are compared with the exact solutions of the problems and the finite element method （FEM）. This comparison illustrates the accuracy as well as the capability of the CVMLPG method.
基金supported by the Natural Science Foundation of Zhejiang Province,China(Grant Nos.LY20A010021,LY19A010002,LY20G030025)the Natural Science Founda-tion of Ningbo City,China(Grant Nos.2021J147,2021J235).
文摘By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.
基金supported by the National Natural Science Foundation of China (Grant No. 11571223)。
文摘This paper presents an interpolating element-free Galerkin(IEFG) method for solving the two-dimensional(2D) elastic large deformation problems. By using the improved interpolating moving least-squares method to form shape function, and using the Galerkin weak form of 2D elastic large deformation problems to obtain the discrete equations, we obtain the formulae of the IEFG method for 2D elastic large deformation problems. As the displacement boundary conditions can be applied directly, the IEFG method can acquire higher computational efficiency and accuracy than the traditional element-free Galerkin(EFG)method, which is based on the moving least-squares approximation and can not apply the displacement boundary conditions directly. To analyze the influences of node distribution, scale parameter of influence domain and the loading step on the numerical solutions of the IEFG method, three numerical examples are proposed. The IEFG method has almost the same high accuracy as the EFG method, and for some 2D elastic large deformation problems the IEFG method even has higher computational accuracy.
基金Project supported by the National Natural Science Foundation of China(No.11971085)the Innovation Research Group Project in Universities of Chongqing of China(No.CXQT19018)+1 种基金the Science and Technology Research Program of Chongqing Municipal Education Commission of China(No.KJZD-M201800501)and the Science and Technology Research Program of Chongqing University of Education of China(No.KY201927C)。
文摘The Burton-Miller boundary integral formulation is solved by a complex variable boundary element-free method(CVBEFM)for the boundary-only meshless analysis of acoustic problems with arbitrary wavenumbers.To regularize both strongly singular and hypersingular integrals and to avoid the computation of the solid angle and its normal derivative,a weakly singular Burton-Miller formulation is derived by considering the normal derivative of the solid angle and adopting the singularity subtraction procedures.To facilitate the implementation of the CVBEFM and the approximation of gradients of the boundary variables,a stabilized complex variable moving least-square approximation is selected in the meshless discretization procedure.The results show the accuracy and efficiency of the present CVBEFM and reveal that the method can produce satisfactory results for all wavenumbers,even for extremely large wavenumbers such as k=10000.
