复合材料层板边缘效应的准相似元解法

用准相似单元法研究了复合材料层板直边和孔边的自由边缘效应，得到了一些具有工程意义的结果。准相似单元法是相似单元法的推广，利用这一方法可以求解三维的应力集中问题，避免了求单刚时的重复计算。由于采用了按相似层求解的技术，应用准相似单元法求解应力集中问题时对计算机容量的要求大为降低     

应力高梯度问题的无网格分析

基于移动最小二乘法的无网格计算，采用线性基函数即可得到C^1连续位移场，使得应力，应变场在整个求解域内保持连续；节点之间脱离了单元的约束，对求解域进行离散和加密节点时变得十分灵活，因此适合分析应力高梯度问题。本文简要介绍了无网格方法的基本原理，给出了确定节点影响域大小的方法，应用无网格方法对带有V型缺口的受拉方板J23-10曲柄压力机机身进行了受力分析，得到的应力集中部位的计算结果与实际值更为接近。     

SH波散射与界面圆孔附近的动应力集中

建立了求解含有界面圆孔的二种不同弹性组合介质中ＳＨ波的散射和界面圆孔附近的动应力集中问题的Ｇｒｅｅｎ函数法给出了一个具有半圆形缺口的弹性半空间水平表面上任意一点承受时间谐和的出平面线源荷载作用时位移函数的基本解取基本解作为Ｇｒｅｅｎ函数，建立起问题的定解积分方程最后给出了界面圆孔的动应力集中的算例和结果，并讨论了不同介质参数的组合对动应力集中的影响     

弹塑性力学问题的虚单元法

具有有限差分法特征的虚单元法，可视为是有限元法向任意多边形单元的扩展。在材料细观力学性能表征、非均质材料力学分析等非线性问题方面，传统的弹塑性有限元法具有网格数目多、效率低下等不足之处，而虚单元法使网格划分更加灵活，为材料的弹塑性力学分析等非线性问题提供了新的思路。基于增量法弹塑性力学原理和双线性投影算子，建立了弹塑性力学问题的虚单元法求解技术，提出了弹塑性力学问题虚单元法的应力更新方案，研究了弹性力学问题虚单元法的精度和收敛性，讨论了虚单元法求解弹塑性力学问题的网格依赖性。同时，开展了任意多边形和凹多边形单元的数值试验研究，结果表明，虚单元法无须分割多边形，仅需节点自由度便可求得单元刚度矩阵和应力等效荷载，程序实现简单，计算精度高，改善了传统有限元的网格依赖性和塑性区的网格奇异性。     

无限大板开孔弹性波的散射及动应力集中

采用弹性平板理论及复变函数理论，对含孔无限大平板弹性波的散射及动应力集中问题进行了分析研究，建立了求解平板开孔动应力集中问题的复变函数方法。若同时采用映射变换，就为求解平板开任意形状孔的动应力集中问题提供了一种规范而有效的方法。为说明问题，本文给出了平板开圆孔及椭圆孔动应力集中因子的数值结果。     

自适应一致性高阶无单元伽辽金法

近来提出的一致性高阶无单元伽辽金法通过导数修正技术大幅度减少了所需积分点数目,并能够精确地通过线性和二次分片试验,显著改善标准无单元伽辽金法的计算效率、精度和收敛性.本文在此基础之上,充分利用无单元法易于在局部区域添加节点的优势,发展了一致性高阶无单元伽辽金法的h型自适应分析方法.根据应变能密度梯度该方法自适应地确定需节点加密的区域,基于背景积分网格的局部多层细化要求生成新的计算节点,同时考虑了节点分布由密到疏渐进过渡的情形.采用相邻两次计算的应变能的相对误差作为自适应过程的停止准则,将所发展自适应无网格法应用于由几何外形、边界外载和体力等因素造成的应力集中问题的计算分析.数值结果表明,所发展方法能够自适应地对高应力梯度区域进行节点加密,自动给出合理的计算节点分布.与已有的标准无网格法的自适应分析相比,所发展方法在计算效率、精度和应力场光滑性等方面均展现出显著优势.与采用节点均匀分布的一致性高阶无单元伽辽金法相比,它大幅度地减少了计算节点数目,有效提高了一致性高阶无单元伽辽金法在分析应力集中等存在局部高梯度问题时的计算效率和求解精度.     

SH波入射时半空间界面裂纹与圆形衬砌的相互作用

