R-函数理论及准Green函数方法在正交各向异性板振动问题中的应用

首次将R-函数理论及准Green函数方法应用于求解固支正交各向异性薄板的自由振动问题。首先引入参数变换,将正交各向异性薄板的自由振动微分方程转化为双调和算子的边值问题,并应用R-函数理论,以解析函数形式描述复杂边界形状;利用问题的基本解和边界方程构造了一个准Green函数,该函数满足了问题的齐次边界条件;通过R-函数理论构造适当的边界方程,消除了积分方程核的奇异性;再采用Green公式将其化为第二类Fredholm积分方程。数值算例表明:该方法减少了理论计算量,精度较高。本文还证明了其优越性和正确性,是一种新型的数学方法。     

含角应力奇异扇形板的振动特性

将富里叶-贝塞尔级数引入积分方程方法，给出了一种研究含角应力奇异直边简支扇形板结构振动特性简捷，高效的积分方程方法。数值计算结果表明该方法不仅具有较高的精度和稳定性，而且还为分析更为复杂的扇形板结构的振动问题提供了可靠的前提。     

用边缘荷载格林函数方法计算含裂纹Reissner型有限板的弯曲

本文在文献[1]所得到的受边缘荷载格林函数基本解方法基础上,利用叠加原理,通过边界积分方程的方法,分析了含裂纹Reissner型板的弯曲断裂问题。计算表明方法正确,便于应用。     

一种求解含振子及弹性支承圆板振动的新方法

将富里叶－贝塞尔级数引入积分方程[1],推导出一种研究含振子及弹性支承圆板振动特性的新方法,根据积分方程和富里叶－贝塞尔级数理论,首先用第一类贝塞尔函数构造圆板的格林函数,然后由叠加原理将圆板的自由振动问题转化为积分方程的特征值问题;进面将积分方程形式的特征值问题转化为无穷阶矩阵的标准特征值问题,计算时根据精度的要求,截取无穷阶矩阵的标准特征值为有限阶矩阵的标准特征值问题,采用Q－R算法,计算实践表明,本方法不仅具有运算简捷,精度高,适用性强的特点,而且能从整体上对系统的动态性加以研究,从而为这类系统的优化设计提供有;力的 工具。     

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

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

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

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

矩形薄板弯曲问题的二维广义有限积分变换法

运用二维广义有限积分变换解法，本文推导出不同边界条件下矩形薄板弯曲问题的解析解．在推导过程中，选取满足边界条件的梁振型函数为广义积分变换的积分核，由此构造出广义有限积分变换对，通过对薄板弯曲问题的控制方程进行二维广义积分变换，可以将控制方程转换为易于求解的线性代数方程组．该方法无需预先选取位移函数，无需进行繁琐的叠加过程，求解过程思路清晰，说明该方法更加正确合理．最后通过计算实例对比，验证了该方法的合理性及所推导公式的正确性．     

和界面接触的刚性线夹杂对SH波的散射

利用积分变换方法，得出了两相材料中单位简谐力的格林函数。根据简谐集中力的格林函数得出了和界面接触的刚性线的散射场。利用无穷积分的性质，把和界面接触刚性线的散射场分解为奇异部分和有界部分。通过分解后的散射场建立了和界面接触剐性线在SH波作用下的Cauchy型奇异积分方程。根据所得奇异积分方程和刚性线的散射场得到了刚性线端点的奇异性阶数及奇性应力。应用刚性线端点的奇性应力定义了刚性线端点的应力奇异因子。对所得Cauchy型奇异积分方程的数值求解，可得刚性线端点的应力奇异因子。     

寻求含孔洞薄板弯曲基本解的一般方法

1 引言用边界元法求解含孔薄板弯曲问题时,若采用通常的基本解作为积分方程的核函数,孔洞也必须作为离散的边界来处理.如果找到一种基本解,使之对孔洞边界条件自然满足,就可以避免沿孔洞边界的积分,而只剩下沿边界的积分,使数值处理大为简化.为此,本文运用的复变函数方法,给出了寻求含自由孔洞薄板弯曲基本解的一般方法.这一工作在理论上和应用上都是很有意义的.2 基本方程不含孔洞无限大薄板弯曲的基本解为v~0(P,Q)=(1/(8πD))r~2lnr (1)表示无限大板中P 点作用一单位力在Q 点引起的挠度,r 数与Q 两点间的距离,D 是弯     

层状地基任意形状刚性基础动力响应求解

提出了基于积分变换、对偶方程与精细积分算法求解多层地基任意形状刚性基础的动力刚度问题. 首先在频率波数域内圆柱坐标体系中利用圆形微元的对称与反对称特性建立多层地基中格林影响函数的波动方程,然后将应力和位移关系表示成对偶形式进行精细积分求解以提高计算精度和稳定性. 再将任意形状刚性基础与地基的交界面离散化为一系列圆形微元,利用格林影响函数建立其平动与转动动力刚度的矩阵方程. 该求解方法高效、准确并且计算稳定,适于任意复杂多层地基任意形状基础动力刚度的计算.      

