含孔边裂纹各向异性有限板的二维问题

应用边界元法，并结合一种特殊的复变函数基本解，研究了含孔边裂纹的各向异性有限板的应力强度因子修正系数的计算。这种特解边界元法避免了在裂纹附近划分单元，可以获得较好的数值计算精度。通过算例分析，结果表明本方法精度高，运算量小，具有一定的工程应用价值。     

采用解析积分的常单元边界元法

本文采用一种在局部坐标系中的解析式来计算常边界单元及内点应力的方法,推导了常单元的积分解析式,并对若干典型例子进行了分析和讨论。算例结果表明,采用本文的方法,能获得较高的精度,尤其在计算应力集中区的应力及靠近边界的内点应力时,效果更为显著。     

处理动边界问题的无网格/直角网格混合算法

把无网格/直角网格混合算法发展用于涉及动边界的非定常流动计算问题.采用无网格/直角网格混合算法,计算区域整体采用直角网格,只在物面附近嵌入局部无网格.对于动边界问题,本文让主体直角网格保持固定,用原混合算法中的局部无网格区跟随运动边界,使得动边界物面附近有不变的点云结构,避免了点云重构;无网格区与直角网格交界面随边界运动跟踪确定,由于主体基于直角网格,方便了交界面的跟踪.用发展的混合算法,结合非定常Euler方程双时间推进方法,先对翼型作俯仰振荡的非定常跨声速流动进行了数值模拟,计算结果与现有文献和实验数据进行了比较与验证;结合刚体运动方程,成功模拟了二维外挂物投放过程,算例展示出发展的混合算法适合处理大位移运动边界问题.     

三维弹塑性边界元1/r~3积分奇异性的处理

本文将文[2]的基本思想与方法,推广应用于处理三维弹塑性边界元法塑性单元中的1/r~3积分奇异性问题,推导了有关的具体公式,并从理论上证明了该方法处理1/r~3奇异性的有效性。针对奇异性单元。本文提出旋转与坐标转换的分割技术,可大大减少工作量。并使计算机编程更简洁。数值算例研究表明,文中采用的方法是十分有效的。     

胶合搭接连接的界面应力

采用边界元法计算了分析了具有理想界面和非理想界面的双板条胶合搭接连接的界面应力。经计算得到了完整胶合界面时界面应力的分布，得出在胶合界面两端出现了应力集中现象；胶合界面发生开裂后的界面应力分布，理想界面裂纹尖端的复应力强度因子和非理想界面裂纹前洞的应变能密度。     

一种采用解析积分的直线单元的边界元法

本文提出一种采用解析的直线单元计算受有离心力的二维问题的边界元法,推导了相应直线单元的积分解析式。算例表明,本文方法能获得较高的数值精度,尤其在计算靠近边界的内点应力时,效果更好。     

新电位移边界条件下的压电介质断裂问题

应用边界元法，并结合一种特殊的复变函数基本解，研究了压电介质的断裂问题。由于该基本解完全遵循本构方程没有对电位移边界作假设，所以结论同非穿透性裂纹假设的结果有较大差异。通过计算表明，外加力载荷不但会引起裂纹尖端应力奇异性，也会引起电场奇异性，同样，外加电压也会引起裂纹尖端应力奇异性。     

表面开口矩形缺陷兰姆波散射的边界元模拟

本文介绍了边界元法与兰姆波本征模式函数相结合的混合边界元法理论模型,用混合边界元法计算了铝板中A0模式兰姆波在表面开口矩形缺陷上的散射.通过数值计算,给出了不同频率下的模式转换特性和散射系数,并给出板中兰姆波传播和裂纹对兰姆波散射的波场快照,本文的研究工作对板状结构兰姆波无损检测技术的研究具有参考价值.     

基于比例边界有限元的结构强度分析方法研究

航天结构经历复杂严酷的载荷环境，且诸多大传载关键结构存在小尺寸倒角。为有效处理航天结构小尺寸倒角强度评估问题，提高结构应力应变求解精度和计算效率。本文基于比例边界有限元方法对典型小尺寸倒角的航天结构边界进行离散，构建强度计算模型，针对典型结构不同倒角尺寸开展强度分析，得到不同尺寸倒角对局部应力的影响规律，给出局部区域的应力分布，并与光测试验和常规有限元分析结果进行对比。结果表明，与常规有限元方法相比，比例边界有限元方法可以有效减少计算自由度，提高计算效率和计算精度，由于其半解析的特点，可以方便的得到内部域的应力场，为小尺寸倒角航天结构强度评估、寿命预测及优化设计提供有效的分析手段。     

