功能梯度碳纳米管增强复合材料板弯曲和模态的广义有限差分法

由复合材料构成的板结构一直以来受到很大关注,其中功能梯度碳纳米管增强复合材料(functionally graded carbon nanotube-reinforced composite,FG-CNTRC)具有异常优越的力学性能,使得诸多学者展开了对功能梯度碳纳米管增强复合材料板结构力学行为的研究.本文以FG-CN...     

含肿瘤皮肤组织传热分析的广义有限差分法

生物传热分析在低温外科手术、肿瘤热疗、病热诊断等临床医学治疗和诊断中有着广泛的应用.由于健康皮肤组织内肿瘤的存在使得肿瘤附近区域的温度会明显升高,这一特性常被用于检测皮肤组织内的肿瘤生长,因此有必要开展生物传热数值分析的研究.本文以含肿瘤的皮肤组织为研究对象,将一种新型区域型无网格配点法——广义有限差分法应用于能描述含肿瘤皮肤组织传热过程的Pennes方程求解.广义有限差分法利用泰勒展开式与移动最小二乘法将计算区域内的每个离散点上的物理量导数表示成其与邻近点物理量及权重系数的线性组合,进而构建得到仅含各离散点未知物理量的线性方程组.该方法不仅具有无需划分网格、避免数值积分等无网格配点法的优点,同时还克服了大多数无网格配点法中插值矩阵高度病态的问题,为此类方法在大规模工程数值计算中的应用提供了可能性.本文首先介绍了模拟含肿瘤皮肤组织传热过程的广义有限差分法离散模型,随后通过不含肿瘤与含规则形状肿瘤的基准算例,检验广义有限差分法的计算精度与收敛性;在此基础上,通过数值模拟研究不同肿瘤形状及肿瘤位置分布对皮肤组织内温度分布的影响.     

边界粒子法结合正则化技术求解Robin反问题

运用一种边界型无网格算法——边界粒子法求解Robin反问题,结合Tikhonov正则化技术消除反问题的不适定性。该方法仅需边界测量数据,计算精度高,特别适用于反问题的求解。数值算例显示该方法在求解Robin反问题上具有很好的稳定性和收敛性。     

基于时间依赖基本解的奇异边界法模拟二维狄利克雷边界标量波方程

奇异边界法是一个半解析边界配点强格式方法,具有无数值积分和无网格、编程容易以及数学简单等优点。本文首次将时间依赖基本解运用于奇异边界法,计算模拟二维标量波方程;结合确定源点强度因子的反插值技术,提出了二维狄利克雷边界标量波方程源点强度因子的一个经验公式;引进了解决波方程基本解G奇异性的一种无奇异积分处理方法。数值实验证明,基于时间依赖基本解的奇异边界法可精确高效地模拟二维狄利克雷边界标量波方程,在计算效率、精度、稳定性和适应性等方面有明显优势。     

奇异边界法求解多连通域的瞬态热传导问题

提出一种基于奇异边界法结合双重互易法的数值模型来求解瞬态热传导问题。奇异边界法属于配点型边界无网格方法,相对于网格方法,其具有无需划分网格,只需边界配点的优势。运用差分格式来处理热传导方程中的时间变量,将原热传导方程化为非齐次修正Helmholtz方程。修正Helmholtz方程的解由齐次解和特解两部分组成,齐次解通过奇异边界法求出,特解由双重互易法求出,源项由径向基函数近似。通过数值算例检验了本文数值模型的精度及有效性;算例结果表明,该数值模型计算精度较高,误差基本都在1%以内,具有很好的稳定性,能有效地应用于求解多连通域的瞬态热传导问题。     

物理信息依赖核函数配点法的研究进展

