A new complex variable meshless method for transient heat conduction problems

In this paper, based on the improved complex variable moving least-square （ICVMLS） approximation, a new complex variable meshless method （CVMM） for two-dimensional （2D） transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.

A new complex variable element-free Galerkin method for two-dimensional potential problems

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.

An improved complex variable element-free Galerkin method for two-dimensional elasticity problems

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.

Analysis of elastoplasticity problems using an improved complex variable element-free Galerkin 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 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.

The complex variable reproducing kernel particle method for two-dimensional elastodynamics

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.

Complex variable element-free Galerkin method for viscoelasticity problems

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.

Combining the complex variable reproducing kernel particle method and the finite element method for solving transient heat conduction problems

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.

The complex variable meshless local Petrov-Galerkin method of solving two-dimensional potential problems

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.

New complex variable meshless method for advection-diffusion problems

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.

A complex variable meshless local Petrov-Galerkin method for transient heat conduction problems

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.

Analysis of variable coefficient advection-diffusion problems via complex variable reproducing kernel particle method

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.

An Unsteady Two-Dimensional Complex Variable Boundary Element Method

The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology in this work is suitable for modeling diffusion problems with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The underpinning of the modeling approach is to decompose the global initial-boundary value problem into a steady-state component and a transient component. The steady-state component is governed by the Laplace PDE and is modeled using the Complex Variable Boundary Element Method. The transient component is governed by the linear diffusion PDE and is modeled by a linear combination of basis functions that are the products of a two-dimensional Fourier sine series and an exponential function. The global approximation function is the sum of the approximate solutions from the two components. The boundary conditions of the steady-state problem are specified to match the boundary conditions from the global problem so that the CVBEM approximation function satisfies the global boundary conditions. Consequently, the boundary conditions of the transient problem are specified to be continuously zero. The initial condition of the transient component is specified as the difference between the initial condition of the global initial-boundary value problem and the CVBEM approximation of the steady-state solution. Therefore, when the approximate solutions from the two components are summed, the resulting global approximation function approximately satisfies the global initial condition. In this work, it will be demonstrated that the coupled global approximation function satisfies the governing diffusion PDE. Lastly, a procedure for developing streamlines at arbitrary model time is discussed.

A Conceptual Numerical Model of the Wave Equation Using the Complex Variable Boundary Element Method

In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.

围缘扁钢加强的开圆孔无限平板应力解析算法

为了研究围缘扁钢加强的开圆孔无限平板应力解析计算方法,近似认为围缘扁钢区域的平均应力处于广义平面应力状态。运用复变函数论的方法,将应力函数表示为两个待定的解析函数;列出围缘扁钢的边界条件、围缘扁钢与平板之间的应力连接条件和位移连接条件,求出应力函数,根据求出的应力函数便可求出应力。为了描述围缘扁钢减缓孔口应力集中的效果,提出围缘扁钢加强系数的定义。算例计算表明:本文提出的解析算法计算结果与有限元法计算结果吻合良好,围缘扁钢加强的开圆孔无限平板应力普遍小于开圆孔无限平板应力,围缘扁钢加强效果明显;相同体积的加强构件,围缘扁钢的加强效果优于环形加厚板。

磁电弹性材料中含有带四条裂纹的正方形孔的反平面断裂问题

针对反平面剪切力作用下,磁电弹性材料中含有带四条裂纹的正方形孔口的断裂问题进行了深入地探索。基于线弹性断裂理论和Riemann-Schwarz解析延拓定理,利用复变函数方法和留数定理,通过构造合适的数值保角映射函数,将解析函数边值问题转化为柯西积分方程组进行解答,获得了磁电非渗透边界条件下裂纹尖端断裂力学参数的显式表达式。通过与已有结果的对比,验证了方法的有效性。利用数值算例描述了孔洞尺寸、四条裂纹长度和机–电–磁载荷等因素对裂纹扩展参数的影响规律。结果表明:水平右裂纹对孔口裂纹扩展影响更显著;垂直裂纹影响水平裂纹的扩展趋势;在磁电非渗透边界条件下,随着机械载荷的增大,裂纹尖端的应力强度因子变大,而电载荷和磁载荷对裂纹的扩展与机械载荷的大小和方向密切相关。该分析方法为求解复杂多连通域的智能复合材料问题提供了一条有效途径,研究成果为含裂纹的复合材料或结构的优化设计提供了科学依据。