最小二乘和奇异值分解法在边界元法中的应用

采用最小二乘与奇异值分解结合的方法，给出求解系数矩阵不满秩的线性代数方程组的数值方法，进而将此方法应用于边界无法中，处理给定外力的第一边值问题，特别地用于处理给定外力的三维裂纹问题。此外，本文还给出求解三维有限体裂纹问题的超奇异积分方程组，并使用有限部积分与边界元法为其建立了数值法。最后计算了若于典型例子的应力强度因子，数值结果与现有文献的相比，符合很好。     

接触问题数值分析方法的研究现状与发展

接触问题涉及到应力集中、边界非线性、甚至材料或几何非线性，问题很复杂但又十分重要。随着电子计算机的发展，有限元与边界元理论与数值方法也取得了很大进展。目前，这两种方法已成为工程接触问题十分有效的研究和数值分析手段。本文综述了接触问题中的有限元法与边界元法及其工程应用方面所取得的进展，着重围绕小变形弹性接触、小变形弹塑性接触以及大变形接触等方面的内容。此外，文中还对今后研究的发展方向作了展望。二十多年来，在小变形弹性接触问题中的有限元法和边界元法发展的比较成熟，并在工程中获得广泛应用；小变形弹塑性接触问题中的有限元法发展的也比较成熟。然而，它们在如下方面还必须开展进一步的研究，例如，不同材料非线性模式的接触、不同蠕变模式的接触、大变形弹塑性接触、基于实验力学的反迭代分析、以及复杂的工程构件的接触分析等。     

机翼非平行边界层稳定性研究

从Navier-Stocks(N-S)方程导出曲线坐标系下的抛物化稳定性方程(Parabolicstabilityequation,PSE),研究机翼非平行的可压缩边界层稳定性问题。发展了求解PSE的高效数值方法:引进法向变换,使得在临界层与壁面之间的扰动量变化最快的区域有更多法向网格点;采用包含边界邻域在内的完全四阶精度的法向差分格式,这对方程精确离散至关重要;以及全局法和局部法相结合的数值方法及其新的迭代公式,能大大加速收敛并得到更精确的特征值。算例分析研究了扰动增长因子和形状函数等演化曲线。     

双重互易边界元法在求解N—S方程中的应用

提出了用双互易边界元法求解Ｎ－Ｓ方程的新思路。用Ｐｅａｃｅｍａｎ－Ｒａｃｈｆｏｒｄ算子分理解将时间相依的Ｎ－Ｓ方程分解为线性和非线性子问题，线性问题用共轭梯度法解除压强－速度的耦合；非线性问题进行局部线化。对所得的Ｐａｓｓｉｏｎ方程及相关类型方程采用双重互易的边界元法求解，消除了传统Ｎ－Ｓ方程边界元解法的区域积分问题。     

填充各向异性介质二维开口腔体散射特性分析

应用非重叠型区域分解法(DDM)结合有限元法(FEM)和边界元法(BEM)分析了填充多层各向异性介质的二维开口腔体横电波(TE)散射特性。对腔体外的区域采用BEM法分析,将腔体内的每层介质作为一个子域,用FEM法分析,各子域间通过传输条件进行耦合。分别计算了腔体中填充各向同性和各向异性介质时的雷达散射截面,数值结果表明了该方法的有效性。采用这种技术,大大地减少了对计算机内存的需求。     

带压力梯度的非平行流边界层稳定性研究

发展高精度的高阶谱方法,通过抛物化稳定性方程及其局部法,研究在压力梯度作用下的非平行流边界层稳定性问题。与经典的平行流边界层稳定性结果相比,显示了在某些条件下非平行性对稳定性的关键作用。提出了压力梯度对临界雷诺数的影响研究,探讨了不同压力场对流动稳定性作用的规律性。     

一类具有混合边界条件的奇摄动拟线性边值问题

研究一类具有混合边界条件的奇摄动二阶拟线性边值问题,在构造形式渐近解的基础上,利用微分不等式理论证明了解的存在性,并得出了解的任意阶的一致有效展开式。     

含边界和损伤因素的纺织复合材料细观胞元模型

纺织复合材料的细观力学分析通常以胞元模型为基础,胞元边界条件的合理施加是获得精确分析结果的关键之一。本文以平面机织复合材料为例,讨论了细观胞元的选取和周期边界条件的施加方法,在此基础上建立了二维细观有限元模型,通过与全厚度模型分析结果的比较,研究了周期胞元模型的合理性,特别是边界效应和局部损伤等非周期因素对分析结果的影响,给出了这些因素下胞元边界条件的处理方法。