矩形域内裂纹对SH波散射问题的边界元分析

运用边界元方法研究矩形域各向同性弹性体含裂纹对ＳＨ波散射问题，利用分离变量求出方程在方程在矩形域的特征值，在正交基上将Ｄｉｒａｃ函数分解，并和基本函数作比较，得出格林函数的显式，导出以裂纹张开位移为未知娄物积分方程，计算出散射位移，算例结果表明，利用边界元法解答有限域的散射问题具有良好的精度且计算速度快。     

夹层杆系结构的位移计算

对于夹层杆件的微元件进行了受力分析，建立了位移基本方程，并对各种载荷作用下的杆系结构进行了位移计算方法的讨论。     

矩形域内裂纹对SH波散射问题的边界元分析

运用边界元方法研究矩形域各向同性弹性体含裂纹对ＳＨ波散射问题，利用分离变量求出方程在矩形域的特征值，在正交基上将Ｄｉｒａｃ函数分解，并和基本函数作比较，得出格林函数的显式．导出以裂纹张开位移为未知数的积分方程，计算出散射位移．算例结果表明，利用边界元法解答有限域的散射问题具有良好的精度且计算速度快．     

空腹夹层板的拟夹层板分析法

本文采用拟夹层板分析法来建立空腹夹层板的基本方程，并引入二个位移函数Ｆ和ｆ来表示三个广义位移。由于假定了计算中面的存在，将基本方程简化的关于Ｆ和ｆ的联立方程组，使问题得到圆满解决。     

基于μ-U方程的饱和土壤固结分析

在给出饱和土壤介质固结有限元分析基本方程的基础上，讨论了饱和土 壤固结分析问题。采取μ－Ｕ形式方程对饱和土固结问题进行分析。利用两点 插值格式与加权余量法，建立了以土骨架位移与孔隙水位移为未知量的固结 问题求解方程；由此克服了以往以μ－ｐ形式方程求解该问题的一些不足。已 在微机程序ＤＩＡＳＳ中实现。文末附有算例。     

层状地基上梁的边界元边界元耦合解法

分别对地基接触面和梁进行离散,假定地基反力的分布情况,并确定梁单元节点和反力未知量;将无限长EulerBernoulli梁的基本解作为梁边界单元法的核函数,然后把Euler-Bernoulli梁边界积分方程应用到各节点,建立起基础梁的边界积分方程组;将层状地基的基本解作为地基边界积分方程的核函数,通过边界积分方程建立起梁各节点竖向位移与地基反力未知量的沉降-反力柔度矩阵;最后,根据地基与梁接触面的位移协调条件,建立起层状地基与EulerBernoulli梁共同作用问题总的边界元-边界元耦合方程组.根据该理论,编制了相应的程序,通过与现有文献对比验证该理论的正确性,并分析了分层地基特性对基础梁的影响.研究结果表明:相比有限元-边界元耦合法,边界元-边界元耦合法的效率更高.     

多锥结合壳流固耦合稳态动力响应分析

采用结构超级元和流体元分析了流场中多锥结合壳受简谐振力作用下的动力响应，利用广义变分原理推导流固耦合系统运动方程。用模态综合法求解系统稳态动位移和动应力，对算例进行了讨论。     

真空预压地基最终应力和位移的解析解

基于非固结渗流问题的基本方程,导出真空预压地基最终应力和位移的解析解,分析地基在渗流体力和压差面力联合作用下的变形规律,认为地基竖向与水平向渗透系数之比ξ对真空预压的效果有重要影响．ξ值越大,真空度、应力沿深度衰减得越慢,竖向沉降和向内侧的水平位移越大,因此设置竖向排水体可有效增加加固深度,提高加固效果．     

弹性力学问题的第二类边界积分方程——边界元法及其二维基本解

本文导出了以单位集中不连续位移（第二类边界）产生的场为基本解的两种边界积 分方程——直接法和间接法边界积分方程，按 Ｈａｄａｍａｒｄ主值存在，进行高阶奇异积 分并给出了二维问题的两个基本解，与通常的第一类直接法和间接法相比，第二类方法 对处理裂纹等不连续结构面较为有利。     

直接边界元法及其在弹性力学问题中的应用

革通过弹性力学问题的基本解将域内微分方程变换成边界上的积分方程，然后在边界上离散；由已知边界位移和边界应力直接求出未知边界位移和边界应力，并得出据以计算整个问题域的位移场和应力场。最后运用此方法求解一个弹性力学问题并与有限无法的计算结果进行了比较。     