三维二十面体准晶弹性半平面周期接触问题研究

借助平面弹性复变方法,研究了周期刚性压头作用下三维二十面体准晶弹性半平面的周期有限摩擦和周期粘结接触问题.从应力、位移分量的复变函数形式表达式出发,利用半平面Hilbert核积分公式、Plemelj公式,得到两类周期接触问题的解答.针对周期有限摩擦接触问题,得到了三类周期刚性压头(直水平、直倾斜和圆形基底)作用在准晶体半平面上接触应力的封闭解.针对半平面周期粘结接触问题,求得接触边界上作用周期尖劈形位移时接触应力的显式解答.当不考虑相位子场的作用时,文中求得的理论结果可退化到正交各向异性材料平面弹性周期接触问题的对应结果.数值算例说明准晶弹性常数对接触应力分布规律及大小的影响.本文结论可为分析准晶材料压痕实验结果及准晶材料性能提供一定理论依据.

非圆形隧洞衬砌支护下应力与位移的新解法

以往在利用平面弹性复变函数方法进行衬砌支护作用下非圆形隧洞的力学分析时,一般将洞室外域保角映射到象平面的单位圆外域,并认为可以利用同一个映射函数将衬砌截面映射到象平面的圆环域内。当衬砌厚度很薄时,这种方法是可行的,而当衬砌相对较厚时,利用同一映射函数所得到的衬砌截面形状发生严重失真。为解决此问题,本文引入两套映射函数,分别将洞室外域映射到象平面的单位圆的圆外域,衬砌截面映射到圆环域。这种方法可以适应于不同的衬砌厚度,利用这种方式获得的衬砌截面,在衬砌较厚时仍能保证衬砌厚度相对均匀。在求解解析函数过程中不再采用传统的幂级数解法,而是使用边界配点法,将两套映射函数结合,获得了该问题的应力解和位移解。本文以马蹄形隧洞为例,求解结果与ANSYS软件所得数值解吻合很好,分析了衬砌厚度对隧洞应力和位移的影响。根据本文解答可以方便、快捷地进行相近条件下隧道的初步设计。

基于复变函数理论的无中导连拱隧道应力与位移解析解

无中导连拱隧道由于后行洞衬砌搭接于先行洞衬砌拱腰之上,围岩与衬砌受力情况复杂,开展相应的围岩与衬砌应力和位移分析具有重要的工程意义。基于复变函数理论和Schwarz交替法,提出了深埋条件下考虑衬砌的无中导连拱隧道围岩与衬砌应力和变形的计算方法。针对无中导连拱隧道两洞具有共用边界的特点,在利用柯西积分法进行求解时,将单洞衬砌与围岩接触的边界积分条件转化为仅存在单洞时的完整边界积分与共用边界积分之差,对于Schwarz交替法中求解“附加面力”的过程亦进行类似处理。给出共用边界应力与位移的计算方法。通过算例,将上述解析结果和数值仿真结果进行对比,证明了方法的可行性和准确性。

基于子空间辨识方法的建筑群热负荷预测研究

根据建筑热力模型和热负荷计算模型建立了状态空间模型,分别利用未指定状态变量和指定状态变量的子空间辨识方法构建了建筑群热负荷预测模型,并以哈尔滨市某小区换热站为例,验证和比较了2种方法的热负荷预测结果及训练策略和预测时长对预测精度的影响。结果表明:指定状态变量比未指定状态变量的子空间辨识方法的热负荷预测平均绝对百分比误差降低了10.00%;未指定状态变量的子空间辨识方法的中期预测效果优于短期预测,平均绝对百分比误差降低了13.34%。

Analytical investigations of in situ stress inversion from borehole breakout geometries

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.