各向异性柱体扭转的充分必要的边界积分方程

本文列出各向异性柱体扭转问题无量纲化后的有关方程的边界条件推导和验证了基本解，并指出一些书中基本解列式有误＾［９－１１］，列出了充分必要的边界积分方 程，进行了数据计算，并与习用的边界积分方程所得结果进行了比较，表明在退化值附近，习用的边界积分方程所得的边界剪应力会出现巨大的误差，扭转刚度的误差则要小得多，而充要的边界积分方程计算的结果则如终保持良好的精度，再次显示了它的优点。     

弹性力学中一种新的边界轮廓法

利用基本解的特性,将面力积分方程化成仅含有Ｃａｕｃｈｙ主值积分的形式,基于这种边界积分方程,提出了一种新的边界轮廓法,对于三维问题,该方法只须计算沿边界单元界线的线积分,对二维问题,则只需计算边界单元两点的热函数之差,无须进行数值积分计算,实例计算说明该方法是有效的。     

A NECESSARY AND SUFFICIENT FORMULATION FOR TORSION PROBLEMS OF ANISOTROPIC BODIES

Non-dimensionalized equations and boundary conditions are presentedfor the torsion problem of an anisotropic body.The error of the fundamental solutioncited in some boundary element books is pointed out after an examination of thefundamental solution.Furthermore,a necessary and sufficient boundary integralequation is given for the problem and compared with the conventional boundaryintegral equation.Numerical results show that great errors of the boundary shearstresses obtained by the conventional boundary integral equation appear with a smallerror of torsion stiffness.Meanwhile,the necessary and sufficient boundary integralequation always gives accurate results.     

轴流压气机端壁附面层的理论预测

从Ｎ－Ｓ方程出发，推导了轴流压气机端壁边界层动量积分方程，并提出了一种适用范围广泛的叶片力亏损模型，结合边界层的卷吸方程和速度分布表达式，可以方便的预测压气机环壁边界层的发展，与传统的叶片力亏损模型相比，该方法的计算准确性和适用范围均得到提高。     

Dirichlet and Neumann boundary conditions for the pressure poisson equation of incompressible flow

In a recent paper Gresho and Sani showed that Dirichlet and Neumann boundary conditions for the pressure Poisson equation give the same solution. The purpose of this paper is to confirm this (for one case at least) by numerically solving the pressure equation with Dirichlet and Neumann boundary conditions for the inviscid stagnation point flow problem. The Dirichlet boundary condition is obtained by integrating the tangential component of the momentum equation along the boundary. The Neumann boundary condition is obtained by applying the normal component of the momentum equation at the boundary. In this work solutions for the Neumann problem exist only if a compatibility condition is satisfied. A consistent finite difference procedure which satisfies this condition on non-staggered grids is used for the solution of the pressure equation with Neumann conditions. Two test cases are computed. In the first case the velocity field is given from the analytical solution and the pressure is recovered from the solution of the associated Poisson equation. The computed results are identical for both Dirichlet and Neumann boundary conditions. However, the Dirichlet problem converges faster than the Neumann case. In the second test case the velocity field is computed from the momentum equations, which are solved iteratively with the pressure Poisson equation. In this case the Neumann problem converges faster than the Dirichlet problem.     

简支梯形底扁球壳弯曲问题的准格林函数方法

以简支梯形底扁球壳的弯曲问题为例,详细阐明了准格林函数方法的思想.即利用问题的基本解和边界方程构造一个准格林函数,这个函数满足了问题的齐次边界条件,采用格林公式将简支扁球壳弯曲问题的控制微分方程化为两个互相耦合的第二类Fredholm积分方程.边界方程有多种选择,在选定一种边界方程的基础上,可以通过建立一个新的边界方程...     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

Application of boundary integral equations for analyzing the dynamics of elastic,viscoelastic, and poroelastic bodies

