A new approach is proposed to solve the elastic-plastic fields near the major-axis line of an elliptical hole. The complex variable method is used to determine the elastic fields near the major-axis line of the ellipt...A new approach is proposed to solve the elastic-plastic fields near the major-axis line of an elliptical hole. The complex variable method is used to determine the elastic fields near the major-axis line of the elliptical hole. Then, by using the line field analysis method, the exact and new solutions of the stresses, strains in the plastic zone, the size of the plastic region and the unit normal vector of the elastic-plastic boundary near the major-axis line of the elliptical hole are obtained for an anti-plane elliptical hole in a perfectly elastic-plastic solid. The usual small scale yielding assumptions are not adopted in the analysis. The present method is simple, easy and efficient. The influences of applied mechanical loading on the size of plastic zone are discussed.展开更多
The complex variable method for the plane elasticity theory of icosahedral quasicrystals is developed.Based on the general solution obtained previously,complex representations of stress and displacement components of ...The complex variable method for the plane elasticity theory of icosahedral quasicrystals is developed.Based on the general solution obtained previously,complex representations of stress and displacement components of phonon and phason fields in the quasicrystals are given.With the help of conformal transformation,an analytic solution for the elliptic notch problem of the material is presented.The solution of the Griffith crack problem can be observed as a special case of the results.The stress intensity factor and energy release rate of the crack are also obtained.展开更多
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, 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,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equi...In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency.展开更多
The complex variable reproducing kernel particle method (CVRKPM) of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape...The complex variable reproducing kernel particle method (CVRKPM) of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection-diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method (RKPM) and the element- free Galerkin (EFG) method.展开更多
In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE metho...In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE 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 reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is present...On the basis of the reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.展开更多
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.展开更多
The element-free Galerkin(EFG)method,which constructs shape functions via moving least squares(MLS)approximation,represents a fundamental and widely studied meshless method in numerical computation.Although it achieve...The element-free Galerkin(EFG)method,which constructs shape functions via moving least squares(MLS)approximation,represents a fundamental and widely studied meshless method in numerical computation.Although it achieves high computational accuracy,the shape functions are more complex than those in the conventional finite element method(FEM),resulting in great computational requirements.Therefore,improving the computational efficiency of the EFG method represents an important research direction.This paper systematically reviews significant contributions fromdomestic and international scholars in advancing the EFGmethod.Including the improved element-free Galerkin(IEFG)method,various interpolating EFG methods,four distinct complex variable EFG methods,and a series of dimension splitting meshless methods.In the numerical examples,the effectiveness and efficiency of the three methods are validated by analyzing the solutions of the IEFG method for 3D steadystate anisotropic heat conduction,3D elastoplasticity,and large deformation problems,as well as the performance of two-dimensional splitting meshless methods in solving the 3D Helmholtz equation.展开更多
Grouting reinforcement serves as a critical measure for the stability design and disaster prevention of deep-buried tunnels.However,the grouting reinforcement effect on the unloading stress path of surrounding rock in...Grouting reinforcement serves as a critical measure for the stability design and disaster prevention of deep-buried tunnels.However,the grouting reinforcement effect on the unloading stress path of surrounding rock in high tectonic stress zones and its mechanisms for disaster prevention and control remain unclear.This paper proposes a functional gradient bearing concept for grouting reinforcement of fractured surrounding rock in deep-buried tunnels.The functionally graded reinforced zone is designed,exhibiting a radial mechanical property gradient where stiffness increases and strength decreases from shallow to deep regions of the surrounding rock.A semi-analytical method based on the complex variable theory is specificallydeveloped for elastic-plastic mechanical analysis of reinforced surrounding rock.Subsequently,the effects of rock strength and stiffness on the morphological characteristics of the plastic zone around deep-buried circular tunnels in high tectonic stress zones are investigated.The study shows that:(1)The annular grouting reinforcement design has limited effectiveness in preventing the failure of surrounding rock;(2)Grouting reinforcement significantlyenhances the bearing capacity of tunnels in weak and fractured rock masses,but excessive grouting to control the plastic zone is not recommended;(3)The bearing capacity of the shallow surrounding rock and the stress regulating function of the deeper surrounding rock can be fully utilized if the functional gradient bearing concept is considered,achieving coordinated bearing between shallow and deep surrounding rock.The research provides new insights for the gradient grouting reinforcement and stress control optimization design of surrounding rock in deep-buried circular tunnels.展开更多
Based on the nondestructive test data of operating railway tunnels in China, this paper summarizes the basic characteristics of the complex contact behavior between the rock mass and lining structure. The contact mode...Based on the nondestructive test data of operating railway tunnels in China, this paper summarizes the basic characteristics of the complex contact behavior between the rock mass and lining structure. The contact modes are classified into dense contact, local non-contact, and loose contact. Subsequently, the corresponding mechanical model for each contact mode is developed according to its mechanical characteristics using the complex variable method. In the proposed mechanical model, a special algorithm is introduced to detect whether the local non-contact zone is re-contacted. Besides, a novel conformal mapping method is also proposed to accurately calculate the mechanical response of the concrete lining. Finally, the accuracy of the proposed method is verified by comparing it with the finite element method(FEM). Several parameter investigations are conducted to analyze the effects of different contact modes on the rock-lining interaction. The results show that:(i) the height of the local noncontact area does not have a significant effect on the contact stress distribution if no re-contact occurs;(ii) backfill grouting can reduce the local stress concentration caused by poor contact modes;and(iii) reducing the friction coefficient of the interface can lead to a more uniform distribution of internal forces in the concrete lining.展开更多
On the basis of reproducing kernel particle method(RKPM),using complex variable theory,the complex variable reproducing kernel particle method(CVRKPM)is discussed in this paper.The advantage of the CVRKPM is that the ...On the basis of reproducing kernel particle method(RKPM),using complex variable theory,the complex variable reproducing kernel particle method(CVRKPM)is discussed in this paper.The advantage of the CVRKPM is that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is formed.Then the CVRKPM is applied to solve two-dimensional elasto-plasticity problems.The Galerkin weak form is employed to obtain the discretized system equation,the penalty method is used to apply the essential boundary conditions.And then,the CVRKPM for two-dimensional elasto-plasticity problems is formed,the corresponding formulae are obtained,and the Newton-Raphson method is used in the numerical implementation.Three numerical examples are given to show that this method in this paper is effective for elasto-plasticity analysis.展开更多
Based on the moving least-square(MLS)approximation,the complex variable moving least-square approximation(CVMLS)is discussed in this paper.The complex variable moving least-square approximation cannot form ill-conditi...Based on the moving least-square(MLS)approximation,the complex variable moving least-square approximation(CVMLS)is discussed in this paper.The complex variable moving least-square approximation cannot form ill-conditioned equations,and has greater precision and computational efficiency.Using the analytical solution near the tip of a crack,the trial functions in the complex variable moving least-square approxi-mation are extended,and the corresponding approximation function is obtained.And from the minimum potential energy principle,a complex variable meshless method for fracture problems is presented,and the formulae of the complex variable meshless method are obtained.The complex variable meshless method in this paper has greater precision and computational efficiency than the conventional meshless method.Some examples are given.展开更多
A stochastic simulation of fluid flow in porous media using a complex variable expression method (SFCM) is presented in this paper. Hydraulic conductivity is considered as a random variable and is then expressed in ...A stochastic simulation of fluid flow in porous media using a complex variable expression method (SFCM) is presented in this paper. Hydraulic conductivity is considered as a random variable and is then expressed in complex variable form, the real part of which is a deterministic value and the imaginary part is a variable value. The stochastic seepage flow is simulated with the SFCM and is compared with the results calculated with the Monte Carlo stochastic finite element method. In using the Monte Carlo method to simulate the stochastic seepage flow field, the hydraulic conductivity is assumed in three different probability distributions using random sampling method. The obtained seepage flow field is examined through skewness analysis, and the skewed distribution probability density function is given. The head mode value and the head comprehensive standard deviation are used to represent the statistics of calculation results obtained by the Monte Carlo method. The stochastic seepage flow field simulated by the SFCM is confirmed to be similar to that given by the Monte Carlo method from numerical aspects. The range of coefficient of variation of hydraulic conductivity in SFCM is larger than used previously in stochastic seepage flow field simulations, and the computation time is short. The results proved that the SFCM is a convenient calculating method for solving the complex problems.展开更多
The elastic interaction between a screw dislocation and an elliptical inhomogeneity with interfacial cracks is studied. The screw dislocation may be located outside or inside the inhomogeneity. An efficient complex va...The elastic interaction between a screw dislocation and an elliptical inhomogeneity with interfacial cracks is studied. The screw dislocation may be located outside or inside the inhomogeneity. An efficient complex variable method for the complex multiply connected region is developed, and the general solutions to the problem are derived. As illustrative examples, solutions in explicit series form for complex potentials are presented in the case of one or two interfacial cracks. Image forces on the dislocation are calculated by using the Peach-Koehler formula. The influence of crack geometries and material properties on the image forces is evaluated and discussed. It is shown that the interfacial crack has a significant effect on the equilibrium position of the dislocation near an elliptical-arc interface. The main results indicate, when the length of the crack goes up to a critical value, the presence of the interfacial crack can change the interaction mechanism between a screw dislocation and an elliptical inclusion. The present solutions can include a number of previously known results as special cases.展开更多
Using the complex variable method and conformal mapping,scat- tering of flexural waves and dynamic stress concentrations in Mindlin's thick plates with a cutout have been studied.The general solution of the stress...Using the complex variable method and conformal mapping,scat- tering of flexural waves and dynamic stress concentrations in Mindlin's thick plates with a cutout have been studied.The general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of cutouts is obtained. Applying the orthogonal function expansion technique,the dynamic stress problem can be reduced into the solution of a set of infinite algebraic equations.As examples, numerical results for the dynamic stress concentration factor in Mindlin's plates with a circular,elliptic cutout are graphically presented in sequence.展开更多
This study aims to investigate the feasibility of deriving in situ horizontal stresses from the breakout width and depth using the analytical method.Twenty-three breakout data with different borehole sizes were collec...This study aims to investigate the feasibility of deriving in situ horizontal stresses from the breakout width and depth using the analytical method.Twenty-three breakout data with different borehole sizes were collected and three failure criteria were studied.Based on the Kirsch equations,relatively accurate major horizontal stress(sH)estimations from known minor horizontal stress(sh)were achieved with percentage errors ranging from 0.33%to 44.08%using the breakout width.The Mogi-Coulomb failure criterion(average error:13.1%)outperformed modified Wiebols-Cook(average error:19.09%)and modified Lade(average error:18.09%)failure criteria.However,none of the tested constitutive models could yield reasonable sh predictions from known sH using the same approach due to the analytical expression of the redistributed stress and the nature of the constitutive models.In consideration of this issue,the horizontal stress ratio(sH/sh)is suggested as an alternative input,which could estimate both sH and sh with the same level of accuracy.Moreover,the estimation accuracies for both large-scale and laboratory-scale breakouts are comparable,suggesting the applicability of this approach across different breakout sizes.For breakout depth,conformal mapping and complex variable method were used to calculate the stress concentration around the breakout tip,allowing the expression of redistributed stresses using binomials composed of sH and sh.Nevertheless,analysis of the breakout depth stabilisation mechanism indicates that additional parameters are required to utilise normalised breakout depth for stress estimation compared to breakout width.These parameters are challenging to obtain,especially under field conditions,meaning utilising normalised breakout depth analytically in practical applications faces significant challenges and remains infeasible at this stage.Nonetheless,the normalised breakout depth should still be considered a critical input for any empirical and statistical stress estimation method given its significant correlation with horizontal stresses.The outcome of this paper is expected to contribute valuable insights into the breakout stabilisation mechanisms and estimation of in situ stress magnitudes based on borehole breakout geometries.展开更多
The antiplane problem of circular arc rigid line inclusions under antiplane concentrated force and longitudinal shear loading was dealt with. By using Riemann-Schwarz's symmetry principle integrated with the singu...The antiplane problem of circular arc rigid line inclusions under antiplane concentrated force and longitudinal shear loading was dealt with. By using Riemann-Schwarz's symmetry principle integrated with the singularity analysis of complex functions, the general solution of the problem and the closed form solutions for some important practical problems were presented. The stress distribution in the immediate vicinity of circular arc rigid line end was examined in detail. The results show that the singular stress fields near the rigid inclusion tip possess a square-root singularity similar to that for the corresponding crack problem under antiplane shear loading, but no oscillatory character. Furthermore, the stresses are found to depend on geometrical dimension, loading conditions and materials parameters. Some practical results concluded are in agreement with the previous solutions.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 10932001 and 11072015)the Scientific Research Key Program of Beijing Municipal Commission of Education (Grant No. KZ201010005003)the PhD Innovative Foundation of Beihang University (Grant No. 300351)
摘要A new approach is proposed to solve the elastic-plastic fields near the major-axis line of an elliptical hole. The complex variable method is used to determine the elastic fields near the major-axis line of the elliptical hole. Then, by using the line field analysis method, the exact and new solutions of the stresses, strains in the plastic zone, the size of the plastic region and the unit normal vector of the elastic-plastic boundary near the major-axis line of the elliptical hole are obtained for an anti-plane elliptical hole in a perfectly elastic-plastic solid. The usual small scale yielding assumptions are not adopted in the analysis. The present method is simple, easy and efficient. The influences of applied mechanical loading on the size of plastic zone are discussed.
基金the National Natural Science Foundation of China(Grant Nos.10372016 and 10761005)the Natural Science Foundation of Inner Mongolia of China(Grant No.200607010104)the Natural Science Foundation of Inner Mongolia Normal University(Grant No.QN07034)
摘要The complex variable method for the plane elasticity theory of icosahedral quasicrystals is developed.Based on the general solution obtained previously,complex representations of stress and displacement components of phonon and phason fields in the quasicrystals are given.With the help of conformal transformation,an analytic solution for the elliptic notch problem of the material is presented.The solution of the Griffith crack problem can be observed as a special case of the results.The stress intensity factor and energy release rate of the crack are also obtained.
基金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.
基金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 for Graduate Student of Shanghai University,China (Grant No. SHUCX120125)
摘要In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency.
基金supported by the National Natural Science Foundation of China (Grant No. 11171208)the Leading Academic Discipline Project of Shanghai City,China (Grant No. S30106)
摘要The complex variable reproducing kernel particle method (CVRKPM) of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection-diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method (RKPM) and the element- free Galerkin (EFG) method.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Special Fund for Basic Scientific Research of Central Colleges of Chang’an University, China (Grant No. CHD2011JC080)
摘要In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE 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.10871124)the Innovation Program of Shanghai Municipal Education Commission,China (Grant No.09ZZ99)
摘要On the basis of the reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.
基金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 National Natural Science Foundation of China(Grant No.12271341).
摘要The element-free Galerkin(EFG)method,which constructs shape functions via moving least squares(MLS)approximation,represents a fundamental and widely studied meshless method in numerical computation.Although it achieves high computational accuracy,the shape functions are more complex than those in the conventional finite element method(FEM),resulting in great computational requirements.Therefore,improving the computational efficiency of the EFG method represents an important research direction.This paper systematically reviews significant contributions fromdomestic and international scholars in advancing the EFGmethod.Including the improved element-free Galerkin(IEFG)method,various interpolating EFG methods,four distinct complex variable EFG methods,and a series of dimension splitting meshless methods.In the numerical examples,the effectiveness and efficiency of the three methods are validated by analyzing the solutions of the IEFG method for 3D steadystate anisotropic heat conduction,3D elastoplasticity,and large deformation problems,as well as the performance of two-dimensional splitting meshless methods in solving the 3D Helmholtz equation.
基金funded by China National Postdoctoral Program for Innovative Talents(Grant No.BX20230267)National Natural Science Foundation of China(Grant Nos.42402274,42307206).
摘要Grouting reinforcement serves as a critical measure for the stability design and disaster prevention of deep-buried tunnels.However,the grouting reinforcement effect on the unloading stress path of surrounding rock in high tectonic stress zones and its mechanisms for disaster prevention and control remain unclear.This paper proposes a functional gradient bearing concept for grouting reinforcement of fractured surrounding rock in deep-buried tunnels.The functionally graded reinforced zone is designed,exhibiting a radial mechanical property gradient where stiffness increases and strength decreases from shallow to deep regions of the surrounding rock.A semi-analytical method based on the complex variable theory is specificallydeveloped for elastic-plastic mechanical analysis of reinforced surrounding rock.Subsequently,the effects of rock strength and stiffness on the morphological characteristics of the plastic zone around deep-buried circular tunnels in high tectonic stress zones are investigated.The study shows that:(1)The annular grouting reinforcement design has limited effectiveness in preventing the failure of surrounding rock;(2)Grouting reinforcement significantlyenhances the bearing capacity of tunnels in weak and fractured rock masses,but excessive grouting to control the plastic zone is not recommended;(3)The bearing capacity of the shallow surrounding rock and the stress regulating function of the deeper surrounding rock can be fully utilized if the functional gradient bearing concept is considered,achieving coordinated bearing between shallow and deep surrounding rock.The research provides new insights for the gradient grouting reinforcement and stress control optimization design of surrounding rock in deep-buried circular tunnels.
基金supported by the National Natural Science Foundation of China (Grant Nos. 51738002 and 52108376)Fundamental Research Funds for the Central Universities (Grant No. 2021CZ111)
摘要Based on the nondestructive test data of operating railway tunnels in China, this paper summarizes the basic characteristics of the complex contact behavior between the rock mass and lining structure. The contact modes are classified into dense contact, local non-contact, and loose contact. Subsequently, the corresponding mechanical model for each contact mode is developed according to its mechanical characteristics using the complex variable method. In the proposed mechanical model, a special algorithm is introduced to detect whether the local non-contact zone is re-contacted. Besides, a novel conformal mapping method is also proposed to accurately calculate the mechanical response of the concrete lining. Finally, the accuracy of the proposed method is verified by comparing it with the finite element method(FEM). Several parameter investigations are conducted to analyze the effects of different contact modes on the rock-lining interaction. The results show that:(i) the height of the local noncontact area does not have a significant effect on the contact stress distribution if no re-contact occurs;(ii) backfill grouting can reduce the local stress concentration caused by poor contact modes;and(iii) reducing the friction coefficient of the interface can lead to a more uniform distribution of internal forces in the concrete lining.
基金supported by the National Natural Science Foundation of China(Grant Nos.10571118 and 10871124)Innovation Program of Shanghai Municipal Education Commission(Grant No.09ZZ99)
摘要On the basis of reproducing kernel particle method(RKPM),using complex variable theory,the complex variable reproducing kernel particle method(CVRKPM)is discussed in this paper.The advantage of the CVRKPM is that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is formed.Then the CVRKPM is applied to solve two-dimensional elasto-plasticity problems.The Galerkin weak form is employed to obtain the discretized system equation,the penalty method is used to apply the essential boundary conditions.And then,the CVRKPM for two-dimensional elasto-plasticity problems is formed,the corresponding formulae are obtained,and the Newton-Raphson method is used in the numerical implementation.Three numerical examples are given to show that this method in this paper is effective for elasto-plasticity analysis.
摘要Based on the moving least-square(MLS)approximation,the complex variable moving least-square approximation(CVMLS)is discussed in this paper.The complex variable moving least-square approximation cannot form ill-conditioned equations,and has greater precision and computational efficiency.Using the analytical solution near the tip of a crack,the trial functions in the complex variable moving least-square approxi-mation are extended,and the corresponding approximation function is obtained.And from the minimum potential energy principle,a complex variable meshless method for fracture problems is presented,and the formulae of the complex variable meshless method are obtained.The complex variable meshless method in this paper has greater precision and computational efficiency than the conventional meshless method.Some examples are given.
基金supported by the National Natural Science Foundation of China(GrantNos.51079039,51009053)
摘要A stochastic simulation of fluid flow in porous media using a complex variable expression method (SFCM) is presented in this paper. Hydraulic conductivity is considered as a random variable and is then expressed in complex variable form, the real part of which is a deterministic value and the imaginary part is a variable value. The stochastic seepage flow is simulated with the SFCM and is compared with the results calculated with the Monte Carlo stochastic finite element method. In using the Monte Carlo method to simulate the stochastic seepage flow field, the hydraulic conductivity is assumed in three different probability distributions using random sampling method. The obtained seepage flow field is examined through skewness analysis, and the skewed distribution probability density function is given. The head mode value and the head comprehensive standard deviation are used to represent the statistics of calculation results obtained by the Monte Carlo method. The stochastic seepage flow field simulated by the SFCM is confirmed to be similar to that given by the Monte Carlo method from numerical aspects. The range of coefficient of variation of hydraulic conductivity in SFCM is larger than used previously in stochastic seepage flow field simulations, and the computation time is short. The results proved that the SFCM is a convenient calculating method for solving the complex problems.
基金The project supported by the National Natural Science Foundation of China(10272009 and 10472030)the Natural Science Foundation of Hunan Province(02JJY2014)
摘要The elastic interaction between a screw dislocation and an elliptical inhomogeneity with interfacial cracks is studied. The screw dislocation may be located outside or inside the inhomogeneity. An efficient complex variable method for the complex multiply connected region is developed, and the general solutions to the problem are derived. As illustrative examples, solutions in explicit series form for complex potentials are presented in the case of one or two interfacial cracks. Image forces on the dislocation are calculated by using the Peach-Koehler formula. The influence of crack geometries and material properties on the image forces is evaluated and discussed. It is shown that the interfacial crack has a significant effect on the equilibrium position of the dislocation near an elliptical-arc interface. The main results indicate, when the length of the crack goes up to a critical value, the presence of the interfacial crack can change the interaction mechanism between a screw dislocation and an elliptical inclusion. The present solutions can include a number of previously known results as special cases.
基金The project supported by the National Natural Science Foundation of China
摘要Using the complex variable method and conformal mapping,scat- tering of flexural waves and dynamic stress concentrations in Mindlin's thick plates with a cutout have been studied.The general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of cutouts is obtained. Applying the orthogonal function expansion technique,the dynamic stress problem can be reduced into the solution of a set of infinite algebraic equations.As examples, numerical results for the dynamic stress concentration factor in Mindlin's plates with a circular,elliptic cutout are graphically presented in sequence.
基金funded by the Australian Coal Industry’s Research Program(ACARP,Grant No.C26063).
摘要This study aims to investigate the feasibility of deriving in situ horizontal stresses from the breakout width and depth using the analytical method.Twenty-three breakout data with different borehole sizes were collected and three failure criteria were studied.Based on the Kirsch equations,relatively accurate major horizontal stress(sH)estimations from known minor horizontal stress(sh)were achieved with percentage errors ranging from 0.33%to 44.08%using the breakout width.The Mogi-Coulomb failure criterion(average error:13.1%)outperformed modified Wiebols-Cook(average error:19.09%)and modified Lade(average error:18.09%)failure criteria.However,none of the tested constitutive models could yield reasonable sh predictions from known sH using the same approach due to the analytical expression of the redistributed stress and the nature of the constitutive models.In consideration of this issue,the horizontal stress ratio(sH/sh)is suggested as an alternative input,which could estimate both sH and sh with the same level of accuracy.Moreover,the estimation accuracies for both large-scale and laboratory-scale breakouts are comparable,suggesting the applicability of this approach across different breakout sizes.For breakout depth,conformal mapping and complex variable method were used to calculate the stress concentration around the breakout tip,allowing the expression of redistributed stresses using binomials composed of sH and sh.Nevertheless,analysis of the breakout depth stabilisation mechanism indicates that additional parameters are required to utilise normalised breakout depth for stress estimation compared to breakout width.These parameters are challenging to obtain,especially under field conditions,meaning utilising normalised breakout depth analytically in practical applications faces significant challenges and remains infeasible at this stage.Nonetheless,the normalised breakout depth should still be considered a critical input for any empirical and statistical stress estimation method given its significant correlation with horizontal stresses.The outcome of this paper is expected to contribute valuable insights into the breakout stabilisation mechanisms and estimation of in situ stress magnitudes based on borehole breakout geometries.
摘要The antiplane problem of circular arc rigid line inclusions under antiplane concentrated force and longitudinal shear loading was dealt with. By using Riemann-Schwarz's symmetry principle integrated with the singularity analysis of complex functions, the general solution of the problem and the closed form solutions for some important practical problems were presented. The stress distribution in the immediate vicinity of circular arc rigid line end was examined in detail. The results show that the singular stress fields near the rigid inclusion tip possess a square-root singularity similar to that for the corresponding crack problem under antiplane shear loading, but no oscillatory character. Furthermore, the stresses are found to depend on geometrical dimension, loading conditions and materials parameters. Some practical results concluded are in agreement with the previous solutions.