无网格局部Petrov-Galerkin方法在弹塑性断裂力学问题中的应用

采用无网格局部Petroy-Galerkin方法来分析弹塑性断裂力学问题.这种无网格方法采用移动最小二乘法(MLS)来构造近似试函数和采用Heaviside函数作为加权残值法中的权函数,由于近似函数不满足KroneckerDelta条件,因此采用直接插值法来施加本质边界条件.如果不考虑体力,所形成的整体刚度矩阵只包含局部边界积分,而不包含局部域积分和奇异积分.采用增量Newton-Raphson迭代法来求解弹塑性增量形式的局部Petrov-Galerkin方程.数值算例结果表明,该文方法对于弹塑性断裂力学问题的求解是可行的和有效的,并且所得到的结果具有较好的精度.     

非均质中厚板的无网格LRPIM动力学分析

用局部加权残值法建立了非均质中厚板的局部径向点插值离散系统方程,采用无网格局部径向点插值法分析了非均质中厚板的自由振动和强迫振动问题。用径向基函数耦合多项式基函数来近似试函数,用四次样条函数做为加权残值法中的权函数。所构造的形函数具有Kronecker delta性质,可以很方便地施加本质边界条件。该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行。在计算过程中,取积分中的高斯点的材料参数来模拟问题域材料特性的变化。计算结果表明,利用该方法计算非均质中厚板的自由振动和强迫振动问题可以得到具有较高精度的解。     

各向异性势问题充要的随机边界积分方程

本文针对各向异性势问题提出了一类充分必要的随机边界积分方程。数值计算结果表明在退化尺度附近，充要的随机边界积分方程较习用的随机边界积分方程有较大的优越性。     

各向异性势问题充要的随机边界积分方程

本文针对各向异性势问题提出了一类充分必要的随机边界积分方程，数值计算结果表明在退化尺度附近，充要的随机边界积分程较习用的随机边界积分方程有较大的优越性。     

复杂载荷三维裂纹分析双重边界元法

提出可用于高温、高转速状态下的热动力机械三维含裂构件热弹性分析方法——双重边界元法.首先建立了考虑温度及离心载荷的双重边界积分方程组,并对边界积分方程组的选取及适用范围进行了讨论。然后提出角非快调元模型离散技术。接着提出超奇异积分方程分析去除奇异性方法及数值积分技术.数值算例表明计算结果与有关权函数解十分吻合,说明了用双重边界元法计算复杂载荷条件下三维应力强度因子的有效性.还讨论了有关热应力强度因子权函数解的适用范围.     

《无奇异边界元法》内容简介

孙焕纯等著《无奇异边界元法》一书共有上下两篇 ,上篇阐述虚边界元法的理论、方法与应用。虚边界法有三种 :一般配点法 ;最小二乘配点法 (超额配点法 ) ;最小二乘二重积分法。*分别对弹性空间、弹性平面、薄板、薄壳问题给出了一个从弹性空间方程出发的统一的数值解法 ,抛弃了板、壳理论关于变形和应力的一切假设 ,又对位势问题、弹性平面问题等给出了边界积分方程离散化求解的系数阵元素的解析计算式。下篇针对传统边界元直接法与间接法的边界积分方程的充要性问题进行了论述 ,并对位势、弹性平面和薄板等问题建立了充要积分方程 ;其次是…     

平面弹性力学充要的随机边界元

将平面弹性力学确定性的充分必要的边界积分方程推广到含材料常数随机的不确定问题中去，给出了位移的均值以及偏差的充分必要的边界积分方程。数值计算结果表明，和确定性的积分方程一样，习用的随机边界积分方程在退化尺度附近，无论是均值还是偏差都存在巨大的误差，而充要的随机边界积分方程则始终保持良好的精度     

扁球壳的边界元几何非线性分析

本文从壳体位移的三个微分方程出发，采用付立叶积分变换的基本解，利用加权残值法推导了几何非线性边界积分方程。这种基本解的壳体边界元法类似于板的非线性边界元法，各种变量物理意义明确，能方便地处理各种复杂边界条件及有开口情况。文末算例说明本文方法的可行性、收敛性和精确性，并与二变量边界单元法或有限元结果相比较，吻合较好。     

扁球壳的边界元几何非线性分析

