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

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

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

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

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

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

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

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

三维位势问题边界元法中几乎奇异积分的正则化

采用一种半解析正则化算法,计算了三维位势问题边界元法中近边界点的几乎强奇异和几乎超奇异面积分.该算法适用于三角形线性等参元.对高次单元将其细分为几个三节点三角形单元即可应用该算法.由于几乎奇异性,与内点邻近的单元上的积分,采用半解析正则化积分算法计算;而远处单元的积分仍保持常规高斯积分.对三维热传导算例,计算了近边界点的温度和热流.数值结果证明了该算法的有效性和精确性.     

三维边界元法高阶元几乎奇异积分半解析法

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.      

功能梯度矩形板问题三维理论的半解析解

针对功能梯度材料矩形板问题,基于三维弹性理论,将位移和应力分量作为基本变量,通过双三角级数将其控制微分方程转化为常微分方程组的边值问题。采用插值矩阵法直接对常微分方程组边值问题进行求解,得到了功能梯度材料矩形板三维位移、应力场的半解析解。通过算例给出了材料参数按指数形式和幂函数形式变化情况下的功能梯度板的弯曲问题。对比有限元法和状态空间法,结果表明：本文提出的状态空间与插值矩阵法相结合的半解析法能有效地分析材料参数按任意形式连续变化的功能梯度矩形板问题,且具有良好的精度,精度可达10-4量级,能够满足工程需要;与其他方法相比,本文方法具有实施便捷、计算量小等优点,根据其力学场分析结果可设计出满足各种不同需求的功能梯度材料。     

复合材料切口应力奇性指数计算

提出一种计算广义平面应交状态下复合材料切口应力奇性指数的新方法.在切口尖端的位移幂级数渐近展开式被引入正交各向异性材料的物理方程后,将用位移表示的应力分量代入切口端部柱状邻域的线弹性理论控制方程,切口应力奇性指数的计算被转化为常微分方程组特征值的求解.采用插值矩阵法求解该常微分方程组,可一次性地获取切口尖端多阶应力奇性指数.本法适合平面和反平面应力场耦合或解耦的情形,并可退化计算裂纹或各向同性材料切口的应力奇性指数.算例表明,所提方法对分析复合材料切口应力奇性指数是一种准确有效的手段.     

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

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