在过去几十年里, 尽管有限元等传统计算方法已被成功用于众多科学与工程领域, 但是其在数值模拟无限域波传播、大尺寸比结构、工程反演和移动边界问题时仍面临计算量大、计算效率低、网格生成困难等计算难题. 本文介绍一类基于物理信息依赖核函数的无网格配点法及其在上述难点问题中的应用. 物理信息依赖核函数配点法的关键在于构建能反映问题微分控制方程物理信息的基函数. 基于这些物理信息依赖核函数, 该方法无需/仅需少量配点对所求微分控制方程进行离散, 即可有效提高计算效率. 本文首先介绍满足常见齐次微分方程的基本解、调和函数、径向Trefftz函数以及T完备函数等典型物理信息依赖核函数. 接着依次介绍非齐次、非均质、非稳态以及隐式微分方程构造物理信息依赖核函数的方法. 随后, 根据所求问题特点, 选用全域配点或局部配点技术, 建立相应的物理信息依赖核函数配点法. 最后, 通过几个典型算例验证所提物理信息依赖核函数配点法的有效性.      

基于二元连分式的散乱数据递推插值格式

本文首先基于新的非张量积型偏逆差商递推算法,分别构造奇数与偶数个插值节点上的二元连分式散乱数据插值格式,进而得到被插函数与二元连分式间的恒等式.接着,利用连分式三项递推关系式,提出特征定理来研究插值连分式的分子分母次数.然后,数值算例表明新的递推格式可行有效,同时,通过比较二元Thiele型插值连分式的分子分母次数,发现新的二元插值连分式的分子分母次数较低,这主要归功于节省了冗余的插值节点. 最后,计算此有理函数插值所需要的四则运算次数少于计算径向基函数插值.     

广义有限差分法求解Kirchhoff和Winkler薄板弯曲问题

本文将广义有限差分法用于数值计算Kirchhoff板和Winkler板的弯曲问题.广义有限差分法是基于最小二乘原理的一种区域型无网格方法. 相比于传统的网格类数值解法,广义有限差分法无需网格生成且无需数值积分.通过数值实验结果表明,广义有限差分法可以有效地求解两类薄板在不同横向荷载作用下的弯曲问题.     

奇异边界法分析含水下障碍物水域中的水波传播问题

研究奇异边界法模拟水波在含水下障碍物水域的传播过程.奇异边界法是一种最近提出的新型边界配点方法,具有无网格和无数值积分、数学简单、编程容易等优点.首先研究了奇异边界法分析典型水波算例的精度及效率,并与边界元法的计算结果进行比较,然后通过数值模拟讨论分析了水下障碍物位置、尺寸及形状等因素对水波传播的影响.发现奇异边界法的计算精度较高,且与边界元法的计算结果吻合较好;数值结果显示水下障碍物的不同高宽比对水波的传播影响明显:障碍物无量纲高度越大对水波的屏障作用越明显;障碍物无量纲宽度增加对水波的屏障作用先增强后变弱.在高宽比一定时,斜率变化对水波的屏障作用不明显;含吸收边界水下障碍物可以得到较低的传递系数和较高的反射系数, 对水波的屏障作用更为明显.     

A Novel Localized Meshless Collocation Method for Numerical Simulation of Flume Dynamic Characteristics北大核心CSCD

局部边界节点法是一种基于非奇异半解析基函数和移动最小二乘原理的新型无网格配点技术,该方法把每个节点处的未知变量表示为该点对应的局部子域内节点处物理量的线性组合,该文基于局部边界节点法对数值波浪水槽进行了研究。首先,通过基准算例确定了Laplace算子非奇异半解析基函数的合理形状参数值。进一步,基于合理的参数选取,用较少的离散节点即可成功模拟波浪传播行为,将得到的数值结果与其他文献数值结果比较,可以发现局部边界节点法用更少的局部点即可得到较好的数值结果。最后,以保护近海岸建筑物为目标,模拟了水下防波堤对波浪传播的影响。结果表明,当波浪与梯形防波堤发生作用后,波峰变得比较陡峭,而波谷变得相对比较平坦,为近海岸防波堤的相关研究和设计提供了数值参考。