利用复变函数和Green函数法研究了双相介质半空间界面裂纹及界面附近圆形衬砌对SH 波的 散射与动应力集中问题。首先,采用映像思想构造满足自由边界条件的散射波表达式,进而求解所需的 Green函数;其次,采用裂纹切割技术构造裂纹,并根据连续性条件建立了求解该问题的无穷代数方程组; 最后,给出了不同入射波数时界面裂纹与衬砌的相互作用。结果表明,裂纹的存在显著放大了衬砌界面的动 应力集中。     

SH波作用下界面任意形状孔洞附近的动应力集中

采用Green函数和复变函数法求解了平面SH波在界面任意形状孔洞上的散射问题.首先,取含有任意形状凹陷的弹性半空间,在其水平表面上任意一点承受时间谐和的反平面线源荷载作用时的位移场作为Green函数.然后,按契合方式构造出界面任意形状孔洞对SH波的散射模型,利用所得Green函数按界面位移连续条件建立求解问题的定解积分方程组,求解界面孔附近的动应力集中系数.最后,给出了界面上椭圆孔和方孔边缘动应力集中系数的数值结果,并讨论了不同介质参数和孔洞形状对孔附近动应力集中系数的影响.     

弹性力学的杂交自然单元法

为了解决自然单元法无法直接求解节点应力以及应力解精度不高的问题, 将应力杂交的思想引入自然单元法中, 与弹性问题的Hellinger-Reissner变分原理结合, 提出了弹性问题的杂交自然单元法. 算例表明: 杂交自然单元法的计算结果与解析解吻合, 证明该方法是可行的; 在求解应力方面, 杂交自然单元法比自然单元法有着更高的计算精度, 而且可以直接求解出节点的应力.     

重力坝问题的非结构分层四边形高阶亚参元分析

首先详细论述了一种通过增加棱内"虚节点"和面内"虚节点"方式得到的四边形升阶谱单元及相应的阶谱函数。然后,结合矩阵的一维稀疏存储技术(CSR格式),设计了一种分层四边形高阶亚参元方法,并编制了相应的FORTRAN语言程序进行了实现。该方法具有很好的承袭性,即在生成高阶单元特性矩阵时始终能够承袭低阶单元的特性矩阵,这样可减少计算量,大大提高了编程和计算效率。最后,通过将所设计的方法应用于重力坝/腹拱坝问题的求解,验证了本文方法具有与相应常规有限元方法一样高的计算精度,且当使用分层Q_8元和分层Q_(12)元时,随着网格的加密,坝踵及坝趾区产生应力集中的区域明显缩小,坝踵-坝趾线段上的第一主应力和y方向应力也基本重合,能很好地解决坝踵及坝趾区的应力集中,可为数值求解重力坝实际问题提供一种高效分析方法。     

夹层板的边界单元法

本文提出了Ｒｅｉｓｓｎｅｒ型和Ｈｏｆｆ型夹层板的直接边界单元法。从夹层板弯曲问题的基本方程出发，利用偏微分算子导出了夹层板偏微分方程的基本解，由此建立起问题的边界积分方程组，采用边界单元法求解，所得结果具有较高精度。     

求解弹性力学问题的有理单元法

传统的位移有限元法采用多项式形式的位移试函数,对于边数大于4的多边形单元,构造满足单元间协调性要求的多项式形式位移插值函数是一件困难的工作。本文利用逆距离权插值的思想并考虑到单元节点的分布,建立了边数大于4多边形单元上的有理函数形式的形函数。利用有理试函数,采用Galerkin法推导出求解平面弹性力学问题的有理单元法。采用有理单元法求解弹性力学问题,求解区域根据需要可以划分为任意多边形单元,极大地提高了网格划分的灵活性。有理单元法不依赖等参变换,不同单元的形函数表达形式统一,方便计算程序的编写。     

高速公路沉降的非线性等阶径向点插值法解

无单元法一个突出的优点在于其只需要结点信息而不需要单元信息。先介绍等阶径向点插值法这种新型无单元的形函数构造思路，接着给出了它非线性求解平面比奥固结问题的主要方程，然后对一软基高速公路的断面沉降进行了计算，并与非线性有限元法结果进行了对比。可以看出该法不但计算精度高，而且在解路堤分级施工的这类移动边界问题的沉降时，比有限元法更方便，具有较好的应用前景。     

阶梯式Timoshenko梁自由振动的DCE解