Two approaches (classical and nonclassical) of the boundary integral equation method for solving three-dimensional dynamical boundary value problems of elasticity, viscoelasticity, and poroelasticity are considered. The boundary integral equation model is used for porous materials. The Kelvin–Voigt model and the weakly singular hereditary Abel kernel are used to describe the viscoelastic properties. Both approaches permit solving the dynamic problems exactly not only in the isotropic but also in the anisotropic case. The boundary integral equation solution scheme is constructed on the basis of the boundary element technique. The numerical results obtained by the classical and nonclassical approaches are compared.     

薄板弯曲问题的虚边界元-最小二乘法

本文提出求解任意形状的薄板弯曲问题的虚边界元-最小二乘法。本法首先利用薄板弯曲平衡方程的格林函数和离开实际边界上分布的未知的横向荷载和法向弯矩函数建立满足实际边界条件的积分方程;然后采用最小二乘法和沿虚边界分段离散化的待定的分布横向荷载和法向弯矩函数得到求上述积分方程离散化数值解的线性代数方程组。导出了一系列的数值积分的公式,并求解了许多例题,数值结果说明本法完全避免了奇异积分及其复杂的处理方法和耗时的运算,而且在边界及其附近区域解的精度比普通边界元(以后简称边界元)法大大地提高了。     

非惯性坐标系中的边界层方程

本文用匹配渐近展开法导出二维翼型在非惯性坐标系中的一阶和二阶边界层方程,消去压强项后,一阶边界层方程与经典边界层方稃相同;惯性力的作用在二阶边界层中才出现。     

一类平衡大气边界层边界条件研究

基于标准k-ε湍流模型,首先利用湍流粘度方程和剪切应力在整个边界层内恒定的假设,推导出一类耗散率表达式,并根据常用的湍动能入口剖面方程以及平均风速剖面方程,计算获得相应的耗散率方程;然后在输运方程中添加自定义源项,通过已经确定的平均速度方程、湍动能方程、耗散率方程计算得到相应输运方程的自定义源项表达式,并进行空风洞数值模拟,从而得到了一类满足平衡大气边界层的来流边界条件．通过将这种边界条件与由湍流平衡条件得到的边界条件进行比较,表明本方法获得的边界条件更适用．并且,本方法无需考虑修正壁面函数和修正湍流模型常数,因而计算更为简单,可为平衡大气边界层的研究提供一种新的思路．     

Green quasifunction method for vibration of simply-supported thin polygonic plates on Pasternak foundation

A new numerical method—Green quasifunction is proposed.The idea of Green quasifunction method is clarified in detail by considering a vibration problem of simply-supported thin polygonic plates on Pasternak foundation.A Green quasifunction is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the problem.The mode shape differential equation of the vibration problem of simply-supported thin plates on Pasternak foundation is reduced to two simultaneous Fredholm integral equations of the second kind by Green formula.There are multiple choices for the normalized boundary equation.Based on a chosen normalized boundary equation,a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome.Finally,natural frequency is obtained by the condition that there exists a nontrivial solution in the numerically discrete algebraic equations derived from the integral equations.Numerical results show high accuracy of the Green quasifunction method.     

用面力边界积分方程求解断裂力学J积分

证明面力边界积分方程被积函数的散度等于零，应用Stokes公式，对平面线弹性问题，将面力边界积分的求解转化为边界点的位移势函数的点值计算。应用边界积分方程的求解结果，推导出J积分亦可表示为边界点的积分势函数的点值计算，无需进行数值积分，实例计算说明该方法的有效性。     

Quasi-Green’s function method for free vibration of simply-supported trapezoidal shallow spherical shell

The idea of quasi-Green’s function method is clarified by considering a free vibration problem of the simply-supported trapezoidal shallow spherical shell. A quasi-Green’s function is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the prob-lem. The mode shape differential equations of the free vibration problem of a simply-supported trapezoidal shallow spherical shell are reduced to two simultaneous Fredholm integral equations of the second kind by the Green formula. There are multiple choices for the normalized boundary equation. Based on a chosen normalized boundary equa-tion, a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome. Finally, natural frequency is obtained by the condition that there exists a nontrivial solution to the numerically discrete algebraic equations derived from the integral equations. Numerical results show high accuracy of the quasi-Green’s function method.     

三维势流场的比例边界有限元求解方法

