边界元法分析薄形层合结构

薄形层合结构由于几何形状的特殊性,其力学分析是数值计算的难点.边界元法分析层合结构具有较大的优势,但对于薄形层合结构,边界源点和对边上的积分单元距离很近,边界积分方程中存在几乎奇异积分,常规的数值积分方法已经失效.文章引入一种半解析化方法,计算薄形层合结构边界元法中的几乎奇异积分,使边界元法能成功分析三维薄形层合结构的层间界面应力和各层内点力学参量.     

正交各向异性位势问题边界元法中几乎奇异积分的解析算法

几乎奇异积分的计算困难阻碍了边界元法的工程应用。本文针对二维正交各向异性位势问题边界元法中近边界点的几乎奇异积分,采用分部积分法,导出一种直接的解析计算公式。该解析公式可以精确计算线性单元上的几乎奇异积分。对二次单元,可将其细分为几个线性元,采用该解析公式近似计算其边界积分。当内点离当前积分单元较远时,仍保持常规高斯数值积分模式;而当内点离其较近时,因常规高斯积分结果失效,则采用该解析积分取代高斯数值积分。数值算例证明了该算法的有效性和精确性。二次元计算结果比线性元计算结果更精确。     

插值矩阵法分析正交各向异性板切口应力奇异性

基于切口尖端附近区域位移场渐近展开,提出了分析正交各向异性复合材料板切口奇异性的新方法.将位移场的渐近展开式的典型项代入弹性板的基本方程,得到关于正交各向异性板切口奇异性指数的一组非线性常微分方程的特征值问题;再采用变量代换法,将非线性特征问题转化为线性特征问题,用插值矩阵法求解获得的正交各向异性板切口若干阶应力奇异性指数和相应特征函数.该法可由相应的特征角函数对板切口的平面应力和反平面奇异特征值加以区分,并将计算结果与现有结果对照,表明了该文方法的有效性.     

二维位势边界元法高阶单元几乎奇异积分半解析算法

准确计算几乎奇异积分是边界元法难题之一。目前,对于一般的高阶单元的几乎奇异积分尚缺乏通用高效的计算方法。本文在单元局部坐标系中表征了二维高阶单元的几何特征,提出了源点相对高阶单元的接近度概念。针对二维位势边界元法的3节点二次等参单元,构造出与单元积分核具有相同几乎奇异性的近似奇异核函数。从二维位势几乎奇异积分单元积分核中扣除近似奇异核函数,把几乎奇异积分项转换为规则积分和奇异积分两部分之和,规则积分部分用常规Gauss数值积分计算,奇异积分部分由导出的解析公式计算,从而建立了二维位势问题高阶单元几乎强奇异和超奇异积分的半解析算法。算例结果表明了本文半解析算法的有效性和计算精度。     

二维正交各向异性位势问题的高阶单元快速多极边界元法

研制了一种适用于二维正交各向异性位势问题的高阶单元(线性单元和二次单元)快速多极边界元法. 在快速多极边界元法中, 源点对于远场区域的积分采用快速多极展开式计算, 而对于近场区域的积分则直接进行计算. 高阶单元的使用使得近场积分, 尤其是奇异积分和几乎奇异积分的计算更加复杂. 通过引入复数表达对其进行简化, 若边界采用线性单元插值, 近场积分可直接解析计算; 若采用二次单元插值, 则给出一个半解析算法计算近场积分. 高阶单元奇异积分和几乎奇异积分计算难题的解决, 使得高阶单元快速多极边界元法不仅能够计算一般结构, 也能被应用于超薄体结构, 拓宽了高阶单元快速多极边界元法的适用范围. 数值算例表明, 若计算精度一定, 高阶单元快速多极边界元法较常值单元快速多极边界元法使用的单元数量显著减少, 且高阶单元快速多极边界元法计算时间与自由度数量成线性关系, 其计算效率仍处于$O(N)$量级, 因此高阶单元快速多极边界元法可更加高效求解大规模问题.      

非线性多点边值问题的插值矩阵法

建立了插值矩阵法的基本理论，用于解非线性混合阶常微分方程组多点边值问题，制作了常微分方程组求解器ＩＶＭＭＳ，可以支持计算力学中的有限元线法。     

三维声场边界元法几乎奇异积分问题分析

以三维声场问题为例,提出一种准确计算高阶单元几乎奇异积分的半解析算法.首先分析高阶单元几何特征,构造近似几何量,然后应用扣除法,将奇异积分核函数分解为规则核函数与近似几何量表达的奇异核函数.规则核函数积分采用常规Gauss数值积分计算,奇异核函数积分采用半解析算法计算.给出三维声场内问题和外问题经典算例,计算了近边界点的声压,结果证明本文半解析算法的有效性和准确性.     

功能梯度变截面梁自由振动和稳定性研究

基于忽略了梁截面剪切变形和转动惯量效应的Euler-Bernoulli梁理论,研究了轴向力作用下轴向功能梯度变截面梁的横向自由振动问题,将轴向功能梯度Euler-Bernoulli梁自由振动固有频率和临界荷载的计算转化为变系数常微分方程特征值问题。运用插值矩阵法可一次性计算出轴向功能梯度变截面梁各阶振动固有频率和临界荷载,分析了轴向荷载对轴向功能梯度Euler-Bernoulli梁自由振动固有频率的影响,即轴向压力使梁的第1阶固有频率降低,轴向拉力使梁的第1阶固有频率增大。在简支-简支梁(H-H)边界条件下、不同截面宽锥度系数c_b和截面高锥度系数c_h,且区间划分点数n为40时,本文计算结果与已有文献计算结果之间的最大相对误差不超过0.00768%;在简支-简支梁(H-H)、固端-自由梁(C-F)、固端-固端梁(C-C)这三种不同边界条件下,不同c_b和c_h,且n为40时,最大相对误差不超过0.101%,说明了本文方法的有效性和良好的计算精度。     

正交各向异性材料切口尖端热弹奇性特征分析

研究正交各向异性平面V形切口,计算其热弹奇性特征.通过引入切口尖端物理场的渐近级数展开式,将应力和热流平衡方程转化为关于奇性指数的特征常微分方程组,采用插值矩阵法求解,获取切口尖端的热流、应力奇性指数和对应的特征角函数.算例表明,该法精度高适应性强.