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

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

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

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

胶合搭接连接的界面应力

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

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

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

平面斜边界问题的边界元分析

本文对处理斜边界问题的边界元法的方程进行了推导和分析,通过引入局部坐标系,建立局部坐标与总坐标间的变换关系,使边界元法在处理结构斜边界问题中更为简捷方便。通过几个算例及某发动机叶片燕尾形榫头的应力分析,证明结果是令人满意的,尤其在计算边界应力时,边界元法比有限元法更为精确。     

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

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

三维疲劳裂纹扩展仿真的双重边界元法

运用双重边界元方法求解三维疲劳裂纹前沿的应力应变场,基于Forman理论、最小应变能密度法和Elber模型,采用增量步裂纹扩展分析方法,并根据裂纹几何形状的改变,对裂纹面进行网格重划分和迭代计算,模拟了三维裂纹的扩展和预测裂纹扩展寿命.获得了裂纹前沿各点的扩展长度、扩展方向和应力强度因子等特征量.扭力轴表面裂纹扩展的应用实例表明该方法正确合理.     

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

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

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

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

含共线分布多裂纹板的剩余强度

采用各向异性体平面弹性理论中的复势方法，应用保角映射技术和Faber级数展开，导出在任意载荷作用下多裂纹板应力场的级数解，引入当量屈服应力修正裂尖塑性区，并利用Swift韧带屈服准则建立舍共线分布多裂纹结构的剩余强度分析模型，计算结果与试验结果吻合较好。文中还讨论了各裂纹参数对结构剩余强度的影响。计算结果表明，本文以当量屈服应力代替材料屈服应力，避免了复杂的弹塑性分析，思路简单有效。便于工程应用。     

一种计算复合材料层合板面外应力的新型线性分析方法

复合材料层合板壳结构的面外应力对于层间分层的起始和扩展具有重要意义。首先简要地给出了经典的非协调位移有限元列式,然后基于修正的H-R变分原理导出了单元面外应力变量与位移变量的关系。面外应力的有限元线性系统为面外应力边界条件的引入提供了方便。同时,通过该线性系统可直接获得节点上唯一且连续的面外应力结果。与非协调位移有限元的结果比较,本文的方法可明显地提高数值结果的精度,且计算效率及资源要求均优于非协调辛元法和非协调广义混合元法。     

环形激波聚焦流场特性的间断有限元方法(英文)

针对环形激波聚焦时产生的高温、高压特性,采用间断有限元方法对环形激波在同轴圆柱形激波管内的聚焦流场进行了数值模拟。为了验证采用方法的有效性,首先对经典激波管问题进行了求解,计算结果与解析解吻合很好。其次,对入射激波马赫数为2.0到5.0的环形激波聚焦流场进行了计算,获得了环形激波在圆柱形激波管内聚焦过程的气动特性。模拟结果表明,采用间断有限元方法能够有效地捕捉激波聚焦过程形成的绕射、反射和聚焦等主要流动特征,而且在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力。最后比较了网格尺度对聚焦点峰值压力的影响。     

高精度微分求积曲梁单元的建立与应用

首先由能量原理导出曲梁弯曲问题的控制微分方程，在此基础上应用微分求积法原理分别给出了曲梁内点和端点的微分求积方程，由此形成曲梁单元的刚度方程，从而建立了微分求积曲梁单元，并给出了曲梁结构刚度方程的边界条件。通过算例分析，得到了微分求积单元法结构离散时应使单元数量少的原则和求解精度与单元长度比基本无关的性质。与有限元方法的结果比较表明，本文导出的曲梁单元是一种具有很高求解精度的单元。》     

An Efficient Method for Local Buckling Analysis of Stiffened Panels

The local buckling of stiffened panels is one of possible failure modes and concerned by engineers in the preliminary design of lightweight structures. In practice,a simplified model,i.e.,a rectangular plate with elastically restrained along its unloaded edges,is established and the Ritz method is usually employed for solutions. To use the Ritz method,however,the loaded edges of the plate are usually assumed to be simply supported. An empirical correction factor has to be used to account for clamped loaded edges. Here,a simple and efficient method,called the quadrature element method(QEM),is presented for obtaining accurate buckling behavior of rectangular plates with any combinations of boundary conditions, including the elastically restrained conditions. Different from the conventional high order finite element method(FEM),non-uniformly distributed nodes are used,and thus the method can achieve an exponential rate of convergence. Formulations are worked out in detail. A computer program is developed. Improvement of solution accuracy can be easily achieved by changing the number of element nodes in the computer program. Several numerical examples are given. Results are compared with either existing solutions or finite element data for verifications. It is shown that high solution accuracy is achieved. In addition,the proposed method and developed computer program can allow quick analysis of local buckling of stiffened panels and thus is suitable for optimization routines in the preliminary design stage.     

一种适合于分析复合材料层板的三维杂交应力超单元研究

本文基于三维弹性理论和Hellinger-Reissner变分原理,推导了一种适合于分析复合材料层板的三维杂交应力超单元。该单元不同于常规的三维单元,它将复合材料层板沿厚度方向划分成一个或几个单元,每个单元中所包含的层数可视精度和计算费用的要求而定。在应用应力杂交模型导出三维超单元时,本文首次提出一种场变量假设方法,即在超单元内独立假设位移场,而应力场则分层独立假设,因而单元能较好地模拟层合板的力学性能。考虑横向剪切变形和局部扭曲变形的影响,保证了单元分析的精度。同时由于单元的自由度数不随层数的增加而增加,避免了常规三维单元在分析复合材料层板时自由度数随层数的增加而增加的弊端。通过对几个典型的复合材料层板的静、动力分析及与试验的比较。表明了单元能有效地用于各种层板的分析。     

随机有限元方程一般式的两种推导方法

８０年代初发展起来的随限元法是解决随机问题强有力的数值分析工具，如随机动力学、随机场及结构强度的可靠性等等问题，本文从有限元基本方程出发，分别采用对有限元方程按Ｔａｙｌｏｒ级数展开和对有限元方程求偏导数等方法。推导了随机有限元基本方程的一般式，并用数学归纳法证明了上述两种方法推导结果的与等效性，本文研究表明，Ｔａｙｌｏｒ级数展开法推导，数学意义明确，但推导过程较为复杂；而按求偏导数法，数学意义不够     

基于超分辨率重构方法的湍流流场重构

从低分辨率流场数据中获取精细流场信息具有重要的研究意义。基于卷积神经网络的超分辨率重构方法是近年来发展的一种较为有效的精细流场重构方法。本文采用高效亚像素卷积神经网络（Efficient Sub-Pixel Convolutional Neural Network,ESPCN）,对Rayleigh–Bénard（RB）对流的数值模拟数据和湍流边界层（Turbulent Boundary Layer,TBL）的实验测量数据进行了超分辨率重构,并与双三次插值方法（Bicubic Interpolation）的重构结果进行对比。对比结果表明:在较小的下采样比下,ESPCN方法和Bicubic方法的重构精度相当;在较大的下采样比下,ESPCN方法的重构精度明显优于Bicubic方法。此外,ESPCN方法对数据梯度较大区域的超分辨率重构效果优于Bicubic方法。