加权残数配点法解正交各向异性板的积分方程

本文推导了一般各向异性板弯曲的积分方程,运用加权残数配点法求解了正交各向异性板弯曲的积发方程,本文将部分配点取在边界上,另一部分配点取在域外,只用关于找度的基本积分方程,而不用关于转角的补充积分方程,简化了方程求解和计算程序,由于正交各向异性板没有争析形式的、实用的基本解,本文提出了两种新的近似基本解;加权双三角级数;广义各向同性板解析形式的基本解和加权双三角级数的叠加,算例表明,本文提出的解法和近似基本解适用于各类边界条件的正交各向异性板,具有简单、可靠、精度高等优点。     

A new boundary integral equation for plane harmonic functions

In this paper, we propose a new boundary integral equation for plane harmonic functions. As a new approach, the equation is derived from the conservation integrals. Every variable in the integral equation has direct engineering interest. When this integral equation is applied to the Dirichlet problem, one will get an integral equation of the second kind, so that the algebraic equation system in the boundary element method has diagonal dominance. Finally, this equation is applied to elastic torsion problems of cylinders of different sections, and satisfactary numerical results are obtained.     

非均匀弹性材料中反平面的运动裂纹

用解析方法研究了非均匀弹性材料中反平面运动裂纹问题。首先采用余弦变换求解非均匀材料的基本方程，然后根据混合边值条件建立裂纹运动的对偶积分方程，再把对偶积分方程化为第二类Fredholm积分方程。给出了数值算例，计算结果表明材料的非均匀性对动应力强度因子有较大的影响。     

New integral equation for plane elasticity crack problems

In the usual formulation of singular equation approach for crack problems in plane elasticity [1,2], if one changes the right-hand term of the integral equation from tractions to resultant forces, a new integral equation can be obtained and is presented in this paper. The newly obtained integral equation has a log singular kernel. Interpolation equation for the dislocation functions (the undetermined functions in the integral equations) is proposed. Numerical examination is used to demonstrate the efficiency of the present technique, and a number of numerical examples are given.     

圆板弹塑性弯曲的简单样条积分方程法

提出了圆板弹塑性弯曲的简单样条积分方程法.以径向转角作为未知量建立积分方程并结合B样条函数进行求解.这一方法具有域积分容易处理、精度高和计算简单的优点.计算结果表明本文解与文[2]解吻合良好.     

Convolution of the impact three-dimensional elasto-dynamics and dynamic stress intensity factor for an elliptic cracks

This paper presents a formulation for three-dimensional elastodynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack is determined by solving a Fredholm integral equation of the first kind. The results of this paper are very close to those given by the two-dimensional dual integral equation method. The project supported by the National Natural Science Foundation of China (K19672007)     

裂纹端部细短纤维的应力分析

基于裂纹端部存在与其裂纹面相垂直的二相细短纤维分析模型，采用叠加原理推导了求解纤维表面应力分布函数的积分方程，通过简化得到了该方程的解析表达显式，该积分方程的特征值方程是纤维几何参数，材料常数以及纤维相对于裂纹位置的相关函数，当材料参数不满足特征方程时，积分方程将具有唯一解；并借助数值方法，给出了纤维剪应力分布算例，和纤维对应力强度因子的影响。     

非均匀介质散射问题的体积分方程数值解法

将非均匀介质视为某一均匀背景介质中的扰动，可建立用均匀背景介质格林函数作基本解的体积分方程．给出了配置法求解体积分方程的数值方法，首先解得扰动域内各点以速度扰动为权的波场函数，然后回代计算得到观测面上各接收点的散射波场．与边界元法和Ｂｏｒｎ近似法计算结果比较表明该方法具有很高的精度，可得到穿过非     

Dynamic response of elastic layer on stiff foundation under time harmonic surface vertical concentrated load

In this paper, the line-load integral equation method proposed in reference [1] is first used for solving the elastodynamic problems. A set of one-dimensional regular integral equation is derived for calculating the dynamic response of elastic layer on stiff foundation under time harmonic surface vertical concentrated load. And the numerical solution of the integral equation is obtained.     

An analysis on entrance region effect of the laminar radial flow between two parallel disks

In this paper, B. B. Golubef method^[1] is used for calculating the radial diffuse flow between two parallel disks for the first step. The momentum integral equation together with the energy integral equation is derived from the boundary layer momentum equation, and the expression of secondary approximation explicit function in which the channel length of entrance region varies with the boundary layer thickness can be obtained by using Picard iteration[2] in the solution of the energy integral equation. Therefore, this has made it possible to analyze directly and analytically the coefficients of the entrance region effect. In particular, when the outer diameter of disk is smaller than the entrance region length, the advantage of this method can be prominently manifest.Only because the energy integral equation is employed, the terms in the pressure loss coefficient can be independently derived theoretically. The computable value of the pressure loss coefficient presented in this paper is nearer to the testing value than that in ref. [3] when the entrance correction Reynolds number Re<100. Therefore the results in this paper within Re<100 are both reliable and simple.