本文基于微分容积法和区域叠加技术提出了微分容积单元法（Differential Cubature Element method，以下简称DCE方法），并用之求解阶梯式变截面Timoshenko梁的自由振动问题。根据梁的变截面情况将其划分为几个单元，在每个单元内应用微分容积法将梁的控制微分方程和边界约束方程离散成为一组关于该单元内配点位移的线性代数方程组，将这些方程组写在一起并在各单元之间应用连续性条件和平衡条件得到一组关于整个域内各点位移的齐次线性代数方程组，这是一广义特征值问题，由子空间迭代法求解该特征问题便可求得系统的自振动频率。数值算例表明，本方法能稳定收敛、并有较高的数值精度和计算效率。     

轴对称弹性体扭转问题的无网格自然邻接点Petrov-Galerkin法

本文将无网格自然邻接点Petrov-Galerkin 法应用于轴对称弹性体扭转问题的求解．无网格自然邻接点Petrov-Galerkin 法采用自然邻接点插值构造试函数，并且采用三角形线性单元的形函数作为加权残值法的加权函数．自然邻接点插值构造的试函数满足Kronecker delta 函数性质，因此本质边界条件的施加十分方便．由于几何形状和边界条件的轴对称特点，原来的空间问题简化为二维问题求解，因此计算时只需要横截面上离散节点的信息．数值算例结果表明，所提出的方法对求解轴对称弹性体扭转问题是行之有效的．     

TSP问题的单元划分法

提出一种利用单元划分法求解ＴＳＰ问题的新求法，该法计算量小，且计算结果理想     

动态断裂力学的无限相似边界元法

对弹性动力学的相似边界元法进行了进一步研究,推导了相应的计算公式,并在此基础上提出了动态断裂力学的无限相似边界元法.与传统的边界元法相比,相似边界元法由于只需在少数单元上进行数值积分,大大减少了计算量.对动态断裂力学问题,无限相似边界元法由于在裂纹尖端的边界上设置了逼近于裂纹尖端的无限个相似边界单元,可直接得到裂纹尖端具有奇异性的应力,而不需要设置奇异单元,从而突破了奇异单元对应力奇异性阶次的局限.另外,还讨论了无限相似边界元法得到的无限阶的线性代数方程组的求解方法.     

SH波入射时垂直半空间中双相介质界面附近圆孔的动力分析

利用复变函数和Green函数法研究了垂直半空间中双相介质界面附近圆孔对SH波的散射与动应力集中问题。该问题的解答采用镜像法,首先构造出含有圆孔的直角平面区域出平面问题的Green函数,然后利用契合技术,并根据界面处位移连续性条件将解答归结为具有弱奇异性的第一类Fredholm积分方程组的求解,结合散射波的衰减特性,直接离散该方程组,把积分方程组转化为线性代数方程组可得到该问题的数值结果。最后,通过算例分析了圆孔的动应力集中情况。结果表明,与全空间中界面附近圆孔对SH波的散射相比,由于垂直半空间自由边界的存在,孔边动应力集中系数明显增大;另外,入射波由硬介质(波速大)进入到软介质(波速小)时,与均匀介质相比,孔边动应力集中更显著,最不利的参数组合,孔边动应力集中系数几乎提高了一倍,入射波由软介质进入到硬介质时,情况相反。     

极限下限分析的正交基无单元Galerkin法

基于极限分析的下限定理，建立了用正交基无单元Galerkin法进行理想弹塑性结构极 限分析的整套求解算法．下限分析所需的虚拟弹性应力场可由正交基无单元Galerkin法直接 得到，所需的自平衡应力场由一组带有待定系数的自平衡应力场基矢量的线性组合进行模 拟．这些自平衡应力场基矢量可由弹塑性增量分析中的平衡迭代得到．通过对自平衡应力场 子空间的不断修正，整个问题的求解将化为一系列非线性数学规划子问题，并通过复合形法 进行求解．算例表明该方法有效地克服了维数障碍问题，使计算效率得到了充分的提高，是 切实可行的．     

边坡弹塑性稳定性分析的无网格法

将无单元伽辽金法(EFGM) 推广到求解弹-- 黏塑性问题, 自行编制了相应的计算程序, 并应用于土质边坡的弹塑性稳定分析. 通过无单元伽辽金法计算结果与有限元法的计算结果的对比分析, 可以发现用黏弹塑性问题的无单元伽辽金法计算程序去求解弹塑性问题是方便可行的, 本文编制的计算程序稳定性好, 收敛速度快.