比例边界有限元法(SBFEM)是线性偏微分方程的一种新的数值求解方法。该方法只对计算域边界利用Galerkin方法进行数值离散,相对于有限元方法(FEM)减少了一个空间坐标的维数,而在减少的空间坐标方向利用解析方法进行求解;相对于边界元法(BEM),比例边界有限元方法不需要基本解,避免了奇异积分的计算,所以它结合了有限元和边界元方法的优点。本文建立了利用比例边界有限元法求解三维Laplace方程的数值模型并用于计算三维物体周围的水流场,将计算结果与解析解和边界元方法进行了对比,结果表明此方法可以很好地模拟水流场,且具有较高的计算精度。     

Quasi-Green＇s function method for free vibration of simply-supported trapezoidal shallow spherical shell

The idea of quasi-Green＇s function method is clarified by considering a free vibration problem of the simply-supported trapezoidal shallow spherical shell. A quasi- Green＇s function is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the prob- lem. The mode shape differential equations of the free vibration problem of a simply- supported trapezoidal shallow spherical shell are reduced to two simultaneous Fredholm integral equations of the second kind by the Green formula. There are multiple choices for the normalized boundary equation. Based on a chosen normalized boundary equa- tion, a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome. Finally, natural frequency is obtained by the condition that there exists a nontrivial solution to the numerically discrete algebraic equations derived from the integral equations. Numerical results show high accuracy of the quasi-Green＇s function method.     

Quasi-Green’s function method for free vibration of clamped thin plates on Winkler foundation

The quasi-Green’s function method is used to solve the free vibration problem of clamped thin plates on the Winkler foundation. Quasi-Green’s function is established by the fundamental solution and the boundary equation of the problem. The function satisfies the homogeneous boundary condition of the problem. The mode-shape differential equation of the free vibration problem of clamped thin plates on the Winkler foundation is reduced to the Fredholm integral equation of the second kind by Green’s formula. The irregularity of the kernel of the integral equation is overcome by choosing a suitable form of the normalized boundary equation. The numerical results show the high accuracy of the proposed method.     

基于比例边界有限元法的高阶透射边界应用

基于比例边界有限元法和连分式展开推导了无限域弹性动力分析的求解方程,实现了一种局部的高阶透射边界. 采用改进的连分式法求解无限域的动力刚度矩阵,克服了原连分式算法可能会造成矩阵运算病态的问题. 该局部高阶透射边界在时域里表示为一阶常微分方程组,其稳定性取决于其系数矩阵的广义特征值问题. 如果出现虚假模态,采用移谱法来校正系数矩阵以消除虚假模态. 通过两个算例验证了该高阶透射边界的精确性、鲁棒性.     

一类Signorini问题的边界变分不等式

本文讨论了一类简化的Signorini问题。首先将原问题和一个边值问题建立联系，其次将原问题的解分解为不带不等边界条件的变分方程的解和一个变分不等式的解。然后利用边值问题的边界积分方程将变分不等式等价地化解为边界变分不等式。这样原求区域上的第一类椭圆变分不等式问题化解为求一个区域上的变分方程和一个边界变分不等式。最后说明了边界变分不等式解的存在唯一性。文末计算了柱面和半无限刚性基础的摩擦接触问题。结论表明文中方法具有较好的精度。     

Multiple reciprocity method with two series of sequences of high-order fundamental solution for thin plate bending

IntroductionIt’swell_knownthatthecomplicatedfundamentalsolution[1,2 ]forHelmholtzequationΔu(x) +k2 u(x) =0　　(x∈Ω:boundedopenregioninR2 )isu (x,y) =-iH(2 )0 (k x-y ) 4,thusit’snotconvenientfornumericalcomputation .IfapplyingthesimplefundamentalsolutionofLaplaceequationu 0 (x ,y) =-ln｜x-y｜(2π) ,theexpressionforthesolutionofequationintheclosedregion Ωisc(y)u(y) + ∫Γu(x) u 0 (x,y) nx -u 0 (x ,y) u(x) n dsx =-k2∫Ωu(x)u 0 (x,y)dΩx.Astherightsideappearstheregionalintegrationinclu…