本文从壳体位移的三个微分方程出发，采用付立叶积分变换的基本解，利用加权残值法推导了几何非线性边界积分方程。这种基本解的壳体边界元法类似于板的非线性边界元法，各种变量物理意义明确，能方便地处理各种复杂边界条件及有开口情况。文末算例说明本文方法的可行性、收敛性和精确性，并与二变量边界单元法或有限元结果相比较，吻合较好。     

不连续Reissner矩形板的自由振动

用分区加权残值法研究Ressner矩形板在几何形状,边界条件等有突变时的自由振动问题,将研究对象按照结构几何形状和边界条件的具体情况划分为若干子域,在每个子域内用不同的试函数代入该域的内的控制方程到内部残值,并代入板的边界条件和各子域的协调条件得到边界残值和连续性残值,然后用最小二乘法消除残值,得到特征方程,文中讨论了该方法的收敛性和计算精度,求解了开孔矩形板的固有频率,并与已有结果进行了比较,结果表明：？该方法收敛性好,精度较高,适用范围广。     

调和函数边界积分方程的充要条件

本文以调和函数的边值问题为例,探讨了边界积分方程的充要条件.文中首次提出了超定问题的概念,并建立了超定问题有解的一个充要条件,它也就是直接变量边界积分方程的一个充要条件.文中首次阐明了边界积分方程与变分原理的内在的联系,还指出了间接变量与直接变量两类边界积分方程之间存在着一一对应的关系.文中的慨念、思路和论点不难用于其它有变分原理的问题的边界积分方程.     

扁球壳轴对称非线性弯曲和稳定的权余解

本文从扁球壳的积分方程组出发,通过新定义的残差表达式,用权余法详细地研究了扁球壳轴对称非线性弯曲和稳定问题.通过数值计算可以看出,本方法应用方便,精确可靠.     

加权残数法中试函数选择问题的探讨

本文应用泛函分析和复变函数逼近论的方法初步探讨了加权残数法中试函数选择的原则,并对完备性进行了较深入的讨论.在此基础上,针对组合结构问题,提出了“基因子”概念,并推演出了一些“基因子”及其试函数的一般组构形式.最后,以四边固支矩形加筋板为例,进行了比较计算,取得了积极的效果.     

Necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions

A necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions is established in the present paper. A super-determined problem of harmonic functions is proposed for the first time. Then a necessary and sufficient condition for the existence of solution of the super-determined problem is proved. At the same time, it is a necessary and sufficient condition for the correct formulation of boundary integral equations with direct unknown quantities. A relation between boundary integral equations and variational principles is discovered for the first time. And a one-to-one correspondence between boundary integral equations with direct and indirect unknown quantities is indicated. The concept and route of the present paper can be applied to other boundary value problems possessing variational principles.     

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

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

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.     

A boundary integral equation formulation for the Reissner's plates resting on two-parameter foundation

In this paper a boundary integral equation formulation for the Reissner's plates resting on a two-parameter foundation is established. With the aid of the Hormander Operator method, the equations of the corresponding fundamental solutions are converted into a sixth order partial differential equation with a scale function as an unknown. In order to reduce the equation further, two auxiliary functions are introduced. They satisfy a second and a fourth order equation respectively. The expressions of the auxiliary functions can be derived easily. The fundamental solutions of the Reissnei's plates on the two-parameter foundation arc expressed by a linear combination of the auxiliary functions and their derivatives. The boundary integral equations are formulated by the use of the weighted residual procedure. The fundamental solutions obtained are taken as the kernel functions of the boundary integral equations. A few examples are studied. The numerical results show high accuacy and efficiency of the present formulation.This work was supported by the National Natural Science Foundation of China.     

比例边界等几何分析方法Ⅰ:波导本征问题

提出比例边界等几何方法 (scaled boundary isogeometric analysis, SBIGA), 并用以求解波导本征值问题. 在比例边界等几何坐标变换的基础上, 利用加权余量法将控制偏微分方程进行离散处理, 半弱化为关于边界控制点变量的二阶常微分方程, 即 TE 波或 TM 波波导的比例边界等几何分析的频域方程以及波导动刚度方程, 同时利用连分式求解波导动刚度矩阵. 通过引入辅助变量进一步得出波导本征方程. 该方法只需在求解域的边界上进行等几何离散, 使问题降低一维, 计算工作量大为节约, 并且由于边界的等几何离散, 使得解的精度更高, 进一步节省求解自由度. 以矩形和 L 形波导的本征问题分析为例, 通过与解析解和其他数值方法比较, 结果表明该方法具有精度高、计算工作量小的优点.