波动方程的无网格-精细积分方法

将无网格法和精细积分用于波动方程的计算.在空间上用无网格法进行离散,用修正变分原理处理本征边界条件;在时间域上用精细积分法求解动力学方程,然后给出两个波动方程的算例.数值结果表明此方法是稳定、精确的.     

基于径向积分边界单元法的内场声学分析

由于Helmholtz方程的基本解是频率的函数,传统边界单元法在处理声场特征值问题时具有天生的缺陷.采用Laplace方程基本解生成积分方程,通过径向积分法将在此过程中产生的域积分项转化为边界积分.此方法克服了传统边界单元法系数矩阵对频率的依赖,同时克服了特解积分法和双重互易法对特解的依赖,将内场声学特征值问题转化为广义特征值问题.最后通过内场声辐射分析和声学特征值分析验证了算法的有效性.     

两类发展方程广义差分法的误差估计

讨论二维区域上两类数学物理方程一次元格式的广义差分法。关于双曲型积分微分方程和Sobolev方程,证明了最优H1,L2和最大模误差估计,其收敛阶与线性有限元方法一致。此外,还获得了近似解的超收敛结果。     

边界单元法三维应力分析计算程序

以三维弹性力学问题的基本解 (Kelvin解 )为基础开发了边界单元法三维弹性应力分析计算程序 ,并对其进行了验证。结果表明 ,该程序可用来求解三维弹性应力问题 ,尤其适用于三维应力集中问题。程序中采用动态分配内存 ,自动选点积分 ,并利用分块解法求解方程组 ,有效地节省了计算机资源 ,扩大了求解问题的规模。     

An element by element spectral element method for elastic wave modeling

The spectral element method which combines the advantages of spectral method with those of finite element method, provides an efficient tool in simulating elastic wave equation in complex medium. Based on weak form of elastodynamic equations, mathematical formulations for Legendre spectral element method are presented. The wave field on an element is discretized using high?order Lagrange interpolation, and integration over the element is accomplished based upon the Gauss?Lobatto?Legendre integration rule. This results in a diagonal mass matrix which leads to a greatly simplified algorithm. In addition, the element by element technique is introduced in our method to reduce the memory sizes and improve the computation efficiency. Finally, some numerical examples are presented to demonstrate the spectral accuracy and the efficiency. Because of combinations of the finite element scheme and spectral algorithms, this method can be used for complex models, including free surface boundaries and strong heterogeneity.     

二维Laplace方程Dirichlet问题直接边界积分方程的Galerkin解法

用Green公式和基本解推导得出的直接边界积分方程来求解二维Laplace方程的Dirichlet问题.对直接边界积分方程大都采用配点法求解,还未见有实际用Galerkin边界元来解的报道.对Laplace方程的直接边界积分方程进行变分后,利用Galerkin方法,同时采用线性单元变分对方程进行了求解.该方法需要在边界上计算重积分,推出了第一重积分的解析计算公式,对无奇异性的外层积分则采用高斯数值积分.数值实验表明该方法是可行有效的.     

多股电缆线电场的边界元方法

用边界元方法．研究了多股电缆线电场问题．在控制方程的基础上，利用其基本解，通过加权余量法建立内点和边界的积分方程．对边界剖分进行二次元处理．得到离散化的计算格式．给出了多股并行电缆线不同轴心距下所形成的电场特性．     

采用磁致伸缩模型进行磁弹性有限元分析的基本理论

基于磁致伸缩模型，推导了用于磁弹性有限元分析的基本公式和理论，该模型有易于理解，推导简便，便于有限元实现的特点。     

弹性力学基本方程弱形式

弹性力学基本方程正确的弱形式将是有限元分析的基础.如直接从基本方程出发,由于其整个方程可以有正负号的差异,往往得不到正确的弱形式.因此从泛函分析的角度出发,基于共轭空间的概念和泛函分析的基本定理准确地给出了弹性力学基本方程的弱形式;给出了连续介质在位移或物理常数间断面上的条件.将三维空间的弹性力学动力学方程,理解为定义在四维空间域上的运动方程,导出了弹性力学动力学方程的弱形式,在此基础上推导出了与Hamilton变分原理同样的结果.     

二阶变系数常微分方程的恰当因子解法

针对二阶变系数常微分方程求解的两个基本方法 ,给出了求对应的齐次方程特解的方法和求恰当因子的方法